{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "\n", "# Photometric redshifts\n", "\n", "This example simulates galaxies with a simple photometric redshift model.\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Setup\n", "The simplest galaxies-only GLASS simulation, sampling galaxies using some\n", "redshift distribution. Then add a model for photometric redshifts with\n", "Gaussian errors.\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "import matplotlib.pyplot as plt\n", "import numpy as np\n", "\n", "# almost all GLASS functionality is available from the `glass` namespace\n", "import glass\n", "\n", "# how many arcmin2 over the entire sphere\n", "ARCMIN2_SPHERE = 60**6 // 100 / np.pi\n", "\n", "# creating a numpy random number generator for sampling\n", "rng = np.random.default_rng(seed=42)\n", "\n", "# galaxy density\n", "n_arcmin2 = 1e-4\n", "\n", "# photometric redshift error at redshift 0\n", "phz_sigma_0 = 0.05\n", "\n", "# parametric galaxy redshift distribution\n", "z = np.linspace(0, 3, 301)\n", "dndz = n_arcmin2 * glass.smail_nz(z, 1.0, 2.2, 1.5)\n", "\n", "# compute the over galaxy number density on the sphere\n", "ngal = np.trapezoid(dndz, z)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Simulation\n", "Simulate true and photometric redshifts.\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "# sample the number (not density) of galaxies from the Poisson distribution\n", "n = rng.poisson(ngal * ARCMIN2_SPHERE)\n", "\n", "# sample n true redshifts\n", "ztrue = glass.redshifts_from_nz(n, z, dndz, rng=rng, warn=False)\n", "\n", "# sample n photometric redshifts\n", "zphot = glass.gaussian_phz(ztrue, phz_sigma_0, rng=rng)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plots\n", "Make a couple of typical photometric redshift plots.\n", "\n", "First the $z$-vs-$z$ plot across the entire sample. The simple\n", "Gaussian error model only has the diagonal but no catastrophic outliers.\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAdAAAAHGCAYAAADE0tqGAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAnOBJREFUeJzt3QmwbXdZ5v8VwiSNgoDMJJCBzBMhkASVoEDEFCXV1S2mbQkKtNio0OnSJrStf6A0dCEi3VIJNgVppWkQuxK6mWIMBsSEIQNkZgiQgDI4A4JEL+dfn2We48pmn3vPPffec/Y553mqdu291/qt31rrDPu73/f3DvstLS0tDVVVVVVV7ZbusnvDq6qqqqqiArSqqqqq1qACtKqqqqrWoAK0qqqqqtagArSqqqqq1qACtKqqqqrWoAK0qqqqqtagArSqqqqq1qACtKqqqqrWoAK0qqqqqtaguw4LqPPOO298fPaznx3fH3XUUcMv//IvD0972tPmjr/ggguGn/zJn7zTtnvc4x7D3//936/6nN/61reGP/uzPxu+8zu/c9hvv/328A6qqqqqzSoVbr/61a8OD33oQ4e73OUumwugD3/4w4dXvOIVw6GHHjreyP/8n/9z+JEf+ZHhmmuuGWE6T9/1Xd81fPzjH19+v7sQBM9HPOIRe3ztVVVV1dbQ5z73uZFHmwqgT3/60+/0/ld/9VdHi/SDH/zgigAFzAc/+MFrPifLMz8wMK6qqqq2p77yla+MBlW4sKkAOtWOHTuGt73tbcPf/d3fDaeccsqK4772ta8NBx544OiKfcxjHjP82q/92oqwpW9+85vjI2KuE3gWoFVVVdV+u/BkLmwQ0XXXXTfc+973Htcyn//85w8XXnjhcOSRR84de9hhhw1veMMbhre//e3Dm970phGip5566vD5z39+xfnPPffc4T73uc/yo+7bqqqqane036L2A7399tuH2267bfjbv/3b4fd///eH17/+9cP73ve+FSE61T/8wz8MRxxxxHDmmWcOL3/5y1dlgcZkd75aoFVVVdtXX/nKV0bDalc8WFgX7t3vfvfhkEMOGV+feOKJw0c+8pHhNa95zfC6171ul8fe7W53G0444YThU5/61IpjWLYeVVVVVbUWLawLd1bcslOLcVfrplzAD3nIQ/b5dVVVVVXbUwtpgZ5zzjljzucBBxwwBve8+c1vHi677LLh4osvHvc/61nPGh72sIeN65j0spe9bDj55JNHi/Vv/uZvhle+8pXDrbfeOjz3uc/d4DupqqqqtqoWEqBf/vKXR0h+4QtfGP3Qxx577AjPpzzlKeN+a6PT5Na//uu/Hp73vOcNX/ziF4fv/u7vHl2+l19++arWS6uqqqpqSwURLeqicVVVVbW1tVoebJo10KqqqqpaJBWgVVVVVbUGFaBVVVVVtQYVoFVVVVW1BhWgVVVVVbUGFaBVVVXVwmnHjh1jHQDPi6oCtKqqqloo7dixY6w8B6CeFxWiBWhVVVW1UPr6178+VpUjz94vogrQqqqqaqF0r3vda7jvfe87vvbs/SKqAK2qqqoWRjt27Bgtzrve9a7Dd37nd45ds/bff/9hEVWAVlVVVQu39vmP//iPo+W5O/Bc78CjArSqqqraFGufO3YCyNtvv334i7/4i/G49Qo8KkCrqqqqhV/73LFCZK5n0ATPP/3TPx1B+pd/+ZfrEni0kO3Mqqqqqu2n/ffff1zznLf2CYjgSQBp3wMe8IDR1asFJt373vceO6joF70egUcFaFVVVbUw2n///UeAzirroeD5Hd/xHcM3vvGNsRc0i9N7zx7gqS/0egQeFaBVVVXVprFO999//xGengHz7ne/++jGfeADHzg+g+96Re12DbSqqqracO1YRQQtMD7oQQ8a7n//+4/jwFPDa9bpesOTCtCqqqpq05Tu23///cccUfBkgVoDtc0a6XqX/StAq6qqqk1Vuu+bd4ASOL/whS+Mj9UeuzfVNdCqqqpqQ3WvOwKEAFD6irXOVCSaV0xhOv4hD3nIuI31mmPXSwVoVVVVtVDpKxSX7jR4aN74abrKepf9qwu3qqqqWpj0lf3vWM/cHZfu9Nj1VAFaVVVVbapuLDsWpF9oAVpVVVUtvEv3q5MUl0XpF1qAVlVVVQtdkeibM9bmovQLLUCrqqqqhdXX51ibK9XMXe92Zo3CraqqqjZVisvUQg00FVf4yle+MhZXmBe5uy9UC7SqqqpaOO24A4y0UoeWBBOB65e+9KXhtttuW9d2ZgVoVVVVtVDaMRNlS/PSVOLeVVz+z//8z8fyftqZeW47s6qqqmrb6euT3p8ACZ47a3EGuPe73/2Gr33ta2Nrs7Yzq6qqqralVlr3pNkSf/ZljMe80n/7SgVoVVVVtVDafydRtvNK/MU65bpdT3UNtKqqqlo47T+nPN+iFFCICtCqqqpqU+heC1JAISpAq6qqqk3t2t0oFaBVVVXVXtOOVVQD2pOKQRvVeWWeCtCqqqpqr2jHnC4ps7DcG51U1rtk30oqQKuqqqq9AqmvzwT5BJJTWK4lEGh6LYvSyoyaxlJVVVXdSSuli+xu/ubtt99+J1hyve4sx3M112LOQNccmW898z+jWqBVVVXVnbS7VuKOFerWep5GzdpmLoXfV6ptO2v1ev/5z39+3JZrypxZC90oa7QAraqq2mbalXt2d9JFduykbq0HWNJ+++03/OM//uM4zvOsxTjPNTu1Nr1OSb/AeQrujcgLXUiAnnfeecOxxx47fNd3fdf4OOWUU4Z3v/vdOz3mbW9723D44YcP97znPYdjjjlmeNe73rVu11tVVbVZtJo1xN1JF/n6TqxVc4MlaTWmS8p0XKzKrJfOrp/+9V//9fDFL35xGaDGB8yB6UbmhS4kQB/+8IcPr3jFK4arrrpquPLKK4cf+IEfGH7kR35kuOGGG+aOv/zyy4czzzxzeM5znjNcc801wzOe8Yzxcf3116/7tVdVVW0F9+xq00XutRNrdXouAEypvak7F1QD1pTky3yOuc997jMWiVcsflpQPgFJK7mD10P7LS0tLQ2bQH54r3zlK0dIzuqZz3zm8Hd/93fDO97xjuVtJ5988nD88ccP559//qrm9+3IL0orHFZvVVXVVtS0h2ZAtqfguf3224e/+Iu/GB7wgAfcqR7t7LnAznuQBb8vfOELIxzp3ve+9/DABz5wOSCIjAVXc+qw4jlu26x97q17WAsPFj4K1w+LexYguXLn6YorrhjOPvvsO207/fTTh4suumjFef1i4q/PD6yqqmqra+qeBbRpZ5O1aMfETes5LtbZcwVyAaxzPuQhDxm+/OUvj+/Bc3odWQONqzZty3xuOyZz7azd2bZ04dJ11103fiPxQ3/+858/XHjhhcORRx45dywf+YMe9KA7bfM+vvN5Ovfcc8dvGHk84hGP2Ov3UFVVtYiKlZegnj2JYP36LlzCO3MFOycL8v73v/+3wdM1mStWp7n/7M/+bLRI9fykWK9ZH11vLSxADzvssOGjH/3o8KEPfWj4mZ/5meGss84abrzxxr02/znnnDOa53l87nOf22tzV1VVLbr2VmeTe62hwPsUkGA4r+sKUAIka5N72PtvfOMbyxYs6DKyjJkGJa2nFtaF64d0yCGHjK9PPPHE4SMf+cjwmte8Znjd6173bWMf/OAHD1/60pfutM1721dSmq9WVVVtR+1uQYOVtP8qInZnm2B7PU0/mXXBmgdY//zP/3xc+/yrv/qr0crMXLyGBKixQm3nEl5PV+7CWqCz+ta3vnWnNcuprI1eeumld9p2ySWXrLhmWlVVtd21Hp1Ndtzheg0w4yreldVqHDiyMOWPCiLNcpssDUt0rhsws93r9U5jWUgLlHv1aU972nDAAQeMP/Q3v/nNw2WXXTZcfPHF4/5nPetZw8Me9rBxHZNe+MIXDk984hOHV73qVcMZZ5wxvOUtbxnTX377t397g++kqqpqcZX1yZ1Ziytt3zEpxsBjmMjYRNqCsjXWlQJ+dhbINGsdZwwl/cUYD67c6THDdgeoHzhICnH2zUJRBfB8ylOeMu6/7bbbhrvc5Z+N51NPPXWE7C/90i8NL3nJS4ZDDz10jMA9+uijN/AuqqqqtkYN3MDTZzJLz/Zv3hENa32Sm5XL9dZbb112rxrP7Wq7ZxYl6zGu4gQyzTvfStG7864vlmyuf0+jirdkHui+VvNAq6ra7gKmrE1SrEWWoCAen48+J8Hsa1/72rg2mW1ACVqe7QfMwDRl92aDhVY63+5e397Ob10tDzbNGmhVVVW1b7XaiFrb73Of+4zrk4HVNLgHeEA0btbMOwu03Y3gXW3Vo/Wqi7uQLtyqqqpqddqbbsuVAoumhQxS8GDHHSkjYMmtm8pAcbV6dl07Kwe4UlEHmndP0+L0nqf79lZU8e6oAK2qqtpmfTt3N7Ao65WJdM065YPuiIadWoJZl1wtyHLdCsdz+7qXlPqbvafdrXq0r1WAVlVVbVLtKp9yb1q0K4H1Oyfb8jrbVwMy5wNPwaHWTwUlWXf0OpWKMu+u7nfeNe5LdQ20qqpqk2otVYD2tNVZxv3NHW3IVmqHFisy+1fqQZpoXS7fP/3TPx3+xb/4F+PxQGrbtEzf3r7fPVUt0Kqqqk2qve22XI1FGwv1L+9oQRZYrtQcO+uWqbuba875vJbL6RiWp8YhwGmc9VVijeb+1ttNuzMVoFVVVZtYe9NtGYsu+Z7zLLzZNmQ75pTQmwWxMdY3XasmH/JBVRmaAhZEUyTBNueYV4Bhvd20O1MBWlVVVe1SsTxBbl4bsqmmEbHJAZ3mjQaetk/zQz24ZhPtO9vzc6r1LpowTwVoVVXVNtNK8InlGMBRADgb7Xv/nZTQyxj7kpupOMFnP/vZ4W53u9vw+c9/ftzPukzu6GwkcUBKs+7afRF9vBY1iKiqqmoLaKUgnV0FCgnS8TrBOqDHWsw647QPJ3nvuPve975ziyNMldZj1jDvec97DgcffPA4p1KsWevUcSXW7GwBhJV6iW5E0YR5qgVaVVW1yTSvsPtqLbLp+iQoJtIV5BK0E5gmWChVhqxjJqBnV8p5XMe0rJ/trNeUA/ye7/meEZIruWrnaSOKJsxTAVpVVbWJNA+Wu5MPOoVPatam32bWKvVSFsRjXxTIZtt3f/d3L1uY81zCuS7HJfXEOTX7MM71gqp5psesxhW7KNG4BWhVVdUm0jxYTmHFuhPVmoo+K61Ppnye99yoUkjMp1UkSLIMBQvFdZuuKmApX5PAbzZnNG5gaSvJ7/RszPRaksc5LRe4O4FM5tvIACIqQKuqqvaR9kWk6LxemVNYKfA+m3M5D6Jxm7I4k0qiNJ/X4Dm1DHPeVAyK+/X2228fPve5z43PAo8yF8vSNuIWBvV55f1252cz7bYSS3gjA4ioAK2qqtoH2leRorPuy+laY8riRbFQUxVotlC744GPC/d+97vfWAXINec8iYaNWKkPeMADhi996UvL50oBeZG1hxxyyOj+nVqq5s06aM5Ju/rZzH75yH2ah8XsuvdW+cK1qgCtqqraJHVq55XKC3ymFqntn/70p4eDDjpo2d057aMZcKXRNRh5jvXqtWOm65xknzG2ffGLXxz+/u//fhzn8YhHPGLZ8pTDKbLWg7WYoghcwokUTgAT8Bo/2xh7FrDTgClQdo0bGUBEBWhVVdU+0L6MFJ1XKi9rmlmTlCaStclYlUBGsd5cI/dqAO9YViPosfKIxRnrMNA297+4o2at/M7DDz98dPuayxzm47oFxltvvXWcA4ydPwBMwFLK9nmf86z05SOdXjy7V7DeSIA2D7SqqmofaF9Gis7mQSaghoUIWAoWyLsEQ2NTlAA4QS55nrm2FEtIH0/HJfgnOZZxqXLjyuf8zjtcw8cee+wIz4Awbt8DDjhg3AauAMrinLp2XQMXcPJXp7mcqykan36jOyt6v69VC7SqqmofaXfrtq426GiedRurjXUXOB1xxBEj8CgpK7H2Zq8x8ExRd67a1KedWrzOB4j3uMM1nJZj6c5CXLfmBOxYs9OCCCJ9bQPZpLlMLcl5Xz6m15CI42m5wI1YBy1Aq6qqNkhTYNJqg47mAWZaiF2xduuRQAiKWXcEHQFDs6X6pnMGmq4lcxs7z6X69TvWYHMvOYdjlpaWRlcvuRZQ95w5XZfXeQSS0y8QKxWozznMuVLR+/VQAVpVVbUBml3HnLpLV2NVTQFjrmkqy9QFOi3YPg3AmbptKeDKnNPgoVmL91vf+tZY3o8VyX0bsElpcZz1UZG6trM2HevZOmjuGcSdM5Zp4DntBDObRjO9hqyhbmRx+QK0qqpqAzQbKBPoTd2ysxbqSoCYTWUBqbhjvQ6osm4pAMc+sHJcrMl5vTqzvjotvgCUn/rUp0YIisblpjXfQx/60PFYzbCncDYuwUh5zAZXuQ7rpKDIBZzc0XlWd34u064usbTXMze0QURVVVUboNlAmWlAz2yuJFgEkvOCZlYKupkGGwUusSxByjrpzTffPKa8mNPYnMNr+6WZxDIOaFmwD37wg5fL/gkqYmECqchbliYJZHK+FJ0P+BI1vLugi9XtMb1f2oji8rVAq6qqFihKNy7U6brjbPrJbHGEleaaLfE3DdQBQfODnLVK8DvwwAPHsc6R/p3ctSDo+LhbU9bvUY961Fi9yHsP50q0r3HcuNY957VEm7Wkp67c9BhdyS07e79pjbbexeUL0KqqqgWM0p2u+VkTpETBsuDmBRxN55pdF03dWzIOLJN7mb6cxgRG8jjTgSWKFRz38AMf+MARnNN7iTuVO/bRj370cklA95HCCvPWeAPVrH+m0lEChmat1dn73Yji8nXhVlVVLWD/zqmVBSgJ8AE5a4zT/pzzXJZx307XRadzm0tA0T/8wz8sA5PFSM4LZInmBcpcR9yxD3/4w+fCP9eddJfANmufgBhAzt77FIrgKV/UePe6K7dsjl3Puri1QKuqqha0Tu6slQU8WQNNQfVZ1+x0LCWq1ZjZDi2OBahEs8atOy+9ZVo/dyVLb8fE5TovBSVF6APDWMZJmUnAUtzO4A3s1ldzPxsRbbuSCtCqqqp1FpjszJ252ibVq4043VkPUaASRDRdt4zlmBq0s8dO02e+egdwUwlp3peCwG7awDtu3mlD7+k9gbv9SX+ZLaaw0Z1YqACtqqpaR826MxMwsxpN10Vjea4EkAASGMGaQCrF2wNBFh5gEmDlWrxOhGsgGbdx3Lnf/OY3x7mBL/VtU64vkcWUer1x56Z4feA5r8NK3MVTS3NfFuhfiwrQqqqqdRJ4pBVYSutNi6ivRllDXMmFOrUIA1vATCuwwJu8t9aoqwrX6l3ucpcRstZYAXEaCWusfYCbog9/9Vd/Nb5Pb1CWoteOnXUrx/U8LSloPKsVmGc7rMwLsJpXwnAjXboFaFVV1Too7sfp+qVAnNWmXEyLIkwLDEQABnyAZm4VgliWsSiBGyCtK1oXZcGCl5J7Ulm8D8zBLMUWEsnrkTliPd/9jgbeCd4B4liXWdsM2FaynlN8fras3zzNpq/QRrp0C9Cqqqp10M7WL1djRc2rXBR5zZL8xCc+MUbVajGmWtBjHvOYcZ/xQGnupKcADwgCrbJ8YAh+LMjkdAae0d3udrdxnMdDHvKQEZKOMWfaok0t3Fmw7SwAaWcpPSuNW6lG73qpAK2qqloHrWSB7SowJnCdNs5OvmXWJB0LYmD28Y9/fASfQvLWFVmTlHMB7BSKxFKdppUAsW0p0RdrE3iB2nnudcf9JNjIObl0jQHwWLsrFX7Iva3Fgpz3M9mI5toFaFVV1TpoJQtsZ4ExAYztQAWKgafgHdbiNChJEXdl9eRPWl8EOuMTvBO37W233TZWEOJyZQ0DMYsyrtqU4ourOUUWuH5dg3ki53WNauKC82c+85k75WMedNBByykq0wbgcWfvTgH9nTUT34iI3AK0qqpqnZS1wKklNi8wZrYYAqixJgncwCMVgBLck+ActWmnxQuMBUuADLCSrkLK+BnDPWtupfsSLBSLMp1PWJZJvaHALoXib7rppuERj3jEeA5zxaoN8CipK+kQ45EKS6uxIBcpErcAraqqWqNWs3a5mp6f8yxTxyW9g2UHSF7HSkxEq/EsQudIjVpuW2MDSdtZjkkvyTaWaqxZco1x21obDdxSTk+3layh3muyBpt7cA6wZuXG4g0U8yUhoM89BICrtSB39oVjvVWAVlVVrUGrWb8DiWk915V6fiYwJpGuQCSiVnAP2LDwAhxgzFolGAFeLNnkZn7yk58cn1mMdNVVV42F342bluVL+gngUVqExTWca/NwHudm4aawQZRuMfaZa1rgPvun7dDs40ZOqsy09+iutKtgpPVUAVpVVbUG7cqVmGLorDwWGysvKRvzXJaBLfjZJ5AHyByb1BXu0RQeoLhuFUeI5RjYgZNIXOuc0lQcw5JN+ozjXC+L0nguXGMJRO2fFmBIesw83WNO15d5AVJxX5sn6S4pXr87Wm3E7r5WAVpVVbUG7cqVOK3nai0x1X5yLEts6v5leQKesaBrv6AgUHvkIx85rmEK4gn8YrmBXNy6CeThRs22wCmu2Gk+amDL4iUAj4VqTKJtA1L34tqA7xGPeMTyWuZKXV/mfclYjQt2WhBivQvE744K0KqqqjVoXlL/tFj7tP4ra1DBgQT1zKsbaxuQGQ+ArMZpmojtHtyn5ooFx0o1D3etAKAUhE+7MfMAI8uQCzYpK+AWK07krPfGgnBcw1zIrgXoXC942g7y9773vce5ppC0HUitla4UIDXPBTu7Tux1UnSm8yyaCtCqqqo1Kq7EabpJIlADvqxtTgNr0rR61v0LRI4Fm0S/EsuUS9c+ka6JtAWWQw89dHThAmiCjcDHvK6JxZqcTj1A0xx7GtyjYLtCC7feeus4VmQu6AKu185hO1gHqn/1V391p8LyrOMEPu2q0ffUBRt4pmvMNEUn+21fBJftpugHeu655w4nnXTS+APzDeoZz3jGmBy8M11wwQXjt7Lpw7epqqqqfa1YTMDISgOfpHAkFQUQgMNnEwUInkFH0JAx3qczCpACHssycDnggANGizaF1p3TeBA8/PDDx7XK9P5Mab9UDuJ2nbpNuY2TVsL1CpKuz3qoa06Ki3Hm9HnMSk5Q1P6TiGHXbryxgJ5gqcByJQsynWkc7znA9AXCY9qKbaXeqRulhbRA3/e+9w0veMELRoj61vaSl7xkeOpTnzrceOON45rASvKLnYI2f6hVVVW7q90pUu4DHlRYZeAWC9OD9QhQoJR1x0S4+owCG/CI9Wi/eVh9KeYOeo73zHoETMBLIXbHSDsx9pBDDhmrAF1zzTVjUQXzpOvJ1CIGeeA2L8g6fpo/ah4WKWiyXCmpMsB897vf/U6F690365TmFZNfreKyTbGGVDlalBZmCw/Q97znPd9mXfolCsX+/u///hWP88fItVFVVbUn2p0ScwmisT8WHEixpoAGoADLl38WVaw34ANOUbfTqFcu16mSt8mAOPLII5cDfJI6kmpB3K1gaby5WaquJakpaWgNvIkOTgF410jGBdiONfe0OXYAC7bfuKNzSqoKgb5zpKvKatNTYp1mzTZfWKat1BalcMKmAOisfDui1HRcSf4YfVPi7lBE+dd+7deGo446au5Yv/C4OSj5UlVVVbtT7SZjA5dsA09WJVglijbBOADnc43LE3TBx5KTtUXwANqcPxYY+cyyPwFAYGPtEPh8hjl/AJQats6TDjBZd/U+0cGgBYIsvsAw1uNsP07bA83bb799ebxrcK8s30B9tT1OY3H6+cyz9hepcMKmWAOdCgxf9KIXDU94whOGo48+esVxhx122PCGN7xhePvb3z686U1vGo879dRTR/fISuus8bF7WBuoqqqiNHOeQmRemkWAZUzSRkAFmEixgMA3BQ/SAzTwBJ4jjjhi/AwDWkCbnheIWYss2FtuuWX49Kc/Pe5PwBJwOo9zc7nmuhgUOZfnWIzOCZaAz0rltZtakrMVk5LyMoWdgCMCdkBzza7PlwAWtOW03XGz7myddJEKJ8xqv6X8JBZUP/MzPzO8+93vHj7wgQ98m2tjZ8of5Zlnnjm8/OUvX5UFCqK+rSUnqqqq7auV1kBno0YTyCP9gwsWTEDPmOREqgKU5tEACGbG86pxoZ544ol3cmVOz2vuQBLsplWIzAPWwOlcQGjdEyR9viV/c9rrM8UMaHrORBHny0DGZNy8FmKUKORE4Zov1vasNrL59e7Iz5thtSseLLQL92d/9meHd7zjHcP73//+3YIn8emfcMIJw6c+9am5+9O4taqqap7mpVr44AePrF8CVOrQpmSfdBPwSLAOKw88A6jkcwInT5no2RR4D4ycI9ZWurBYd8z5sy3XBtqsVOd1XR4MggTipMQfj1yCnKYdVeb1GjXXPLdpLNQvTL5A5OcVYAfKs188NrL59b7QQgKUUfxzP/dzw4UXXjhcdtll47e33ZVf1nXXXTf88A//8D65xqqqtofm5SlO9wECiw/UFCQQh8EDZl2SdRgIp1KPL/fgBqKCfliRWc/0WiQr0IFL8klTds81kLnBiNUJwIAJjCxfn59g7j23aiJYwdD+RAoH2DurDrRat+n+dwAza63TKkKL2EVlSwNUCsub3/zmcT3TD5irg5jUaeHzrGc9a/zDsZZJL3vZy4aTTz55DL32y3nlK185ujWe+9znbui9VFW1uTVrcQIMiy9l+jyDFlj5TAK9gAro0oEkwTzWCVmfiTi1Peun5mfBkmdromAbOBkLzs5LIBqrNW5XUb32gzlQ5jMz5877KShXWmecB7iA8N53RPZmXPqOrhTss8jBQFsKoOedd974fNppp91p+xvf+Mbh2c9+9vjaeoCIr8g3ruc973kjbP0BW1O4/PLLx7DvqqqqvaWp+5HFGGgAiBJ2yfGUVmI/1yu4pQSe4CHjABnkfI6xOM2Rkn7Tknuxco1xPAtWbmkaZqc0oPmBycPnYVJqElBkLAF3gDuF2GoLtE9B+J13wDvu2p0F+yxyMNCWDSJatEXjqqq2tws3aSWs0uRgJgUkVYdSWSgtvpKCBx4pMu+h8AvYpbqP9VKABRfWZzqn2Od8cdEal1zTuJWBNO3JQN3xrjPrksDq3MBqjtn6tLv7M/nqHWuZ084y84KtFj1gaMsGEVVVVW20ZvMU48KMBeo5KSAgx0pMLVpg4RUDEyBNz07pKEkzYU2CHbeuBxABsg/vLF+xGgPKVCEyvw/6gDGWKxkLoranZm3aoE0DfFLFKPmnq4Veyvh9dYU1za0YMDRPBWhVVVtGa7F6VnPM1L05dWGmKHssVNBLhSAAASyWjLHOY2z6doIdeHKzOtZ21iiLMwUYAM7xYGgtNBBVlYilysJMSp5rEYmbTirTdmW2J/DI+xSqjwWdgKTdKZl3r52sac6L6l3kriprVQFaVdW2LL/nQ35366xOj5uu5QFD8j8BicvWfuub1iwBMRZhUlY8Utv7+OOPHy1ExeCltiRYKKBNOzJjpOZZD7Vf4XZrrNyyce/GpZrAJdcTa9UxcTdTXL6xRKPVRMnuv5M1zdl1UtqK1ujCVyKqqqpajbL+OK3SsyvQxsrb2TFZ7wOaHAe6s9ZUgAikgCadRbAQiHo4Ph1ZzJl1Ua5dx3kPtvYbZx7nsc15ZRWkW0kaclsLTdF5lYxS1N18U1et6FuPBF6mbi+r1LlYommjNg0IWmsFof0ncJ02x97Z72UzqhZoVVVbQqtNk5i6FxNwkwLqs8dMYZvgn6llFetq2gCalQdI5gUneexJVUmZPmPBMSXv7Ac+a562szBvvvnmcax93Lz6gHLvWq90DikzaW9mvkTfkvM7n+tO5G7ct67LOeNiTmPuFJdJ55YUyF+rtbj/pDDEztJbNrMK0KqqtoRWmyaR3pLABEjep4D67DGBbXIwk385BUFSQjwDILAJBuJS5Y7NtaX9F6AZa1/yQLNO6DkuW/BlwbJkFVewn7UKply9zs8ijVU5BXiChezL2qt5zA2YU7BR3LquOYUiQBzwpyX51rLGvP8WTF+JCtCqqraMdpXLGIsya39xV650TMrUARXY5H1AkL6aKR5vPVJwD3CBVRpSJ61E0BAAcaemqwqrMO3BwJFVev31148u1qx5mscYa56pIpQuVQKE3Af42pbOLOYBU8exQFM2MPBkYcYCt29apSiu5qkbeE8ia/dfZY7pZlPXQKuq2jZKtCxAAAIw7mxNblqiLkFAU7CAjopBKVYAWMAFhsDKIgVSQPvsZz+7bMUCom1g570x056YgoqcC/xE37ICXau1TRAHYucgYwDwoIMOGqGbDixJhwGurImCY1qtAW6+BLBwc32A7NpiMefns9o15u2kWqBVVW0bJZ8zxQVSTm9n7smUqEukaly3cQM7JpV+gCoFFQJJr50TTLlgsyYZKzXrj8AuipaVxwUsqCcQJpC0Dup6zJti8NZAU5jBceY3lpXrfOZ3n9OiB0lZcU85jtwHIDtu6qrOlwjnnhZOmKf1LqCwkQUbCtCqqra0Zj9gPU+7lMy6J3NMoks9uFbj8p22AEuEKWABYlJLUic3MPXM+gOdwNMcXgOcwCCQvvLKK0cgmo9rNy7YlOZLWzL7wDvrlmBpzdI4EM61xaKduk+nrtcUgzC3a/EFwFppCkakvynYTpt4r+S+3bHOBRQ2umBDAVpV1ZbTrvI8Z3tVToOFWKdJWWEFGs/iy/qidBIWZQKD0maMFRe3KHhau3QesLQdnMAKID2r551iCFyonj1sd7zrEOQE3o5xDgADuhtuuGG5OTZ4mp/VCJ7GOWfWShMoNLsO6Z5dR9JrQNl9ZUwiZ/0cEn3sZ7RS4+uN6Liy0R1eCtCqqraUVpt6Mh2f3ErNroHEeGuBxALjFrWdO9X7lPVj3QFrmmqDmffgF+sWWFmOrEOBRGA1DdwBPsE95DwB4RFHHDHehzSYAMtx5nae1M8VqZuKQ+YHY/uSrrLaoggpRD/bNHteqs9KbtN7rXPHlY3u8NIgoqqqtpSmwS758CcfsFnvnAdbz+m9mWCflLlLTmdcutYOAxj7WHysxLhYYymCpGuJS9YzqzYNtb3O9XoEiln31J7RdaSTinVSbl2A5Sp2ndqXBebOAbwe6bqSWr2xsCPXaIy5Z928aahN7sl1Jg2Fpj+zHZM51ztlZaNTZGqBVlW1pTRbq5YLlPUHIivleQJloAUIqT+bdJHAL2udcfVmPBcvkDmX48AR4FIBCHwBFpScz7U8+tGPHq1Frlavbb/22muXIe08XLopdpDKRc7HbWy8MfZ77Zpz77EYPScnNM8BzbQE4RSCK4EpX0Sm1unfzLFu1ztlZSNTZGqBVlW1pRSLiwANhECQu5R1N1UsLXAAqORzpqTdtDIPVykYgmiKJ4BWSu9lbudmfYKe83p98MEHL7ttwTapJ9y14AnQ1kMpEbTeAyyrM8USlO0T5ON457zmmmvG/dZNr7vuutEiTCF7z7lO+1MTdzb9JOk5sSanpffmrXdOrdP7zrHqt5NqgVZVtaWUIgGUIgcgA6CU1l1T2IIL2KXmbYJ+plV5Uhg+ZfvABeTSvgxIQMi5zWmbYwDRuKSwgLKcS8BOrVxrr86RggryPLlzWcUgBYDTlJsEK6WAgjVZEE4wUCobuUb3ZY3VfSWFJU24p0E4qWQUKzk/y1mIbrTbdJFUgFZVtaU0hQJgcbXKtwzAphaT14Ft3KYpwQckqSaU3EyWofVH+0HOeJamtcp0RjEGIK0bghcLjys5LlkWqIf3omhT2CHAAkPrnK7VfucHPWAk0HKNV1111Qgxxwp4SplB75Myk3q2qXO7UsuyNAd3vwm6SuF61z8Lyq1aWWh3VYBWVbVl10Cn64asx0Am1haQJCoWtGwHO3MAW9yfABJYHnXUUaNlC2KsWZYqaKcLCkilSLzX6RkaFy6oGpcyfcaZz9xxF4M9oAE3UFk/pVi/IA3aruHEE09cLlhvjrhtQdqzbQFm5piuX+bnlbVWx4C3n1Ws2T1JD9mxgYUO9rUK0KqqtuwaaKJHk5uZAuwAkbqv6VwCekropSZt8iONA1dzyctMugrImc/rpLEkyhZsQdJ75wAm12BN03Gx+myPK5VlmOjZNNKeFkSIVZiUGpaqMc5nW9ZcjQsIA+Q0zo7lOK0mZJsvAmANlikgn1J+e5IesmODCx3saxWgVVVtGSUAJh/S3K0+tMGGwMoYIIuFCi6JvE1JPOO4RblpAcRr44zPuihYCQCiVAECTecMNONGTUSuNc+kn3gPnOAH3pQm2iCXdmPk3MalcEMsS9awcZ/5zGfGceAX8KXDivOlC0vWPmfXL6frmtmf+9wT6H19gwsd7GsVoFVVbSlwstKSdpFgmryO5Zf+meDD1algO/glEpa7F7Qy5rjjjhthlPN87GMfG6GX7irGxy2ckn/OCaYg6Jk1ax3VuqZjnJcCWS3KbHf96YgSgJovbcjsB0oP1wKOzgXgvgQkQjZWo2PA11wpADEPiLPrmntjnfNeG1zoYF+rAK2qalMrbsIpOG0DjLhcfYDbpkk1K8+aoZJ8nuPupFigwMj9ypVrP6sPAADQmmgidlP4AFwTJes4MMu6Ynp/ukbwJG5XVnAsP+Azt/N5OA7A4/o1d/qWUnI4PYxzTazhuK4DPvedgvG5v/XO0bzHFo7YbR5oVVULq1h8cSem9dY0nzOVfBKlah/3pfeAErdh1gBTnCCWaWCb/d6zSBN09NjHPna5qlDcumB9yy23jGkmicyl9OI03jjHg9+0qwoQK7knCMhaKQh77XrAUwRvAoDizk3kbazc5HO6fscDfAA7Wx2Isi66EQDb/w5LdqvBk2qBVlW1kJoNQAEXLlYgCjjS5DpuQrACwkSaJo0j1YJScYgVartI26yLinSNBQdexjqHc7I49f0Eb9YekMqtdB7bvHasa0igUiJyFThwjDk/+MEPju5g1+1+0sHFOYxnFceCjbUJfFkrNca1cPHmfARQs9G1eUz7jG7nogf7QgVoVVULqdkAFGuCyZcEJNGsKaQ+dRM6DiRZfwAZi1Kk6Y033jhClbs0li0LkgWXJtWA5TjAAV0wTO4k6zCVgliZzgdqCfIBtqyjmlP/TqkqzqdYAtctKKsmBIDeJwqYFWt7CuC7NpAFa89A7Z4cn8pHWSflqp5G106LJaQH6lZMI9loFaBVVS2kZgNQwAacRJxaS/Q+yprgtBen4BrPIJPenMZkDdO4pH2AU1I5ACkVg1iHAAeAjk2gD6jF7ZuOLInAdV2HHXbY6OJ1raD68Y9/fJzbdXAvg6/117iWHZPIWl8Mkj5jXud2Hbalb2c6w4C7fRSIziuWsJUiXxdJBWhVVQuplQJQAIlFlsbRcVECW9ybCcKxbhhLFnjAyFiQ4UYFLlBzHAsR5IDL/KDH0mQVSlMBWNfhNSvUXOZmrQIbIE4L0HPPOo/encBv3TMuX9apcenMAoazvUpdl/MAZAraJ7o3ViiLOAUYpjmeWzl1ZJHUIKKqqhZWswEoLCsgSVrKdN0PALlc08cTsFKUgFUXwCVwJ27NVA0yBgzJMUBnLfKP//iPR6CCpPOnQIH3rNnkfxrPDUyifYHY9YGn84Bpisi7XuONA3kA9Qj8wBG8HXPQQQeNj1wrd3TWdr02BmwDyRZ7Xz/VAq2qal21J6XdZt26wAZaKaoOPkAChmAKToALcumiwnJl0aV4AeAmwjYNsL02BzdsolsFA4EVcIIwyKoc5BpYpc4LenJE417mUgZxAj3nBj0WcurPJk/T9bg+x7CWwdzPKaUFze0+XB9gm9f1CmByvnlFEbZi6sgiqQCtqmrdtFJpt9VCdVqmL11OBAZx26bJdNYTk1/JsgPDVBgCJ9Bj5YGW6+EmTR9Pr1mnoMZtC0SONW/cwOrPJqApzbABLQFKaY+mNGCCeY4++ugRvMCY2rrGZr3VORKZa5uoYEFS5L5STzcBVK4tFubsz6zF3tdHBWhVVeumeaXd0gYsLlmgSvPmnbUq4wIFFttAElBBKjmZQGO78SzEWJ6sUlZkXKdgZC7jUvcWgAAwnUhSRAH8vGbdsgxzXqBLcI81U3MDoXvJPNy37tcY9+x+WZSZ28O85gJPEE8xeQ/XzzK2L27nQnJjVYBWVbVumrpgU281634sMJGzBHDzLNMpgI0HLUqUbooFAOX1118/WnTmSEswz6DlWOdIagx4cpGCmzm95hq1j8s1a6kJXAI2gOZKjVWsEhBQgrH3sWDN53j3az9AJt3Ewz2xogHXvK4xjbqTowrWKSqf/Ndq41WAVlW1borbNoEt4JEgIdBLQXfbjQvsAqlUIAISlp61xHRAScF4ENYrU4qHY4BJiod8TOcFL9utnTrGOQJY8wIcKxTMANjrdGexzzNLNO5h0HUN5nXdzpPKRfZz49LUdZ3jQdw1CXDKvYGj68raq7GBKSA7h/dd29x4FaBVVa2r8uEfS9IzWLD4UoYPzIAtjaLBEFhSW5alaQzApBhB1hbNlxxKkE29WvPKyzQHCzORruZLi7C0F3N99nERgyAIsyzjkjUnQNrOBcw963qSk5koYRam60nULTet64gl7NzG5NxJr3Efzu1aXaf7z7it2hpsM6oArapq3bVSl464clPMAFTSxsv4VOgJIL0GtBRGAEjAswaZaj3m52oFYwUOWIdpgG1O++J6ZanaZxvYcdsCJRcroLoe66yu21okeHq2dmmb6+LmNRfIGu941q753KdrdX+Ok3rDvZtG31kf9cXAOPfA+s3PJpG7e5LfuZUbXK+3CtCqqtZFsx/cSbVgdcVlC2QgEWAGkAQcSTUBFhanY1lngASWLLtEs4JgiiakATUL0thE3nKJEtesbdZCQdF1AWCKwDsfSIKg3E0QZiFOIZw1WMUNnvjEJ47X4mG/GrdZ702fUA+vzRELOaUIATLNsF2/LwFpxWbbWuvabvUG1+utArSqqn2ulT64s84Z12zyILlGU3CdgAbEUmc2BRNs46YFL1Ze2ntxqRpjDmMAByxj1bISwRKcPLLuCVwABZ7eg1qqF7l2lmL6hrIsRfea27ZpWzMP53VdxjmnnFJWrEeKKaTlWdy2WesEUWu57sO5nNuYqRt7LeDb6g2u11sFaFVVG5a+Ap6iWVOMHYg8p6j6wQcfvNzdBDBUF7IfXBL0A5TGghLrT9WeuIRBywMMwcz5AAz4pkFCyf8ETvAGY5aoIB/XbZw1WjKPOQNarxNNC27q5QbEgJfygKkOFMsWHO1L/dxpIQTv3VsgHAvUtXAzr7Ux9VZvcL3eaim/qqr2eR/PlJeL5RRXZSDDmgQTYGNxxeWZyNRYdcnPZDFyyVrTTEuzuEspFp41UddxzTXXjO3IwBZkjTWHdUYwTkcVwHRs1mGBjpzLedMVJnmZ5nJ+92IsqKeIvOOB3/zm9N786cKSMoQs7tkep7mnvDZ22hs07u3pz3g1apWivataoFVVrYu7NhGmQBH3qffcm3IoQS0QZcVl7S/9PWMtOR7MWJJyPY2zJul8oOh9Oqmkn2eebUsVIM8sS7DL+iOgJrf06quvHreDq2Afz+YBWu5g7lzXD47kHq2XAp/57XdO92ftMgXup6k6gJ5AKc/Jf6XkjcZaTNTxzn7Gq1ELMOw9FaBVVe1VrbTOlm4poMDCDEgIENNIGoC8DmwBNW5eEDOn/eZgvXG12gZ8om9TaOB7v/d7R2vQcVzB4MkiJMBhPXLdAtuVV145HHPMMSMYnY/VeNRRR40gdJ50Ugnk0hEl0bppbQaEqWPrWFaqcwHrIYccsgxUD0pJQfefkn5x9U6tReOmObEB8OzPuFpfLaQL99xzzx1OOumk8Q+C2+YZz3jG2FpoV3rb2942fnP1zdQ/w7ve9a51ud6qqv5JCQTKh3lqtU7bjiWqNaknwAFm9nFtAl2iZ1mY2oEpTqA4ggbV1ic/+clPjmBjyV577bUjsBzHFWouc4PsCSecMDzlKU8ZU0bSKgzwACk1cV2P/VlzBCSWp/lcA0j6TMk65ZFHHjkWRzj22GNHmBprH5AH4s7N6kw6inVL86UikXMl9zNRxtPI4yhBRQma8vMyb362059xtf5aSIC+733vG17wghcMH/zgB4dLLrlk/MN76lOfuhyBN0+XX375cOaZZw7Pec5zxvUO0PXwD1hV1d7X7Bpc3IqxjHyox7WY1AxWGtgBV16nYTU4sLISmRu3baoJBZBgMs23BCQWovf2AQwrMm5XwTop0QdmYOf4WJpJPWGdArNt7svnjXsCvbhOfUE3l4fr9QyOyv2BKbhyR3vt2X7HWENlgYJsSviRuQE4AE01o6myVpy+pcYl/Yc8dy1zY7TfUiovL7D8QftDA9bv//7vnzvmmc985vgH/453vGN528knnzx+Uzz//PN3eQ7/ePkDTZujqqrmK7CcRnNOXbc0dSsazzLjigUAr7lLk9PJCkuBBI80xfY/HTeqjyqfA87p3B/96EfH+ROlmobYYMxiVB6PZQlm/qedmxUH7KzYrB86h88Kr7l8ndc5A15rpCkCD4aHHnroeKxo39yfawLgdIPxxcD1Zh73Joo3a5z5OeX42Z/lFIj5WSc/Npbwzo6p9kyr5cGmWAN1ExRXxzxdccUVw9lnn32nbaeffvpw0UUXzR3vjy/fKvMDq6pqz9JSVkqRmKZmJF3FNkBkAcaychz3p/lZVkCkdRjLDZxi/YEVqw+IuXHNw6o1hqWaykA+/IDHZweLE1y5bNPo2nOKEjjOOOdkrbL6nMe1+awQoOTLvA/WAC3WsfeszHzoSsNxPKg7p3t3Xsp68Ly1znkgzP5p67IU4CfXYVu+kFTb3IU7lX+mF73oRcMTnvCEsYTWSvLHmaoikff5o523zuqPPQ//jFVVrU5JS6HpB/tKIEgVoqSHxPK0tgl2tqcQQdy53LPgpJG16j/2ARjLLmt+KWsHogKFwDCWIjdwcji5aJ0/QUbgmUbcZL/rsA88uVXzJds2D1YtaxJUzau8n/XXBCalylLSbcyX63TewBWgKe3bppGxKwFwdv80Lcj9xkrdnZSWas+18BaotVDrmB/4wAf26rznnHPOnSxWf+yFaFWtTivBcl6KROAZa47FmWAj/3f+vxMsA54JrgHRFHwHTpaWMQCUijxx8ZrXs/3mM96xIMmq9PBl2pxgnDVU5yPn4+51jb54m8s1uRbgBlTbzc8NK7CJxclqTYuzFMFPv9L0/ATUrFMak5zOpM7syc8/lYymLeIajbt+WmiA/uzP/uy4pvn+979/XKjfmeLimSoh7vOUWpRVVa1Nq80nBAxrmB5x74JR1ivBzzqjtUWgScoGt6fXXK/mAEOgYwWKxgU3AAXhWJ6sR1BLnVrj0/IM+FxDWo8letYzi9UxsQiNMw94mxOMuZLJNdkO5Ao5mCsdX5JeYrz5gD5VhyLwtH1Pg3/itp1dC622uQvXNzrwvPDCC4f3vve9YzDArnTKKacMl1566Z22ieC1vaqq9YvKBZFpdG6KxAMKCw+AWHgAwhoEK5AEUZajZ9GwqflqfCJyY70CkoCgtPzigv3IRz4yBhYBHyuQlZnrybpmrD/eJiD1LIKXi9W1JeLWs32+IJgPqK1nOj4uZl/qvXauaZuxuLMFDTmn5xSESApKXK57qlYW2ljddVHdtm9+85uHt7/97eMfRtYxrSH4Q6VnPetZ4x+mtUx64QtfOHZAeNWrXjWcccYZw1ve8pYxOfq3f/u3N/Reqmo7dFehWEL+R+NWTFNo76e9MZMX6T14skJTrcdxLDrpaKxG5wEp0EwgofP6Yp1SenGjAqqxUkjAmPUnoAfoRNj6DEnPzkTHWod1HYGlz5VUKjKXa+IWNsa9sJRdVyxU5xfhC/SzLuoEI2WtMmvAe1LPdlatLLRxWkiAnnfeeePzaaeddqftb3zjG4dnP/vZ42t/sP7Io1NPPXWE7i/90i8NL3nJS8Y/chG4Ows8qqpqbZotJRerE8BYawmkSREFkAEQH/T+b0GTm9O+WIzgBrLgB1jmsM7IXQtUieQ11rkA0LxcuykwkNxMIHM8sE67njgHoKaMn30+SwDZtQE0gMaaFi2cewRgS0Kg7N7NIzI3pQkFPU1LD4K2+/c+aTquyfO0nm17c25eLSRAV5Oaetlll33btn/9r//1+Kiqan3TWOKmBYtYpHlOZxKAAB/rdrFMAU30LC9TAnZS69bx4Jbo29SZtV/ADxcs6KWhNRduqhlxv4IbUDsvuKVgA6WuresF0lyPZ/AGbG7dI444Yjw/2CuEAP7O5dxgaI4Ug7fd8axO12Vt1c/El4Dp/Hujnm21GFpIgFZVtdiazfnMOiQYpCKQR/pkgox1SOCRlgJuAKTQiTlYd1kbBUvgAaIE66TTCreuMp1Ke5oX1MANmFmy5kjhBPN6zUJN6b4Uimc9upaA0PW7Xufggk2wkXsyh23p5Zm80+ScpnSh145hpebLw3RNNMUPplZ769lubhWgVVWtuYQfOCSaNBGr5HXad1E6pHCzgllalSl9xyLlcgURsANdD8eYnzWZtUVWJHCZy7NjWJyCkXJOQLWEk/qxYMjKTZS+a3UdAWFK4hkTaPsSAOTOA+Cu0Xj7WM3Ax4Xsi4Dt5kzZQGN9AdiV1Z5czkbQbl4VoFVVfZvAjQsznUeirNklmhSAkuPpGC5SVhsLEBTAC1hALUUNEvDDauSWBUjgYV06FliAEqzM4XjQcgzYpOUY8AGk+bhJzc8KjCWn0IHrMg8Xq2O4ZIE1BeRjPTseoPPsnljKxBXr+t2Tn0e6wqS8X9Y2Uyh+JTfsSpWaGkG7eVWAVlV1JwFhIJO2XYEosLDIAIUCD+Ah47gwgRSAYwl6HTdlihjEyrSfS5b1CY6eAZM719zyL0EwKSNAlCIIAo4cA5SZ07lE4KcYPNjH0k2ZQPOmApJoXmB0X+7VuZODmmhaqTU6uwAfaLM8p6X48jrW+bygoJVSTuq23bxayDzQqqo2Tin6Tp69j4AhNV89gCgRsIDjdYDAGgzcUjQhtWvBSq6n8YAJyI5PoI2xtnmw8gQV2c+SBB7WpGsRJJRzOW8Ci7KG6hgwVPidpcu6ND/QZd3TNTlfUk6SdgLOjgVlADanazfOsdPSebGO46Zdqaxe1kSTj1ptbtUCraptrHlpFHHbgic357QbRQDArRnrDCC5Y+P+ZGECGvCIaGWxsdCyrkhAZ35zgZP3zpMAJNZj0kNE1yZaFuhYg4ofpH8niDo3aKWcX1ylAAjA7sc++ZeuJfVsvWf9Og9Iu177vbZmympN/qbzpT6v174EOK/zObc5nde5VgoKmhd5S01l2ZwqQKtqm2qlNIq4bbli476dBW2KIgjeAbAUMQAe+wEJjIAVUJPiIcfS2iIQggsrFGiBLBDNmiqrNmXxjGXNpn1ZgonsM6/x5nHNwJyG3KxLEbfJH3Uv4G/NFcxZv0cdddR4zSoNpV6tuVmagM56dZzzp+pQ1ktjTfpZrKZQwrxAovQBbSrL5lMBWlXbVPNaksVi8qFunW+296fxQJdAGMBKuTzuXEFDsdZADURZtM4FZF4Dk3FJOTFWGgu4AlN6aToPGBnv/AAZiAFzCjE4nzVSxybYJ8UUHAOIrExQY9mCKQsWRBNRC7axqK2fgpi10biQA1AP15U81fx8/Nz8HGYLJewqkMi1ez3vd1AtvgrQqtqmmhcVOht9G8sylpWyd5TWX6ywuCxZcKBpn7EsQWATWatYO2DZJ7AILNNODOS4eQELnFinKdAOJmBpHGgF3s4JfI5NUJB95nUfzgPuyUNNkfkENiUyGExZoo4FWu5h9w6oLGeBQ4kANibpOtNG2Kli5Oe5KzfsbCBRClA0lWVzqgCtqm2q2Q/ztBdL9C3IJOXEAxDB0HPSNUAG2EAGELNWSYJwuHhZewngAdH08AQiEPPaeaSxRMYAGmsza48eCUaKxZkgIdcNfM7vXNyu5rU/cLvqqqvG9VYWLrdzomXda6xdVrcvBamG5jxp1k0Aa3yaWqcQQmrprsb9Olu7tsXgtzFA//f//t/DmWee+W3bFXP/sR/7sT2dvqqqfajphzlLM40bQBKgPFJHlhUHUiDCFcpiBFBKeTogYr0BJYs0FlaAwxo0N4szwT/Js2RpOk/q5YKYeViO5gEr18NSM8a1ACxwUwozUEoKpl42+AJx7i9F4LmSzWm+rMN6dj9JzXFc6tamlGA6t+yNQggtBr+NAfrTP/3TcwH67//9vy9Aq2oTaTb6FmQSDWv9EkBAkcXJFeqRll8EnKCSdUqWIsvUvCw4+8BGyzFjABU0jREwlC4mQO4YlYFSP9c1mB9EXZfrSYH55JgCMYFZ2pyRMdy0jrcmC46uOxG1KXKftV6v06kFTLPW6XXyYd1LCuXn3ms9bj+tGaCXX375+Owb2xVXXHGnAvDcPNw/VVVtHiX6NlZf1j0pLsyUuTMObLwGS+Mdzz067cOZzwE1a82V5tYACmpcpgBpDtsSEJQUGJYpOAG5MV7bnzq3qW5ke1y7WZ81b/pyWvcU/XvyySePoHbtLGqwjes6UbvWaxMUNJu+Y3x6k5JrynPq3a5F7cqyzQD64z/+4+Ozf4h/82/+zfJ2f8xyp/7bf/tve+cKq6paNyVFhaUZEKWDCbiBln2gxOLzoX/YYYctW4pAAqqBTGrhJr8TXI1NlR+QsoYJvqxQ8GYdppUYUBoHeoBrrHGAzQoUKASGjrfW6jwg6lpdu/HGmTfWZoKCwMpcrhGg03ZM+zTbuHATIRsr2/HO5z7d494oBt+uLNsQoNZB6HnPe97wP/7H/9ib11RV1V7Urqyb6f6U6ktd20TkWq8ECdvUk00loLh9wYdLlNLiy3Pae1mnjJUKVoCWwJvkfurMknVPoLU+aR9wJlAJaGL5xsIUrOQ9uEXmBVTzADJ3MNi6Nz2CARk8bWN5Otb1feITnxitZdeZ806LHcQqTS3g5MPuyRroztKJqi2+BgqeXC/vec97xn8gLpkf+qEfGv/Yq6raWM1aN2DhfWA63U/JrUwAEZCwHJNK4jhAjPXJwiN5nKCSdmYEtqy+tBrzOnmUzpPoV0CdRrumPu3HPvax5UII9mcu121uRRgIcMk8rtnnEDA6BgxdF7jlfgI8D/eR964lJQH9TIz3xSCBUPmZmQ9w4962fU+sxpWKzFeLrz0GqLXQpz/96aMbxx+6AAQBRP/v//2/4dRTT907V1lV1ZoEDiJGfeD74E9uZ1yFsX5sY2UCFsuLVZfuIo5NVCzwARPAsObkhaYQQvpbAoH9J5100ni+BOYkj9J6pOPA0LYA3XlT8g88WX+u3TnTNcU1JYDJnDxh3Km+uIvAdY8pZGC747Nem0IHIO4Lvi8T9rtG926f+3Yfro0ANkCbulddSyKLV5u+spJWKjJfbQOAvvCFLxzOO++84Ud/9EeXt73tbW8bfv7nf36s6FFV1cYooEzz5pSai7WTD2zQS9F36SXAyBIM3CjNoFlo3LPeSwnxWtCNCFVro+YFLoC64YYbRusQmKwrgpz3goNYd6AJxrFqWZWsx0T/pnB9UlGSDgOsKeoey9a1XXvtteP5gREQU2t3WhkIjF0DACaiNtZwfk7m9mybn8PU/e0hxiOVmfaWxdhUls2p/Zam4bNrEBeOb4PTb03+YK2N5FvcZhA3lX96HyTT4tlVtVnlAz4WFiCwplLEwAc/0IAFq83/L+B88IMfHP9vU+oOOATzAIb/DWNSnxZozZcm1KxKr/0f2cd9SmDmGNuB0XnT9cT8wOGa0pkl9WZ9NEmnAWf7FUdIEXfWqrm4al2n8yfwxzWmjF/K+7kf66DGutd8mbAv67yzy07OFYs9oMznXKNmt7ZWy4M9tkCPOeaY4fzzzx9e8IIXLG973eteNxZorqpq463P9Owk8JpanlkbZZF5z+IClfTs9OwYIPHMQrQ/vTyTRuLZWEBhSYI16AAiINsPVM6RfE0AE+DjGNZuzuVaAA7wgJGA0nUDdNJZjMmzOUAyHV0sJ7FMc0+5VkBlwaa4PMUinX6pSLnAWYs9VmItxmqvABQ8rYG+6lWvGv8ZfPvzB/t//+//7U+4qjZIgWOKoGcd1PvAEywohQi4YNOqC8xYeY4FLfvjqrSPG9P/eQJ64gaO5QdirFbWp/H2c7fGLQyIQBy3qcjeG2+8cbQqs+4KoMBmTnBzneZ2XvO7Vsd7TnEHIOUqds3OC4asWBB3/YFn+o7GsswjXzryJSMwbXBPtU8AauFevcsPfehDy1G4j3/848c/4KqqNkZxLQIngSHrC6zyDKYgw0q0npngGf/D6a0JctPiBrEybQNDIAtAHZOqP2CZakSpMwvM5rHPdQG4dVI5pdPyfLmGrMVyMacoQ/JGuYeBTUBSiifkWnz2sEQFNhLXciJpXVus8hRHSDPsrIvOszQb3FPts2Ly+ZbpH4YFmm720wILVVXte6U4AaW7ii+2gOL/MiXqgNE44JBv6X84BQhinQFZOo/4/wY1ENMaLOktXMOAmPVN1qJjQTkFC1JDN2uK3oOnACQAcw3md52pViQf0xjnzHFgmipCzilvlCWbbijGBYIszWyL+3e69hm3rHO4Dz8bXxhcJ/dv1zardQHoK17xiuGlL33pcOyxxy7XjCR/7AVoVa2fEtiSJH8AADSWGUACGvAkgpTlBx6OYSUCWlyZqTaUNVHHmTevATJgBWfzOy5WquOSvpK+n8R6BGXnMs5r1x3wppVZon59KY9FyErljiXPgCfAKBHBrgFgATqQJF8aAlXn9/NIM2v3KZ/UeV0DeDouNXBbIajapwD9jd/4jTFdpUFDVbWxShWhKfBS0g48ffiDCwvSs5QVY0EHzFhnARh4sOYSgZtgfeuMlGAeMs8UxMCW9JLUm/Vsn2cws9/1mCOu5BRoN5d9zpniDwDsIQDJttwPuQdwzHv7EqiUUoD2T3ucpuqS67BdHitr17U6PvtbIajapwD1zye8vKqqjbE4AYa7NXVfs24JYFlztGYIMBo9JPkfbEEqwT2pvEPAxuJilaVlGRdrqvnYDzbcth7OCTy2pzYtcCfgx7yuCciMz3pl1h4JXFmW3//93z8CLQDPfbhugYrmsZbpHlL0HgRZmlnrJPcH+InkNT75nlkHjWuZAQCagqNsy/poKwRVex2gqXlJv/iLvzg8+9nPHs4555zlUPnIN7qqqvYdPK1fAgT3pSo6PuTB0MPrpGekZmvK61ljBLtU1ElbrmmHkkTtOofjACiu0UTRcp3GOmU5xiXr/OYHJOdIFGws1ETA5rzep/l1SvAlCClRwMAGnMZxsYIsS9N4kDeXL/NpDg7ofi5AHus2mpblE+1rjPNN1z5bIajaJwD1x5tvltPqQ1Ol20JVVXtfYOjDX/ALkMVtCkrgBlpxbWadEBy8ZlUqOMASTOBO3LtpLh1L0lhgNA9YOQ8oAqnjs6YKjLmugNuYWJ+pg+t8irkn8ChdWpwv7clYjq7FMaxq1xZLEQhB1TGxZKe1e1MUwjWkwpD9ibqNMsbPa2cBQ833rPY6QPPPUlXVxiqgCVyyFgiQAEYpoA4EtgeYSt8lj5OrE8QE0aQTCsAJ9vHsHKxAZftYgQAGdMn/9IU51iegpkVZirA7fwrHO06dXLBzHLcymDneeV0HGLoWAm7uWV8MgDTBRyDsulnKXL2BZK6B3PN0vXSqwrFaiDSWqqrWV/ngT5EBEGRd2g5qgWcCYbwHFvuBSeRprEPAAjXQS51aAEq3FVAzN4EdcHHdsvA8EvWa4gZpsm2elLtLu7BUGjKnsWBojrRKS+GEHOf89rF6E4yUHNB0dvEMrK7J+6S+OD/Num+raqEA+kd/9EfDW9/61vGf0DdcheV/4Ad+YG9MXVXVHMWiSlBP1ibzUKgAjEDJ+p5iAsCT4ugicb0Gn0DQa+BMVxVzmyepKcaY2/+5eeP29ZpF6HWCfsDSHLmeRMWyQEERDEHW2mhK/4EfazHrpEDsGs2ZOryxbL0HY6BMKT73l/qlgJ/1y7XkdLbWbbUuAP31X//1MZXlOc95zvCYxzxmzA/7iZ/4ieFFL3rR8Au/8At7On1VVTv5cAceJfBUBEtRdlYkiAAT6CTyNF1GgBHIVAEK0BLNmpZlqQmbou5ctralwljSUlLOzzjPYAdoLEyWr6jaRMWaj8sVpIEza6dZb82aaK6HO9p8QOjc7ilu3qyLutZpoYT0M407d56LdjUNxpv/Wa1bHigLNGWz6N/+2387nHbaaQVoVe1lTT/cU3pPYA6gKabufQoe2M+tCaSglH6ctgNtrETPwARSxqeIQYBkG1A6B+hMcyg9ANJ5ldXjamWtgqFjAFrqSVJVyLUAP2+V6zAORFmQ+nuCMcjHEk1Aomsyz/d+7/eOY1jV3LWu3bW6hqTzsJznAXQ1cGz+Z7VuAPVPlfWRyB91vnVWVbXnitWUZ/roRz86pqNw0X74wx8erc1AhqXHavQaGEEtgTjpnBKLML1CjQMjwGX5AWWiWs1lm+NZhoKL4jZ2PQlUYhWaw9wBtOPS49M4QM0XbmNAV2CS69fdybEqBD3qUY8avxSwlBMkZT7jkr4CsEnV8XMA3eTDzq57msN1xzpdCY7N/6xWq3/qVLsHYmX+q3/1r4Y//uM/Hr89vv/97x9+7Md+bMwP9c04j6qq1qZZqyn5mGkFZtkEHFPdJ2uJ3Lrg5Euu/alzK1fyaU972nhsKgYlLURVIc8p4m6/86Z2rn0J4MkaJLdr1mHtAznzxUULhrFoHZNqQ8Dk+uKqZZGmHCDXr+sWtJS1VccoL+jhS7qfA9ibw3UGsvnyPtuj2DUkvzU9PucFFzX/s1o3C/Tnf/7nx+eLL774Ttvf/e53Dz/3cz83vm5OaFWtTvPW5+JStE9hAABgdSU/M5ad4uqia2PVAaOAH+BIsYUE5xgPssDkC27q0JozlmO2gd60lZn9tseFmx6clP6gaZjNmnR9rivbfdG2PaX9EoTkekCb1ele033FPQI7wMYK9hwlbcdz0nVmcz9ne5/uam2zKS7VXrVA/RGfeeaZ4x/7VKl6srNH4VlVu2dp+sBPHdesO3JzBkxcmCDCnWkJhbsT6BwPaOnKEiuRlQoYxnOJsiidI7VoA0CQS6RurilAsR0wAdG83rP4sp6Z3pzkepOrSXHtOgYofZ5kHZVFzO1qjPlAFbzch3tLlaGs1ebLg/OxrKNUIMqXkFicuR4WJ7leXzBqWVbrZoH6Z06qyjve8Y7lHKtIUvbll18+PP/5z9/ji6qq7aippSkIJm28wMo+AE0vTmuX1vOMs804wT4Ca8CRixXMfOEFPK+Te+lx0003jfMBq2Ozjpk+oEAnICdF3eP6jEWa1BMQBzrnIJ8LQAdygSpYATVIsi5zP4Gk+4jLmHs2Vi2wsjx9CU+OqGtynWCbZtcgm6jhVEAyZrrGWbdsteEu3P/8n//z8MY3vnF4ylOeMrpo842OrrvuuuEFL3hBAVpVa1TctlytSQ1JRC2AWRcEzlh1AGhMihsAHXA6nluWdQY+rDVAAhTvWXQAlnKcyQMFznwxjvs2+Z6JtDUnWIFh1lfjJnYO5/cMhM7B3Zoi71mzZH2m5Zkx5lDcwf25X2u7zpUm2clHZX2aO5b2FIhTi52ANLmgUd2y1YYGER1xxBFjsJB/GIUS0tqoqqo9V6wkbsw0ZphGs4MDaLH6QAh4UmjAeKABQ9Zg3LveSy8J4EA3XUiMt42lCMQsSs/gCHBgmFQZMHMuFjI4ObdnFmyK1BtjLtfGTZviDElhAThWYure+vxIykuqIaXsXtzO6ewCxJ4FD7lG55wWf/dzy9pw2pnF0u0SUrUwUbjWWkDUP6ZcT26jqqr2XFm7m1bPARhQ8GCFxdpTLAEgQQvkbr755nE7gLG8zMOiAxbHAKX/2bguWXYAlzqyKYqQaF3XYJvXyR+NNZcoW1av6wPYlA9MIYXU1gW8uI+TM2ofOCZa2HbAc05jFYZgVZvPPqDXaUb+py8K0xzNrNUmfcfYpNc4v7HTtdCq2vAoXC4W6SpPfvKTx959l1566V69qKrarvBMOUwAZaVRLDKWYYKAAGgKMQDi/nQcy4sLFWD8r8Y6TJuxFFkg58yaJ+s0pfymaR8p0+dcsXpZuYmijYVpnLmBy/kCuwQgaT+WNV3gtATkfhN1a4yloIy3DYADePM5xnXN5mjO5m4mpYZaDKFauDxQLhUQ9Ufpm6Ew9b0l8z796U8f/3n8I1100UU7HX/ZZZctRw9OH7WOq0VSImPnWUO2szAByHPcnQEdSIEncAm6MY5VByb+zlls4CNFROBNrDTWYZpee836cyxwOc68IOqRIvAJKkqtW0rOJPgBbqJswQ3Y088zgT6+BLg2luMpp5wyvk+rNJaxACFz6cV57LHHjsB3r9zN5mUZa3um4AKL1heBaTBQgqmi2SAhz4nRWCnfs6rWzQJVfJobaCprHcr4/dAP/dDwspe9bPwH2hvyz37ccccNP/VTPzX8y3/5L4e1XuNsg++q2iitVEIuUE0AjEfWDQHHexYmqAFXSu2BpzHcnSldB0h59r/oC2jOGauQJRdYZlyKIHikNZr9GRPgA1ncpmm8najcuGdj0bpm5w+8AjTABU1zuS7vPexzbSzrVBnyv5wApuSUxgWdKFvP+VnOBgk16rZaGID6ZjhP/sgvueSSsf7t3nLlqpLisbuKi6eqFk3z6qumpizXa9JEwDMwSKk9kBNV6xivsyZpPnOYD4SyDulZrAJLNO5P8yaX0pi4OUE61wGYKQofi812X5QD9hR7T9s0APN/B6JpiWabc7oGj6SysKBBFdC4bc2XtBtzsEop0OTmpdnSe5SfZUoNzsvrbNRttSnamfknvfDCC4eNlkosvnFz/fx//9//NzzhCU9YcaxxabpLvp1X1b7SvPqqWfMEkPwNplBJrFEgYVWCnZJ4nsE0x5gjbtXUn2WZcumCI+AZb17Q8d7/KxCDjusCu6xfpnBD8kVTNtB1xMp0vO2xaIEQANOPNNAybppDKlUlEM8aZ9zHSaVhmeZLRO5/JWsTPBN5m7SWWprVempLNNS2vnL++ecPj33sY8d/pNe//vVjhPCHPvShscXaPJ177rnDS1/60nW/1mp7al4iP1D52wVISs/MqUUKNKxQoAU0Lk4CrlQQAlDroikebwnEOViDtgGU8eT8qZkbS9T5KC5cXiXwA8QUpk+RBe+tRybVBficy3Gu3xxg7HzmsSZruzHmcZx5FMJnIRvrnsCTa5cSCRwr27Gzbtz8DFOWb2eWaFXtK+23lBpeCyofIKzbZzzjGbt13BOf+MQxWOF3f/d3V22B+mBIInhVrYeyBkqp4cpiTCsxYBCgd/XVV4/wYEmCZf5O/Q0DsP8TQLM/+ZHmM4d9ntN8W5APqNmWerMpuxkLkqvVPFnf5Gp1rgQBca8CX9yrWa8UXOhY2x2b7YKJnBswE/zDqnRuXyLsTyF4cHSuRCRPO9Dk5+T+sp5snliigXUhWu2J8CDN2XfGgy1hgc7T4x73uOEDH/jAivv9k7VNUbXR8kGfdfuU6PMAtgTzACDw+GfOmmkgk4LuWSsFUMck/9Fr24xP6kkstkAzhd1ZkaDqXM4NhiDIdezDBKASZORDBazSuSVNul2fACfrpo5PZSFgZNE6ly8AIGuOaSF5YDRuNko5pQrj/k4louTMUizRpqxU66ktC1AuIt9eq2rRu61kOwhyeQIi0IAE1ybgABxAgah9sfDShzMgBF4gBiRgYm3alr6gKQMYKxD4ADD74hJOHqVz+jaeNdGsXabgvfcChZwj1q/XclhBNJWScl1gylq0z72AsPNM3bMJuHJtsc5jUWftONvTfWUK2H4xrrY1QH0L5jKKfKgAon867qNzzjlnXD/5nd/5nXH/b/7mb47/xEcdddT4T2wN9L3vfe/wB3/wBxt4F1X1zwp0sq6XbcCQggWAAjiACWaJaE2pPA9juGpBggWXovPgA0yUUneglNZlsdDizk1PTMcZk/XPFEaw3f8d8AVqCTxyXa7RdqDkauWCdSwg2g/orlH5T3Mlgti5wd25AdK+bGOtJu1lZw2tV9rflJVqvbWQAL3yyiuHJz3pScvvzz777PH5rLPOGi644ILxH1ZIfOSf/z/+x/84QtU/l8TsP/zDP7zTHFW10TmgCRby95pqQEDmtS+H4AOQQMAaVXEIEMCIJenvO9GpSWehdC8h1mMAGCszbcZsB2Ggyzpp8j9Zq3HdOn/WE9NLNGklgOncChzIQY0L1vVJTUnwki+0jnFvXru/BALlflOr1ry+FNgWAE4Drmg2h3ZejmfdttV6a+GDiBZt0biqdkeAxFpLWgp4+PLHcvKv58sgmKUwgXJ6qS0LmLH2WFvJmUxbM1BJ+zCuVmuOqWiUAgcJDKJpAXbnBh7wcy4gj5vUnFnrTB4pUFkS4QEyd5p7C3hy7WlrltQxEbXGK/wAtmCcvNUUeE8uaaxjEkg0a0U6LlY7dY2z2tfa9kFEVbXRmtaTBSjgI7BJQFBcpVlf9E8LOunFya3pn9nxxiYQKIUPvAe8NJYGZNuNSWEEj5QGJAAMJGNBZo2TlWhO1mACjVJBKPVonduSivmU4EvUsPVNFiqgeu8eUkghwEyAkGfBRmnbRrM9PKNduXSraqNUgFbVXg4IihLskpSUwENUq2NAKQUEuGxZjKADrqAJUL79pqtIonIT1JNC79P+nSkqQIm4jUuWEpGba8k8KRZvblA0HtDBLZBmDWc8axfoUwvXszGpRhTr1r2kOH5A6tnPwLyAG4imYIOfQ6zVFE1oWb5qEVWAVtVerG2bfakQlGR/cGGpfexjHxuh4rgTTjhhhEsKInDRerYeCH6Oyz5KdK55E4kbwABjIJ21Uc9xG6fZdYCb4CHgnNbBJRBngSaVxrH2Jf8TOAFP4JA5nZuL1j2yOEHYa2uiKVyf6FtiSSYAyX7Ado4UgpgWSFipzm1VLYIK0KraS7Vt46K0z/ogCAGBbaADoFdcccVyEff3ve99Y/lJ40FT+b3kVcb6Ax+Qc7x9yfNMdG0qCnkGzJQCpAQLmQvsnTdl/1JEwX7u36SxOIabVlUva6nyNdOkGsy4odMlyTUDJrdq0mpYn8b62RjrHriozctqzZcBrwPPVCcyh3uxRgzAzemsFl0FaFWtQSuty3lvbQ/wuGVZWWDAZQkMrLhPfOITI/yASIAQgOoklGCZrFFOy/RZj0yeJVAmGjaF2+N+jZuUvAbJ1JFNkYKsjYKa+Z0rxevjijVONS/1a53fucEN0AHStZvjkEMOGS1nljSYsnADWV8Y3Dfr088izbFjYbqeuHDzcN588eh6Z7XoKkCrag2aty4HMKy2wC5rglrzgUjgKFLVdpABT25Q4AKzdEEBkKw/xl0c6CRiliXJeuNyZQ3Opq2kJq3rcF7H0rTSj7nsB1NzJa0kVYXcn9fmc2+JmjVejicoct9y16o0ZD/r0rWIxA20fXHgpgVvQJ8tuZefZ6CZR9c7q0VWAVpVa9R0XS6VhFhjLDkg8R64jAENVmnyPbkub7rppvHYBA1xaQJecjyNjzsXeAAlwAIw1pv3KemXKj0Uyw3I4h4FaHCMZQpusaTB3DmMBWOws15rHBC7h6xvgiE4Aqr7MadrB3LnAUH3S+ZlkYKu66WVgoGmP8/AvqoWWQVoVe0FsS6B0Ad/Su+lBydQpmMJiy8F04EHpFh43LupNZvqPOBmm+dYhaw7QAWpjE2t3HRGyXpo1koTkJNiCuYCU2uf3hsLyMBpn9cpdiBdBQQT9KMSkeOAFpwvv/zy8cuA10CZCF6Qd6+sceCdthzrmma1VVSAVtVeEkgkPxPYWGsCg7KWCSpgxzrzABTvuXPBK7VoA70UODDOXKknm5SUlOqzP+uX4Jx6tc7JCk4JP8clotYjwUe2X3fddctAVfhAhG0CjczBBZ0Waq6JtZlC9e7RsVzBcT0nHYWAtUFB1VZUAVpVe0GgACLplOKZVckyBZN0XLFOyEJlVd5yyy3LEEuNXKBjASbAJ6koWQ813rjkgqa3p/HGZB01jatZkXHjsiDjwgVzgHac8Z6955pNv8+ktqRfKGCm1Zl7FGkc9/L0elOI3lj3nTQZ0b1Z462qraACtKr2Uj9P0AFIa4NAktZgKSoAcCAFJimYYFyAl8IDrLpYjSzLdFBh3QFnrE5KFG7K7QEZ0DqO1ceCdExKBSZy13Nqz4JqApC4iBNQBHgiilmtLGTgjJVrv2v2JQF03Rv3rmclCH2ZSEpPitknb7WqtooK0KpaowIIwMyaIVgl7YMlp+A6aAi8ARGWmZSVADT5l3Gppq9m8iVT4i/1Z2dLV8dV6xzmYy06BuiIFei6EiyU0nxZj01EL3CCntQTqSmOcU/gB6TpsAKWBx100HgfLFqgTsm+1LklxzpPAop8sWhKSrXVVIBW1RqV5tHWMEEk1YNYayBoH3mf9nzWE0EV9AAv65GUnp5ZwwTWwDGdUpKCAkZZL2Ux5pi0JktgUaxTSkWiNM9OFxbn8dpzLEfnFRDkOYFJoJpat4ceeujwyU9+cjnAKVGzvgC4X8+BdkvwVVtVBWhVraHurfdgCZQq7bDkwIv7khXnvSICLLBEuNqXlJKUrUtuZPpgAmzWHlO3NgFAce0mSCn9O2OlTjud2J7WYiAIlFNXrTHpNWqclBnbgDGlAl1zLEnbPHPTBuLu0bz5mcStnGCiRNw2aKjaqipAq2o3696CifcszM997nPLNWpBCxi5Ue2f1p8FEa7dFI4PIBMwFCCmqHpq3sZaTC9N24Epa58BoTmANcXeA/SAL9sI9FPnFmBdv0f6euqn67150/7M9QOq+2elclt7zXXrvkAzLt12Tam2iwrQqtqFgJOblhIIwzIDQiCxH6iAkwWYqkGgBkqg6risbcZSc0wKrxNApll21kUT5Qp0zuncxifgKBZp3LWglVQS4Et6CdAZ75qtd8Yqtc81m5+FDIJAmR6fthsH9uYg8yvdR2Ap5cW41Oaty7baLipAq2onirsyRc5TkUcaSLqm2O45eZeAClLyI+OOtQYKYImsnVqcIJRC6iCUgvABIoFv1hkT2cvCBVUCZS5acEuh+XRaicXqvWCepNY4Hgidx7WCqevzWvnB1MQF4cAw1uw0HWXazq1dU6rtpAK0qmbWOQEHkLKmKN/ROmfWMFmHIlET6EPg6OEYcGWxAqq1RdtTuSfgjLvVNuOSK2nctCA8pYpR3McAKjLWvElNSXPuWL3GgqQatSANkqlrG7dwrGH7WZBxQQNmiiH4IjALzqllWVhW21kFaLUttKvm11nnZA2CR6r9pIA7iHqAk6o9iiBwcXLbgiT4Jf/SXNYSyba4brMmGegm2jWBRdNrAUX7stYJbtPnBPek72jgnMIKJGo24ARc6TPuX3qNtdtp0Xmg9SUh6SruK51REmxknrplq+qfVYBW27r59Wx/T8AEF+5Z64HgBFqKCHCtGpN6r8Zw2YJuwGweQTUgDLjA5njnB9CsWUaJljX3tJcnOXd6dsaizBosUBovL9M1mN+5wdB47lyWJFiCMYhrTcbCtD/pLKxlr8HYXK7ZfacNWxpyp+hCVVX/rAK02rbNr6cKAFOQXYQt+AEQ+DgeRLhx5T8mKIj1OY3OBVXPXLiBdVqbgdwsJJ0jJfimAtlYolMrFCATRWsMS9GxSV1J/VprnSndZ7viCO7dsQrEs5B9KbDtyCOPHAFqvuSTpuRf3NixcFN2sKqqArTaxs2vpwrsAEjAT9JNAIrlBlws01hkQMnl6TgQBbisXToPC3HaBYUCxKnA0f4UTEh927xPgfk0xGYxpoG2+WIluyZromDM8nTPgp3MlyhflmbSWVjJ+bIAjtNiCDT7BSM/v+Z1VtU/qwCttmXz63kBQ7ZrEm0dkJUGKqxLsGF5siod5zkpJCJwk3JCYJdAnZS1o6SjTAEKygFm5piuS8YaTGuyWH/ZRizarHPGXZxzeO8avHefxvlSAKQBJiDnddJoUrPXnH4u9jW3s6q+XQVotS2U9IpYiLGowIE1yeUJFuAh8pR7EyhZlx/5yEeWU06sIYIImLBAY2mCkvlBCpBZd46JtQhslACigClpKwlaStm+BBslOjdwNa/r4KZNjVwAdB4uZvBOEQduWtYnaJJ5beeKZVVz68Y69UXBNfk5TGE578tHVVX/pAK02nbBRKmiA1rpiJKWY7YDFtfozTffPEaugibAJPUk1po+lwAWKy9wSc1aoErAkDGxUinpKyDp2FioicZNEYVE1SaK17nlbgr28RokgfPwww8frWFjvDe3LwHuTQSuca4BMB1nDOvTmNx7rMtZWDa3s6rmqwCttl0wUcrlxXpjtbFCs2YJjGCkpm3cpWDGkkvfTFB1fHI8bWMdJo8zpfhSdYjidqVpJG6Cc+LiTa1bSr5monBdq2tL+gqrkpXMIpWK4v5AF2R9CUg/UHJfiiSApnl9kfCzECiVNVcqLKtqdSpAq20ZTMQqS5NrRd8FC4EPQAJSgOp9upx4HasTgFK8IG7bQM/4BAbFjTurFKMHzIydpqhkHTSFDuwDSXOlKpHtYAeW3M7HHHPMuJ/r+ZRTThnTbcg47ubMCfrO4WfgGoDfz2a6FlxV1a5VgFbbRrPreZSOKqDDqgSmPAsesj3rkwADQHHnJoLW9rQNyxpnWpMlB3Se4p6NlZkyf0CZtJdcI7Al8tY6ZqoR2X7UUUeN7lhANx+gJqI2UbPmSIF4VqnXLFFWqnkct7M82aqqvl0FaLWtlNxGcGJ1csumQILX4JLgHFaofNAUKHCMtcS4WG3Luuh0/TIC0bhPV1JcxGkhBnCxNKduXgAUHezaU6WIxZlgonRmES3MJeu6Xb8x7hksWdy2OwdAp0Yu7SpPtqqqb1cBWm07ASJ4cmXedNNNI6gAJPVq01waLOOqNSbRtPYnCClRs0knmdUUgispc1N6f6YdWSzTrI06P4gCovfSbljMHq5TNSH35rr19uSe9UiUsehbmpblSw1ezyzZpqpU1epUgFbbpgYuJfcTpMDT6/TUTANr4EyD6XQu8ZzatVy+FLCl/u08Jc90NUqqS9zMAJcCDalCxD1rH8uTtcxitA1AuWO5dw844IDlvqC5LlBlgWbdd/rzcf3pBDMNJqqqaucqQKstK1AAw6kbNa7KBOGkKEJ6baZJtvdJARH1KkgHfI13bNZFWY5xo04hmtSPFKPflSs3zatZgIAI5gDptTmALbA3zvWAZxpyp7OK6056itfmdX+qE82zLBOZnAIN5phWJKqqamUVoNWWFGABoahZoEzlIA+vweJBD3rQcs3bpKOIVmWtAWOsv6wHpsJQXLxRLNKpAs1YoAGU7SnVl9zRKbDA2LgE8ijHZ7yAJmugce96BkWWpvuMxZkUFfcDwLZZ7zQ+c0670uTZPbuOum+ravX658zuqtoCAgwwkO/IDUuKv2efOrfpj/mJT3xitObARQDRtddeOwbeAFDcmrH4WLIpNjCF51TTwgPE5RpLNeX0UkAha5op5+fBkrQvrmTABN+4beVtsi6NAU3jABB0p+3LzBsXtPEAmrXNdKVJC7UUwQ/wd+aOrqrqzqoFWm259U5VeqwJeg+SAKEoAisuBQg8g411Q3mg3KHAAiDTnM5UKorlllSSlc6f5/TmnFYfMk8CjkAvheJB0DV5ZlVmHGA7nlUcAKcrjEhbY1y3Y7ia09HFNZs/AVHmzevU551G205h2ijcqlq9CtBqyyhrlEABQtyzKWTAYgNQVpwxABPXJeszcOSOjUWW8neJhF2pKMI8ARJQTvM7088zfTa9zz4WIjnG6wCTS9a1g2Xq54K+16zbRz/60aObN+ud5k1x+ljCKQYfMILwtHTfarrVVFX17SpAqy0jIODupFQLAoQAMmudIlWBzXauTqAEIDBNpC2BUFyaWcNcrVJYIbmiCRIC7rhcwRpMQc044EwNXOuziZpNBxiWJqAmJxU4paqYmwVKAWkqDaV4/rRg/Lxaty0YX1W7rwK02lJpKgEgyChlFxcnKMSdC2xcvAAGmooPGAdKKYhgnTTrgl7PCxTamUAQqFPhB9BiQXoGKvOCqnGx/FI/1z4F4u0zl/xO2+NytT5qrCApEblcuq4daAPT/Exm25HNA2QLxlfV7qsArTatknaR1BQC0+RTTkEBqKBkvTOWnnVRAGK1Jf8RMGMxplhCGmGvpihCrNe0RksrsmlR+FibKQVoeyznfBkAQfB0L4ofpAi+7VmzJPvtY0UnaCg5pNOep7Uuq2rvqwCtNp2lGThY04x7leVFSfFIj03WGUiKsAVKILn++utHiInUTdWeRMCmNi0opch70ld2JdfGJWuNEuS4i8m5UhSBuGFZlFlv1REGRK3PgiaLmfXstTHug/VoTu9TwCHFFVijoodZ1SoTBa75cjFNV6mqau+pAK02naWZqjlxtQIMC857FmYgBqLACk6O4cb0AGKBOcn9TC6m47JWmrzK5F0mmChK7doo9WvBDBBBO1V/nCOBPYDnXOZTXg/Undd1OadjwZWlaTvLMkUVVBgCUMccffTRo0XpGC7ozJ2iDwkWokbWVtU2ygN9//vfPzz96U8fv6n7pn3RRRft8pjLLrtseMxjHjN+qBxyyCHDBRdcsC7XWu17gQ0IABorTXQtMKbOa4JgCMQCVVABl+uuu278O3Icq9NzaseCDGDa77goBRNAEsCSuxlN4UnGpIOKuVmHYApwrL+U2AP6BPOksIOSgrm/tC0DUtax+2ONpvUZkAKk984n2Oiggw4aoWsu50iUbgDuOeuhVVVtcYD6ln3ccccNr33ta1c1Xp7cGWecMTzpSU8aPvrRjw4vetGLhuc+97nDxRdfvM+vtdr3CgxSlABcvGZtgR5QWc9kfQLYjTfeOELMOCBkxYGUvyuvgcnxIBi4eZ2uK8n5TJeUlOmb5nROle4ocbd6xN1KaWtmXtfri54WZEBp+7HHHjuey5cAgDTO8cApOMgY21nA7tX7rNGSdU/bgdW9ZF/XPqtqG7pwn/a0p42P1er8888fQ/pf9apXje+tA33gAx8YXv3qVw+nn376PrzSaj0UGKRvpUIJBHbW/liKnlmnLDmVh44//vgxTYUFZz8L0b5UF+K+TW6l96BkTDqupCLQagKHArysn3qk6pDzE5j6EhAAJrUksDvyyCNHqxr0E7QEiMCaTi9Z50wuaeYASfPw2MxeV922VbXNALq7uuKKK4YnP/nJd9oGnCzRleQDb9olwwdTtdiBQ8mJBI+AiTXJqgQ+VuiVV165XG2HZyJBRkDiuAQXxeWbgKTANmX6bF8NPMFyWhDBMWAHkmAY+CWFxbwsZ1/4vOZu5YYFW8/+DlPByP3YRr4UTK1K12r+umarauO0JQAqGT4fNFE+jBIEMqtzzz13eOlLX7qOV1mtNXBousYJjqATSxRUgcRrEARTv3fPqf7DrQmexmatEdgokEzE7PT8u1LqyMZF6nz+1lyfa/LwdxigO3+6oqSgg+MPO+ywEZof//jHx+NZxnH7TisHEbduclNZs3XNVtXGaUsAdC0655xzhrPPPnv5vQ9d603VYigBPgR6seACVO5MgBHtCpbJgwQsBeMToRq4Wm/0O54WiKdpz8zdlfODYOraOpfzpjKQ93l2vcn9dA9JWQHEBEmlA0w6qeTaXbf7SpcWUG1aSlVtvLYEQK0VyQmcynsfUPOsT/JB1Jqfi6vZ+qzAZ42QWzMWKAssUbcsSr/vG264YbkBNQs26SG8FHGxTl2za4Fnih94BjP1aLNWCXoBZKoMgeFjH/vY8ViwBMQELsUFnLXc1OENMN2bR6Jx/SwKzqpaDG0JgJ5yyinDu971rjttu+SSS8bt1ebTvAo6rDNQ+shHPrKcrpS+nDfffPMINGueQJQUFtABK5CKGzVFE1ZqSbaSUokokbjpnuI6HvnIR47eCy5Y8BMwlCLygaV0GikswGgN1LEsUtC1jsu1a7x7UsLPl8JpGT5grsu2qhZLCwlQH3zccJEPRukpSSbnfvWB9Du/8zvj/uc///nDb/3Wbw2/+Iu/OPzUT/3U8N73vnf4vd/7veGd73znBt5FtTfWPhMko47tVVddNa5zyhMWaW0/TwPrT37n1VdfPVqa3LoJELPP6xR3X6umHVlYjtYp44IFTCk0AAnUrEQABGnP9idS2HjjUsNWMQUQdd3mt1abgg7T/M26bKtq8bSQABVJKaczylrlWWedNRZI4MrzLT4S0QiW/+E//IfhNa95zfih9PrXv74pLFtk7dN78GH1gZIvUml8HavP34R1ULBKUfhExM4WPVhJaTg9T85BKboAZp7junWNLMmTTjppDApyHVyv5pvuB0/HSrPRScW9WOdkXaf7StZAA03P02jkqqoWQwsJ0NNOO22nKQTzqgw55pprrtnHV1ath9t22tg5OZvpmQlgXKaU7ioqDYFTIAVa4MZKnFYP2pUclw4ss0USuH2zRulv09yxRsE8tW5B3BqtnEwAtY2VbIxgIyUE3YdCCo5PYXjuWTI2XViSYjMbjVyIVtViaCEBWm1PeLIiWWjpaJLAIQ8uWhABkNSzBSXb06IsxQwALrVqV2t9UurrgqX5vGc9Ome2x/o0P1dr3LrWNq1TsiqtZ7oG712ncVy33LTKTDrWvaVzCquTC9h6qW3G57pnLfLWtK2qxVEBWm24Yjla+wZLAAlUQcV6OMuO9WZ78idT/CCl81K0HdDSlmw1xRCiWJ7OEUs0/TkTfJQG2Yng5nq1Lm9fApsAMVHAtgl6cu3AmrxV12l/1jktO8QaBdHMPxuN3MjxqlocFaDVhghoAJC7M0qLMtu5ND0DTdI8uHLBVjF4YAEvAAUVa47pqpJatrsDz6li0cYCTdoKsKfKEIgDHth5D24g6D3Ikmt2nSxPa5r2syzdJ4s0XxJYlI6burEpRecTjey6uhZaVYujArTaEHiyLIEFpKwNAgjgpHC7IDHwAkkWpoIExgOrY0WxxgI1B9CmR2banaXa0O4qrckCy9S4jXs4jbJZkmDJLevZeNej7B6gcueynO0DUO5p1+hLgrmmlial5OC8dU/Q7FpoVS2WCtBq3QV80pAATooSayy5nvYBJkDZBqqsNvAALO5cx3GJghOlFyYBXFp/7a6y9jp9DWZcr7E6yfWA7BOe8IQRmFJqrIFyQ7te3VXSnxTwWa7ux/5UI/JIA+9ZTdc9HTONxKWuhVbVYqgArdZdoAJCWowdeOCBI5iAk1hxrEvPIMJtCzSsNhYo8AJn+nUGnnG1zta43ZWANjBLxG2ChoCPtQlUIKniUFqlibIFfkAES2u33MgpGQi4ooXtB0HXCLZZX42VmfSYqeKidVx6m8Y6TV3croVW1carAK3WXeAJfLGqgMcDMNL82phYmaw7Jfy8B1FwSrUfAiPgS7/Q3dG0ulCsVtvMmaLtwBWgAyJ3bIKJANB9cDHHCjXGa8+uy1ws07Qqc53WQ1eC4NRFG3hmPbT9PatqcVSAVusuAEiBdSANKEHSM1CwSFmTLDvQZK2xRoEs/TKTUkJrLQhPYM0qTt1ZcAR2UHQt9lnPVLCDdZmgJeX7WKGsTnJP4AqqrNe4Yb233fWrPpTAoJ1B0D4QTj3fWJ0FZ1UtjgrQat0FCiwrMFT83XsVhbhHWVjgBZ6eYxUaYy3S9gA0lmKaYa8k+zN2nms3cwEb8LEagRTsQDBBTkAYy9P7VA1Kcft0X3FsXK3JR2WZgi0Yz3PbztO0/2fhWVWLpwK0Wvfm2B7coyn6roiAqkLJ4wQKuZUica0bgho3q9dZpzQuxQZ2Bk8wSxQt8MUqnCodUFK5yHvjWZuOBVJrnp4B07UBovsI0BMolP6zUyszzcDXkn6Stc+qqhZPBWi1rpWGWGjey+VUHB50ktaR9UGAAro0xBZUBEYAaluCcVYj83GfGp880dTHZZl6sDADOuBM020P8LrllltGKzOlBMGTwNP25KmyQGN9Oi6wLASramuqAK32uQBJwfSkp3hvjRPIkvvJOky6hv3J5yT7QMtzLL7VKpG5gdk0x5PFCaRcwwCXfqHeszQp1uVBBx10p2hhrlhroO4pReDTLLtRslW1PVSAVvtUSdfgDgVL1mMACm5x44Il0MXNCmaKKASA9jt2d5U1VPPEfZyuK4AocMl70EsRBwAVJZvKR1y3jhc45LpcP0sajFPkIAUQqOuVVbU9VIBWe1XTcnQgBXxcsClwkFSOFEUARWPAKlYp+MRqTNDP7jbAnl5P0mZYnOZMUJGIX7BkObo2uZ5yOl23MbZ5pNxgri2Nr1PLNlrJVTtd/y1Yq2rrqACt9rq1yUKzFpgye57BEWxSWci6IqgAUQomJPcyaSJe22b/agojOG4abZs80+R5xk2cDir2AyWIctNKU3FtQBr3MsvZfhG0gT8QG787P5OW4KuqracCtNprSgk6Vhurk4CQNQla6WoimMi2wJar1nOaU5sHZOKyDRR3pnRKcawo2pTJy1okkIO66/AssCgF2vXfVM+WSzedXRzrOXVsHZ9UFtfFUgXg1f5MqCX4qmprqQCt9priogQ70AEa1htrM5GtJCAINLlQAQ04PYMn4CYYZ57bFoCNn5WxsTazlslKNFfOkaIIoOm6BAHFhZuqQSmCIOiJxZlj0tsz/TtX645tO7Kq2roqQKu9vvaZAgKsOUBUaShrnECini33KEs06SSA6RE4rtQIO7Cd59KNqzYNsdOf07U4v+OSJnPwwQeP50jBBBZzito7Xo1eQHcv9ieXc3ddsC2GUFVbVwVotceat84HToAp11OhhHQqAT6l+VifqUaUoghglyLxa9G0OHvyRtPqDFxdjwLvLEy5p0ceeeR4LjmprFHXxNoEO5apfVy1CYhaKwCbB1pVW1MFaLVHShPsFCpgYRLLUyUhrlAAA9eAFkSNA7iUuvNsDTLH70y7CihiVbIWpZlM1zRF0wIkCxQM03Q7dWvVnk23FMBLofgUYKiqqpqqAK12qZXSMLKde9SaJ1hlzRFUr7nmmhGeQJp8T2AFzmknFvBk7e1uD8+sh05durFozQeU8jUT2HTccceN59WvkzXqurhtwdazgCHXz307rSRUVVU1TwVoteY0DNsEA3mAKMsvTbGtKbIoP/KRj4ygZOVpnp1C8MbHdbvSeueuBJaszOR5Zv0UPAEVqDW9BnfXo6MKkKZ2rsAh21UZOuyww0brN4FDhWdVVbtSAVqtOQ0jtWtTYSh5m17rrpIUFdvBE1hB1fusee6pgDMdVIDe+UHee9fpOgQMOa+1TtBM0XjWKRnrWpPbuVKj66qqqqkK0GrVaRgpV5f1ThJ4o9gAyLLmWJaB6k033TTCKfvBzmuu1wQM7anMZV0zVYbMyyq13mlNUxcV5xXE5FqU7gNRQAXJrIWmkIPgoaaaVFW1Gv1T3H9VraC4bdOCTEQtGHoNQGkm7X0sO1G2oMR1yl0KqMnPTK7n7pbmc655blVzpqIRMEo5AU7RtsBqfTPVjNxDonKNBUsyJjmjsY73pEF3VVXbQ7VAq1UFD3ktb5OVmWIDxIqUAgKyAoe4RVmlxqU/JvCZJ6/XsuYZeDpfAoZS+MB7IHVe12UdE9gf/ehHL5fiS4BSLE9ABXOQT+pMeoe2YlBVVatRAVrNBSalh6dcyNSrJe5ZwGHhWeMEpXQ5+exnPzuugXKVxhJN6goBWaoVzdPO9ilwEBcy2AFoImhdH7dtCr+nu4prY/EmytZ1sD6JpUruxZi0UPMloG7cqqp2pQK0mhtty1pMriZIssbUi7WuyeIEx2uvvXbM8wQfa42xMu0LCNO1hAAvhd5XKsnHOkxU7VS2pfi7B8CxNEHR+cEydW5TKJ716Zzcua6HtTmtbUuxks3JgjZnKwZVVbUaFaDV3GjbWSswRQVYeKxS77lpQZS7NGuRyvSBX7qupAVY2opFs/DkenVc1lpBMOA1d0rs2Q6SmTfuYYAHP1axXp4p48d6pqx/5lqS5wmwjpMTqhsLi7XwrKpqNSpAq7lFz6dpIVyaWQdlvX3iE58Y4QmkrFNWqm2su/T/BMN0SPE6OZ8rCRjTuJpFmOpAATmwOXfcrKDHwsx5ABNoH/vYx47zkMhbx4Aot63rBd90UgkoPTsm66NVVVWrUQFajYo1Rmk9BkqpBcvKA8tPfvKTy6X5POwDKgE5gnisk3LRep/qQqzFXUW1pkoRwAFZ3L+UICFwZV1mnRYcuV1Bz2ugB1YPgUMZb63THAl+mnXR2s4FXVVVtTsqQKtRYAWACahhBbLaUhg+bl1gBTBl8ozRaSU1ZdOAmkVIWVsEx53VrwW6dG8xP/ja5pg8WLLgCoBK8bk28BQcBNregyBrGNCN5VoGStdEwEoJbGrFoaqq9kQF6DZXciRBJcE/XLBgk8LuxtgGcKDDElVVCBiBDahE3Xo/7emZYKRdFX9P/87kYLqWANB7+52HZcrVetRRR43rlZ/61KeWLVfXBZ7myTpqXMHmSQH7pKbMK01YVVW1OypAt7FiUQIigU0KwwMQ2KSSUMZYR3SMaFzAisXo2bEsx6x5AhOYAdxKMgbcpqBNwJE5gDwNr7llv+/7vm/c5pyPf/zjx3OynO1nDYsQztqpgCfzphCER74UNN+zqqo9VQG6jZXIW9YeYHpwiYKmqFSuUACS23nzzTcv54KCF0CxQkEJvLhvzWMMiCZNxX5gjHU5K/sAMWunEdeviFoWZ1JRFIN3bmudSVdJkBMgAqwvAY5Jt5c0waYECDl26qKuqqpaiwrQbaxp5C0IpUgBcKafJ8tSQXYQ5QoFRZYeeAGScd5nnTN1bu0DqViWSUUJTAkYk17C6mXpgizL0TxctqDIQtRVJRYxSAMlgHPtgmrSUgLeuJDjok3QU9zKdd9WVbWnKkC3scAD2Lhik/+YbipctYoksCy5b4EULAEU3MArrcE8c6Oy+FLUPWOzhsnq85zG2SxO0GNZppNL8jqNBVTwBO1TTz11OOaYY8ZjQC8dXUA/JfrAFzRdhzECiqau2emXhViehWdVVXuiAnSbCijBBOBAx3sA4vJkTYIY0EwDggDplltuWS6kkOAjx6fHZlzBU7E6M795vea2jevYumUCl0DTOPAFSBG3HoDK6rQfsJPu4vqlqQCoXFT1b1mgs4rF6R4Kz6qq9oYK0G0cPKQ/Jghy2QoKAkwABShQSgeTQE+VIaCTPgJAoJkAoiiW5GzeJ0sR2GxPCzLWJZcrMCYK17k8bJcqY002HV8SmZuAoEAwPUDt2xkc8yWgqqpqb6gA3cbBQyw+z9JXAFORBJYcdynXLTgaK4CItZem2LZxxYIiCzDVgaJZeHLrgmaCiECOlZv1zUTRxvIF1xSo9941Ze0SbK2RTiHpejwED3kuKKuqWg8VoNtQqSwEeixO4LnxxhuXy/AB40knnTS6arlFBQmBYgrMA2sq+nhOasjs2ioIpv0Y6DmnY1N6z9pn6usCMjet7Vy4gD6FpIAm1qhrTO3blcoQNrK2qqphuzfUfu1rXzsGmbBg5Px9+MMfXnHsBRdcsNyRIw/HVfNblrEmWXyJvGUFXn/99SNAwfTKK68cg4i8T25nUl64SmMhpvzfVFkzTTCQgB5QY2U6xu8GNNNqzPYjjjhi7PZy3HHHjeUDBQyxWO1nDWcNNHNM1fXNqqo2Qgtrgb71rW8dzj777OH8888f4fmbv/mbw+mnnz6mVAgumScfyPZHidCs/lmsNGBUOYglaP0yLclAktWpfuynP/3p5SIGgJvyfGCbBtQAB5K+qKT4Adh6n2hc0EudWRas34ltyflMIYcjjzxytIT9bhNUlNxO8MxcgD0PkHXbVlW13lpYgP7Gb/zG8LznPW/4yZ/8yfE9kL7zne8c3vCGNwwvfvGL5x7jwzmdOKr5SlF4a5npZgKgAoNSWchaaIoOsARTNMHPFwSjFEnwbJ64TgNTlqrX9qVmre0HH3zwCFiQjiUat6xjWMMpjuA9SCcwyfW1Y0pVVYughQQo6+aqq64azjnnnOVtLJsnP/nJwxVXXLHicQBw4IEHjhbSYx7zmOHXfu3Xxrqp8+SDeBo9moLjW9lty3oDM25QFqifKWixELlJWZwAx/oTWJSgH78PP1PwNAdoevbzNpc57PdwnhSCBz1wZL2aE/gAMzmh1itjxaaykLGHH374eE5zO0+K2Xd9s6qqRdJCApS140NfIMlU3osInafDDjtstE7lDLJ0fv3Xf31MwL/hhhtGl+Sszj333OGlL33psNUVeEpZASnWprVMMPTMGgWuVOsBVs9cutOuKCn4nrZiiYhN82uRsyAMeB5gmQbYwGcfq9PvQiF4sLXWCoj2p6G2332gbO6saXZ9s6qqRdNCAnQtOuWUU8ZHBJ4CU173utcNL3/5y79tPOvWGuvUArX2tpUEeCJowYuL1hcLkEozbGAFSlYiuAImsNnG6gNCc4BWrEOwY3l6T1mT9PMzR45LIQbQkz/KbRvXK7dsUlZsk+8pcAiEgTRBSEDu9bxApaqqqo3WQn4yJc/Ph/5U3q92jdOH/QknnDC2vJonH+xb2R0YyxM8uT9BCPysL/oZSk8BO+5bYGV5+pknMpccA5SsTj9PFqSfme1JJwE5++PWBbzUtk2RePsB23lF2ia/1Hnivk1QkWOn3WEUebAtIK8VWlXVomgh01h8YJ544onDpZdeuryNZeT91MrcmXzgXnfddaNlsx36ec4WLwAo0EpfTCBLrqRHyuqxHIGOVQm2KbgekAIWiHpOHVz7p707U/c2x4Jtmlk7zntjwJAL2fkdL/jIuKzFxiXsd5bfm+tMkQbHOV9VVdUiaCEtUOJePeuss4bHPvaxw+Me97gxjcWHcKJyn/WsZ42RndYy6WUve9lw8sknjxaOD9pXvvKV47rfc5/73GGrKlGp85pDg5NIVuAEUhYhy8/PRFECihUIouDGOk3aSlKAzMfyZF0m19L4vAdhAAZJ52PJes8t7NpA1nkB0RjzW6fmtlU+0HVxnbtOx6VDS+QcDSKqqmoRtbAAfeYznzmu1f3yL//yGBF6/PHHD+95z3uWA4tAIOtw5MNf2oux1tJYsJdffvmYX7hVFRcqpTl0+mB6BiWuWlBjIfoZgRhr07YUY+cyZX2miTV4Gk+2gVzcrCzGuHSNBdcUSgDHuNhdV4KDQNocgocEEYGnAhmAGcvZ61ilubdUGGoQUVVVi6j9lqZdjLex4soEHWDZTBZorLNpygfY+AKhohBg2k/uz5cPgT2p1sRC9cUjQT5gmv6e5vFzYSX6uaTgfJpr228O6UKOty3RuIANmOAHuF4rEQjG5rB9Z/cQi7PQrKpqEXmwsBZotWtNS9gBT9J/SEUm5fqATiQuyBnHDQ5mhx566BiUFSuehQia4GXOrDt6n+IJngNCVqt9ab4NyCzLRNJy55o7Llmu+IAw0bmz95AavbNWdSsMVVW1iCpAN7kACKBYcaCXqkFAx90t19I3KNBLezLAA0/jWaeOSUcVVimQpdG27Y5hTSqrl6hccwCzZ1HTrFgWq3xc501VKPDz2rasr85albmH6b6ueVZVtegqQLeAgCeRrcCYqkDWP1MRSAEK65fcq9Y/uSYIJJO7CYTGsy5jGZoHALldvRYQpNQfeCYi1lxgZ5xj4+o1h0AvykrBPGsyKTcgLNioheGrqtoMWsg0lmr1SjeV5M2m76YIV+uaKa1nPyswReMBTUoJt2wAys0bizMABUJKZSGgdqz3goJA99GPfvQ4NzcwyHLrWkMw1jXsKgUlFZHA3HPq4haeVVUtsmqBbmJNA3DAh4UJiCAWAIFj2oGx8GIpfuxjH1uOrE1OaJpj2+94YEsKCug6h6hm83LPGg+e1j2Ntx+IpRI5Ntama+JCrju2qqqtpAJ0C3RWYeUJ4gE27lQWIJglIhYgP/vZz44Qs46Z4BzHgxur1BiWJ/hKS0kdWuAlY1iy5rXfWPvN6WF+qSm2g6vzmsv5dhVNG9j7QiD4aLbfZ1VV1SKqAN3ESuTttddeO8JHezBWIWBx5yY1BWStWRrPJcvCFIUrMChFFADTcYAJqtJWWJ+O8ZycTXDLOiVQmg+IAdCzbazZtCZLusrO3LEJIjJv3MdVVVWLrgJ0E5fvAyzAtP74oQ99aIShfYCYrijAyWrkguW2Nda4WJusykTAJsIWFFmT1kTBU1CQmsIBb2AHtGmmnXQYOabcwI4B48B1V6kobYhdVdVmUwG6Sdc9U3CdZQha1iEBTvUmYnmCGNdqUlUAz36gAlbHpNMKaKZKkcCgtC2TlgLSLFpW6THHHLNc5J17Nu7fVCayD3TT2aapKFVVbVUVoJtMKTTATQqIQAaOAGgdNLAEMWugaVQtWAjcvAY91mEsTlYpa5XLlVs3AUBHH330OKfzgLH94C36NvODcNqS5bj0X233lKqqtrIK0AVTciKzFjj7Ps8sUK5YEARGFqKxKZuXaj6CeQT9CBxirQKd41Kiz9onAAKkuQCSa5fVad5Ym+DpwTol4ARQ88r1TNkr412fYKDCs6qqrawCdEG7qxC4xbKbWnNcpYohCBICTCX7brjhhuXjWYdcp4CYrihx0To2pfiA0rZ0Tck6JvcrILJOWarOGxeu+VwDS9UYc9hnO5du+3ZWVbVdVIAuoHsWhNIhBUBZdVyzsepE3V599dWjWzbrk45L02vgTY9OViFATpuIJ42FHM+adHzakwEpqzXrrPbHYnUdrgdMnQckU3XIHKlzW3hWVbXVVYAukOKeVSGIgEqaCmgBGzgJClKij/UpbQVgWZXglbXOWJ+sxETKmtd81jq5XDM+DbPjkgXOI444Yvn83LYpw2c+czk3axPQaRo9W3hWVbVdVIAukOL6BCbWqACepKywFIGV9ci61G2FhRo4sjaTmsL16niuXNtYirYDIisyVYMcB7bpFQqWrFjvze341NUF1rhuAZbLNu3PoqahVFW1nVSALpgCtbQUE/CTfEpiDQJb3LXe289S9ZyCBwFi1jlBlwXr4bjkfQKsbYDIAj322GNHgLsGsDbWmilYAqT558Gzqqpqu6kAXSBNI25ZnJTUE+uYSUEBMS5X29KKDPRsT0cVMGVFmi9jAZW7lxuXC9Z+wLadlQuM9nEXJ8cTPM3ptetKY+6qqqrtrgJ0ASNwQQ3IwFPBA4BjkQaIxqY3J0uT5chSTIAQeIIoi1GwD4B6b35zc/ECtP32gWNK7xnn4VwptMDiJcd43zXOqqqqAnThInAJBClRrrYDWurOemYlpvUY4E3BmtzPQFjZPXMBLVcv2NpnPOgaC56ArLsKwLJmBSBx1eY6Er3btc6qqqoCdOEicEEqLcesf3qwEOV96vEJfMZxraoIpEYtEAKedBbBQSzNrJkeeOCB434gNI9jFEYATtscx0pVBB6QuWtZmkl58fA+1+X8gN6i71VVbXcVoAsWgZuoWdDjvgVDwBORy9JM38+kmHjOOqgHKxMoKVWJANDap5J93qfcngAj5wTQRNlOlUChFHoHT3PNFnaoqqrajipAF0RJV0mZPiAEUoDUCNt+rlluVQ/7uHJBEWTB7Pjjjx9TXTzSNYVbNh1bHJPqRixOViqQ2rerdmNpsD11M9edW1XVdlYBumAdVrxWqIC7Nu7U1KdNkfjUr2WNsgpZnAcccMAYPGQcqInKZXV6jvuVRQuGom3B2Hop63O1KSmzbuZ2WamqajurAF2AtJU8g9M111wz1rYFK5akfUAIeKw/blyWJ0gCoTVS21mUxtufAKK4WHVVAWfjARlw0/fT+NSv3ZXitnXuum+rqtruKkDXWYFlgnPSAsz7tCNLNSAP1uEnP/nJcR0U9KSWpJCB4wA0FuVBBx00Qo0l6llT7ETSHn744SNUuYOVA2R5Gre7rtg2vq6qqvonFaAblOsJotP0EFG2atxqeC3NBBiNB7wUhBc0xHpUy9baJaByw4Irdy9QZi3VWDDm8k2BheRwes8F6xrqiq2qqlqb/ql/VbUuAqxpjmeAl1xL65VeC/wBRu5VIAVPlqcH9yzAcvECIaA6Rq6n99y93oMqgLJyWaFxEXPbskxjddYVW1VVtTbVAl1HTYNwgIxLleXJdctSTOcU4+R0skZjHQoaOuSQQ5aLJ9x2223LvUFZokBqH3cu2cYSXSlAqK7YqqqqPVMBuo6aDcLJNtYiEIIn9yxYeh0LMtWEuGGNtZ81Kqo2XVbA0xgWrWPB2DhgjqW72mChqqqqatcqQNdR3LRK6WU9kiWaIB+wVDXImifrMwXfRc9y2zrGemUKwsvdNP5jH/vYuP0JT3jCOCbNsx2X0n/JG3WeVhCqqqraOypA10lpTcZdC3LSTqK0FWMtegAiSALgCSecMFqUtgMkqxIQrXeaUyQtaxZQAdL+WJqejQVn8p6VWtdtVVXVnqsAXSexPMETGEXbAiLYZY0SNK1rgiJ51mkFNLlouXaBFARZnixUIJ5nUWZ9MyBljZIgJeOrqqqqPVejcNdJQOghsIdrllX48Y9/fAScnExWKOiBZNZDQfKWW24Zo3HTKBtwbU8e6K5K8KVhtkfdt1VVVXtPtUDXSak/C4633nrrspV58803j5YhaxTkWKrAqJyftVDgBFpyfHJJV1t+DzDbBLuqqmrvqwBdRyVfk2s2VYgEBYFkonFZngSiyf9UJN5xLXxQVVW1OKoLdx8rlYFSws9rUbUgKFcTGMGTNSoVBTS5d0EVXFmnxrfwQVVV1WKpFuheVkCZYB1AlK5imyIHLEjvFUUQOQuWtgGokn0ACaRAq/l1onMFEs2Lnp2er2CtqqpaPxWg+6jWbaxKrz1b05RyEmtSIBDwRbalGPw0BUXAkeONnQXo7PlqnVZVVa2fCtC9qLhoiZXJKrTuCZ6ChHRMYUlyzUplYX3aLqhIcXiAlKJCYGitlFU6m34y2wYt52uD66qqqvVT10D3okAuEa+ewUwQEMuTNekZNLUTU4ZPRK40FXVsI8d5rJSmMmt1BpiOaY5nVVXV+qkW6F7UvIbTyfNkeXLFegCm/TfddNNycQTFDqZAjDXruFi19k2tXM+g2eCiqqqq9ddCW6Cvfe1rx9ZbIPP4xz9++PCHP7zT8W9729vGxtHGH3PMMcO73vWuYb0VCE5hxppMLVtroAcffPDoxj3xxBNHl63XXLnT9JRZazbW5Twrd1cFFaqqqqptBNC3vvWtw9lnnz38yq/8ynD11VcPxx133HD66acvl6Wb1eWXXz6ceeaZw3Oe85zhmmuuGZ7xjGeMj+uvv35YBAWsQCkal2sXMBVPYIGm84rXK1mzO9teVVVVra/2W1InbgHF4jzppJOG3/qt3xrfa9vFFfpzP/dzw4tf/OJvG//MZz5zdI2+4x3vWN528sknj0UIzj///F2eT6F3uZeCdtJTc18pwT++DCSQyCOQbSBQVVXVxmm1PFhIC9Q64VVXXTU8+clPXt6mVqz3V1xxxdxjbJ+OJxbrSuMF4vghTR/rpaxjcu16WOdMyb0GAlVVVW0OLSRA1YNlpSk0MJX3UkLmyfbdGX/uueeO3zDyYN2ul7KOCZrWQK2P1iVbVVW1ubSQAF0PnXPOOaN5noe0kvXS7DrmajqrVFVVVYulhUxjYZGBiW4kU3k/bUQ9le27Mx64NrIo+zRlpaqqqtp8WkgLlEUmxePSSy9d3iaIyPtTTjll7jG2T8fTJZdcsuL4qqqqqtpyFihJYTnrrLOGxz72scPjHve44Td/8zfHKNuf/MmfHPc/61nPGh72sIeNa5n0whe+cHjiE584vOpVrxrOOOOM4S1vectw5ZVXDr/927+9wXdSVVVVbUUtLEClpfz5n//58Mu//MtjIJB0lPe85z3LgUK33XbbGJkbnXrqqcOb3/zm4Zd+6ZeGl7zkJcOhhx46XHTRRcPRRx+9gXdRVVVVbVUtbB7oems980CrqqqqxdWmzgOtqqqqqkVXAVpVVVVVa1ABWlVVVVVrUAFaVVVVVWtQAVpVVVVVa1ABWlVVVVVrUAFaVVVVVWtQAVpVVVVVa1ABWlVVVVVbqZTfeisFmdazsXZVVVW1eAoHdlWorwC9Q1/96lfH5/VsrF1VVVUtNheU9FtJrYU7aZf2Z3/2Z2OPzv3222+3v60Ar6bcW72Obu9162o73W/vdevqK3vhfmERPB/60IfeqWnJrGqB3iE/pIc//OF7NIdf1nb4A6Xe69bVdrrf3uvW1Xft4f3uzPKMGkRUVVVVVWtQAVpVVVVVa1ABuhd0j3vcY/iVX/mV8Xmrq/e6dbWd7rf3unV1j3W83wYRVVVVVdUaVAu0qqqqqtagArSqqqqq1qACtKqqqqrWoAK0qqqqqtagAnSVeu1rXzs88pGPHO55z3sOj3/844cPf/jDOx3/tre9bTj88MPH8cccc8zwrne9a9iK93rBBReMlZumD8dtBr3//e8fnv70p4/VRlz3RRddtMtjLrvssuExj3nMGOF3yCGHjPe/Fe/Vfc7+Xj2++MUvDouuc889dzjppJPGqmIPfOADh2c84xnDxz/+8V0etxn/Z9dyr5v5f/a8884bjj322OUiCaeccsrw7ne/e8N+rwXoKvTWt751OPvss8fQ6Kuvvno47rjjhtNPP3348pe/PHf85ZdfPpx55pnDc57znOGaa64Z/6g9rr/++mGr3Sv5Q/7CF76w/Lj11luHzaC/+7u/G+/PF4bV6DOf+cxwxhlnDE960pOGj370o8OLXvSi4bnPfe5w8cUXD1vtXiMfxtPfrQ/pRdf73ve+4QUveMHwwQ9+cLjkkkuGf/iHfxie+tSnjj+DlbRZ/2fXcq+b+X/24Q9/+PCKV7xiuOqqq4Yrr7xy+IEf+IHhR37kR4YbbrhhY36v0liqnetxj3vc0gte8ILl9zt27Fh66EMfunTuuefOHf+jP/qjS2ecccadtj3+8Y9f+umf/umlrXavb3zjG5fuc5/7LG12+Ve48MILdzrmF3/xF5eOOuqoO2175jOfuXT66acvbbV7/aM/+qNx3F//9V8vbXZ9+ctfHu/lfe9734pjNvP/7O7e61b5n42++7u/e+n1r3/90kb8XmuB7kK33377+G3nyU9+8p3q5np/xRVXzD3G9ul4YsWtNH4z3yt97WtfGw488MCxgPPOvg1udm3W3+ue6Pjjjx8e8pCHDE95ylOGP/mTPxk2o/72b/92fL7f/e635X+3q7nXrfI/u2PHjuEtb3nLaG1z5W7E77UA3YX+4i/+YvxFPehBD7rTdu9XWg+yfXfGb+Z7Peyww4Y3vOENw9vf/vbhTW9609jV5tRTTx0+//nPD1tNK/1edX/4xje+MWwlgeb5558//J//83/Ghw/a0047bXTrbyb5e+Rqf8ITnjAcffTRK47brP+za7nXzf4/e9111w33vve9xziE5z//+cOFF144HHnkkRvye203lmqP5Jvf9Nuff8QjjjhieN3rXje8/OUv39Brq9YuH7Ie09/rLbfcMrz61a8efvd3f3fYLLI+aL3rAx/4wLDVtdp73ez/s4cddtgYg8Da/v3f//3hrLPOGteCV4LovlQt0F3oAQ94wLD//vsPX/rSl+603fsHP/jBc4+xfXfGb+Z7ndXd7na34YQTThg+9alPDVtNK/1eBWR8x3d8x7DV9bjHPW5T/V5/9md/dnjHO94x/NEf/dEuWxVu1v/ZtdzrZv+fvfvd7z5GwJ944oljFLLguNe85jUb8nstQFfxy/KLuvTSS5e3cXl4v5Lf3fbpeBIht9L4zXyvs+IC5mLhAtxq2qy/170l3/o3w+9VnBSgcO29973vHR71qEdt2d/tWu51q/3Pfutb3xq++c1vbszvda+EIm1xveUtb1m6xz3usXTBBRcs3XjjjUv/7t/9u6X73ve+S1/84hfH/T/xEz+x9OIXv3h5/J/8yZ8s3fWud1369V//9aWbbrpp6Vd+5VeW7na3uy1dd911S1vtXl/60pcuXXzxxUu33HLL0lVXXbX0Yz/2Y0v3vOc9l2644YalRddXv/rVpWuuuWZ8+Ff4jd/4jfH1rbfeOu53n+43+vSnP710r3vda+kXfuEXxt/ra1/72qX9999/6T3vec/SVrvXV7/61UsXXXTR0ic/+cnx7/aFL3zh0l3ucpelP/zDP1xadP3Mz/zMGGV62WWXLX3hC19Yfnz9619fHrNV/mfXcq+b+X/2xS9+8Rhh/JnPfGbp2muvHd/vt99+S3/wB3+wIb/XAnSV+u///b8vHXDAAUt3v/vdx1SPD37wg8v7nvjEJy6dddZZdxr/e7/3e0uPfvSjx/FSH975zncubcV7fdGLXrQ89kEPetDSD//wDy9dffXVS5tBSdWYfeT+PLvf2WOOP/748X4POuigMSVgK97rf/2v/3Xp4IMPHj9Y73e/+y2ddtppS+9973uXNoPm3afH9He1Vf5n13Kvm/l/9qd+6qeWDjzwwPHav+d7vmfpB3/wB5fhuRG/17Yzq6qqqqo1qGugVVVVVbUGFaBVVVVVtQYVoFVVVVW1BhWgVVVVVbUGFaBVVVVVtQYVoFVVVVW1BhWgVVVVVbUGFaBVVVVVtQYVoFVVVVW1BhWgVVVVVbUGFaBVtQX1iU98Ythvv/3mPu573/tu9OVV1ZZQG2pX1RaUHqWzja9vvvnm4Vd/9VeHpz71qRt2XVW1ldRi8lW1DfTpT396+L7v+77hoIMOGi6++OLhXve610ZfUlVtehWgVbXF9bnPfW6E5wMf+MDhD//wD0frtKqqPVfXQKtqC+uLX/zi8IM/+IPDfe5zn+E973nPMjwf+chHDh/4wAc2+vKqalOrAK2qLaq//Mu/HJ785CcP+++//3DJJZcM97vf/VZ97D/+4z/u02urqq2gArSqtqD+9m//dgwW+sY3vjG6bblvo+c+97nDbbfdNu6/973vPfyv//W/xu0idH/rt35reNSjHjU86UlPWt72+c9/fvnY0047bXjTm940vjbHGWecMdz//vcfjjjiiNHCrartpAK0qraYvva1rw1Pe9rThi9/+cvDpZdeOjzsYQ+70/7Xv/71wwEHHDD8wR/8wTj2x3/8x5f3sVQ/9rGPjft2pm9961vD05/+9OH0008fvvSlLw1veMMbhp/4iZ8YX1fVdlHTWKpqiwnIrrjiiuG//Jf/Mq5zTtc6jz766OH4449f8dgXv/jFqwoy+vCHPzxatz//8z8/vj/llFOGJz7xicO73/3u4dnPfvZeupOqWmwVoFW1hSSonsuWXv7yl3/b/te+9rU7BejDH/7wVZ2H+/Yzn/nMnYoyWDc98cQT13TdVbUZVYBW1RaSNcuvfvWrqxq3mu3yRVmaUVy03MLWPa+99to9vuaq2qzqGmhVbUMJKvrsZz+7y3HHHXfc8Ja3vGXYsWPH8Du/8zvDpz71qXH74x//+HEd9Lzzzhtuv/328fHHf/zHo2VaVdtFBWhVbUP9p//0n8b1Ti7YN7/5zSuOe/WrXz1G6UqBueqqq4ZTTz113H7Xu951eOc73zlWNWKNPvShDx3LBIJqVW0XtRJRVVVVVa1BtUCrqqqqag0qQKuqqqpqDSpAq6qqqmoNKkCrqqqqag0qQKuqqqpqDSpAq6qqqmoNKkCrqqqqag0qQKuqqqpqDSpAq6qqqmoNKkCrqqqqag0qQKuqqqpq2H39/y99+73bXtHKAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(5, 5))\n", "plt.plot(ztrue, zphot, \"+k\", ms=3, alpha=0.1)\n", "plt.xlabel(r\"$z_{\\rm true}$\", size=12)\n", "plt.ylabel(r\"$z_{\\rm phot}$\", size=12)\n", "plt.show()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Now define a number of photometric redshift bins. They are chosen by the\n", ":func:`~glass.equal_dens_zbins` function to produce the same\n", "number of galaxies in each bin.\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [ "nbins = 5\n", "zbins = glass.equal_dens_zbins(z, dndz, nbins)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "After the photometric bins are defined, make histograms of the *true* redshift\n", "distribution $n(z)$ using the *photometric* redshifts for binning. Use\n", "the :func:`~glass.tomo_nz_gausserr()` function to also plot the\n", "expected tomographic redshift distributions with the same model.\n", "\n" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjsAAAGxCAYAAACEFXd4AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjAsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvlHJYcgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAnFRJREFUeJztvQd8m/W1//+JtZcl75XY2XtPEvZM2avzx2qht7dc4BZoacu9lN5Cb2m5beEWKJSWMm7LnxbKhhJICAEySMgge9ux471k7Wn/X+drP47s2I4kS3oeSefN60GW9Fj66oksfZ5zPuecMT09PT1gGIZhGIbJUHLkXgDDMAzDMEwyYbHDMAzDMExGw2KHYRiGYZiMhsUOwzAMwzAZDYsdhmEYhmEyGhY7DMMwDMNkNCx2GIZhGIbJaNRyL0AJdHd3o6GhARaLBWPGjJF7OQzDMAzDRAG1CnQ6nSgvL0dOzvDxGxY7gBA648aNk3sZDMMwDMPEQV1dHcaOHTvs/Sx2ABHRkQ5Wbm6u3MthGIZhGCYKHA6HCFZI3+PDwWIH6E9dkdBhscMwDMMw6cWpLChsUGYYhmEYJqNhscMwDMMwTEbDYodhGIZhmIyGxQ7DMAzDMBkNix2GYRiGYTIaFjsMwzAMw2Q0LHYYhmEYhsloWOwwDMMwDJPRsNhhGIZhGCajkVXsjB8/XnQ9HLzddttt4n6fzyd+LigogNlsxrXXXovm5uYBj1FbW4tLL70URqMRxcXFuOeeexAKhWR6RQzDMAzDKA1Zxc6WLVvQ2NjYv33wwQfi9q985Svi8q677sJbb72Fl19+GevWrRMDO6+55pr+3w+Hw0LoBAIBbNiwAc8//zyee+453H///bK9JoZhGIZhlMWYHpqPrhDuvPNOvP322zh06JAY7lVUVIQXX3wRX/7yl8X9+/fvx4wZM7Bx40acdtpp+Oc//4nLLrtMiKCSkhKxz1NPPYUf/ehHaG1thVarHfJ5/H6/2AYPEuvq6uLZWAzDMAyTJtD3t9VqPeX3t2I8OxSd+ctf/oKbb75ZpLK2bt2KYDCICy64oH+f6dOno7KyUogdgi7nzJnTL3SIlStXihe/Z8+eYZ/roYceEgdH2kjoMAzDMAyTmShG7Lz++uuw2+345je/Ka43NTWJyIzNZhuwHwkbuk/aJ1LoSPdL9w3HvffeK1SgtNXV1SXhFTEMwzAMowTUUAjPPPMMLr74YpSXlyf9uXQ6ndgYJhWQYf6VV15BfX09ysrKcOWVV8JkMsm9LIZhmKxBEZGdY8eOYfXq1fj2t7/df1tpaalIbVG0JxKqxqL7pH0GV2dJ16V9GEZOyBJHUcvDhw+LlOvx48fx5ptvitsZhmGYLBI7zz77rCgbp8oqiUWLFkGj0WDNmjX9tx04cECUmi9fvlxcp8tdu3ahpaWlfx+q6CKT0syZM1P8KhjmZPbt24fdu3eLKsKLLroIV1xxhXgf79ixQ+6lMQzDZA2yi53u7m4hdm666Sao1SeyamQcvuWWW3D33Xdj7dq1wrD8rW99SwgcqsQi6MuDRM0NN9yAL774AqtWrcJ9990nevNwmopRAps3bxamekl8k8l+1qxZ+PTTTzm6wzAMky1ih9JXFK2hKqzBPPLII6K0nJoJnnXWWSI19eqrr/bfr1KpRKk6XZIIuv7663HjjTfigQceSPGrYJiTofYHNTU1WLJkyYDbFy9ejPb2dpHSYhiGYbKsz47S6/QZJhbeffdd0QKBopMkyCXoT+7RRx/FlClThJhnGIZhsqTPDsNkEiRoqAkmmZIjhQ5BfaTmzZsnvDzUS4phGIZJLix2GCYJtLW1iTOOyZMnD3n/7Nmzxew3qkRkGIZhkguLHYZJAkeOHBGG+6qqqiHvp1EoNNyWPD0MwzBMcmGxwzBJEjtUhUXtE4aCUlkTJkxAdXV1ytfGMAyTbbDYYZgkdEymiM2kSZNG3G/8+PFiiC2lsxiGYZjkwWKHYRIMjYUg4/HEiRNH3I8iO2RkptYLDMMwTPJgscMwSRA7lL4aPKR2MHl5eaJkklNZDMMwWTIIlGEygrUPoWHdLpR5/MhZ96sRdx1z7r3C18PNBRmGYZILR3YYJsHUtzlQXhBdc0qagt7U1CTGpjAMwzDJgcUOwyQQrz+ITqcXFYXRix3y99D4CIZhGCY5sNhhmATS0OYQl7FEdojGxsakrothGCabYc8OwySQhnYHdFo18nMNUe2v1+uRn58vxM7cuXMhB7/f8fuY9v+3+f+WtLUwDMMkA47sMEwCaWh3iqgONQ2MForuUL8dhmEYJjmw2GGYBNLc6UJpvjmm35FMytRzh2EYhkk8LHYYJkGQ0ZjMyUVWU8xix+/3o6OjI2lrYxiGyWZY7DBMgqCKKorOFNliEzvFxcXisrW1NUkrYxiGyW5Y7DBMgpDESqxih6afk1GZxQ7DMExyYLHDMAmCxIrFqINeO/Sk8+EgM3NhYSHa2togF9TU0OP0sG+IYZiMhEvPGSaBYifWqI5EUVERmpubIQfVO6tx6PNDCAaCqJhagTlnzYFKrZJlLQzDMMmAIzsMk0ixE6M5OVLsUGQn1ZGVmpoa7N24FyUTSjBzxUw0HmnE9tXbU7oGhmGYZMNih2ESQCgUEtVUxXmxlZ1Hip1AIICuri6kCqoAe+2115Bflo85Z8/BhLkTMO/ceWiuaUZ7A4+vYBgmc2CxwzAJgIQO+V7ijeyQZ4dIpW9n165dcDgcmHfOPOTk9H4UlE0qQ25hLg58doD9OwzDZAwsdhgmAUiDPAtyjXH9vs1mg0ajSWlF1tatWzF16lQYI9ZMZulpS6ehs7kTHY3c94dhmMyAxQ7DJCiyo9PpYNTHVok1uCIrVWKHxlPQPK5FixaddF/RuCIhgBoO8QgLhmEyAxY7DJMgsUMDPWOZiTWYgoKC/ghRstm+fTtyc3MxefLkk+6j11A2sQxN1U3oDnenZD0MwzDJhMUOwyRQ7IwG+v3Ozk4kG/LiHDx4EDNnzuz36gyGvDsBXwBt9fL1/mEYhkkULHYYRkFihwzDNGMrmZAJmqq+Jk2aNOw+ZFI2WU2iFJ1hGCbdYbHDMAkoOyeRkpeXN6rHkX4/2dGdI0eOQK1WY/z48cPuQ6ks6r3TWtfKVVkMw6Q9LHYYZpTY7XYhCBIR2SGSPf388OHDqKysFNVfI1FQXgC/xw+33Z3U9TAMwyQbFjsMM0okcTJasUMDQUmAJDOyQ1GoY8eOjZjCkqBmgxTh4QaDDMOkOyx2GCYBYofSQhaLZVSPQ8KCBFMyIzv19fXCEzRx4sRT7qvWqGErtrHYYRgm7WGxwzAKKDuP9O0kW+xQ9KikpCSq/fPL84XYYd8OwzDpDE89Z5hRQuIkLnPy2odOuim/9hD21bYAY4+fvP+590b90L/f8fshb9+2fhv8IT+e2vlUVI9TWFGII9uPwNXpgiV/dJErhmEYueDIDsMkwKA82kosifxcA7rcPoST1MzP3mwXqalosZXYgDEQ4yMYhmHSFRY7DDMKKL1DPWusVmtCHi/fYkR3d48QPInG5/HB6/LCWhz9Wsm3Y8mzoKslddPYGYZhEg2LHYYZBR6PRxh+aZBnIrCZ9eLS7kq82JEES15xbFEoa5EV9lZ7wtfDMAyTKljsMMwoU1hEoiI7VpNeGJ3tLi8Sjb3FDp1RB32foIoWSns5250Ih8IJXxPDMEwqYLHDMKOAUlhEoiI7KlUOLEZdciI7rV0iShNr1Rj9DqXrHO2OhK+JYRgmFbDYYZhRRna0Wi0MBkPCHpNSWcmI7Dg7nHFVVNHv0MBQEksMwzDpiOxih/p+XH/99SgoKBBfGHPmzMHnn3/efz+dUd5///0oKysT919wwQU4dOjQSaW/1113HXJzc8UZ9i233AKXyyXDq2GyDcmcnIgeOxI2kx72BBuUaYK5z+1DbkFuzL+rUqtgKbCw2GEYJm2RVexQW/zTTz9dNDn75z//ib179+I3v/nNgDLehx9+GL/73e/w1FNP4bPPPoPJZMLKlSvh8534MiChs2fPHnzwwQd4++238fHHH+M73/mOTK+KybbITqJSWBJWsx5dCU5jUVSHiLdXjrXQymKHYZi0Rdamgr/61a8wbtw4PPvss/23TZgwYUBU59FHH8V9992HK6+8Utz2wgsviO6vr7/+Or7+9a9j3759eO+997BlyxYsXrxY7PPYY4/hkksuwa9//WuUl5fL8MqYbIrsjB07NqGPaTMb4PD4Ra8d8vAkSuxQKspkNcX1+xTZOX7gOLqT1P+HYRgmYyM7b775phAoX/nKV1BcXIwFCxbgj3/8Y//91dXVaGpqEqkrCUoZLFu2DBs3bhTX6ZLOrCWhQ9D+9MFOkaCh8Pv9cDgcAzaGUUpkh9JYwhDs8SfsMUnsmGwm5MQpnigi1N3dDXcXT0BnGCb9kFXsHD16FE8++SSmTJmCVatW4dZbb8W///u/4/nnnxf3k9AhBs/xoevSfXRJQikSGspIs4qkfQbz0EMPCdEkbRRdYphYIdFM6dRElZ1HRnaIRJqU4zUnS1BjQelxGIZh0g1ZxQ6dKS5cuBC/+MUvRFSHfDb/8i//Ivw5yeTee+8V6Qdpq6urS+rzMZndYyfxnh2duExUF2WKErk6RjfbSmvQih49LHYYhklHZBU7VGE1c+bMAbfNmDEDtbW14ufS0lJx2dzcPGAfui7dR5ctLS0D7g+FQqJCS9pnMDqdTlRuRW4ME2+PnURHdtQqFcyGxPXa8bl8CAaCwnczGkgsUXNBhmGYdENWsUOVWAcOHBhw28GDB1FVVdVvVibBsmbNmv77yV9DXpzly5eL63RJZ9hbt27t3+fDDz8UUSPy9jBMsnA6e02/ZrM54Y/d22snMWLHZXcNSEXFC5WtOztZ7DAMk37IKnbuuusubNq0SaSxDh8+jBdffBFPP/00brvtNnE/9S6588478fOf/1yYmXft2oUbb7xRVFhdddVV/ZGgL33pSyL9tXnzZqxfvx633367qNTiSiwmmZDwJqFDgifR0NiIRKWxyFRMa4x1TMRgzPlmeBwe4VViGIZJJ2QtPV+yZAlee+014aF54IEHRCSHSs2pb47ED3/4Q7jdbuHnoQjOGWecIUrN9foTH9x//etfhcA5//zzxYf6tddeK3rzMEyyxU6yUqAkdpoS5I9x290wWo2jFmW5+b2vtbW1NeHl9gzDMBkrdojLLrtMbMNB0R0SQrQNB1VeUVSIYTJF7OSadCKyQ+bi0XZnpshOvP11IqHSdaKtrY3FDsMwaYXs4yIYJl1Jqtgx6hAKd8PrDybEsyMJldGg1qhhMBuE2GEYhkknWOwwzCjEjsUyOtPvSGks8RyjbCwYDoXhdXkTEtkhSDSx2GEYJt1gscMwcUAmXdqSl8bqFTujnZFFhmL0IGFix2wzo729PSGPxTAMkypY7DBMHEgjRpIldkx6DXJyxow6skPmZPF4CUhjEeY8s+hhRa0dGIZh0gUWOwyjQLFDlVO5xtGXn7u6XNBoNdAZersyjxaKEIXDYXR2dibk8RiGYVIBix2GGYXYSZZnR6rIcoxS7FBkh6I6o63oGqoii2EYJl1gscMwcXZPNplMYuhssqCKrNGmscizY7D0DhZNBHqTHlqtlsUOwzBpBYsdhlFY2XkiuyiT2EmUOZmgCFFhYSGblBmGSStkbyrIMIngkQ8OxrT/XRdOjX7ntQ+ddJNj0w7kYgywtgPJgiqyKLITb2NBKjv3eXww5hoTuq6CggKO7DAMk1ZwZIdh4sDh9sNiTIzpd6Q0VjjcDY8vvsaC1F+Hys4TmcaSxA5VZDEMw6QLLHYYJg5c3gAsRm1Sn0NqLBhvKkv02AESHtnJy8uDy+VCIBBI6OMyDMMkCxY7DBMj1GPG7QvAnKBy7uGQIkckrOLB6/SK9BeZihMJzaIjuPycYZh0gcUOw8QIiQ/y0ViSLHaosSCJFafXP6pKrNFOOx9O7HAqi2GYdIHFDsPEiBRpSXYai0SK2aCF0xO/2El0CoswGo3Q6XQsdhiGSRu4GotJe8g7cnTX52g5Xg2jxYop80+DwZy8snBJfCTboCw9R9xix+mBrdiW8DVRtImiOyx2GIZJFziyw6Q1NLrgb3/7G47s3AydwYT2xjp8+uZf4bIn74uY0kr0hW/UaZBsKFUWTxqL0myeruREdiSTMosdhmHSBRY7TFrz9ttvo6amBgvPuxwLzrkEZ155A/QmC3ZvXJO0YZWUxqL0UqK9MImM7AT9QYSCIRgtyRE7HNlhGCadYLHDpC11dXXYvn07LrnkEhSUjRO3qbVazF5+PrramlF7YGdSnpfERypSWJIvyOkJxF12nugeO5Fih7pIh0KhpDw+wzBMImGxw6QllKZZtWoVSktLsXDhwgH35RWXoWzCNNTu/0Lsl2gorUSRnVRAaSwqc6fmgrHgc/X25klmZIeOrd1uT8rjMwzDJBIWO0xacuDAARw/fhwrV64ccpTCuKmz4XE6hIcnKZGdJJedS1AEiUQFCZ5Yuyer1Cpo9JqkeXYITmUxDJMOsNhh0pLNmzdj3LhxmDBhwpD324pKYbHl4/ih3Unqnpw6sUPE6tshsaM36+OaqRXVuiwWqFQqjuwwDJMWcOk5k3a0trbi6NGjuPbaa4fdh77kK6bMwoGtnyLg80KrT4x3hUzPKRU7fRGkWCuyqHuywZwcvw5B5myr1Yqurq6EPWbrY49HvW/RHbcn7HkZhsl8OLLDpB1btmyByWTCjBkzRtyvtHIyerp7EprKcvuCIq1k1qfGs2PUa6BS5cRsUqbITjLFDmGz2TiywzBMWsCRHSatoOqfnTt3YsmSJVCrR3776k1mWPIK0NZQi7IJUxPy/K6+CEuqIjsUoSJhFU8aq2R8CZItdpqbm5FqQuEw1q9fL1oOaDQaLFq0CJMmTUr5OhiGSR9Y7DBpxeHDh+Hz+TB37tyo9i8sq0RD9QERjUmEf0WKsKSqGqu/104MaaxgMIiAN5D0yA6lscgonkpcPh/+tvkzdL7/PsYXFsId8GPLs89h3rhxuGj27GH/jTntxTDZDYsdRrE88sHBk2774uP34eoC/rKDJm6feup2YUUVqvduh7OzDbn5RaNek1QVZUpRGiuexoLU/yaZPXYiIztut1uM69Bqk388guEwXtu2Ff5QCDedcQaKLBYhYnfW1WHV7t0w6XQ4Y2piIngMw2QW7Nlh0oZQMCDmX8WSkrIVl0GlVqOt4VjC0lgGXa+PJlWIkRExiB3JR5MKzw6RSJPySKzbvx9tTheuWbRYCB2CIjnzKitx1rSp2HD4MI62tKRkLQzDpBcsdpi0ofV4DcKhEMrGRy92VCo18orL0dnSmDCDcipTWP1dlPsmrUeDEB9jyLOkT4nYSYVJud3lwo7aWqyYMgWlVutJ9y+bOAlVhQX4cN8+hJM0JoRhmPSFxQ6TNrTUHRWpqFgnmlsLS9DV1pSQbsoU2UllCktKY3l8gahHM5DY0Rl0yEly9Il67VAJeirEzscHDsCi12NRVdWQ91OE59zpM9DpcWP7scRE8RiGyRxY7DBpQXc4LKqqiscN3URwJKjBYMDng8fZlZ6Rnb5eOy6XK2qxk6wxEUP12km22Gl1OnGouRmnT50CtUo17H7FubmYM3YcNh45Iiq2GIZhJFjsMGlBZ0sDggE/isdNjPl3rQW9Jdhd7c1pG9khnE5n1GKHuienglT02tlRe0yYj2eUlZ9y36UTJ8IbCGBfQ0NS18QwTHrBYodJmxSWwWSBJa8w5t+l7slGixVdrU1p6tmJTezQfsn266RK7ARCIeytb8DccWOhyjn1x1W+yYRJxcXYeuxYUobAMgyTnrDYYRQPfWm11tegaOz4uHvlUCrL3tY86lER5J1JdWRHr1VDTV2UoxQ7VHquN2aG2KEITSAcwrxxlVH/zsKqKrQ4HDjeeerWBAzDZAcsdhjF43HaxQRz6pkTL9bCUjg7WhEOR2fyHS6qQ5hTNPFcggSe6LUThdjx+/2i700q01hSr51ksLehAeMLi5BriL6MnpoN2oxG7KmvT8qaGIZJP1jsMIqnrb5WmGHzSyrifgxrQbGIzLjtHQloKKhBqolW7EgNBXUpGmeRzF47Lr8fxzs7ML20NGZxOKO8DAebGtmozDCMgMUOo3jaGmtFrxy1Jv70kdmaLy6d9va4H4OmnROmFHt2+hsLRiF2pH2S3VAwFb12Djc3YwzGYFJJ7DO+yMzsC4ZwtLU14etiGCb9YLHDKBpKO3U0HUdhRfSejaFQa7UwWnLh6mwfdWQn1QZlJUd2ktlr50BTI8YV5MMYxyiKQosFJdZc7G3gVBbDMCx2GIVjb2kUXZMLy+P360hYbIVw2ttGFdnRCbPw8L1elCB2jEYjVOrUrDFZvXaofLyuvQPTSsvifozpZeWobm0TM7UYhsluWOwwiqa96bgoHTfbCkb9WGZbPpyjiex4AzCnuBJLwmLQimnvpzICkyCiaEsqSUZF1rH2dnT39GBiUfzDWycXFwuhU9se/785wzCZAYsdRtFQCiu/tCLukvNIqEeP3+sRW9wNBWVIYcXSa4ciO7m5sY3TUKTYaWtDgdkUUxXWUD13aCPvD8Mw2Y2sYue//uu/xJdY5DZ9+vT+++lM9rbbbkNBQQHMZjOuvfZaNA/64KqtrcWll14qQvfFxcW45557op4hxCh/yjl1PS4oHZuQxzPn9UaHXHGalOVoKDhY7JxqZASJoXQXO9RXqaatTZSQjwb6PJlYXIwjrS3cYJBhshzZIzuzZs1CY2Nj//bpp5/233fXXXfhrbfewssvv4x169ahoaEB11xzTf/94XBYCB0K7W/YsAHPP/88nnvuOdx///0yvRomkdCk8p7uHuSNouQ8EqPFBpVKFXdFlhyjIiQkkXUqsUORHTnSWNRrJxjs7UM0Wjo9HnR5vaK/zmihVJbL5xefHQzDZC9q2RegVqN0iD4a1LfjmWeewYsvvojzzjtP3Pbss89ixowZ2LRpE0477TS8//772Lt3L1avXo2SkhLMnz8fDz74IH70ox+JqJF2mCoOarxG2+AKFkZ5KSy90QRTbl7CzLQma35aRnZ0GjXU6u4RxQ6JfxIdckR2CIruFI3CYyNR09YqRkOMy+9tFzAaxublQadW48iRI6ioSIxoZhgm/ZA9snPo0CGUl5dj4sSJuO6660Raiti6das4U7zgggv696UUV2VlJTZu3Ciu0+WcOXOE0JFYuXKlEC979uwZ9jkfeughUUEibePGjUvqa2Tio6O5XkR1EuHXkTBZ8+BxxJ5yEQ0JZRgVIUHHgFK5I4kduo/SNXJEdohEpbKoCqvcZoNWrU6IwK0sKMDRo0cTsjaGYdITWcXOsmXLRNrpvffew5NPPonq6mqceeaZwnfQ1NQkIjPSB6kECRu6j6DLSKEj3S/dNxz33nuviBxJW11dXVJeHxM/5NdytLckzK8jYcq1wdUV+8wkrz8khIRckR3CZDKNKHYk83KqIzskrkiMJSJCSseYZlpV5CUmmkdUFRaIv/FkjbRgGEb5yJrGuvjii/t/njt3rhA/VVVV+Pvf/w7DKKowToVOpxMbo1yO9U2tTpRfR8KYa0PA5xViSq/Xx+TXIeSK7BCniuxIYkNEdhpTty6KntBzJmJkhN3jgdvvF+mnRFFVUIiNfp+IGk+ePDlhj8swTPogexorEoriTJ06FYcPHxY+HjoTGxwap2osyeNDl4Ors6TrQ/mAmPShpqYGBpMFRos1oY9r7vP/tLW1pcUQ0MFihzw5I4kd8sAl80RhOCialIjIDkV1KGtZnkCxQ+XntD5OZTFM9hKz2PF6vfB4PAPOwB999FFhFh4tdNZKRsKysjIsWrQIGo0Ga9as6b//wIED4uxs+fLl4jpd7tq1Cy0tLf37fPDBB+KDbebMmaNeDyMflNJMtF9HiuwQ7TE2mjsR2Un9ENBoIztSQ8FEH7NoIO9bIiI79Z2dKDRboNck7jjT8ZgwYQKLHYbJYmIWO1deeSVeeOEF8TNFXSj19Jvf/EbcTr6bWPjBD34gSsrpLJ5Kx6+++mpRGvyNb3xDfHjecsstuPvuu7F27VphWP7Wt74lBA5VYhEXXXSREDU33HADvvjiC6xatQr33Xef6M3Daar0hQQ1ReiomWCioWGiVOEVT2RHq1GLTW6xM1zPGDkaCkrQ32siIjv1nR0Ym5+4qI7E+PHjxXuK0pcMw2QfMX9yb9u2DY888oj4+ZVXXhGG4O3bt+Mf//iH6G9z6623Rv1Yx48fF8KGzrKpZPWMM84QZeVS+So9D/kBqJkglYpTpdXvf//7/t8nYfT222+L5yQRRAbOm266CQ888ECsL4tRoF8n0eZkCSpljzWy01uJJV9URxI7VBVGYpCaaMrVUPD3O078DUpUN1Vj//79yNmeMyCy9G/z/y3qx/UEAmh3uXHapMT7asgLSO8pigxTqpxhmOwiZrFDKSyptJVSV9TkjwQJRVvoSyoWXnrppRHvJwPpE088IbaRPsTefffdmJ6XUTb0hUSRApiT88VNqazY01gBWf06ktgRa3G5hhQ7FFmhNg5yYDAb0B3uRsAXgC7O49TUlwajsvNEk5eXJ44fix2GyU5iTmNRNcPrr78uSjkpbUSpJIJ8M3KF0JnMgr6QqJ9SsqBeOyR2YhkhoJTIDjGUb4deixyjIiT05t7KNp8r/jRRk90uvDq2IYTcaKFoE50YxXpCxjBMloodSlWR14Zy4EuXLu03C1OUZ8GCBclYI5NFUCNJau2fVLGTaxPz02Ix1JJBWUmRncGQF4WOXaobCg4WO16XN+7HaHJ0odRqTZrBmt5T9N5K1FgLhmEyWOx8+ctfFmfen3/+uYjsSJx//vn9Xh6GiZf6+nrhS0mm2DH0pcc6OztjMijLHdmh6kQy3g8ldiRzsFyRHUpdUTrb544vskORqUZ7r9hJFhTZoZEa9B5jGCa7iKvPDvWwoTNIKvMmsySxZMmSARPLGSYeSEjTF3oiZiyNJHYoehCt2KEvYkpjyR3ZGan8XOqeLFdkh44nRXfijezQ+qmZYJkteWKnuLhYvLeoMIJhmOwiZrFDXgeK4pDJ75JLLhGTygkqE//+97+fjDUyWSZ2aFYZRQmShUqlFqIgWrHjC4QQDnfLHtkZSewM6J4sEyR24vXsSNGW0tzkiR16T9EwUBY7DJN9xPyNctddd4lwOn0pRVaEfO1rXxMzrhgmXih9Rcb3ZKawIqtzohU7ShgVEU1kh1ovUDsGuTCYDHFHdshLY9brYEly9+exY8cKsROLOZ1hmCwUO2RE/tWvfiU+NCKZMmUKVzowo6K1tVX0U0qF2MnPz49a7ChhVEQ0kR25qyFFZCdOzw5FiEuSGNWRoMgOHb9EdHtmGCaDxQ7N5hmqx0dHRwd3LWZGBUULpVSDkiI7bm/vtGyTQblpLGlUhJxQrx1KY1GELlaou3FxbvLXL52kcSqLYbKLmMXOmWee2T8uQjIm0ofbww8/jHPPPTfR62OyTOxQUzxKk6ZC7FCDTIoknQqPPwiVKgc6GUdFRIodWjdVFUVCAkhusUORHUoP+T2nPqaDT6BIrBVZkh+ZolQf/duz2GGY7CLmT28SNWRQptJzmkr+wx/+EHv27BGRnfXr1ydnlUzWiJ1Zs2al5LnoC4+g6A5VF44EVWIZdRpZBmwO12uHBEJk2orEjnSfXOhNfY0F3T4R5YmWpqYmcVmcojSc5NthGCZ7iDmyM3v2bBw8eFDMsaLhn/ShSyMjaD7WpEmTkrNKJuMhDwVtqfDrDBY7p8LjCwqxowQoMkFEprJEabzbLbvYkQROrCZlSmFptVrkJaFz8lBQmpQE1uDoGMMwmUtccXmaW/Sf//mfiV8Nk9VRHYLKzlMB+c7oCzYasSMiOwqoxBquizL1uqIvbrnFjkangUqjirn8nIQHDRQe44m/+3IslJWViQ7aNPmenpdhmMwnKrGzc+dOEdEh8yj9PBJz585N1NqYLBM7BQUF/ZGLZEMpKYruUPo1Gs+OEnrsENLxoUiOhCR85BY7orGgKfbGgiR2RESvugapQEpbigowFjsMkxVEJXbmz58vPpCoAyn9TB9qQ/WpoNs5NMwocfjnUNhstqhKkCmyU2RNjQg7FdRHh6JSkZEd6edUCcVTVmTFUH4uRVhozl6qxA5VjZKwJrFDn2cMw2Q+UYmd6urq/vb99DPDJBIaYtnS0oLTTjstpc9L6dho3s/k2VFC2flw5edKiewQFNlxdZ5cGj8cJHSomjPVERZKZUnd3xmGyXzU0Q7QG+rnwXBXUiYepI62ckV26LmHq7Si+zx+qsZShmdnOLFD/iPa5MZgMaC1rjUmsUMUFhYieomUGLFDhRYktJI5moRhGGUQ81/5N7/5zQF+AYmamhqcddZZiVoXk2UpLErBUFfjVEd2qH2CNMx2uKhTd3ePYjw7w4kdJUR1pMiO3+tHd7g76q7Z1B/IkOQxEUOJHfq3j8azxTBMFoqdL774QpiQN27c2H/b888/j3nz5omzM4aJ16+T6j42FNkhRvLtUAM/QinVWEOJHSWUnQ8oP+/p7bUTrdhJ5oT7kcSONJOLYZjMJ2axs3nzZtFX55xzzsF//Md/4Ktf/Spuv/12/PrXv8Zrr72WnFUyGQulEWjidapKzgdHdgi73T7sPlIUkyM7yem1Q2JHjpMkiiRRNR77dhgmO4i5zw618v+f//kfURHy4IMPQq1WY926dVi+fHlyVshkNNRQLhgMnjRYNhVQ6ozev1FFdhTm2aExF5SGIZ8OiR2qLlICNDKCiCayQ5Wb7e3tvZVYMsAmZYbJHmKO7NAX0/e//30x+fzee+8VIociPe+++25yVshkvDmZyqmltEIqobQZRXdOFdmh/Qw6+ediDTUyQmmRHbVGDY1WE1Vkh/wyFNmTI40VKXa4sIJhMp+YP8EXL14sznY/+ugjUSpMHxQ0L4sEz80334zf//73yVkpk5HU1dWJJm+pGP4ZT68deq8bdBpFVexEdlGm+Vi0RqWIHSm6E00XZUphEXKKHYqQURftVJvjGYZJLTnxiJ0dO3b090Shs94f/ehHwrD88ccfJ2ONTIZHduRIYUlEE9lRylysocQOCR064VCU2DHpo0pjkdihdLhczRClaCKnshgm84k5svPMM88MefuCBQuwdevWRKyJyRJISFAqQ06xQ5Gd/fv3D3s/iQmTgiqxJHMtRZpI7CipoWCkSbmrrUtRlVitjz0+5O05X+zEvrY2FE+b3n9b0R23p2RNDMOkjlEZEagHCZkkB7diZ5hooCosQo5KrMjIDgkayew7GLrPqKBKLCmaKlVkKVHsUGSn+VhzVA0F5RS6RIk1F81dDlnXwDCMAtNYdDZOpeY0J4vCz1S+GbkxTCx+HfqSlkrA5UB6bofDMez7XWmRHWKw2FHCXKxIz47fM3JjQTImk9iRy68jUZJrRbPDwSZlhslwYo7s/PCHP8TatWvx5JNP4oYbbsATTzwhztD/8Ic/4Je//GVyVslktF8n1c0EhxM7Q/V7EZEdhXl2Bosd8r1QRZuSIjuEz+MbNn3U6Xaja/0GjAkE0bp5C+Si0GyGNxCAOxCAmaPSDJOxxCx23nrrLbzwwguiqeC3vvUtnHnmmZg8ebKYmfXXv/4V1113XXJWymRkM0G5R4zQqILhIjt0ti8iO/kKieysfaj/R9PRvWjpdMN1LBfmps4B9wnyrLKLHa+Tys97fx5Me19EisSGnBTl5orLVocDZpmjTAzDKCiNRYbSiRMnip+p7FWaLXPGGWdwNRYTNWROJZ+MnH4dgpoKUgpoKLFDPaVCoZAyIzsGLVxeP1zegPhZSUTTWLDN5YJOrYZZP7QYShU2gwEalUqsh2GYzCVmsUNCp7q6Wvw8ffp0/P3vf++P+EizhhgmGr8OVRTJ0UxwMCTahxI7UtM+ZYodHdy+AJxeP8wK8xRRU0HaRhI7FNkpMJtlTWES9B6kdbQ62aTMMJlMzGKHUlc0DJT48Y9/LDw7er0ed911F+65555krJHJUL9OSUnJkBVQShE70qgIk8IiJwQJnFC4G3aXVwgfpXGqxoKS2FECRRYLWp1OuZfBMIySPDskaiQuuOAC0aOE+uuQb4emoTNMtGJnwoQJUAIkdijSlF6RnV4B1uEksaM8MdbfWNCIIb1QHW4XpisgqieJnX2NDcJHpqRO2QzDJI5RD/whYzJtDBMtXq9XlB3LbU6ONrKjtD47BEVzwt3dcHmU59mRxI6zwzmk2HFSf65QWDGRnUKLRUTJOj0exayJYRgZxM7vfve7qB/w3//930ezHiZLojqE3A3lJKT5UmRIjpzRRZEdapKpVk5Vdz8kcALBMAKhkGLTWC21LcAQBU4dfRGzPIX0BqLIDkGpLBY7DJPFYueRRx6J6sHIbMhih4lG7FBvGKU0oSSxQzidzgEDIcWoCPGF3PvlLDe/t+8ckAra5muB2x3Ay4HDsNgHVTXlnQm5Izt+r19En1SDUkN2jwc5NHHeYIASMOl0Ymsj345CUmsMw8ggdqTqK4YZLY98cBCfr96CMTkqPLr6kKzrkHA7OrHxSDtCb29HQdmJUvjdGw7B1eXEv89Rno+DTixyVDkIh7qhM4w6G52cXjs9gMvng9VoPCmyYzUaThJBclJooYosNikzTKainE8bJiugiIS9tRl5Rco5g9Ybe1MXPs/ACI7f54VGJ28fmJHIGQOEwz3QaJWXZ6NhoJI/ZzB2jxs2ozJSWBJFlly0uVjsMEymoo43DfHmm2+itrb2pEGgv/3tbxO1NiYDcdk7EAoGYC0sgVJQqTXQ6HTwewY2lgv6fTBbKdU28D2uFMaMyRERHrl71YzUWHAosdPp9qCqoABKoshsxufV1cIDxTBM5hGz2FmzZg2uuOIK0VyQys5nz56Nmpoacca+cOHC5KySyRi62prEl7O1QDliR4ru+AaJnYDPA01xOa0aikR5Gqcfaiqo1qhPEjtU3k2enfmVlVAS0tgI8u1UyL0YhmESTsxprHvvvRc/+MEPsGvXLtFM8B//+IfoUXL22WfjK1/5SuJXyGQU9rYmmG35UCugmeCpxY4PWr0yTLRD0UOmGAUP66boDnl2IiHxQ6blvEE+Hrnp7ebcO8aCYZjMI2axs2/fPtx44439c4WoZwpNYH7ggQfwq1/9Ku6F0MR0OuO/8847+2/z+Xy47bbbUFBQIJ7j2muvRXNz84Dfo1TapZdeKqp7iouLRRdnmmfEKJOu1ibYFOTXGSB23Ce+6MLhkEi3KVrsCK3Tg+5uZSoeMik7fDQM9ATUy0ZJZecSNB8rz2jisREMk6HELHaoFFfy6dBcoyNHjvTfR43i4mHLli34wx/+cFIHZurWTDO3Xn75Zaxbtw4NDQ245ppr+u8Ph8NC6NB6NmzYgOeffx7PPfcc7r///rjWwSQXEq+urk7YikqhNHSmgZGdYF9EQtFiJ9wDlWoMAn5linvRWHBQZKfT7VZU2XkkPDaCYTKXmMXOaaedhk8//VT8fMkll+D73/8+/vu//xs333yzuC9WXC4XrrvuOvzxj38c0Helq6sLzzzzjDA8n3feeVi0aBGeffZZIWo2bdok9nn//fexd+9e/OUvf8H8+fNx8cUX48EHHxTzugYbpyPx+/2iY27kxiSf+vp64e1Skjk5MrIT8HlFREfy6xAanfK+lAk6jpQOUqtz4PcpU+xQRdZJYsfjEWXnShzLQJ2UW52cxmKYTCTmTxwSH8uWLRM//+xnP8P555+Pv/3tbxg/frwQJ7FCaSqKztCcrUho3hZ1tI28naasV1ZWYuPGjeI6Xc6ZM0cMlJRYuXKlEC979uwZ9jkfeughWK3W/m3cuBO9VZjkQVV8VPVkylVGM8Ghys/9feXnAX/vl7ROoZGdYCCMnJwxUClY7JBnx+3vbSwYGdmhdJESKTSb4Q0E+meiMQyTxdVYVIUVmdJ66qmn4n7yl156Cdu2bRNprME0NTWJidg2m23A7SRs6D5pn0ihI90v3TeSyfruu+/uv07iiAVPasSOrbBUmaXS/b12XDBarAj4e70mSu2zQwKHhI6SIzuUxnL3QAie3L60FfXYGV9YCCUijYpobW3t65zNMEymIFvrVarg+t73vocPPvhAVHWlEpp3RBuT2rQLiR1rkTILeyPFDkEpLZVKJXrwKBHy6ZD3RW/UKjqNRTESSmWR2OktO/cqzpwsQeuiY0pihyLVDMNkcRqLfDU0P2jwRhVTFRUVogSdvDWngtJULS0tojcPVXXRRiZkGjpKP1OEhnw3drt9wO9RNVZpaa/BlS4HV2dJ16V9GGXQ3t4uKvcosqNEqBRerdH2NxYM+rzQ6o2KjEIRfm+vwDFatP0/K3JkRERjQUdf2bnSuidL0PgKEjwkdhiGyXKxQ5VOZC4knw15dmijn+k28t9MnToVt956qzAcjwR5fahXz44dO/q3xYsXC7Oy9DNNoKYmhhIHDhwQpebLly8X1+mSHoNEkwRFimiw48yZM2N9aUwSoaiOaCaoQHOyhMF0ovycIjtKTWFFprEMRo1iIztqrRpatQoOr7ffr0MorcfOYN8Oix2GyTxiTmNRJdbPf/5zfPe73x1wO5WOU3UUNRmkEnKK0PzLv/zLsI9jsVhE9+VIKE9OESLp9ltuuUV4ayhyRALmjjvuEAJHqvq66KKLhKi54YYb8PDDDwufzn333SdEF6eplCd2ioqK4Ncq999FJxoLnjAoK7nsnASOTq+G3qCBzxuEEiFxa9GfaCzY6VFu2Xmkb+donC00GIbJoMjOqlWrTqqckiI1dJ9Ukn706NFRL+6RRx7BZZddJpoJnnXWWSI19eqrr/bfT56Kt99+W1ySCLr++utFw0NqcMgoC/JojR07FkpGH9Frh0rP00Hs0KbUyA5h1p/otWN3e2AzGhVZdh4pdpxOp0i5MgyTxZEdirJQoz9q+BcJ3Ub3EVS6SZGbWPnoo48GXCfjMvXMoW04qqqq8O6778b8XEzqoL5GlGqklgVH2qFYyKTcVn+sP7Kj5JRbr9jRCLETULDYydUb0N43goF67JDYUTKiIisYEA1SuUKTYbJY7PzkJz8Rnpy1a9di6dKl4jYqHSfBIZWhk2+GjMoMQ1Dna6rGEpGd9oGGc6WJHb/Xg+5wuDeyo9CGgpLYySswCrFDPXfC4W6oVMqLmFAaq6YvLUSenYlFRVAy+SYTxnQFhW+HxQ7DZLHYIR8O+WQef/zx/pTStGnTRCXVihUrxHXqqswwkX4d8lCRZwewK9qzQ6LM63Eh6PcrO43lDUFnUItNXPeFYDQpa7gqITw7fp+YV0fTzm0mZUd21CqViFCzSZlhMou4+uycfvrpYmOYWPw6Si3jltD3lUS77R3iUqnVWDT40++XPDuafvGjVLFDA0uburrQ3dOj2LLzSEiUs9hhmMxCeXFvJiObCSrdnEzoDL1fxK6uXrGj1MiOGPzZ09NvUCaUalImgzJR39cvS+meHYLFDsNkHix2mKTS2dkJj8eTFmKHIjlU2SdFdpQqdiRhQ0JHo1VhTM4YxYodiuwQTQ4HKLCXm+Ju6fGKHRpETMZ6hmEyAxY7TFKhqA6RDmKH0mw6owlup13RE88jxQ4NA9XqqPxcmb129BoNNCoVWh0OmHV64YlROr3eMoiKLIZhskjs7Ny5U8y1YZh4/DqFhYUwKLiR3OBUlsfRJaIlGoU2QJTEjrYvhaXkXjskIHMNerQ6nYpuJhgJvV9p3ZzKYpgsEzsLFizoP8uhqec054hhoiFd/DoSFNnxuLpE2blSDdUkbDRadX+puZLFjuTb6XC7YU0Dvw5BY2psNhuLHYbJNrFDf/jV1dXi55qaGo7yMFFBg1xpMGtaiR2DCV6XQ7F+nciycwn6WanDQCXfDvXYSZfIDsEmZYbJwtJzGtdATQLLysrE2S4N6SQj51AkYkwEkxk0NjYKYZxOzdmosSANA1V6Q0GpCougn9ube2d6KRGDRgOHz5s2kR1J7Ozdu1fuZTAMk0qx8/TTT+Oaa67B4cOH8e///u+isWA84yCY7PPraLXafsNnuqSxAn4vVJre/jVKhMzINABUgnrtkACiMn8lpt7GjMlBIBSCJY2G89J71m63i+gkvYcZhsmSpoJf+tKXxOXWrVvxve99j8UO088jHxwc8vbta7cgFFThf9ccRrqgN5jQHQ4pUjRIkLCx5hkGRHa6w90IBbtFKbry6BGNBbVpUIkVKXZIPJI/kSLaDMNkWen5s88+2y90yHwqlRYzTCT0RWFva4K1qBTpBEV2wqEgehTsSwsMkcYilGpSDnd3I2fMGNFBOV2giiyCfTsMk6VihzwYDzzwAKxWq5g4ThsZmB988EE2LjP9+NxOMVTTlmZiR6s3IhwKQalfyzTwkwZ/DhA7/fOxlNlrJxTuhk6thjsQQLpAs9zoM47FDsNk6Wys//zP/8QzzzyDX/7yl/3zsT799FP813/9F3w+H/77v/87Getk0ozO1kZxaS0oQbpFpHJyVIqN7PQ3FDQMEdlRaEWWJxCAUaeD0+dDOsEVWQyTxWLn+eefx5/+9CdcccUV/bfNnTsXFRUV+Ld/+zcWO4ygq7UJRosVOkP6VOAQwT5zcnd3GEpEEjTSAFBCrc6BSp2j2DQWVWLlGY1pKXYOHhzaj8YwTIansTo6OjB9+vSTbqfb6D6GIextzWmXwiICPi/UGq3w7Sh9VMSAMRcKbizY5fEg32xOS7FDn2mhkDKPK8MwSRQ78+bNw+OPP37S7XQb3ccwJBQcHS1pLnaU+QVHgqZ3lMXAyialih1fMAhfMITCNBU7UkUWwzBZlsZ6+OGHcemll2L16tVYvny5uG3jxo2ip8q7776bjDUyaYajoxU93T2wFaav2AkFlWmmJROyVts7ADQSqdeOEqM6RHFuLpodXUgnpIqslpYWlJSkl/eMYZhRRnaokzLlsa+++mrRdIs2ajh44MABnHnmmbE+HJOB2FuboFKrYc7r/bJIJ4J+nxgZEfB6xFm9IrsnR5iTlR7Z6fJ6xWWp1QqXz59WFZs0vNZsNvP0c4bJxsgOUV5ezkZkZli6qL9OQTFycmLW0rLj93lgsOSKL2WK8ijNYD14VITS52PZPR5oVCoUWSyizw5VZtFg0HSBK7IYJjNIv28jRvnNBFubYE3DFJYU2TGYc8XPfq9b8Q0FB0R2/CF0dysrGuXw0kwsAyx9Q0DTzbdDqSyO7DBM+sNih0koPo8LPo87Lc3JBEVzTBabYsXOsJEduq2nBwF/SHGRHavBKCafp6PYocgOGZTTKf3GMMzJsNhhEt5fh0hnsWOwWEU5t9/jVlzUbDixo9Upc2QEeXZsRoOYfK5W5YhIT7pFdsLhMDo7O+VeCsMwqRI79GFbW1srOiUzzHCdk42WXGHyTdc0lt5ghFZvEFEqJREKdSMc6h7Ws6M0sUOfFyRucg1GIR4pupOOkR2CfTsMk0UGZfrwmjx5Mvbs2YMpU6Ykb1VM2tLV1py2fh16fwd8Hmj0BuiNZsVFdoYaFaHkYaB+rx/BcFhEdggSO640EDutjz0+4D0R2roNh+xdKJg06aR9i+64PcWrYxgm6ZEdqq4hkcNNtpihCIdDcLS3pGV/HakZIqUsKKpD0899CvPsDDUqQkKlyhGNBpUkdryO3pQVRXYIi96QdpEdikgVmM1odykryscwTJI9OzQA9J577sHu3btj/VUmw3G0twojp604ff06hFZnEGk4xUZ2hkhjEVq9ssrPPc7ehoK2vkqsdExjEflmE4sdhsm2Pjs33ngjPB6PGA2h1WpF461IeD5WdvfXoWaCFlv6NRMkAv7eL2KtTt8b2VGYZ4fETo4qB2rN0OcoSmss6HF4YNZooNNoBogdSg1RxCRdKDBbcLCpKe3WzTDMKMTOo48+GuuvMFmCXWomqBo4tyldIL8OodUbhWeHIj3d4XCcrTeTMyqCBM1wX7i9IyOUM8DU6/RirPFEU0YSO9RY0E2NBXU6pAs01ysQCguhljvo5I5hmPQg5o/xm266KTkrYdIee0sTyiZMRbpClViERq/vryajjspQyPfycGXnEnSfw+5VlNixRoiD/l47Xm9aiZ18U+97gVJZLHYYJj2Jq8/OkSNHcN999+Eb3/iGGJJH/POf/xRVWkx24nNTM0FX2vbXkSI7NARUpVJDb+wTOwry7ZAfZ6hKLAm6L+ALQ0meHWtEZEcaE5Fuvh0SbNQjiH07DJNFYmfdunWYM2cOPvvsM7z66qtw9X0AfPHFF/jpT3+ajDUyaYC9tVFcWgvTdzp0wOcTlViEFNlRkm8nmshOMBBCOCx/t18yqg+O7Bi1WiEa0k3sUBVqntGEDrdyhC/DMElOY/34xz/Gz3/+c9x9992wWCz9t5933nl4/PET/SmY7PPrGMwW4XVJV8ijI4kdjU4PlUqlrMhOFGJHmp8ld7LF5+41IkdGdtK1saDk2xltZGfzW0ej3rfL0YmyeTpxMqnRaERzw7KyMvEzwzApEDu7du3Ciy++eNLtxcXFPDAvi6Hhn+naXycyjSWJHfpipoospczHEg0P/SR2NKcUOz6v/GKHojpEZGSHSFexk28241hHe9L/jWuOH8EXez9Ha3szKg7lQa1Wi95PdB9Vv1JUfcWKFSgoKEjqWhgG2S52bDYbGhsbMWHChAG3b9++HRUVFYlcG5NmzQTLxqevOVkqPbfknfgSoVQWDTVVAh5fED3dPSNGdvSGXiGkhIosKjsfSuyQb8fpTT+xQ40FPf4APIGASMclGp/fh3Wb3kdtfTUqSsfh/DMuwZU3nyVae1BKsLm5GQcOHMDWrVvFZ+3SpUtx/vnnc6SHYZIldr7+9a/jRz/6EV5++WVx9kt/iOvXr8cPfvAD0YOHyT6cHW3ifWBNY3Nyf2RHN7b/uojsKMSz4/L6xeVIYkerUyEnZ4yI7EABkR2dUQf1oDYEFNlp6LQj3ZAqsjpcLhjz8xP62F1OO1ate0MIngvPugxVFRPFZ6vUw4w8Q5TCou30008XfknyTh4+fBhf+cpXUFKSvj45hlGsQfkXv/gFpk+fjnHjxol88syZM3HWWWeJ0CpVaDHZ6dchf0tuXno2ExzKs0OQ/0gpIyNc3oC4HKkaS6TeDNRrJ6SIyI4x94RfRyKysWC6iR1qb5ToiiyXx4l3P3xN/HzlRV/D+LGTRmxcSJGcM844A9/5zndEiuuZZ54RER+GYRIc2aG88R//+Ef85Cc/ESMjSPAsWLCAB4NmeSVWbho3E+yfixUKDRA7J0ZGyO2AAVy+PrGjG/lPliI/Pq8C0lhOD4wWI7Y0bRlwe52zC3WOOnxStxEG7cAUzJLSJVAqFKGiiqxEip1gMID3PnqD3Dq45LxrYDaeKPg4FWRYvvnmm/Haa6/hpZdewpVXXon58+cnbG0Mk2nE3Ru2srJSRHeIeFuoP/nkk2KrqakR12fNmoX7778fF198sbju8/nw/e9/X/wx+/1+rFy5Er///e8HhG1ra2tx6623Yu3atTCbzaLp4UMPPSTOepjU0NXajNKqyUhnpLlYGl2E2DGaEAoG4A+EoNOqZY/s0KBPlXrkYCz5dpQwH4siO4UVhUBHb/pNwtRnsPb4gyeJnbSYkZWg8nOKbH26ZS1cLgeuXPm1IYVONNVbVfpFOB5y4PXXXxePSSeeDMOcTFyf4BQ6feSRR3Do0CFxnaI6d955J7797W/H9Dhjx44Vg0Xp9+kP9fnnnxdnKGTAI+Fz11134Z133hH+IKvVittvvx3XXHON8AgRVKVw6aWXorS0FBs2bBDGafINUaiX0m1M8nE4HPC6nbAVlyGdOK326QHXWzudONK1C8vbgygM95bPN3Z2oaNrD5ze8xQgdvxi0OepoDRXR6u8qbdwKAy/xw+DhYRj54D7jNoTYqcg+kCGYkzK+xoaEvJYh2r243DNfpyzfCXyrPFXVpGf54wl56E1Zx/efPNNGI1GTJs2LSFrZJis9uxQ5OV73/seLr/8ciFCaKOfSZjQfbFAv3fJJZcIsTN16lT893//t4jObNq0CV1dXUJU/fa3vxU9fBYtWoRnn31WiBq6n3j//fexd+9e/OUvfxEhXIoIPfjgg3jiiScQCPSG/Znkcvz48bRvJkj4Ar3REIPuRLTBpO+tunF6BkYn5IrsjGROltAb5J98LpWdD+XZ0WnU4gvaE5A/1RYr+SYzHF4fAqHRHV+vz4NN2z7G5PHTMWXC9FGviyLr9DlKXkr6PK6rqxv1YzIMsl3sUNqJPDuUKrriiivERj8//fTTIsUULxSloXSV2+3G8uXLRYllMBjEBRdc0L8P/TFT+mzjxo3iOl1S34nItBaluijaMNLoCkqJ0T6RGxO/2DGY0ruZIOH193756iJSK0bFiZ1Tp33IoEz9eCi6Iqdfh+iN7Jz8xUzRHSrjTsfGgsRoOyl/tv0TcXnawjORKEhAXnvttaL9B/VBa21tTdhjM0xWih0SIIsXLz7pdoq8hOI446EmhRTN0el0+O53vysMd1Th1dTUJMzQ1NcnEhI2dB9Bl4PLLqXr0j5DQeKM0mLSJnmPmNihs8h0nocl4Q8EodWoxDgDCbVaJW5z9lVCyZ3GijayQ/j7StXliuyQqNGbemdhDcaoI7GTjpGdEwNB46W5rRGHqvdj6fzTYdCfHPkaDeRTpHmF1NmeThzJ88gwTJxi54YbbhDRncFQZOe6666L9eFEfnnHjh2idwQZjclgTKmpZHLvvfeKNJm0cdg3/mgc+aTSvb+OFNnRD2GYpVSWw+1TRmRnhLJzCX1f9Mfv9ssa2aGoDkUbMkns6DQaUTofr9ghX+KWL9Yjz1aAqRNmIhno9Xp87WtfExFyybTMMEyUBmWagyVBZ2x/+tOfhF/mtNNOE7eRUKGqqHiaClL0ZvLkyf3RoS1btuB///d/xR8s+W7sdvuA6A51EiVDMkGXmzdvHvB4dL9033BQFIk2ZnRQ9Iyieek+JoLwBkjsnPznYDTo4JQxSkKEwmEhxqJLY/WNjJBRoFEllsE8fLk+iZ22vg7L6QZVZHW44xM7xxuPobG5HivPvmJYIZgIaJQEFXJQOuvTTz/FmWcmLl3GMBktdqg6KhISJcSRI0fEZWFhodhG8slEC3XiJU8NPQdVVa1Zs0bkoglqnkWiijw9BF2SqbmlpUXM5iI++OAD5ObmilQYk1woIkahc0N+EdIdKi/XR5iTJYx6DZweedNY7r6+OdGksag8PUeVA5/HJ2saK7cgd9j7jTqtMChT1CHethVy+nZq4pgBSK91++7NKCksw7jy8Ug2VPBx9tln48MPP0R5eTkmTZqU9OdkmLQXO9TDJlnpJKqgItOx0+kUZyIfffQRVq1aJbw0t9xyi4gq5efnCwFzxx13CIEjRZQuuugiIWootfbwww+LSAN1cb7ttts4cpMiczK1sHekcTNBCYqc5A9RPURpLKent7pI9lERp2goSAivDFVkyZnGcnhQMn746jwyKIuTmmBoyNSh0iuyth+rRbi7G6oYojPHjh0Tfh2K6qRK4J1zzjnib5R8kGQRMPV5jhgmG0leLDUKKCJDqS/y7dBQO0phkdC58MILxf3Uy+eyyy4TkR0aSUGpqVdffbX/92lEwdtvvy0uSQRdf/314vEeeOABGV9V9kAfpNQrKRPwDePZoSgEpbHk9D5EMypicEWWXJGdYCCIoD84ZNl5ZBqLSMfyc+q1093Tg84YK7IonZRvK0xJVEeCRNXVV18t3rtvvPEG+3eYrCbmTmnk8H/sscdEtIfECp2hRbJt27aoH4v66JzKbEc9c2gbjqqqKrz77rtRPyeTGCgSR34qqmTbW4+0xzdcGsugRTjcLaaOmwyJn3YdDW5fQHxxaaOI7BByRnakaec0KuJUYsftCyLfnH5ihyCTcqEluq6IVAZOQzvnzliatKjOSN2WJ1gWYtVbb8LfaMCsqfOw9PKJSVkDw2SU2KHUEpmTv/zlL2Pp0uT98TLp0UxQRHbqhy/zTwdCoTCCofCQ4wv6Gwt6/bKJHYrskECgiebRQEZmp0yRHamhoCF3eIOyXqMWKSB3GvbaMWq10Gs0MfXaoSIKaq8xsVKe+YGVFRMwa9o80d+nrLgCAIsdJvuIWexQ2ogiKaeffnpyVsSkjTmZfFXkpQKa0r4SixiyGiuisWBpvkU2sWM2aBFtEkJEdprki+yo1CroDMN75kRjQZ0mLcUOrZ1Mym1Rlp9TJPyLL74QaXaVSz5v29L5Z6Ch+Tg+2vQ+Lvx/i0Tqn2GyiZg9O9Shk5pWMdlNRvl1+kZFDJXGMui04gvOIWMXZZfPD/MI4mGoYaDknQkHw7JEdqjHzqkivma9VqQG05F8sznq8vOdO3eK9gxDNWJNJWqVGmcvuxAd9nbhH2KYbCNmsfOb3/wGP/rRj0R1AZO9zQQbGhoyR+z0NbiLnIslQakjiqrIOTKCSt9pDdHS32tHhlSW1FDwVJDx25WGkR3Jt9Phckdl+CUPIxVgKOEEsaigBPNnLMK6detG7DDPMJlIzGKHzlAoNDtx4kTxB0xl4ZEbkz3NBDNlzIYvcPJcrEgsRp3MYscv1hBLZEeuxoIU2RnJnBwZ2SHjdTpSYDIhGA7D4R25JQF1F6e/lQULFkApLJi9DEVFRaK7Mp20MEy2ELNnh2av1NfX4xe/+IWYQ8UG5exMYVHOf6Qu1ekW2Rk8FysSi4xdlCl6QH12esWOJ4bIjh/+FAs0WitFdsZOPXXEjzw71GeHRIMmzfwjkRVZvb3fh4bG4JAxWeoQrwTo7/aqq64Sw5w/+eQT0YuHYbKBmMXOhg0bxLTxefPmJWdFjOKhLtbUlZW6J2cCYhTDCM3tSGgcb+2CXH6iULhbREKiRa3OESbhVEd2Ar6A8AmN1GNHQno95NuxmtJL7OQaDEKgtY9QkUVRE/LrLFy4MKmjIeKBGoHSCImPP/4Y06dPz5iTFoYZiZj/CumPw3uK8C2TudDZO1ViUdfrTIEEhWGISiyJXBnTWFL35FjSWBRt1Rl1KY/s9JedR+nZIdK1IotmZI00EPTo0aPic3Lu3LlQItSklUbsUHdlTmcx2UDMYueXv/wlvv/974uxDu3t7XA4HAM2JrOhKfH075wpfh3JszNUJZYECQ3yl1BzwVQjzeWKxaBM6E36lEd2+hsKRhHZoTQWiYZ0FDvS2IiRxA7NCaR5gdLMPqVB6awrr7xSNDykaD3DZDox5yG+9KUviUsa7xCJNNSPzxIyG4rqEJkkdiiNlTeCqVaKqlC/G6tZL19kJ4Zh2xTZSXU1Fvl1NFoNNCMIR4mcMWNE9Rt1UU5HqNdOdWvrkMNM6TNw//79im+6SumsFStWiOosmjFI09IZJlOJWewkaygokz5ihz4UM2mooEhjjRTZ6YuqODy+lIsdpzcAnVYNjVoVc2TH0ZbaSKvX4R2xc/JgzDpt2kZ2yKTsCwbhdruFCTmSI0eOiIrVWbNmQenQZPS9e/firbfewk033aRoccYwKRU79MfBZLc5OZP8Ov1DQEcQO7mmXoEjh29HVGLF0FBQQm/UyxLZiabsPDKV5UrT8vP8PrHf1tZ2kthRegorEo1GI4Ytv/DCC9i+fbswVDNMJhKz2CEH/6mMb0xm4vf70dzcLMLzmYI0F2uoURESdB+VpVOURek9diR0Jp2ojKJOypRaSpVnp3RC9JU9VJHV0hXb9HClkGcyiVQceV7Gjz8xyZz6T1EK67TTTkubKAn1TKNeQDTzcMqUKYpogMgwsoudofoyRP5Rs2cns/vrkEchk/w60qiIoYaARr6/5WosSAblXFMckZ2+aBRNP0+F2Onu7obP5YuqEityyCrNJQt3d4vBoOkErTfPZBSRncFVWHRSkA4prEguuugiHDx4EP/85z/x1a9+Ve7lMEzCifkTprOzc8DW0tKC9957D0uWLBFnBkxm+3UMBoMI0WcK/UNAT2GqJbHjkKEjcbxpLDIoE6mqyCJRRYInljSWSaftbUTYN64jHX07FNkZnMKiDsXpkMKKhP6uL774YuHfocgUwyDbIzs06XowF154IbRaLe6++25s3bo1UWtjFOjXoahOuoTno8Hf90WrP0X0Q/TakSON1TfxPFYMpt4Ii9ftTZlfRzxvLJGdPoFJJuV4BJ0Sys9rIsSOVIVFE87TEYpGUSPEd999FxMmTIBOl37/JgwzHAmLHdPoiAMHDiTq4RiFQWftlMbKNHOyFFWIJrKT6jRWIBgSWzyeHZVGBa1eK1JLqRQ7MUV2+roou9M0slNoscDpdIrKK4KGI1MKiwZ/piN0EnPppZeK17N69Wq5l8Mw8kZ2SPlHQmFoGnhHzQbnz5+fyLUxCoLSlYFAIKP8OoTXHxhxLpac87HibSgooTfrUyZ2qHsypc5IZMXie6GIWroOBKVeO/W7m/H+i5+jtKgcG7euQ2dtALVbPKgbcxTpCEXuqYcaeXfmzJmTcSc3TPYSs9ghQUNnACRyIqHqgz//+c+JXBujsBQWdV2lmViZhNcfgjGKuVMUXfEHQmKjvjdKHRUxOJWVsjSWI7ay88hUljuNy8/H5OSg094mxE5dYw0qKyakfZqX/Je7du0SvXf+9V//NWNm4DHZTczv4urq6gHXacgdGfL0+tQ2W2NSb06mjqvUlyPTIjsjNRSM9OwQFN1JldiRIkmjiex0NHYgVZGdWPw6kamsdE1jqVUqWM1WdDo60OXoRJfDjqXzz0C6Q5/pl19+Of7whz/g008/5cnoTEYQ86d2VVVVclbCKD6yQy3lMw2vLxiV2JGiK+TbKbSmpns0jaegzsk6jVrxaSyK7OSX5cf8e1SR1eGUZ6J8IsizFaDT3o7ahhohEipKMiPNSx7MM844A5988on4u0+36jKGGUxcn6Jr1qwRG/k4yLgaCaeyMg8a/EkDQDMxf0+l53lRDK6MFDupbigYb1qE0ljUVDDZjQXDobDo1hxvZMcTCJ6UFk8XbLn52Hd4l/g3KisZC40mviicEqEGsVSK/uabb+Lmm28WYo5hskbs/OxnP8MDDzyAxYsXi7RGuuenmeiiOkSmmZMJry+6NBZFWMhM60ix2Ik3hSVFdgiK7mjykyd2vC4v0BPdtPOh5mPRCZPU7yjdyLcVwO12wuNx4fQl5yEd2PxW9ObpSuN8bKtdjS1btmDZsmVJXRfDKErsPPXUU3juuedwww03JGdFjCL9Ovn5+SfNAEp3qHPvqYaARkKdjFMZ2aG5UaPpP2MwG/obC1ryLUn16xDxiB2aj0Wkq28nz1oAh8shTLzjyk+MjcgUyHi9pGSJiORTSb3NZpN7SQwTFzHHJan8eMWKFfE9G5O2YicTozr9oyKiFDskPNIqskMjI8b0RV6S7NehCK80oiIWaD4Wka4DQXPNNjjdXVCr1LBaMlMIXHDBBaIA5Z133knbdCPDxCx2vv3tb+PFF19MzmoYxUHitqmpKTP9On1fsIYoSs8Jq0mPrhSOjCCDcrxl50SOKgc6gy7pIyOooSD5deLxdPSmB9X9ZfbpBr1mf8APkzFzh2dSJ2WajH7o0CFRks4wWZHGou6aTz/9tOiwOXfu3JNKkX/7298mcn2MzNTX1wtPRUZGdvpSJyMNAY3Eatbh4PGBgx+TRSgchtcfHFVkR0plSWkmpfXYkTDrdWkb2emwt4v+U5nWkmEwU6dOxezZs8UcxEmTJsFkSk1FIsPI2kFZ6pS8e/fuAfexWTkzeOSDg/0/H965GcfqXPjLjk6MGWNHNo6KkMg16kUEgoQI9VhJdlSHGO3MKEotJT2y4/DAVhR/CodSWekqdhqaa2E2mMUJAaV4MvkzkAaFPv7441i1ahWuueYauZfDMMkVO2vXro31V5g0xt7SCFtRaUZ+iFNDQWoQeKpREZFpLMLh9iM/DjNuLEhG6NGksaSKrLYkRqPE1PIuD8omlY1K7LR0uZGO1DfVYWx5lUj3erxumIyZZeIfXL1VYZiNf778Prpbragsn3DSvksvn5ji1TFMdHDjBGZY6GzV3toEW1H8X2RKHxURrTl5gNhJgUlZiuyMOo1FIyNc3qQZS4P+IELBUFyVWJFix+MPIBgOI92q+RpbjmPK+Oniut2Rmm7VckKvtaKsEuu3rEUgmJ4+KyY7YbHDDIu7qwOhYAB5xZkqdqLrsRNZek6kwqRMkR2VKiem9Q0X2QkHwwj1VZ4lI4VFjMazI6Xqujy9j5UuNNjtCIVCmDx+hvDtkH8n06EI7xlLzoPP78XnOzfKvRyGiRoWO8ywdLY0YEzOGFgLS5CJ9I6K0MZUOURDQ7tSMIKB5mJRxGO06UOp106yys+pEosYbWSH6PKmZmhpoqhpa4Vep0dRQQls1nx0dmW+2CFyzVYsnrscew/uRHNrg9zLYZio4HG2zLB0tjTCWlAClTozK02iHRURidWkS0lkZ7Rl5xJS7xvqopxbkJuAlQHjX/6s/2d3XTOK6low5c1t8Y+10GqgysmBPc0iO8fa2lFeOU287rzcAti7Mj+NJTFr6nwcOXYQH29eg6u/9A3RZ4hhlAxHdpghIY9HZ3MD8orLkalEOypisG8nFZ6d0TYUlNBR6i2JjQVFl2d9/PO7CPpdmpFl96SPSdkXDKKxy47y0nH9A0E7utqzpuke9Rc6c+kF6HLasWPPFrmXwzCnhOU4MyRelwM+jwt5JZkpdvpHRehPLXY2Hj2RnmhwBtHQ1oWJEbdFsnxiQULWR4JqQmleQr6U9MbklZ+T2JHSUKPBYiCxkz5prLoOEjZARWll/4ysYDAAp9sh0jzZQEFeIRbMWoLte7Zg/NhJKMznyeiMcuHIDjO8X2fMmIytxOofFRHjNHCKQLhT0O2Xyttz4xi/MJxvJ2lixxtISASKGgumUxqLUlh5JmO/sCmwFYnLDntqmk4qhfkzlyDfWoCPNr2PUDg5JniGSQQsdpgh6WyuhyWvAFpdYr5w031UhITJqEMgGIY/mLwP9kAwBF8giNwEeHakiqxkdFHu7umB2x+AKQaT93BQdKjL60mbNNCx9jZUFRT2XzcaTMKsnG1ih6rQzj7tQnQ5OrF99wkvF8MoDRY7zJB0ZLpfp697sjR1O1qkKIa7rw9OMpA8QQkTO0nqokwdqEmcJCKNRY8RCnfD7Vd+7xaXz4d2lxuVBSdSlhQFzbcVZp3YIQryirBw9jJ8sXerGC/DMEqExQ5zEuTV8Ti7kFdSgUwXO9GOipAw6XsFSDJTWZTCIhKZxkpGY0GXr3ediRI7RGcapLJqO3qrrsbl5w+4ncROe2crspG5MxYJ0fP666+L3kMMozRY7DAnQVVYRKY2E4wcFUElz7FAfXZycsYkdUq3w9MbhbEYRy8ipMhOd7gbQV+vwEukOVmqpEqEZ4dIB9/O8Y4OFJjNMOl0J4kdp8shjMrZhpTO6ujo4JFCjCKRVew89NBDWLJkCSwWC4qLi3HVVVfhwIEDJ01Zv+2221BQUACz2Yxrr70Wzc3NA/apra3FpZdeCqPRKB7nnnvu4bOLUZqTTbk26AyZO9k41lEREiR0jDptctNYbr8QEIkaNmqwJKexIJmT6RjGKhiHguaTkeBJhy7KVIk1OKojiR2KnlEJejZCr//cc8/Fhg0bUFdXJ/dyGEY5YmfdunVCyGzatAkffPABgsEgLrroIrjdJ/pt3HXXXXjrrbfw8ssvi/0bGhoGTNwNh8NC6NAgPvoje/755/Hcc8/h/vvvl+lVZYY5OVNLzuMdFRGJ2Zjciizy7EijKRJlUE6K2ElQ2bmEzWCE3atssePy+4VfZyixQ712KNKVjb4diRUrVqCiokKks+jznGGUgqxi57333sM3v/lNzJo1C/PmzRMihaI0W7duFfd3dXXhmWeewW9/+1ucd955WLRoEZ599lkhakggEe+//z727t2Lv/zlL5g/fz4uvvhiPPjgg3jiiSeEABoKv98Ph8MxYGN68Xg8cNo7kFecuX6deEZFRGLU65Kexso1Jq4KTmfQIUeVk/CKLPLsJFLsWI1G2N0exaewiKHEDnURtubmZbXYob5OFKGnz+7Vq1fLvRyGUaZnh/5AiPy+DxISPXR2cMEFF/TvM336dFRWVmLjxt4hdHQ5Z84clJScmN+0cuVKIWD27NkzbPrMarX2b+PG9XZBZXpTgkR+pkd2AsH4IzsGbf9U8uT12ElcZIeiDZTKkuZYJbp7cqLIMyo/skMprHyTCWb90GK016ScvWKHKCwsFJ/Zn332GQ4fPiz3chhGWWKnu7sbd955J04//XTMnj1b3NbU1AStVgubzTZgXxI2dJ+0T6TQke6X7huKe++9VwgraeP88gmOHTsGg8kCgzkxc5QyaVSEhNmgg8frT1pPGJHGSlDZuQRNJU9kZCcY6u01lOjIjscfgF/B6Y+6jo4hozoSUvl5uvQLShbLli3DpEmTRDor0pbAMMh2sUPend27d+Oll15K+nPpdDrk5uYO2JgTYifT/Trd4XDUoyKGwmjQItzd01++nkhC4TA8vkDCys4lKLKTSLHjTGDZuYTNaFT09HPqAdTmdGHsCGKnwFYoqrFcbieyGYomUjqLTmLJc5nt4o+RH0WIndtvvx1vv/22KFkcO3Zs/+2lpaXCd2O32wfsT9VYdJ+0z+DqLOm6tA8THeRlamxszOhmgkTA741rVERkZIdIhkm5v8dOgiM7CRc7fWk86VgkUuwotfxc8utUjhTZyevtqpzNvh0JqrK94oorsH//fmzbtk3u5TBZjqxih9Q+CZ3XXnsNH374ISZMmDDgfjIkazQarFmzpv82Kk0nX8ny5cvFdbrctWsXWlpa+vehyi6K1sycOTOFryb9OX78uPg3yfTITqDPFxLrqAgJ6Qs+Gb6dE92T9QlPYwUDQQQTFI1yev3QqFXQqRNTHk8YtVroNWp0KDTtcbyzQwgyi6G3lH8oTAYzdFk4NmI4yGNJn+NUjNLWxseEyVKxQ6krqqJ68cUXxVkAeWxo8/aFsck8fMstt+Duu+8WUR8yLH/rW98SAue0004T+1CpOomaG264AV988QVWrVqF++67Tzw2pauY2FJYJpMJptzRT9tWMn6fJ65RERLUjJD6wiQnspPYhoKDe+0kyqRMaaxcg06kKxIFPVaeyYQOtwtKpLZ9ZL/OibERBWi3Z2cn5aGgghE6+Xz11VdFqxCGyTqx8+STTwqD8DnnnIOysrL+7W9/+1v/Po888gguu+wy0UzwrLPOEqkp+qOJ7NxJKTC6JBF0/fXX48Ybb8QDDzwg06tKb7FTVVWV0C8wJeL3ekAvMdZRERKia3CSKrIoskPGaa1GndDH7W8smKBUFkV2LAlMYUnkm8zocCkvsuMNBNDmco7o15HI1hlZw0FFJvT5TSeyH330kdzLYbKUxH6ixkg0pjW9Xi965tA2HPQF/e677yZ4ddkFlfhTGosiZfUZ3nbI73GPuvMvpbJcfWMdEl52nmC/TjJ67ZBnp8RqQaKhsu4jLS3is0FJoptSWPRxNZJfR6LAVoS9B3ciEPRDq+HoMlFeXi56pZElgaq0xo8fL/eSmCxDVrHDKAcqv6cQM/mmNnwx0BCeafi9buTF2VBQwmLUo61rdBGI39t3nnTb523H0I0e9Axxn1J67VDZuS8QTIooI7HjCwbhCQROmj0lJ3XtHbAaDKI8/lQU5pcIsUb9dsoyvDlnrN2Vqe/OP/7xD3z3u98VKXOGyapqLEZ+qqurxYdPUVERMh2K7NBAz9FgNurgTEJkx+cJwmCML70WVa8dhzdhZeeWBJadS+SbzeKyU0aT8u6WopO2z44F0a0af9LtQ2HLzRNp9baOE0UTTG93ZUpnUTk6FaVwOTqTSljsMP1ih6I6SkodJAufl8TO6ASFxaiDPxBCIJjYgbM+Twj6ZImdXGNCIjvJKDuP7KJMb0ElVWQFQgF0uO0otUV3IkBCp7eTMoudwVAhytVXXy0iPOvXr5d7OUwWwWKHEf11aMBqtuTRExXZIVx9peKJoDvcDb8/BL0heZEdj8Mz6jPqZJSdS9Ckd6vBiHaXciqymrt6OyKXWKOPehbmF6OVIztDMnnyZJx55pmi3Yg0noZhkg17dhhRhUWh5cF9jjIRep1+nxvG4tF7dqQv/nxrYrwHPl+IXPvQG5PzZ2m0GhGmMQ9eP/Sj6OOTjLLzSKj8vNOjnMhOk70VJp0RFr0pJrGz//Bu0U1Zo0l8uk+pbH7raFT7mbqrMG5cLV555RXh3zFG4YVimNHAkR0GNTU1og+GNIA1kwn6vejp7hl1ZEdUc6ly4Ozri5Movw6RtMhObu8XiqfLo8iyc4kC6rWjoPJziuyU2IpiEneFecUiGsQl6CP7d0KhEPt3mJTAkR0mq/w6lMIiRit26FiRb8eVwMaC/WIniZ4dglJZ+WX5iis7jzQpbzt2DOHu7lG1B0gEgVAQ7a5OTC0bOurp3rxlyNu13d0INzai7tN1MFdMFreZli5J6lrTDTrBIv8ONZX95JNPRB81hkkWHNnJcjwej2j2lQ0pLMmcnAixIxl0nQn07Pi8Iag1Kmg0iffCEGqNGjqjblQmZRo5kayyc4k8kxHdPT2KmJHV4ojdr0OQSMszWYVQYoZnypQpOPvss0WHfDItM0yy4MhOlkN+HfowzyZzsug5E+cQ0EgsJj1aOhI33drrCSYtqjPApDyKNJb0u8koO5coMPWWn9PYiIK+UnS5aLa3waDVI9cQ+zryzXlodbYnZV2Z5O8x9IzDGEcuHv35H3HVyq8j12w9ad+ll0+UYYVMJsGRnSyHUljk1bHZbMgGqKGgzmBETs7oU3bkW6FeO4nyG/i8QegNyT3/IJMypbHixe1wJ63sXIKaCWrVKnS45Y/sNHW1ipLzeFK8BWYbutxOhMKJbU+QadCxPWf5Smg1Wqz+9B2EQokZVsswkbDYyXIkv062QJEdnSEx1VNUfh4IhhPWa8eXisjOKHvtUGSH5nbRdPJkIYZpihlZ8pafB8MhtDs7Y05hRUZ2unu60enO8PkrCUCv0+PCMy9Dl6MT6z9fy4ZlJuGw2MliXC4XWltbsyaFJXl2dMbEiB1KY0nDOxOBxx2A0aRNehrL7/EjFKdAc3e5k5rCksg3m2Ttoky0dLUJsVIat9ixImdMDvt2oqQgrwhnLDkfB4/uw77Du+ReDpNhsNjJ8qgOkXWRHX1ienpIJl1HAsrPw6FuBHzJ654sYerrCRRvKkuInSSmsCTyjSa0yyx2mrtaodfqYTXGV3mmylHBZsplsRMDUyZMx+xp87Fh6zo0NNXJvRwmg2Cxk+X9dWgWlllmE2gq8Xmc0JsSUzat02qg06oT0muHzMmE0ZTkNJZ1dL12SCSZ9SkQO2YzvIGAGAgqZ3+dUmvhqFoyFJjz0OHM7MG6iWbZgjNRXjIWq9e/iy4+dkyCYLGTxWSbXyccDiHg80HfV+2TCHJN+oREdiiFRRiSnMbS6rXQaDVwdbniKjunFFgyy84juygTcqWyyFTc5uyI268jUWDJQ6e7i03KMTYcPG/FxcLH8/7HbyEQTFx7ByZ7YbGTpXR1daGjoyOrxI7P3fsFrzcmTuzQ2IhETD8nc/KYnDFJ654sQVEKk80Etz12ESH9TkrSWCaTGAgq14ysVkeHaGpYYi0c1eMUWfKF76fdxRGKWCChc9FZl8PjdePD9e+JMS8MMxpY7GRxVIe++KqqqpAt+D19YkehkR0SOokoiY/Gt0Pem3jFjjUFYkejUiHPaEKbM3F9jGItOddptKIx4GjIM1uhVqnR6uB+O7Fiy83HeSu+hOONx7B69Wq5l8OkOdxUMEs5cuQIysrKsmoAn7cvsiNKz+2JEztubwChUBjqUUwB97qDMCTSnFz9ybB3mbwtaK1uB6rDvTdMODOqh3R2kt9JLyaep4JCixmtLpnEjr1VVGGNdoQKVWMVWvJESoyJnXHl44WHZ8OGDSgoKMCiRYvkXhKTpnBkJwuhHhZHjx7FpEmTkE1QZEej00GdwCnUUvk5DcccrUE52X4dCZNZKyq/Av7YfCSuThfMttSZ2QvNFrQ5XTL5ddpFM8FEUGQpEGkxJj6oOmvp0qV45513cOjQIbmXw6QpLHayEJqF5Xa7s07s+DyuhPp1pMgOMdpUltcdgCHJlVgSJktvGsrjiq3SyWV3wZyXQrFjscDt96e8IuuEXydBYic3H26/By63PFGqdIeia1/60pfEHK2XX34ZjY2Nci+JSUNY7GRpCkur1WLs2LHIJpIhdow6DdSqnFGJne5wtxgCmtA01ggYLb0RJLczENMaqVw9pWKnryVCqn07jfYW6DW6Uft1JIpyC8RlSzt/SY+mQuvaa69FYWGhmJJOBRYMEwssdrJU7FDXZLU6uyxbPnfieuxEnnVSKms0YoeEDnp6kt49WYKmqusMGrid0afeyNBM6c9Uih0qP6fp4W0p9u00j2Ie1lDQIFGz3oSWtqaEPF62Qido/+///T+oVCr89a9/hc83+sIAJntgsZNlBAIB1NbWZl0KK1mRnURUZEk9dvQpSmMRJosWrhjEDvl1iFSKHRI6VIKeSt9OMBxOSH+dwRTnFrDYSQDUAPW6666Dw+HASy+9hFCI+xcx0cFiJ8s4duwYwuFw1omdcCiY8IaCEjazAQ6XN+7fF96ZMWNgTFEaizBbdDGlsUjsUENC2lIJ+XZSmcaq7+wUfp1EmZMjfTttnS3cXDABUNd3ivAcP34cr7zyCvfgYaIiu/IYWcwjHxwUl/u2fIyWpgBe2NaOMWOyp0LE5+ntEZOUyI5ZD5fXL8rP4xU75NfJUaXu3MOcq8PxGju6u3uiOuMhczI1I0xUamcotjRtOem29mAz9jS2YHKjasBzLyldkpQ11HW0C7+OzZib0McttBQg3NKKDnsbigtKE/rY2UhlZSW++tWviujOW2+9hSuuuCKp700m/eHITpbR3liHgvJxWffBQH4dIhmRHavJQJabuFNZYtq5ObURE7NVJ0zHVAUWbWTHkpdYv1M05JkMCIbCcPlSU5FV196RUL+ORL7ZJrwmza1sUk4UU6dOxVVXXYXt27dz00HmlHBkJ4ugcQkuewcmz12KbPTrJCuyQ2ksoitOsUPpJFt+72OkCnNub8m8y+FH7xSq4aE0AUV2KqZWINXk9R3bTrc36WMqAqEQGru6UFU6Kyn+o8K8YrS0s28nHja/dXSYe8yoNM/F3//8Nqq3dWL+zMXi1qWXT0zp+hjlw5GdLKKtsVacseaXjUO24XU5oDMYoVIn3hej1/VOP7fH6dsRPXZSHNnRG9RQa1RwOk5tUva5fCIKlEpzsoRBq4Feq0any5c6v06CzckSxYWlbFJOUtPBBbOXYsuO9dhz8Au5l8MoFI7sZBFt9TWwFpZAq+s9q88mvC4nDObE+jASYVKmLsbBQDhlZecSJHrJt+Pq8kc1JoKQI40lpbIospNs6jo6YNRpYTUm53WWFldg1/7tcLocsCTxvZiNLJpzGkKhIDZ8/pHoybMUHNlhBsKRnSyhm0pqG+pQVDEe2Uhvj53kRSbIpBxPZMfjDorLVHt2CBI70fTacXe6odKooDfLI5IplZUqsVOZX5A0P1tpUYV47MbW+qQ8fjZDx5VmaM2cOhfrt6zFjh075F4SozBY7GQJnS0NCAUDKKrIninnkXjdDhhMyTubJpNyl8srGu/FgqdPbBhT2GMnUuw4u3ynXDNFdmgmllymdorseHwB+IOh5Pp17HaMy89P2nPodXrk2wrR2Hw8ac+RzdD7c8WiczBt0iy88cYb2Llzp9xLYhQEi50soa3+GPRGEyz5yfEjKBky2HrdlMayJDWNFQiG4fH1RmpiqcTSaNXQ6lKfUaaKrHCoW3hylDQAdCixQyQzunO8sxPdPT2oLOgd7ZAsSovL0djCYieZgueMJedh3rx5eO2117B79265l8QoBBY7WUJrwzEUlldmXck54fe60dPdk1TPjtXS+4Xc2tXbzydaPK4gjObUR3UIi7VvYnvH8E37KOrjtrtlMSdLWIw6qFQ56BhF48Zo+uuYdDrRsTmZlBWPFZ4d2pjkQJ9x1Hdn7ty5+Mc//oEvvmDTMsNiJyuw2+2i5LwwS/06VIlFJDONRSMjVDlj0NI3ViFa3C6/LH4dghoZUkWWo2P4L16/x49gICir2MmhCkKzAR3OJIqddvLr5Cf9ZKCsuLd8n307yYVMyldeeSUWLFiA119/HVu3bpV7SYzMsNjJAg4dOoQxOWNQWFaJbIQqsYhkprGoj0qu2RBzZId67Jgsye0fM+IQU5sezvbhIztS1MeSL08llkS+2YgOlycpj+0PBtHU1YVxSU5hEXqdgX07KRQ8l19+ORYvXiy6LG/evFnuJTEywqXnWcDBgweRX1wBtVaeCIISzMlavSEpPXYiybMY0NIZvdgJBcPwe4NiTpVc5Fr16BghjeVod4hKLGOuEXJSYDHiQH2rMCnrNOqEV2GRX6cqBWKHKCupQG19dUqeK9shQX/JJZdAo9Hg3XffRTAYxOmnny73shgZ4MhOhkN/3NXV1SjM0ioswudyJNWvI5FnMcYU2ZEGcdIEcrmgyA4ZkMPDzPWiqA/115Hb60VpLCIZvp2atjbYjEbkJdmvM9i34+obYcIkF3rvXnjhhTjrrLPwwQcfYM2aNTFXTTLpD4udDKempgahUChr++sQXrcrqSksCZvFKEqk3d7o5jhJPW7kSmMRuTad+OAnwTNcZCe3UP4GeLkGHTRqFdqdnqSInfGFhUgVkm+ngauyUip4zjvvPFx00UX45JNP8Oabb/K09CyD01gZzoEDB5CXl4ewNQ/ZitfVhdz8yUl/HkpjES12FyYYTt2vxe0KQKtXQ6NVQS56K7LcwqRsLbIOuI9GRFAlVtWsKkV8WVFzwUT7dro8HnS43Thr2jSkikjfztQJM1L2vNnEcLO01CjF5PwleOfvq7Hzk2qcf/rFWHF16v7tmSyN7Hz88cfCQFZeXi4+zMg1Hwmdcd5///0oKyuDwWDABRdcIMy2kXR0dOC6665Dbm4ubDYbbrnlFrhcsVXEZCp0/EjsTJ8+XfY0hJydo8mgbMq1Jf25ck0GUSLdao8ulSWGcMoY1SHUfX6coUzKNPyTzn7lNidLFFqMaE9wRdax9nbQn0ay++sM5dvhfjvyQALzorMuQ0NzHf750Rvw+ZI/d43JcrHjdrtF86cnnnhiyPsffvhh/O53v8NTTz2Fzz77DCaTCStXrhzw5iShs2fPHpGLffvtt4WA+s53vpPCV6Fc6uvr4XQ6hdjJVsicTKLPaBkYtUgGOVTxlmtES5Rih9JYcvp1JHILcuFoO7n8XLpNOWLHBI8/ALcvujRhtCmsMqsNek1qex1VlFQK306XozOlz8v0Ulk+AZecezU6u9rxzDPPiPYcTGYjq9i5+OKL8fOf/xxXX331SffRF9Sjjz6K++67T/RLoAZRL7zwAhoaGvojQPv27cN7772HP/3pT1i2bBnOOOMMPPbYY3jppZfEftnO/v37hUAcNy77ppxLeJxd4tJoSX5khyjJM6M5il47olmfMyBrJZYEeXJI2Aw2bVLZucFigEYnT9PDwZCQJNoS5NuhqFWq/ToS5SVjRWl0bUNNyp+b6aWkqBxXXPgV4Wn84x//KE4OmcxFsQZlqiBqamoSqSsJq9UqRM3GjRvFdbqk1BX1UZCg/elDhCJBw+H3++FwOAZsmSp2pk2bJo5HtuJxdEGlUkFnTE2lTWm+RYidU5kfKTpBpedyp7EI8upQ40CPY6CIIAFEUR+lYNBqYNbr0OqIrZcRsbul6KRt3VENau0a+MdMG3B7KtBotCgrGYu6RhY7cmLLzce3v/1t5Ofn49lnn8XevXvlXhKTbQZlEjpESUnJgNvpunQfXRYXFw+4X61WizeutM9QPPTQQ/jZz36GTKatrU1sVH2QDZxW+/SQt/ccrYY61IXldX9MmdgJBEPodPpQYB2+N01bV6+wUEIay1bUG/Xqau2CydorCinK09XWhQlzJ0BJUHSnLQ6xMxQNnc3QqjUotCRv+OdIjCsbj807PkUwGBDih5EHin7fdNNNImPw97//XZSpr1ixImt9jpmKYsVOMrn33ntx991391+nyE6mpXooxafVajFhgrK+rFKNw+0ToxxSRWl+71iFpg7niGKHTMzk8THJNCoiEq1BK9JVJHbKJ5eL27xOL4L+IKyFyfc6xUJRrgnHWu0IhsPQqEZXxdZob0aprVh0v04F7s1bBlwv8Djhq6vD4dX/RGVB73GPxLR0SUrWle1IlVvjtAvQovbhucf/jvXv7BQDRdWDGpEuvXyiTKtkRoti8xulpaXisrm5ecDtdF26jy5bWloG3E/5V6rQkvYZCp1OJ6q3IrdMTGFNnjxZdA7NZrpc3pSKHaNeK56vqXPkhnFUnk4prByVMv4EKZVlbz1h0iThQ+QWKetvozDXJKJONN5hNARCQbQ62lFmGxgZTiVWowW5BjPqOxplWwNzAorkLJ63AucsX4mjtYfw1upXuPFjBqHYyA5FJEiwULfL+fPn90dgyItz6623iuvLly8XLnoa8rZo0SJx24cffij8EuTtyVboOJHZLpuPARHu7obL60+p2JGiO40jzJsiaKyEuW/quNyMf/kzuOuasaeuBVW+MeJDv7O6AeUtdkx/R1kTo20mvWguWN/ZgXH5+aNKYdH7Y2z+8CdFqaAivxR17Y1CwHHaRBlMmTAdedZ8fPDJ23ht1f+HC8+8DKVFJ0femPRCVrFD/XAOHz48wJS8Y8cO4bmprKzEnXfeKaq1pkyZIsTPT37yE9GT56qrrhL7z5gxA1/60pfwL//yL6I8nUYj3H777fj6178u9st0Hvng4JC31x7YiX3VHTDUAmuah94nG3B7/Oju7km92MmzYNuh4asB6YuNxkpYxslvTo4ctBkMheHw+mE16tHp8oomfkqDJqBTKqu2vQOnTYr/ceraG5BnssJikG+aOzE2vwz76g+jy+OEzaSsKFo2U5hfjKtWfh2rP30Xb6/5B5YvPAszp8yVe1nMKJA1hv75559jwYIFYiPIR0M/UyNB4oc//CHuuOMO0TdnyZIlQhxRqblef+LL669//avoI3P++eeLgW9Ufv7000ObVbOFpppDKCgdB61OGZEDOf06hMWY2uNQVmARESXahqvEorESFqtyxA4N2pTKukmMtbu8KFCg2CFKrGbUd3aKyEw8dPd0o76jSQgNuSmxFkKtUuM4p7IUh0FvFL14Zk6Zgw2ff4QPN7wnKnmZ9ETWyM4555wz4kA2Cus+8MADYhsOigK9+OKLSVph+uHzuNDZ0ojZK85HtmN3eaFW5cBs0KVc7BANbU5MHSJ6I01GV0oai9CqVbCa9KLSicQEVZRRtEeJFNvMaG51C99ORV7sY1DaHB3wBf0YVyC/2CGhU2YrEuJr9jgeW6A0qG3FikXniDTWx5+tESfSX/nKV0b0hDLKRBnuSCZhNB07LERi8TiuGrA7vbBZDKLqKZWQaDDptTjeZ/IdDKWwSIQZTfJXYkVC6SHqYSOVdhf0zfpSGjQBXadWo66jPa7fJ4+MXqNDYa48JeeDqcgvQ3NXG/zBxHWGZhLLxMqpuHrl10XBBzWxpawET05PL1jsZGAKq7CiChqtclIkcmF3eWCVIRVDYnNskRX1Q4xgIFo6XSi0mlIuwqKpdOpy+9Bod8Ji0EGvVWYlH/l2KKJDvp14oJTR2IIy5IxRxscflZ33oEf4iBjlYs3NEw0IqWCGRhO9+uqrnNZKI5Tx184kBK/LAXtrE0qrkj/hW+nQWVeniOzIk4qpKMwVYmeosz/qsFxsS01H51go6jtW1MdG8vAoFRrcSb6dUDgc0+85vS50ursU4deRMOoMKM4twLE2HgyqdKhp7WWXXYYvf/nLYsjyH/7wB9TV1cm9LCYKWOxkEE3HDokcc/FYTmF5/UH4AyHZKopI7PgCQbQPGsFAbRFI7JQpaAyDBEVzyLvT2OEUUR4lQ/OsqLFgXUdHzFEdaiJYnjewM7vcVBWOFeXw1P+HUT6zZ8/Gd7/7XRiNRvz5z3/G2rVrEY5ReDOphcVOBtFYfQhFY8dDrVWWF0Quvw5Bnh05qOhrxjc4lUXih0q8pU7LSoLSb5S6cvr8KFR4ZKfIYhFzsmiQZ6x+nRJrkRgToSQqCytEdRmnstIHKo65+eabRaHNJ598Iqan04geRpmw2MkQ3A47HB2tKK2aIvdSFFOJpcoZA0uKe+xIkGggX85gk7LUbFCq2FIaGlUO/MEQzHplC2YSZhTdiUXsBEIBNHe1KqIKazBmvVGksmpaOSWSTtCQ5bPPPhu33HKL8O9QWmvz5s1sXlYgLHYyhMaag1BrNCKyw1Bkx4NcsyFlc4+GYlyxFXUtg8ROhxN5FoNizb/0Ga3TqNHp7o2MKZkJhUVodTrh8vX2UzoVtW0N6O7pEVEUJTKhuBL1Hc2iLJ5JLyoqKkRai/rEvfvuu/i///s/dHZ2yr0sJgIWOxkAnUU0HNmHksrJUCksPC9r2bnMTfGqSmzCn0P+ocjIjlKjOvQ+cvsDyDMZ0NTpgtKpKiwETVg42toa1f7VrXWiiZ9Jp8wU3fiiseKyppWNyukIlaVTY9vrr78e7e3t+P3vf49NmzYJnx4jPyx2MoDOlgZ4nA5UTJ4h91IUgegA7HAjP1feL7Wqkjyxltrm3gGbYoBlhwtl+coUOw6PX6SwqopsaO5SvtgxarWosOXhcMvAYcFD4Q340NjZgglF46BUDFq9ME5Xt9TKvRRmFNAA5n/7t38TUZ5Vq1YJA3NrlIKcSR4sdjKA+sP7YLTkIq848+eBRYPbGxCVWPlWeSuKbGa9mMt1rE/s2F0+UaFVqlCxQwKHvDBTygrQ7uw1UiudySUlwrcTCIVG3K+2rV5cKjWFJTGxpFI0GKRZWUz6otPpRJTnW9/6Fnw+n5jduG7dOq7YkhHFTj1noiMUDKC59jDGz1rIU5P7aOuLSpBBWE7o34NSWZLYqW/r9e+UK7DsnGjpcov+OmMLrNh2tEGIH/pZ6WLno/37cay9DVNKSkdMYZXaikT0RMlUFVZAp9HicFMN+NQlvdn81lFxuXDsRdhh34y/Pv063v7/PsLpi89FeUlvyjKSpZdzy5BkwpGdNKe59gjCoRAqJk6XeymKob3LDYNOA6MCKorGl+ahscMhIk1kVqbUmskg/7oGQyk2EjfFVpPot0NbfcfQHaCVRL7JhAKzGYeah09lkYGZoiUTipWbwpJQ5aiEUflwcw17PTIEmn+2eN4KMW5Cp9HhnTX/wNoNq+Dx9o5lYVIDR3bSnPoj+5BfOhYGszKjBXKJHUphKSHSNb4kD93dPSK6U9fahXFFCoyU2I/B4Q3C52hDSZkKsAdRofOhtr4dPYXhk4+jrQpKYkpJCXbU1oo+NUNV3+1vbBQjJihqkg5MLZ2A/fWHUdtQjfFjJ8m9HCZBFOQV4fILv4KDR/di845Pcaz+KBbNOQ2zps4TJexMcmGxk8Z4nF3oaKrH3DMulHspihM7k8cWyfLcG4+2nxQxsQd68NrmIzh8vBXL50w4sY8y5lAKGrt6+xIVmXtnqpVbDdjf5ECnN4h8o/IiUVuatvT/7FN5cbSzBm/u+3DItNuqA4dQnrcYWrXyXsdQ5JttKLTkY++hnSx2FJqaihc6cZg2aRaqxk7C5zs34LPtn+BQ9T6sWHwOObYStk7mZFhOpjHHD++FWqNFSSV/IEpQmTcZlOU2Jw8YClpsw4FjzSLCU5KnTHNyQ5dXCB1qKkgUW/Ti5wb7wHEXSoRGgtCkeZrpNRiHxyemuKdDCiuSGRWTUd9YC7sjvmGnjLLR6/Q4Y8l5uOKir4rPiLc+eBkvv/wy7PaT38NMYmCxk6aEQiEcP7QbFZOmc2+dQVEdJZiTIyGx09zhFB4MuQaTjkQo3I0Wpw9l1hN9iSjKU2bVo65T+WKHGF+Uh7r2LjEvK5IjzR3QqFVisng6QT139HoD9hz8Qu6lMEmkuKAUV638Os4+7ULU1tbi8ccfx+rVq3maehJgsZOm7N69GwGfD5XT58q9FEXR0uGETqsWJd9KobzQCq8/gDE5Y5CTI7+PaDAtLj/C3T0idRXJuDwTOtwBOH3KH05JvYHC4W4cbz9hqqZuydUtnZhQnCdMoukEGZVnTJ6DQ0f3weePrkM0k55QZGfqxJm44447cPrpp+Ozzz7D7373O2zdupVN6gmExU4aQj4Qmr9SWF4JU26e3MtRFM0dDpTkWxRhTpaQTLPhcI9iU1hGrQpWw8AIYYWNxm2MQW0aRHeoeqzYZsbhphOeqYYOh0hrTipVkDkqBsi42t3TjX2Hdsq9FCYFaLVanHvuubj99tsxadIkvPXWW2LW1tGjo/MJMb2k1+kOI6ivr0dDQwOqZqyQeymKgjwxLZ0uzJ8qX9XNmzmHT7rNYfeiw+SHK9yGf4QOQKNVidvHQd5xFpJwru/0oNxmOEkgkmenwmZEbYcbs8oUWEU2iCmlBVi//5jw6eQa9TjU2I58s1Fs6YhBbxRm1t0HtmPO9AVQc7o6K7BarbjmmmuwbNky0YH5hRdeEOLn/PPPR3l5eqVjlQSLnTSEwpz5+fkIlSurBFhuOhxu0fW3JF9ZZfj2VjfMeQYxZbOjyYWSSuUIB4fdB5c/hMq8oT1OVflGfHK4FQ5fELl6ZX/Zjiu0QqtRC5/O5NICNHY6sWzKyc3b0om50xdh3+HdYiPBw2TXcFHqwLx//358+OGHePrppzFz5kycd955OLox+h5Y3KywF05jpRkulwt79+7FkiVLFJWqUYpfh9IuhTYzlIS9xY3CMgusBUa0NShrDEDTcQc06hxRfTUU5OOh+4+0Kn9WFqULyZ9zpKkD+463QqtWoaoovdO8FnMupoyfjh17tyAYDMi9HCbF0Gf8jBkzcOutt+LKK68UUf0nnngC6zZ9AJdbWZ8lSocjO2kGeXWoARUNmdv0CQ8MjKSpwymEjrqvfFoJhEPdcHT6MGlOsUiz1expRSgYhlrTm8pSgtiZYjMKkTgUdCwnFJhwtM2FuRW2YfdTCtPKC7HveAs+P3ocy6dWKuq9EC8LZy/D4Zr92H3wCyyYtUTu5TCy9fCxYtmES7Hv0C7sqNki3hMzp8zF/FlLRMqTGRkWO2kEDZSjFBZFdfR65VQbKQExUby9S7ZmgiOlsHq6e2ArMiFHNQbVu1vQetyBsgnyRxycXT64unwYlzfyB+WkIgsONjuFkflU+yrBqEyCrNnuxtSyAqQz7s29jRNJrk3o1uHz99/AuE7fkPO9TEtZBGUDVFVI6cxpk2Zi94Ed2LlvG/Yf2Y0ZU+Zi7vSFMBqU03JDaaT/aU8WQUKHpuauWMHG5KH663h8QYwtll9ERNLe6ILRooXBrIXOoEFeiRnNtb0DQeXmeLUdGq36pJLzwVAH5XyTFodalB82p5ERoXAPjFo1OlyZU7I9t3KGuNxxbI/cS2EUgFajExG/r1/xTcyZtgD7D+/G3956Dhu3roPbo/yUsxxwZCdNoCZTmzZtwqJFi2A2K8uTogSOt9hF87iSAuV0KKa0FRmSSyfY+m8jc/K+zfVwd/kAm0HWtdUfs6O8ygpV56lTU1NLcrHpaBvsngBsJ16ObHR65g15+9HmOgS7yzC+eBI+O6yCyTBX+B4Maf4nQ9Gc+VUzseXoF5hWNkmMlGAYvc4ghozOnr4Qew7uENGefYd3iSq+eTMXw2xUzueh3HBkJ428OoFAQDSdYoYWO+VF1iEHQcqFo82DYCAszMkS+aVmaHUqNB2TN7rT1uSC3xvE2CjTaePzTaIXz75m5U5CD3eHsb/+KMYWlGLhpFmwu51o6Bh+Gnq6Ma18EmzGXGw4uFX032GYyPETNFT061d8CwtmLcWRYwfxtzefw6dbPkRnZ6fcy1MEyvlmYIaFRM7GjRuxcOFC5OYqq6xaCQQDfrR0OsVYBiXR1uiE3qiBydo7XJOgDsolVTaRyiIhJBd11Z2wWPWw5kXn/SIfzLSSXNS0u+HxK7Mq6GhTHfwBP2aMnYzC3HyU2Aqw9/jhjOlCS0J++ZSFaHd1Yl/9yf2cGEan1WHB7KUivbVo7nJU1x4S3Zj//ve/4/jx48hmWOykSVSHzMkc1RmatvpjIi2jJLHTHe4WZeaF5Sd3cy6bYBOm5WNH5Bny6HEHRBVW1eT8mNoXTC6yQJ0zBrvrWqA0AqEgDjQcRVVxBSx9Js3ZlVPh9LpFaitTKLYWiiGh26p3w+5WbpSNkd/TM3/mYnzjyptx6aWXoqmpCX/605/w7LPP4sCBA6KgI9tgz47Ccbvd+OSTT0QFlk0JZgkF0lhzEBV5ZliMyqlQa29yIegPo6Tq5AaCZFQuHpeLmoPtmDC1AKoUl0cfO9QOjUaFivGxvZ+06hzRSXlHYztmVBSJyielsO/4EWFOnl4xqf82mylXTDvfW3cYS6oWDlnFlI4sHD8bDZ3NWLd/Ey5bcL6Yo8UwQ0Fdt7sb87G06hLUaqqxc/s2bHj/97Dm2jBn+kLRw2lwZ+5MbULIkR2F89FHH4mz77PPPlvupSiSgN8nIjtKKzlvPtaF3AIDjJahBUHF5Hz4fSHUHkltPj0YDIvnrJyUF1evn6nFFhi0GuyoboRS6HI7cbSpFjPGToJRN1DQzBw3WaQOPzu8PWPOZqn8+Kzpy+DwuLApg14XkzxycnIwfuwkXHHhV8SWZy3A+i1r8f+98Wd8tv1TOFzKqBBNJhzZURCPfHBwwHVHRys2vrsKUxcsxx/WZ3e+dTiajx1GD3owoVw5PVV8nqDorzN5fumw+5AIGjchD4f3tGDseFv/vKxkQ9EkSvlVTYnveFGTvvkTyrBh/zExaLNc5tEcZNTdVr0HJr0Rk0tPHp+i02gxf8JM7K05isrCOkwsrkQmQNVYp01ZgPUHPke+yYbFWCr3kpg0oaSoHBcWlaPLaRdDZqlPz6792zC2rEo0KVzcPV6Io0wj815RhkBna3s/+0hMNa+aPl/u5SiWhqP7UVA6Dka9Fkqh/nAH1FoViipGFgJTZhcjHO7GkX2tKVlXwB/C0f1twqtjMMY/52p8kQ1leRZsPnwcwbB8JmviUEMNOl1dWDRp1rAf0FSdRemsTYe2w+HNnB4kU0onYObYKdh8ZAdqjh+RezlMmmG12HDawrNw3VW34Myl58Pr82DVujeFoZmsE11dmRXt4ciOQjl+aDfsrU1YuvIa5Kg4Jz8U9rZmdLY0YsE5l1CMB0og4Auh6Zgd48iLox75XIIEx6QZRTi0pwVllVZYaVhoghi/+uQ009baDthanTjPmAv9EPdHC6VVl0wei3e2HcDWIw04beo4yEGHqwt7jx/BlLLxKLCMXEJPVUxvb/8Qa/duwCXzz4NGlRkffYsnzoXH78WH6/+JL519JcpL5fm3YNIXtVoj+vLQ1tLehJCtBR9//LEYPjphwgTMnz9fzOfSaJQ9CPhUcGRHgXicXTiwdT3GTpmJ/JIKuZejWGr2boPRYkXR2AlQUlSHPCLlE6PrXzNpeiHMuTrs3FwvKriSRYfbjwPNDswut0GfgLlcZE4mwXOkqR3VLanv4+EPBvDZwR2wGS3Cl3MqtGotzp25HC6vGx/v2yTMzJlAzpgcnDl9KcqKx+K9dW+gtqFa7iUxaUxxQSmuuuoq/OAHP8AVV1whOva/+uqr+PWvf4033ngDNTU1advKITNObzIIeiPtWv8BNDo9pi86U+7lKBa3w47mY0cwY+lZiskvd7m8aDjaibFT8qM2/+aocjB36VhsXH0E+75owqyF5QlfV3dPDz6raYfVoMX0kgR4bOzHxMVEbQ+azWF8tnMXzNNKUDTU5HTbyT6a0UJCZdPBHaKJ4LKp86KuRsozWXHOzOX4cM96rD/4Oc6YtliIhXSHXv+FZ12GNZ++i/c/fhtnLDkX0yfNlntZTBqj0+nEsGnaqCnhF198gR07dmD79u0wmUwi0jNz5kyMH58+/h4WOwrj8I5NsLc1YcmFV0OtVY4PRWkc3LYBOoMRFZN6ZwYpwWO1YVc1tAY1xsZo/rXlGzBzQRl2b21Ars2AcVFGhaLli+Od6PQEcNGMsoROLad01rLxBXD7Q/joUAsumF6KPGNy37Mk3NbvPwa724wzZiyCURdb6q8ivxRnTF+KT/ZtRnd3WERFMqF02791O5brSqAKN2DNqy+gvmIylkycO+Rr46GhTKzT142oxPKJ40Saq7ruMNa8/hne+OsHonNz1dhJGD9uEsqLx55Uxq6kknYWOwqiufYIju7eimmLVnD6agTaGmrFsZp75kVQjfDHlUoO1ragvsWOSctLTunVGYrKyflw2H3Y9Xk91JoclI07uT9PPNR2uLG30YEF4/JQaE58XxwST2dNKcKHB5qxen8Tzp1aEvfzDDfvSoIiOZsP7URTpxHLps49pU8nkgZXff/POkMOZo6fiM2HvkCDo1FEhyJ78JSbK9K2w/Jpkxciz2TDliM70NLVhuVTFqEoN1/upTEZwJgxY1BSWCa2ZfPPQFtHixA+R2sP4sCRPVCpVCgtrsC4sipR2WXLja1pabJhsaMQ6urqsPOTVSgdPxnjZy6UezmK7quz97O1yC+tQNn4qVACrXaXiOrMGF+C9tL4ep7Qh8KsReUIhbqxfWMdQsHuUUd4Gru82HC0DZX5JswoTV6JuE6twvnTSvHRoV7Bs3R8ASYWJnbypi/gx2eHdqDT5cBp0xagLG90fZXK84tx5swlwvezZudGLJw4E+X5JUh36H00vXwSiiz5YobWuzs+FFVbCyfMhl6jnCaQTPq/z4oKSsS2ZN4K2B0dON54DHWNx7Dliw3YtO0TmE0WVJRWCgFEUR+5YbGjABoaGvDiiy/CWliCOadfqCg1rDQ/EwlCmoW1+IKrFHGcOp0efPDZPuTnGrFs9gS8i5NDwMNRZ/eedFvBtAJ0+kP45OMalB/vwviZRf2RonExTEmv6/Rg/ZFWlOTqsWJiYdKPFXVXPm9qCbbUdmDj0TY0dHmxuDIfiehZ3GRvw7Yju0U/pTNnLo4pojMSBRYbzpuzHNuO7sHGAztQkV+M2VXTgDSfkE7QMbp04XlihMb2mt2oaT0uRBCNmugdpsEwiYE+W6hJIW3UlTkUCqKxpR51jTVoaD4uoj5EjWcebr75ZvnW2cPtN+FwOGC1WkVfgVQP2qyursZLL72E4uJitJcuhUbLZ18Sp9U+PcCU+vH2wzha34aVp81UxBys5g4HVm8+AINOg4tXzBKXb+aMfkAj/Uk2VttRvadFjJaYOKcYecUmVOYZT/m71Lfn8N5WtL11FJX5RiyfUCgaAaYKWjsNC91a1wEq2pg5ZSqmlRdCo1bFnMaikuo9dYdQ29qIYmsBFk+enZSRD7Tm4+1N2FV7UAwSnVU+A7PGTkFRrnIaVY4Gb8CH3XUHcLDxqPA8zTj9PDEmgM6408VcyqQvPr9XiJ/pp5dg3ryRU9XJ/P7OGLHzxBNP4H/+53/EwDM6oI899hiWLl2qWLFDh50mma9evVr0Mvja176GJ9bVpOS5003sOD0+IXRaOpw4e+EUTKwolHVdwVAYOw/X44tD9SjOM+OCJdOh1/V6hxIhdiQ8Tj8Of9GMrjYPcvMNWDivFCUVuUN6guj91FzvxIFdzXA7/TjLpxJzrOSKfvmCYexusONQqwuqMWMwvsAkxFexRY+cYdbUqb1CvA67p3f8Q21bI7RqNWaOm4LxRRVJfy2hcBjHWuvR2mEXzQcLLXmoKhyLysJyWI3ydopOVLn+gcYjqNGE4XI7YDAYMX7sZJSXjBWbPkazN8PEQrIMylkldv72t7/hxhtvxFNPPYVly5bh0Ucfxcsvvyymu1LERGli5/jx43jvvffEJU0yP//888UZ1uBxEdnOrP2/w77qJuyraYZOo8bZCyejrDAxxt1YkASMzx1Ac50DTTV2hAJhUWJeOa0QYxJY4TQY+vPsbHHj+KEOwBUQQqewxIzicgty8wzoCXejvdWN4zV2uB1+5BebRPn63M9T3/tmKDyBEA62OFHd7obHH4JOk4NyqxEFJi0KTDpRvUXRhhanD4c884Wp1uFxw6DTY3Jppeh8rFGnNtteairD8fZGHGk+hvrOZoTCITFJvchSgEJLvkgR5RrMwgOjhFRqrBiXLO6tqqk9hGP1R+Fw9nbKtVnzUZBXhAJbIfJsBbCYcsU2UoUNw0QLi50EQAKHpoI//vjj/d6OcePG4Y477sCPf/xjWcXOb98/gHAoCI/DLrr9Nh07JC4ttnzMWHZO3FVXI/2zDXtf3+1D3U9+iMG3S9d7hmkiNaC51ODH7BlmHT096MaJ36PHDocCCAX8CPi8cHfZ4ezqgKujBeZja6BR52B6ZQnmTC4b8UtPPFP/0438lh68LvqyJULhbhG1CQTDCITC8AdCcLi8WOupgcvug9cTEmKjqNyM8sn50NPIhTj/ek71Z0d393T3IBzq7t8M3T1orOlCU4MD9jYvvJ6A2IfWQQKootKGojKzqOaq2t4JjSpHpLA0OWNA38kUUZEue7fefPvg7+uBV+P/Mqfj2hZaLiImtHW4HTje2YqWrg7YvS74AgEEwyFxLDRqDczGXCEkyHxcklcEvUYLdY5KVBlRCTWdEKj7LsW66b/+y9610mXkffEQWY1FQqfB3oImewtaHR3ocNlFVZjUqNCiN0Gv1Yn0GokfutSoNKJDM1WnaHLUUKtUYngn9fTp/TeQLunfoPe1SLcPtWR6LQOuj1JgDS49d7mdaGg5jpa2RnR0tqHD3oZgKNh/v15vgNlogdFggk6rExEgnVYPnU4vXhe9TrVKA7WafqbXqxY/977Ovn+rAVvkaz5xfajXlY5ikhkaFjujJBAIwGg04pVXXhGdHyVuuukm2O120fVxMH6/X2wSdJAqKytFRVQixc7FF1+Mz7Z8fuKbWPqAU6nEH3gqSZd/ZunDj77QdGMCKfWbDAWlYHyqMHLUOUJEqNX0byfPWnTqHBFForlbdFxyVGNEt2Z6a5FQo0quUDCMnnAPdPYTX1aJJNbXHsLJkTgSL0KwRHyN023dY1RCIJGYSFyH4zEx761TaUf8Owp194o3WicJfhpGSlu47zLlf2tjEnsUaP20dff/3Pua+jfpdnFGM+iXUzyGg7VQ+lBYXCCaEyZD7FBwg77vSfRkbDVWW1ubaGldUjKwbJSu79+/f8jfeeihh/Czn/3spNvpgDHKwSP3ApgE4JB7AQzDKABXjXNEMTJanE5nZoudeLj33ntx991391+ns7OOjg4UFBQkNGwqKc5ER4wyET5WscHHK3r4WEUPH6vo4WOljGNFUUYSOuXlI4/aSXuxU1hYKHLGzc0Dp17T9dLS0mHnftAWic2WvFJm+sflP4bo4GMVG3y8ooePVfTwsYoePlbyH6toIkZp32RBq9Vi0aJFWLNmzYBIDV1fvny5rGtjGIZhGEZ+0j6yQ1BKigzJixcvFr11qPTc7XbjW9/6ltxLYxiGYRhGZjJC7FBDvtbWVtx///2iqeD8+fNFH5vBpuVUQ6myn/70pyelzJiT4WMVG3y8ooePVfTwsYoePlbpdazSvvScYRiGYRgmoz07DMMwDMMwI8Fih2EYhmGYjIbFDsMwDMMwGQ2LHYZhGIZhMhoWO6PkiSeewPjx46HX68VA0s2bN4+4P01jnz59uth/zpw5ePfdd5EtxHKsnnvuuZMGCNLvZQMff/wxLr/8ctERlF7366+/fsrf+eijj7Bw4UJR7TB58mRx/LKBWI8VHaeTB1OOEVWcmQ6NyaGByRaLBcXFxWKW4IEDB075e9n4mRXPscrWz6wnn3wSc+fO7W8YSP3t/vnPfyruPcViZxT87W9/Ez1+qKRu27ZtmDdvHlauXImWlpYh99+wYQO+8Y1v4JZbbsH27dvFHxBtu3fvRqYT67Ei6A+nsbGxfzt27BiyAeoRRceHxGE0VFdX49JLL8W5556LHTt24M4778S3v/1trFq1CplOrMdKgr64It9b9IWW6axbtw633XYbNm3ahA8++ADBYBAXXXSROIbDka2fWfEcq2z9zBo7dix++ctfYuvWrfj8889x3nnn4corr8SePXuU9Z6i0nMmPpYuXdpz22239V8Ph8M95eXlPQ899NCQ+3/1q1/tufTSSwfctmzZsp5//dd/7cl0Yj1Wzz77bI/Vau3JduhP9LXXXhtxnx/+8Ic9s2bNGnDb1772tZ6VK1f2ZBPRHKu1a9eK/To7O3uynZaWFnEs1q1bN+w+2fyZFeux4s+sE+Tl5fX86U9/6lHSe4ojO3ESCASEkr3gggv6b8vJyRHXN27cOOTv0O2R+xMU3Rhu/2w+VoTL5UJVVZUYIDfSmUK2k63vq9FAjUfLyspw4YUXYv369chGurq6xGV+fv6w+/B7K/pjRWT7Z1Y4HMZLL70kImDDjWuS6z3FYidO2traxD/s4C7NdH24/D/dHsv+2Xyspk2bhj//+c9444038Je//EXMO1uxYgWOHz+eolWnD8O9r2jSsNfrlW1dSoQEzlNPPYV//OMfYqMvpXPOOUekVrMJ+nuidOfpp5+O2bNnD7tftn5mxXOssvkza9euXTCbzcIz+N3vfhevvfYaZs6cqaj3VEaMi2AyDzoriDwzoA+NGTNm4A9/+AMefPBBWdfGpC/0hURb5PvqyJEjeOSRR/B///d/yBbIj0IeiU8//VTupWTMscrmz6xp06YJvyBFwF555RUxq5J8T8MJHjngyE6cFBYWQqVSobm5ecDtdL20tHTI36HbY9k/m4/VYDQaDRYsWIDDhw8naZXpy3DvKzJLGgwG2daVLtDw4Gx6X91+++14++23sXbtWmEuHYls/cyK51hl82eWVqsVVaCLFi0SlWxUNPC///u/inpPsdgZxT8u/cOuWbOm/zYKW9L14XKVdHvk/gQ5/YfbP5uP1WAoDUahUkpDMAPJ1vdVoqAz0mx4X5GHm768KcXw4YcfYsKECaf8nWx9b8VzrAaTzZ9Z3d3d8Pv9ynpPJdX+nOG89NJLPTqdrue5557r2bt3b893vvOdHpvN1tPU1CTuv+GGG3p+/OMf9++/fv36HrVa3fPrX/+6Z9++fT0//elPezQaTc+uXbt6Mp1Yj9XPfvaznlWrVvUcOXKkZ+vWrT1f//rXe/R6fc+ePXt6Mh2n09mzfft2sdGf6G9/+1vx87Fjx8T9dJzoeEkcPXq0x2g09txzzz3iffXEE0/0qFSqnvfee68n04n1WD3yyCM9r7/+es+hQ4fE3933vve9npycnJ7Vq1f3ZDq33nqrqBb66KOPehobG/s3j8fTvw9/ZsV/rLL1M+vHP/6xqFKrrq7u2blzp7g+ZsyYnvfff19R7ykWO6Pkscce66msrOzRarWivHrTpk3995199tk9N91004D9//73v/dMnTpV7E/lwu+8805PthDLsbrzzjv79y0pKem55JJLerZt29aTDUjl0YM36fjQJR2vwb8zf/58cbwmTpwoymCzgViP1a9+9aueSZMmiS+h/Pz8nnPOOafnww8/7MkGhjpOtEW+V/gzK/5jla2fWTfffHNPVVWVeN1FRUU9559/fr/QUdJ7agz9L7mxI4ZhGIZhGPlgzw7DMAzDMBkNix2GYRiGYTIaFjsMwzAMw2Q0LHYYhmEYhsloWOwwDMMwDJPRsNhhGIZhGCajYbHDMAzDMExGw2KHYRiGYZiMhsUOwzAMwzAZDYsdhmGyjnPOOQd33nlnQvYdfD81pf/Od76D/Px8jBkzRgwaZRhGXtQyPz/DMGkAfaHPnz8fjz76qNxLURyvvvoqNBpN//X33nsPzz33HD766CNMnDgRhYWFfPwYRmZY7DAMkxACgQC0Wm3GPdepoAhOJEeOHEFZWRlWrFgh25oYhhkIp7EYhhmRb37zm1i3bh3+93//V6RlaKupqRHRittvv12kcCh6sXLlSrH/+PHjT4pgUFTjv/7rv8TP3d3deOihhzBhwgQYDAbMmzcPr7zyyohrGO65onkst9uNG2+8EWazWYiQ3/zmNwPup/3nzJkjfr+goAAXXHCB+J1I6Hl++MMfCmFTWlra/1oGp7HoWN1xxx2ora0Vx4mOxXDHbyh+8Ytf9O8TuXFEiGFGB0d2GIYZEfqSPnjwIGbPno0HHnhA3FZUVCQun3/+edx6661Yv3591I9H4uQvf/kLnnrqKUyZMgUff/wxrr/+evGYZ5999rC/N9RzRfNY99xzjxAbb7zxBoqLi/Ef//Ef2LZtmxBgjY2N+MY3voGHH34YV199NZxOJz755BPhuxn83HfffTc+++wzbNy4UQiY008/HRdeeOFJx2rSpEl4+umnsWXLFqhUKhGBGu74DYaE0s0339x//f7778f777+PL3/5y1EfX4ZhTobFDsMwI2K1WsUXttFoFFGNSEhgkFCIFr/fL6IXq1evxvLly8Vt5Gv59NNP8Yc//GFEsTP4uaJ5LJfLhWeeeUYIovPPP79fuIwdO1b8TGInFArhmmuuQVVVlbiNojyDmTt3Ln7605/2r+Pxxx/HmjVrThI7dKwsFosQOZHHarjjNxj6XdqIn/zkJ0LokPdHWi/DMPHBYodhmLhZtGhRTPsfPnwYHo/nJJFAHpwFCxbE9FzRPBb5Z+j6smXL+u+nVNS0adPEz5T2IhFEAodSYxdddJGIouTl5Z0kdiKhdFhLSwuSBUV0/u///k8IHUqFMQwzOljsMAwTNyaT6aTbcnJyTkoDBYNBcUmRFuKdd95BRUXFgH10Ol1MzzWax5KgCMwHH3yADRs2iCjKY489hv/8z/8U6SryAUlEVlsR5KMhH08yoAjSCy+8wEKHYRIIG5QZhjkllIYJh8NR7Ut+FEoPSTgcDlRXV4ufZ86cKYQIGXgnT548YBs3blxMa4rmscg/Q0KFxItEZ2en8NBEChfy3/zsZz/D9u3bxWt97bXXIMfxI6FDaTYWOgyTWDiywzDMKaEvXhIMVEVEVU2Dy60jOe+880Sfmcsvvxw2m02kZCiCQpAf5Qc/+AHuuusuERk544wz0NXVJUzHubm5uOmmm6JeUzSPRWu95ZZbhEmZKq3IoEyRG4o+EfSayHtD6Su6j663trZixowZSPbxk9Yg8fOf/xxPPvkk3nzzTej1ejQ1NYnbKaUWbaSKYZihYbHDMMwpIVFB4oGiKV6vtz9SMxT33nuvuP+yyy4Tht0HH3xwwP50naI/VEl19OhRIYgWLlwoqqRiJZrH+p//+R+R8iLxRQLp+9//vhBFBIkiquCi0m6KQJFJmUrTL774YiT7+EVGbijtR+ukNUhma4nNmzdjyZIlCV0Pw2QbY3oGJ9cZhmEYhmEyCPbsMAzDMAyT0bDYYRiGYRgmo2GxwzAMwzBMRsNih2EYhmGYjIbFDsMwDMMwGQ2LHYZhGIZhMhoWOwzDMAzDZDQsdhiGYRiGyWhY7DAMwzAMk9Gw2GEYhmEYJqNhscMwDMMwDDKZ/x/rkGv5iaqmtAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "tomo_nz = glass.tomo_nz_gausserr(z, dndz, phz_sigma_0, zbins)\n", "tomo_nz *= ARCMIN2_SPHERE * (z[-1] - z[0]) / 40\n", "\n", "for (z1, z2), nz in zip(zbins, tomo_nz, strict=False):\n", " plt.hist(\n", " ztrue[(z1 <= zphot) & (zphot < z2)],\n", " bins=40,\n", " range=(z[0], z[-1]),\n", " histtype=\"stepfilled\",\n", " alpha=0.5,\n", " )\n", " plt.plot(z, nz, \"-k\", lw=1, alpha=0.5)\n", "plt.xlabel(\"true redshift $z$\")\n", "plt.ylabel(\"number of galaxies\")\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "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.13.2" } }, "nbformat": 4, "nbformat_minor": 0 }