{ "cells": [ { "cell_type": "markdown", "id": "e9036402-5b70-4886-aabe-22f0c1da094c", "metadata": {}, "source": [ "# Visualizing the results" ] }, { "cell_type": "markdown", "id": "355cfcd0-b9a2-4c21-9213-b74ef304e9f1", "metadata": {}, "source": [ "There are multiple ways to visualize the results we [retrieved](tutorial_1.ipynb) from the data. One way is just to plot multiple sample paths as you already saw just after the MCMC sampling. Another way is to calcurate and present the Highest Density Intervals (HDI) to quantify the posterior probability distribution, which is demonstrated in this notebook." ] }, { "cell_type": "markdown", "id": "5f905d60-47e9-41ee-951f-6ccc07d8e3a4", "metadata": {}, "source": [ "Let's import some helper functions in FRAP." ] }, { "cell_type": "code", "execution_count": 1, "id": "52a41d8b-ff38-47f1-a39f-b75cb210b196", "metadata": {}, "outputs": [], "source": [ "from frap import visualizer" ] }, { "cell_type": "markdown", "id": "cfd5b8c4-7064-4d23-9bda-5c3c0282afee", "metadata": {}, "source": [ "We first perform the Kernel Density Estimation at each radius using the MCMC sampling results. Here, `nskip` is a number to thin-out the sample just to reduce the calcuration cost." ] }, { "cell_type": "code", "execution_count": 2, "id": "6603ec3a-b6f2-4b31-82c6-936cb5bd7b56", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Starting KDE calculation...\n", "Processing key: log10_Sigma_d\n", "Processing key: log10_a_max\n" ] } ], "source": [ "kde_results = visualizer.calc_kde( 'mcmc_reuslts.pkl', nskip=100 )" ] }, { "cell_type": "markdown", "id": "d7d41f1c-f373-4f49-ba67-34161324f266", "metadata": {}, "source": [ "Once you have the results, you can calcurate the HDI with a specific source distance (for conversion of the radial unit)." ] }, { "cell_type": "code", "execution_count": 3, "id": "7fb787cb-14fc-4eb3-a1d9-1ef9a5324dbd", "metadata": {}, "outputs": [], "source": [ "hdi_data = visualizer.calc_hdi(kde_results, D=114.9)" ] }, { "cell_type": "markdown", "id": "92db9341-3700-4aab-9f65-d85c05fd20be", "metadata": {}, "source": [ "Then we can make a plot like this." ] }, { "cell_type": "code", "execution_count": 4, "id": "de144a67-684e-470e-b312-3a86e9146303", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAt4AAAK4CAYAAABDOfoYAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjksIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvJkbTWQAAAAlwSFlzAAAPYQAAD2EBqD+naQAArPFJREFUeJzt3QmYFNXVN/DDvoMQVFZBcAnIaBDRAEIgKAE+iRHyAgIBTdQYjXF78dEYAf38NCrE6IsaTRQwEkgQ14RVxA0RRFFBSVS2EHBUhAADAgLzPef63smdmqruWm5V3Vv1/z1P69DLTHd1d9W/Tp26t1p5eXk5AQAAAABArKrH++sBAAAAAADBGwAAAAAgIah4AwAAAAAkAMEbAAAAACABCN4AAAAAAAlA8AYAAAAASACCNwAAAABAAhC8AQAAAAASUDOJPwLhHDlyhLZt20aNGjWiatWqYTECAAAAGIbnotyzZw+1atWKqlcvXNNG8DYYh+62bdum/TQAAAAAoIgtW7ZQmzZtCt4HwdtgXOmWb2Tjxo3TfjoAAAAA4LB7925RKJW5rRAEb4PJ9hIO3QjeAAAAAOby0xaMkysBAAAAABKA4A0AAAAAkAAEbwAAAACABCB4AwAAAAAkAMEbAAAAACABCN4AAAAAAAlA8AYAAAAASACCNwAAAABAAhC8AQAAAAASgOANAAAAAJAABG8AAAAAgAQgeAMAAAAAJADBGwAAAAAgAQjeAAAAAAAJQPAGAAAAAEgAgjcAAAAAQAIQvBPy9ttv03XXXUennXYaNW7cmI4++mjq06cPPfPMM0k9BQAAAABIEYJ3Qu6++26aMWMGnXnmmXTPPffQzTffTAcOHKALLriAJkyYkNTTAAAAAICUVCsvLy9P64/nybJly6hbt25Ut27diusOHz5MZ599Nq1atYo+/fRTatasWaXH7N69m5o0aUK7du0SVXIAAAAAMEuQvIaKd0J69epVKXSzGjVq0NChQ+nQoUP04YcfJvVUAAAAACAFCN4p27Ztm/g/93wDAAAAQHblKnjv27eP5s+fT7fffruoNLdr146qVasmLpMmTfL1O/bs2SPuW1JSQg0bNhSHFrp3705TpkyhgwcPBno+W7dupWnTptFZZ51FHTt2DPmqAAAAAMAGNSlHVq5cSYMHDw79+M2bN1Pfvn1p06ZN4t/169cXJ0hyjzZfZs6cSUuWLKGmTZv62gngEyv58Y888kjo5wQAAAAAdshVxZtxKO7fvz+NHz+eZs2aRS1atPD1OO7DHjJkiAjdLVu2pMWLF9PevXtFgJ49ezY1atSIVq9eTWPGjCn6u7gyzhV3HmKQw/qpp56q4ZUBAAAAgMlyVfHu3bs37dixo9J1N954o6/H8lCAa9asET/PnTuXevToIX6uXr06jRgxgo4cOUKjRo2iefPmiao3h3s3X331FQ0fPpwWLVok2kw4gAMAAABA9uWq4s2jiITFwZv169evInSrRo4cSccff7z4+fHHH3f9HTx8IIfzZ599lh588EEaN25c6OcDAAAAAHbJVfAOi9tJeBxuNmjQINf78AmaAwcOFD9zNduJK+IctJ988km699576fLLL4/5WQMAAACASXLVahLWunXrRHBmXbp08byfvK20tFS0tKgT4nBPOfdzc7W8efPm9MQTT1R6bM+ePalDhw6xvQYAAAAASBeCd4Cxtlnr1q0976fexo9Rg/dbb70l/r98+XJxceJ+b6/gzTMiqerUqSMuAAAAAGAPBG+fY3dLPISgF/U29THspZdeCvcOEVHbtm0r/XvixIm+xx0HAAAAADMgeFtgy5Yt1Lhx44p/o9oNAAAAYB8Ebx94jG71REsv6m3qY6Li0K0GbwAAAACwD0Y18aFVq1aVpnn3ot6mPgYAAAAAAMHbh06dOomJctjatWs97ydv49kw1RMrAQAAAAAQvH3gkyZ79eolfl6wYIHrfcrLy2nhwoXi5wEDBuCTBQAAAACVIHj7JGeZXLp0Ka1YsaLK7XPmzKENGzaIn8eOHev31wIAAABATuQueO/cuZO2b99ecZET4/CJker1ZWVlVYJ3SUmJqGwPGzaMlixZIq7nx3PovvTSSytmtuzfv38KrwwAAAAATFatnJNkjrRv3542b95c9H4ctKdPn17puk2bNlG/fv3E/2ULCgfv/fv3i3937dpVBPKmTZtqea48cU6TJk1o165dGNUEAAAAwEBB8lruKt5RQ/t7771HEyZMENPDV6tWjWrVqkXdunWjyZMn0xtvvKEtdAMAAABAtuSu4m0TVLwBAAAAzIaKNwAAAACAYdBqAgAAAACQAARvAAAAAIAEIHgDAAAAACQAwRsAAAAAIAEI3gAAAAAACUDwBgAAAABIAII3AAAAAEACELwBAAAAABKA4A0AAAAAkAAEbwAAAACABCB4AwAAAAAkAMEbAAAAACABCN4AAAAAAAlA8AYAAAAASACCNwAAAABAAhC8AQAAAAASgOANAAAAAJAABG8AAAAAgAQgeAMAAAAAJADBGwAAAAAgAQjeAAAAAAAJqJnEHwHIi9fX76z4uWfHpqk+FwAAADALKt4AMYRut38DAABAviF4A2jgFbIRvgEAAEBC8AYAAAAASAB6vAEiklXt+R9tr3T9oBObY9kCAABABVS8ATSG7oZ1aoqLeh0AAAAAQ8UbQJNnl35U6d/n9zuxInxjhBMAAABAxRsgIg7XMnS3aFSn4uIM4gAAAJBvqHgDhOQcsYTDttPTSz4U/0fFGwAAAFDxBohY7eZwzaG7dM+BShe+rnr1auJ+GFYQAAAAELwBIpLhmo0e2ElcGIdvteoNAAAA+YbgDRACV7BltVsN2ZIM37L9BCOcAAAAAII3gKZqN49iUnbgUMVFBnK0mwAAAIDIDVgMAMHIfm21hYRDt5w0R06cI69z3hcAAADyCcEbIARuHeFKtjqSCQduOXqJDN9HjpRj+QIAAICA4A0QAbeSyMq2DN1q+L6g/0lff9GqV6PLfvc6lnXOjoxgNBsAAFAheAOEPKnSWe1WOcftlvdFEMsP3tHiyy0Lkp9ICZ8zAAAzYQIdgBC4gs3VbtlK4jZBjnodZrHMFzVs807a/x14YqJ/U8LETQAAZkHwBghY7VZxK4nfUCXbTdbe83+wzDP+OZEn0/Lng3/uMv5vsbzvzmEt5Qg6vEPI1wX5fAIAQPzQapKQw4cP01133UUnnHAC1alThzp27Ei33347HTr09bBzYA9ZvfZz4qSsOKr3RRtAtsle/kcu71mlBSnO0M34fAPnaDpptLoAAIA7BO+EXHXVVXTjjTfS2WefTQ888AD169ePbrnlFvrpT3+a1FMAzbia6CdYqSdZQn7wThdfOIDHscPlDN0y6PNF/k2G8A0AYA4E7wSsWbOGfve739EVV1xB06dPp0suuYT+8Ic/0HXXXUePPfYYrVq1KomnAZqri2F6aLkVALNYZpesLqvBV35GdI5q42xn4TYWGfTlha+TO3x8XxxpAQBIH4J3AmbPnk3l5eV0zTXXVLpe/ptvh+QDknrxG0rk2N3cOuK3jcDZboLJdLJLvrdxntTIn1UZ4osddVH7u7HDBwCQPgTvBHBF+xvf+AadeGLlk5zatm1LrVq1QsU7QRxa+EQ3Dkjcq80X/lkO++YWwJ0nVfJoJkFbR5ztJqg+ZpfbZ0OeWBm131oN3c5Jm7zI6juq3gAA6ctN8N63bx/Nnz9fnNA4dOhQateuHVWrVk1cJk2a5Ot37NmzR9y3pKSEGjZsSE2aNKHu3bvTlClT6ODBg56P27ZtG7Vu3dr1Nr5+69atoV8XBKMGaLcT0ZzUgMwhXR27G0O1gdtnpdAoIlGOdjgr3Ryo/XwG+T5yZwCTOAEApCs3wwmuXLmSBg8eHPrxmzdvpr59+9KmTZvEv+vXr08HDhwQ1Wq+zJw5k5YsWUJNmzZ1Df1c8XZTt25d+vLLL0M/LwjfFysDElcNOZCo/ddqoOHr5EgmcuzuoKNV8O/j38OPxbCC2eS3lYM/h2F22tTfH3SIQP68ys8+V90xxCAAQDpyU/FmHIr79+9P48ePp1mzZlGLFi18PY6H/BsyZIgI3S1btqTFixfT3r17RaDm/uxGjRrR6tWracyYMa6PlyHdzf79+6levXqRXhf449UXK0eekGMfy7YTvqgnVMpqd9gRSpxhHe0m2dupK/TZkC0fYXqt+bMoP4fqiZt+qaOroOUEACA9uQnevXv3ph07dtALL7xAd999N40cOVKMp+3HjBkzxMgkbO7cuXTOOeeIn6tXr04jRoyghx9+WPx73rx5ourtxH3cXu0kfL1XGwroo4ZcrvY5K47q4XiuSMsKN/9fnlBZbKZKPzCsYDbJMF3oSIj8zAQNvmro5s9P2M+e+hnHiZYAAOnITfCuUaNG6Mdy8GY89naPHj2q3M4h/vjjjxc/P/7441Vu79atG33xxRf00UeVT6zasmWL6P/m2yG5SU28cCDnYCLDtQzfal+337G7i41ugmEFs9nCVCwUq8HXT/h2tkdFnZBHPh5VbwCAdOQmeIfF7STLli0TPw8aNMj1PnyC5sCBA8XPixYtqnI7V8X5Pr/97W8rXS//zcEdklEsGMnRRzgcc+CWlW4d1W75+yWEH/s5T3gMEnzVkXScIZz/zbcFHcHEb9VbnmcAAADJys3JlWGtW7eOjhw5In7u0qWL5/3kbaWlpaKlpVmzZhW3nXbaaXTZZZfRgw8+KHrDue1l+fLl9Oijj9K4cePEyCgQHzmEm59gpAYbDkccUJgc/UTHFOD8PGQ1HezlHGbSzwmLMvjKKrY6ygn/Lvn5Us8tkJVuXaPo4ERLAID0IHgXwa0gUqFebPU2fowavNnUqVPFEIY8YyWf2Ml937feeivddNNNRd+k3bt3V/o396b77U+H/4SboBPeOOmqOGJ0E7vJ6rRaMQ5ywiMHdDmSDpM7YfLkXiedoVtSwz8AACQHwdvH2N3q6CRe1NvUx1Qs6Jo1Rcj2E7SdeKId1cSJE32PPZ53MiQFPSnN7b5xVBzlc8SY4HYdPZFHQzgshznhUY4yolbM5WdC/b1xhG71JE/5N00dXlBtwcF3BAByF7xvu+02bX94woQJ2n5X1vFJmI0bN674N6rdekebSIPabsLPEaHC7NAnP0cycKvhOGxo9XO+QZyfC1n15kvcf6sYvyeaMnxXACA3wZurrHySYBTl5eXid9gSvHmMbvVESy/qbepjdODQrQZv8M/vaBNu4trAO9tNwNy+beZsyZCBW0elWH7G+O/K8JtEwFSr3sxt4qi4grT6GtW/Lam97W7UXniEcADIdPDu06ePZ/B+5513aNeuXRX9zm3atKkYp/pf//qX+Pmoo44SJxrahHuxJX4tp556quv91HG61cdAekyeoEYGB656m364P48jlKg7RPJntf1Dd+hTf1dSYVI90VNWvZOqXDvDdsM6NWnmgnWVlrfzBOTRAzuJ/5cdOFQRvtGmBQCZDt4vvfSS6/Xct/zyyy/ThRdeKKriJ55YOUB8/PHH4kRCnladx8G+4447yBadOnUSE+XwyCZr1671HFKQb2M8G6bzxEpIh9y4h5npL0moeqdPTlLj9l7IYSTj6rlOkzzfQA4vyN+VoK9Phu2gk/KoYVtSx8xX8XCe8r4cwNXwDQCQq5MreSZHngnyiiuuECN3uDnhhBPoj3/8IzVp0oTuuusuOuOMM2jo0KFkAz5pslevXvTqq6/SggULxHTzbu0zCxcuFD8PGDAghWcJuttM4uZsNzGhzzaPZPXVebKkU1bfG+fwhkFaTpyVa68dl2K8wrbbfdQALneI0HoCALkK3hy2uf3EzygbfJ+HHnpIPMaW4M14rG0O3kuXLqUVK1bQWWedVen2OXPm0IYNG8TPY8eOTelZgi1tJm7tJupMllkMeCaTgVEdq132IefhvVCr3n53VtXQLVtC+PF+QnQU6u/nEM5Vc7SeAECuZq587733RCW7efPih/z4Ptzn/e6771Iadu7cSdu3b6+4yIlx+MRI9fqysrIqwbukpERUtocNG0ZLliwR1/PjOXRfeuml4t/chtK/f3/tz3vFBvNDpKlMbzNRK3cceoIergc9Pd3yPVAr23kI3ZKcrVVWrAvtuLqFbjnLaxRyhlh58YOr3xy+mXxObjOBAgCYolo5p8kIGjZsSAcOHBChln8uhANt06ZNxXB4znCbhPbt29PmzZuL3o+D9vTp0ytdt2nTJurXr5/4v2xB4eC9f/9+8e+uXbuKQM6vTxeeOId3ahat3kTnfqudtt+bB13G/038f+09/4dsCH/OE/hwomVyfd0s78tchmkO0s6edsmrrcQtcPsNzjqq3pI8YsGcvd952okCgOTJvMaDjBQbhS5yq8nJJ58sRjTh9pEbb7yx4H35PocPHxaPsQ2Hdq7uT548mZ566inauHEj1apVi0455RRxUulVV11FtWvXjuVvL16/gxo0aoyNR0bJHls5PjTLS5uDCSdTcmDDSXpfU6vebjNbFmorcYZtOQqJTlzh5r8j/7ba+83PzevES3yfAMAUkYP3RRddRFdffTX96le/EpXv66+/vkrlm1s5OLDyBDzcD37xxRdTGmS1Oiwen5tHZ+ELmMu2w8xqr7cz8CB8x0eeTJnVEyeDUF+/uhMoP5McaOVJjV4VZ3m/OMnnoYZv+ZzkiZdq+Fa/X5iABwAy0WrC7RaDBw+mRYsWiVBdr149MWoJj+Utx7detWoVffnll6JH+txzz6X58+eLIfrA36GLb171F/rh4G8hIASsZtrSPuDWMxvnmNF5p/Z22/IZSYI6LKAM3zLQFgrcatjmwMvkSY+6yX5ut50At+flfA5qFRzfKQCwstWEA/Rzzz0n2ky4lYSr26+88krFRDsy19eoUYOuvPJKMZwgQnfQZYzZDYOQFWOb2gf4uao9tOoIEza9Dlt2ctRqN1SdUZJ3SHgnUA3dxarbasiVPzt7xKMsb368/L2FKt9MrX6rz0mtgqP9BACsrHirPvnkE3ryySdFhfuzzz4T1x1zzDGiAs6jgWBGx3B7UJ2vnkNDB56GinfGTqz0qjaqM/aplW9U6KLDCZX+PosyoLpNcuNW4ZaS2JHxem7Fqt+Fdgrw3QKApCreWoN3FNySwideHnfccWk/FSODd7Va9ULNKpfXNgIbWwjcWk4Ywrf+EWR4meL75O+z6GzzUMN2lKCtVtjd/l0MPz+19aTQrJcM4RsATAjexjRac1W8Q4cOaT8No2GM5+wvIxlk5NBoHCbQaqQ3SModGezEFv8syosM3MVCNy9T9eK8zfmz877Ox7v9HpXaeuI1hKEM5Gp13Dn2t40nZQOAnYwJ3syQ4rvRQ3xBdsmAIQMNv+8yTOR9cp0ooUg+Vh17Gr3dhTnDrlxezjCu3t8tIHsFax3UHQLJT/iWAdwtfAMA5Cp4Q3GoyhQmT0i0rc3ELXxzVdZ5smWe3n85A6F8zVFmJES1W1/4dt4niSMHXn9HDd+ylaRQ+HZWv91mvQQAiBOCt0XcJrSA/8jaRpNDhdpykqfqnPpe8mtWL0ECOKrd0RQK1mm06rg9H3WHoFj4Zl7hGwAgCVjjWIaDB3pTvZdNFsj3V4ZG2XIiq95ZH+HEOZ40ky0icjpzJ7flof4e+Vj0dodjyudNDgHodiJmsaEGVc4Jd9ThDjHMIADECRVvCwz5zgni/+jz9kdt0bAZhwvZcqIGzqzsYPidxKVVk7oVAUrOpujseXe2o7j9Hobebru5naDJnCeB+q1884XDN1pOACApCN6QCRy0sjrhjNrvLQNk1tpq3MiZE916dGX1X7afuC0Teb1cZlnZIYP/8NN24ocavqU8fMcAIHloNbFMHloN4D/UQ+pydBv5Gcga5yRCbrMlSvI2vp9X+4la6ZazVOJ7kz1ebScySLtNea9SP2f8GHWGTbSdAIBuqHhbxGvDAZXbL7IWrmTLiSqr7SZqX7dauZQXSfbvynHOOYBzaJIX/rdaMUdvd7Y5207UlhO/lW91nO+sfr8AIH0I3hY4t2Oz3PT4RsFBK8vtBFluN1EnuJGBWVYh5cQtfOFRXuRIL85KphyjWQYo+Tu8KuKQXUHG+HZ+htz6vbP0XQOAdCF4W0SObAHe/d1ZJaveMkRm9fU6W0xkyHZO3MLXy9kK1eq3epGhm3dYbB3XHfT1e6ufl2JHFN2GGET4BgAdELwtkeVKri6yjzfLsvw5cGsxYc731M/oFfLfvLyy/pmAwtSqt5+2E3V8b+fRRT/hWx1ZB9VyADD25Mo2bdpQ3bp1034a1sBJP5Tr1hs+8iFft+097bLNxK3FRIZmr9fIj+PqN1fKneFbVsttXz4Q/kRL/vzwZ0Q92ZJ5jfHtdkKverJl2Mq3+hh8HgHyzZiK95tvvkkbN25M+2kYjzcMWW0zCENu0PLQgiM32FlrN3HbceLQXCh0y+ucrSfyoraoQL45+739TC2vVr1l5VvXDj6q4AD5pj14l5eX044dO2jLli30z3/+0/MCeqqdUDmAZr2alLUgKXecvIYPLPR+OsO3ehKmvD7rnweoyu09d35vgrScMOcJlzqgZxwgn7S1mvz1r3+l+++/n5YvX0779u0reN9q1arRoUP/6buDws7q0JReKf0UozPkuP/Z+TrlcHkcBGwNl+pIJl7V7rzukEA8orSc8OP48VHCt/NzKsO3rd9hAEip4n3DDTfQ+eefTy+88ALt3btXVL0LXY4cOaLjz+aKXGGrG4O8V0zU3uC8DBnntoG2/XPAOxFu/IQRvo9XKwrCTH4FqXoXGmZQrXozdaz4oBfmFdpt/w4DQILBe8GCBTR58mSqWbOm+P/7778vrj/66KPp448/ptdee40mTpxIzZo1o+bNm9Pzzz+PXu6Q8lLRDStPVU/b+7ydQUNWGGVvdlBq0EbgBh2jnKhje+ughm+3AI7eb4B8iBy8H374YdE6csstt9B1111HnTp9vSKrUaMGdejQgXr27CmC9zvvvENNmjShn/zkJ1SnDmZgjIoDF/q888n2WSydU8M7RenNRuiGQp8Fr53zYmN76wzffgM4gjhANkUO3itXrhT/v/TSSytdzy0lzuECp06dSp999hndddddUf9srslpsvNOnV48b0cDkpjFMs6Nvxo4MMMkxMWr5SjI2N5uPeBR+Wk/kRDEAbIlcvD+4osvqH79+nTsscdWXMfVbrcTLM8991wxVvff/va3qH8Wct4b6DaMYF6qnUm8TufnyutzFjScq/dzzg6IyW4gCUGr3hJXvXVVvqOMloIJegByPqpJ48aNq4xQwi0lO3fuFCdaNmjQoOL66tWri17wrVu3Rv2zAFXkcVIhOa67rv72QgG6WLj2G75lyPAKMXl7DyF+8jOlfkbVEU646l0oVMuqN4dzvl+Q4Qidv9v5WOdoKWG+yxgdBSBHFe/WrVvT7t27af/+/RXXnXTS14fBly1bVum+H330EZWVlYnwDRCV88TCvAY2Oaxg1KMfSR49kYHHbdg2AFOr3mrPtzOoy+uct8mf+bFych6v0VKiyusRUACbRA7ep556qujnXr16daWWEr7ul7/8JZWWlorrPv/8c9EHzidinnHGGVH/bO6ooVJuHOSManmjjv/MGzIe4SNPI5owfr38unWF1iQ22PJkMv7cOoNHmLG7AYJy7pwH6fWWZIBmzqAtb3O7qI+Xj3UTdZZM9IQDZDx4Dxw4UITsZ555puK6K6+8ko466igRxo877jhRFW/ZsiW9+uqr4vbx48dH/bO5xMHEbbi1vFY5uNKb52qpCSeU+g0J8j5qVc/tpLW8HrWA9KeSl/yuUwqFa7+PLVT51jlFPQCYI/KxrR/84Ac0bdo0atr0PxvMY445RpxAeeGFF4rp4T/55BNxPfd781jfHNYhGrnSDjvucZaYEEDTpg7LFzS8em2Yg2z4nfdVq9fO0K2GDQ45eZn8CMzA3w+vz3yxXu84uPWMq33fOo4E5fH8F4DMBu969erRuHHjqlzfo0cPWr9+vZhCfsuWLeKEy7PPPlucjAmge5rxvOGNqNpuo7Pyry5fr75T59Tbztuc75EzdKvVQYxmAmnhUOu2PvGaSl43+d0tFr6dzzkMnIAJYIZYz3LkYQU5bEO85Io5TxUNr2nGIRhn5c+tJcRNoduLPdYr0OTp8wtmVb2DjHCiW6EdZ+d3yRnEw46Agu8agMU93pAuHWfC20qeWJln6uuXJ9uGGVM7ic+UM8zkuT8f7B/hRDc/YV/2gDtnvwx6FBCzYgKkB8EbrG0zkRvGKFOM24xfN7dpuAWEMCdUxdnC49ZiwjCaCdg8wokuhYYpjDOE48RLAEuDN49q8thjj9GAAQOoRYsWVKdOHdFm4nXBON7R8EZBhi11JZ2XlahakcWJlV9TR0jws+EN22ISNBzI+xYb9SGPO05gDhOq3uroKEFnyXQL4X7lZbsBYIrIx5R5QpzBgweLyXI4gENy8jT5iLpxkEEuL689CHlCVtjHenELAYVm41NvdwvceO/A1BMtnb3eSZ1o6feky2J0j4gCAIYF70mTJtFrr70mKtmjRo2i733ve3Tssceiqp3gWfiytzfrK1mMZqLnJCqvarcbNVx7hQ8ZENwUCiwYChNMHFrQ7UTLtMO35DeEBw3fOOESwKLgPWfOHDEb5X333UdXXHGFnmcFgcefzTL1Nef5ZNJCQwo6yQ2uc4Ma5ITKQhVrVZRAkvWdRbBrvepW2EhjbG+371bQCjgq3wAZ7fH+7LPPRHX7kksu0fOMIJI8BfO8j2jithx4w+wM0eoU0k5ewd1v6A4DbSZgOrd2rbQ/t2FPwGQ4WgiQoeDNU8HXr1+fateurecZARTg7O/O64gmbiObOEOy342tM6jHGboltWqX5/cPzON2FEZ+Xk0I32oA1xm+81S0AbA6eHNP9+7du+nvf/+7nmcEgaV1GDQJcmOgnvQkYUQTb35actyWaRKhW1YU0WYCpnDb+XNWvTl8qyNK2VT9Tit8FzraBpAG9TOp87O5YsPO5IL3TTfdRM2aNaNf/OIX9NVXX0X9dZn09ttv03XXXUennXYaNW7cmI4++mjq06cPPfPMM5F/txqQ5Eo1Kyu5rLyONBXqAfcSZ+jm0JL0GMkAQRXaKZTh27Tqd7EAnnTbCdbfEPfn6/WAwdnrvknvHEYO3scddxzNmzePPvroI+rWrRvNmDGD3n//ffrnP/9Z8JInd999t1guZ555Jt1zzz10880304EDB+iCCy6gCRMmhPqdzvCSh5MOs1zZj2M5BZ32PY3lizYTMEWQz6LaehL0kmb7iZ/thI4A4vY7EMRBh9ddQrKfz5au++igJa2dfPLJNGTIEJo6dSr9+Mc/Lnp/HgXl0KFwYw3b6KqrrqLp06dT3bp1K1139tln05133knXXHONOGoQVBpDXKVBbYnIw+uNwjm+udtQk24tJurj4yKfV9hxxgHSHNe7UPFDBl4/R3PiaufyO/Z3Xoaf1QnDLZrh9QLBuNB7FDVQ637/Iwfv7du3U9++fWnduq9XJn4m0cnbRDu9evWqch2Pez506FB644036MMPP6Rvf/vbWv+m7SuKQl+UtA/x2qjY4WV50mrcOBCgvxtsHdfbS5D2qTjHBneO/R22rSvK9iNsODKR+lrkzzY9f5u9HiIsq+9R2LBdqB1F13sfudXk1ltvpQ8++IDq1atHt9xyi5hMh9tONm7cWPACRNu2bROLgXu+g+BKhbNqmLc2DDn5ClaCX38eeHl47ZCo00mr00onCTtLYCNZFdZ9lCbOEzX9tp5E6fWO0hNrS8tJoQAG6Sx7v+J6j3T93sjB+/nnnxetI4899pgI4T179qSOHTtSu3btCl7C2rdvH82fP59uv/12UTHm38V/ny88i6Yfe/bsEfctKSmhhg0bUpMmTah79+40ZcoUOnjwICVh69atNG3aNDrrrLPE8vLLLWjKlaw8hGg79cMtD/W6bUBwqNR7xyPIjlhSO23O6ht2msBEzs9l3OuZuHZKC1XTo/R6u1WBiz3GNsVeR1ZeJ6Rz/kJNHRPo8Bjew4YNoySsXLmSBg8eHPrxmzdvFq0xmzZtEv/mMcj5RMdVq1aJy8yZM2nJkiXUtGlT1xYZvq8f1atX9xzbnHce+MRK/l2PPPIIxcW2w3rFoL/b/3IKujFPoreboc0EbOWn5SSouKek599ZqOe7WK93XCekmbxt8vt60HqS7vK3WeSKd6tWrahWrVqiZzkpHIr79+9P48ePp1mzZlGLFi18PY5P6OSTQDl088Q/ixcvpr1794ogPHv2bGrUqBGtXr2axowZ4xnauaXGz+X00093/R1cUedKPQ8xyCH/1FNPjbQssiYPXzqTpFXtBrCt6u3W4mfLd8Pte56XWS3dxm32GoouSl+xCbxeG8ZS1yvqex559/373/8+3XfffaJafMYZZ1DcevfuTTt27Kh03Y033ujrsTyk35o1a8TPc+fOpR49elRUp0eMGEFHjhyhUaNGieERuerN4V7VvHlz0R7ih1vFnMc5Hz58OC1atEj8Hg7g4K3QiAJQXLHRDSSM2w1QldsJWnKkE1Zo3aQGdL/rsDir3oXO/+DnmsYoJ0lUvU0aQs6EfnSTjzTY+l6Eec6RU82vfvUrmjNnDl1++eX0wgsv0FFHHUVxilJZ5+DN+vXrVxG6VSNHjhRjbPPJn48//niV4M394BdddFGov3348GER6p999ll66KGHaNy4caSD11ThciVq+hfNzwdYrdbICVgwJF30dpO4q93Ov6+2mdjymQRwkuG72DrI635uQVxtOWFxjHbitSMuw3ca4to+6QhwhY4EOHdS0tzOBn2tNmUCnYK8n3F/xiIH77Vr19Idd9xBV199NXXu3JkuvfRSMVEMt20UwjM3JonbSZYtWyZ+HjRokOt9+ATNgQMHimDMVWlduJLOQfvJJ5+ke++9V+yk6CBX1s4Vq42h1Pkhdn5J0N9d3JEj/ofpjDqWcLHqnAwQaDGBLPK7oXaOn1+oV1x+V+Ls+Y7S750nhY5qeB0hiGsIRqZ7fGoTw/frMVa7/byfLKnPf+TgzScqcmBlu3btEqONmDiBDo8zzgGYdenSxfN+8rbS0lLR0hJmYhsn7kXnfm6usnO7yhNPPFHpdh4JpkOHDpH/TtZ4HaKVlVPTVhxp4uXx9JIPq2ywC21ow2zY1Sp2sYo6QjfYLsp4wE5yoy6r4IUm5okjfBeaYCfNlhOdCr1XbhVP52tV7+P1/uheVkHaYXR+HtXfk+Vt6fwirWHy+iQDuJYG2qAT4qQxgY4cM5u1bt3a837qbfwYHcH7rbfeEv9fvny5uDhxv3ecwdvEvVvb+7pMdEH/k+jZpR8V3dCGbTGJUsXGaCZgK51hx9mCklb4jhIo5bZEZwDUsX3yE7qdy9stjPvpyfdaVkFfS5g2kTxmhLD8nI+RRgCPPKoJV5HDXJLGY3dLPISgF/U29TFRvPTSS2Jnw+tSrG989+7dlS4HDxwoeoa9zWeqq889qRkVbVZohem17IIuU52tI1lcwQPonphHftfiGOfba+fb7/jephVLwoTuqJOLRR0VxqZlaOPfnh8gdKvUz0FcOSpy8Ib4tW3bVkzyIy9//N29Be/v9kEz7Ute7HlhNBP9G1p5CbsjEzZ023jOAUDcO4tBw7euAB73rJZhRNk+hQ3dOvj9vbYM7WfycwtCx/seNHwvXl95tL1CchO81ZM9+URLL+ptxU4QTcqWLVtE/7y8/Ojyaz3vm+Wp4+WIJhBuCml1Kumkl7vNfaMAYcI337/YY/yGb93V76izWsYhrjG0o74euU31u6PiFbKLhVr+Hc5LHgL465r/js6drbgq37kZJJkn+lGna/eauIZvc3tMmho3biwuUt/OjV0/rGFmLEyT22twtplAtqDNBPL2OS/WI66ODV5IHH3fxUY5UZ9jGF6vK4mTEqPMA+Hc9jjDt7rM3IZjDBom3cKiCSe7mtr7Pd/jsxnHEY44htuM/Oxuu+22QPevW7euGOv7lFNOoe7du3tOq65bp06dxEQ53F/OQyB6DSnItzGeDVPHiZVJcFsR84cl6okfaVG/NM7XhbaF4nRumKPuyOH9giwpdGJhoXWrnxM0/UxJr3O872KjnMjnxIKcbOYM287XVChQ+tlG+Qm1YSqUzrDttWzdlpnb9taPQqOoqCdwsrQCeJy54fUAOyhunyt1+ajXx6HQ+xv08xb5GU6aNKliOMGgONhee+21YuZJDsVx4pMme/XqRa+++iotWLBADPHnxCc6Lly4UPw8YMAAsoE6lndeoG0h2c8D2kwAKgsTRAqFbz8jnTi/jzoCeLGjpG4VWD8KvYZio6cUCnpBQneQABZkXoNCOyxB+HmezpE2shS+Xw8Rut12TpLgZ8SfBnX8T+4YOe3yRDh84ZP+5Cgdbdq0obPOOktc+Gd5PVe6ecr3kpISqlWrFn3xxRd0yy230IUXXkhJkLNFLl26lFasWFHldp6Bc8OGDeLnsWPHUhaZcuJEsTaTYmyo3KdBZ4VZR7VbXUnhPYO8K/QdCBqqnL3fYb+vcvK1YtxG/vC6zs/vKrTOd9s+BNl2+Q1k6gnnQXZevO7rdxsWdOcg7lE2ks4Or4f4XWkPuCA/5873IExLU3UdQ+XxBDD//ve/RbD9+OOPafPmzfT666+LC/+8fv16uvjii2nnzp30ne98h9555x0RurlaznhGx6efftr33+Tfs3379oqLHJ6QT4xUry8rK6v0OH5+HPp5J2DYsGG0ZMkScT0/nkM3z7rJuA3FOV28rdy+qKaEbzfyA+wcfcOm3vW0xFENwUgmAMmG76A7z7oCuDryUTFuQ+/pHo5PnuAX5EQ/PyHI+TqjHC1Ql1Wx1+M8aTLMMHfO35P0iZg6ssPrIfrf0w7dTlGXe7XyiLPZzJ07l4YPH06/+MUvxHTohVx33XV03333iZD9/e9/X1zHLR9TpkwR/37mmWd8/c327duLQF8MB+3p06dXum7Tpk3Ur18/8X/ZgsLBe//+/eLfXbt2FYG8adP0K6o8ZjcfSeCRTNSTK+WHV34g1ZDqHG+50MQlaVUgvb546hfMLXjza8KslYWXKy9DnkRHR493lNFM3D53qHgD+FsPsqgnBoZdB6jBPYlRpHRNrlVsuakhWfc5MOpykjtOcZ/0p/6tpAsyftflUYP6fANDt5f9e/fQpCGnueY17RXvqVOnih5vbhkp5le/+pX4//33319x3VVXXSX+v2rVKkoCh/b33nuPJkyYIKaH5+fObS/dunWjyZMn0xtvvGFE6A4ryAoljcq3ydV28L+yVy9utzkhdAP442eIwUKiVr/Vtgu1OhznuURRK7aFgq2zuq17JlA/lWmv5xaV2wRAzkscVfEg09yHNd/iSQCLifxJ4BDLvdt+RgDh+/B9V69eXXHdcccdJ8bL5tYQv2S1Oiz+e7feequ45L2nV3450g5GWf6S2TiyiVe1u1A1R8LJrwDJnWzp5wTMMOsC52N0nFAY53BtXqGbxT37sdsoJ2F+h0rXco7rBE2vEy51FtcaxrCz4rUDmeQcIZEr3gcOHBAtEYUmpZH27t0r7suPcapXr17Up5I7fr5AbicDJFWJdvbmFfsbapuJGwxN50/cKxDn547/rV4AQN/Jll5HkfzQOfGOrBbHVQUPW3zxakdIKnTr+P1uFXndRxjiqIA7t+m6csT8mFpMnMvZ6+hO3CK/shNPPFFUvR966CG6/vrrC973d7/7HR0+fFg8RuJ+mD179tAJJ5wQ9alARGGr38W+bGG+jF4nViLYxU/3iaxpH00BsJVc38nqt5OfcKJz7G/n42UVPGphIGzV2ys8JhW6JblDErTo4dV37gzfOospzgp41O1qXDNP6lboM+E8Ry6OIztaK94/+tGPxCghPBb3nXfeSV9++WWV+/B1v/71r8V9uKeaHyNxTzXjfmsIhz8gzrCUxF6brmll/c7aBv7oODLg1WaCHR8ASnzH1HlUyVkNL1YVj2PaebfKYbFLkKnXiynWO51U6I663S3Ud64eZdBNHYXGlFbP+TH0wwcZwUb3OQCxBW8ezYTH5uZKNp88ecwxx4hRQ0aPHk1jxoyh7373u3TsscfSzTffLO5z9tlni8dIctSRc889N+pTyZVCK2u3D06QL1YSJ044BfmioYLqDcEYwD5h1mlu7V1+AniaCoXvMGGrUItJ0oK2iAR9nnG9LlPC93zNoTvqSbVxfo4iB++aNWvS/Pnz6bLLLhPVbO7jfvnll2n27Nk0a9YsMc43j6fNt/E42fPmzROPkX7/+9+Lcbn58RCPMB/ktEYfydssnKZBmwmAfbwq4V6ijPcdd0D1EwBNaTEJE77DjCEedx9y2uF7fgyhO0oFO64+e0nLPO08Fjb3b2/cuJF+85vfiEo3T7nOF/6Zr+MZIR9++GFq0KBBpcc2bNhQjFVdo4b/6TYhXXGH8rRWmllvPwryWCe0mQDES9eRPGcAj2vCnTjW8V6zA7qdeGdSi0nQ0BYmEAZp2bExfDeMIXRHEefnSOtpo23btqVrrrlG568EjYIOJeQ1XJDu5wQAAIWHGAyq2JCE6pCDafE6iU2G7yCT0Oh8HerOSNiKqXP4RZ3VeLffb9LQjmmPYNJCY2iO40RLLRVvMI/XdLZxDwkYRqHnhqnik4M2E4D06Sx2mNx64iccBZ2EJmrgUpdD1JNRne0hOp6f1++3tSA2X/Pf0b0s4lrGCN45E/aDrmsEkzDPRZ0qHuLdYKPNBMCM77K8JN16kjSvUOM2C2Oh36EjdMtlIZeHjvDtbBExPRj6affJcouJk9/f95cXPiS/Ar3i2267Tfy/efPmdMUVV1S6Liiesh308HtoKI3DSEG+dDpWnqB3BsuwMPIMgN7vUtTCh5/WkygzXUZpmYhCR+gs9BzSWC5h206YzrYInTNdxmWmx1jocfwdr2Ub9DNYrZwH4fapevXqYnSSk08+mT744INK1wXFQwtCYTzLJ594ypMMNW7cuMrtckXMX45nl37kOrmB84PiNuV3Wpy9XW7BW74GteKNUOeNPxPqYVm/OzPOw6uqMEcb8B4B6OcM32Gm7PbTtpF0EcRtW+VX1OdaaN3n/DvMxPDttRMRNYTHeaR5vobe7qTfE/Wz4gz8h/bvpaX/73zPvKYK9Kr79OkjQvZxxx1X5TpIl/OD4DWblilV77BfOgQ6/5WtoLxCd1B4jwDiP/nS63tWrELuZ/2gbk+SCjVRZ36Mws/fddvGmkidhVFHFdz0qneLBHeEnEcXwv7tQMmHx+T2cx2Yz7Qvkw0rtLwy6XMCkHd+d2wLBXAZvgsVP3RPNR/XKB1Rq91B/p5cJqZXvXUFcFmo050X5mvoH0/rfdDxN/WO4QJWSLvq7fWlc2szgeiK9SYGXc6oaAPYwyuA+w3fLIkAHjR862oxCftYG8K3jgAeV/huGKHNxPZCXWyjmhw8eJCee+45uueee+j++++n1157La4/BRaOoe33Syf7uyEcvyvYtKeSBoB4ee0w+1m/qusHOeReHEMQJh1mw6z3bF1XRpl8R+fkOvMj/g5beu0LCbzLsWfPHnr66afFzyNGjKA6daq++FWrVtGwYcPoX//6V6XrzzrrLHrqqaeoRYsWUZ4zBORWQUir6h32S4d2h+S5fT5Q7QbIjmIjnRQLnGolXGcYKlaZjaviKdd5fgpDNlW9dYyCorPy3TBEtTup0UuMrHgvWbKELrroIvrtb3/rGro/++wzGjx4sAjdPGCKelmxYgV9//vf1/XcwQcTP6DOL53th41s4FWZKlaxwg4PQLbHBy80xncxznGvdVTBnZVZdQIaeV3UcbELjfjl/DlLVW8p7PKLWvkO+7iZSpXbxEwTe/B+9dVXxf9HjRrlevtdd91F27d/vXDHjRtHy5Yto3fffZeuvfZaEb7feustevLJJ6M+b9BQfUxykHxW6G9l4ctk6mei2EbC9o0IAEQTJXzHNRGPc+IZGcB1TJRTbLIhvvhZFrafh+Sc4TrO8O139tEstpZEDt4rV64UwwcOHDjQ9faZM2eK24cMGULTpk2jHj16UElJCU2ZMkUEcQ7fc+fO1fHcQdOKIMnwrWuWKijOueFwfjZs32gAgL6ddB3hm8UxC6bu2R8LVbv9XO/2O2wWJXz7yQ9hQvdM5UhHlkJ3qOD9ySefUM2aNalz585Vbnv//fdFqwn7xS9+UeX2q6++Wvx/9erV4Z4tuJJ752FWBEkF4aDh3m2ljf5i/5wtIs7Pht9JIwAgP3S1ltm4XnG+dj/Lgl+n7QWMsKFWTi3vN4CHyRotMtJaEjl4f/rpp2JWHp6x0q0azmrXrk1nn312ldu7dOkiquHbtm0L+3whocNJcXD74hV6XjauvE2n9mAWWr4YSQYg+6KMdGJjKA3zfNByEiyA6ziiPtOSsdITC9481TtPZe6G+7dZp06dRPh24kp506ZN6csvvwzzXCEmcfd6FxsrNstfMJOoQdvPTo2z4oMjDgDZp6vlRDItfPttM/G7LLJUIIpSoHNWv9UqeNAWk6wLHLyPOeYYOnToEK1fv77KbcuXLxcV7e7du3s+vqysjBo0aBD8mYKV4TvNscKh8AlQAADFAmfUEy5NUeykyjC3ZUnYMb7dqt/OS9jnklWBg/fpp58u/v/II49Uuv6jjz6id955R/z8ne98x/WxmzdvFhPrtGnTJtyzhUj8fKHiCMpR+sjR8gAAEB+vI1lydA8d62FTqt5RdgSKVb1NeY1RZD3wWhu8L7zwQjEyyb333itmpfzHP/4hxvb+r//6L3E9V7N5RBM3r7zySkWvN5j3hZJ7p87DRGHIxxcK3X73rPNScYhjQxpmg4mdHQCQoobvOEY5SZrfbZDNr9GEdo+ZGe/tDh28OWD36dNHtJvceOONYnSTAQMG0Jo1a0SbyXXXXUeNGjVyfeyf//xncR+3Ey/BHGEPD0lBAnsevmRpibLDgp0dAPBqPXFeTG858QrEQXYmio3tnYUdDIZtsoHBmz377LN03nnnVZqVkl1yySU0YcIE18dwK8qCBQvEzzyzJdghSOXbeTJF2PBu+4ori3BiJUC+v9/qzrj8OWg1PM11u1f411lkyEr4Zs4ZQ5P4e3kRKhk1adKEnnvuOfr4448r+rr5hMp27dp5PqZWrVoisPP/O3ToEP4ZQ+QPt9/qgwzOvFKVgdprJRV2VipTKyQAAOAvfPP6n7cThdb/vE5PI1zpDsGmvs44qt5J70C0yMkR8Egp6YQTThAXP9q3by8uoL9S8fr6neJnOYmOV58UXxfmi6QGcCe1Gh50yKAgXzJUXAEAzCTDtx+8DUo6YOkq5CTxOp3b6LTDaJBiXZS/kSehWk3APEn15HqN0Ymp4LMBJ1YC5JOO4kax9YdJRzKjrOv8vs4ghS51cjPn70irdUXHEIPFzMzotPCFIHhD4mN0+oUQmA6cWAkAYdcbJq23ObAWCvtu67piOyB+149+dzKcswmrcy2oP6cdvm3/GyZB8M6htA/rBPn7CIHRFDsT3w+0+QCA3/WNzWNeB1nX6Tih1K3C7Sbt8B1XbpiZk+EDnRC8c8aUD7nX8zB1hQwAkHU6drJNKZZE2ZbornrLqrbzIu/jpzqeZviOo+VkZs76ulUI3mAck/oAbYUqNQCEXXfISxTFKsJJBMgg25Kgr9fv0US3YO1sKfHLhPBt+u+0AYJ3TuV5bxPc8UbElGoVAKSDR8mSI2WFDd/F1iNxF1eKBVM/gVl38SJs2Hb7PWnhoKwjO8zMef5A8M6hNPcy8/6FS4tJJzwBgLmcgTPOo2dxVm6LBVRdRYY01q1p98lH2Y7PzOEoJk4I3hllcq90nr9waUAVGwBMk4WWwrTXrba1nCB0fw3BO4OysEJDhdaM5Y1ecYB8ku0maY6olFQYles5tc3Ghqp3moJWvRG6/wPBOwPUcBRkBZB020fQoYNkNQHhLx/VGwAwg+3tJoVCaaFzWcK8zjyuN8OOcoKj3V9D8M6QICsAE78AxSY7AACAdIQN38XG9DaZ6UWfNHu9g2QInNtVGYI3JAJfvHQVOuyLth4A0N1+Idc7fugKjzqKN2Gr3nldj/K23W37Lq9Hi0lVCN45l2QgNrHKnmVBNiB5PFwKAMWpY3q7taBEHfPbtKp3oR0NUyvgaVa9C7WdyNtN2vaXFpjMKCk1E/1rIPz973+n0047jQ4ePEjz58+ngQMHprJk+MuQxAcO1W5zcFWmYR187QEgHK/wydd7BVY/6x3eFkUJaH62ZW793UHCdLHKfxrrV95xSXsbK7OE+jxMCtvFdvaSXn6oeKfgZz/7GdWqVcuY3q8kPnQmfwnzwq2qXejwqKnVHQAwR7Ew6udomq4ZGXVUz6NOGpRWy0naQwirlW1Tt/elHstI5qWkliGCd8L++Mc/0ooVK+i///u/yQQmfkHy2iuX1vJFmwkA6KCj5SRM+InzxPwwLXtJb8NMatcxMVOonyuvZaVr588PBO8E/fvf/xaB+8Ybb6T27duTSeKqevsdQtC54kQYjIe6YcAODgBEVajFxOYAGfcJpWBO6E7684fgnSAO3A0bNqQbbriBTBLXHqqOMI92B/3UDQM2EgAQ93o66KgfQaqOfu8b5O9HPWE0DWlPI2976FbFvRytC9779u0TJyTefvvtNHToUGrXrh1Vq1ZNXCZNmuTrd+zZs0fct6SkRAThJk2aUPfu3WnKlCnihMc4cHvJ73//e7rvvvuobt262n9/1JUEh+84qt6mHnbKC6/PBW8IEboBIKl1jl9hDvn7DVVJrfPyPLygiUYHCN1JVL2tG95g5cqVNHjw4NCP37x5M/Xt25c2bdok/l2/fn06cOAArVq1SlxmzpxJS5YsoaZNq648ysvLxX39qF69OtWuXVv8fPjwYXFCJT/v8847j+IW5cxqDt8mHe4DAIB8kSN1FBvpBJOuuS8TFLwoctU6zuVoXcWbcSju378/jR8/nmbNmkUtWrTw9bhDhw7RkCFDROhu2bIlLV68mPbu3Suq6LNnz6ZGjRrR6tWracyYMZ6hvV69er4up59+esXjpk6dSuvWrRPV7ripe/RBDz3p/JAFnR4ezGLbYVYAML8fOkgVuFjlW0c7QLH1nHxdcczaGQcUzaq2l4RZJnGfaGldxbt37960Y8eOKr3TfsyYMYPWrFkjfp47dy716NGjojo9YsQIOnLkCI0aNYrmzZsnqt4c7lXNmzenadOm+fpbsmK+a9cuuuWWW2js2LHi78hK+/bt28X/P/30U3Fd27ZtqUaNGmSCqFXvtMcUBQCA5PgJplwUmv/R19s9XZXvJIImh++ws1kGfb26oOpNkT8bcY6Pbl3wjhJOOXizfv36VYRu1ciRI+nmm2+mjRs30uOPP14leHM/+EUXXRTob+7cuVP0lD/yyCPi4iR/35YtW6hNmzaUNnUg/DAf3DDTw6qHC9EXBwAAhcJ30Eqk28Q5QUZsiTJqS9KT6pgwoU6aSjVWqeWy1H303rrgHRa3kyxbtkz8PGjQINf78AmaPIvkQw89RIsWLdLyd4855hh6+umnq1z/4osv0v/8z//QhAkTqGvXrqKaboqo4TvqhxQn/QEAZFOYIKqGb/nvpMhqd5jwnWbVO89Ga/586D6CkJvgzT3W3ErCunTp4nk/eVtpaaloaWnWrFmkv8snb/7gBz9wHdObceU9rSnjdYfvPO9lZwn6uwEg6jrELahGCaLoXw4mj+0mpTH0ZMdxBMHKkyvD2LZtW8XPrVu39ryfepv6mDTt3r270sXvyCpRyS8tf+jUi64Wk2IQALEcAcAOuibRSUrQ7UuUqeRxkmVyRsd0NETnlPK5qXhzn7Vahfai3qY+Rjfu7fbbL84nXqomTpzoe8zyqD1KzsfJKrif+wIAAKRJZ+h1C99BRnRJstc7b0pjnPRGhnk/Q1z6kZuKt834xEseHUVebrrpptQqxOoHjn9WLwAAkE9+tz9pVIDjPG/I74guScvjTJajY+79V4cZjLJscxO8eYxu9URLL+pt6mPS1Lhx40qXOnXqpLpi0xm0MQECAADkQRqjduUhfJcm+BrVscHD/t3cBO9WrVpV/Lx161bP+6m3qY8ByAP01QMA6Bel6s2BPUxoz9MJqaMTfq1RwndugnenTp3EBDZs7dq1nveTt/FsmFFHNAFKZKxVKA6BGgDyxius6l4fxtlmo94f81yYVdEPG/ZzE7z5pMlevXqJnxcsWOB6n/Lyclq4cKH4ecCAAYk+PwAAAJv5PdHQGUDjDJQmFnOKvV5Z4ZbFKPUSdFnlodd7dIqV/TDLNzfBm40bN078f+nSpbRixYoqt8+ZM4c2bNggfuYp3gEAACC4oFXlvFRz5Y6A2+tVW0pk0M7zssoqK4M3T8O+ffv2ioucGIdPjFSvLysrqxK8S0pKRGV72LBhtGTJEnE9P55D96WXXloxs6VzunibmbrHa+JzyjO0owBAkrIUIIOsP9XwrV7kbYWq9IWCe962t6WGvCbOWJ+WHcx28OYp1o8++uiKCw+3x+65555K1//85z+v9LiaNWvSc889R+3btxcnUZ5zzjnUoEEDcRk+fLiYnIZ/98yZM1N6ZfnjdYgIIVA/LFMAMKndRIpjfOu4Qz2/ziiTBjkDdrHA7XxsEFk+yXK0ha/NyuAdBYfu9957jyZMmCCmh69WrRrVqlWLunXrRpMnT6Y33niDmjaNfzzsJMJVlqoJAACQvZ17uZ2KI3y7BdQ4CxBhfneQwO0UZBtv6pHvrBh+zkm+72vlNEqbNm2K9Hgen/vWW28Vl6ziL/L8j7an/TTAMLxhcKvSoBoOAEmzeTvltS41fdnpmHnRBKUW70TkruINAJUhdANAliR1tFeuO9UAbvL61Ma2jCy+HgTvHDF5DxFjeCfH5A0DANgvSCU4rpmWTRxGMA5hlp3JWSAPELxzwrQ9Q0wVn374lhcAgLjWMyYF0biej/P3JvW6w+xcRJ3u3ASlFj93huANxkEYBACwkwnr77xUu7NSiMvba0DwBgAAgNTCd1ztJnnZAQm77GyvHNsKwTvjnF9IfNEAACBOYUf7UCeSCfv4vFW7w75eWyvGpRnYWUDwziC5p+38Qtr6RQMAgGyT41nL6nfQAG5ixdyEtpssje1d+r/P1fYsg+ANiXN+0fNYpQAAyKqoY1wHnRZdnW690HMCdzaF79GWh26G4A2pyMKXBwAACgsbeP1Wv/2EbrB/lJNSw59fEAjeAAAAYKRCgdqG0J1UpT1Kq43p4bs0Iy0mEoJ3Dpn65WI4HAgAkA3OdpMo63dZ+VYvpofupOhYBqaG2tKMhW6G4J0D6p5wlj68AACQfbLtpNh1XlDQsfNky1LNodttx81rhy5ONWP/C5AqXjHN/2i7Me+CSV9qAABI7gTLKCdcMlS3C+PQ2LBOTS3b6RaN6kT+PVGfg+7QrX6GOBepIVv9bMnMpGNZukHFO6NM3sNH1R0AIB90tpvYKonXrGunxIR+79KYQ7f8Wb1EGVUnKATvHDLpcBKGEgQAyI+oVW+IX5rhuzSB0O1HnOEbwRsAAAC0M6G6bcJzSJKuoJhm+B6tIXTrOAE3rvCN4A2J4S8w2kwAAAD0090Dn3T4LtX8d4KcgFvod+gO3wjeOWHidLp5r0wAAGQd1ut2L4ekwnepxhYT3S2sundoELxzwO1DY0Kftw07AwAAADbRvW0NG775/n4vzOQj4uosqlGXL4YTzNlwTmnx+sJieCgAgGxLa1tkWpU5ieUQ1xDCHIpnLljne6jBtMJ0meZqt8o5FGHY4QZR8YbEmLw3CwAAkBVxHFH2W/m2oYKdZt83gneOOD8kJrSbAAAAgL4qfpxHktXwbVrbSFmM1W6d4RutJjlR6PBT3LNUIdwDAORb0u0mprWZpLEcdM1k6eQWqrkNJasV7mK5ipdzkKWMinfOJfVFcf4d516iqStJAADQj9f5WO+HX3bFJH3+VNqhuyylwRrCLGcE74xLe8VWqNqNEysBAPK9PUp7G5V1eRg9rCziRDk67D1w2Pd90WqSM26HnuTZynG1m6S9JwwAAPaF77AtGaaH+SDtJs7XEmSZqK0QcbScmGRQiqGb//Yz75T5vj8q3jlS7IOJXmwAADCF6QE6blFCd16OLJcZUtE/t2Mz3/dF8IbYZqfCFPEAAJBk+LYlrCf9PE0JqFlrMQkDwTuHvL6ASbWE4MRKAADQeSJm1k7W1Plaoo47baIyS0M3Q/DOGT8fUh1V72K/w8YvCwAAmBdGbQ3cST7vLIbvQZbmCATvHAjy5XYOjh8FTqoEAACdZGVbvWQR93OrPd06xv5Ww7fNAbwswYly4oDgnVOFvnQcmKP0fKO3GwAAoLhiOw7OAB4VB1abQ2uZxTsMEoJ3Dvn90rlNDVtMsfvYvqcKAACQBbaF2DKL+7pVCN5QkKx++6mAy9uCtJhk9TAhAACAqdtB23q+yzISuhmCd44F/cIVCt9hQjcAAACkw5bwXZah0M0QvHO6Rx32A+zWfiJ7uouFbtO/3AAAAGlI6+iv6eG7LGOhm2V7DlGIRZSqdpa+PAAAAGlMJa+TqVPLl6UQunXMFloMKt45l+ZeLvq7AQAA0t8umlYUKzMgdHtdFxWCd454tZskEb5NPYwFAAAAFGmM76CPKfvfv+N1MSF0+7ktDHOOKUAq5CGmpP4WAAAA6KWGw7DtEXIbLdtOmFvriVfIDhq+BxmQCfyGap1tQAjeEDu3LyPaTAAAAMKHvLi2o24BvNh9i3EW+AalHLrDLDtd4btaeXl5eeTfAr589tlndNttt9Hzzz9PpaWl1KxZMzr99NPp/vvvp44dO1a5/+7du6lJkya0a9cuaty4sbal7PbB4S9FXCdVuE2ag+ANAADgrljA8zPjJbjTkT+cy3fvnt00oGt7X3kNFe+ErF+/nvr06UO1atWiH//4x9S2bVv64osvaOXKlbRjxw7X4J1Gy4nu8I1qNwAAQLLV1bRGSDFdT01HCeTvCbOMEbwTMnr0aDrmmGPolVdeoUaNGpEN/d66Qnjah5QAAACyAEeLzVpuYXZwELwTsHTpUlqxYgU999xzInTv37+fqlWrRnXq1KE0FPqgqCFZhvCwAdytxQQAAADiD4uoelMiOyr8Nxa/s9v3/TGcYAIWLFgg/n/UUUfRd77zHapfvz7Vq1ePzjzzTHr11VfJVFGGG/R6DPbWAQAAIE6cNUzNG9YF73379tH8+fPp9ttvp6FDh1K7du1E9ZgvkyZN8vU79uzZI+5bUlJCDRs2FCcwdu/enaZMmUIHDx7U/pw//PBD8f9hw4aJpvvZs2fTgw8+SJ988gmdc8459NZbb1HS/H4gOXzzxc+Ym2mOwQkAAACVmRo+sxa4z+rQNLutJnwy4uDBg0M/fvPmzdS3b1/atGmT+DdXnw8cOECrVq0Sl5kzZ9KSJUuoadOqC5EHgOH7+lG9enWqXbt2RdBn3/zmN0W7Ce8ksP79+1Pnzp3p1ltvFdebrFCIVoccKnS/vK0AAAAA0pbllpOeFuYK6yrejEMxh9bx48fTrFmzqEWLFr4ed+jQIRoyZIgI3S1btqTFixfT3r17RRWdq9Dcf7169WoaM2aMZ2jnFhE/Fx4mUOJ/s7Fjx1aEbnbiiSdSz5496eWXXyabP7CyKo4KNwAAgHlsDKhZfU3WVbx79+4tht9T3Xjjjb4eO2PGDFqzZo34ee7cudSjR4+K6vSIESPoyJEjNGrUKJo3b56oenO4VzVv3pymTZvm62+pFfPWrVuL/7vtIPAOAI90cvjwYapRowZldU84jS8IH52488476aabbkrtRNYswfLEMjUdPqNYnibD51MPNU/YuEwzMYFO+/btRTV64sSJBfu8eRxtPpmxX79+9OKLL1a5nRcFj6e9ceNGUZ3moK7Do48+Spdccgk99NBDdPnll1e6rVevXvSPf/yDtm/fntgEOm7iDt9pBO8kl18eYHlimZoOn1EsT5OZ8Pm0seWkZ4H8YMIyDfo8rGw1CYPbSZYtWyZ+HjRokOt9uA1k4MCB4udFixZp+9vnn3++6CV/5JFH6Kuvvqq4nttali9fXvE3s8rWw0EAAABZktb2OOxJjz0zmB+sazUJa926daKVhHXp0sXzfvI2ntKdW1p4WveouEXljjvuoGuuuUZU3bmdhWet5KniuSWFp5FPW1wtJ1n80gAAANjKa7vsNwM4H1/occ779kQmyE/w3rZtW5WeazfqbfwYHcGbXX311SKA85CFfFIon3DJPeTcm9ShQwfXx8guoK1bt4rDGBL3McXRy9Tl6Bq0YsNOrcPrqM87afJvp/kcsgTLE8vUdPiMYnmazPTPJ2cAP5zPXz6O84NzWL24X+tuQ5ap/Pt+urdzE7zlkH6M2z68qLepj9E1bTxf/JJ/n4cchPDatm2LxacRlqd+WKZYnibD5xPL03RtDdnOc27jXu9CchO8bdSqVStav3491apVq9IwhHFVvAEAAAAgGK50c+jm3FZMboI3j9GtnmjpRb1NfUwaeJhDrzYUAAAAADBDsUp37kY1UfdCuGfai3qbnz0XAAAAAAA/chO8O3XqJCrIbO3atZ73k7fxZDe6TqwEAAAAAMhN8OaTJnmyGrZgwQLPHp2FCxeKnwcMGJDo8wMAAACAbMtN8Gbjxo0T/1+6dCmtWLGiyu1z5syhDRs2iJ955koAAAAAgFwH7507d4op1uVFTozDJ0aq15eVlVUJ3iUlJaKyPWzYMFqyZIm4nh/PofvSSy+tmNmSx9gGAAAAANClWrmf0b4N0759e9q8eXPR+3HQnj59eqXrNm3aRP369RP/ly0oHLz3798v/t21a1cRyHlGSQAAAACAXFe8o4b29957jyZMmCCmh+fxsXmc7G7dutHkyZPpjTfeQOgGAAAAAO2srHgDAAAAANgmdxVvAAAAAIA0IHgDAAAAACQAwRsAAAAAIAEI3gAAAAAACUDwBgAAAABIAIJ3Qt5++2267rrr6LTTTqPGjRvT0UcfTX369KFnnnkmqacAAAAAACnCcIIJGTlyJC1evJiGDh1KZ5xxBn355Zc0a9YsWrlyJd1yyy102223JfVUAAAAACAFCN4JWbZsmZikp27duhXXHT58mM4++2xatWoVffrpp9SsWbOkng4AAAAAJAytJgnp1atXpdDNatSoISrghw4dog8//DCppwIAAAAAKUDwTtm2bdvE/7nnGwAAAACyK1fBe9++fTR//ny6/fbbRaW5Xbt2VK1aNXGZNGmSr9+xZ88ecd+SkhJq2LAhNWnShLp3705TpkyhgwcPBno+W7dupWnTptFZZ51FHTt2DPmqAAAAAMAGNSlH+ETGwYMHh3785s2bqW/fvrRp0ybx7/r169OBAwdEjzZfZs6cSUuWLKGmTZv62gm44IILxOMfeeSR0M8JAAAAAOyQq4o341Dcv39/Gj9+vBhVpEWLFr4ex33YQ4YMEaG7ZcuWYoSSvXv3igA9e/ZsatSoEa1evZrGjBlT9HdxZZwr7jzEIIf1U089VcMrAwAAAACT5ari3bt3b9qxY0el62688UZfj50xYwatWbNG/Dx37lzq0aOH+Ll69eo0YsQIOnLkCI0aNYrmzZsnqt4c7t189dVXNHz4cFq0aJFoM+EADgAAAADZl6uKN48iEhYHb9avX7+K0O0cp/v4448XPz/++OOuv4OHD+Rw/uyzz9KDDz5I48aNC/18AAAAAMAuuQreYXE7CY/DzQYNGuR6Hz5Bc+DAgeJnrmY7cUWcg/aTTz5J9957L11++eUxP2sAAAAAMEmuWk3CWrdunQjOrEuXLp73k7eVlpaKlhZ1QhzuKed+bq6WN2/enJ544olKj+3Zsyd16NAhttcAAAAAAOlC8A4w1jZr3bq15/3U2/gxavB+6623xP+XL18uLk7c7+0M3hz2+WTOWrVqiYq6VKdOHXEBAAAAgHSVl5eL4aZbtWolzv0rBMHbB16YEg8h6EW9TX0Me+mllygoDu8Y3xsAAADAfFu2bKE2bdoUvA+Ct8F4iEI2c8Fyqt+gofh56cadVKNWbapZq7b4d4M64U8YdbP3wGHx/3M7/qdaH5fF6yuPMBPmtcjnW8jzL39c5boh3zmh4mf5Ws/qUHz89bis2LDTc9nofo/DLEOv58KP1fn8gjwXSA6/x/K9UX/2ksT6A9JdZwF4bcOCfi5XOLZ/Ntpbtocu6F1SkdsKQfD2QV2QfKKlF/U2Pwu/GNle8k5ZHapT/nVrSeOjW1DDOvG9bXUbEJUdOESvlB6kQSc2j+3vzP9oO9Vp0DDya+Hny2YuWOd5n1r1/vdOigUrP6n4ed4bX7cSXdD/JPF/ft09Oya3QXt9/U5q0KhxxXKRz6l69f+0F/kxemCnSMsw6cd6/b5C7yUkw+2zpL7Xhd53uf7wg79r8jPvV5zrJdus/fzrHaAk11cAvM1yktsw9XPpVwPlsbZT24K9IHj7wD076jTvXhPe8G1uj4mKK0x1YwzbThyGeePJG8Q4NnJyQ6trB0INai0aVe19L91zoFKgcAa7Vk3qivs8u/QjOnKkvOL6pDdmvFx4mfDz49Dt9lq88POXrytsAE+an4AdZBnooH5WiuHlzN+ToLw+93HucPhZjuprV59L0M+T3++1XMf4eZy6nNXHIIQDpB+6IRgEbx86deokmuX5ZMe1a9d6DinItzGeDVM9sTKqO8YOpmrVq9M3up5H3zj9vECPDRvC4gjffjeycQVuxq/p/H4nVlzPYVvej3+H/Jmfa5LBWw3d8rkEIe9vegB3hsukg7WfgO213JwhW/47yPeD32evsK5+LsNy+17xMlc/4168vkNuOwQ6PltB1gFu93UG92LvRZD72hiGUPWGJD5nEF21cj4VM8fat29PmzdvpokTJ9KkSZM879enTx969dVX6bvf/a6YmdKJF+MJJ5xAGzZsoLFjx1ZMuBPF7t27qUmTJtTv5mepZohj+s5wEXZjqQaFMBsrnYHbT3BzBm5n0JGvQX1eHMCdISiJlhNekcnnIZ+DjjDqttORBrfQlkTY9lO59huw0wxr8rPh5+8WatmQ37tiO6s6jwok/bnzc/RBLoeo6zRTIXxDHBC4i9u7ZzcN6Nqedu3aRY0bF26dQcXbJ578hoP30qVLacWKFXTWWWdVun3OnDkidDMO3iZQN6yychVmY6hurMJWwHVWuP0EBrXCzdyesxrAOWir4fvpJR8mujHTGbrV3xPlfY8qbPVeVxgM+pr9hrFCnwd1A+V2P3m7vK3YBi3Id83rvrLKzt9BuUz8VsHdFLu//L1R2lWSqqCrVfMsBHBUviGOzxTolbuK986dO8XU7dLpp58uhn/hCW5uuOGGiuvr1q1LDRt+PZIIO3TokLjvmjVrxHjdXNHu37+/aD+ZO3cuXXLJJaJCzW0o8+bN0/JcZcW7/jfaiFaTtmedLy5hyQ1ilI1gkEqRrl7uIFU69TU6Q3ehIKRWnWXLSRJVb1nxln83Djred9NaSaIGba/qaLHPdLHPgjNYu91e6LYw1N/n9Tuc1XC3Cnhc71eQvnk3SXxuw3wenIHdpFYWVL5BB4TueCreuQvesrXET4V7+vTpla7jyWz69esn/i/H7ebgvX//fvHvrl27ijaUpk31BLWorSZxhrBCATyO1pIgJ4f5Dd1u4Zsr3XI0ET7R8pHLexZ9vI7Q7adtxus+poTvOCvcOtqmorQXmBZk/FQ2i4Vw53czySMUNrRMRTlx1u2xSQdx0z6zYA8E7uDQahJjaH/vvfdo8uTJ9NRTT9HGjRvFrJKnnHIKXXjhhXTVVVdR7dpfj69tKvUEwijcDtW63R5FmCCgbpz9hG6VvL9a9Y5rZBcn52t0CxyyPSBMMNL1vnuJK7TpCF5ZCtxBnhffx20DKofwk+0nzhF/wn7G4pDmScNR1mHOx3qtJ93oWt+g7QTCfm4gXrmreNtEVrw7Xz2HatSpr/3EOxNHvQgb5MJWu4tVveNsN/FqMyn0/kQJuLrf9zjbFKIEbl3VRlNDd1BhK9/MlACu+4TxtMjPpleoj6NSnpXPMcQPoTs8tJpkLHhPev5dqtugkdaqIm/ATN1o8esM+hrV18Mbr6DVbmf4lu0fMngH/V1+3LLgo4qQr1b2WKH3JurnIOp7H3cwi7KDEGaIv7yEE7eNaqHzMExsPbGtgBCVrtFX8vD5hmgQupML3tUj/i1IwNQrLqDfXDyAOux/U/xbR8sAb0xNnCEwbOiO0pdZCAfwoDPrBeE2UU6xICFvj/I5CPPe82PUMGZS6Ob3Xe5wIXT7J5eV2/dGx+csTvLzp34us4Z3iOKcqdjUAIgQmPwyh+QgeFvg5w8+TddNW0Q9fjBW+8bQlA1W1I1nlN5uSb0/n1gZZ6XPbUUXdNZEHYHFpMDtN3TLkK1ewgTuvFUCvV6r7eE7DwGccQEgbBHAlmClPk8E8OSXOSQjX7vSGaHrJCh5wl1a4zy7PZ+g4ggEHES4BcT5u5M4Wcn5PhTrCY16sqXfoBLnTkiQwK3zBLQ8he5iJ1wWYuJJlyaOWx8nHXMpmA4BEMvcL5OG7gwDFW9L6apEmbAhjVqlcuvtjhrGLuh/EiXJ+T7KSq76b6eo4UJWC/1cTAndqGxTrDscXq1aple+wx7NsTWAh6l82xpsbX3eprNxuc5XjvqobVhxtoPGAcHbAnsPHI4leKnS2lDFdQKX7mpmHF9sdXZMr9k21aDpDOOS6WEoSuhWX3PYnSr5WeD/q5c8c75+ddkifJvN1rAB5rA1dLud96B+H2z5TiB4W+DP4/+LfnfZIHrpycoT+ki6qt5phe+0Rufw+9zUqeTj+P1unCHTq9Jr6yH1ICdRhj1hEgG7+DJyLudiUPk2Q9gTLk0NXMWel6nP20Y2Lsv5RWbBtq36jeBtgT88tYSeXfomnTFkTJVqlK7glUbLic6gr3M0ExlI+ATLuPDKQY4VHvR16B65JWlB2kv8hEG3ijb447WsCn3G1PBt8tGWtAsKSbAhZIA5shi6VbaEbwRvi8jKn1v4tm2UkzhaTHSeYOH1u+JYcTnfu0Kvw6vqbXIAiiN0O0M2wnZ46rJztjZ54fdPXkz+7MlzFOSIJ1kK4bLKFzRk2Bi+IDrb3vf5/9s6EnRITRuG30TwhtSq3lH+XqGNva4QxidYxhEq3FaAcsZNvzsPNle9g4xeAukIuhNrcvh2G3YwS8KEb4CsVLnT+E7InQL1snj9Dt+PR/C2VJwn2cW5YZJVJx0hP+n+ZlM2brYNnRTm81ms2g3pt5x4tZ6YHMKdY35nqQoeZP1kW/UTorHh/Z6vhFhdlWud22y356degkDwtsAlQ/vT6O99m+b+8Q/i33GeZBdn1dv0KagLBRHdh6+cI5roYHLg0dFigv7tZPlpOXFrPTH9s+g2VKbt4duGw+t5CI8mMn25zfcIszpbsaKMeOL3+TWoU8P370TwtuTkypkL36BhP7oksb8Z14YojtCtBgObqqFyEpuwsw2qr9uG0U2ihm6w4+iKTRVwlpXwDWBL6J7vUT2Og/p7g4TvOJ9f9neTM+CsDk1p7eeHXYOK88OgY3Y5P4EwKJ0bNbdhBONqv+CRTVo1qSue//n9Tozlb3j1d2clbPr5LPmpqsqNSFaWiy2zWvLnUp7k5Jf6/ZSzXrox5eiXabP4hiGre37XhUnMxOv3eUD2lmuxkNswwaM0zplfmfN74vZ843qOCN6WkhtDt6mddZAjAejYCMXZYhLniXhy6njdG6ygh7yCTPNt6rTexSbJYcUCgwkhIesKfdbcdvSjvPemTUOfhfANkGbgdm7bTGuBaugSwJ23JcGspQJawpdOujZCpmxcTeE2hndUOne8dNExyRECt5k7+ro+r85qeJrrijiO9hXi9X1F8AdbmB623aT9HM1fQlAwfIetQiW1EdIdBJ3Ph1973MPP6dwQu+08uT1/Z9h0e+/DtACYBkMH5qPq7cUtYDrXGWlN7hVn+FVfo/P1yao7C/ocgrabQLYkWZxTA7fN26A04ORKi3mtXHUGRdlyYlKLiboxivPESv59PJZ3EkMKZnFD6fdzWGhKeFS70+G23IOMchKFc4QU20+2dA5ZqK4X3daNSY62knYfMNgduuM8KTLLELwt0L17d+rcuTM98MADRYOIaYcok6hWJRVada9g4jqkbdIoEqZ9HiGapMK3JGfGTPozrWu95QzcQYoReR9txW1bhx2FdOkeZzuvELwt8Oabb9IHH3xAV155ZSp/P0zVO68bCx3BNGiV18RhBf2OZIKJcuyTRvhOa4dSxwQ7znHDg4T6MOE7CzNZynUgjnj5k8QOCarc+iB4W8i5MopzFkuV35W/TRPlhH39UVZ0PFJKmGVT6PC/ifyMZALmKhR6kv7cpRG+nS0f8vvvd9ZLXTP0xl35RhUZikGVWy8Eb8vFOYulyu8GJM7Q7dzoJhHeePnyWN5x7URkMYAG6e32gkqXHZL8/Don5kmK+t136892C+C6Q3Ke206wLkh3xwmhWz806YD2s/3jrHQnNXFOXNwOAcfxGtIeHxnV7uyTo+rEObKS1+cq6QDq9V2S16sjkRR7TJTnYNL5G3FB0DYrdKOXWz9UvC2VxsqpWNUlzWqMTSvrQmN463gdafZ5o9qdLcU+j0n3e6tMCqFR+riD8rOeDdrnbUq7idfnzab1e1bYfp6AyRC8MyKpDV+hnkf1doiH1wbIpJYVVLvzRQ3fbpesnXCZJqxfQe4oxbmzhPaSeCF4Z4BXn3dcGyW3cWbjrvJkYQNbbEUZtqpjW7uNbc837/x8LuVY7M4LiyuA5zV8ZxWq2v6gp9t+CN4QWhKHVdNuoZCT6Kgbdx2H4HRMpW4aTA8PTm4BXGcIz2P4jjKpmc2hG8EcoTsrELzBSmm1Vth0oolpYaTYuN2QXXFWwbO2A6uLLX3eCNRmvT82beNsheBt2cyV8B+2h7iwQdTPhsrGMIINcD7eGzWAZ3lHE+IJlKacCJrFCXIQupOB4G3xzJV+JtIBs1aGYSfPMV0WW2cg/h0jVL3zLY8h2kQYwSRZCN4Z4VVBykIVqNhrSKpSGmdvpY7XYPKOF9pMIImhCLOwvstrnze4Q4tJ9iB4Z1iWKpBZei1Za7vREXbQZmIH3e+TbDvR0fONdUS2q5l5rI6jxSSbELwBfOBp4+PcAAZdwfoNQElVABF6IAqdJ13mpertR9Ce3TyG27yHbkgegjdAAtVkuRJ1m7WSQ3QcFd8kwjCq3fkT19EJdeSTsJPw5GkH0PZ2kyjBEjsI+uCEyuRh3Jgc4HCUxRP6wAzFwg76uyHKzi5X5fgzhIAAoA+q3elBxdtyxUY2sb0C5FZRRZADMK/qHdeRmzBtKHlqNylW9Q46nreJ8lbhxgmV2YbgnSG2j2vtxcSdhygbMg4FOnrG0z4hETNV5psM2s7AHVf49jsKionri7jgSCaEYfuOmO0QvAF8klU09ZB32MoET0Mvg4TO6obbEY80q384OgFpDUGYp6q3TnmrLucVWrfSg+ANVkuq6stB2XRZPeIB9onze+nnc56nqnfWw77XbVndQYj7daHanT6cXAmQEJ610m1UExv5aTMpVpVMu1UG4uX2/uo+uoOq3X/6vAt9H2Wftyk751kNzbZI83sz03FOQh53klHxBggYIKIM4ZW3nkxTNvRgBl0nYOJzld91Cthf7W7RqE7F59bmITHDQvC2QPfu3alz5870wAMPuN7uZ0NmY7+jW1U17WnR49jI2Vb59fNZSvt9guyT431nbb0H+a2MJ/F6Tap2t1DCd54COIK3Bd5880364IMP6Morrwz1+KwdykG1KzjdAcTPZwrvE3jROfSgV/jO2novqiDDCuYt8EJ6xasWjup3HgI4gncOq0CQLN6IpXmIz8QAYluVH8z8HPjZuctD1dv2WSzdwj7Cv168DTKp2u2UpwCO4A0Qkd8NhNd08VkLsNjxgyDcxgLX9ZkzcacT9EE4t4ufVs0WOQjgCN4APumY9Mb24Ou3eog2EwgjTPjGZy178hio43zNaZ9UGSY8t1ACeNYgeIORTDs8rGvjHvfrSiKEoIoIJsKRluJMmD6+2JjdeQzdSUirzUSG7ighembGqt4I3jliWpiNGvBMbLMwObja9v5DPumueqc9eytAXnFgjlq5bvG/j89S+EbwzglUKc0hp1FPcsdBx/vvdxjBYlV3W3eYIDlhPyN5r3rrDCeoPCcni20mcQTlmRkJ3wjeCSkrK6NJkybReeedRy1atKBq1arRRRddlNSfzwTTNqppH7JNA3bgwFR57/XOaj9s1mVx7G4d7SVOWZpwB8E7Idu3b6dbb72V3n77bTrjjDO0//68VBHT3Lhe0P+kiqpv0BUZTxcfl7y895AvOqveaDfJbtEAlXkz3+c4dgRbpBC+5cgqxf7mX17wv41H8E5Iy5Yt6V//+hdt27aNnnzyydxVh4OwsR/Tz8o/6xUpP20mAHHC56+4NMdyhuSlVe2OS4si4dtvUA5Tudf12hC8E1KnTh1q3bp1In8rCxsf01oanNU3v1/AtKoxune8dO0MoToPaX5ebNypDyJrJ6FlXdaq9XG0mLjxCsLOkzl1BHD5u9z+ZtiQj+CdM1nf8CTB9Mq1145X1MPtpu0MAQRtNwEAu0O35AzXMnSrt+uuVLsFevl3jm1Y2/fvyVXw3rdvH82fP59uv/12Gjp0KLVr106c5MgXPvHRjz179oj7lpSUUMOGDalJkybUvXt3mjJlCh08eJBMhg0PxMXm1ibIliwc8cuqrFV4TV4WSfZ3Jx26JWfQdhP2OXmFdRm0owyTmKtmr5UrV9LgwYNDP37z5s3Ut29f2rRpk/h3/fr16cCBA7Rq1SpxmTlzJi1ZsoSaNsXJbpCdE5kwWyWk3W6iM6Tw59n0o1ZJrZeK7aTwci/W7iPfG7SR5bu/O63vVAuf09C7Belixci4XlOuKt6MQ3H//v1p/PjxNGvWLDG0nx+HDh2iIUOGiNDNJ0ouXryY9u7dK6ros2fPpkaNGtHq1atpzJgxZAobq5C8UXT7Mtj4WlTVq1cjU4RpN8HREkhTmFCX53aTYofWs3aCJSrp6XG2eJiqhUulOq3zIXIVvHv37k07duygF154ge6++24aOXKkOOnRjxkzZtCaNWvEz3PnzqVzzjlH/Fy9enUaMWIEPfzww+Lf8+bNE1XvtGXxcKvpr8mGlX8awQOVMEj6c2T6uiJOaYUgG9Z/psjKsrL5RN4WBb4nce9M5Cp416hRI/RjOXizfv36UY8eParcziH++OOPFz8//vjjlDbb2xtMdeRIea5ObPU6AhEGAjjo+AzJSxQ2fhdtCJBZCZRpbbPlRcfvSooN1e5C3EZGiVuugndY3E6ybNky8fOgQYNc78MnaA4cOFD8vGjRIkoDgk28glbR3DZC5/c7MXPvle1tQJAveWk3SQvCd3g6W4DibieyudrtttPgHKkkTtlq9IrJunXr6MiRI+LnLl26eN5P3lZaWipaWpo1a1bp9qlTp9K///1v0S/O3nvvPTHCCuvTp4+4JCELJxeZElz9nKCU1mFwDsReK18/n4EgVcE8H9oHc0+6LPQdgGQhkGdHWqOYxEE98TKp14OKtw8826RUaBIc9Tb1MdLkyZPplltuEVPHMz4Zk//NlxdffJGSgGqPPiZt0J07IoWCcJDPQNTPiyk7SJA/xXYG0W7iryXBKzAjSOdzGMEshW4pytCAYSB4+xy7W+IhBL2ot6mPkXhElPLyctdLoXHEd+/eXenCQxhmkc5+4jgV2mDLlSo2SlWXCYBuYXfsbFjPxD2DZZyFgzS+8zavZ3SH5biLQlkK3WlA8LZA27ZtxUQ98nLnnXf6fiz6b/W6oP9JgR/z9JIPybYdCFt2ggAAssCkI6gQLwRvH3iMbvVESy/qbepjotqyZQvt2rWr4nLTTTflqueWdx5sej1ulRfTKgSFQnWQQ/C2vTeQv6p3oeID2k2SYXM1OunlwAEco5JlG4K3D61atar4eevWrZ73U29THxNV48aNK138jj0OEBaq3ZAFus51gMoQpOMR9Qh1Ev3dphWRbITg7UOnTp3ERDls7dq1nveTt/FsmM4RTaLo3r07de7cmR544AFtvxOAgwdX/JwXAAAVgnZyoTtq+EbLivnQVOQDnzTZq1cvevXVV2nBggViunknPkFy4cKF4ucBAwZofZPefPNNUenWGbhM23O1IfDxIW3bDgEWG04NVT/Iy9CCWR5eFbJzhCbsNsa2bVOeoeLt07hx48T/ly5dSitWrKhy+5w5c2jDhg3i57Fjx+p8j3LDphDoZ9SAtFeG6L0GyM46J4vrKN1srczrqlKj2m2H3AXvnTt30vbt2ysucmIcPjFSvb6srKxK8C4pKRGV7WHDhtGSJUvE9fx4Dt2XXnppxcyW/fv3T+GVgamqV69GWVTsxMooVUgAnTC6U/SgFva7jHUAlgPkPHh37dqVjj766IoLjxjC7rnnnkrX//znP6/0uJo1a9Jzzz1H7du3FydRnnPOOdSgQQNxGT58uBhfm3/3zJkztT9n9HibzU/FyIZWGoAsjm7i58gPvp/xBGmE7uA7hmGOVCZx1MK09lSb5S54R8Ghm6d5nzBhgpgevlq1alSrVi3q1q2bmJXyjTfeoKZN9c/Uxz3eH3zwAV155ZVF74uZAs2tGMnD2niPAMyR5XYT3RC69ZOhOWprINpM7JG7kyt59sgoeHxunvJdTvuehRPs0paXahPGvAaArEJrmT44PyfbUPEGo6tObmEVFWN3WC4A4WW1AODnBEtdsA4CKA7B2wLo8TbLkSPlaT8F62CDDGl/zgqdYJnVdhP05Ga/up1Ufzfog+BtgSA93kEnT4FgcAjQ/ygROLEKTIHvbXLhLe2dbBvXOzb0d2MnTh8Eb4AM4xV6XEOpIcyAqdIOfzaJI7Rh+du/MwDxQfAGCMHP0QKuIj295MNcVgqw4QXb4Cgg2Ia3MXFXu9Fmoh+Cd4bZUJHkAGtbf+UF/U+ivE+eA2AyVBjT3dHG8s+WPBaP4oTgbQGcXGn+lMxZmnYZIKv8tF3h3BewAU6qtBeCd05PrrSBLeNemzxGugnQdgIm8LMuse3oW9ACAWQLTqq0E4J3zqG6AwAAACpMER8fBO8cVxjTru4g9Gf3SAR6PMHmIyx5XDdlqV0u6+ufuN8rHCGJF4I3pCrt8J/GRj2NjUJcQwq6QWsJ2My2dZIOaJezT1zvmQzdOKEyPgjeAAmO3MKhNOlgakOfPIBp8lj1DgM72tmD0B0vBG8LYFQT+w7/OcdXTbLiDADe/HwX5Y4ywndhCN3JH+HMUktQXiF452RUk0IbG2xckjnsl0TlOc0NofzbWe+vBPN5fQ+CfAfz2HISBr7v2WozQbU7fgjeOZfWxqXYxDmoEJup0PvitgFGRQxsh8JEcQjfAP4heIOxbO9NdjszPAsbKNvfFwC/0HKipzUC7RHmw0gmyUHwBtC8wlI3MjhsB2C3LLSc+J1EJ0hAlkez/E4Xn3T4zkKRI+lliO1VMhC8Mw7VyXgcOVLuuqKSK0avHrw0Wy/QvgN5ovu7lvWWkziHFMRwhViW8B8I3mAVE3qGvXZm/Gxc0qrC6NgBCxrcTXivAHR8lrNQ9U4TdvrNhjaTZCF4WyCJ4QSTrOZkuXLE4TvL1R0cQQGbuO38RfkMZ3ndFUWhnWysM+yANpPkIHjnZDhB06o5ea0goQoMYKe8rrP84CN5WeypNglOUM0OBG8w8rBkVqsk2DgB2B2+ueqto/Itf4/XxcaAh/VbvLJ8NDVP8C4ChHBB/5Po2aUf4fCcB1T2Icvhm3tiORyHPTwvg7VXFT3q7/fCv7dQ5Z6Dnd9+bA7Z+J7bD5PmJA8V74xSV4hcPeaLCSe4FJs4J+uwoQIw4zsXZX0YdnxvtZpdaD0Yx/jhJvTw5q1dAkcAwA2CN0AMY+SaeJZ4lJ2vLLf/QP7o+Cz7LSA420f4cX4ea0OBwm+wDLK8dU64g+ALJko8eG/fnq89Xtt6F4G0V5jkmN95gao+5Emh9aYzbIcJ03laL6uBO2/V8ULiWhYmFojyIFLwHjRoUKAgvXjxYjr11FOj/EkAq4Kl3GhyxSepQIrgC5AMPy0hUSrXeZyyPutDsoYV1zIxoQUpbyIF74ULF4ogzYG6kEOHDtENN9wggvqnn34a5U+CxfLa380nYjIc9gTIHq9wrCssJ73O5IAXpN1D13qN/54zXKLqDVkUKXh37NiRSktLRaAeP368CNhO69evp549e9KUKVPoyJEjNHz48Ch/MpeSmEDHFFnuJba1Em3CSbkAJn+2neHbzwmUUBiq3l9Dm0n2RAreq1evptGjR4tA/Zvf/IZ69Oghgrb0+OOP0+mnn06rVq2i+vXr06OPPkqzZs3S8bxzJe4JdCAZNle8/e4M2bpzAdmmewZLP+Fbd+i2vd2k0PIuFC5R9UabSdZECt4NGzakP/7xjzR9+nRq0KABvfXWW9S1a1d65JFHaMyYMXTxxRfTnj17xHV8G/8bzBbXyt32jUaQE1VMPmHFlGElAbIoykmUxX6vLn7WT7LPWnfoLRS+3SrcqHrHw+RtVB5oGdVk7Nix9Pbbb9MZZ5xBZWVl9LOf/ayisn3NNdfQG2+8QSed9HWfK5gr7sOiWT7s6jxBxe2EFVSDASDNAoaJJ9L5CfeYTVM/Ez8LeaFtOMETTjiBzjvvPPFzefnXw6eVlJTQL3/5S6pVq5auPwMRoNoJQaE6DmCGLBcuClW281z1RptNNmkJ3jxSyYABA+jWW28V/+bKd82aNWnNmjV02mmn0QsvvKDjz0DO2VAxzuIhvKye7Ar5EscMlpCcvIbQPO94ZFXk4D1//nwxpOCSJUtEZXvy5Mm0cuVKevXVV6l9+/b0ySef0MCBA8Vwgm6jnoB58tCPHeehu7wewrNhxwjA9p3KLK2f/YZpG8OnqSfTZ7E4lKvgfe2114r2ks8//1wMLfj666/TddddJ24788wz6Z133qELL7xQjHrCwwnyqCcff/yxrucOlhzO9LuhsHEoQedrMyF0+1nh667yIXQDZKvdJKkjAUFCdRar3mm8JhO2U3kWKXjfd999op/7Rz/6kRhakIcOVDVq1IhmzpxJ06ZNqxj1pFu3blGfM1goi/2Jpr6mYiHY784NDsEDgCmiVL1NrT7nbScDNARvDtY8nOCMGTNEsPYybtw4MeoJh24e9QQgC4oF0yzM1GnbEQiAMLCTmf1AamP41t1igzaTDARvDtM8gY7fUU+WL19O119/fZQ/CQnJUh9h3gJp0hsYtJmAzUz+Lse1fjYhgIUJ0Db2epsGbSYZmDI+CB7p5O67747yJyGBQKSzSpuFqq8OCKcA6cvC9zDq+tSk4JV0kLax6g3Zo20cbwDwriwlucL3Ey6Kjelu44muAJD9XuM4ZtQESJLW3U0+0XLnzp20d+/eikl03Bx33HE6/2zmde/enWrUqEFXXnmluETFoQqH7JJx5Ei5CLBZqLQBZBnWi5BlJrQXgcbg/de//pXuv/9+0cO9b9++gvetVq0axvMO6M0336TGjRuTDhwC/VYLuE0kymHJIH2IWa2wXtD/pIqKN8I3gJmCrBchOh3Lmn+H7dsMr+XA18dRHDOpzSjPIrea8MQ4559/vpidUla6C114TG8wn66+bPR3f83E0O3VbhJkhAcTXxdA1vF61fYT4KMEy7wNLQjZEil4L1iwQMxUySdN8v/ff/99cf3RRx8tJsp57bXXaOLEidSsWTNq3rw5Pf/887Rx40Zdzx3AeE8v+ZBsZHslCSAMDCtYGVoSAQwL3g8//LBoHbnlllvEjJWdOn1dJeV+5A4dOlDPnj1F8OYZLJs0aUI/+clPqE4dHOrIA9urMcWg0guQLdjZtC90oz0Iche8V65cKf5/6aWXVrreeWJlmzZtaOrUqfTZZ5/RXXfdFeVPgkUBOi9tJm7LSF4ne7xN5Ww3QcUPsi5LO81R1s9+T7bTPYqIrt9lw46BKfi9Rn93RoL3F198QfXr16djjz224jqudrudYHnuuedS3bp16W9/+1uUPwkJMik4m7qxNGkZOZdX2GWGyh+A+aKse9IOYVkPzegjh9iCN4+0UatWrUrXcUsJTwvPJ1pW+kPVq4te8K1bt0b5k2CBrLeZZHHFz5XuoNVuU3eGAMIqNr49mCdMBR3BGKwN3q1bt6bdu3fT/v37K6476aSvD60vW7as0n0/+ugjEcg5fEP2g3SQagw2dOmHDfUCAGADmyvnSfWnY/zujAXvU089VfRzr169ulJLCV/3y1/+kkpLS8V1n3/+uegD5xMxzzjjDMqjw4cPi/72E044QZxg2rFjR7r99tuNH9M8qVYKBD69UI0GyIekjjDqGnvb5rBs67Yv7dYi0Bi8Bw4cKEL2M888U3Edz6x41FFHiTDOM1RyVbxly5b06quvitvHjx9PeXTVVVfRjTfeSGeffTY98MAD1K9fPzEazE9/+tNUnk9cVWa0mZgRgnEoFSD7O6ZJFUYQlgEMCd4/+MEPaNq0adSrV6+K64455hhxAmXbtm1FNfeTTz4Rk+bwSZgPPvigCOt5s2bNGvrd735HV1xxBU2fPp0uueQS+sMf/iCGYHzsscdo1apVxu9hBwnUujcGWdlIJg3hGyA49HnbB8MKgk0iHfOpV68ejRs3rsr1PXr0oPXr14sp5Lds2SJOuORKr65pz20ze/ZscWTgmmuuqXQ9//s3v/mNuN3kFhwO0ugTc5fX/nTsDAGACbgan9f1cDHYbmd0yngvPKwgh+0LL7yQBg8enNvQzbii/Y1vfINOPPHEStfzUYFWrVolXvGOS97aTIodOThypDy13vUowwkCQPZ3qrnvN2ww4wqzjVVmdZ1o6hFB3X3w6O/OUfA2DY8tPn/+fHFC49ChQ6ldu3biZE++TJo0ydfv2LNnj7hvSUkJNWzYUFTyu3fvTlOmTKGDBw96Pm7btm2i190NX2/LEIt+gnXQNpNiGzcbwiO/ZnXZ5G0HBCBrbDrZ27n+iYs6kY5Ny8dENu60gD65Ob2YZ9nkyntYmzdvpr59+9KmTZvEv7ln/cCBA6JazZeZM2fSkiVLqGnTpq6hnyvebnhSoS+//JLiwsFVx559sXYTXvGH7e3GStwuNuwMAUB+yJ2CYtuSJKrcfv8GWmTyK3LFm3uX+QTBAQMGUIsWLcRQedxm4nVJcxxvDsX9+/cXI6vMmjVLPF8/+CTRIUOGiNDNI7QsXrxYTBDEgZr7sxs1aiRGcRkzZozr42VId8NjoHOvfBamRwcAyNIOZNbbTaJAwQQgnEgpmCfE4SoyT5bDAdxkvXv3ph07dlS6jof382PGjBliZBI2d+5ccfKonI1zxIgRYtSWUaNG0bx580TVm8O9ivu433rrLdffzW0mzt5vUxWqeps6dXqaLuj/9WRSadJ1xAMgj8ESLQEQl7h36nhbjf7uDAZv7nd+7bXXRCWbg+f3vvc9OvbYY42cnZKfY1gcvBmPvS1Dt2rkyJF0880308aNG+nxxx+vEry7detGixYtErN3qiGbR3zh/m9+vE1Q4QaAPIVvG8ax5vWyyUHLlJ0Yk4oROGqQT5HWJnPmzBEnJ953331ijOos4nYSruizQYMGud6HlwGPT/7QQw+JgO3EVfFf//rX9Nvf/lZMniPxv5lNwRuV7XxDfzfktTppcvhOashXv73UhR5v6nMzgS07eZBij/dnn30mqts8IUxWrVu3TrSSsC5dunjeT95WWlpapaXltNNOo8suu0xMIHTRRRfRo48+KpYZj+HN46DzyChJwyQRekOoHFkARwMAskUGuaz2e+d5rGdTKt+QL5GCN59oyCcO1q5dm7KKW0EkryEBnbepj5GmTp1Kd9xxB7366qvi6AD3gt966630+9//vuhz2L17d6WL14matsnCUILFjgrY+hoAIPvh2+TWlLjkIWzneWcq88Gbe7o5CP7973+nrOKxuyXeyfCi3qY+RuIjAzfddJOY0ZODM/eDT5gwgWrVqlX0OfBEOzxmuLzceeedlBU2HxYEALsF2THOaviGbMrjTlUugjcHyWbNmtEvfvEL+uqrr/Q9K6iET8LctWtXxYWXO5hZ6UYPPEB2IXzbe1Kl6bCc8iNS8D7uuOPEEHo8WgeP3MGjf7z//vv0z3/+s+DFJjxGt3qipRf1NvUxOjRu3LjShcdKB3M4K2B8mC/N6eLjgrYZgMrh27Tqd1LnmPgNifJ+SZ00qM6umca6LGobC06szIfIE+icfPLJYnKZtWvX0o9//GM69dRT6fjjj/e8dOjQgWzCY3BLhaZ2V29TH6MDn3zZuXPnSiOigBmyFq4BwN/33ln9TjuEJ3W0zW84lAGY749AWXm5xLk80N9tvkjv/vbt28U06jzyB/MziY7pE+04derUSUyUwyOb8M6F15CCfBvj2TC5/UanN998U1S6wR4mTJ4DAMmM8+0nfOcxfObxNfuhfk7iaDFBf7fZIn0reFSODz74QJxYeP311xs9gU5Y/Np69eolRiNZsGCBmG7ebWdi4cKF4ucBAwak8CzBJFx5Srv6BQDmHPVSw7mJ4ZSrpHk6P8WEkU2cM6Oa8DmAZER6p59//nkxecxjjz1Gw4cPp6zisbY5eC9dupRWrFhBZ511VpWJhDZs2CB+Hjt2bErPEiC+mdrQ3w1ZlcRMhl7hXA3kaQUvro7mef4Bfu/TXL/pfN/RZpKTCXR4DO9hw4aRDXbu3CnaY+RFTozDJ0aq15eVlVUJ3iUlJaKyza+Vx+Bm/HgO3Zdeeqn4N7ehOKeLN7nHO81JdPjvoj8aAPLMrVfcVrwTkfTJlIWeS16hzSTjwZtPIuRxqGvUqEE26Nq1Kx199NEVFx6mj91zzz2Vrv/5z39e6XHcOvPcc89R+/btxUmU55xzDjVo0EBcuNLPY5nz7545c2Ysz5t7vLml58orr6S8yEKFNQuvAQCSC+BRw7ecQTfpcKueSJl26E0r8JvQvgI5CN7f//73ae/evbRq1SrKOg7d7733npj0hqeH5xYb3ungYRQnT55Mb7zxBjVtiqAFAADh2VD5dgu38jpbRjFJozDitrMC+RPp2/GrX/1KtFpcfvnl9MILL9BRRx1FJtu0aVOkx/P43HxCKV8AAADiHCXFhgALwfB7ivc33yJ9q3kIvTvuuIOuvvpq0YPMvc5nnnlm0Qlk+vTpE+XP5g73eHM7D7ea5KndxAZcNXFWLrLWv46WGYDkyHWH6RVR2VaSpXWdzfjESvR35yB48xje3HLBeCrz22+/vehj+P6HDpl/KM0kGMcbAMDukU3CVr9Nr3qb+Bz97BTofL9N++xAxmeu5JE+glzkSCKQXzb0MAaVxdcEAOmSJ1wGXb8EPcGSK6VhhqKzpZ87Tn5DtzyCgW0FRAreHKLDXMAsaawIsnR4MkuvBQDsXr/kaSKctAWtdKs7KXnfYcmzyBVvsBtCYz6E7dNGfzeAGVApBS+YOMcuRgTvli1bZmqaeVsm0AEAAPNlZZIdU4Xp0Q7ymCSmhseJlfYwJu1y/ze4w8mV9pAbRlSKAexi4gmWbidbQvrveZjHoM0EjKp4A2QpdGelfQc7DwBgMx07Khyy1aBt8s4Z2AHBGxKVtTGuZUDN2msCAPMEmVY+zNTxWeoV9tPSUSxEy9tRhACdELxzwLaVhm3PFwDAJGFGNkGPcHFhqt1xtwdh4hz7GNPjDWC7LFW9sfMDYC6uemM4OnvWhXivQIWKtwXiHtUkyOFLsBfCNID935Es7eDbVN1Gbzfogoq3BbIyqklewr0NG+8sP38AABPWhRiFBtyg4g2JQrUGACC6vBQyorZ4hA2/auiOUu2Os80E/d12QvAGAACw6IiPnwIGn2AZZmQTAIgXgjeAxo21DRvtQmx//gAQTZaGFEyaWl1Hmwl4QfAG0MT20Gr78wfIkzhOiseQgnqhzQSMDd6YLj77PYPoRzQDwjUAgH8YzQQyGbxLS0vp8OHDaT+N3A4nmNRJjzixEgBALxQ19EDABquC97vvvkuXXXaZCIc87F2NGjU8LzVrYgTDMMMJfvDBB3TllVfqeLsAqkAlHMC+74OfYkbeT7CMMrJJ0DAu3w/MVgmxBu+pU6eKiuyjjz5Kf//736msrEy0jhS6AIA5ocKWkAFgK/6OxfU9K9TrHWbqeNADs1VCLMF7xYoVdPXVV4s2kSuuuILmzZsnrm/WrBm98MIL9MQTT9BFF11EtWvXpubNm9Of/vQnevHFF6P8Scg4hEAAyBIb12kY2QStJ2Bo8L7//vtFBZvD9//8z//QwIEDxfUctL/73e/SqFGj6LHHHqM33niDqlWrRrfccgudfvrpup47WLLyRw+iuRU4G0MBQBJ0fDeSPLqkaz2LkU2iSaLNBHIcvJctWyYCNQdvlbOd5Fvf+pYI5uvXr6d77rknyp8ES+HESgAArF/9hkt5sVHcbSbYOcpx8P7000+pTp061K5du//8wurVaf/+/VXue8EFF1CtWrXoqaeeivInAQAAEhGlQo2jSeHIsM3hki9pBPCwI5xg0hz/Svcc8LxkXaTgXb9+fXFRNWrUiHbv3k0HDlReeBy6+b6bN2+O8icBAAASEzRAF2vhSqPdxJYwwwFbBm5J/hwlfPsd2UTHkIJxT5qTBaUun0c+EVieDGzL5zWV4N26dWsRsg8d+s+XvWPHjhVD4Km2bdtGu3btwqgmOZsJDQDAdjZUr73a+cKMbCIrzabQEb7jllS12/Y2k9L/DdUyaKuBW70+y9XvSMG7U6dOYkSTNWvWVFzXt29fEa5vu+22ipaTgwcP0i9+8Qvxc0lJSdTnnDtJTKATFwR5AMjLCcl+A7oNQT5pxUJ1UuEbE+nEb7SPncEsV78jHRMZMGAAzZkzh55//nnq2rWruI4neeGAuGTJEmrTpg2dfPLJ9OGHH9KOHTvEiZg///nPdT333OCjBzwxka1wYiUAZAmCczwtJsXwfUwOYhi7W6/RAzuJzwa/57ZX+rVVvIcNG0YTJ06kVq1aVVx3/PHHi/G6udebw/by5cvpiy++EKH7hhtuoNGjR+t43gAAAJDDo4ymtZwk0Wbid+fEZBygRwdsfcpi5TtSxfuoo44SwdttBJPvfOc7YkKdLVu2UJMmTUR1/IQTTojy5wAAADJTNS/W1uDnPs6ji7aNrBE0REepevOyiesILKrd8Rmdscp35CnjvfDslWPGjKGbbrpJzGqJ0J0/eai8AADE2c6SZltL3JVldejAsI/1C8E4XVEr1qMzVPmOLXiDedJYgQepLqBvEgDyxGudF8e6MGhgibuyGCV0m1T1RJuJf2FG2InyeOfY4KaMFY7gDVWgUg0AkAw/08oHCeJuw8JGDTwmhW4ToZpemM6gO7rIUINquHYbG9yE4QrjndcUrGNjjyAAgM38BOug/d6m0hm6+XeZtlMBlamBV3fPtxfn33L728V+B1ODuc6dRFS8IRaomgMA6JV0O15cE+noCDFhf4dthSXTRnBJO3SrnBXrMKOmFKqcx9VbjuANscH43QAAYEqI1N0Sgtkq0wvdo5VArLaWRP0dbrfrfg0I3gAAABmrerv1ead9UpnuQ/YmTG0fZ3932q8tqrjbgEY7ppoP+zucveDFquBRIXiDK7SKAADYSefRRl3h1vYQ6YRqt7ukT1ocrYRmHb/L7We3f0eB4A3aV9oc2tFmAgAQD9uGXo17FJOgoV5XaMZoJpU5Rw6x0egiz13HTgWCtwW6d+9OnTt3pgceeCDtpwIAAGBM6A76e20JyzYdIVCr3LYGbj8KvbZPyw6SX3Z8AnPuzTffpMaNG6f9NAAAwBBhhxcMO+12lKH7sjJed9JsWm5ZDtxOcidDvj9Bq+CoeOdM3Icow/aG23boFADANlGDc9AqbJJV2yT/Freq4KTKr5lwwm6SvE7EHH7OSb5/ByreoB36uwEAzJgMTVdA5PAdJmQlUbUN89x42Zi8rbKh2p2H9hIvzte8f+8e8gsVb/A97TAAAJgjjSOFfivLafQo+/2bUXZG4h7NxKbe7ryG7qgQvAEAAEBbFTbuEyrdJPm34j5B04ZqN4SH4A3aYBhBAACzRZ1Ip9C43nw9X/g+NoRH06aPt6naHWamSPgagjd4QrsJAIDd4mgZlOFbBkWTAnec7SZxBvU0jhKElbcTKnVD8AYAAMhgn3ecJw/KkG1SYEziOcTZZmLCMvQL1e7wELwTUFZWRpMmTaLzzjuPWrRoQdWqVaOLLrqIsgRtJgAA+WNClTsKv1XsJKrdkA8I3gnYvn073XrrrfT222/TGWeckcSftArG8AYAAB0K9aBHrV7HUe026YiBH+jtjg7BOwEtW7akf/3rX7Rt2zZ68sknyaag66fPG8MOAgDYVbxAn+5/wnSxanZc1W7bQjfogeCdgDp16lDr1q0py0yeiAAAIK/ciid56M8N2r7B4dp5kdfHUe1G6M4vzFwJvvCK223Fg2o3AACYJOhMls5tm9yuxRm6baxy4yhJTive+/bto/nz59Ptt99OQ4cOpXbt2omTFfnCJzD6sWfPHnHfkpISatiwITVp0oS6d+9OU6ZMoYMHD8b+GrJWzUa1GwAATBP2pEUO2urFhOdkijwcLYmbdRXvlStX0uDBg0M/fvPmzdS3b1/atGmT+Hf9+vXpwIEDtGrVKnGZOXMmLVmyhJo2rdozV15eLu7rR/Xq1al27dqU5ao3qt0AAOb0eb++fmfaT8Paqnfc1MBtY7UbclzxZhyK+/fvT+PHj6dZs2aJIfr8OHToEA0ZMkSEbj7hcfHixbR3715RRZ89ezY1atSIVq9eTWPGjPEM7fXq1fN1Of300ylLZFVbhm35f1S7AQDsZFIwzTK1nxuhG6yrePfu3Zt27NhR6bobb7zR12NnzJhBa9asET/PnTuXevToUVGdHjFiBB05coRGjRpF8+bNE1VvDveq5s2b07Rp03z9LbeKue04ZHPPm85KN4YSBABI/ogltwzY3vbgd2jBNNsjcBIlWB+8a9SoEfqxHLxZv379KkK3auTIkXTzzTfTxo0b6fHHH68SvLkfPGsT3wSFCjcAgF3tJrJoAslC6IbMtJqEwe0ky5YtEz8PGjTI9T58gubAgQPFz4sWLUr0+QEAAIB+aVT2sxa6MXFOjiveYa1bt060krAuXbp43k/eVlpaKlpamjVrpuXvT506lf7973+LPnP23nvviZFZWJ8+fcQlSTgRBwAAsi7pkyyzFrhBv9wEb541Uio0mY16Gz9GV/CePHmyODlT4pM4+cImTpxYMHjv3r27yoQ8fAEAAAhbVJF93nkIiUn0eiN0gx+5aTXhsbslHkLQi3qb+pioeCQVHo7Q7VJs/PG2bduKscbl5c4776QswImVAADpzWCZF0nuWORhJwaiyU3F22Zbtmyhxo0bV/wb1W4AAAAzqt5ZHx0Gw07qlZvgzWN0qydaelFvUx+TJg7davAGAAAAc3q9s17pxoyV+uSm1aRVq1YVP2/dutXzfupt6mMAAACy2MaXp4qm7up0XnrkQZ/cBO9OnTqJiXLY2rVrPe8nb+PZMHWdWBlV9+7dqXPnzvTAAw+k/VQAAMBizj7vPFUyZUDWEb75d+QhdOdppywpuQnefNJkr169xM8LFixwvQ+f6Lhw4ULx84ABA8gUb775Jn3wwQd05ZVXav29OLkRACA/MAGanvCdt9FL8rRzloTcBG82btw48f+lS5fSihUrqtw+Z84c2rBhg/h57NixiT+/PEHoBwAA28J3nkI3qt3xsDJ479y5k7Zv315xkRPj8ImR6vVlZWVVgndJSYmobA8bNoyWLFkirufHc+i+9NJLK2a2dE4Xnya0mgAAQJzyFrLU8F0ogMvb1fvlIXRLqHbrV62cU6hl2rdvX2kyGi8ctKdPn15lPO1+/fqJ/8sWFA7e+/fvF//u2rWrCORNmwY7OSUOPHEOj9u9a9eu2EY18TvRgm6oeAMAxMdr3T7/o+3UsE7VAc3y0K9cbKfDGTLzGLQlTBEfzP69e2jSkNN85bXcDCeohnaerp1nknzqqado48aNVKtWLTrllFPowgsvpKuuuopq166d9tPMNIRuAID417Ne4ZtPsHQL33kfatCt8p3X0A3xsfKbJ6vVYfH43Lfeequ4AAAA5OkES656Q2V5DNhuUOmOn5U93qAPqs8AAACASncyELwtgJMrAQAgTtzfjOAFOJkyfla2muQNj+ONKeMBAAAgDmgxSQ4q3pAotLYAAJg3gyUAJAPBGwAAIEeFDsxgCSq0GCULwRtQhQYAAAEhLJ/Q250cBG8L4ORKAACIG8IXQPxwcqUFsnJyJfq7AQDMgYl0AEc4koeKNwAAQM6gzxskHOlIFoI3AAAAVEAVNB/yOIRg2YFDqY/og+ANAtpAAACyJ+i6PW9BLK9s3bkq+9/gXOzi9Th5tCfIY3VDj7clJ1fWqFGDrrzySnGxEYI9AIBdwaxFozppPw2IMXTbtJNVpoTmYuZ/tN01QKuPLfR7+PEN68QXjxG8LZCVkysBAMD8Eyw5kM1csC615wTxsS10lykB2u95CXw/Ds9BH+f2d+MI4AjeAAAAOeQMKJAPtoXuQSGCc5STh+Vj4/puoMcbYm8HQZsJAIB9bO0DBvvfz7IIoTuO56ETgjcAAABYWRWFbLWYlDlOgkyT/Pu6wzeCNwAAQIYVO+pYKFjYVCWFwkwM3WUuI4pw4E07dMcZvhG8AQAAcqpQwJFBDeEbdCrzGLrPpMCt0v2ccHKlBWweThD93QAA9pIjnGB4QXuZMlFOmBFKTKJrpBMEbwskOZwgB+XX1+9M5G8BAIC5wwpKCN/2MuVohSl92yaMdIJWEwAAgBzzE4ZMqJhCOKa8d7aGbqeoM1wieENs0GYCAJCt9bEpFVSw472SITUroXuQ0oceNnwjeAMAAEDRIFHoZEu+Tr1kgfM1JfW6dPyttIcPNGlYwDio4TtoAEePN1SBPm8AgHzxO4ul2u/tdhtzu71FozpkMq+gqwZX5+vS/Zqcz8HrOXn9Xef9kw7dbqOUZNkgpe9774HDvh9Xrby8vDzG5wUR7N69m5o0aUK7du1K7ORKKeoJlmgzAQAwS7H1ugzeUUdtcOLAalr4LrTjEOQ16XhdQarTXssyzZFLbB+tRIe9e3bTgK7tfeU1BG+DIXgDAEDS4Vt38DYpfKuBW0dQjRLAwz4X57JMs60ky+0kQSB4Z0SawTtq1RsVbwAA86QVvJ0hNckQnkQLhvO1OV+fn1aWKH8vzbaSvIfuoMEbPd4AAADga0zvKNz6peMO30lVg52/363PXedzSHOIQFS5o0GriQUV75NOOimVmSvDVrxR7QYAMFeaVW9nOI0zeKc9skcWIXS7Q8U7Y5KcuVKF0U0AACBOcYwSoruPG76G0K0HWk1AK1S7AQCyMbRg3FVvr6H6ogRwVLnjgdCtD4I3FISqNwAAxM1tDPBiJyd63Y4qdzxwEqUeCN6gDardAADZKajEdZJlIYUm4Sk0mY3zdtADlW79ELwBAAAg1EyWcSkWohGy44fQHY/qMf1eyBBUsgEAsgXr9eKh0+2SN2gv0Q8Vb9ByaBIrcQCA7Emj3SQNbqHaGTr5CIC8X9aXSR53MpKS7U8OJAKhGwAgewWVtNtNTGunkPfJegBHi0m80GoCkQI2QjcAQLZltfoZNmDy/eVjstSGor4OtJjEBzNXWjBz5a5du1KZQAcAALLPzwgnSc1mmRQOmDrDpe3LB4E7GsxcCQAAANZNqpOEuKrTtvXDq8sBFe7koNUEAAAgx4K0DNreUhFXZVdtPTGV2wgtatsMJAPB2wLdu3enzp070wMPPJD2UwEAgJyyPaDF3U7h7Ps2OWgjcKfHnmMiOfbmm2+ixxsAAMDwnQfnyCdJt544A7/tO0txHsXxc25DHBC8AQAAcs7vNPI29jKzpCvQsic+iWWFXu1wbVPyvkkHcLSaAAAAgC9eFVQTWiuitpgUC218e5BgF2ffN3q1q74vtgxvjOANAAAAoYOLDJYmhm9doTvsfXWHb69+7bzqqSFsJx3Y7TpWBAAAAKlym9HSxFkuw4RuvyEsSGuOrmWDsbbjC8xB3s+oMIGOwTCBDgAAJC1sAJHh0oT+7zgq3VGXVZRx0BG6k6lQh/3sB5lAB60mAAAAEDnYBG15iGu69aRCt/wdcY+DjtD9tST6uJNoO0l/tzQH3n77bXriiSdoyZIltHHjRqpTpw516tSJrrvuOvrBD36Q9tMDAABIlBom4xj9I4nQ7fx9xaql6lCDfl5r3gN3z5ROlnT7uzrbUFDxTsDdd99NM2bMoDPPPJPuueceuvnmm+nAgQN0wQUX0IQJE5J4CgAAAIlUvYtVdZ2BUucJiH5+hwmjX/h9rXkM3T0NHKFE5/NBj3cCli1bRt26daO6detWXHf48GE6++yzadWqVfTpp59Ss2bNqjwOPd4AAJCmMJU+9URCZ2W3UBU3Sg90sd9tWo+wsx/eWfHPa7W7p2GB2+97HKTHG8E7RVz9vuGGG2j58uX07W9/u8rtCN4AAGBzAHcL3oXCZNjwbULoDhO++bV6Vb9tC91RRgbpaUHgLiRIXkOPd4q2bdsm/n/00Uen+TQAAABiCVVhTyYME75NCap+l5U6zKDa/63+2wZhhmPMM+t6vPft20fz58+n22+/nYYOHUrt2rWjatWqicukSZN8/Y49e/aI+5aUlFDDhg3FXkr37t1pypQpdPDgQUrC1q1badq0aXTWWWdRx44dE/mbAAAASXELj8UCZdB+7yCjoiQZCv3+LecEODZNiGNiL7YNrKt4r1y5kgYPHhz68Zs3b6a+ffvSpk2bxL/r168vTnTkXmu+zJw5U4w+0rRp1Q9TeXm5uK8f1atXp9q1a3vuPPCJlfy7HnnkkdCvBQAAwOSqd5gQ6XfCGZPaS9KelCVJCNs5q3gzDsX9+/en8ePH06xZs6hFixa+Hnfo0CEaMmSICN0tW7akxYsX0969e0UQnj17NjVq1IhWr15NY8aM8Qzt9erV83U5/fTTXX8HV9S5Us9DDHLIP/XUUyMtCwAAgKyFrmKjo9hy8mHWQmrWXk8arKt49+7dm3bs2FHpuhtvvNHXY3lIvzVr1oif586dSz169KioTo8YMYKOHDlCo0aNonnz5omqN4d7VfPmzUV7iB9uFfOvvvqKhg8fTosWLRK/hwM4AAAAFO/3dgZxP6HbhKCYlcq3CcsyC6yreNeoUSP0Yzl4s379+lWEbtXIkSPp+OOPFz8//vjjVW7nfvCLLrrI1+X888+v9FgePpBD/bPPPksPPvggjRs3LvTrAH+4lYd7+f22BwGWZ9LwGcXytPHzmWTV29nHLXugbQnd6nPp1qY+LfrjfXTQwm2SScvS9nVoJoYTbN++vWgDmThxoucJltxOwq0kXNXmCW24TcXNFVdcQQ899JBoX/nkk0+0PD/+m2PHjhWtJffeey9dc801vh6H4QSjwfLTC8tTPyxTLE9bP582VHBNC4t+vu8mLlfTlqOJ61AMJ+hi3bp1IgCzLl26eC48eVtpaaloaXGb2CYoDvkcurnKzu0qPH28qmfPntShQ4fIfwcAACDJMGZiULSZSW0pJgdum1nX4x11zGzWunVrz/upt/FjdATvt956S/yfJ8rhixP3eyN4AwAA6GFzaEx7p8bmZWeD3ARvHrtb4iEEvai3qY+J4qWXXgr1ONkFxGN+82EMqU6dOuIChcllpi47CA/LUz8sUyxP2z+fXY7+z3lXKzZUDopndWha5bqk7N4d/nwwU77vvGyTXn78ntm0zdxtyHZe/n0/3du5Cd42ksG/c+fOaT8Vq7Vt2zbtp5ApWJ5YpqbDZxTL02T4fGZ3mXJu457zQnITvPnESvVESy/qbepj0tCqVStav3491apVS8zMKaHiDQAAAGAGrnRz6ObcVkxugre6MLh1w2viGr7N7TFp4PHF0fsNAAAAYLZilW5rx/EOq1OnTiLIsrVr13reT97GwwnqOLESAAAAACBXwZtPmuzVq5f4ecGCBZ6HChYuXCh+HjBgQKLPDwAAAACyLTfBm8nZIpcuXUorVqyocvucOXNow4YN4mee8AYAAAAAINfBe+fOnbR9+/aKi5wYh0+MVK8vKyurErxLSkpEZXvYsGG0ZMkScT0/nkP3pZdeKv49aNAg6t+/fwqvDAr54osvxJjnY8aMESO9NGjQQJxo2qZNG/rBD35ATz/9dNEFyCc/8Oym/Dlo2LCh6Mnq3r07TZkyhQ4ePIg3gIh+/etfi5N55QXLM/zwUnfddZeYIOvoo4+u+Kz269dPfAb//e9/4zPq0+LFi2n48OHUrl07qlu3LtWrV0+c/zJ69Gh6+eWX8Rn9X7wNnD9/Pt1+++00dOhQsbzk99hrVmenTz/9lK6//no6+eSTxXLmlsvevXvTH/7wB19DpfGAAD/96U/p+OOPF+8Vf/a/973v0dy5cylPy5PPF3vwwQfpv/7rv+iEE04Qy5IvvFwuvPBCevHFFxN5P7L2+XS6/PLLK34Hz2JuxfIst1C7du146RS9jBs3rspjN27cWN6+ffuK+9SvX7+8bt26Ff/u2rVr+Y4dO1J5XVBYzZo1K72//L41aNCg0nWDBg0q37t3r+vjN23aVOW9r1OnDt57xd///vdK34dCqwgsT28vvvhi+bHHHluxDGvXrl1+1FFHVVquq1evxjIt4siRI+U//elPKy23evXqiYt63bXXXovPaHl5+dKlSz23hxMnTiy2uMtXrVpV/o1vfKPiMQ0bNqy03v3e975XfuDAAc/H/+1vfxPrVXn/xo0bl1evXr3i3xdffLF4T7O+PP/5z3+WV6tWrdL9ebk4P7c//vGPyw8dOhTb+5G1z6fbelZdzpwNCzFleeYueLPdu3eXT5gwobxLly4iuDVq1Ki8W7du5ZMnT7bqQ5w3/J6eeeaZ5Q8++GD5+vXrK+1M/eQnP6l438eMGVPlsV999VV5SUmJuL1ly5blixcvFtcfPny4fPbs2eIzwLcNHjy4PK94WfTs2VMshx49ehQM3lie3l577bWKDezQoUPL33zzzYqwwTuFK1euLL/55pvLN2zYgGVaxGOPPVbxOfzhD39Y/uGHH1baSTz//PMrbn/qqadyvzw52DRt2rS8f//+5ePHjy+fNWtWeYsWLXwFm3//+98V9/3mN78pPreMt4lTp04tr1WrlrjtZz/7mevj+fMsCyG9evUq/8c//iGu37Nnj9jeyvfprrvuKs/68uRtEt+HHzdjxozyrVu3Vqxj33///Uqf21/96lexvB9Z+3w68bq0Y8eOYjmcccYZRYO3ScvTyuAN+cR7t4WolTGuOKj+8Ic/VNz2+uuvV3nsn/70p4rbX3jhhfI8+u1vfyte/+jRo8VKsFDwxvL03hh06NBBLLerrroq0PLHMq2qb9++YlmecMIJYmfP6eDBgxXLe+TIkblfnm7VU1moKhZsOADKIwrOnUJ2xx13iNtr1KhREapVXPDg2znc7Ny5s8rtl112WUUV3JajymGXJ4e8t956y/N23hEfOHBgRdX1yy+/1P5+ZO3z6XTNNdeIx3ERg4usxYK3ScsTwRsygyuJXtWv3r17i+v79evnuSI8/vjjxX3Gjh1bnjeyWsWH4T777LOiwRvL093vfve7ivDhtjEtBMu0qpNPPlksz2HDhnkuNz6qwPc577zzsDxd+A02xx13XEU7iBuuXHNI5PtwBVtVVlZWcZTn1ltvLVgF5gsfybBV2KDo9Je//KViebz99tta34+sL8/ly5eLFqaTTjpJrGf9BG+TlqeVJ1cCuOETeaTDhw9XOqFj2bJlFSfOuuETMwYOHCh+XrRoUe4WMJ9YvHfvXvrNb34jToYqBMvT2+OPPy7+zydUqZ/HYrBM3ckJxN599106dOhQldu/+uoreuedd8TPZ5xxBpZnSP/4xz/on//8Z8F1JJ+Mziehua0jX3vtNfryyy8LPp5PfOP5NNwen0de2ysd70eWHThwgH784x+LEyEfeeQRX+tZ05YngjdkxksvvVTxM49aIq1bt65i5JsuXbp4Pl7eVlpaSjt27KC8+P3vfy9G+DnnnHN8DaOJ5em9QVi1apX4uVu3bmJFf9lll1Hbtm2pdu3adOyxx9KQIUPob3/7G5apTz/72c/E/z/++GMxEgT/X92Y8kgnPARsx44d6dprr8VnNCR1Ujk/68gPPvgg0uPff/99yju5veJ1w0knnaT1/ciy2267TWyDfvKTn9B3vvMdX48xbXkieEMm8NBsd955p/iZ91p5qCBp27ZtFT+3bt3a83eot6mPyTIe8mr8+PFiWKWHH37Y12OwPN1t2rSpYkhKDoO8Euedms8++0wMfcn//+tf/0rnnXeeOMKgDl2FZeqOd1TuvfdeEU6efPJJOvHEE8VkaHz55je/KcILh/OVK1dS48aNsTxDCvr546Ey1eF65eObNm0q1iXFHp+X9auXjRs30u9+9zvx84gRIyp9dnW8H1m1evVquvvuu0UR45577vH9ONOWJ4I3WI+r2T/60Y/ok08+EYedpk6dWmXsbok32F7U29THZBmPt7tr1y4xhqo8rF8Mlqf3/AISj1Nbq1YtMT8Ar8D5ts2bN4sWFMZjxnKgxDIt7pprrqGnnnqKjjnmGPFvbmmQbQ28o8PLlz/D+IyGF/U7LX8u9Fj19rysX93wZ5fXA9xe1rx5czFvghPWsVVxqxm3mPD/77//fjrqqKPIL9OWJ4I3WO/qq68WlUT2wAMP0Kmnnpr2U7LCE088IdoevvWtb9F1112X9tOxnmxnkj8/+uij9MMf/lAEcHbcccfR7Nmz6bTTThP/vuOOO1z7luE/OJxwRZCPEvDy497Lzz//XFz4Z55I649//COdeeaZ9N5772HRgdH4+z5q1Ch66623xHph5syZ1KpVq7SflhV+/etfi/M5eF3ALWY2Q/AGq/33f/93RYWbK4i8R+zUqFGjShtyL+pt6mOyiGfv4kpijRo1RDtEzZo1fT8Wy7P4cuGWCJ5N1al69eriMytnYuUNMJapN26D+stf/iJax1599VU699xzRZWQL/zzK6+8IvpjeabiK6+8Ep/RkKJ+p+XPhR6r3p719asbPoGSZ1p95plnxPr2T3/6Ew0YMMD1vljHVsY91//3//5fcQIkzwYalGnLE8EbrHXDDTeIqd7Z5MmTRZB0o1YUuKfZi3pb1qsQN954owh+fPIf98ry4Xr1InuVmfM6LE93au8gL1MvXKWVuP0Ey9QdH+rlUQsYh2q30Qu4n/jnP/95xcga3EeP5Rlc0O809yRzCHI+nluqZBtQocdnff3qFrrHjBkjdiK52MFHG/loWFzvR9ZceeWVYvtz8803i/MInNsreeSQz5uR1/GIR6YuTwRvsLYSJk+u4JMtrr/+es/78hBWXGl0nt3sJG9r0aIFNWvWjLJ+cg976KGHxJ698yJPVGXyOt7RYVie7vgzU+jEHUk9qZKHscQydffhhx9WbFB51BIvfHTB+bnGZzQYdaQHP+tIdecxzONPOeUUylulm9vMZOjm9qk434+s2fi/3+ubbrrJdXvFLTuMR5KS13HbqanLE8EbrMOH6rnCLUM3h/BC+ISJXr16iZ8XLFjgGYYWLlwofvY6/AdYnsXIzw4Pd+VFHarq+OOPx2fUg9xZVo8MeLVNOQ8P4zsfDLfrcA99oXUkj/PP7T5u68izzz67YjQTr8fzeyi/F3lZx3Lo5p7uP//5zxWhe+TIkbG/H2D48ox1eh4Aza6//vqK2b4mT54ceDruatWqlb/xxhtVbv/zn/+c+ynjVX6njMfyrOyVV16pWG5PP/10leV2+PDh8lNPPVXc3rp1a/FvLFN3+/btq5gN8fTTT3edMp6noO7Zs6e4T9OmTStNSY3PaLgp4+vXry9mmXS66667fE0Z37JlSzFlutPPfvYzcXujRo2smTI+yvLkz+KIESPEfWvWrFk+e/bsQH8n6vuRp5lAxwWYMt6E5YngDdYYP358Raj5zW9+E+ixvNEuKSmpCDwvvPCCuJ6DD0/d27hxY3HboEGDYnr22QreWJ7efvjDH4rl9o1vfKP8ySefrAiMmzdvLh8+fHjFcp0+fTqWaRFXXXVVxfIaOHBg+XvvvSe+s3x59913ywcMGFBxu3Oq8rx+RjnUfv755xWXtm3bitfJ60/1ep4iW8VhuUWLFuK+nTt3Ll+1apW4/sCBA+UPPvhgee3atcVtHKDdbNiwobxBgwbiPr179y7/8MMPK6aT5/eGd9L5Ng44WV+eHLpHjhxZEbr58xZU1Pcja5/PqMHbpOWJ4A1W4NAiN7DVq1cvP/bYYwte7rnnniq/g/dy27dvX/F7eM+3bt26Ff/u2rWr1ZWYJIM3w/J0x0GjT58+FcuvTp06ohor/12ouoNlWrXqzYFbXXa8PPmiXnfhhRdWqnbneXnKCmKxC4cVJw4jvMMo78PV6Vq1alX8m3d09u/f7/m3//a3v4llLO/fpEkTUUGU/7744ovLjxw5Up715fnyyy9XXM/Lr9j2yqsaHvX9yNrnM0rwNml5IniDFXgD6ufLWizY7N69u3zChAnlXbp0EdUZ/uJ169ZNtK3wni/4D95Ynt64qvr73/9eBPBmzZqJlTtXXbkKtmzZMizTADiozZkzp/z8888vb9OmjahMcfDmStmwYcPK//rXv2J5agw2paWl5ddee235iSeeKHZSjjrqqPKzzz5bfJ7V1igvH3/8cfmll14qdnj4fWrevHn5ueeeK47+5GV5Ll26NND2atq0aZ5/P+r7YZp2KQZvU5ZnNf5PvF3kAAAAAACAUU0AAAAAABKA4A0AAAAAkAAEbwAAAACABCB4AwAAAAAkAMEbAAAAACABCN4AAAAAAAlA8AYAAAAASACCNwAAAABAAhC8AQAAAAASgOANAAAAAJAABG8AAAAAgAQgeAMAgBa/+tWvqFq1anTXXXcZtUQHDhwonteLL76Y9lMBgJyrVl5eXp72kwAAALv961//opNOOokaNmxIGzdupAYNGpAp3njjDerRowd17dqVVq1aRdWro+YEAOnA2gcAACK7+eab6csvv6QbbrjBqNDNvv3tb9P3vvc9Wr16NT3xxBNpPx0AyDFUvAEAIJKtW7dSu3btqEaNGvTJJ59Qs2bNjFuizzzzDF1wwQV0yimn0Nq1a9N+OgCQU6h4AwBAJL///e/p8OHDNHjwYCNDN5PP7f3336dly5al/XQAIKcQvAEAcuq///u/xUmH5513HvHpPn/605+of//+9I1vfENc///+3/8r+jv4cY8++qj4edSoUZ7327lzp7jf8OHDqaSkRITgunXriko5P477sL20b99ePJ/p06d73ueiiy4S9+H/u6lduzYNGzZM/PzII48UfV0AAHFA8AYAyCnueWYdO3YUI3+MHj2aXn75ZRGmOcSefvrpRX8Ht23wiZWsd+/enve777776JJLLqE5c+bQunXrKq7/5z//SbNmzaKePXvS/fffT3Hq06eP+P/ChQtj/TsAAF4QvAEAcuqdd94R/3/sscdo5cqVohK8a9cu2rFjB33++ecVQbWQV155Rfy/bdu21KJFC8/7tWrViiZOnChGFdm3b5/4G3wy5oYNG+jqq68W97nuuusqdgbicNZZZ4n/f/rpp/T3v/89tr8DAOClpuctAACQWZs3bxbhl3F/Nvc9n3rqqRW3c7uJHytWrBD/P+200wre77LLLqtyHVfVjz/+ePrtb39Lhw4dogceeEBc/vCHP1AcTjzxRDHcYVlZGS1fvpy++c1vxvJ3AAC8oOINAJBDamX5zjvvrBS6g9i2bZv4/9FHHx3p+fyf//N/xP9fe+01ipPcoZDPGwAgSah4AwDkOHg3adKEfvazn4X+PdySwvyMZsJtJQ8++CAtXbqU1q9fT3v27KEjR45Uuo/sF48LP0+u9svnDQCQJARvAIAcB2+uNPOIH2Ht379f/L9OnToF7/f000/ThRdeSAcOHKi4rnHjxmJkE245OXjwoBj5ZO/evRSnevXqVXreAABJQqsJAECOT6zs16+fltYNDs1evvjiCzHMH4fu7373u/TSSy+JEyz5RE4+0bG0tFSMdpIE2dfut4cdAEAnVLwBAHKGg/CWLVvEz36GDCxE9nbLQOtm3rx5tHv3bmratCk9//zzVL9+/Sr34fDtpWbNmkWr1Bzi/ZDPM2pPOgBAGKh4AwDktM2EW0y6dOkS6Xd17ty5on/biwz5J598smvoZi+88ILn4zmwq7/HifvEeZjCYrinfPv27eLnTp06Fb0/AIBuCN4AADkN3qecckqk/m4mx/p+9913K/Vvq/gETvbhhx+6Vq257YVnzfQihyrkPnGe3MdpxowZvk7K5HDOIZ0r6L169Sp6fwAA3RC8AQByGryjtpkwDrAcZPnkSNk37jRgwACqXr26aPPg2TG3bt0qrufH/OUvfxG3N2rUyPNv8EmZjGe85PHAuVWGcfvKvffeS5dffrmvUVXkmOP8unk8bwCApCF4AwDkjM7gzSOTyDG4n3vuOc+Ja8aPHy9+fuqpp6hNmzZ01FFHifA7YsQI8f9C08X379+ffvSjH4mfeXKd5s2bi/YTvvBslz/96U9pyJAhRZ+rfH6jRo0K9VoBAKJC8AYAyBEeTYRbPnQFb8bBl3G7iFsrCPv1r39Njz/+OJ155pliSL+vvvqKTjjhBPrlL38pdgR4SvlCpk+fTvfddx9961vfEo/nlhGutnPFvFBol7gHnWer5MeOHTs25CsFAIimWrnXWhIAAMAHDsEnnXSSmBTn5Zdfruj7Nsltt91GEydOpIsvvpgee+yxtJ8OAOQUgjcAAEQ2a9Ys0cIxaNAgMXygSXhSnvbt24tRTf7xj39Qu3bt0n5KAJBTaDUBAIDIRo4cKdpI5s+fTytXrjRqiU6dOlUMI/iLX/wCoRsAUoUJdAAAIDKe9v3hhx+mZ555hj7//HOjlmiDBg1o0qRJdM0116T9VAAg59BqAgAAAACQALSaAAAAAAAkAMEbAAAAACABCN4AAAAAAAlA8AYAAAAASACCNwAAAABAAhC8AQAAAAASgOANAAAAAJAABG8AAAAAgAQgeAMAAAAAJADBGwAAAACA4vf/AQm4kY/tN67RAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "axes = visualizer.plot_kde(hdi_data)" ] }, { "cell_type": "markdown", "id": "a09acd68-cdc1-405a-b14d-1a4e19bb53dd", "metadata": {}, "source": [ "The colors correspond to the 68.3, 95.4, and 99.7% (1, 2, and 3$\\sigma$) HDIs. You can directly access to this by e.g.,`hdi_data['log10_Sigma_d']['hdi_map']` that have integer values of 0,1,2, and 3 for >=3, 3, 2, and <=1 $\\sigma$ regions, respectively. Here, `hdi_data['log10_Sigma_d']['X_fine']` and `[Y_fine]` have the grid values. It is also possible to directly edit the plot using the returned variable (`axes`)." ] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.13.12" } }, "nbformat": 4, "nbformat_minor": 5 }