{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Tutorial 6: NearestKernel" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "This tutorial demonstrates how to use Monte Carlo data directly to create a 2D kernel that can be used for PSF modeling in SPECT. In this case, we'll use the same Ac225 data as from tutorial 5." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "import torch\n", "from spectpsftoolbox.operator2d import NearestKernelOperator\n", "from spectpsftoolbox.simind_io import get_projections, get_source_detector_distances\n", "import pytomography\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's open the headerpaths for all the PSF kernels and load the projection data. The data was acquired at a 128x128 matrix size and 0.48 cm resolution. The point source was shifed upwards 0.24cm in x and y so that it corresponds to the center of a pixel. We thus index [1:,1:] after loading the data; that gives us an odd kernel size" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "headerpaths = [f'/disk1/ac225/440keV_PSF128_shift/point_position{i}_tot_w1.h00' for i in range(538,1638)]\n", "respaths = [f'/disk1/ac225/440keV_PSF128_shift/point_position{i}.res' for i in range(538,1638)]\n", "projectionss_data = get_projections(headerpaths).to(pytomography.device)[:,1:,1:]\n", "distances = get_source_detector_distances(respaths).to(pytomography.device)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We'll normalize the projection data so thatthe PSF sums to 1:" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "projectionss_data /= projectionss_data[0:100].sum().reshape(-1,1,1) / 100" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now define the nearest kernel operator, which takes the PSF kernel closest to the corresponding image plane. We simulated over 1000 positions, so the rounding will be very close:" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "psf_operator = NearestKernelOperator(\n", " psf_data = projectionss_data,\n", " distances = distances,\n", " dr0 = (0.48,0.48),\n", ")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let's test the predictions at some various distances and pixel spacing" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "# Input size\n", "Nx = Ny = 182\n", "Nz = 128\n", "distances_test = torch.linspace(0,55,Nx).to(pytomography.device)\n", "dr = (0.48, 0.48)\n", "\n", "# Manually create kernel size\n", "Nk = Nz - 1 # size of the kernel\n", "x_k = torch.arange(-(Nk-1)/2, (Nk+1)/2, 1).to(pytomography.device) * dr[0]\n", "xv_k, yv_k = torch.meshgrid(x_k, x_k, indexing='xy')\n", "\n", "# Create input with two point sources at x=64 index\n", "input = torch.zeros(Nx, Ny, Nz).to(pytomography.device)\n", "input[64,91,64] = 1\n", "input[64,100,80] = 1\n", "\n", "# In this case, the kernel size is as large as the input\n", "output = psf_operator(input, xv_k, yv_k, distances_test)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Lets look at the output:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAVwAAAD7CAYAAAAmeCzOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8fJSN1AAAACXBIWXMAAA9hAAAPYQGoP6dpAADEJElEQVR4nOz9yY8t2Zbeif3W3tacc7y5TTSvz0x2mURJoorkqAiwQIEFNQNNChBQA4000P+hP0FTAQI0lkYFCNKoBAIFqEpFskiVUBLJJJMv870XL1509153P401e28Ndmt27Lj7vXFfvPC4ZwVuuPs51tu2ZWt/61vfEuec42xnO9vZzvZ7N/WHPoCzne1sZ/tQ7Oxwz3a2s53tO7Kzwz3b2c52tu/Izg73bGc729m+Izs73LOd7Wxn+47s7HDPdrazne07srPDPdvZzna278jODvdsZzvb2b4jOzvcs53tbGf7jqx67IIij170bGc729k+OHNufHCZc4R7trM9SZM/9AGc7R3sHLZ+7y0+WFHyQgO2+Puh5b/vtnS8f8hzEHwcYv4A+77P4jVR4S/B4fDHOf3Ojw94OmPgw7FzhPskLD445YMlxb8leyoRkOPYMSx9dp+973O1Dy/yvbGncp/PBucI9wnY1PFIesB0iHAspyPCpxbtvou9b4fzfb9W/mWQj1LxtF4QH7adHe6TsBzZuvCoySTatUfLZfu+TY3/kPa+X0DyHrf10PbV7BsB8Z9lgdX7oKazfR/s7HCfiEmKaL25xQfrfSFE30Vk/EOJvt+H073vWvjtyzySL1hDgg3j4alfyx++nR3ukzCHW4xUVUie5ERKXNrbfXDD3ErH8S4P7ts6nvuSfu/TcZSOaglymSfr3maZ+yCct3mhnFombkeYz2ycG4rlngqk8PueEXz/7exwn5wdD9h59OOdsGH64N9n37fEy9s8mI89x6dsSzOXp+Jk5/ZhO90zS+HJWJwyTpkJy9DC0nrfxh7j0N7XQzRnZMyPY+nzU1Px+N2pa7DEjpjvY2mZU1bu59ted8FTAMHj8OGflI+sIj/CT8GJnWGPs8N9sibcH+VE6thTtD/kQ/m+930fde8PsZ2z/SHt7HD/IPZtHp55lGCnUa4oEBVghvL2PmafbxNB/r7tVET6NlDD+9z3u1oZhb6rCSI1Iu3sc0tmoZyd8VOwM4b7B7H3McVXxV/Fw+bsbNr5tg/iU5ry/T4SbO/7/N8daxXR4Mr7Z/39Bc5VZU/Tzg73SZriiCZUmisf8vLBfAzeW9pjMu0PlRqfsm9DC3tK0dy7OUJBe8GocKrOjfm+ispXwMUE6dmegp0d7pMzj936x3h5qprr7C3HsEIslCgdwSlHu/Td3JYe9u+CY/tYWtnb0ML+0FFiTJKFajJnkTRbUWHm4qPcfH/hD3/cZ3usnTHcP4i9jwRIFC4pHZ6d/e5YntI+9rbf54i/rZ3CYx9zbR7DKY4/FffjqN9VtPzweYlo/y86XjfiQlQropLzzaI156z/U7NzhPsHsfeV0JHw/6yvEPFbceAw/uEVlR5cbw/hipEeNYcLTldCnT7Gh+w+HYiHHP5DVK9TL5z71n/Ifn+4sXNmUrLrIQTr7517qrzbs5V2drg/JCunnxFOkHefxHho4nFLfneR1jxZ9Bh7n9oJ79Nk9nP+tUqO9lhD42xP0c4O90lbOcUESUmVYpHAWkjTVFEh0fKQOr17xDJ+uWN7SBug/O4xBQnl54aMdd633VOFE/Nl7jvGuOwp7HtpvYcSj/NlZi9EUUV3FQUS7mnsJnCOdJ+0nTHcJ23Th2/Ox51MTyVHuw9Xp/2+7X2L7MDUoc2d7dtEhaeW/TaO7v7rLcV/9+7z7GyfvIlz7lFP37mn2XdhmneTU5yKpniSfIVIFbBbW/wcJstO7V0d8akIcB7dlZHpfdGinulBzLdXOux5pn4uIPNQZHoqMv621+IUe6M8nojB6wl8EM3h0kvTubi9t6Xg3Xc8Z3uf9pieZmcv+r2y9xHBqEAbsmF7NtGL/Kt1yVk9ZO+C0ZbO5dT+TlG2LO69RMEPJc1k9vvbwiOn7KHrWmDsk13Fx9EeHd3xMb1Nxd05Mv6+2Nnhfq/s2071iwjSjd7pRgxXGqIDjp95B1w+jPfhoPc5njkeOy8rXlp+Hq2eisDmcoh56u0Wo985VnoKO116IZSsjPIc7ML69x3j/Lvjz6W4Vw6HOIuuLsFZrOvDNypETbbgVr9Lh4e3jYj/0JDTD9fODvd7ae866GNGu6hEE89BndDCpPKcTmffYi8PUbRKhxmTW49JSpXdKtzC58sOZkr8v++4liwuU8IRpRhQ+f1SmfSS4y63N1/m2BJcUH6WpqT+nKOzzYUO8338PuzsbH+fdk6aPTlbcignJPqcdxhyJOkXfooiR4FvwxZY2veSneLnqtlnj1mu/DY623Ldch3FaUdfFkOU68VtnVq/3P60sMIXKzwU1c+tLNMNEXsodMj3a15VduoF8z4Kac72Xdg5afaDtOwIssOdYrcTgH+xVHQpcltKtM2TUg9NvUs7BQksLXefQ3loCN+XPFuigJ3a/hyGuO/6POx8o6MVtfJrLyVd3MiyVsL7SvKd7X3ZY5Jm5wj3B2mnH77kfEvVKVkWw3l3kr1wHHmVDvxURL1k7xp5l/aA2M+9EffSvtRsm0uR9mkTJLwMw9+pbNdSRrLTmcnSNTtHtU/NzhHuk7OlSHA5WSUISm8QVErEKGlwWKztKRNoQEE9MixHhad4r0s2pT09jEHeV0E2Z1acOve3pX49Zv+nEm2PjcRnMpqpsCHitNP9Sfje2n6hb9k5iv0+2znCPdvRAx0tRbqLpb/3VWkd7aH4vYQuplPraVv3+T6/zTCcR9GnvnsXW4I85toSS8ssWQlDKFRwumUgk8RpzgUOP1g7h61Pzu6b7kZ2QEE5ciP2CHP0JkfUsMj/1GFrjvsjz+PjyF2E8zaXObVTRzxNhD1QcXUU5c/pW/ftY2mb86TiUlR/XwFFeRxLtDLSeko1vigFhQpMkcSNCCJDzvUF5DOn7T0Guz3jut9XOzvcJ2NLDxssP1SeneCpR7PPRSEOCFVoETt0rqBApaqnBXL+yWOLR1Nm85ecxtS+XYFDSA6yJML90HZP8Vkfii5PO7OyZf2p9ZyzCBZRVfDLCkkvOkXiTrP08prb2zrU+xz2Q2PqbO/Dzg73SVn5wJQP9ryv1Zyj6gn0In4qG6vQJCye+J4TO+Us4/415TR5Woqr0+fHxzd/mO/DdO9zfjOusZvu27en8eyLJc7r9PgeOhaIM4epc585rdSJQQV2gQtLSTou53qQykMK4h1wWaAyRSgU4mzhdJew64dsft7F9YECJ86zgels42zv084O9w9m71LcUGKHD5Ws5ml2dqje4YryU1frRqwbg+4CWBuXjw//0hR8yXGVsMGpooZ43NP1M3wxd2Knk2fRKaRo2hmOnGeBgy472zmkcF/Z8WP+niYdj8RoAmbu2+bkQhSH9ZGtFEmXWJTivq3+wanzLrUZptdtOkM56y+8bzs73D+IzZM8j3G8p4j4SzanQcUICrRqU3JmNHusG70TRuFkxDkQIf1Me0tY76kpc/4uR0mnWAEFthqTRhMHsBRBymQdSfzUEEVKE45z6iQeVOBKn5061ulLYooDl456GhW6EA2XSUmRCq1WReeGgv7lwLo486h8JJzYIsDEET4EeyzBBiXGvIRfPxY+esjO+PF9dna435mVEe3bTA0fg63NnYRJ6K0EnDA++IPZ+oSN8rdeB0cFhGSODVFv4WRd+ZAuR57TqWicemeM1R/9NGJ28fgmtrT9ue5AWQIbv8i4p/9l+t0ytryU7JrbMdtiitOe2l441lB4EtXbIqRjkl5CPL6iGCXBEQ9F5vPjPEVhE7Ii2anrcx/d7b59vE0gcLYzLewPZm9D/o926kG4r1S1XAas7T2UYP3DrVSVIi4VnMJxZVp5zG9jMQlXcVxc4Y8vQwPl/uaww3R6vtT9YAIxpONXR8s/rphj7lCnbIfj7UyXL7+TAkpQkq81sW1OcrYlq+MhyOhtbHYu6eVb/DuiBj5mbJ6LLt7FzhHud2bv+uY/VV1URiLLpaUpweMMSq0zfhgiLE2DiKKuLlBSoaTGuoHb/a+KbSoQQo+0U7Sq8rOcpEqOZ4JNGmIUnM9IUvSVsdwyKcfCsjYnw+ZasgE+ib/n6DtsQ+qw/lhcpxjhRedzXLDg0rGfoptJWlYVxQ0iikqtcM5ibI+SikpvsLFJpBtDAq28jvNE5EN2ik7nnarDgRsWkmXLlMGH9/E24zne6/eFST/dKPoc4X6v7b7Musz+Pl5m3tnhVDNCP+3Vs1JSjv5+nN0T+bjpw70MNeTvl4skHmnFuc47YSwdaxSfeXwU/JBlDYsMK8yFhDyOO4lwTxY9vEsRx8K1O9lM9LGshKfr7L4Pdo5wvxc2744LGYMzs8/mEck84TP9XAhRre1BVIh5qlThlFgCzrAfXmPtGKKz4AwCCb+clk+ixdQpWKVOwX7deTFELMrweOIEAih4wLkM2DLvhJCTUApRlUcQbT9JlB1JHk6cbTivEN3FSDd9XVDJ5tuQ4mVx/KIoCkbC8akwm6hUU0A2NXV1wWj2jOaAsYcpJW/yQiqlM+edMk7ZdBz48xuZ9qYrH3n/ovXX722iz7d1uu+L7fD0nf05wv2922OiklMP0bzK6LHrHZufUvrEjI3i5CfMZ9KrFJn5+n+dpqNHEWCRVJMHncM0ui2XSVoD5YsjUtQmkWEREUqOSucvg9IS/irz/U/P49T65efLyxw/SrLwmTE91sbrfzydPx3tl8vflzw9RcvLXODJ9o4i6jM2+/u0c4T7rextubRzGGAu2s09fy/Z0sNZPkAhWgzYo4MAFFqsBaUrlAqFEDHfFDBdY/rAaJg6BUcsO/V43qSoIFG8BnLkXbIWCqrUBFvNWXSRykdkzky3nU45YLeJuhYLCGKG/3i9SXIt7Pf4Si4VZKiTf8/pYROK24IJGusGRns42k8Us3ElcyHsJePHSzOg+x1/jPz9ta1C6fA4cbI++n0XuKI8Dk4cy9nmdna438reFfN6mwz0qQfBcjzYywhITx3Nqa3YMeOLaJSq/V5T4gccKiR4yqnp6YjbT81jwicyFVReY4KtlnzVXBBQTtuliGaB8LdK2gNlpLaEC2c4ZHbdjzDe8uUy5xSXZz5N3E2+m1HAlKpS1C+i0NKkaxkx9RjtZnZCxumncMYpW4pqsx5yOt1JMUipl3EfD/khOzvat7Gzw/3ObamC574p4lI0PP/+xFSzpHe5KdUoOivrBj/1ddOiiLiMosKJwprc7mV6BGV1RHSMldcKsGXEplKENWUMaKbOszj2mUOU4rvo2Jw7TDHmgMMuZvqdzLZ32rnMI+TJlY3RImMovc2CQdHZ6oDdKvEvMB9Za7TSOGsTx3nKKy4LHU5F06ciyuNxIIljnbs2x+uXXihuPEKsl67G6e/O9jZ2drjvzR47KN+WX7lcCJAdySkMrkxEFVQxqQGLNTuc8pGWUk3iiPqiiBqvoTtMJAOtU7kNTAqICqyVEJGmQ6nCcZjwYJcVZrNsfRmpJocwn/7P+bFwDKMQounixbFIaVPp+HIEqY/Oh5RcrMgFCdFhhvN3LsEhElgJSurJS6R8wfm/gzwjtoBi4ktxKUHnj6I4yXjl0jlN8HOZXhcpIIVIRZvyfZdePm72c8nOzvht7Jw0+07tXfDeEyT05LDuw/GCcygiulxWanICzU0jXz8Vrj3EIDVatYGn2+TCiKNE1sIhFpFUTNrlYy8pXyWmeILCFeGEBWpVOrcj55yJ/YsdLeYFHkfHlZeTe843UcpSgvEUBcxgC1hGinOaHveSzaly5XJFwUdgRwBH97V8CZ0WqHkXl3B2I4+1c4T73uyxzvQx0UK53IkIopBTzHSmESZRC5TTagc4u5/tM/NAXdRbkBqFwlWb0CDRR2ej2WPsLsg6Ru2CaaLHQwQ2LSNCwHNjQQSkqLToNpGONV2Cwjm7pWWLafKsTDX/HaJr0ZOEVnQ8JdMhisX4/avJ8qk4YbKPfMy6ukzbrPQqRbfxb4elH24TO8SVGG6k7FHAGNKmYz/dKaOIyst1IUAttrhGZVLRztYvk2+Rdvg2gUGk+6W9L6x/joKjnV9NT9JO05emdh+Jf/nWl1GRkholFTrCDSkhNIve5kmdeYHFLMJ8t4KKU8d+z7YW2tfkdZaj5IfsWJBGJTjGiwDppFORk5G+ui/KMKaX26nI/qH9JztNcTu6xieKXmYLPeoYzvbudo5wj+y7eBsv7WP+0MyXOe0McplmrpbKmwkPf8AxU/NCZ0MmvUkOIdpou1CO2hK5uLH0t61e4DD04w0xCpscdZm5d5EB0YQEky+lzQuXSb0igiwpTPPlmEEVcyqVFNubT6PD/jMuHf+umEbFBa84RrrSLPYfi8wKES8I5FkJKmhUxKISm6PqhNEeO7983QZSo8mEjefzy9e4wG3jTCBG6ChU6gZsQaZ6DRkvfhc64n22tP5SxPthOvezw120P8Rg8HzU499h2dmeLirwW5hSofLDXPztLJYxToYxtgMFmhZFjtjAMxgAzLyyq3AK0/2GZI/LQudTbHGW7Jokq2JSaorxzjHX6OimbYJm21vg0y63TY1T7ykcs2yRllYlRyqoxPLIDj1W3XkFtml0m+0ka+KIdlYeT0GbS1BEvr5Z7S2e87KT/5Cd3x/Czg73yL6LwXcfPjvPDD8m4p6W9E62Nod+07fGf+dGrGoQLIOxWDtS6wtENJWscoSorzGqYz98g7M5i66CKAsBm8SZwnH4aax1hcjMEdc1Oq8ml7mGSC1Gl9lZTHFVlyLoikkfsPTZOMN04/4WosuyIECmzR2PnFU4vkptEFEY26NVS60vAkabk2NesKbDuj4d82Q7pYUXSmYTlCW3Ml3OTTnOkxeYG6csihIzp0yYPWQlRW2Jq/uuz8qH6+DPDvc7t3lEcer3pfXiMksDfv4ATfG9JYywrI5yzqJCdGYCpGDc4D9D+0SaKHrVpOV9x4h+NrWe70TN47YFmzm0kxjvNFpdlJGcfJZZA/GYc6KxdMJx2dnfZQlx4dBENUTHKKKo9QVa++uipEKrNYPdY+0wYSV4pzcenV9WVcvRfb5fOq8bqvZiwnSynhCSlBmKmUAIs/2VZSEP29tAD+eI+T47O9zfu52CCeJnj7Fj+tRx1FFub5a5ThFRfNjCsYhCqSZlzWPyx1jvRK20OBRaQFFTSUujLxnFO+TR7OnHmGGPUeF0yixFhOXmGG555LNqsvjZfP3EXy2c/PK1YhKl+hLaqZPNugoqXNk5/JA/S9BF4C1LSpTVVHqFiMZhqWVNrTYYN+Bkim+7UCQRm3imwoeE9UZedR4rnjetcK4LG4nwip78jbNo7btJ2MB6cBIZEY+JZku7b5ze9/o86zA8ZGeH+53YYwjkcyuryEpHsJRci8sVU/m0N5WI/UkcxnnYwTmDtQeiQzH2gHUVTVUXW9cJxzVuoFZrtHisUqsWQTPaPcYe0n7LpJTXO1CpjHUaRebfp8tMne4chpjySYuSY8nsgGPc1maGRWQNxEhwFtXm/c6utFSIalKhiFYtWjWsqueTdXp7Fyh0PcZ2KWmWXwCZzpZpchFv12RNYJWi3il8UB6Ud9weuukRl5Nz3vHaBSzYL1Fs5Ohcl5db+vux350Nzg73O7C3HYTz0k638NmSzaPgzN883RKGwD3131s3Lu6hzLBXrEDAKeOdSbXBjSY4XJJje6xNo86lCH5h2SLBlp100U1hYjniTuR/8VKStphy50qsUwmmSP+qEmMj/qxl4x24GAa7Z7TdxNlOepcRseUMLeRilOhgl5KfD0SPEfdN0Xg4h4LRcGLFdH7TZc7O8/dhZ4f7ndopqtfS4J4ric3J5eVfpUhNfHA0U5v2ECu1W8tCBfDNJaOTqtUaFFTSogOsADBywIlFxCeMKrWZCNw4LCKZluS1B1SBM9rJsaUpv2R8tJwOT6lgY3IontJVJrjGtM3c3qa8kjl55OGU3NX4KKmFh1z8/qvcbNNZLANKXaGlYnRdwrqXJBldiD6PLZ+nXy7APqXoTCzBnSfGCtNqg0iV9HXjrKJUXju2TCGcOuT4gj9VwPDQ2H0MhrtU4PNh2H0h09merE0TZqXN+37NxWr8Mh5rjM5OiUKhk8MRFIraV1RFnq4qSn7vOa4YZZ4ufsjLTI/rFMwQnYs6+n16PY73lyhdDxRixHOe789HvHqy3NzKjg6T6HmpMGEpopflsuT5OR2VEb81bvtYewxOe8ZyT9k5wv2D2NtgZmWWOn4/XaaUQXyosizzYw1ZU7aIOmerK6lp5JKaFZoKw4iguJZP6WRHb+8SngkgVjHYOzLuWFLIqsRdzQULMYLNbAmv8WBDciofS6R6zQsY/HrjZDtxvbkTLAs8XHlsofgh4Z5lok8qVPEyWjcfU6s1tVofRbXRwcaEWj/eZLrbgjkXi0cylOCcQUooaFbUUHb/TWyRoOaWaXIj1nZLe4xXkymFcCl6fgxt8dT277MPK6ot7exw34u97TTqsdvLpblu8h3MnVm0RdHuU5YSVIVcYHiAxSqsGpMDsZhJpGcDxLHS1/TchWn2CIwTZ3BUqIDFyXGUenQFJFdwTa+MysmxyXmfYD9Eh56i2YqsQTvFhWPknbnLuZxZSU7IqeLvcj/GDThnioKHciq/xH5g8n1J4Us6uOX5O4vElj3xeGNRxcRhL13XeWI1H3l0vvdTvcrfH1dW/iE71lN2driL9j7Ke+d0sJLGNf9Oz75bxurylr0YiwudFaY4XIzEytLN3FHBFfQwv2BBkg9mnMWpEWUrBrunk7vkL2paFJoDdwiK5/IzdtUbtvIl1o0Y26FDSelod35fRWRqXQ9uySmUUWvGeI+crngR75iMio50MV0W9hnxY986KJYxZyy3dLbiVLo9EjDcGMFXeoVWLY26REsdeMp+2cHu6MxNvoYhcebPf/rymPNjU+uiKDoTSrVLpxtfpDESL18A2INPeIbPImPkmIN9XxS7ZHMucDzmuTNdGt/l32eLdna4i/a+B8l9WeLyu3sw1zkLISV31BHlZ/p3UWvPAuSwgHdGxxAz7VYN3sGJxYb9r7hEuwpNzSAdtVoz6HVRWTWtqIqRaU52TTHkaTcC3yRyEokmSCRP7XE2aoons8HBQmYFTLDpgrdann/CWcv9FTCLiKJWaxp9SSOXKDQjHRaDQmPEY9rx3KclvGp2TP4aHAnYzDm25W0K980zSaokkBMVx47x4/DieBSW+9B4n2/jPkpZSU08O9u5nR3ue7Fvi1uVuNqcJhU63YYOsx7BnW/rGFqIf+eMtz1++BLeV2T8C9aCsQeMWWP1kCAEwwDUXNkXqPCfYaDjjkZdQAV2GNIpJYeVeKhj5sQW3FRLPzmuCEvkUt+ZU4nQRwFzuMKRpe0Hx+2PZZxEmtPy3ZKdwORz57zTbvQla3nBiksADnKHYcBhPYNDtSGy7Y4adcbI1HOAo0zjEnOhPK58LzI/eMTNEoM+Ki/kJwuq3WmWwkNjVvH4ziTlOnNa4tlKOzvcb2XvilOVjnEOM5yISCbkdQ8dTPHd4FyL5QRZzH5Pt5kVs/KDHXmpDdYNdObWV1Ipy0ZeUNPS4IsjxvBQVtKy1i9olNeGHe2eQ++j4qRGNon4gqMNGgOZvVA4wcLZxmUlUctyIiuey+TUCkeftgVTLm55H1LSzjtepZrgsH1TzXX1kkYuWXEZXjrQssFi6Nl7x+uMx7vtiLX9xOHGYzSTZGHQ4M0YRjgG0jKJdREoe0qtkqhQZJKIqADTRCcfHLVdpqLl/c8j1aVS8YcgrhIDLu2M4S7Z2eG+sz2EU92Pwy4vtxRVlDjcQ9uL00jDKQ2FxRbnUj6EpKl85Jwa0zOqzjsg/RGVq9DokE4b8KSxGk0NAoPeIaIY1DZDC46QVIvHEZ1ojDobjnQPZjhtZAxAdpwlDDC1+FKKVWUFFrxAA8vON0MXSvnChkqtadQlDRsqVzFIh8NSu6AORhf2mGGJspINjiPwdO0pEp3zJFlynhmOKBN286IK6/J5usXqsiznOB8V07F8ajwvRcpxa3No4b7l77P3kT/5/trZ4b6TxUFx34A65Yhltsx8naXlIkG9mEZOvg/OuOy5NaeJFdVZc31ZiPhhWEfl5ItTjVcHG8DojlF/yiA1O3aI85DCym3YuEs66Rjo6OQKFKndurHdxCl65+CjLyUNKO9IrBsZzS5Nv/O5hsSYmhZSxG15LYSx0CjIGKxfxouozx3z5O/g/LRaoVRFpdcoqTxmqy5ZyzMMA7fyDRUtNS0H2WIY2NlX9PaOfrxlNPtZ26L8YsnnH6vaFBxBGVmTIjrbWOHmy6lDcYrdHyUOI6UuYcGzl66Pipei2lNJtfjdUo7hbXHfJVuKgn+Yjjba2eG+sz00MO77/jGDMdrUyZ5OkN2zv9LBnkxKTylceQcxYvOSg4aBkYERQyWgQrRb4ae1Gs1WtYwcPA6rpp2Bjw5tRvFKGOUsSTb/fek7Jxw1t3TYTOEq9j+HK6IpVQWNCM+jbdQltazR4VExWFTYr2FgdB3GeczW2nGmDnbPdc17pJxZzKUk5y8WmWgtnLAyIfjWBRBzp3tqwJwhg3exs8N9K3tMOeOcJzv/PK8v4fK7R5ftxhJfXfy9tF5RnZQi2yWBl3ycUxHvGYOgeAAP9g1GBmppqdyGDS1RguWZWlEp4WB2GBmCILdGq4Zh3CUVq3zIZeY+J9ZiIixG2nMcFkCrBq3a4OSGRBUrYYuIaboUJRbwQoAalFRUap0SXbEfmXUjlax5Jj9O23tmP2JNyyt5xV629G7H6DpG22HtMDnWyOf1tDQfyWdBmYjTTmlwDl/wMEl6hTbmlWwSCyJKQlo3YG1ZSh22FR38QqQ7xf1PJV/TDSrGwxLl610j1Hlk/Zjlfxg0s7PDfe92Kokwt2WMdWpLkdFSxFGyG8qP1fLvxfZLjuh95aDO+WTQYPeI0oEi5ohSjwpBi1CJ0LKi5ZJDcGQjfaoyW8ZbiwiuKA9eKpVd/ExUcjyT6jN8hDfJ6MdzxgKh6i288MreZIhvpmkx1C7oSIQXnSUWgwze6bkhzQJOnd/xBT2OrhfPscSsJVTrkVv3JKHziB+fuIf3j7NTpcPR7jun70oX4ek7Wzg73Le0x97wpcjz+OE63ao6RL9LhQop2itqzyTIKU4I9QUdDJh22Z1ig5aYnc8Rl182ZPpdaK0jln68xeqRTn2EpqJjxYqGJnBHrYNn7oraNYyqo5M7zPi1x0OrqxRJRkcxd4Y6dJCIEVv5InAJZ21xeLWyKB9pnC+yqPXGH4eNFW8FPBGdVrheUc3L2D4pf1V6TSUNjb5EUdOzZ8UlH7nndIzs6RikY+DAaDsPKQS+8mgOCaqY8o2zQ0xskFJQJ3SbENGU/dUERaWadLwAbf0cawePF9sDxuyKUTPHqI/H1fGsaD4uHaBJHTqOkm9xlhUtcoaXkr3fpvz32yz//bWzw30v9hh+YrQiOZUSYWUblTjZcxO2Qdk8sHwYJpaWKfe19Lv/21d95Tr/UqLw+AxtqiQTUXTuDi012r1MkIIWoVZCKxXGWVZcYsWgpfIEipg9d34qn3VivRN0YjIsMKtGK5ednEWoykqVV5EuFc4pRq1xG85ZlK7D+cap9xTWKPUjNu6K2tUMs+hVoRMdK287N730DIp83Y95vsVLrUjwpW8LGcjynKOS27KKWoBsZvdviv2f4mwvRLHJ2d7HHDgFQ5xtyc4O96QtDbL7Bl6MFpYTGvOE11wMfIpTUWynIMELodSzKqbF86RMQTWaVVOl7URqUnTgIbKKUoSxOAGmUW7kmBp6OntDpduItqJE0EqolNAohbU1G3cJAnv1CnEa5yqsHnMSTUgOJWbto1OzqTigKN9VBbm/4KFGLDZvxyf4tDQpGgZSqW26ElInaCBtNzj/Dc9YuQuu3SUjlp3s0K4KV1CnBJYqXgLGHnwpMFW6X1mjIs8g5tV1czW0CINUep2+j052GLeUzIQ4U0EUxoyJqTLv+pvHRsHRZp4/gGl+oHTOj8Nsc1LvhxOVvk87O9x77W0SAvOpVml2YUs+gl0SMSlFRjJHc9obK7Z7AYr3QFFYEKPg+NXs4c4lobnCrNSZzc6hSpStRNy3MNg9B7nhtdxw4TZsbMvKalBCLZ4xsLLeYQzqJ3TujoN94x2gtoU4txd7EZcdKJD26ZxN0o9xaj2JfMkwQ4pSnUZLk9bzgjIWrZa0eL3jratN0kmopKVxa2pXo4Nja13LiMHISMeOzt0x2h2j61PUqdUqX/PEtIj6tA1lFZySKiT4pi9MFV58WrWT6DwWTGjVYt3gXx7hXqWKu1ggERtiLuK5/sV+smT8yKbJ1eXt+Tvh/39Ke5fJch+qnR3uSfu2A+O+CNkufFdGwTOO7CwqiVJ96fOCQJ+WS1FTietCvOVRvrBsR142H4xJLFU4XIf1nRIYsXbAqI6dukGhaF2dXitaCRVCaytwawauQKCT2yAE02BtjmD9sU95uoKe4LwxEp5qLujENojJrqilEJNL5bQ/SkhGGCJuQ6kq6SQoarTUVE5ToQNcorBoRkz41zG4Hcb56xDPpaySS8cVIt18rbOEooQX55Iimj+HotXRDOsudXZzS/TCOZY6CqLCyz3js8f47lKS7D4625zBEO0Uc6b8zC38/mHY2eH+3ixiX6cywHOu4zKmFjPSc5k+t1AfX5b2SkrMTCOYqJ617IynySUCPSpHUE168PvxFoBO31FLS+dWwIpawWihEuF51bI3FTt21LSs1DMETQ9Y7SUgU8+vWUILoKmuwjFPM/CxAKDSawazxRgf/UWMOa6jVJ0KGMppuDVDav4YHXkUpdnwjNo1XLoLNAotCuuM5yDImMp6tdSBoRBeCjN82RRt21NpMwU8wBgw83BPUKmceB6Zl447Xq9KrRjNLty7uPxshqKKqBqLmyTD4vhZmmGVOYY4cyvH6FJuYK4oFu2hWSF8SI73vlfY2b61nSqRLH9GWyrFDLZQGba8DJzuEFDsZ0aVOtpUnA7Pky8ziMK6kdH56jKLwznPUhABEWGlFWutWbsVjVtT0fpEW+wHpuacYHLnCHLi6NQx6tBxYqKcVfw8ihwldpPQwWm3aKl8RMsKTU3j1qzchjoUc8SEoAtkMCPHIjdL+8l/LF/nSUQfkmbJ2cbIfGm98p7MxsVcfewIG37Q7lu2HJsPbTN+/1hO7odj5wj3pL0L7vS4dXLCYj7Y5gkw5zEx55AiGZaTH0VUUZTsTiKbhNlOE0ansuY+yvSt0409oPWGSq0wtk84LkBdecWs3fg1UinWcsHWrKlHYV0pNkq4rAXrNC/HF3x9uOKvrOaNaFCgpU3Js8R+QCe+b4xUowOeJ5usG+nGG09XK3BeHXqulUmmCCHU+sIvo73Gba3W6exrWbNx17x0z2ikog00tyHMGDrp6GTPQJeKHeZYsrU9ld4E/LXGYTxVrJzepzGQe5lVapNeGh5PvkjRbaXXKBT74dVU+jHQz2Kbe89Djth7FuKJTjgmshyaUj1MiiRZUpYLiTfHUtVcGZUuJXnvi1Yfo5Nw1lI426Ns7jxPvbnLLrqnEhEPbaPg196zTG5Pc0wJi4pU80gqcUQTZWrh+CaOb8AyeGfkDJ21tKFVd51YV94JX3YXXl1MYMcrHCZhrMJUA7Z0sKciSIedMBfS8T9gkfolaB/dSsuKS1Zuw0pqGqXQIjgHnXOMWJ8wY/D/3JDghASHuNLRZTioxJ1j9VmEa2R2DyKeHPnDURhIRE9glcT1XZzaL9O7SgbNyTEzyw2cunoPO8O3gQhK2ln52Q/Tzg735Bv1bW/60jTrVHVMiEQoebhCFheZRrqCIKHDrKTIVB85xZy0iXqwvpMrhcrWlNlwf6IkSyiO2TGqzHcF71wO5gY0XMozaqNZaY2WcEYCaw3XjWK0Gy7Ghq295lfKsHUDIr45pZPslJSqUFQF1cvvazR7X4obhM6t8+XDCsVgs2jMHFIouax+f167VkvNlXzCC/sJaxoaqXhe11TKY5e9cfRmZM+BndzQ40t5e3vHaPaMZp+2r5VGUYVCii7DJajkQD0Wm1kKcehFmMPDCR4qSS8j8NFu4AGfrmRTx6yVYLmktxQFitvu0t9+NlVGtaeejblA+qmod26nMNtTv//w7OxwF9+w72LRSd5HnZkvr4vv51zZvJxLSGK5yeMIdBrdBf6sKFx5m4sHLsIN+cEO7cJDBFWW2Zb4qNcaGFDkBJd1lk4ObFlzaWsapTEOlECr4cIJrtVoEdQorLhkkD1OG6wa6c3dpMIMyFn4gipWliJr8c4WgiOF0Lq9iI6LhJVSHjuupKVSLbWsad2aNQ1rqWmVplL+JTEGPNricBJKed3o2Qmhfc5cqGaOcZfnMG9ead1IpVfpxRUTgT6aNcTOGvFnKuOdsTryCybzs4+TWnEMSQE9nWIivE10Oh/rp/jpP2wn+jZ2drjA+xkQZTZ3aZvH0WTJhYzTvPl0b3F6F5Mh5QMvZYQT8cHSocYHNCxSONvpVDhXTE1k/8qHK1DD/K9Nci4H2aKp+chuGK2jt6QItxJYaUGH/W26SwbVUcuakQ7njKdZFZStiM9GB9vWz9KxaKrASR1Db7EQRSaHm7vpWudbBFXiy2RX6tkRjLDRFatKqBUoEUaTiFeJ7jZyYLB7RnPweGoopLDiy4gVOQkYr0nWAy6vXyzW8EyJ0RwSUyJeB39tC2ZF6lZcbCZVt8UXimc5+J5meV+EWVPpdO/vGnJqZha/p1gmMnHi8o/RwP1wHfDZ4T7ayij0vkE1X2bOvS0zx/mz0w9AmLofqfcfZ6tFlV0bTh1XTqQlO4p6mNCUJNDDyql6Scr3XSFu2KgXAOzsACNcjIpahI12WA29FW4HQQk855LKVnytfsfItJ13dJZeICdEgqInEaXFogI9TqFSJOgry2w6LvBT/nL6XklLw5or+5xrNlzoikZ78R0tggMOo2O0jjpwcXt2XhWsuA6xKWVMkJVwRmQbOGtTkUmeOWT2hbF9UZBhfFmxvmS0XYBJhhzdkl+C4Q/ieJriu4UynMS29KHIRmq/TnqxWaZjcjpWjq2Eux5KcJ2aOX64Ue/Z4T7KptGne6sBU0YMS8kBYTqI59nsXHFW/j3fR45Ufa+sOdd22icrbPseahiFU43sgXLX5eeWMWnlAnSMKCscjGZTaRrlYzUtjlr51u8XukYZ4Q01/eSYdIgWFYYgRRjYBIPdp2X8NTlBW3M5qg0XjSoUPQBoalq35pI1G1WzqrzKmRIPgVgHvbUMznqdX8iJMjsQq6kqvZ4UXERM1+9SZ3pd4ZeckLoAg49StS6q5ETRyCVWLLCfONt4X/z2I5sgrgeRmeKPz8NQfrsVuEMYQbH82S/v7ThnEA73hJXBw0OVZWcr7exwH2WxbHG50eD0jX3MCMjfuYXvy+XmljFUnCFqHrgCy4vrTiqwQoFCpFxB1mGNegk5k573ETVb/Xe++0JyvDH5U3JdBYw1ns+qGnbmawa1o5I/o3ZXNEqoBA5WMA4Oxk/ZP1opNlVNZyqa7o94Yw/8O733XRPMXcAic1LOC4JrTKwiw3NoK1kxutAi3PmoN9LJ0pVxI1hfMNDoC7Rac+VesnEbLlXDSntGwnWjuKhgO0JnHJUSOuP4TH7HrXzNYXyVEnPRWcbuC8O4C9cvR9Kxy0XCm1WTEmPRan2BloqN/hgRhQkyj30oG46yjwCNvkolvSXYMb/3cdwkpxnuZRpPwUn774+7A09tKZFcOtPH5CuWxG/m2/1wIt6zw/3W9hCl5b5Be6o8Mn62ZDGhlMWkIyRQVmylozuKZu6xSFWa0Il9BKulOVp8XmjguamGQfcMGKzzyafeCtaBcdAoTxdTCJWC3lbY3iewjAyI7EO5QY74BI2a0MS8JKPMI2NXADmBSlVGjioUOKzdijU1K61otFc4W2lYae9wrfNY9+AMO3XD4PaYSducvL2ogwuZuTG/PvE7300iY96eKdFSSRvazfuXR+cOWAoYIZy7osawrAhWJm1l7hTTy1WO1y1nORNe8fvKa5yttLPDfSubv6nvw7DmjricfnmvVj6g7pSDDc500u58vkgsjgg19f5nSMyoBiWx08E467rguwMoabLzkjwV1apJ2XiPDeZpfOwcq3UTCgx2XgBGKnZyi6D54tCyNzW10lxUcF07NtrRKsdvD5rtKDjnHddPDn+NO7nhpv4ytK0ZWFXPk0MzAasFT5My+NY2qTuEqElkLCisDBg30uprGtnwjB+xchteqguuqoqfbDSV8iwKLf6udMZxN1h+Y15zp27o3B0OQ63Wnl1gio4S5BLrOPVPVLRAbcuFBQZrSTMPpXL7ntF1GBlYc00vezp3m1r2xO0NZpsSZ1PHXyW+dbq3s+q8iaRnHFPhXs+Xe/jFHMe6njj2rBA2fzbKdUqLOZEPS1ns7HDf2WT2+32DRs2WO8a47mUnOHv0EE0foMcN2KVuCceloCQHFqPI5c4LU1J/nCobNzJwoJc9BzvSGk1vNRt8ZBv/bSp/HTojDFZxyQrrLAfx8oMRD55Hsaljw1HlVnx55Qy/ktpHteKVtzb2kkvWrJSm1UITHK3go2/rYLCOwWbdhFIDYX6dYjXY/Pv5NQoHNIu2KzQ1Cp04yNN183mmZJwr4YOMuz5U7DFXjktayjFZ+ohikcfZ20ID72u/T8fODvedbCl6nVNpciQwHVhlVJuX9w9F6Nbq/LTQhf8mkUSsGgpJMhE1iVrnlU1lRVLZ9+r0mWVtAuuGwHdtMkeW3OoFwFrPImira3pzx2C2bNWXjKrjgmsqo+hNhXWe3+qdGrysLc8qUKJptHAzrNCD4uB2fmpth8QyqMR3cdDiK8NaufRRYaj68st4dkJv7xKUsNLXrOVFYkH8WD3jutE0SlhXHkYwDnoLvYHOOvajo7c2kKi8gzeRKRCUwTxrIjM3ouCOKJV6rJX3Y87PrfWalX6RoIT4UhnoCnGcCqsauuENoz1krD3BQypFy9b1uFBEMp0FKXIeoHCw5dibiCDNx/DSDK6MasuIeD77esjxfjhRbWlnh3vSHvO2LlkGS4kAmA7CuaOzs89s4khK+C873SWbJ0x8ya6aNV+MUUxU/JqeQYHbRRiCkpcbtlxUMQlF25vQzXcIbbtj8mdwOwYZGJzBBPwWoBJHoxw6nNDz2l8DLVCLYmU37GXDoPYMdofDooNDG90Bh0pOKR/bgJ1FoJW0qTCjYYOmolGKlVZc1t7ZNso7fyv+ThgHnbV0zvh27+wY3C5VtSlVU8GEYSASuL8UjrWESAsIId6TSm28iE9ytiZF9aPrAvNjGv2VxSjzCY3fb3S2FBFrdLD3vWBl5jjvd7Kn7b5o9V1Kd3+Ymgpnh/tWtnTzl/GpXKZbFkPEAT6l0iRZxUIyL2eSczQzT9o423tnXPY0kxkEMJtK5lY0gcEQpAEBrOsnztR/Nk6ccvx+WqDQYWyX2AqRKXCo7uicL4Kwzl+DRjkuK4sOQgJrbalVRR06RVzYFZ27xsiAUwbjBnRwnCOHENHuEtwRS17Hgp2gQ3vzKgjZrNwFK7dh3SjWlfBR653tSjsORtLLwAU62MENHNSWzt3Rm23abq3WKH01KQRJkAMKHTBkI53XyA0vOaVqv65U1LJJka3XI9PB2Vo6e5ebUs6m+akYJejo+qiWNGNJvdKik430sIDtL9m064gja+WW43MWECz+7q/EaTsFZd1HlVzOVzx1Ozvc92LHteF50JZMhDgNmw6mslFk/D4mf2I442yfnGfs0BAflhQpR4ih4G2m9jQhsp1/lzrpuiygXSpSKakSQ6GUNizxxbS9gkUA0LNnKztuhitqJVxUmYr0STuw1oYxOOI/umz4phPUbo0d/fUbpQN22HAtdRCcqaT1jswr1np2AyExaG3g8erUk+wj95IL1bDSQi3e2SrxlXC99VHufnTc9JZXds+d3DFwCNF1fkSMGxhDT7f40imhFiszQZuUrDQoqVipZzx3P/J4bmpFA1YCr1daRkgFD2WRhSpnJqJQwanm+6mKpGng/s6LYAq8No2vYjTkF/2MbZHG9n2c8vkzcIoq6WbfnXK6S7PGp29nh3vSHnuzy7fxEoY710mIdk/RQfg+Y7ABVigwuLSNNFYtKnS8nUINKi03x3sJFU/OjqEH2MpvqeD5iuTeWvFvf1ahEqzIpCcqV/jZux0HueNufEmrFDeDZqWFVilabbioB4xTOIRP2xpB2I+Kg23p7CV38tqX/UZmREgy+ejQR4giigofyRoGRjmkl0JNy8ptuNYtzxtfulspqJULxQ3CYAkcYcd2NNzJHVt5XVCyMuvB2H1qogkZgokt2lWo9EtykWQdh0p86/hn7jkGL2g+MmLEoBmxATqxIVE3VwZToRrN2D7h6qM94GwUMi9pgVnECOzkZX00/oIU48RmRTO+Km0uthS3UQYWZVL4vmfiIX7v0jI/DPvAHe59OO27krFL7GkJWsgDudQizT8jt3UAl1uhZCJ7qCKSGOnaVPZb6p9KaIE+F1FJ+ggxAYZvTa7J2ORgtkjBGR3GbSpVhWlGPjqXWl9M9FpFFKPq2MsNX8prZHjB9bBipX3Dydd9g3XCSvsHtlGO5w0oUdSqod5fo8wvuJNb3siXiXI1Z0woNKZgNazVCxSaipaN86yE503FR63iovbR7UY7BivcjZ53+6Z33A6GzmUJxoN9A3gYwcZoNhQbaGknM4mq0NU1tsO4nkqtqVSbynVXcs3aXVCjsTgG2aNdReUqj9cKbN3XQathnwRywBerpGjXeZzesiSes6yJccRQYJpcjVcy31MDFB2jwW8vLRcTwaWzVWm5/LOMjOFxOgs/TOw22kNh1g/Y4iD5PZchHkUWReS5QCXy5h3vJPKYLHvqQYkRyQL2dQ/1J0kihn+JthQaM5a9xyb7ieujUmLLH31OqBk3sJMbDm5gNzoOBjoj7I2ityqkBz03d6Mdz2vHdSNc1ppLVmzcpU8uxeRTMd216cWVr4GmpqKldWta19JKxboS2pAkaxQpYRcZEwdjGa3DOIsV6xHV4GBj0cXcYnFGxm7V5J8OswetWmq1ocbjtn6uYpP7Kq9ZvF7WjoHtkPUTSg3cY5nG/DKV8EJfpoqp4l9wuvc4tjizOtZaOLbjkvMfpsP8tvYBR7gx6rzve1h2yPNp0jzDG3+PoiFQ0sOylqyPIko8LVeFzTrvTrC36eCX1N48shKmkQ6Q+pLF8t3U9ryInKigrOCydvRcVl1PaGCxKCLuQ6kaE7Vt8VGhiPblrWI56Jav5Ws4QKuvuKiE14PCuJqLamRTjfytqy29VWzHKlRgKYxrqXrFnbumkz2GIVG1oq3cBZqandwk3mxFzQv3kpd6zUdtxScr4bp2HruFBCPEu2csdM6wZ2AIEEalPBXNl9v6e9jqK9BX4ZoZKtVi3JASa7H6TFdtYiKs5JqWDc/sS0ZGfqe+oHYttasxMjLgeC2fM7g9B3PDaPaF6E64j2E24mQK2zjGpCwGRe+zkjGRnO8Mz2WehC364Z0MQuaFOjmKPXbcZaT6w0t+vat9wA73fdl9EXKJoy5/P9G6fWcC+swBp+li3v/SdDzZPGKdcUcjHhl/98ss1+F7vYfgFCXDD4N0bN2e3XjBwWgGJxgnOOdFbS6bgdEqVtpwM1Rsm4qvK8XBaFajd36dxLMpji0el3f1OCyNW3MpLReV5rL2GglXlaeiKRyjk6PMvXMOU7yAowg4+EgYyS9KIuaK9pFukSz0kW3tS5FRtGxo3JoYn1fOC/M4yV0kRtclTYhSDxiXoaAU0c5w/Hyfjj+bNpWcn2/JPAkXIyZhl8p951VqaZ/vMmZ/2LDBfXZ2uIt2CvQvv4f7HOl0vanjnfYls0VUm6Pqaavz8P0MO/VOIw9652Knh7JLbJXYCH6Z0We5i5r+EvsT0SlJFvFZFehO5Xb9UU4bOMrsoW/kMvFNLZaduuOr7golLZ+uFLXy51try6cvb1OWvfnSclltuBtbjNNszTU7t+INrxOMoEI23QWnVbs2AA41L9w1P7toedF6CtjPViPPakOtLNYJt6OGQRNdTaX8/RwZvVOkTbc4Jgeds4mWFgshTGiBs6qeh2V1Ot+RA5aBK/ucmpqBgZaGj7ni1nW8kdfs5JaBA4PdBX3dIW2HUJ4dy3lLBskEaZL8MoxQQ2ZslAnSaeSbaWYZCorLRapZ6oNXJuEerEwrnfAphzrHe909y/6w7OxwF+1tp0Hzt/5j3vxl8mzaameR+3iE4c7YCunz3DVhqeJ33ndrbl7XdfoiidVUnuUQKqFCa5vUoUDK5Y1/fYhNClhOLApFx8hurDkYzcEIg1OYsLLSDlGO61XHYdQ8qxu2jeJ1X4GBvVthxBwVPnhZlwbtKloaNqqmUV74/LJy1MqhxHOBjROGoF7mq8wcd4NhS89Bdhgy77iUpBRRiNMI/jyskM4/iesEBoVnSGw8VSxc50t8lBsb5URucYR0olhNlFe0dsp/hvCCCzS+LDBfYObJic6EhiYv7ln0ejRBK8ZSygmMJyLc2fbieunzpd/L5SI7YSlB/cNUEDs73JP22Jv9mLd5jmYXHWUin0cK0j2qTkfUMDtzohnjm7AJQuSjQ8bbBqrQHGZwznpRmOBArB2Ljrg1o/Oar7FFuU+8KHBDij5NKKhQRVQWCxH2dGxNzc1Qs6kUnVEMVmGdUClLtbJcX+1RyvLT7QbjKm6GiqoXdsOKzvUYmTpcgMvAAFipiquqYqWFq9rxoja0IZLurWJ0sDOKQ6CE3Q2W35lbXqtXbHmVrlU04zxuHEuHEbBSI+llqZKzLfnBDWs0NUa83sLzqmW0jld275kOTmMlQwrGdoxmn6lmE2peVeC5XtHNhqhXyA44RbS6mY6ZJbjBThOsGbvNTU4zLmshsGPmNtdr9uM8Ctto7kvKZSvpZT9sOzvcR9kpnHbJ2c7f6ksJtaktJc+OLTrn4ve0fgkNzKMe/1PLasqXlayydQrbdZikhTsl9GcB8FRNJTW64KeWFivCenZ00nFwLV8dWhyKn640ra7Z7Rq06lh9YuAGuqHiWT1iHLzuK1qtse6CW9PwBQO1q6nIHNnnesVKa9ZauGoUL1qP2660RcTjthYYrOJmFLajsDc+gRZNU2PIRRQivrTCOUvvdmlftdokRoFzNhVDGDfQqkuEmo279iwJahTC1oy+zBlDz8Ag3pFrao8XS2YWaNUgVmEL/eWSk+uX2UzugT+4wrEWMFK+JyeoYoRhFZNmcTwhDzpMdzT2p047f14moZei3cmR/GDt7HAfZUuY0yngv/y7hAnK7Xib1rFHcfGSxzjfv/99zoPNSZv5A5anmX5/x8m1sqtDWa47Xy5CCUr5yNYVoicxOx8jPScx2VQyK2Bwewbp6NzAq87r5W6vFBut2XUNm02PvvDH2BvNs2ZAi+Oj1jd3HK1GdfD1oGlpaKmwOCoU13XFphIua+Gigue140JbWmUZQ3IuavLeDMJu9G10hihU42K3XE8Li9dKUzOwZ3A7tNQo6lQybN1diEYHcAR4xbMTNm7DmppGNA64c10qeChhESURI480NH8tUUUUSi5GsWHmELm5x0ms8r6V46EPDjX3ugtbJpWOMzIVsynhrlN2CjorP58zfpaKIeb2w4MT4OxwH2lL050Sf7ovwZaj1jz9ChqqM1wtZZRL2CDhbFNHG/fpPyuj2diLLGK5Gc+LCZdYqZTOJPI7XexO0FKasf2kh1lpOdr1DlSJrwgDqFWdkk1xmZ49O6m4NStkaPj1vqazwlW9YlWPXN/esPnI8teuv6HfarpDRatf8k3f8K91w7qqONy8YHQWg+NCaWpRXNVeAezTFVxUjo+bkcvKcFWNdNbDFl/sG25HYTvAdnS86gyvx44b9TrNCGLVmj/e8FkodIhwQVYq84LrWrXUsmEtz9K6FYpGNK3y10e7FYOz3LkOi5ehjG3XD+NrRnPIsIGbFqkACbtVMIGCIpbrnWvAg12fxsGk/FemTnY+vpRaIVJhzI5YgRZFlPz1eMgJzgOTNEo4ndeYO+G3cbRPD+c9O9xvZXMHtKx3e5rXOMd0KbiQ08+W8Li5pSk/gbPpMnZb/oxOt1xncryFVsIpsyF5VHb2LSGROF22YhCXPzcMDNLTuZGD0dyNNZvK9z87DBWmg2oD9QuoDyOr/chPdlvqnePLrqK3isvaJ9wOxoQyYS9Kc1HBy8ay0ZZn9chGG1aVwY5gnHCwwt4IB+PojJdh7BgZ6FDhv6mpQP+qOcK6MR4KCpKRlbRs3BWj+GSXRlGJb9+TtuUyCdAGSMKEtjlRCnPJfIk3gblw8pbkY0sMlrczCboZVg5ek6Ggi/kNLwUeS8c8d673Rcjl8k/Leb6LnR3uSXvM2/P0tKmMCibJB+BeOpmoqbOdENgjx9b3KZs3hiwFU9K6aRp6LMuYhFecSsUMqWEkU+J9jI5jpOusL3hoi2hYiecLqKBlUOETTSOdd0q01PjE0p1sMc5wO7RcVv5qGScMdxq9Nsha0/ztZ0it+dHhK9ovvYPe6JaD0dz0ittBs6l8ZPsnF45P2pG/+8nXNJVBBWYCwG9fX3Ewmr0RbgZ41Vk6m6u9Qpc0BC/9aFO065tNampqV3OQHYP4pF3NhrryTISGDa1bs3IbLtyKldTUyjvbRntGxG4Y6QOGCwGqcLtUNKHm9KyQGKv0CiU13fAqwQk4y+h2iFRU+jJDC6mJqKeSGbNLFK/EWkmRa4n7+perMTsMuzwOZwU00w4P82djib2wBI+VJrPv35af+/Qc9Nnh3mtLCkhwfKPf/sY/JiFxaop2ivKVD2fOZCj3mxNlc9rR6WNViRdadi2AXMF2FCUXL5UoR1jTotCoQDuzOA6jY2/EMwiMwhpBtKAuK+TlFaxb1j/5Emf3vNhu2I2a61pjQxJso33Z7lXlo9rrqwN1bRHlsEawRtDKzxis8/86a+mtxbjocDW18y8bf4tnervhv1iFF/UbmkD9ivzfmpqV1Ky1TpGsClGph0B8sYPfZ05ulrONeG3nxSf33hsy/JHHhnd4rpjNHK9cYrn+TO8bO4+zpXXvq9r8sOwDdbiaY1z2MVhSHEzT3mTHXR3UEf7lly5pWtPsrqf7lDSw6uS0PvaxispfshAF++9yQiSyDdLDjEacOnKYkZmgiv1HmUZRmirKPIYIuTd3vo1NjHQlO1iD17J97n7sMV5y0cLaraip+KobcNT87rLmqq5R2lG9qFA/fY7745/C1SW1Eq5//RV//H9/RfOF4W7UvK4VrwdFoxytgh+vej5eHdg8G9CtQ7Uw3gndraZWhqoQPd8bQ+dGhsCK3bhLWtdSU7F3XSi57bHYUFI8MsiBFZe0bh3YCzUbd+W3J1sUyle36YqLOt9n4xzGOTpGOgYOsmOUkVITQqsWawff2SFqMagWrdqkNeypYbnoIDpa4/pUURZfpFO1sNzGXQQclc88TJx8poiVLAVXdvsl5x6yFd2B49+pRPhU3iON4nu+++HaB+pwlyK7Uzd8PjAeXxAxV/o6bTPdhHLvacqXa+YhdhFoJg/MMR4bBMdlirX67WZKVdpXoXUbmQeWYXJp5p0IcjtzQ+zG4MSmbD+QKrIi2T9GdxbvjAbr+LrXfNw3dPuKi/2A6v002VUVsl4hLy7Z/PRrPh1v+dFuTa1qL1wuvrDhsh5Z1yOq9s5WbRSq81jmYDWdUUn7VsIVtDgExcZtiFI0Ppb1J2zFYhjTeVgso2RBGfDQwLV7wcZtWGnNqhJWWlLV3O3gz2/Eeuw6OHBD1Ls1k5eqb+UzTl6MKQp2TPUU0iwjFzDkopf8sr3PckScqWBH8oyJjjYtZlgc124mujTdW1xovtJsmR+u8/1AHe5jb+hS1Zc78R0Ly8W3/3wbx/Sw+fRuir1NH8jITtDK69cO402x3VxpBqBUMxGaGe0BVdz28oHUukm4bYxaRwvG5W4KcwFzL/YdH/whOAWLFd96x4llJzesuGRj/fEqVHBq3pEdjOVXW8V1VXG7bbl+01HdHpB+AGtxFxvEWZq/85Lnl6/5G/vXfLRb8eqQ8eMXqwMXqx59IehLhVxUqEOP3Fp6o7gdNdvRd+WtlcJYh3GGNQ1XqvVtdZzBiWWgw4aXRqRv1ax8kQJdghUAalfzsVyxqjSbSgXtBkldgF91js4aOjr2suPOfY0Nvdi80Pg4gRKcs6l7hipoXUmbIr1sp49uLAGe39eyF9rxbClUrolNhTB+3dgr7zBxoCVFLIooOTdQOub8ZMXxXxY1zGeIcbnpMT2cQHu6WgwfqMOd26m3alniq4plT5mvsvF2XyQco9RMExMH02Z/0UmWiv7ZKVs3Bp55nkJOxMHDlD8up4OA9VzvoOwsa0yPofdC2hKnldOIuHQCubTVR8M+cx8iW0fanw6iMn0QB1coahdExEWoxGvkHix8ub3g+e/2VC921L/+AhkG0CFK/8UnVM8v+PlPX/Oz2y3mzSvsncV2YAdBtEMqwe4s9lXH7ouKb15t+PKw4vWgff8yB4O1OOeoqRDE6+AWjkWh0C5AIqGirg4JQItJCbaNu6SlZqO9uPlFpWhUdrbGghR0AeV8r7U+dHOIL6zRHtLsRKkqNe00tj+amicKX+xsIXW6JyY1tazS/SewVaZjjwA95eRpTsL6MehF1YOTd1kEfz6OsxOej/cIsy0VOdwXxT5mBvn0HG20s8N9cLof39QP0VzyIDhNFp8NpoTJ5UaRvqDAO10/LS1xtIjHedzOD90MG0w2HShhOItlpGLN3ASdynYB+uE2UJqa5HDL5M1c9zZ+FluRe80EExwToQFOTRWghEF8zzTlFHWowFIIIv7VczDCl4eWj99sWH9+Q/XZK8QY+PQFrFe4P/oprl3Beg3dgWq/R774Ct7c4f7iC+x2wN4OmFubnO3nt5d81Ve8GRTGBUfojzQVTgzO4MJn4GEC7bw+QseeOujrRgZDhEfWNKykDjCCog0JPC2+dc8YhoTyni2tJ+jUbdgXS3jMM2LplVozhu4SU1WvXH3m7Igo5f/hJXfGAl4qX7p+iC87sin3279QjT2AG4PsZyFik45DJsk1P95LUycc8X1shXfh4T49OztcYLmWe+ktfN9bOS5fTqsUSrVFBNECCmcPeb8x6nQW5waUahGpsLafCJTEaMGlRFiGCSQpdo0Z28WiQoGDUtVEY1WpalLcEKOrVfPC7ynQj2q1JoqqxCaI0yuU1bEqaRldhw1OWVBseEbtGlZuwyADnexp3Zra1VyypkJxsCN2dNwNwrZW3I2Kr3drVr8daX+2p3qxwfyD/2jxqjuAP/nrk8/q/+KfoP/1Z2w/G9gPNd/0DbuE3zpG61gpjQvc4NFZRhdxW0vrWmL7G4WwDnq7PiLXOLHUrqGm4kp5CciPV9p3HVb59d1bX8nm77IwyMBethzcDYP1ZcKx6aRqcoWhdYPvuCGKWl9kcZtQRqyJOgn+R9llQ6kp9U/JrPABXxLsKxrHIrolUMqKl3YJZYXGlZF6diz/CcfjtHS2pyCA+Mz8sJ1saWeHuxiJyuznEva0RO0qnXRkKJTN/XzbGxPKLEst01hW6bPRDZawjCq5kGHLJzQLjmhbLqh6qYrReNEZL7dYp/Y5sXUNZOHwwe5zx4fwsGppJxqx/gx10YG2BvGEfh24uE1wrjW155Vi09+tVGgRhnAOxvmEVm+FvanYdTW4A6xm/NQHzPzH/wDN/xP9T/8CgM7mrrzgaVqtUqkqbrASgjuHcrE+L8TozntQ7SoqNBav/rWmoZGKdaVTwYUWv71YPmwDOwG8+4kykjGyVcRyaA8JRKhgtL4dvIcXah98isWZAA/MyrrLMVBWqsUxYPEv4cStVk2IYsORhaq0slAiOtTpvij2fwwhHD8LJRRXOt4Px7ku2QfqcEuc1c4+gzwt00QyeblMFJvxb36YDL44sPUmPBC9d6LBwcWa9myRFZAbEBp38HhaVImawQWR4pWaK4Z1la5SQiU1UqwuqKRh1GsERauvfRRl94gKHFLVErUUlNQ8r38BEKbQPvlVyzq1iQGPR0as018BS4VOzknjm4BHnPqChkau0UrQImwqRSVeaKYWuKihVXA7KkYr1NpS/cO/iflH//Dtbm3dYP7xP+L5rkP957/hr24v+bgxbLRCi+IQ2vsIsNJwN8A3ndc7sA72o2VwlsG2ab7ikhP2zuIX6w2XtbCphJWGF43jqnK8bEZe95rb0aufKRG+7noOdqSWhhWXjCrDBLWsg3auLxXem1dUqkVqlfBxE+5trS/8dQ4vwFjlJihGu0+NLCH0pMPf+2i+ZU+EMWL0PNCHppRCLnLxy1RFp+e+YMtYoNBwKPIK+VlQIenWMWE1pE7UMT+QE6fHtsSBv4/98DQc+QfqcKPlt+8iX/ZEOa1Q4FeTgUgYdSpllSMJPf4N4I6m5iqtO6nymvFnJ5zYwgnHKiXrhgkuKygqaVDiK8IETSMbTBBK8S3QvSCLwssFaqm5cM+B4HDFa9k2bk0VCgTy1ZOAH/oTq4OTrWRa1mpx1KJotKDEX+vrxtOnrmsfHbbahRbmjrU2rFc99j/8Hz7uNi6Y/Z/8R1z8+X/O8z/vsTQIMNQwVL68V4BWObR4hdsYldbK6y6M1rse53yH31yVJnzUKi5ruKj8Np7XlovK8LweGKzQWaFWviW7FqEWzcp5DN3IkMaBpg6JtJqBzjtjDNq1SQgodoIopS6dMyHazO2GDNP27YniR8bgsZlLW34XubsR80883oTRlgL1YYwWLJrSykKLjO3OZpEFn/gpMw7excS5x1Rng8gPyTcfC8ocO1e7+BZP3yV6TFBuKgoNjvamVkmH1jlLpX0788HchWMosFzIhQriW5RHPHV0PcO4LXiyHh6oQ9fY0fkS2rroIhtFsV/wEwBey+9C1n1FTUvtGgbpsRiu7AsaalY0KLyj8A4jHhestEIJ9MU8vdHCWitqBVpJatTYFpOGEtscLPzta8vzeqRV3rHtRk2tvA7C3/nxl3zyP+pQ/7v/zQP38WF787/+P/HbL655dVilz7qgv2uc8GaoeD0o3zLdeYjBuuPH3zofEbfK8Ulr2GjLR21Hqy1tNaJwaOX41e0lnx0aPttr7kb4Ym/prWM3+kaVQ8CLnXPsGTAhWeaw3Mkto3hth9ataVlxJ2/o2ZO4wElkXCUo59Z+QW/vaNQlShQHcxM0j4fEpx7tHmP65Axjos7YPs2KYuPKOPZi1+AIPxl7mETAkyo3F1vGHyhLzYFp5wkgae8Wz5VzHW/vdL9fkW1JyztlPyQv+kiT4rfpTZ++rSuQqbJ+trz8ovBLcKAABJqOVm2IMDK+pqRJjjN+liPcOrW3iVhfLRpVVxM6loiiklU44pYoHxh5ohUtFTUbewngq5ycf1i9XmvDEAbKhnYSoTbKO1ctglb+eq1ClLoKzlSLd7grDZVApTw0UCmoxCVcc26ftAPP6oGVNmmZ0OmGy5cd1U83jyIIPWSXf2z4qNsiArU2SDgmY4XDWLE5tChpGawwOKEPDjeKlo8udocQnhWR7EobPr7Y+bJhcbggAbnW3hlfVr4j8UWt0KOPko0CbQqw1ECHJLbE6AzWWSqeewKajKzcJTUtRnwBxiCeE60CaKOpaZW/t414jdxKfPJyVJFZUiH6IlW0uUj7clkbw8+MSLoacTxOf2Z4rdRJ9l+WUevsWZp/F+G4xeTb29j3x9k+1j5AhzuNaiU4ven3eYoVhUCORLXLqNflljZRnb9MSmnVUmnvFON0D6CprlCqogoRqcMmKlB8EBp96aEB5Z1pI5sUWUSnGkWzV1x6Z+pyounaXaWpPsClzZHeSlWsgnygxU+no9UK1loFbNPRaD9NbpR3jNG5rgMNqlYuTKF9u/NaTRXPIkZ6VRsuK8On6x3rZmSz7mhaw/rjkXEr7N40XPypRv7sZ293U09Y8w9+wUt+RftXhvbZiFpB9Uzjekv/jePNb9for5+n5W+HGgd8ujrgEHajDqpmmj+62PNys6euDHU9cv0j7/z6W83YK/qu4rIZeG40xvmov7eKnRKUKAbrkqiNFpAeKqO4rCoPW4yOWjTP64ZXfceXvObSXdKgk/DNTu7SsUb2BAoatUmRslE+OSdOJ9ioYoWSip35mtF1DGz9OCvggjjmrA2lwgmb9RoTumBBJBjLjOHvatFtpi4kkWNe5CyUNBi7S1rLH4J9kA4XKCCEY4EPP0Uqq3Yq5nJ3c0X92KU1J7JyWawNEnyRYxmbMs4HmpC1aGP0GtkBMcESk1a5JTapnUuMYhUqRESOgTGQdHxHhpVorHOBs+lr/WslaMmk/RjVNtpjcsb5z7XAde2d8WD936sCe/XnDReVoVY+6vPXz6PiWhwrbVhpw/WmY7Pp2Xw8oFrQ1xpVG6wZUC8ucJ+8eNc7O7WffkL1yRe0bw7oC0E1oK5qXG/Q+4Gm9W3ajfPMhWeNvyfr2t//Slk21nBRKZ6v91xeHlhdjujWUj2TEJwZZOcYBp2uV4yYwV+ndSU0Tmisvxax5Ng7YH+dn2l/7wdrERE2bkOT2tV4iw05B+kmspcWQxu0fCNjxO87NvE0WJeFc7Rqcc5MIApYqkbLjJfIaBFRkwRcyYKIAUyOfjNdLWp/pO4TxH3rxf3+EO3DdbhATHbl7rfekmxhWeqImji5yGeNzf5qfYGIYhi3YR2VHohYrlnpNVq1NMpnj3u2k6PxuFxuQy4oGtmgpGLDM8BTlAbpityvonZe9vAiFDd4Mr9ikJ5BekanWIdCg0YpL6hinacu4WjE468r7aGDWoVutuLJ+5KcMjyrHY1y3I7KO1DlaJX/zAQHc1kZGmU5mFD+qiy1sqy0QYuj1paLy47V85Hm5zXSBDxCD7RmRD66xH780Xu5w/anP0Y++guqLzvUSiG1Qq4a6Eb07Ui1MrSVYbQOY4W1HtHiWAfHu3KjhwOs4uqiY3010H4MaqORWuGsowoVdnLrG1UO4d9oM1RyUUXqW5GcrT0NbrAecrjWmiE0tdQIl9L6pF3gSYAvJR5C0jMKAkUORUxoivhZjogvtfDFMoYxlGiLaGpZ+xZBJo/zzIA4JT4ToTCPCU+SeeRIN87ynN1BEZiIlDPFovX7IjXzh2kfpsOdQAix08GxvkDuA2UnmKunbvWM4y6vExIPJlTlKOupWeU0LE7V8t8hqaGq5JSd1GipGOw+LaekZiywO+t8rdRKroHchvxL9SWx/HTjrlm7FWN4KHd0aKcZxqCV4NPcIevsnU0vDixsnQtsApIYSx0SYW8GT+syjsRCOFj/GmqCItervqYKEa8OdaWDVWzHipW2NMYw9Jq6M/SfDYgaEO3VvbavGpq/eoX6+a8wP/rJt7rNAHJ7B7sB2zlsZxBtqaoDjA7bgx08FGSdYJ2wHyqsE+76Gi2OpoovCUNVWVTjsAO4O4taOezBsf2NZnvX8tXthl/tNny2r9kaoQs902Ktom9aGfDcyNF1cDuY9LdzpMSaKcZKH0SBaucTmgOdn+ZT43XPhoTPR0ZM7MfmC1L8C2R0XS5kKdgvk2sWo9bQZ60fb8M4DIk1xsQlN9FpTvQ+/Kjzx3HaxWQe8VJ0+/1KiL0v+zAd7oItdT6Inyf2QazuCthuKY8YNQtcESUgHv/KGgMeVzPSZfghjKko8i2isdgEN4zSodyAlSElN+L2ar0uIAbDgTuv1oWlpmXDBvCCLL0bUHgcUKNpqXyNj4MxUJ+iDeGPWvkEWaN8IkmAfQh8IsSgRNK6mwoaBVuj0Nb5oDXACaPzTnd0hpUS+r6i3htM70A5lIZhp9luW158fYf63av3c2O3O9x+xGwDvKEdajXiRofdO4ZeM1jNaD3/tw8MBoegQ3lzUxka8S3cRcD14JTDGYfZw91ty+2+5XXf8rrXvB7El/ZaH9U659/xg4XdGHi9RZLwYAxea00Sg6Hk/Tp8HzSLpU1sVxP+2VRQYWTALTQitQzYELUmWlmw1PuuSOaWn0Pg4QbnmarT4ixwwScmehmZeRNtycFPlch+eE62tA/T4SYeYX47O3pKgty8oaNzltHchYh0NeU3QqbGB4c7pvLdbDGzG220ewSV8F0gtcr2xQvTwWrdiHE58t3bV/Ryx45XxOr+KBSzldfcytdEoZWVXKbj1K4KlLCaymk612FxbMYNFYoRS0vFM92yM46vDiQM+6rWKbLVAhe1SvDD3eiTajFdttbeybzpPUXsosJLKirH7fgJV5Xhqh6olaWpDKNRHEzF+l/1vBx/jf5fvtVdXbZ/+ee8/mcjf/X5xzTaUGtD+6uRftR8td2wHSveDBU3g2ZvFZfaUqv4wnF0RvOCnqtVR3eoGL9S1LXH/b/8+pJurNgONbdDxZddw+86zdcdbIMsYxxSWmA/Or7sD6xUxUZrvhk6dnSMGBSKZ7KhdyOv5RbtNBUVnXQMeI1eQbGF9FI1jBzkjr17xcHcsONronLbvAO0CYI58bv4tykoW3HMlYUUnq/dhITuIXHDvevI1MlSBMc75LHgqsfEdJk3mTfAnNsP0/F+kA63FFnOHxZ/zyCH8nffpqTEn/L6SVjEFe2o4yZjRCukhENUjDKhzbaSOtXMQ46cHb63ViK1h4fJ2A4rIxJI72lfonyjQzekpElMokQFrBWXDHRoanrJPM8YZRnXUluNcZYBEzS+oqMVn2wTXzAgQKWExvif8VEZKu9wv+ksm8rzXj00ITSqZjtqtkbTKstaG0zkxt6tWX0+cPXgnXzY7G+3vHmz5lXX0qgADwR8+euuZW8UeyO8GvxPU6ugsQu1dTShEMI5YRg0xiisMVij+Ga/9r3YjOJmqPimV9wMHjbYDjbNFJSAiO/B1jGmIbWjYys7vHCNpnU1AyMH2SbK10HuMDG5FKb5I4OHBsRgqRmDSLm/vwH/L+hdMSdRJmmdmDSe8sUiwVsTnYW4TohoY3BRfj6FBmYwg1iW2P6P14v+4dgHXPgwrwUvWfqBFjaRx8vVM2WRxLy23ZtN5b9SYLiJLB54iaWoSEw2uCJqiOv4n1kHVUmd6DvgdWzL/ZfyfE11hRb/UEarVEurrn3U4gZW6jrpvebr4SUJo1j2c/cjNm6DwTAyslU31K7hyj3D9y0zXLKmjccLbHSVps2RyxspURe1L+nVynN6P23zMPzFpufnmx3/4X92wP2v/qe45+/GWND/5L/k6//DX/D/+dWn/Gq3wjjPq3X4n731SauDgde9j0Ab7aP0j1rYaF9F9kk78LPLLVocWll6o+lHzS+3F9yNvvPEq1747c7xTTdyM/bsyW3PgdTLTCFsZcdB7tLsY88NUfbRBGfqp/9BtAabi1tsh3UDo+tTQcx++IbRHHyPNdFZNyN0ejC2S5VnpdhNNOdGjD2g9YZKrYJc5DKJv4TYkmZuOLM0dsLzY21fRLrFjBJY7goxL+Wdf/b9tnPhw1tYKbzh9T9HlsQ4HAqSlkLBMzxy4N4xl044rz/dJgSFLuA4ch7TdgQF1meEl7Cwkqrjj0GHRFxIyOG5loOdrje4PTZ0NohOQqEx4qedFsNObhIdyRFFuk3ClL2qmFfbirKLxnhcMlar1YFe5hEJm5xwZ7wmbq184m03am6Hmv7ffEn7X/1zzP/iP3n4Bi7ZzR3OeYjDAnurGG0WtBmsx5Y74yvnhgRkSyjnzTgrBIqcVfSjZj9W3I2K7ajYjsJuhN3oe6X5FNbgYQAXo1KTHEwne3r2yRn3LqqH1Rg3MLpDijJNmPFE3HV0faJzWYbQgscHBhGaSnCA2ExJjNVkbsxNJ+N5RUZODAIiTDAZuyUcUJT4JisdePx79v1kn+Wzdfz8xKV+aPaBOtylmxs/C0PBgX/LKqYCG4b0deJBlgpasagiYlzRceeoNy9cXn7f3GUqc+fxrwxhKJwojD1gbZ/pN+EhSAk7POtCS4OyFVbGFPVUeo01Bw72m9A3q2FUU0giat5OVKykC5GyFxHfqBcYBnrJmPJOblBoVu6CqP2qEGrb+M4KrkUCx3cfgoGVVlTKV3httPCsEW5Hxbpv+OxfXfHJzV9y+el/i/n7f2/hnp02ef0Kt+vRlUNwOOcTfl3o3Osdr0tJrc64pKdgg5P2zAtPaYtsC+dgP1bcjh5C2I7C3Qg3vWM7GnprsDgGGfxLSfnr17FjpONg34TrUqdklnE+ibrWL/zLLDBU4pTfupHBZUgpjRjr2QI6YKT9eItzNpfiBiZBSojZGdecGJXlQh9r+yOI4DhyU2H2546w4qirO9HQnYjVzC1q56Y7l47wh2gfqMOFfEMfusH3AfvHPMKMS03f6KlmnLKOPE/Nprhv3p4lD1yHb09eMioiu2EJV3bYxAuGgO2m6WYUvBlRUT4xSgeGirPBDSkJY+gCfhySe4nTWUwlA07dyZ6KGuu8WLcVh3baF1s4MBKFxwntyhUrK+gaLirHVWW5qke0soy94P67f4fuupO6uIu2Xvs+aD/+kh9/ueNgNG/GltHBwXhaVmQRjNbT4KoQztbK/6sUKYHWjb5TMMDdWHE3akbn4YjfbA3b0fCN2TMwev5zSHT5e2jS9alk5RkF7Hw5rSh8Y0efyIpMAoDY7ywWEUylGH2FWaz0itsvNRDSOAhVZPPPwpYmn+XGkQV8kJLMxfohiPBN1qalu/Pl07In+Lbu6Fmc/56XzN89TYf8gTpcd+L3+5aDY+ccGZalo/Of3YuMTzAti78NI74/lE9bxQHrouxeSoAooq7uPDkXHXGcqsY+Vz7q8cwKa0e0NNRVm7i/voijZjQej9Mqf6dVi5aKMUASrfaprN7tfHVdUe2kqVOnCMNAHTQdBvrgcKPcn6MR30rcOB99UvuKq8vKcV2PXNc9SjnsoDj8v9+g/80N6p/8a/QvruDv/Rnub/5N3/3hhLl2BR+/oP7Zb/nkN3ccxoq/2LbYIpoF72x769hUisYjNmglVOJ88ixEtgfjy2+dE7ZGsx0Vg/VR82+6LVv27NRdeulFypYqnIwK+sGxn5nHZbMWRoQQYnIrwjvHpedeEa4sWoCC0hWKcY7yLhPHGGdQubTXMwdGos7BBAZwdoFT69J3E0huMotTkKLd+QwuOuKFbaZtnIIbnqbT/UAdbmnziv9Tn923ruN4nWl0OxmEqbQxDtjDdLkZhWayHcmSeqn7Q6IC+YRc+ZBaQrQdLE4/y/YtqRxZZQnAfCYWE0s6USkq80m4mkamTs85ixGv7dCyQTsvtuPEcucOoeFO1J91XiRHK64bL9V4oS2frDp++vKG5z8/UD1TSK1xxtF/NiKfv0L9d/81+sU/Q3+0Qv7sZ7hPP8L+B//B0d2xf+OvIZXm+eFfov77V/xmu6FWNa96r2vgHGgttNrLRWZdCHjZWK4qwyerQ8Jxo67C7aC4HX0XCeN86UlLy+AGhkDjqkJBSmzLM7g9UeGrkY2PTt3gPy9nOVgqvc5QjtlPotjI7zbG47ZlYY2xPVE6MXYRgSz1Ob3Xq2KbvgTXj58mj73gSNPYLaleCA59xDSYwAYzR0wakeUL5NTzFT9/qN3607KzwwXyW7R8a556sz60nROCymWUMhmIpdJqoWSWItq87SmX8Xj7RyRzYsVQcYSqjGjKREee9k6J8UF7Fd86RxUPrsMkDDluc77/sj16H0R2NL5rr44JtcBUWGlHqx2XTc/F8576Y4163viqsLuB/lbT7X1n30pbKm350ed/TvVHn6OUwv38p7jL63zsz19A26L+f79k8/WXfPLrAzujqEVhxd8pHTDli8rDGXUQ5bmsvMO9bAacg95orKvorGJvFQfrK8mMdWhRKCdUaAZCdE+LcgojAfbBoBxEdbgK30wyXt8yMeZhGNJQjMLw+Z7k0vPIrklQhMuRbez2MJneh3ukZy9X50ZiE8nF4oS0gUj1Uqm5ZPx7zqvNUEL+ZBkq+HDsA3W4ZWTqOHaScaCcepPOB41Prvm3/hImHARyInOhiApS2+kZ3SxxcWd9qsA/HLZ4QIzdIVJRaV/c4B+8nH1WZGpZyfH127I+WlLxQWISMVs3Ys2Bdf2CWq1z88PggFMUV0TBKy6D/OPax8RFF4iBkRs6LtyalornVc1FpVhpuNCO5/XA5bqjeW7Rn1zCZYv91Q3mjeHuZsPdvuXz7Ybeen7sn/03r/j4l7c8//gv/fn87exwAdx6g/zDv0979f/lT371OfwO/u3dM+5G4XbwGhHe2cNKwUeNYVMZfrres6pGVvXIYBTGKW7His8PNa96YTvCVwfL3WB543YM0tPh9WJrWirncdkDWyyGjbzIjjf0gotO1DFlIQzjNvWdi5rI6X4Zk6CjxK8t8V3JAkrxb1XAT8fc3EOe6STnG9ctOprMo92wbVISr0aUX96larZ4HLoY76eerfkz9XSj2PvsRLh0tocHxEMmy+scRYHCKeL3MseXVFo52WxJ5VkqzJht422/sxOKUl42Vj1FU6FXr2/Bk88rxvFe59Uw4vUCVKCJ+Wo1X5SgdTgXY2E02L3B7OHQV2z7mtux4lVf82VX8dV2w803a9zrO6+bsGD200/h05dcvBh4tuq4ql2oeotUNFhrx0VledkMfNz2vNgcuFp3rJqBOhxP7OZwMD7xth8dB2MZGRmLxNiStc63KPLXIESyyfmZdG3Ke+FZLbnrR7S5Ti0UeG8BB6TlA/w0L1Zgdi/TegvjK68X2qXfs8zcjsf4D9OZPsY+0Aj3sbhs+VkJFyw5yNiptIhY8VO9vNcT+w3TsUwbs9MW1SG6KKPknMCwKdKN2g5lax7P1/TJM0FRaS9SXZZuRgGdpOAfsT9VIyGa6oY3dLyhrZ+F6Mwn0np3h1YttWxoQs8zi58mj8GRDBI0J4riEt9SMTuMmMRqlKXrK7pvFG7cIdWO/e8qtncrPr+95NVQ8/mh5s2g2I5wUa1xCD/9yy+o1p8jf/YGd/1sen3rBvtHP2P9d/6cP2pf86c3V3zZNVinabV3tn/jouPH6z2/+PQ17eVIfRVoYBZuPrN8tV2HqBpuBnjdWb4Zeg5uYJCBCs2F27Bl53m3MgKj7wPnKp65FxzYs5M3E8eXpA8DFUtJxUX7Kb25C8UMviddZ26xdkzXcCKRiJqwEKKDtfaAiIeRzBiphPEeB3Wv1GsvJmcpEmUZs51CCDaN5HxPbUjWTQuFoov20e38+YmwXZxpPmRPnzL2gTrc+yziTI/Jgr49znvc3TR+4QsqFvch+FLkGYab6WEZV00t17lPjGf5s1IbtcRylyIqLV6hLFHJKLVZYzIuTJNlpHJRpUyF1JtX922U0IRIUwQcwmgUw17jjIBy3N22bA8Nb4aa20FzN/oos7fwZtBsuobutwZ1/Rr56qtjhwu4lx+hf/KM+lXHy7ZnsIrnjWdG1OJYaUtbGeq1ob5y6CuVBG5E5c4VvRV649utG+ej9Ng8s6aidnXqBZeuh1gO7Ok4oNBpxAga8VdxSvWbYeiptHfGSBCn8z0jQwlzedE5jBQ/Wxof947p5HTzMglmeCDaPRlshKN/nBN9uo422tnhnrT527d8O88HqmPa9Xd5gMXo9CjZMPs+Y7hlQ7/sCK3r0/f+wcldK8py0jSNjFVqkjV9SzzXYTyJPpT/RuwwTlG1any7dXJrn0YuEwUsRrzgp8s6pMpiFdpAh0KxcZvkcNfUbFTNdaO4CvzbVnmZyF3foG5JNKxX+xW3Q82vg+zhN50E9TH4/KDYm4Y/+fOX/OzuDc//7JeotsX+7OfTi183uD/9Eyrgp//sFY0O0bcTOqO4qAbWzUBzadFXCv2iwR0M45sBZ4VKWQbnsdvd6NgZE+Y0Qutaaio20mCd167t3AErHrM1dNzyJUoqGjYJy43XLyp4NeoS5wxjlPgsKsMSxzZ00K30msFtQ0cRf38Gs82VYqJS111je/83TcFAKYsbyrGm0phZbKIaXv6u+KzEdiUxF04ED0fOfM7GKZ+7px/Rzu3scN/ZyoGzDDF4m0IMeY0iSRYHt8wH33Q7JTG9jGzLSKZ0vDiOstSpU6sUvavQkyio3F5M3pR80tK01Am6QPwReW0FXyQRnb4Vr26lXINGs1E1F5XmqobLCq4qy0pbRBy90dx1DYdRM1jF54cVu1GxNcJulFyOax1745kOX+w3tF+NPPvll8jFGuYOF3A//QmyP/Di55/RfDlSKct+qLgdai7qgUpbVAvSFtdLkarLvFh4fpkOeKe7ogkt4eMdc/56OehlTO3mcdDJXXohxkKHqF8bE5Ke+eFfbrH4IVZ/WZclPsvkWXxZlqyUCY5/FOWWiVRmnF1VDOnyvh9HsaeFZ6JeyX3Ocs5gWPrsh8VqODvck/bQYIHpgFgidRfJgnL6NivpjeplEvC0knCedUqZRL1J1IZZVjodmsXJsSRkdKhRAEdEYUyfqsiiyAlEB+7J+mUmXYpIObVvCY7DEulgTLrNeqfjr2eF4qLSXDeK69pxVTme1YZNiDo7oziMmjdDzcEq/v225mBIUW1nfPHCwTia0J3is32LcMUv/u1vaKpfoz7919g//bPpJXn+AvvXLO3f/ldULw60n71if1vz6s2GdTNQ1yNqJUhTXK8KrFHsRk1ncxcHgI4ehWIlK1Sonot4rKAD6j+kxOLIgLW2aHszBEW3ChNFZkRTSZPGi08vDl74245JF8GPjdzxt3xhLiVV8/2P42aKI8d2OJnJoB63zVQ4EfFkScHEshaDP/KpI31sBLtUav+07OxwF+2h6rNy4EU62PzzKRE8yUHKccQQ69Enpb2L+Fr4Lu4hDnY7Ege1p4z5PlKKrIWQpvyBMhYFcJTU/sEtnHA6U+cl/Qb2vvOr6ODAdZJ+rGWdaGEmKFtRcE/9EWe4paVmLTXrSlhXwoW2XFWW5/WAEsdgVegJpviqr9iOirvRswK2gwuiM/6aau/fGB180SmElle/bHmhbml//mvUxcURtOBefoT8Z/8zqv2ey+2Wq1/9lo//9We4m8F3hKh9g8nxy4HxzjeZ/O1rn2S7GYTOOvajZW9G2tBOfnQ2tCuy7DnQyz61u7FBLL6VyyBp4+UY+0ITQUSjyS9XMys8SdddlwUSXjiorlqvABaqy/x1aVK3hvQyjlGyzZACKLRa+fFBlZzu1MI6MVHrlkWdoJi1Jed8CiMu14t/Ryda5lCW7O1yJt83OzvcB20J0J/9PUl4zavVTg280gJuVeBmOcoom1lOMTUpBu8kCeJGoCkezhmlDO9wNQ2WaWQ7T6S4kBjTQSgnE+1jVj38LRbnQsQW9HfzfrPDrdG0SrPSkoodVtpyWY8MVnEwnl/bW6/E9WbwSly99dipcf4l5tXFxDs6K9yNQqs0r283rL4YaL54jXzxJfz4x6BnOgIvP8p38E//DPW3f436d79Ef/4N9otb7I13todXNV/fbfima3kzKB9lW98CZ8BQh/MaQ4zvdRQGRobQJznPPhrWDCgGdqn7QrqvcQpfFDrMRWEim0GrJkgu+pdkFKpJDtZl8XpL5mNny4k1kbxdd0/eymO0cWZTcm9Pj+sjzi7zvMZ83fkBlFhu3sJTt7PDvddOvWUzrUUmjjW+nZczvZ6/aGaO0zuqRCMriOQxAvHrLrELYuO+GNH4YghVrZIzTg/UbAqpJRdUaNV6mlfQUqj0CqVqWn3lW7wX5bsRZ6ykRVD07MKxWGpZp+7CMaqN1WaX7hnP3BVXuuGiVrxoheeN42UzclWNbGq/3WsHX+w3vBkqOuuj197CYfSaB1okKIz51j/GOTrrGKywNcIvb64wVmj/+Tes179BfvZT3Kc/OnEfw5372c8T5iuvX1H9V/+c/rNf0veaz7cbvugaXg+K/ejpaxdVRWMVB2MYnaVjzLxiFK0L2rUy0nKZhovnKNdU4l9ug9tNCkkqaTBuTA43/j0XA4+zEaX8rKPSq+SInbMMZpupfbH1eVEsM+cKx04lflwWnahV7HoSu+1GZbGWFCBE2KFI+ObE2onE8WTfc/plTFLb2Wc/DDs73He2YtBOsK35m1sdDbJoMWGRplMxQ5OgBlV8PmUfxJ8xKeaV+INIeYHvTrbjPwzL1CnKOkqiSWQneGdbSZsZDUEDN5brji5CE1VytJHgnzocO0XrWjaqZl0p1lqx1r4FT6ssrbZoZVHikiJXeuyc11wAr02rQultpJKpsHzSOzCK276hv1GsugGqtxvi7vkLzH/yH1P/N39F/bVhdD5JZ53H0X1reN9WHmC0Cmv9lR8DhqvEd82wzlIVj5jvotym6xThlzg8jMtiQH5GUYfPx8k9FBSisqaGiNeqwPqZBobJd0pqoupYKW40GUPFWIkOPi5vk8OP49JLj/rEbCmupPKYLpzu23V0WFr2h+N0zw73yOZlv0tLxGm1K7DbKRvBBUcrhUhILJOcRys+YphStZRq8O13DuEh9gIknmfro5W4THSWEQeMBPpIH4rK/+ALHuJUtDwGgLq6QFDUah2i4DqwZXPL7RXX/vjC+VbiW7S3bh1iPMMz95wKxZ4ehVBT8Vyv+GRV86L1IjU/Xxue1YareqRRBiUuCHvXWAcrbdES26z7B04rlRzepoJWeUiikqiJ4DgYxcFo7KCQj6+xL9+h3XrdoP+3/3Ne/l/+H1z8ZuTKeM5upYSrWrgdYtdizWDhVVcRG0PuR0tnM1xgijF0sAOvZe17yVH5FvaFMPmoPKOh1/swUgyD22PckKrSYueOibZCwIk7d0PUw1VSU+lVYDD0qTzYFsUJEUZK7IcALXXjK8DnAQDEZa3c5LDD8cw7PpRJN+TYfXrYY+DYpoFKzosY4lupfO4IV+ipOeKzw31rk+K3aRFDSfXKhQrhu+RMFTIZhVPctOS/Tv5OknkAlZdGXMDQSgesAj9WqxaNZxrEjsFK1ZN9KHLSrJJVmOK2yeFmMRaPWMZps6CTA/GvGEVLhRZF7TRaFI1SXNcV17XwcQsfN4ZP2oHLauS66ai0pa4Cm8EpWqU5iGOlHBdVaIlj/U8tvhT3KrAbWhUExsPd2FSWTTVSbwzub/3JO91hAPeTn6L/4f+AF//XX9IZzXbUVCIM2reNH62Pqr1imO/F5px/OQyF0HesoHPO0RoN1rcvr0VzCD3MUmWYc4yM7FTWWBhkhRWTZjVW50SkLbRz47Acg5i5iKLV10e6xaXYUBwvjiIqFg9j2cCGKKl9iM0v/WQeFsvbvV+3QwQSb20pMVwkl32EnCvZpp/PcyVPw/F+oA73Xeu6PT47IXqLpkwKScRGBXw021OqMM0rgrxDXAWMLbajzsmy2NWhCj2tYnZaqSop/PvPQ4Y6kODjIG/0RSq9dVh6m/UGtLQpihVRVHj9WkFRUbNyF/izsFy7K2o0n8vvAMOlu0pkf4Nv0HOlG1Za0xt/LM8aTauFi0q4rOFl4/hblx0/u9xy2fa09cjmWWh+2GnGUbEeQnt4J/xoJTwzirtR4fDtd8A/4j9a9Xy6PuAcWCfsRz+UV5XhR9d3rH6usH/9b7zFvT0283f/Lr/42b+g+dzguPIt38VrKhgndFaxGxWNqhisd752BnH47hHQGzgYx4uxSqPvbqjprU3C585BZw1vbG5TFKl0HT2+TXqb7ncnHXvZUTkP5+z0DQMdB32DoFjLM3p2dGqqMREjYvDO0NiOnm1wuDo5T2NCSbaalhFHTDjzgEmOWBJkVWFcj3MjlfLl5MZ6cR9Udv7T4/LlyInhIMoLtLsxC+pMnG2eVfqo/fvvdD9Qh/sQ7es+K1gHBW41L0Io9WohswMihhp/93u0UzwWhaiMwypVUVcbn502uR4+OuHYRHISFYeKMCV1moppajY6T7GroH0AHiKoXagmCz1jfblqTUtFLRqN8JH9GICNNNQheo0vnFb5aXejNFrgqvGi3tc1PG8cnzSG503PZdvz7NmeamXRtQvywCGCC1VdtbJcaEOrPLRgnVAr38J8pSwv247rVUdTj2htGcf80nvxox36r798LwSiqz/TuPGOXV/TVoamyqyR20PL7VCzMxsGJ6ldD3h9BiVwO/r2QdsRWu2pcM75UeQ7GevklK1z9FbRDjrNgkbrR8fO+gquVcJ1LdppLC5Vulm8UtkVH2EYOLClokWUSknM2Ci0dLiD7MLvfswM7H3zSpVf4rBKjBdrfefhNG6FSRI2bldcfi78hz7XgORAoqyME5eXi5aPswhq0rM0bdX+FOwDdbjvagtZVwkllDF6nYuCBPZAxL98VY+aDDBbJB5iBKFVQyyN1aql0Ze+Y6vNCZS62qCkTpiroENH1y6V2+rC4SqpeOY+SceuqX3tv4xYHGu3CqyChDZzLSuuqio0VXRsuEhtv71Yt/+pJJZ2wqoSaoE2yB5+0lpeNoafrg98fLHj4qJj8+MRvRbM3uF6sMbhrENCZ9xWW2rVY53QqArrBCWOTWV4sTqwaXs2m57NRwN6DVL5/Zuto/mjBvnTX7yXO67//i+4Gv+SH23v2Fz11JcG7d9z7D6vuLlZ8bpvknOP/dA+aTtq5fir7TpUw6l0Dp5n7HV4AbaB/WCcYCwc6uxwdqNjdI568BrCbYgOR+cQIxhn2NDSiMY4X/n2Ui7oneGX8u+padm464S5D9J7KCKUFVfUDLKh0qvEm46dgrVqEltCSU2lWga7pzc+Yja2SwyWiClPeMIF68HDVlUILAhdRnzXijLYOGb4zAuKIHak8F1PRqzreCp2drjvYPltm1kD4piUVNqir5QU9K4SA9M6Rw0Ra7XOVxVFPNVYr3PQjTchCbZCSY1O2gUxW0yYDtbURYTgaUgtF7ygdi21q1MXWQAjI9r5WPeSFSKCcy7UNxk6N8IItYToVUuYVvvt+4y9b0tTi8/kbyqoJIuKv2wMz+qRVlueX++5+vmAulSIFjQWKw52YI3Q9/682mqkD1HrVT2gxXHZ9mhlaSrD1XXH+pOR6uMKdVlBrcFY+M0B2VS4q8v3cq/d3/oT9K+/ZvPnd8nZ6rW/0fXasOoGruoxwQyx6WTshbapDJVSWDyFzbdn95CUEs9v1bWPjrcjDEDlclPFOtQRXNbThFHlBEVNbVVicVyywuIYnecFX7mXRD3ig+yys8Uy0oXtWQYOjByo8IyUUXVgveiOCmwVgMH6BF6koilVZdihYLmkqDQmzsJ3zhaVa9jUZRgyrhyLeZRqfEPLoB0RnetEhjSugyyEQt9POzvcR1kJzBccXClKIIU8JSJjqiVMADkC8FFs2VfMZ5WH0SJipwM4tMBWymO5tVqnaEPNRpoKxxS1a5V4CtfGXiaK0uhGnPThCFUAFioqUVSiGIIDN+HxPLghcA00WoLDJTZe9JFaJdBo7yBa5X9utGOjLZeVYaUNjTY0F4bqo5xhx4GMFmdhHFVysh4v9Wh5ow2VslytO0Q5nBWaS0P1wneEkIsG2hq6AWk6pNbQ1LwPc598inpxQb25QbegavFvE0A1lrq2rLQvse2tYqUNm8oXcQxW+Z5oyrJSwuAkULcEG/BghdcBHh0cRDCB+hbLh5XERKHfZ5+RKRoNSjSdtQzW0iiNdY6DGxPUEBOZB2CgS+MwQgtAqhKMyVItNVZqJESesX17HFMl3czajOVCdq6xL5uQx3G0shw5OePA6Y15Ds/KCZ2HiQ43O9lw5MSWQh7L+f6X/H7gDvexybP5dyXGFJNbuf597nhjOWWlNx63wg/SbnxDKrnEIsYnMIzrE7G9CvPXOLDH0ELbuiGVgMa6eBWiDy2+3HMM27EYvlGfp4RYPK6aliawDSyOb9wW6ywH2VFLw7W7DMpfPoI4GMNoHSKEKFrYDi44Xa+8pQV2EjrfrsCG9OJKG15c7LGD0H82ECG/4Y0wHjSvX6/Z9zVvusjhdZ7jitdW0OLYDTWCd8SiHM/rA3q3R5oDUgmutxw+hxVb6r/6DL5l0gxAXr/GbTuGg2AGQWmHqn0Hj8Ormu224WA0B6PYG1+WfDBedGdwwuu+orPC7ejHWiVwsL56rlHeob7q/YxBAYOD7WBRAbIZrAvO1//cDhbjHIOzKIRKCYO19M4kuUg/oTeMjJ59gmOQnp6dT5CGF3K40Kk0u2RDlDzt6GjjuBvNIcAN0yCiXKYU3DGmn1RRxuamdsYvjuaj30PYdtjPRNlslhfxW/3W9/q7sA/c4X47W1TLP7mMTZEwHMMOkQuZW54bfIfe6WB2zmIld3X1ouXekRps0GgaAlnd+G2KSlNIKJJrKDRVghc66byYouwhYLZRisW6QO53FhFBKR+t9RLwRATrQAfna5yEogX/e7S+0/AGRIfs+66iO9S82a/YjxVv+jpNtY2LVLCgRwEocR5e2LWs3wxUnUXV4doYYX9bU73pqb+5ecu7ecJev8bdHOj3FUocSvsXjnOw3TbcHdokbLM3GuOEWikGKwxWuDO+2eR2lBStHozv9mucd7jbwTFYH7H2oZOEVr58OZYyi/Pkt4M1weEatChGI/TO0DFignoZ4a4ZMRgMCvFt2wPXVjGtPHQYTMD+VcgblEyGUrgojsUIY8XP4nbSMpRjNgQGcXmJ+4/LHMMRx1oLSxaWOSfNnoq9C/LjyEr4GeDPBHAAlUjj5ZvX2B3GHorvSMvEAeoCPWw0B0R63Jhxr1j9lb4rBlqEG0bXY8wtlVrThB5nxvkqfoDB+Yy0lpqRjk7u2MkNCk3DGuUUG3eNxfKF+l04OsPGXdO6Nl2C3jQIQjWqpCVwpVo2WlMpz1U9tIpK4HbUmP2a7VhTi0Urx8FoRiustWFvNL/atxyMcCioVREHvRl8lPyy8c5WC7weal7eXXDdDDTa0BudGAIvdnt+uv4c/Z++w+2dmfsv/iU3/3zg33z2I7RyaLEY5x3qr3ZrdsZHrwcj3A6wqTQr7RXNRuedq/FFYAzWsR+9czXOi5hbFws7gAE6a3k9dtSiaUVzsL5suMInzV7jE1YNNXsOHGTLKF6RbON8UUon+6RQZvCOFOcd4hiuUVlSPNg9h/4VnbxJU3xZcIKJeqibDE3YLun1QsZny3VEKnAWG+ABNXE7MagYZ8FJoTUhdUGtzNBC7pDydOwDd7jvYjFEDdHq/O0q2QnP+395sz5ynZPHXWgGSI4KJiLPgbNrne98K05TCqD4JIefylk7Yuhy9hewKhPs8/YsXnLFY3tGfATUsmFkSDoJCo0SzRCiZE3tq98KNoPF0lpf/NA4lYoCDlZ4M/ipto9UPUZ7O2oGKzTK0VvhdwcfCR5Mfg1WAS/ejrm7rgqfOSo6o7gdamplM7Iewqhnvz7w/N8cSzS+rQ2/OvD61SXf9PklOTgfvX7RKQ7GU746AzeD8/Qv5cuBjXN0Jpcme2lJy+gcxjl66yGA2nrYZnCWwRnu2FM5zcFVDHg8NnJyd4FTe3CKQToO3CWWy1b8/e7cnXe2bpiwVrwT9Y42SkE6DKPt/FhyU7pYdLixHVNZzBAt8nFt0GlOCnWF6BKUEfU0+oX5c1ImwsqKshl8kBz6U0mXeTs73LcyTS7ddeB6ogMWdMEfVMVbuJDhw1eilb2l4hved8fN8owRw4rRsKeYjVjmVWfeIi7sj80wDvuJw400n3XzcUrUHWWVA383duc1tgsJug07XvnjEYWi5lJ5Pu9Al0p775xmsDUrV2Oc5m7QHAy86RXPG4/L9qERo+/U4CGHwcJXB48PH0w8Fo8T60A70yL0ISM/Wmi1YlMphIpKwYvGeUwUx7OuYTvU/N3//X/JxX/6W8w//kfvdrf/xb/gr/7VFX/x6hl/uatTa/TtGFXDfNR601v2o+W1OSTMuxKFEmFr++Qw4xgY8FP9vRx8S3XXYHHcyjeYEK36ZVWq9BvxcI9xA5aBIWD55UvUO05/b+M9duFFHst9o5JYFCrSukkOT4Kozmj2E8jBho4RWhqM6wMUllkF/iW/905Sb2YON2OuUYNhWg4clisduSjf0y90wnDOTL6fOmL/yVOxs8N9K4uPToynZtOZorjBuX76XVw3OmJi0UQmy6flimxsVGjyn4cpGjaQysmJOjLtxq+X6WbJ6YaH0LqR0ewpxcajqpiPYgotADcirkufRZjizn2d9yWBpYFmdC2101gUg/O4nQ40stuia8Pr3ouIG+d/vu6HhEVG20iDNh4TjewIiE7aO+/IAfaFEXBRgRbF3Vjx6rMN7f/rl8g/vueW3mdffIOEiHl0sBu9s70bXNi/87xZYzlYLxXep0STwTqbOM6lDXQM0gVxcZPEfobQZt1HjF7w0ZaJqGKqXXb4tUVSy2G9wyywVutGnLGI9GhXtN2B1CkiLuf3FfMIeQwKla/lCgmvdO9n4zxBA2UDyrSdsGzxnBx1n8aRpE4XnWzeXmYNPR1xm7PDPbL7bp6nnuT67mhqOiDcQBwMsXghrxshCRO+m1apzakv07/jwCUtT0yswcThZpERJtsazQFjOwZzR5Tqi9FPmfTQqqWuLnDYRHQHaPUzAO7G3xGFbqwasGIwMnKQigu7YuUqeuOTaSvtHdY3veLNALeDjwp746fVnTN8xWtGMfSyT/u6sNe0tLTGKzVEhoQO3RV0KrqA215otfBp4Mjejpq//OYZ4399yx//H//P8A//x28FL+h//t/iPvuGSlt0KFbYG3jTO256w8FYdjY20MxRa4/v4nsrXzPSpUq+GLXG5JSxmaJVq80EYlJoRtfR2zsGs01cV/CzFyU1dXWBtUH4PRTEqFCdOJhQOUYUpM/RprWeTRMdY5L0VM2xMy26QPsIOESuE4H0mNCUtM5EdD+9wC35ufIl8llD2nfIiL+DKyAGNVs3rp+P4HHtfL4fdna472zxBgtTfGmaTPODb149Q/7OmUTrKjHbUsRjbpkInqMc/7PEzbxYtLMZVkAUzkTaT0Oko/mIaphs37qBYSx0VZUvxdybV2HSnKGK+FPw5cHx4bPOl6ZuRx99VuIrp/ajhw56a7m1HR0DO3XHvENELwffTt15DrC1rRfFER/xGhP1LOC61lTKSz5eV45P2p7nbcd6NTD8RU89/gvk6gr3k58eXc+j67u9g69ewWi5uOz5aHvgJ/sVB6PZDp6VEUXIIwXL4oVnejkwSJ+m/0PAwaPZEJH6RKa/hlHfIrcw0gnSERRaNxjTJ2cLpCi2vGfH40AFKGCVGo/GIp1wNDkSLRxcbpke3ENsv7Mwhu/Tt83R6hKLIEINU4ZDfqbC8T/Kvv/822hnh3tkj3lLlm/T+fJxIMXBVP49Xy4ONiDCDXOS+FHUG6CIhN2VAzfTZGLknKPh0AUg7ENJFTLX2enO9XNHk6PNqLk6mn0ivaekSkggqoA5RnPhiA7GR4i1+OTSIbQZ763hTrZ0HDhwlzBLf5aaUUZGvGShpkaconFexNs6hxJfvSV4h6tEWIVii2dNz6YdaFpD/xWMb+7Q//b/htooqr9+jby8wv31n/luEJWGtsU1NfLqDbLdwfYA1tFcj1y96XjZDPyu85iyczC4CAiYUMEVoAcZGULC0mET1qoktzuKjSFthEhsnz4r73kJCRlyRRdkKEBCd2QPY03pg379Jux7DPCEnc56EuumKfbtx2es+opJP3FLjpM8Zo/G6SmMdR6oLD1zj3W2T8vODvedbD69OTU4ZjXgAf/Nb/P59xDVkHIjvmm2Ny6fm/7lQZ5bqk8d58Qi5Sc67Bj1FMuWWekysWbtmDRzM43N0lSXoXS4QaFopaJRilYT2uD4K9ZZx+1guRkG7lxHR89AH/DMIQhze8Frw8jGXYVmi/78hhiFO9+qR6MwwQnUSmiULya4rCzPNwdefrJl88c+t+lGYf+F0O8r3vxTh8gbLtZfsrnuWf/IUf2kRV6skb/1c9zFBvnFj5BnF6ztb/mEO766u+CLruZ20PxWBIOlC6yNjdtgMHQSklZkTPVCe+2K0XVYGejNnU9MFnq2gwsVVip34C3ZJHOqViowsCOODDfEYgTNKhUX4Gxy1iUdcSoKo/LyRYBwXNmVx1C26GzjmKkLBkF0unPq1pKDLas546v61PLff+jglJ0d7nuwuS7usr3vQTKn0kyx4Hxs06l/PprZ327aYr3kYcoRZj01DyU0PqPuvNKVKgRtgMQuGKylcyN7Dr7/l4xp/z5K1imZJGi005Pn1VdQhc6woThDi0KHcuJaHI2ytPVIvbGoqxZ7O4B19PuKu7sV//7NNcYJ1WvHx68O/Oj1HS+2e+ofj1S/6OFig7u69CpbH91QP99x1XZs9IZGaS/eg/hkIaCcYBdeumXEHnuOWSy6SFaWplABnFi+N3EWstSjriwZj/vOzrBkFExhp2zBuYYZ0JwCtth9+pRNIIr7bBlq+yHb2eG+k80AfkaOcael6VKe6kOJf2VhkiiR59wQkm+xc0RuoZ4Ux+I2nYGkzjR9ICaaDjCZWsZ26wnfKzpI+GWnmg5AwhUbfeF5nM54jqfTXLtLNmFqejCW150NzslXVfXWcTAGh+MgO1/RFqyWtYcR8DzjuuytFa5TjH53sqN2DTUVz9SKi6riqhYuK3jRGK7qgdVm8E1mDwZ7ZxnvhJvbFa/3a353aELBguL5ruXj7QV/z3zJp9WW6psbpKlwv/gZrtLIR9fUf7TjZ1+94bd3F3x2qFhXitVYcecqjIzcyF24E7lTcS2+ZHpwvlV8JZ6JYJ0vtyVIE5qCNji6DC1EXePIKoFcVODvp0YrHYTlTWqhnjRlUQsFNjkJm2ZDsyKE5a7SNckZlzOwOP4mCWNTrFuOxzLnIeTZXrFvYNq9d25yz3dPw84O9zuxckqVeYnHS0mxDCeXe2jdOc3MLzATe3ZBfT+ufyJySW1NygTdZL2gu4tvkRgZBCKe+O9lB6E3vvXM3g3s6TEyTgRU5hbb9ZhwXJ74H/9z6WdZQuIf51Bk0GvqrUG9MYx3wrBT3B5a3vQ1O+MbTt4MftqrpOLNfsX11wfaL2+RSsHLF2BG6MeUaKyVpRZCx2GNGhX+VRd5s5qoV1G2hk/9y4BKGqLaXInbzi3yZ8sINn7u8dncJidhw9EphiuYotVZVPwYi2IyficFBXGJxnVk832dOMe3dpxlQu1p2tnhvpMtZUXnDqvEr+LbPGJTU6dbYrYpgogSkAuJNKB4UGPLnrKlSblc4FQmBbKsT5qOdC4ejTqKbEuup9ZNmhrXak3NiogMa5FUsADQG8fX455bdtyqVxNHG3t2CYpWLrEYXIhwFZq9bCmF0QfpF+ENn5jzcpDGCbux4us3F+z2DevXPcNQMYyaf393wdd9xZedYjvCN51nUHRG8e9urhiN4m/8069Z/2RLrYNT/OIN5ssDhxt/3a4qy8crDVT87k7T4cJxKVqXE37RWrnEMLCzr1CiWOkXHs8N4kMlfFNJm4RirBsZh0PozJslEvvhNkW0iTvrpjBEZLAYu0OpFVptvOpWOWYKHH46Hgq4ItG7xsUxGNvdlJZZOXMrneUco2Xh71M2z588rWj37HDfm80HzJQeNv18aik6PYV7pQfq1KCc8ngzlhfoYxR/o0JGOzrtMSkyJawQC7EIApMisahD5ZzXSVVSsZJnrNwFG7dB4VkDCoe1DpTymgHFcQ8cGNw+tVmPU21/Fhoh9/NKV1Zs8Vz5fRsZqZymUYqLyifMqoAdA/RGYzph39d0o6Yzmm/6iteD4uvOF1tsB//iUwhfdhWtWvPiiw3PzY7rv/oaUYL5co95Y1Fa0WrLpvIyiK0WLkM344NsfYcFsXTsGOmOHJkSzxAx7NK1nNxiN022zS1WjcXrnwpVUDhhURq0HGvRwS6V5+ZlbfGCV4muuPjCT5SvqZ3OZ8jC70/LWb4POzvck1ZishF/mhOs55nVYm3RKVrNW8yDsWw4WdaFS0oI5c8gRBlpQ8fUsfL3WJ3mM9sVCFkKT3KbE9+JYkR0haKaCo/IFEYos+axDLSWDRfuOc/cC56pVYAQHNYZlET6lE2UKU3N6Dr25hWX+tNQ31+niFBTUbsWIyMjQwIqfPPLTFszjIjTaAzrynfU3VS+w0QdPM9h1PS25mA0uyCf+PlB8aaHz3cjnfU84IPxznhdKTrbstbP+LRraP77bxAF/Z3COYXSjnU1cF2NXFQVeyM8q2oY4VU4pgN3DG7P6A6ZvxzYCCq8vDpzmztxBKeZ+4UdR+9xOyXvNmnRhtkADlzpz4okWelooSIKeufo1qax4isZXXaNsYVN4YzzLG3JAgZ8knt7jA8f29sUMTw9h312uA9aefNPDZL51EZmzjYmy+I2po42rRWSUH4LJ7AqmUazMGMYzKaNedsx+i1bXccHqmAmhAg4bivro4a2JiHDnkRxxDK4gYOt0ZJ1BAB2dqR3IwfZs5cdh6h0pS59cQVmMgcwDPSyT9vuQqQYo+GWTWJEXLsrnqkVz1vFZeiZtgrqYr1VvO4bdkaxHTVf9orbQfjd3rEdLbfjyOAMewYGU2Op+epQY5yw0Wtuhhr+HDZtz9Vlh64t1cry7OKARfjs0DA6zbPG9yNbmQus+O4YPmpfMbgd1g2hc3Jme1QStTF8x11rh9RBWYvvS+bEFElOO8FzY+NFtwRrJYpf5RFv1yfHOq04ixYjWsi6zrMuuTNqmB+XeiGKzctPl7HHy9wLOTydIoZ3sbPDPWlLb8/HvFELcZtFnGpWiZPWmibW3OKgDMuewHRPJb6mNKIRqBBVIQXbYcJImBVEAJ4iVfSt8rqqoTGheAeGA1U0+9u5PmiRHejZ07ngcGVTnKlPO3k6WFC2wvdj6+wd1g3UakMlLTUttauoabiUlhdNzWXlW/rElung1bx2RnE7am5H4bc74ZvOcjMYemvzcUkQdzGOtlMYp1lrxc7UaLniZdvTNiMr7VC1Y3PR45zw7MawM4qL2gvq1GPtuySE8wAdeuwev6B9CXXQKg6luZoWLZW/vm7AiibyWiOEcOwsszMGDxlEASMlHmbwpd4ev48828xOKGiFztfJZfnDcRIMzMt3S8tiTiotm/MRpxgHNq09/e7pRaxva2eHe9IeA8gvFUDEQWbww7GaRQO+wmZOh5mWQZbOWk3hB8ooFqYUHyafl9PG6FBTb7ViGyp0UDV2qrFbNrTUqg2N/7wzXKsXNGzY2EtaWjSKRjS1Utyano6RnewY6NjKa3q3o7d3fntSU7lVgi0sxotfhwjaYrDOFJl4f+4DXmawtS2NVkms5rqyPK9HXzocui50VvimV3zVwTed5W6w3I0DHSNb9t4J4eiko6OjG0a+Giq0XLA3wkbXWOCTfYOIo94YDyuse142PYMVvmpqjFOs9y2d8y3n4/loaYsmngrLgAvnGWcOlV57BDnMbDpzM+nMrFSFsYEmFiLjqD0bxYuWJRPL+5iF7qc83AAlJLrhbGwdMRLijC2O8zgzug/fnb9w4sv4XaPYp4/9/jDr5/4gduJSSqkudsyTXZ5y3bP9o+ofNftpjx+UR5HQffSbopsicipVxZRUIdpc0bo1XlbGL+W3AZ2XNk+qWKPr8MpXMQPvnVJ0TjEZF5WxEsVp6So4hQ691WoloYeaj26b0NZ7dF6dbG98N4W98Z0S/HENGDGpy4XnTPTsZMcNW24Hy+3guB2F7ajZDzXDUGGGrNuw0oaNtqxDV+KWitrVGQbAppdHCeuUCl7+juVZAxAghjFdpwmTRFVHsNG0a8OC00y/n2IEzBkKx/b21C1b/CvtsXSuyNP9Ydo5wj1pjxloJcA/f2tHRaQpljulx8yTCPG7afLEt1gZi6ldbmAJUMreJSjD9n73qW27j2yN3RFJ8RHbA7AyTh5akYpKr8L+TVKjqqSlUZdcuZe+Cgwfg9eiee127DlgQgVZ1BQQ8YLp0bnEtj8Oy+iCLmto5+4wKPwUHcCJSdzWFZe0rKhQXFSal63XTvB6uY5aWT5qBm7Giq/7iu3olcle9R2v2TIErdlRRpRTEBgRZTXY133PwVS0qmY3VlxVl3RjhQ7O3DrhxfpAoy3f9BWg+dFqhXTCa1b0sks4vBLluy3Ecw4Si/H+ja7Hmm24xtOEmG/EOGZ2yIz1ECPbCBVMqw4pli2rxvJsaF5x5pe1YVvTsZwxWROW1XmMHiXSlpg5c375UqT6GCe7lKx+WnZ2uN/KHooc78NiT/F278FuKbm8s4EdSe4uwxXzyHYp0o0Za5mM/ZJOFJxqUAubWy5EiHHqkNpxG4akmpWkH2fXLLVnkeVEUIyuFZrKVdSuoZHKq4+pfLV8R1wvU9kFzd2D8dVtHSODDAx0WAn0K9HowJuNWKsFDm5AGeF2qGi08E1fsdI1z7sarRxRH1eHVuitcmwqYTVoVvYCI0nxISmDubT9+f1YviblvfHXoOBEu7juFPrJrJJ3mbT69R6eCZWztXexucPkxN9Pz5E+1s4O91vZ0gCaE8gNGffyy6bqraOBdU8FWmQUuCFwaM2kYkkQH72EUtqSrQBMVPYjXUhUVBDrFzuoGtujVUOl19T6gkpaL9vodhzUlpqGxq0wGA6B2lTTsGdLz47B7ZOztW5MHQkAtLTocPzWDYyWhO/awEeNZa+1rKlo2bhLrlnzadtyWavQrBLfFXfw23ozam4G4YuD41VnedV37MnNMU1I41W01PgCDis2ReODdLRuTXP4mN5W1ErT2RWbyrDSY9LHdaFlUKMc142iMxWf7D/GiqXjLrW36e326F45jJfNnCcp3TSSVVKhlb/msXtzxNktY2A2NPj2NqFysIAspiJDpTBNCSMEGmEhJJ6ZCjHqnb8IFc51zGdr0wTafThtGSzMcx/RytnjfJmn65DPDvdb2du87Y+zvm+1pzQtjFFP5upmZbGy4OF0RAu2iCyrSYRUUr78tnVBRbJUgTGgXYWmpgnT8Y6BMfTf6tl5ke3gbP0Db1J23jmLUb6PVlLAwoLLJb3OdVSqRXw9GxU1NRWN0qy0x249lOCLHQ5WMVi4G30zx+3g2I6e+nWQHXvZemwVTRPoZWOAUSKOHKPyAZ/4o4fLruaiUnx5aLmuFZvKT8VN6NtWK4/jriphLTVrt2Evl3TcEYtG0n2flfKmDgu27KAQ8GxLYBcME0ec74tK2yhLeif9whbHwCOi4GL9ueRihLfyAsVqkaHwiJlfudby8qdw56dtZ4f7TlaW68ZuDvGb+JY/ZeXbvSyqyNssM8RT5xzLeIO5WRkwXkdBqYo55neUSHE9qKiBGgjyMbstJMfroQSffTcMbPRHNGyoXUvrWlpqOgb2sveOAsvg9gx2Krxd7jvqxCqpaNRlcDBDghUGu8fYjkp9Qi3rQAdraKlYKc1F7btI1Aoa5ajF8VXvGzrGjhJvesON6biRW+7kNYPbU8saTc3aXeBFcG7ycQWHq4K27Bd8yc5c0+6fheOs+XGreNkoRieMzt/9lbKsta92u6prDv0VA73vTeaGCbc67idaxGijtm3s2OAd8JQGZu1Y8KfDZ25MDIRiyXydnb+3KvJ7YzI0VBqW42Nyj+KYdvNCjDgm4XgmV64Ly8GIm/2cf16ut0QXe7rYbbSzw30n82/ljM9mCtkxDjkfeCUXd55EKCk18ZvgwMsHS2LJ7rggCj1v0sfk9xz95LbT0bnP1aXmegrgW8QY8dNzW0SlFZq9eGaCsWE6bbZJByC2ginNOcsYHLNJrVuC85Hat6Oho+UTVm7DSnlcVZGd7e0o3CC87j1me9M7bgbLV2bHnWzZymtGOqJuwyhdklUsnV/EmhGw1FjxJcTfDCuaTriuNdZptkbxk1XPRhteWx/dv2ws1im2o+JgGg7milv5xt8NO4TIOt7PQo0tRLBzK1+Y8WfplGPEq6TytzUWPBQsiZItUIrgzKsUM9QwOYJi/fnnS3YqSp070rdxmAvFRE/Y0UY7O9x3tscMMpg63PsGzFLEnBkLWQIyTvnDvzTtgyx+47VkJ8mvgAHniqIqdROIhPkq9NYq5QDz+vlh9kSqfXJcFRpxoSMEWZhlMFt/HJV3HJHpUEIWY4EdxuorrVoq8SwKKwM1DS0NTWgmCR5KaBS8Co72de/oDdwEgfOv1Vep2CKK5YwMSKGvW6p6xW64NnzvsHSieO3uWPUVLwYd+poJP2qFVllsuE1XlWWwwptauKk0t6aZdr4IUWW+lp51EYXB50mz2MGhXN+E66RpEdFYBnCE6HUMRSy5qCFWmsW/vS1QEk8mykrMtJh9zSzPwOLP+0pz4/7//+19/ZbbNvLlrQKpj+62k0xmfvsA+/6vtHvOnt2zs5NkJml3S6JIAPtHoYACSEnddttObNacjFsSP0AKuizculVlnZK/vtf6GlsB91VmJ4cVgL92slzjfu0xY5bi5BbtxIDqZa23kr1e08MMtnbunAcM8ZyCJalVe2qHrZxtCBMCTQg04r77B3b0A+7jj+DIONAHjBgwUWqOSCOG+IQxHjIlIHfK5w6xQAGfEEYQce4CfPZPwiFDGlP2vEeHHTraYhf2uKMNembsHOG+Fw93CkCI8t/zGHGYAn4dT3jCEUd6hMcovGzikrf0rkoU8BiztiLogyr1E4vwOOOAwAHbsMWPQ4ezZ5w94X88b/C+7/D3zQTXCeSMgfDQJX6ZHDbYY0N3GHgrD6PEjzP38GGA92cwp6CYVv5S2iAH2MQ7zQ9AsNEyF4qBqdTLsFzu3Cx/u+zxAsirJq2qUHdxsCoZGVU9/+1521ojGkC+XIOkGcmVz/6atgLuZ7VbwHprX0s92ASHxN0mKVi1F5WeVNZsGicaXje3aifkwjZAHTln9NhA+FTxdTgBloBaiFMFtjYVdVGeli5vqUygeJ49OtqiwxYb6ZEgHR1YqoLpFUbY4uYBBww40FNKtpBW47nCFrMZdwJdI1mTa04JEWEA0YiBepzoiOfxPilBgD9Gedj+11a0vz6SSMRSSyEHQh+3whmnh6At+hOTXIwiZ97W3pNINtAmXmjpL1a67QJdTfnMvNVlnrXalup5Vb2XqS0Nql5SHuhcfamWVr3gW9vpsb8dWwH3VXbpaXyNX1oSgrem3kObNdaevU6AkDcLkMZ4bgB4zsFq6T3ABtTEM2L3kNUO+jlzh87tETDiEH/HSFLrYBv3GInxGP8Jlbn5MGDyR0z+VEBuISVZeV0C4zj9O4ORphCLnOoJd+4nPMQfcYbHBh3uO/FwAfVuCY9nCZKdfMAYA57oUcolxhFTlJbwAcLHasDu5P8j1c54n8omIutlnUnOQASO8T/41TGm4PGP0094P23wru/hiPC/j5sctDsHwt5FbFNx8v82/Q3v4zuM7ohT+APH8O98jwDk+rbaAl1TqO19qrS26XuZUgt0bQzpockKZVtOnTwcdymwtsTpm/mTtm9boNety0s8oS5ec1nK2AaUC3hek4xdC6j99W0F3E+yl0hVrm1jP7s2cZcmXQAW9Lyt1zuXBRXvdqklTz5O1oFuSnsYeExJ09pRDw9tETOC4ODT8lfB1tolUb1G6nPqMMqSuks90jZwqTW67HM2hzr5iDFEKbUYR4w8YMSpFPJOYCu8p0/XEZDq3OR04jye2esRYzziQI84xDtsvMPJ99g54ORTg0wW1tSRlIjcOcbJM7oo9X6ZStqv3Fzkh0xbL6K9V7YsZnU/lTrKkjFZBZUEiXDFObw0J9v2OXqyZvucYNNmmOmx+cKcXW0F3FfZNe6pJf8vBQNsmmO9bZEQpR9P4mwLpxbNJFddTz35pemMNRMYIQbzLv3YC3cLpGUvUsZZ8rZ2/d+wde+x4/d5WXwKf8hRU+eGMRyz5yje7TGlpurDwEsfNJOyChRPr5WvAZIU0dEWm7jHHjv83O2x7xgdE4YQ8XiUAuLnUPrC/hqe8ESP+BB/wRgOGPwjmHp0vE1JFn0KioX8AJkyzxwyQHuSOrwd3yHC4zT9Lts4jy3t4aLDr8ceY2BsWeRpAYSOIvYc8EPv4PcOj+OEIY7Y4QFgwHdDXpJPSX6nfKzSCno/2upuMU7QYuBL2X5KIykdk/W9frnjs+1LRjbFOwO3Syud5QLjcs+aeajHlm+32qNss/Rb+L6AeQXcmX3OSWADbdaunO+CHnfZ2khyrV4AUJaf+YddPCNOnR8IsuTv0rJX2+AAQE93WakwxVPiFH1WIYhXWUud2mIrbVGWOrtK6ihoVbAeTsJaqbi56kd9jDj7kP8eSArl+CTDUsDXbLUIaXbJJK3VbfYXILUcVEccEODjUJI0UpbcyT3jgB2epz2IgKeNKBciIjZM6F35ju67DmG8w3M8YKIRjrbwGDCFQWiOMAfDtlobkDzfBZCtg6PmPuevX3n/cm+zDheATZYp2wO2SNJ8dbX0QG+zLF+y6rP7rSqF79So+Tsu/H1tUtwSc7epj3P9bGn+qBIsu2QzyRHkUJfIMzVIqyEovcDp34AYzvlHRIn/Y97kYjVaYKZze2z4AVM84TA9imqAt3ig/wIA/BH+T2r7ooW0p8yDWlVC1p2m4bfBMelCIcEgpj7Xvt3gDltssKVOOv1G4K6TH7tLQb4xBhyClFw88CNOeMq61o428En76zGAyOG++wcIjDHxqR02WaIl4+hAqZPuMD2WYF76Cg78H3gasQ1bDOc97k477DvCgyf8sAEenH6TwN93jPf9Fn88S2LHwA84BY/R/zsDrtU5V/ysJqQknbTwzHWNWvneW9C1wbk2waYERyVZps/75LmCmsLK3UfyXIzVtvbBXzqcyJ63fzMvDZ59W7YCbrZrE+M1dolauFbIRj3flxwXZunXbLFYCQoo3o78CAs4LwTVgOStjqWWAe8TD+lwio8I8Dj7Z9HbhnNq0V3qJKhpDVelK3wY0Ll9CQ416gTJQDugd3vs4r2UPSTGzgmd0DHgA/A8BTxPHs/hjGecMNCAIx4xxkNVb1YkU30qfsP5emaeIeR6bXEZfSDkPmBpbJE9fuN/YYg/4IfzPzAFGd/ggWOQdN+9qhUIeE87+OjxH/zfAoZwIE2+sPy19WRjqpdAmuhQ9NOFt4V5iNVjF1M5YN2eaZHatdwvADtfZ5wuMWCoL92vAuWycfO6pSC+D89WbQXcL2o6wdpJuAS2S09+9Xx99VpNf5yyUbM0tctHKt17SyKFjUorrzkIn6pFywEM8QN8GDD65wwgPgzw4SQpwqYWg0biJ8PzAgBzL6m8EC6X4NC5HWKUtF52Drt4h54cOibsO84ANgE4TOLZfsABB37CGaIEUDmaUgHMfW5LTmD4WGe6AcjUggT8RtODrK+oCQAiNfMjPrhfMNKAP6YfAQifew7SORiQ2gpTlG/53nUYpm2+p/b7kPPbh075OQZI2i5zks4ZCZ9uIfuXv4tdTnAgk004ry6WVkRxBMzqanlFVttyh5Nbkq7vC2yBFXA/g92aREuT0C+832oor2gq9ZO8lDdaXR1Vpbus1QnqDfowwPEWm/5d+pF7+JQ5pv3MxnA0WtCQwRTElVcYETCmWq9AUT3I0vmYu9BmD5R79Cw0xgZ7RAroibFlxnEKOBPBEWdFwoAJZzrhiEfpJpE8biKXQTKEEQMGdCx6WA3waQ8xDVJqFTMiV3jm9CBy1GUwDgBAndSJYOC3+IhxvEd33OOuYxwnwibVeDh5KYK+7xjv4gYP8Wew6xB7eahM/mQ68ZaAWFV8Jn2nVk3i/aG6p0vzYV4ZTNLALRVQp/Re1/IKXaAyL0t9oTmGzU6TmbVabSvgfnG75C1cF3oXbs16yCUwIsGMhR/Q4uva422X2EwdOMvHpJwgQzw+bXyokqY2qp75Wvu+5WuRMthSdlUG3JRqrAVmKDIcExwD4yT3x0ckwJUiiyOdMcUBYzxkSsNR7aWHMCFyCoVpI8aUdisPB9+oAnz+GhSUCSF/PSFOCGCM4YATH9DFDs/TRlQCYNynBbgPkpDhCNgyYzfdY6IRHd9J3Qic8hjr75nNe5Z2SJ6p0gOLbZX0e5yMLMw+YO25XsmdEid3epnOWu1ltgLuZ7eWQliarG2kF7PXNsV32UL5bNaBtT5eW3VK6yp0bi8pvnESACdOkf2uyiDTAJl2G5DxCQMqcqc6aGZThR07USbk5ogjGJ2k89Ideuywj3d4hzvhbomw6yR99zAFPI4T/km/YKBjphLO4TknW9SlJRmd2yV+VigM8bJ9uk6p+aAFZtSqBA2S5IHJn3KQa4pHEDF+6f8nDvwz7vx/xxg6+NiByWGT+GYH4JkAR4Sf4nswCCd6xIBHQ39IwodQMwWEFWRnHLuhd+x3L9ubOZNSwIXjbUE6NJ6s+bsC1hTorbZl2Ad9cRR0m0s0wurtAivgfgV7aQrkkrWcW3voMH+9EBSTUSwF3QwYA9mLBZUKVfpfphXscWYa0jCjGQRotU+atG0hktbiHW3h0MHFDpT89QgB2xiBIQScwoSRzzjjiHMKkikfLC166si9jlHHonV3pYlk6VWh29r9ANHMqkfffubjIEkRGIAI9IExeMY5SBoyADARiCK21GEbt1I/gUrJy/LQU5Bcrn0784QrSdcFmyW23Jg/QEU7vK67Q7tCsx71CrZqK+B+cWslZ0A9IWv95FxetvSDKEkGBAfiHaS27ZhrI+iSNFYdXTnzhiqWZ+rQ8Va4zXiW//w5840qafLhNJN4abdYx7uk8TzPgMKHMwImbPsfcvCsc3ts6QE7POAuPqBL1/NhHHHyDAZhjAG/hicMJIGvIT7hafxnHkPfldbr9oEQ4ggy90eusU9VyTb5AaIPm/qhMmLyJeDHCShz8kYcMcQP+H/8v/Au/g3wP6Mf5bP7XppcdgwpuOM6+GmHe/yE2MmYT9MfmPwxH2/j3qVzHiou1z5EQjybehk6X9B4wsk7JRMgC6cq2JqlX7oKyiqDtkNJsVqX24K37cZwK+nnLQD4U4pHfT1bAfeL2EsnhOVm7XtW01i8DquXrM6WNZtsvJwl7eRcOD9Pya11tvKv5u63P7qmQDow8zb1kkKcgCCKBa1fQGD0UUrVlH2A5zBihMcTPWJUCVg4VN4roy3uUjhmCaRxvo1LdX4rD75JTVaOWe65h+ONAcAJ53jAkbY4xAf03iGgQ8edfEPpnEwAEUnXYerR8R2YnitPNYMsMSg2q5Dk/er2lwJm+j2ILVEFqotO4DnzkhfmV9LY3jyH2fM6z/sWyUW3zvHntBVw/zTWPrHt60sTVCdc6UElP6IB7ZKwlh4Vzq+VcmkbnKDKhTgkVcHZANyEOjUYADE4pQjLgcIc1MzSffJHON5i272Hoy0CPDgy9tjAQbxaJkJAxC/0O070jKf4G87hgOP0W76e3NmWikRtdmdZivBodlebHkvE8P5c2tkYCkK92s7tMPkTpnDEpnuHnvdCSUSPITwCDHTcY4gDttMWPf8IZFAu52IwNtgDDAz8KGqF1FMuNwJtKIQ2qCmgO69Z3FouWk+qVnBFelalhAO1VltNvd0+ccHFQ64TIoC59/pStc6tVPhr9tfxbNVWwP1sZl29lha4NFHa5VHr6V5QH2TR+oI+1/CEgHi/lHS4uQh58vgcS5DGR63POhZAq8BWywU2xVLy1Wor8FJs3KaVAsggtuEHbEnaru+i0AIK9wiiv/A0Ysq1bQeM0yF3ERb9bqFEHG8rD7VtMWS3Db7O9sr/Qoqy6/aRPDQBouO90CKpJ1xJRx4w0CElWRDOPmBgwpTQtmNCT+LBT5gAHNHxFt7tcnZc5rlzDVyRYFXdHBrPu36vPFTVUy7zQa9fr2veGFKunRoQtW3WyzliBl2dz63HvTTPryVFfArw/nVsBdzPbkte61IbnmuZbnVq8OVz6NEa7yN19UVehhYZli0qo/pUBVsBluL12YAOLcqS7KhlmV+VZ1RAAMM5KcW4pQfs8R4/xB+hOVCSIKqVDSJGOqdi4iOmMGAKR/R0D2LOyRWq+RVwEM9ctcP5/NYrhsfkT7kkpFIOGshjdJW3K/pduebJH3P3CnlwlACaox4dOowxYPCEIQQwCA+9w4YZ/bTBgFNSQGzRuT0mf0rXPeUHgtXSVpxtLjzfgGmVTWaDagYwNQiazlb2n9NSy8kU9byde7e48foWALf7fHuguwLuZzMFSuu1trzTJR5sPkGL51HzpiLbqVN5F0vjVam8xjMyv41xKqCln1uPUN5Tz8tywV2OsgtQFamY1foybSpaIcaAKQ440xEDzujgKv52wBkjjTjgD1EDTFKHoCRQnIqnma5NgVcCYyJdY+6Sjrj2fDfduxwcs168Ns8EANuNQbdrayD4cMYQH+HdiOA8IgX8Gh7wEHZ45zboWb67jgjvIRXHHHqMdISnAc5tAI9cg0IruNWWxp0CfXruzJlX3qxdUdSJCnU9j/rhrH/LXDNecWhb9VhrZWFL9jHA+e2BLbAC7me2pUlz671LFcXstgV0b9cd1eP0GWzbouDE8p7N/rKf2zbc8mZonJgCqiK9kgBPtQ8VRUQJ+ujeHp4mCSjl9yLG1KxyigPGcMA4PRt5V2pBnk7RFu8u5SFrgLeAq7xuSeQooDv3zn1ud+NQ9LG2pXlEwJm2cK7HgZ6ACLzDBmQiiR0xNtFliZg+LPJYY6IPGifRKkuI7GtrLWdez42qwFFO8OD8PS2Drh7HrtKA+fx8CUB+m17ra2wF3I+yl/Cz16Qx13hafU+3SRBEWhe3g6ZYKs0w7xQ853qr7g5ggDogTggQsOrcHudJPLopHNO2lk+UcYknq15WyIEexzsJ/jSVwmyWWUSQNFu3wV3/d/R0hwf6GX3cZHVCD4cnHHGgJxzoERMGnMNTSoc91neSWFKLEfJnuTNCNHKuKD2/fDinMQjvqx691Qtrs0Zd5ss3mCiFRMPoPtYjJTDgS4CLHGOgI+78FrvQ4Ye+R4iAR8QZHmPqBZe1wcBMCRDjBKZNBkn9Puz1l3+7FMg0D50WtVXni0bdYgJobTnG+fwGauBdmufX6LFb9m0D8gq4n9U+dvLUoD3P9inb3PZwbwjdl85uZUpQ+dcFozpw1garZqNJjRNzokPsUHSyAWMCo5EGTBgksSGOufQjUECt9Uqt2iIPD7XsCokiabfL26sWd+FcFdWCYICyeNIhTind+AgmhxPO4EgYg/C6Pu0rnSgcOtpJoE+DXIYumF0HWkBu77XyuSq/4/rvxmytjWKWsvq2we9r2Aq4H2UveYK/9ClvgxFL3oRd0qH5fK6trbewsrDJbCt8KlMH5i4VUjlXHKXsk8AsnIHE09ZAUMDW5yQHBQwZIkUgJqdeZVYbvseO3mOLOwQK4OgxpWI0T/R7AiSPU/gje7ZVcW51zGA9PDGpq9uh433mZ1VHrJlwkz+k9kE7aNIHxUatYL6KGFOpxFhWCVk5cEHPO/g/MPEJv/Eex/geYXqPmGB6St/3Hu+xoT36fo9z94RxeobH2SgJRO6mmtxZwkN+kOi4RTkSMef0pfWSePoEFJXCIjWhlIKuxlqFzFvpXy+tFL9dWwH3q9vHTt5rPG+xukdVAZSSxutKEEhFDRVvW+RFNtmh5YGByUjG5g8AmzwQ4ohUfgY9ts14pZuEj2OpN2C8WD0fmwdA6w0GTODYz4JhmlCQ02jbe2UfNulcGqDT92yX4SWdcYij0By8RQgBR34ECLiLd1mFwfq/yAD1YDgQnCguArIXLGOaUOoo29VG4nktV34xEzcgF/vOlIIep6230cYGXgq2Kz/7ElsB96uaCsuXuLElT1jsUqRZ3ovN3yKy4pzaq3rVDQgu61kBAaoQp1RExaa8bopMKYEu8V0+XowTomkZE9Gl6HkZM1OHTfcOAHD2zxidLLm3cS/jyNfColwIT6nkotZs8Jkf1sAcMSOGmteMCKDImIAKpPXYxEUbrO852gh3G1PBGEhigaoVAtL1GR48qseYrk0TFyZ/SkV6znBugw/4J0Z3wB29Rx97bCHfg4sOgXoESDt3AmO3+Qmjf8bxfE5jTAkO5oFSepF10Fq4GUCVn79EI2RZWDtHil1uWgrUab9L/G6L+K0Kx4J4xPcI0CvgflW7pFe0E3ruUcRqmxp42x8MmSw0QJedSC1eAmLwVYCpjYaLlQCZ/uDb4E1ENxurLrldqnrlwxmd26XKYHtssM9Ae6KDKBIgioSzf6pojqAJG2YZr4Gt3OJbzwnkQN8lrlbuW0koyPrjnNBQgmOxAjz18Jv99SGQVgshjohegPQMxqF7xBZ7IEI4agxgMLrYoacdQNJXzafAXIghfa+qWNBxKse6sMIx6hHJOCzB1Lqt+ZKXaudSo9Gt2ufAfH6rMlgLtKutgPunsmsBrnbS2h9N8jyoz96NjUDXKoW0twJLZEzhkLoLbFCbVTbYkUQgnGE7RxTQm0vBVNPqw4BN9w4bJ4Vq1LuN8DjQB4w4JQlYKdCtUiy5Js7H1OQAxzswdZiC5TQDYpBebW1R75YKqKqLpZRZXRFUCgITgLJaZgtyVg/swxkxaDsdxhGPCBTg0OfSkrt4D4c+0yqaRszUCXebwFU68pb2SUUHbDndwtVe66grD2DX0AoWGAtI11TULfrqpRLIa+9/+7YC7p/K7KS+lAZpkyhqaU6MpVtuWbyZYBAYLqXCZh4ygQTzJgNrNGC6rPdUTtBG02U5u5SBptlgamM44sC/40RP2OAujy1Cgk0CtGMKlDX92ayCoEodLvxxlrMFABQKWCUv3bcyLBvNz/dPQNt61kvcbUSpYwuktGQ/VVRGiBMO/jcM9IQzH5L4rYenCYjAfXgPTx4f8C8QMfruXmqI+2nuVRpOtxS6Cag1s9YsB9t+l+0DvtaAl0zFEcu0AJr39P1PAdQ1tXe1z26vmVzXJuQ1DlhMPT7Vm1r5U4QGaXxyhZc8m2UA1qLe1T7J+1PAKl62xzkK8IAATgpcH0epP5uj9AXwdIx2zPmq87Jfs6r0b9EZs1IRKh9LUX+bBTe/laWugeqWrZcMFHVESIV8yrGKx6zbn/0THI9gYmzoQa4dQKCAbdzCRw8mzrTCtaI0bd1hOW9bOvHSZ63VqpP632uUgM61pc+uvX/peN+HrYD7Ra2ecJSXdm3igjWd/HabNkChk7/kaVnPJqeAGtlTRztM4YQQTqg9HUZuoR1LoKhQFSURQ4MwBAInLWmMZ1CSnOm+IUzYdO+w735ET3dw6HGOUvN1jEdESF+xMRxxHj+gbt0T0pLaqiKKbApUL4sL6NfJHnKHap2uesN137B0HN4VRUM+Z6EzgJrbrbP/0jEicvWy3t2DwZjCAGKHSAGM99jGDQYMiBSwp59A7DBMj6jNnEtXEPa70evJwVR93aFOZmBokkzRPmN2jGVKot5qHvT6GDB9SYGbb8tWwP2GrPZklovdzJbECsaw9EDrpV7jlgu9kMdxYXuCdM91ECmUhwC2eLZDphLqsdTjnnm3F7jFuWytfq8kEtggodEZXzj/Sy0nY6S6v8HQPVM8gcAIqeuEJynTo1SDjC+VUlQPe+k6b343N8ZYgZulJNqg2LI+fLXX2wq4n2yvebLX/Fcr0Vk+7qX3IgjKTU5A8mQoBc+0/GLFgcYAj3NyTkrkv5Ue2RTQixXKMohJoC6EAczbJCELKeq+KSUUUZIBztSjpzt02OLgf4M3SoTJH5DTgdNSXuRgAb4pLem4S50nSlJHSzsUqiEgZmqjeIkalWe9X1mtUKfJlrtejsdV/QalH7pyTxEwpRKQ/nQC8wa9u0ff3SNyABg40yE/gBy6FGTs0ZP0YztPH/JKRHhqbSekT0mTqksOIgvzZl7YeVPojjodfInjZSx7u5a2QrX9dY1u+/otPOS/nq2A+8l2aQK2trTNp02uFhDrH6H13uwoakC5drzyQag/r4JMIlsiw6PacznewvEWDOmAq7paYvF2p+TVytJ+KuM3w9AMq1YdoOUMlxIWgBoss/dKgFZXq7Zt7lcF4o23PA8k1tRC20Y8RtE4I2pAT+rrBoyYIsDkEgkk52HqBOIqjjk9KFqNbQxX+NulkotLHvHS+yUQV45lQbc9xy0Vw2rACrhvZK8FzmsBBf18OQBSexLeHIkSv6fyqJAlRCIfErDpeJMkS5que2EKkBbLaZaetvPrzAqHKx1z9+jcDjv3Xj6NAVM4JN5W9h/9c+qRNpTC4LZFEJIagcqyXvqnabJFul6SbrbKx1I1zuLBBZOgUaqnaZPMcwZm8SRFvaHHC41qIprgnn434oV6MG8BWKlcSKuKMYPuFAYQjXCUVAuQehE973PSh46TeSP7c+N1E+cHSEuDLMcI9H5Yrfcl+ZfuZ9UwHyP1uqS4een+346tgPvmdmlpdMsTtnTDte1afm3pM6DmJEOWUM2lVtYjLR6ytcqDSsGavIzNLbgV5DgDVN/dg6nD2T+ByKVlu0OHDbSNj9aZLXehLNk1W60GNKN71doCZIJdVZAL1d9Wv1pdn+Fuc2UzUimbeqyNxjgfu/UObfnEsgrQuggxZaMBwMbdw9EWPo4AAVs8SKH19KDs3D533BCvd0IMBbxzIfIsDbPJDS0fK1en17Rst7zUpQyzS/ax4Pxt2wq4b2aXQPLSUuyaWU+j1T1e49B0Uwssde8r6aw7ZZBa0tnWmWpKIUhJyAigzmCSc+eWOjkRgtHzHgBwmn6H4y02fJ/UCx3GcMxBskwJRLNkjxM0RbgNetXZcPV1zMpQylUngNpAlQYy9gWASR60c3fJIxX+eO4Z1+Oqz297jpWEhVzqMgXRHG/hsEXACESgox4THGL00jmCtvBuSJ2H5Zw+SfCQ5G31VEgKBNqg7u5gg5pLUq+Wp23n16U5d0tzu9ZXaG0F3Dezt5pY7RLQSr2A+cQ3QLsYtS6AVHSlSjug5mNVNtaC8JJ6waogjE6WkydNxPBxBIOzJAoQmZSPU6YS5IqlgldrtVRLlvz53DkxY7lgjqVKovEGJe22PIRa3jlTEiiKCFu0p7yG2c/qZs29M9648uohVVVjdBj9c6YQmEY84l8CvoCAKuk97WXbKrjZfO8AtEBNXOptl+1SEkQLwMB8tXZNvdDaCrRLtgLuZ7e3nng16KnVwa5LgZH0r41szzxcbujkgCXRfNU9wAj9teALQQJkoA6O+uKhReFrfRhSPQfV6p4Xl/uV5QfDpWs0PG9OVFCgrANbZdy1l5896kVuuznWBVnWcnYeUJQMQARnja5UJPMINC3upXI4qZx2az4ZeiVrbet9bDeHcl2vsVuKhJfO+e9DmWBtBdw/hV2bpC/5MXjRD1CP2gtTz29eWKYKJlXJDcsA25oWR5HW21JfV0G1c/vs6cYofcty3dwwwIdz7nwbE0fZjs3Wa8jWeK0xnoyCwkieYoCmGjNtEHBGrqxlLi1UPcSa67PJFQCsp52DeVwCc5ZSsNxybFYSIQ4gSnVu8z3ZIFIPQhAqge+gyR+d28tqIQzw9gFxsfuu/Y6WP4tXKQJrrTd8LZPxY8Dz+wFatRVw/3L28fIb+aHWgFqKU18KtNwQ1uuymWr5lPYUk2OWug22C650iFgC3DI2BXU5ppVnTdW1yM+9XvbbTgh15bNrXqg5f0MlFLssr1L52pxvNvsmDj3ECRQKxcGa6kyMkLCNIyNSD6YekQK8BuVgpXmvsdeA3NJ1tvPjY4/9fdoKuF/VPtYrqIuMlI8MH6nBm3Bu9raRbPUOIcAG1yw3dWlqPR9XouN5G8PhsizlfRjSa6l+5cOQeEgJ5miLmmgkTdEkJSzl/xcQtdyuflZawdvAlY/nrJ4o26pnu0HxTm3kv7VWSgVYb1kSILRmblEuaEBPr0+uKy3zwxk+Tgjpeh1vpXsvgJ730KaVRNLHuO/uwKHDFA75YVWPr/U+S2p3ncBAeRzLWWW6v763pEpYmq/XpI4rEKutgPvV7WOXYAtyH1r2eiJi8TjzEvSyd7RY7KSpqVCDkN3X5aATgCT/QiqzODbe7ZwLzTIrsmOtEw9sa/C5Vnip9oLZp7ImWaPlpGcg3ABdlq21+5dx3LJc6zcyKDA8BkxhgBQYGvOxc8IIUlfkrLm2WtlLAbGF8y6C6EvVCS+1T93/27MVcN/MPheHdSmP3b5vioyTensoXKJ+BgkK1cWk7RF1ed5cS+v1RU37DSAuLcN1Ke14C8kiS4XNwzlzt6qWyPK0WdRdvcOunMdIq0SiNYKoB/Mup/1mHWpKWCjLeulmG8JoVBVdOs6QvOIiZcuXmIFcKqdR5rZNgE051XAqXr8JbBEueMwNMId4Tp0r5L+BHvMYFGhzC3vegLFJ752AeELpadan/cbqdMs98ZaUJ5dAd3YBF95fsu8vffearYD7ZvYlJtFyoKIoBmoKoOIPqzY01wN0hVJYSITInm7xPm0VririnzhbFe/nOrVJkzr3xtWjROUdliQG62HWrXXKdZc6vdWSvhr/MvjUSQx6Li5URWzTXdtt7bUvJ5HMTT3kCQFSEN4G4FxqjUTECEHkdIDw4no9WWnS1g6uOFitjWDrdyhl8Br7VI/3+7YVcP80dmkiX5ukXHlouWRfXm5akG0j1go8S5Hn+kdoPSTKU8ZDuy8Qlb5oQPKwc7S+dM61nq1G19uHhRygVkvEaLxT3ScmXW7rQZJVDpTj1TI2c5w03vK6pigq/tqC24LJPa6rtC0eRw9XJUVAerTFLnceBhK3Sz0mGuBpyHpcW4Rdi9SIZ9t+n9J1onjgdh68JOV2yRjXS4q+xL5Pb3cF3Dext8io+Ri+S6gByumtmHlaxZttxew2WGK9H6Be3jdL/pw04HJmmS7FmfrcTifEEeN0KJ0bKo+UZbmNdtyoQHF2N+IVMFwy6/Et8duaDJGvFZnGmB+rlnfNjo8SELNjXB6XetQpoGnuoY9nYCpBtOgDJhzRuX0qBLRBiEnhkLS5+vDUQFjh4H0aiTdAuwSyr8mEbBNzPta+L6BVWwH3T2UvmYTtD0R+XHXwyCyFr56rgG1J3701QgU7zmCrmWXWq5XWM0OmD+oGlTDA9fJATzGz3yVlQSycqwBh4lsrsLYp0PNqZ+11t/dzBvyLgUVA6x3Un9n99bVQEz5tw9xlyZjqcTnztFJRrXT3lftAzTn0yHNbg1lfw1bAfRP71MDCx52LIF0VpNyfBQRGqY2q/ajaH1jxdEUWdj1NsyzJS4BMjlK4W80gs6US21KQhU5QL7EGrVbIX1EBmSop44lRvDe93gL4etySGj33OlXRgHxNITYPrHydSDTGvMlide+bc8uqIUDrTWbagVz5RhKlEFEquFHKyHO8NU0zPSKErtltfsI4PWP0T+kgNmCn7cxtWvjy91vGuMq6voTdULWv9mn2ubwIkv8az+56hthcwqV7iV3+UV6zwiUW4A1hyoEylX/VnRRaDWk7/mabxmOb87GXIvB2oAsBrbTtUo3b66uDuSe8lD47G8Ls/ibv1PwP2r4HpWAO6UN1YZxx5jEvXSfj+lxcYeBL2erhfjb7GK/gJdWXgHY5bT2tfPboARhPMlsbLHPNazUrFVIvDiUwl5axPpzQuwfsNn/DOD0b+Zdpb548U8vhlnNYzla5x8RJZu1wUS8USZocR+RanD8vZRbn3mcrCyveqiRNxFgrHfK5852ovea26LfubxMN7LbajSPfg+SRVt+fKcrjY0AX9+keynGYevgw4Dj9gra4TklaabXJ1mqvtXSFmH9W8/6rl/sWtgLuZ7e3nqxzTezb2BLYNnYpBTeKNna5G4IcLzYeZgFW47G2+7b8bCVzSwBrHwifcC+u8bHzbgptBTdg7iXeGssFJYhpdaHJIiF2+doyT5721+4Qco/nXTxezo2vnO6XsBVwP8nemuO6fRwbgY45Fbd0Yl0WuQO1x5ooiYrvW9KnAvaHWDLWUnQ/Fa2JCDhPHyr9KNDJEeOEkKpiyWeJazWqADGbXVY8NFrwMPPfeVsNHqbjoF7SzwNeYXEb+Vc41lwzN4OYn21be4D2u7FbLtAITdWvog5JxIIpnDP6A0KYpCAQdzl1WoOWTB0chNoZ4yM08aO2VqHSms1Sq6+13v819hJH4/vjiFfA/SR7qYzmLc+ogLrEES5Jf655LvZHcd0TajlJ5g0og24Hph4hdSvIyoScjmo7RiT+syqWA1AFQnPQbMdi1Ri1zWkKtUslKZeDaZP5bIkjtrxofS12pHY88gCZroy9PjcAyaRLreApiEohVuBvVg7EQCpgXsZ5yVq1SyP/+2T78r+Nv4KtgPuXMquBpAVAes0EvuTRtsdL3k+jUOh4lxUBqg89Twlwwwm2hKBtAZO98awdFg/Vai8qSduF+hA6viqAZlOW7X7581Tc2ygQKsVEdeh5Km9thZ8FUIFoy6XWgcHy/wV0a2lcSc1lIJ7gU2Uxpg6OpQMylFZIh7J93C6XbGwfxHOJYW3fnwf6uY1ifEmJDeBis8HV3tAuLcPsDyECiUqQV3b7pR/apR/Ny7OKCo2gXhSjc3fpKKU4du5qEM4muKanty1fyjXM23Dr2IqXW29Tf9ZaxcVm5UQLqm2g8daDZ0krzLAlKOcJE63HaPfV7Vx1HZf2J1CqWdFJ3eGFmg8ASsv4yotuz6m2gulbW5tqvmSrHuRPY7cAUD0hqr06WOAh898ti83f7Ws5JxmAy+8lj6skOJzg/dPCD14PZ8GoeHcthVBbm8hhx1ZSdbOc6qopWDea36vnr++Jfejc/tksSbPs8S5TJfo/Ocf8umywzL6+fN52vC+dH5/Tvvb5v56tHu6f3uwPtpGFLXpmS57Ma3Pf7TldTuNVz1a6/54Q0RYVTxrSnNPvmuO+VkXQerlt2+56u/lV2GSDMp62PgKgXrAGEC+PpZV6ze2SZ3srcNXu03jxKauPeZcOJ52YHW/A1IOIMUx/oBR0lyptlz3d19BRqyzsJbZ6uN+cWe6v/RG1Xmprr/EqLgBYLlZTQNY2ZJShWHpgSQPcjvllY1jmq9v7cSsY1si+ZvVvL5+/HL9NuFg6d3veW2DVXsvCWFEnj2gxIAFhUS/YBp5VOvFVeykEfL9e6Vvaiz3c1VZbbbXVPs1WD3e11VZb7QvZCrirrbbaal/IVsBdbbXVVvtCtgLuaqutttoXshVwV1tttdW+kK2Au9pqq632hWwF3NVWW221L2Qr4K622mqrfSFbAXe11VZb7QvZ/wcIwbJMZ3OK4QAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(6,3))\n", "plt.imshow(output[64].cpu().T, cmap='magma', vmax=0.001)\n", "plt.axis('off')\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We can save the operator for PyTomography PSF modeling as follows:" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "psf_operator.save('/home/gpuvmadm/PointSpreadFunctionFitter/notebook_testing/psf_kernel_nearest_multi.pkl')" ] } ], "metadata": { "kernelspec": { "display_name": "pytomo_install_test", "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": 2 }