{ "cells": [ { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from poke.writing import read_serial_to_rayfront\n", "from poke.poke_core import Rayfront" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Intro to Poke: The Rayfront\n", "Poke's one and only interface is through the `Rayfront` class, and is in essence the most \"supported\" way to use Poke. \n", "\n", "Poke’s sole interface utilizes the Rayfront object, which contains the totality of the ray information needed for physical optics calculations. Rayfront is a portmanteau of “Ray” and “Wavefront”, encapsulating Poke’s mission to link ray tracing with wave propagation. The Rayfront is first established by initializing it with parameters of the optical system, such as number of rays, wavelength, aperture, and field of view (shown below)." ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## Displaying Footprint Diagrams and Ray OPDs\n", "Footprint diagrams and OPD maps are (one of the many) important tools in any respectable ray tracer's toolkit. The footprint diagram is a simple map of the rays on a given surface, and the OPD map is the optical path experienced by each ray traced. Here we show how these data are accessible with Poke.\n", "\n", "For now, we use an already traced Rayfront (for the EELT!) but this demo will later be updated with the ray tracing included. We begin by using the " ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "[{'surf': 1, 'coating': (1.2+7.115j), 'mode': 'reflect'},\n", " {'surf': 3, 'coating': (1.2+7.115j), 'mode': 'reflect'},\n", " {'surf': 5, 'coating': (1.2+7.115j), 'mode': 'reflect'},\n", " {'surf': 8, 'coating': (1.2+7.115j), 'mode': 'reflect'},\n", " {'surf': 12, 'coating': (1.2+7.115j), 'mode': 'reflect'}]" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Load a rayfront\n", "pth_to_rf = '/Users/jashcraft/Desktop/poke/tests/ELT_rayfront_aspolarized_64rays_0.6um.msgpack'\n", "rf = read_serial_to_rayfront(pth_to_rf)\n", "display(rf.surfaces)" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Great! The Rayfront was successfully loaded. Now we use the data attributes to generate a footprint diagram at M5 (the 4th surface in the surface list above)" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [], "source": [ "x_m4,y_m4 = rf.xData[0,4],rf.yData[0,4]" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Now we just do a scatterplot" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcoAAAG8CAYAAABe5+UAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAAsTAAALEwEAmpwYAABM8UlEQVR4nO29e5xddXX3/15zn8xkJgkJuWlIIhhAbioiPCpFBUVR+yC29kIrFcvj81iUVhFQininiA8qSltarFb9VZ9qUSotFGwRtRDKVQORFHIBQi6ThMwtczsz6/fH3hNOZs45c87MuXz3Pp/36zWvOee7917fz9r7zPnOd6/vWtvcHSGEEELkpqHWAoQQQoiQ0UAphBBCFEADpRBCCFEADZRCCCFEATRQCiGEEAXQQCmEEEIUQANlHWNmW83MS/hZXSVdl5jZJdXoa0q/rzez/WZ2RZnsleyHmV2d59yPmdluM/uFmX3GzI4oYON8M+s3s/Pn7ISoOGa2Osf1flcRx7Wa2XNTjrugwP6dZrZtpv3EdDRQ1jHuvtrdLeu95foB1lRZ2iXxT7U5HOgGVpXJ3iWU6Ie7Xx2f82/GTX8Uv+8EXg18HTgPeNLMPmFmuf6GV8b7r5ylblFF3H1rfI1fn9X8sSIO/SNg+eTr+O/1GwX2/xzl+2zXFRoohYhx9+8BK4AP1FrLVNx9xN23uPvNwMuB24CrgS/n2PcvgBfFv0Xy2Aa83Mzekm8HM2sELo33nREzOxW4CHiwLArrDA2Uohi2AfMp8o8yybj7DnefqLWOQrj7MHA+8AzwJ2Z2do59tlddmCgX18W/C80qfxfYB9w9kzEzawb+Bvi/wIa5iqtHNFCKvEzGToAj3H3A43qHZtYUx98eMbMhM+s1s383s7fmsVPU/pPxOeAI4IgpsZczzOyCKW0nmtlXzexZMxsxs01m9qdmZlk2f2fKMUeZ2SfNbIuZZeK2q83s7qx97s46/tQpx7/ezD5sZpvjPn9tZr9bih/luDbuPgDcGL+9LKvv7LjzN6boOizW/lMz2xnrf9LMrjWzzlz9mNlCM/tavP+wmW2Mr+UfTfFrtZldPqWt08y+YmbbzWw8OzYWn8evx+fvgEWx4TvN7A05NPw62ycze7VFsdoDsf73Z2n9dmzreTO72czmlXJe43P0f+PPx4iZ7TKz75rZsVP2+6ssTVvNbFm83/NmNmhm/2pma0vpO4ubgd3Aa83sdTk0GnA58Pki7V0BtAKfmqUe4e76qfMfwKOPwrT21fG21VltjcCPgVHgj4liYcuAr8b7fniKjZL2j4/ZCmwtoPcb8bGPARcDXUTxxRvi9i8XOOZfgQ8DS4CjgV3A1VPOxd05jr863nYH0RfPYqLY7X8AE8CrSvVjhmsyqfeCAvu8Mt5nDOjMcd2+MWX/d8VaL4v1zwfeHOu8F2icsn878AgwAlwAzCOKe34T2JhPH9Esx4EfAe8BFgCnAkOT+xPNbB4GXgW0Ef1TcQMwDry9wGfxF8APgJfEPtwct/828PfAafHn4X/H7TeUcM6XA08BzwJvIBpcjgHuAQaA0/Nc4+fiz8Xr43P6FqAP+FWJ1/wM4r9DotmkA7fn2O/c+Pw3zPQ5IfqMDwFnFPu50k+O81hrAfqp/U/8h1PoZ3XWvh+K2z6Xw85/xl90J852/7h9K8UNlN+c0m7Af8XbTslzzN9Mab8EeNeUc3F3jj6vjrfdOqX9lLj9SzmOKejHDNdkxi80ogFo8hq9LKt9NbkHyjOBv81h5zfj/c+b0v6JuP0zU9obiW7DzzRQfnxK+2eB34hffx84OcexD5FjgMnyaQQ4LKt9MdHgvw9425RjNgD7Sjjnt8R9vGlK+2FAP9EA2prjGjvwzintN8btJ5XQ/xm8MFB2A72xjVdM2e9+XviHI+/nJP57+Fn2NS/mc6Wf6T+69SoO4sWtdn1//Psfcmz7LtF/uf9rDvuXwvey38TfMP8vfvuHeY75xynHfMndv19Cnz+a8n5j/PuoEmyUi4Gs190z7ezud7n7+3Js+lX8+zVT2ifP4Xen2BnnhfNciKnn+uPu/tP49bvc/YEcx2wAjjOzrjw273P3vVk29wDPE/l/+5R9NwELzWzJTELNbDnRPwx7gTun6J5sWxnvM5UJ4J+ntM3pc+HuvcBfxm8PxirN7CxgKfCdIsy8HzgS+MhsNIgX0EApisbM5hPdygH4dY5dJtteNZv9Z8HTOdqeiH+flOeYZ2bZ1yTPTXk/OViVFAsrE/OzXvcWc4CZvc3M7orjhhNxLPWpePPCrP26gMkY2xNT7ZD73E8l77k2s8PN7Atm9piZDUzG+4A/mKplCjtytPUDPe6emdLeF/8u5tq8kmgG9sTklG4KhT6re9x9bEpbOT4X1wPDwLlmNvl39DHguhz9HYKZrQSuAT7o7vvnoEGggVIUwOP8LnffGjdN/pc/lucPdTD+vWCW+5fKQI62SZv5ZlhDs+wr5/FZX6qWY99K89L49xjRLcCCmNnHiGY+I8BZRLcRs+8cZPsw07Xrn6k/d895ruMZ3kNEaTifA1Zl3cWYzB/Ndz6HS2wvZCubyc/LYJ7thT6rufyc8+fC3XcR5c02AJeb2WnAy4jisjPxNaIQwj/OuKeYkaZaCxCJYn/8u9nMmnN8gXbEv5+f5f6lkmul5qTNomZYCedN8e+fu3u+L3gAzKyFaDYyAZzv7jOd88nz12xmLe4+OmX7/KkHlMAfE93GvN7di7mFWA32x7878myf62d1tnyBKP/x94nuznzZ3Q8UcdxvAsSz9Fz8nZn9Xfx6TdY/wyIHmlGKoom/jB+P3x6TY5fJtvtns392V0VKylVlZPIW1cNF2qgkxfpRMmbWTbSyE6CYwgKLib7s9+QYJNun7uzu/cCT8dt1OezNpcLL6vj3f+fYNk1LlXiA6J+IdXH6xVTyfVYrSjyAfZdoUnMs0UyxmOPyVdk6pOLTlDtGIg8aKEWpTC4w+N0c236H6MvmpjnsD9F/7W2Tb8zsi2b2lTzHk7WfEaUJOC98IdSSYv0oCTNrB/4/onSGr7r7HUUctofoFuESM1s2Zdtr8xwzeQ7fPaX/RqLzPFsm45snTLHbTLSKuOrEtzlvIVrh+qbsbWZ2GNGt6u1MX7RTDa4hihNfq3hjbdBAKUrlRuBW4MNm9sdm1mFmS83sq0S5cpe5+y/nsD9E/90fbmYnmNkqXqhCMpXlZnaxmc03s8XAV4gWZXwpz4rKalOsHzNiZi1mtsbM3gc8SpT/+Emi9JsZiW+dfo0oZvZdMzs2vhbvIErbyMUXiWbml5rZe8ys3cxWEMXInsxzTDF8g+hW54VmdqGZdcfn5++I8ilrxQeIFjbdbGZviM/50cA/xdt/z91Hqi3K3R9z96Pd/TPV7lvE1Do/RT+1++GFHLCpP2fMcFwj0Rf0I0SzlD6ixPtzyrT/CqIiBfuJlut/B5iftf0bsc6jgE/HfowQ3cq7BLCsfc8oxkdeyP3L/rmaF/L3sn+2TtGR/XNBsX7k8f3qPHrHgB6i3NPPkpXbWsQ1vSDrOnyAKB3kANGM95+JChHkPDdEi1e+RlSYYZhokP594H3xvr+Xte8FebRP0xpfu+8TrWIdIVpV+udEM+XJ4+4ucG3OyNPf1Xmu+d1F/k0sIvoHYTNRkYzdRLc+jy3iOl0db5t2DoroN+dnbBafk5x/u3n2dabk2uon94/FJ1GIxGBRabb3oEUINcPMPkH0Zf0b7n5PjeUIUVF061UIkRMzW2lm/xXHDqfyVqIZ7voqyxKi6migFELkoxk4mShmt9aiBwUfZWZ/Hbd/0GsQsxOi2migFInB4qeHEN12BdhiZltrKCnt7AT+FHgxUc3QfqL0iFXAme7+3QLHCpEaFKMUQgghClCXlXkWL17sq1evrrUMIYQQAfHggw/ucfdpRfTrcqBcvXo1DzwQQpqdEEKIUDCzbbnaFaMUQgghCqCBUgghhCiABkohhBCiABoohRBCiAJooBRCCCEKoIFSCCGEKIAGSiGEEKIAGiiFEEKIAmigFEIIIQqggVIIIYQogAZKIYQQogAaKIUQQogCaKAUQgghCqCBUgghhCiABkohhBCiABoohRBCiAJooBRCCCEKoIFSCCGEKIAGSiGEEKIAGiiFEEKIAmigFEIIIQqggVIIIYQogAZKIYQQogAaKIUQQogCaKAUQgghCqCBUgghhCiABkohhBCiABoohRBCiAJooBRCCCEKoIFSCCGEKIAGSiGEEKIAGiiFEEKIAmigFEIIIQqggVIIIYQogAZKIYQQogAaKIUQQogCaKAUQgghCqCBUgghhCiABkohhBCiABoohRBCiAJooBRCCCEKoIFSCCGEKIAGSiGEEKIAGiiFEEKIAjRVu0Mzux44AAwAJwKXuPvOeNulQBewEPg3d781bj8J+ACwBTgc+Ii7Z8ysDbgO2A4cBVzj7puq65EQQog0U/WBEhh09ysBzOwy4OPAxWb2auD17v5WM2sGHjeze4Be4NvAme6+08y+CLwHuBm4BHja3a81s+PjttdV3yUhhBBppeq3XicHyaz+B+LXbwPujfcZAzYCpwNrgfbJWSfwC+Cc+PU5Wcf8CjjRzLoq6oAQQoi6oiIzSjO7A1iaY9NVWbdTFwBvAs6Ltx1ONDhO0he39QD9Odonj8m1rS+HpouAiwBWrVpVkj9CCCHql4oMlO7+5kLbzawbuBF4r7vvi5t3A/OzduuK2/K1Fzoml6abgJsATj75ZC/KESGEEHVP1W+9mtli4GvApe6+xcwmZ5Q/Bk6L92kCjgXuATYDQ2a2LN7vNcBt8evbso45HnjU3afNJoUQQojZUovFPP8W9/sdM4Po1ukP3H29mf2HmX2OaNXrn7n7fgAzOx/4rJltAxqBb8a2vgxcZ2ZXAkcCF1bVEyGEEKnH3OvvLuTJJ5/sDzzwQK1liBSycUcvt2/Yxfb9Q6xc0M7Zxy3lmOXdOduBoveVjentQpQbM3vQ3U+e1q6BUojysHFHLzfds4Xu9mbmtzXRP5yhd2iMM49Zwl0bew5pf2bfASbcOeKwjhn3lY3pNi46fY0GS1F28g2UqswjRJm4fcMuutub6W5vpsHs4Otv3vv0tPY9AyPsGxwtal/ZmG7j9g27an25RR2hgVKIMrF9/xDz2w4N+89va2JX3/C09pHMOKOZiaL2lY3pNrbvH0KIalGLxTxCBMlc42krF7SzpWeAnf0jDAxn6GxrYtn8VpZ2tbFtz+Ah7eMTDjj3bt47476yMd3GmiWdinOKqqEYpRCUJ7549LIOvn3fM3S0NjG/tZH+kXEGRzKcdewS7ny855D2nv4R3Cc4vKt9xn1lY7qN8099Mb/eOag4pygrWsyThQZKMZXr79xE79AY3e3NB9t6h8Z4fEcfxy7vOqT97ieimhZnrDt82r7Lu1rZ2TdC3/AYXW3NLOtqZUffyLT2/YOjYLBgXsuM+8pGbhulXJep+05e6z8966Xl/zCJxJJvoNStVyGI4ovLu9sOaZuMkb16zaJD2kcy4xiWd9/VizsPtk+488izvdPa73hsB4Zx6trDZtxXNvLbyGam6zK1XXFOUSxazCMEsHJBO/3DmUPa+oczLO1qm9be2tRIS1NDUfvKRrg2Vi5oR4hi0K1XkVpKXZxz7e1PsG9wlNHMBC1NDSzqaOFdr1zJ9x/cfkh7a1MDE+6MjfuM+8pGuDY+evY6YPrin1xtimXWB7r1KuqK7MU5y7vb6B0a46Z7thyysGNqe0NUUhEn+udx8v3U9hfafMZ9ZSNcG5t7BqZ9Fq67Y9PBBUHZnw8t/KlvNKMUqaSUxTmltpe6aEQ2km9DC3/qA1XmEXVFKcn/SUiwl43a2tDCn/pGM0qRKIqNO/b0DzMwnJmWqL6zf+Tg72Lan9ozSFODMTbudLY1ceSSDh57ro8DoxkWdrRWzUau9ucHRwBkowo21izp5OzjlqqYQcpRHmUWGiiTSSlFAR5/rpdtewdZ2NE66yT45w+MsrijhaGxCTpbGzFgYGSctuYG9gyOsnBeS1VshJToX682chU4UDGD9KGBMgsNlMmk1KIAw2Pjc06Cb2tp5LiVXTy5e5C+4TFaGhuYmHDWLumomo3QEv3r1UYpsW3FNJOJVr2KxFNqUYAGm3sS/GhmgsWdbSzubDu4762PPseqwzqqZiPERP96tZGNihnUD5pRippTbHHr2zfsyll0PIkxr1JjpSHrk43p7Z1tTSyZ36Z4ZsLQrdcsNFCGQ664o4qOJ0ufbEyPSx+xaB7HruhWPDNhaKDMQgNlOOSKO6roeLL0yUbuuLSKsycPxShFkOSKO6roeLL0yUbuuHQ2imcmG80oRdXIFZ/JFXcMJdYUig3lYqbbRqEcTVFdVJlH1JTJWGTv0NghNTTntRgPP7OfvqExOloa6BsaY9/gGHsHRg5pe/iZ/Ry3onPavvna02Rj/dZ94E7/cIamBhgZzbB+8z4y4xP0DWdqrk825mZjXovl/NvYuKO31n+2IkYzSlEVCuVAhhxrCsWGcjHTbUOxyzBQjFLUlJlyIEONNYViQ7mY6beRjWKXYaFbr6Iq6MHItbURuj7Z0IOlQ0a3XsWs0YORk2MjdH1psxHlVTaxs2+YfYNjNDUY575iOfdtfr7o66UHS1cf3XoVZUUPRk6WjdD1pcnG/LZmTlmzgFse3sHY+ASHdbSwvLuN+7fsZ3AkU5QNPVg6LDSjFLNCD0ZOlo3Q9aXJhh4snVyUHiLKih6MnCwboetLkw09WDp9aEYpiiLkByOHkjgeso3Q9YVkY15LEy9b0cWTPYMMDGdobjQyE85LFncEe65VtKA8aEYpZk2uYgFP7h5g/Za9SqRPiI3Q9YVio384w1hmnPWb9zEyGn0u+4czGD7nz3slfVTRgsqiGaWYkXyFy0N4MHJoieOh2ghdXyg2lnW1srlnkIYGY3R8gq62Zo48vIMN2/vK8nmv9oOlFbssDa16FbMmX+HyEB6MHGLieIg2QtcXio3J9necuOLgSlSAkczzZfm8V/vB0opdlgfNKMUhJK1wueKcilHqPOW3odhlaShGKWYkiYXLFedUjFLnSQXXK41mlOIgSS1crjinfNR5UsH1cqAYpZiRpBYuV5xTPuo8qeB6JdGMsk4pNhYZUrxFcSn5mMYcyFr4qNhlbhSjFAcpJRYZUrxFcSn5WKg9qTmQtfBRscvS0IyyDiklFhlavEVxKfmYxhzIWvio2OV0FKMUByklFhlivEVxKfmYxhzIWviYjWKX+dGt1zqklIcoh/6AWz3UWD4mRV/oPuph0fnRrdeUM9eHKC/qaOHUtQu55aEdZCacRR3NLOtqY3B0nN6h0SAekluOhxrXg48h65OPtfdxsn3TrsG6XeCjxTx1SL5FO5t7Bop+CO3gSIb7Nj/PumWdLOpoZu/gKJt2D3DKmgV0t7cUZWPmh+TW1kY9+Bi6PvlYe9uDIxm+/+B2LfDJQdVnlGZ2PXAAGABOBC5x951m9lbgt4DHgBOAH7j7j+JjfggsyDLzW+7eY2ZtwHXAduAo4Bp33zSThnqZUZbj4cpJfZCtfEyWPvkYro/1tMAnpBnloLt/3N0/DzwMfDxufzFwlbtfB1wK/L2ZTep7xN3PyPrpidsvAZ6ObV0P3Fw9N8KnHA9XTuqDbOVjsvTJx3B91AKfGscozewKoMvdr5jSvhx40N1XxO9vBp4CjGhg/Fbc/jPgY+7+s/h9H/Aid+/L0ddFwEUAq1ateuW2bdsq51iNqMTDlZXsng4fQ9dXLz6GXhAhl76XrejiD047oi6KE+SbUVZkoDSzO4ClOTZd5e63xvssAG4BznP3fVOOvx74ubv/IH5/irvfH7/+HvBjd/+WmT0BvNvdH4m3PQuc4e5PFtKXxluvk/HI7vZm5rc10T+c4fHnetm2d5CFHa3Mb22kf2ScwZEMZx27hDsf76GjtWnG9p7+EdwnOLyrfdY2Kmk7FBvSV3sboevbf2CURfOaGck4na2NODAwMs68lgZ6+kdq/neaT9/izhY6Wps44rCOg98tvUNjXHT6mtQNllUdKIsQ0w38JfBxd98yZdtHgH3u/vU8x74feLW7/1EpM8ps0jhQVvLhykp2T76PoeurBx9DL4hQSB/UR+wymIIDZrYY+BJwqbtvN7PzsmaOVwL/7e7fM7MzgF8Bo8CfuvunYhNHAZMzxtuA04CfmdnxwKMzDZJppZIPV1aye/J9DF1fPfgYekGEQvqMF95D/cUua7Hq9SGiAXrydmu/u7/dzD4IXAk8HrevBM4CdgLfIVoNa3H7n7j7ATNrJ1r1ugM4Evhcva56/dg//ZLHnusL8uHFy+an/+HK0ld7G9JXXX1rlnTWzYxSBQdSwMYdvVx3xya27Bmks7URI4ottDU3sGdwlIXzWmoa+3j+wCiLO1oYGpsIUl89xMdC11cPPqZN3+VvWcc5J6ys9ddfWdFAmUXaBsrJ+OTY+HiwDy9O+8OVpa/2NqSvuvpWL66fGaWKoqeAyfhkgzUH+/DitD9cWfpqb0P6qquvnmKUGigTyNR8yZZGY9uewWnxhqVdbSW1d7Y2cc+mnkPiiK1NjRwYzXDv5r2ztj0+4YDPyUYlbc/GRv9w5pAVxtmFprPbW5sap12/fPtW0kbo+srpY6U+I/NaGunpHz4kx7CztSmYz3C19a1Z0pmzlnTaUkZAt14TR6XyJUOPI4YU3zn/1Bdz3+bnVUy7zDai89zEzr5h9g2O0dRgnPuK5UWf68ni9t++7xnlQFZB3/mnvphf7xw85Lso6fmVuvWaEm7fsIvu9uaD/0l3tzczkplgwbwWutub6RuO8pvWLe1kw3MDvGLVgoPxhkLtPuE0NTZw2qoFB+OInW1NTEw4p65ZVJSNQrYx5qSvkrZLtfH03qFZFKv2MhS8nr2N0PXNb2vmlDULuOXhHYyNT3BYRwvLu9u4f8t+BkcyRfu4be9QxT5/m3sG6Wx7Icfw+Bd1s2F7X1n+9pKo7ye/3nNIvdjJ37dv2JXYgTIfGigTRiXzJUOOI4YU3/nJxt2csmYRx6184cugd2iMb977NMcu7zqkPV+h6Vz7VtJGEvT95Nd7OHXtYUUV6s5nY/2WvZx5zFLlQFZRXzZpza/UY7YSRq6HLtfLQ2VDeaCu40EUqw69mHYt9BmW+s9f6PrS+PBnxSgDptiHLocSZ6qXhysv727LWWw+5KIKaUt2z2djflsTz/UOJ+qznQR9pcSPP3r2OoBELvIJ6TFboghKeegyhPHg13p5uPIRh7Xz8DP76Rsao6Olgb6hMdZv3Qfu9A9naGqAkdEM6zfvIzM+Qd9w5pB9H35mP8et6JxmI1/7vsEx9g6MzMlGJW2HpG/VYe2J+2yHrm9+WzOnrl3Ipt0D7B0cZVFHM+uWdeaNH2/uGcj53ZXkB0BrRhkopTx0uZ4fKlsrfUkrqpC2ZPdCNkL/7CRNX7lsJ6GIula9Joxci3YmYzNTA+gjmfGcRYtz7VtJG/WmT4uhwtWXTYifnSTpK5ftJC/y0YwyUK6/cxNbegYSFwcLXV89xOnke/L0hfxA53LZTkIRdcUoE8ZLl3YkLg4Wur56idPJ92Tp6x/OMJYZZ/3mfYyMRn87/cMZDGf9lr0111cu2y9d2lHrr9VZoxlloFx/5ya27hlIVBwsdH31FKeT78nRl4QHOpfDdhKKqCtGmTC27x+aFtsKPQ4Wur56i9PJ92ToS0oxg7naVowyYYQ2o8yVL3n7hl3TYpSKI4ZrO3R99ex7vcQAQ/d9zZLOg99toeZXKkYZKPnyJee1mOKICbIdur569r1eYoCh+z6vxRKbX6kZZY0plC+pOGJybIeur959r4cYYBJ8Dz2/UjHKQJkpX1JxxGTYDl2ffE9/DDApvmeTlPxK3XqtMbmKnIdU4Fi206FPvsv3UH1PQhF13XqtMRt39E4rch5SAfC02dZ5rb+i3vI9jGLp+Wx/9Ox1wSzo0WKegAm5AHiabOu81mdRb/keRrH0QrZDRzPKGpNrMU9IhYzTZDt0fUm1Hbo++R627SQs5tGMssZs3z8U9INw02Q7dH1JtR26Pvketu0kLObRjLKKVLuwQBITpFUQIXm2Q9dXLwn9SbUdUiECzShrTLULCyQ1QVoFEZJnO3R99ZLQn1TbSShEoBlllah2YYGkJkirIELybIeur54S+pNqO5TYpQoO1JhqFxZIaoK0CiIkz3bo+uotoT+ptrMJLXapGWWVyPcg5iTG+hRHlO0k6Uuq7UrHP0O3XYsHPStGWWNyPYg5qbE+xRFlO0n6kmq7kvHPJNgO6UHPmlFWiVwPYk5yrE9xRNlOir4k265U/DMJtmvxoGfFKGtMrgcxJznWpziibCdFX9JtVyL+mQTbIcUodeu1Sqj4uWyn1Xbo+mQ7mbZDKpauW68VIFdhAUDFz2U7lbZD1yfbybT90bPXAVS1EIEW81SJfIUFNvcMqEi3bKfSduj6ZDuZtjf3DARTiEAzyjJTqLBAyIWJZVu206pPttNlu5KFCDSjrBK5ipwnoTCxbMt2WvXJdrps12KRj2aUZaYWhQVkW7ZrnRgfsj7ZTpftShYi0IyySlS7sIBsy3atbYeuT7bTZbsWhQg0oywz1S4sINuyXWvboeuT7XTZrmQhAhUcqBLVLiwg27Jda9uh65PtdNmuRYxSA2WZWbmgPWeMsrO1iXs29RxSRLy1qZEDoxnu3bz3kH2XdrWxbc/gNBv52mVbtmtpe3zCmdfSSE//8CFFsDtbm0rSJ9uyXYztNUs6Z/4iLjO69Vpmbvvldq751yfoaG1ifmsj/SPjPH9glMUdLQyNTdDZ2ogBAyPjtDU3sGdwlIXzWg7uOziS4axjl3Dn4z2H2MjXLtuyXWvb+w+MsmheMyMZp7O1EY9tz2tpoKd/hIUdrbIt22Wzfflb1nHOCSsr8v2d79Zr1QdKM7seOAAMACcCl7j7TjNbDdwO7Ix3fdDdPxwfcxLwAWALcDjwEXfPmFkbcB2wHTgKuMbdN82kodoxyiQXKJdt2U5DgW3ZTo/teolRDrr7lQBmdhnwceDieNs17v6N7J3NzIBvA2fGA+oXgfcANwOXAE+7+7Vmdnzc9rqqeJGHXDHKJBcol23ZLrY95ALbsp0e27WIUVY9PWRykMzqfyDr/dvN7FIz+7SZHRu3rQXa3X1ypvkL4Jz49TnAvbHdXwEnmllX5dTPTK7i5yEVGpZt2S637bT5I9th265FsfSK3Ho1szuApTk2XeXut8b7LABuAc5z931m1gGsdvfHzGwpcB/wcuAY4C/d/aT4uDOBq939tWb2BPBud38k3vYscIa7P5lD00XARQCrVq165bZt2+bsZ7HFz0MqNCzbsl1u22nzR7bDtl3JYulVvfXq7m+eQUw3cCPwXnffFx8zCDwWv95lZruIYpjPAvOzDu8CdsevdxfYNlXTTcBNEMUoS3RpGpPFz7vbmw8p2HvmMUsKFAP2ChYalm3Zro3ttPkj22Hb3twzwF0be6Z99150+pqKPVmkFot5FgNfAi5z9+1mdp67/8DM/hB42N1/ZWbNRAt3XgM8DfyKQ2OUj7v7zWZ2OTCRFaO80d1njFGWYzGPip/Ltmyn0x/ZTqbtchRLD6mE3b8BJwDfMbO7gQvi9meBK83sCuCvgD93920ejeTnA581s6uARuCb8TFfBo4wsyuBDwMXVssJFT+XbdlOpz+ynUzblVzkozzKWZKv+PlTewZparBDkrofe66PA6OZORcDrqTtpBZIlu3a206bP7Kd2/a8liZetqLrkGIBmQnnJYs7grBdjmLpIc0oU0Gu4ufrt+4Dd/qHMzQ1wMhohvWb95EZn6BvODOnYsCVtJ3kAsmyXXvbafNHtqe39w9nGMuMs37zPkZGo++g/uEMhrN+y94gbFeyWLpmlLMkTYUFklwgWbZrbztt/sh2/RYiCKngQCpIU2GBJBdIlu3a206bP7KtQgRT0YxyluSKUSY1/qDYqmwrRpkO25WKIybhXClGGSC5YpRJjT8otirbSetTtqsbR0zCuVKMssxUMkaZxPiDYquyreuQfNuVjiOGfq4UowyQfDHKJMYfFFuVbV2H5NuuRhwx5HNVyRilbr3OEhU/r73ttPmTVNtp8yepttPmT0jF0nXrdZbc9svt3PCTp8hMOIs6mlnW1cbg6Di9Q6NzLhx86tqF3PLQjorYDr3gceh9ynYYfSbZdvQg4iZ29g2zb3CMpgbj3Fcs577Nz+s6zOFcffTsdXOu9arFPGVk445e7trYw7plnSzqaGbv4Cibdg9wypoFdLe3ALMvEDw4kuG+zc9XxPbMRYmTZTtt/iTVdtr8qaTt+W3NnLp2IZt2D7B3cJRFHc2sW9bJ/Vv2MziSSZw/oZ2rSqEZ5SwopSB6SMWHk1zwOJQ+ZTuMPpNqO23+hHSualIU3cxWzfCzZE6KEkwpBdFDKj6c5ILHofQp22H0mVTbafMnpHNVk4IDZrYZ2Arkm9Pudvd3V0hXRSnHjDJXQfTQE3jTlvyfNn+Sajtt/tR74n5Sz1WtCg7c7O5vcPfX5/oBHpqTogSTq9hAEhJ405b8nzZ/kmo7bf7Ue+J+Us9VTQoOmNk8dz+Q98AZtodMOWaUU4sNJCWBN23J/2nzJ4m20+hPPSfuJ/Vc1aTgQL5B0Mze6+5fT+ogWQ5yFRtISgJv2pL/0+ZPEm2n1Z96TdxP6rmqacEBM/ukmW0zs81mtgX4QsXUJIRcxQbSmMAbep/yJwzb8qd+bYfkT00LDpjZPwPv8HhHM/uf7v7DiimqAnO99bpxRy/X3v7EnJOD01RYIG3Jzmm7Prr2YfcZUuJ+Uv25+I0v4ZwTVs5pbJhLwYGHgbas9/WXeJmDuSbTpq2wQNqSndN2fSppW/6EbTttRQ7y+XPXxh427uilEhQzo7wYuAbYRZQq0uXuhxU8KHBCKDgQUqJuUvtMqm35E7Zt+RO27ULtcy06MJcZ5W8CK9x9rbuvAf5s1ipSQjkKDoSUqJvUPpNqW/6EbVv+hG27UHulFvQUM6P8tLv/edb733D3n1ZETZUIoeBA2hLAlfwfdp9p9CeURHf5E4Y/5Sg6MJcZ5VlmtsXM/t3M/gP4+qxVpIRyFBxIWwK4kv/D7jNt/oSU6C5/wvCnkkUHiplRfhe4LKvpfdkzzCQSQsGBNCaAK/k/3D7T5k9Iie7yJwx/ylF0oOSCA5O4++9MaUr0IFkOylFwII0J4Er+D7fPtPkTUqK7/AnDn0oWHShUwu4id78p74FmF7r7zRVRVWGqHaOsduwubfHCNPqjGGXy+kyjP7WIi1ayz0rFKAsNlD8BCi3aWeHu75+1ohoy14Hytl9u55p/fSJOem2kf2ScwZEMZx27hDsf7zmk/fkDoyzuaGFobILO1kYMGBgZp625gT2Doyyc1zKjjVLba9FnT/8I7hMc3tUufwLsM23+6BzO3fb+A6MsmtfMSMbpbG3E47+reS0N9PSPsLCjNXF9Xv6WdXMqOjCbgfI9M9h8zt3vnLWiGlLNGGWtYndpixemyR/FKJPZZ9r8qVUcsZJ9Vj1G6e7fnHVvKaeUGGWtYndpixemyR/FKJPZZ9r8qVUcsZJ9VipGWUx6iJhCKUXR01aQOKQiyEn1R+cwmX2mzZ809lmpwugzpoekkWoWRU9b8ea0FQuvlz7T5k9IxbiTWlw8jX1+9Ox1HLO8e9bf7XMpOJBtZMGsFaSMEIoDq1h4Ms9h2opu10ufaSsuntY+K0ExBQduBP4eeCXwEeAH7v6RiimqAtUsiq5ix+ozhD7T5k+99Jk2fyrdZy2Lom9z9/uAPwBeBvTOWkVKKKUouoodq88Q+kybP/XSZ9r8qXSflVrMM2NlHmCxmb0OeMrdD1gFp7dJYeWC9pwFBzpbm7hnU88hSfGtTY0cGM1w7+a9h+y7tKuNbXsGp9kotb3afY5POOAV80d9VqbPtPlTL32OTzjzWhrp6R8+JEG/s7Wpov4ktc81Szrzf3HPgWJmlFuArwDXmtnbgBdXREmCyFUUvRZFwdNWiFx9qii6+qx9cfEk91mzoujTDjA72t1/XRE1VaISBQdUWCC5Sdf10Gfa/KmnPtNWFKCSfVa94ICZvc3df2xmV03ZdDpw5qyVpIBcBQdUWCDZSddp7zNt/tRbn2kqClDJPmtRcOBV8e+XA9uyfvZXREmCyFVwIMlJuuoz/X2mzR/1qT6DKjhgZi9292ey3r/G3X9RETVVohxF0W/4yVNVS7hPY2Jwvfepog3qU32Wv89aFhw4+NBmM2sGPjxrFSlg445e7trYU9WE+7QmBtdrnyraoD7VZ+X6rATFzCh/BPwK+A7wd0C7u59YMUVVYC4zylKKDSQlSVd9KolcfarPNPRZy4IDFwH7gEeJ0kReO2sVKaCUYgNJSdJVn0oiV5/qMw19VmoxTzEzym3ARuB9wIXAie7+zoqoqRJznVHmKjZQyaecL5vfylN7BmlqsEMKCzz2XB8HRjMVe+K6+qxMn5V8ar36rE2f81qaeNmKrkOS5TMTzksWd6jPKva5Zkll0kOKmVH+HHiLuz/r7p8Ets5aRSTkejP7rJldYWbfNbNlcfsFZvagmd0d/2w2swvibT/Mar/bzJbE7W1m9tXY1tfNbPZnqEhyFRuodEJzvRQzqJc+05YUX+99JjlBP2191qzggJl1u3tv1vtXuft/zbpDs8+4+5Xx68uAF7n7xWZ2GrDd3Z+Ot/0I+H13HzCzq9396hy2Lgcm3P1aMzseuNHdXzeThrnOKKcWG6hGQnO9FDOolz7TlhRfz30mOUE/bX1WveBAFvPN7OPAYsCA44FphoplcpCMaQAG4vZ7s8SeBGxy94G46cVm9rG4/6fd/Vtx+znAx+Ljf2VmJ5pZl7v3zVbfTOQqNlCNhOZ6KWZQL32mMSm+XvtMcoJ+2vqsZYzy28APgLOB7wO/7e5/PMMxdwBLc2y6yt1vjfdZANwCnOfu+6YcfzPwKXffFr8/xd3vj19/D/ixu3/LzJ4A3u3uj8TbngXOcPcnc2i6iGhhEqtWrXrltm3bCvqdj0rHKGsRowsp1qM+09Vn6DGtNJ7zeukz13V+2YouPvfOE2b13Q5zi1E+4u63AJvd/U7g6ZkOcPc3u/tJOX4mB8lu4EbgvTkGyaVA6+QgGdu7P2uX/wDeEL/eDczP2tYVt+XSdJO7n+zuJy9ZsmRGp/NRyRhlLeJlIcV61Ge6+kxCTCtt57xe+sx3nXf1jbBxR+/MX+QlUsyM8h+Ay4GLgUeIBrc3FDyosL3FwJeAy9x9u5md5+4/yNr+SeBOd/95/H4+8Kfu/qn4/ReBfe7+2bTFKGsRLwsp1qM+09VnEmJaaTvn9dJnvuvc3Ng4p1zKucQovwR0An8JfAH46qwUvMC/xf1+J362ZT/RrV3MrBU42d0/kbX/GHCimX2KKEa5EPjzeNuXgevM7ErgSKL0lYpSyRhlLeJlIcV61Ge6+kxCTCtt57xe+sx3nSsVp5zx1qu7r3f3x9z9KXd/p7v/01w6dPdXuPsJ7n5G/PP2rG0j7n7OlP2H3f08d7/K3f/c3d/r7gfibUPu/gF3/4y7X+Dum+airRhyFURPctHhkLSoz3T1GZIW9ZmuPqtdGL3k51Gmgbncet24o5drb38iiALASS3QXas+Qznn9dJnSFrqqc+O1ibmtzaxs2+YfYNjNDUY575iOfdtfj41fVaqMPqsF/OY2Ttm1WOKCakAcBILdNeiz5DOeb30GZKWeulzflszp65dyKbdA+wdHGVRRzPrlnVy/5b9DI5kUtNntQujF7OY537gbuDv3X1DRVRUmVCLoieh6HBS+wxJS730GZKWeukzJC218r8Si3mKSQ95D3A1cLqZfc3M3j7D/qmmkkXRk1B0OKl9hqSlXvoMSUu99BmSllr5X4nFPMXMKI8D/hs4H/gTYDuwBfiZu/+/siuqAiEURVdhger2GZKWeukzJC0hJcXXU/GDWvg/l8Loc5lRfht4HDiBqCrP29z9YuC4WSlJOOUoOKDCAtXvMyQt9dJnSFpCSoqvl+IHtfK/EoXRi5lRfg94n7v3Z7W1AF+MB8zEUeuCAyosUJsE5VC01EufIWkJKSm+Xoof1Mr/uRRGn0vBgd9z9/HsBncfJarUU3eUo+CACgvULkE5BC310mdIWkJKiq+X4ge18r9WBQfGZ9qnnihHwYGQEnTrPUG5XvzXOZf/9eK/Cg6UiWoWHAglsT4kLfVUWCEU/1VwQP6HUhSg0v7XpOCAmE4SE+tD0lJJ2/I/nD5D0lLv/odUFKAa/pcbzShLpJSCA/WecBySFvkv/+vZ/5C0VLrPWhUcEFmUUnCg3hOOQ9Ii/+V/PfsfkpZK91mTggNppFIzyhCSfENKOK73wgohaZH/9Z38H5KWcl4LzSgD5ezjltI7NEbv0BgT7gdfv/HoxUEk+YaUcFzvhRVC0lLv/td78n9IWsp1Ld549OKc38VnH7e07N/7mlHOgo07erl9wy627x9i5YJ2zj5uKbdv2DWtEIGSvOu7sEJIWurdfyX/h6WlHNdi9eLOg9+92d/Fs13xCnMrOCCmcMzy7mkX4+afb51WLEBJ3vVdWCEkLfXuv5L/w9My12uxff9Qzu/iSqBbr2UiVyGCekryDaVPaQmnT2kJp880aqlEYYF86NZrmbjtl9u54SdP1TyZvd4T60PTomshLaEn/ydVy1wKC+RDi3kqyMYdvdy1sSeIZPZ6T6wPSYuuhbTkag8p+T/pWqqFZpRloJQiBElP5pWWsPuUFmmpFy1zSQPJh2aUFaSUIgRJT+aVlrD7lBZpqRctlSgskA/NKMvA9XduYkvPQFUThUNK5peWcJK5pUVa6kXLmiWzf+5kPjSjrCAvXdpR9UThkJL5pSWcZG5pkZZ60fLSpR1V+47XjLIMXH/npmnFBuotmV9awknmlhZpqQctqxdXb0apggNlYPv+oZok0IeUzC8tYSVzS4u0pF1LNWOUuvVaBnIVGwgtOVdapEVapCVNWlRwoMKU+9brxh293HTPFrrbm5nf1kT/cIbeoTHOPGZJqhKlpUVapCWchPt603LEYR2HfL9edPoaFRxIEscs7+ai09fQ3d7Mjt5hutubuej0Naxd0pmqRGlpkRZpCSvhvl60vOuVK6d9v1ajxuskmlFWkEoWIkhCQrC0SIu0pEtjrbRUorhALjSjrAGVLESQhIRgaZEWaUmXxlppqebCnVxoRllBylWIIJQE+tC0JDVRWlrC0jKvpYmXrejiyZ5BBoYzNDcamQnnJYs7qqpF5yt/eyWKC+RCM8oaUI5CBCEl0IekJcmJ0tISjpb+4QxjmXHWb97HyGj0ue4fzmA467fsDUKjzld1iwvkQjPKClKOQgShJdCHoiXJidLSEo6WZV2tbO4ZpKHBGB2foKutmSMP72DD9j6Gx8aD0KjzVZniArlQwYEaUI5CBKEl0IeiJcmJ0tISjpbJ9necuOKQRzeNZJ6nwcLQqPNV3eICudCt1wpSjkIESUj8rYWWJGiUlvC1JEFjSFpqpbGaxQVyoVuvFaTUQgSnrl3ILQ/tIDPhLOpoZllXG4Oj4/QOjQaTEB1KcrbOV7K1hJLMrvMV3vm6a2PPtO/MauVNajFPDSilEMHgSIb7Nj/PumWdLOpoZu/gKJt2D3DKmgV0t7ccsq+Ss3W+kqwlpGR2na+wztfaJZ05vzOrWVwgF5pR1oBchQhCSioOLdk4dC1J0BiSliRoDElLEjSWS0u1CgvkQzPKgMhViCCkpOLQko1D15IEjSFpSYLGkLQkQWO5tNR60U4+NKOsAR/7p1/y2HN9wSbtq7BA+jSGpCUJGkPSkgSN5SpEUK3CAvnQjDIQNu7oZVffSNBJ+yoskD6NIWlJgsaQtCRBY7kKEdS6sEA+NKOsMpPxybHx8aCT9lVYIF0aQ9KSBI0haUmKxnIUIqhWYYF8qOBAIGzfP8Ty7jYarDnopH0VFkiXxpC0JEFjSFqSpHGuhQhCjVFW/darmX3IzP7WzD5qZj8ys9Oytl1qZp82s6+a2Tuy2k8ys78xs4+Z2ZfMrClub4v3vcLMvm5mtftXpEjKUYQg6cnD0ljfWpKgMSQt9aSx1oUF8lH1W69m9lHgBncfMrNzgf/j7meZ2auBT7j7W82sGXgceBXQC/wKONPdd5rZF4HH3f1mM7scmHD3a83seOBGd3/dTBpqeeu11CIE9Z60L43hJYTXg8akJPknWeOvdw7WrLBAPoJZzOPu17r75Pz6SKIBEeBtwL3xPmPARuB0YC3Q7u474/1+AZwTvz4n65hfASeaWVeufs3sIjN7wMwe6OnpKbNXxVNKEQIl7UtjaAnh9aIxCUn+Sdb4652DnHnMkuAKC+SjIjNKM7sDWJpj01XufquZLQOuAF4OvNPd95jZXwMb3f1LsY1vA3cDjwF/6e4nxe1nAle7+2vN7Ang3e7+SLztWeAMd3+ykL5ap4fkIlcRgqQkCUujNEqjNCatuEAuqjqjdPc3u/tJOX5ujbfvdPcPAVcD/xIfthuYn2WmK27L117omMSRqwhBUpKEpVEapVEa01JcIBdVX/VqZpe6+xfit1uIbq0C/Bj4RLxPE3AscA9RjHLIzJbFt19fA9wWH3MbcBrwszhG+ai791XHk/KyckE7W3oGpiXgLu1qY9uewYq0j0844Ny7eW/V+pRGaZRGaZwsLpAUalFwYJWZfdHMrgA+D7wPwN3XA/9hZp8DbgD+zN33e3Rv+Hzgs2Z2FdAIfDO29WXgCDO7EvgwcGGVfSkbL13aUdcJy9IojdJYPxpDLi6QCxUcCITr79zE1j0DdZ2wLI3SKI31oTGE4gK5UMGBwNm+f6jqyfwhJixLozRKY/o1hlxcIBeq9RoItShEkMaEZWmURmkMX2PIxQVyoVuvgVDpQgShJ8QnRWPoifLSqCT/EDXetbEnuOICuQim4IDITSULESQhIT4JGpOQKC+NSvIPTePaJZ05v9tCGyQLoRll4JSjEEESEpCToDHJ2qVR2mulMcTCAvnQjDKhlKMQQRISkJOgMcnapVHaa6UxSYt28qEZZeAUmlHmekL4U3sGaWowxsadzrYmjlzSwWPP9XFgNBPEk9CT8LT2NGpPisZ5LU28bEUXT/YMMjCcobnRyEw4L1ncEYRGnd/C7ZpRippw9nFL6R0ao3dojAn3g6/fePTiaUm867fuA3f6h6MnjI+MZli/eR+Z8Qn6hjPBJBonOUk6qdqToLF/OMNYZpz1m/cxMhp9hvuHMxjO+i17g9Co85u//Y1HL875XXX2cbnKficLzSgTwMYdvdy+YRfb9w+xckE7Zx+3lNs37JpWoGD/4ChtLY0ct7KLJ3cP0jc8RktjAxMTztolHcEkGic5STqp2pOgcVlXK5t7BmloMEbHJ+hqa+bIwzvYsL2P4bHxoLXr/EYFBCa/m7K/qxK1aEcFB5LLMcu7p33Ybv751mkFCu54bAejmQkWd7axuLMNiBJ7b330uaoXM0hrknRStSdB42T7O05ccXD1JMBI5nkaLGztOr9RAYFc31VpQLdeE0quAgVJSTROcpJ0UrUnQWOStSdBY6W1J6mAQKno1mtCue2X27nhJ08lLjk/yYUFkqw9Ccn8SdYeYpJ/tbV/9Ox1iZ9NajFPiti4o5e7NvYkMjk/yYUFkqw9JC1p1B5akn+ttKcVzSgTSDmKENRb0rO0S7u0V1Z7ktJA8qEZZYooRxGCekt6lnZpl/bKak9DYYF8aEaZQEotQhBa0jOEnZgt7WG1S3t42jWjFMFTShGC0JKek5CYLe1htUt7WNrTXFggH5pRJpRiixCElvSclMRsaQ+nXdrD0p6GwgL5UMGBlFFsEYJaJzcnOTFb2sNol/awtKe5sEA+dOs1ReQqQpDW5GZpl/aQNNaT9jQXFsiHbr2miI07ernpni3TniR+5jFLlJxfJ4nv0h5WexoKERxxWMch3ydJe+hyKWgxTx1wzPLunE8SX7ukU8n5ZWgPSYu0h6UxrYUI3vXKldO+T9I6SBZCM8o6oBYFCpKcaJ1Gn6Q9rPak+JSGlI9S0IyyjqlFgYIkJ1qn0SdpD6s9KT6luYhAKWjVax2wckH7tBnlZAB/257BacnHpbZ3tjZxz6YexsadzrYmjlzSQWtTIwdGM9y7ee+cbNeqPW0+jU844InVPq+lkZ7+YZ7sGWRgOENzo9HZ2hSMxjRcj/7hzLTviHpcuJMLzSjrgEoWKFi/dR+40z8cPTF9ZDTD+s37yIxP0DecCSJxWj4lO/G9fzjDWGac9Zv3MTIaXZP+4QyGs37L3iA0Jv161GMRgVJQjLJOqFSBgv2Do7S1NHLcyi6e3D1I3/AYLY0NTEw4a5d0BJE4LZ+Snfi+rKuVzT2DNDQYo+MTdLU1c+ThHWzY3sfw2HgifQrteqS5iEApqOBAnVOpAgV3PLaD0cwEizvbWNzZdnDfWx99LvjiB/XmU1IT3yfb33HiikMe5zSSeZ4GS6ZPoV2PeiwiUAq69VrHlKNAQZITquVTWBrryafQtCsWWRjdeq1jylGgIMnJ4PIpLI315FMtCxHctbFn2t98veZHTkXpIWIa5ShQkKsttMTpekpwl09haQytEMHaJZ05/+Y1SBZGM0oxjVIKFCQlcTqNyeDyST6pgEB50YxSFE0pBQqSkjidxmRw+SSfVECgOmjVq5hGvgIFaUvCT2NhgSQluMunyhZXUAGB8qEZpZhGrgIF2/YO0tHSmKok/DQWFkhKgrt8qmxxBRUQKC+KUYqcTC1QsKd/mOamRsbGx1OVhJ+2wgJJSXCXT5UtrqACArNDBQdESUxNPv7IPz7Kos4mGqw5VUn4aSsskJQEd/lU2eIKKiBQXnTrVRRFOYoThNYeWtK3fJVP5fJVscjyoluvoigKFSfIlcBcStECJazL16Ql7dfK1yMO6yjq70+5kbNDt17FnJgsTpAd83j3q17EMcu7Wbukc1o7hJPcPXPCutdci3ydXdL+KWsWcMvDOxgbn+CwjhaWd7dVPGm/Vr6+8xUr2LRrsKi/Pw2S5UUzSlERSilakNbk7pDa5WtYGlUsIExUcEBUlVKKFqQ1uTukdvkalkYVC0gWuvUqKkK+ogX5EqFr0Z7GwgL1lJxfTl9D+UwWatcCndpR9VuvZvYh4HhgE/Aa4Bp3v9fM3gr8FvAYcALwA3f/UXzMD4EFWWZ+y917zKwNuA7YDhwV29o0kwbdeq085XgySTnbT127kFse2kFmwlnU0cyyrjZ29g2zbd8BFs5rYX5rI/0j4wyOZDjr2CXc+XhPvFAkHe09/SO4T3B4V3vNtYTm6/mnvji4xT9aoFMbQlrM0wpc7O5DZnYu8CngLODFwFXu/oyZLQU2mdlCd58AHnH3q3PYugR42t2vNbPjgZuB11XFC1GQfIt/oPoLIQZHMty3+XnWLetkR+8wewdH6RvOsPaweSzvbjuYrN3d3sy6pZ1seG6AV6xakKp2n3Aw6G5vrrmW0Hx9eu9QMIt2Dj7hQwt0gqKmi3nM7FLgRe7+oSnty4EH3X1F/P5m4CnAiAbGb8XtPwM+5u4/i9/3xfb6CvWrGWXtqMUin3yLO9Zv2cuZxyw9JLF7suDA1ITvpLdPJue/6WXLaq4lNF9/snE3p6xZpEU7orqLeczsDjN7JMfPO+Lty8zsy8DbgU/nMPFR4OKs93/t7p9z988CbzOzP4jbDwf6s/bri9tyabrIzB4wswd6enrm6qKYJbVY5JNvcYdhwSSUh5awnuT2Un11XIt2REFqPaN8A1Fc8ZSsto8A+9z963mOeT/wanf/I80ok8dsZpTL5rdOW4Cxs3+k6PbnB0cAWNjResi+89uaeK53OKjYlAoOVN/X5d1tDAxn5vQZm6ldM8pkEEx6SHy7dZItwNqsbVcCz7j7183sDDM7zMzmm9lVWcccBTwZv74NOC0+9njg0ZkGSVFbcj2ZpHdojPectipn+xuPXlyxp0ysOqw9qNhUJdtD0hKar0cc1l7Rp4Tke5JHvs+8nvARHrVY9XoDMArsAU4EvuvuPzSzDwJXAo/Hu64kWuSzE/gO0WpYi9v/xN0PmFk70arXHcCRwOe06jV8pj6ZZPKpBrnab9+wi617Bir2lImQEsqVhF87Xyv5NJBCT/LI97cgakMwq17d/eI87V8BvpLnsPPyHDMEfKBM0kSVyPdUg1ztN/9865yf5DHTUyaymYwdpa19JDOOYUXtm/T22fpaqaeBFHqSh57wkQxUmUcETTmeWlJPC1l0DsLzVYUCko9qvYqgKfWpJUcv65hWWGBwdJzeodG6WMiixTzVeapIvifmqFBAsglmMY8QpTBZuKC7vZkdvcN0tzdz0elrOOeEldPazzxmCb/eOci6ZZ0s6mhm7+Aom3YPcMqaBXS3twC1X2iSlAUuSW6f39bMqWsXsmn3AHsHR1nU0cy6ZZ1le6rI2iWdRX8mNUimA80oRWpIwhNLkrTAJW3t5ToHSt9IL5pRitSThCeW6IkayT8HKghQf+jpISI1zOaJJaE88SLEJ2qkrX18wpnX0khP/zBP9gwyMJyhudHobG3SUzxEQTSjFKmhFsUMktKer+hCSBor3d4/nGEsM876zfsYGc3Q1BANfB0tjWzbO6iCACIvilGKVFHtYgZJaS9UdCEUjdVo39wzSEODMTo+QVdbM0ce3kFzYyNjmXEWz29TQYA6J5iCA0JUkmoXM0hK+0xFF0LQWK32XE8Q2dGbyblARwUBBOjWq6hjylHMICntoSXhh9auuKMohG69irql1GIGhdpDSbZf1NHCqWsX1lXRhVKvlYoCiHzo1qsQU5gsZpDrSfJr8zxhPlc7hJNsPziS4b7Nz7NuWSc7eofZOzhK33CGc1++nPu37GfPwEjNNZazfe2STi4q4VoVahciH5pRCjFHQip0UG+FBZT8L8qJCg4IUSFCKnRQb4UFlPwvqoEGSiHmSEiLgupt0Y4W4YhqoIFSiDlSaqGDSrYv7mxlUUdLEFqq0a7kf1ENFKMUogyUUuig0u1AMFqq0S5EucgXo9RAKYQQQqDFPEIIIcSs0EAphBBCFEADpRBCCFEADZRCCCFEATRQCiGEEAXQQCmEEEIUQAOlEEIIUQANlEIIIUQBNFAKIYQQBdBAKYQQQhRAA6UQQghRAA2UQgghRAE0UAohhBAF0EAphBBCFEADpRBCCFEADZRCCCFEATRQCiGEEAXQQCmEEEIUQAOlEEIIUQANlEIIIUQBNFAKIYQQBdBAKYQQQhRAA6UQQghRAA2UQgghRAE0UAohhBAF0EAphBBCFKCp2h2a2YeA44FNwGuAa9z9XjNbDdwO7Ix3fdDdPxwfcxLwAWALcDjwEXfPmFkbcB2wHTgqtrWpiu4IIYRIOVUfKIFW4GJ3HzKzc4FPAWfF265x929k72xmBnwbONPdd5rZF4H3ADcDlwBPu/u1ZnZ83Pa66rghhBCiHqj6rVd3v9bdh+K3RwKPZ21+u5ldamafNrNj47a1QLu7T840fwGcE78+B7g3tvsr4EQz68rVr5ldZGYPmNkDPT095XRJCCFEiqnIjNLM7gCW5th0lbvfambLgCuAlwPvjLf1xNsfM7OlwH1m9nKiW639WTb64jYKbOub2rG73wTcBHDyySf7bH0TQghRX1RkoHT3N8+wfSfwITN7A/AvwCnuPgg8Fm/fZWa7gBOBZ4H5WYd3Abvj17sLbBNCCCHmTNVvvZrZpVlvtxDdWsXM/jCOM2JmzcCLgK3AZmAonoVCtADotvj1bcBp8THHA4+6+7TZpBBCCDFbarGYZ1W8IGcP0YzxfXH7s8CVZvYIUezyz919G4CZnQ981sy2AY3AN+NjvgxcZ2ZXxsdcWDUvhBBC1AXmXn/hupNPPtkfeOCBWssQQggREGb2oLufPLVdBQeEEEKIAmigFEIIIQqggVIIIYQogAZKIYQQogAaKIUQQogCaKAUQgghCqCBUgghhCiABkohhBCiABoohRBCiAJooBRCCCEKoIFSCCGEKIAGSiGEEKIAGiiFEEKIAmigFEIIIQqggVIIIYQogAZKIYQQogAaKIUQQogCmLvXWkPVMbMeYFutdczAYmBPrUWUiTT5AunyJ02+gPwJmST4coS7L5naWJcDZRIwswfc/eRa6ygHafIF0uVPmnwB+RMySfZFt16FEEKIAmigFEIIIQqggTJcbqq1gDKSJl8gXf6kyReQPyGTWF8UoxRCCCEKoBmlEEIIUQANlEIIIUQBmmotoN4xs9uAjqymE4AVwDLgdmBn3P6gu384PuYk4APAFuBw4CPunqmW5nyY2dXAGVlNn3X3O+NtlwJdwELg39z91rj9JAL0BcDMrgcOAAPAicAl7r7TzFaTsGuTCzM7E3gnsBtwd/9kjSUVxMxeAnwGeAh4EbDX3T81m89dKJjZfcBw/Hbc3d9oZouAa4DNwFHAx9x9V7x/sP7Efxc/AZ6Jm7qAXwJbSej1OYi766eGP8C7s16vBf46fr0auCDH/gZsAJbF778IXFhrP2ItV+dpfzXwL/HrZuC/gQUh+xLr+UzW68uAG5J6bXJonQc8CbTG738AvLHWumbQ/CrgN7PePw68stTPXa39mKJxmnbgr4Dfjl+/HfhWEvwBDgPOzHr/SeC1Sb4+kz+69Vpj3P17WW8/CNyQ9f7tZnapmX3azI6N29YC7e4+OZv5BXBOFaQWhZl93Mw+YmaXmdm8uPltwL0A7j4GbAROJ3Bf3P3KrLcNRDPLSRJ3baZwGrDN3Ufi9yFrBcDd/8vdf5TV1AAMQsmfu5A4PtZ8tZlNnv9ziHVz6HUJ2h933+vudwGYWStwsrv/PH6f1OsD6NZrVTCzO4ClOTZd5S/cguwCVrn7hnhbT7z9MTNbCtxnZi8nup3Xn2WjL26rCoV8Af4R2Orug2b2f4gG/QtjfRuz9p3U3EMNfYGir80C4E3AefG2IK9NiSRJ6zTM7FzgDnf/tZmV+rkLib9w9/vNrBG4x8z6OfTa9AELzayJZPgzye8B/xC/TvL1ATRQVgV3f3MRu10IfD3rmEHgsfj1LjPbRRQnexaYn3VcF1GMqSoU6QvAvwOXxq93k1tzvvaqMZM/ZtYN3Ai81933xccEeW1KpObnfraY2euB1wOXALj7Y1mbi/ncBYO73x//HjeznxH5Nal7P5Hm5909Y2bB+5PFbwG/Ccm+PpPo1msAmFkD8Gbgtqy2PzSz4+PXzUSLF7YSBfiHzGxZvOtrso+rJWb2hay3RxHFwAB+THSrj/g/42OBewjYFwAzWwx8DbjU3beY2Xlxe+KuTQ7uBY6Ib5FB2FoPEt+efDPwIWCZmZ02i89dEJjZ0WZ2YVbTpPbbiHVz6HUJ2p9J4n9k/jO+nTqb74XgUMGBADCz/wm8yN2/mtX2BuB/AY8ARwI/d/e/i7edBFxM9ASURQSystLMPk+0SGQ3cDzR7clN8bZLiVa2LQT+1Q9d9RqcLwBm9hDRXZd9cVO/u789idcmF2Z2FvAuolvJYx7+qtdXAj8FHoibOoj+kVlHiZ+7EDCzFUT6HyKaTTUDf0a00O0viD5DLwEu90NXvQbpzyRm9g/Axe6+J35f8vdCaGigFEIIIQqgW69CCCFEATRQCiGEEAXQQCmEEEIUQAOlEEIIUQANlEIIIUQBNFAKUWeY2afM7B1F7nuBmf27mb2/BPsdZnalmW2JC2ULkWiUHiJEnWFm5kX+4ZvZBcBqd796Fv3cTVQ8fmupxwoREiphJ0TAmNkVRHV0zyZ6UsaZRI/72pS1z6lE5dweAk4CLieqm3kj0dMbTgFuIiqxdw3wFTN7xN2vjmtvvgTYAxzh7gVnjmb2PaKnp9wJ/A/gh8AS4OXAQ+5+VRncFiIoNFAKETDu/vm4vNf5wChwnrsPTdltELgiLrP3TqKqKJea2flEz838fWCDu38MwMx+SDTYAVwEfMjdf2pm/6MISZcRVce5iqiazHNEhawPEJXx00ApUocGSiHC57NEz168IccgCTAE/ImZ7QHWAC0A7j4R1xLdQlS0PRcXAB+N63F+C/jPIvRscfcJYL+Z7Xb3AQAzmyjBJyESgxbzCBE+5xLdMv2Ama3Nsf0LwKPu/nng+1O2XUg0GH49qwB6Ni92998jemrFh81sUflkC5EONKMUImDM7L3A/yZ6yO0bgVvM7IPu/tOs3b4NXGxma4AXAyeY2auB3yGaYV5PdMv0e2b2ReDtRM84PBZ4h5m9AnDgnyYfJVaA9xE9deT1wBFAd/xsSOLX73P3vy2D60IEg1a9CiHyolWvQujWqxCiMM8BK8zsj4s9IM6jvJzo8V25YqpCJArNKIUQQogCaEYphBBCFEADpRBCCFEADZRCCCFEATRQCiGEEAXQQCmEEEIU4P8HCYu9kc/cxXAAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=[7,7])\n", "plt.scatter(x_m4,y_m4,marker='o',alpha=0.5)\n", "plt.title('Footprint Diagram on M4')\n", "plt.axis('equal')\n", "plt.ylabel('y axis [m]')\n", "plt.xlabel('x axis [m]')\n", "plt.show()" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Easy enough! Now let's plot the OPD v.s. the Entrance Pupil coordinates" ] }, { "cell_type": "code", "execution_count": 23, "metadata": {}, "outputs": [], "source": [ "# Grab the OPD\n", "opd = np.copy(rf.opd[0,-1])\n", "opd -= np.mean(opd)\n", "\n", "# Grab the EP Coordinates\n", "x_ep,y_ep = rf.xData[0,0],rf.yData[0,0]" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "Now it's just another scatterplot. Here the dimensions are in nm" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAb8AAAGvCAYAAADLxdvpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy88F64QAAAACXBIWXMAAAsTAAALEwEAmpwYAAB6dElEQVR4nO2debwkZXnvv09V9zlnzuwwAwMMMOwRRbZRBBRFvaICcb+JV1y5ahLjjSLEaLxeFyQuGHdvxJiExCQaMYkELkIMKIqIDAqyKfvuwLDMfuac013P/aOququr3+6u7q7eTj9fPs2cfvv91fvWW0/12/UuzyOqimEYhmGME96gK2AYhmEY/cY6P8MwDGPssM7PMAzDGDus8zMMwzDGDuv8DMMwjLHDOj/DMAxj7LDOzxhrRGR3EfmIiNwgIptFZJeI3Cci/yAiJznyf0REtMGrJCKPiMg/i8gRKV0jzTYRuSvSvE5E/P6dvWGML9b5GWOLiJwC3An8N+DPgXXACuB04HHgShH5hogUY42qfkRVBbgwSnqrqkqUtjfwv4HTgOtF5AUJnQAHRG/vj957wH7AW4EtwD8BN4rI03pxvoZhVCkMugKGMQhE5LnAxcCVwOmqWkp8fDPwXhG5Cfhbwvvkza2OqaqPAd8Qkd2ATwNfAZ7eJL8CTwE/Bn4sIv8CXAL8RESOV9U7Ojo5wzBaYk9+xtghIhPAPwA+8Aepjq+Cqv4d8EPgTSJyWhtFXBn9e7iI7J5VpKpXAh8AdgP+vo3yDMNoE+v8jHHk9wiHOK9W1ftb5I2HN9/X0xpV+RqwGThORJ7fpzINY+ywzs8YR14W/fvzDHmvi/49UUQWZzz+C6N/b1XVJ9qpmKruAn4SvX1pO1rDMLJjnZ8xjhwa/ftQhrxxniLVBStORGS1iJwJfBiYBf64w/rdHf17cId6wzBaYJ2fMY4sj/6dyZB3Z+LvZY7P/zbetgBsBM4FLgWepao/7LB+21P1NAwjZ2y1pzGObIn+XZQh73Ti782Oz98aLYzJk6XRv1ua5jIMo2Psyc8YR+ItBGsz5I3zzAL39aQ29cTDsrbVwTB6hHV+xjhyWfTvszPkPS7694equrNpzhyIFtWcEL29vNflGca4Yp2fMY58G7gfOElE9m+R903Rv5/qbZUqvItwbvFnqnp1n8o0jLHDOj9j7FDVOeAMoAz8lYg4575F5C3AycCXVfWqXtcrcrf2MeBJqp2uYRg9wDo/YyxR1Z8Av0s4rPlDETlFRJaLyKSIPENE/hL4OvAl4L29qoeIrBCRk0Tk64SrRO8Cnqeqd/aqTMMwQEL3goYxnkTux/4XoTPqg4FJ4LfAj4CvqurPU/k/Avwfx6EuVNW3NCmn0Y22A3iMcMP9d4F/VdVye2dhGEa7WOdnGIZhjB027GkYhmGMHdb5GYZhGGOHdX6GYRjG2GGdn2EYhjF2LFjfnqtWrdJ169YNuhqGYRg94YYbbnhcVVfncaxTTnm2PvFEPq5kb7jhjstVdejDcS3Yzm/dunVs2LBh0NUwDMPoCSLSKhBzZp54YgvX/fxruRyr4J+8KpcD9ZgF2/kZhmEYGVEgCAZdi75inZ9hGMbYo2PX+dmCF8MwDGPssCc/wzAMY+ye/KzzMwzDGHcUGDNXlzbsaRiGYYwd9uRnGIYx9vR3wYuIHAa8HpgBng98hDCc1yeBe4BDgA+q6qNR/nMIgzyvBK5Q1Yu7rYN1foZhGEbfOj8R8YG/BE5X1UBE/h4oAecBP1DVfxGR04HzgTeKyHHAyar6chEpAreJyNWqurmbetiwp2EYhpEnq0RkQ+L1jtTnzwIEeLeIfAA4HXgcOBW4NspzTfQewlib1wKo6jxwO3BSt5W0Jz/DMIxxJ99N7o+r6vomn+8PHA+8XlW3iMg3gTlgD2BblGcrsFJEClH67Qn91iitK6zzMwzDGHv6Oue3Ffi1qsbORH8CvAB4DFgKbCac33tKVUsiEqfHLIvydoUNexqGYRj95Dpg92juD8InwTuASwmfCAFOjN4DXBKnR0+ChwNXd1sJe/IzDMMYd/ro21NVnxSR9wOfF5FNwGrgY8Ai4FMicihwEHB2lP86EblKRM4jXO15VreLXcA6P8MwDANFtH9bHVT134B/SyXPAG9vkP8zedfBhj0NwzCMscOe/AzDMAzz7WkYhmGMGQoE5tvTMAzDMBY09uRnGIYx9oxfMFvr/AzDMMadPm51GBZs2NMwDMMYO+zJzzAMw4A+7vMbBqzzMwzDGHvGb87Phj0NwzCMscOe/AzDMMadMdznZ52fYRjG2GPDnoZhGIax4LEnP8MwDGPsnvys8zMMwxh3FMQ6P8MwDGO8UNDxWvBic36GYRjG2GFPfoZhGIbN+RmGYRhjhjm2NgzDMIyFjz35GYZhjD1qHl4MwzCMMWMMhz277vxE5CDgXOAXwFrgCVX9mIjsBnwSuAc4BPigqj4aac4BlgErgStU9eIo/SjgXcC9wB7A2apaEpEp4Hzg4ehYn1TVO7qtu2EYhjGe5PHktxvwLVX9HoCI3CYilwJvB36gqv8iIqcTdl5vFJHjgJNV9eUiUgRuE5GrgS3AN4EXq+pGEfks8GbgG8B7gAdU9dMickSU9rwc6m4YhmHA2D35db3gRVWvjzu+xDF3AKcC10Zp10TvAU6L01V1HrgdOAk4EFikqhsdmlMTmpuBI0VkWbd1NwzDMKCyyT2P14iQ62pPEXkVcLmq/ppw2HJb9NFWYKWIFFLp8Wd7NEmnxWfJ8t8hIhtEZMOmTZtyOCPDMAxjIZJb5yciJwMnA++Nkh4DlkZ/LwOeUtVSKj3+7LEm6eljpT+roKoXqOp6VV2/evXq7k7IMAxjXIgXvOTxGhFy6fxE5FTgFOBPgDUicjxwKXB8lOXE6D3AJXF69CR4OHA14cKYGRFZ49BcmtAcAdykqlvzqLthGIZBuNUhj9eIkMdqz2OBbwMbgKuAxcBXgA8CnxKRQ4GDgLMBVPU6EblKRM4jXO15lqpujo51BvAJEbkf8IELo2K+AJwvIh8CDgbO7LbehmEYxviSx4KXG1R1iaq+IHo9S1X/TlWfVNW3q+q5qvrWeJtDpPmMqn5QVf8w3uYQpd+oqmeq6sdU9T3RMCmqOqOq74qO9Rbb5jDe/OIXv+DEE59HoTDFihWreP/7P8BXvvJV9tlnf3x/ksMOO5xvfevbvPWtZzI9vYyJiWlOP/2VfO97F3PccSfg+5OsXLmaD33ow3zuc19gr73W4vuTHH74EXznOxfxxje+hUWLljI5Oc2rXvUavve9i1m//jh8f5Ldd9+Tj33s45x//l+y55574/uTHHHEUXz3u//K61//BqamljA5uZjXvvb3+Ld/+3eOPno9vj/J6tVr+MQn/oK/+ItPsXr1Xvj+JEcddSz/+q//zute9/tMTi5mamoJv//7/4PvfvdfeeYzj8b3J9ljj7359KfP5+Mf/wSrVu2J709y7LHP5t///d959atfy+TkNFNTSzjjjDfzne9cxNOf/kx8f5I1a/bhC1/4IjpCCxCMQaJjN+wpC/XmWL9+vW7YsGHQ1TBy5u677+aoo45l+/YdlbRCoYCqUi6XK2kigu/7lEolADzPI0jdmIVCOPAR52mkU9WaTqSRzvO8Sh183ycIgtx0IsL8/HxbuunpaT74wT/jz//8Aw3b0xhdROQGVV2fx7HW/85avf4b787jUHjP/bPc6tVLzLenMVJ85jN/ya5dszVppVKppuMDUNWaTibd8cW6ZJ5GuvQPxEa6ZB3K5XKuumTHl1W3c+dOPvWpTzM3N4dhGLVY52eMFL/85S/rOhCjMaVSGdv2Y2TCFrwYRv/YtGkTP/nJT1m5ciXPe96J+L7PTTf9ijvvvJsjjng6hx12KHNzc/zwh1czMzPD4Ycfzi9+YR1gVnzf54477uK6667nOc95NnvvvTfbtm3jqqt+xMTEBCef/HwmJye56667uemmmznwwHUcffRRg6620XcUdHTm6/LAOj9jYPzFeZ/m3HM/ycREEVVYvHiavfbai9/85g4KhQLz83Mcc8wx3HbbbZUhvtnZWVvEkRERoVAo8spXvA4RYW5ujhe/+GSuvOpHFAvxrS8cffRRXHfdzykWiwRBmcMPfxqXff9iVq5cOdD6G0YvsQUvxkD4r/+6ile+8nXs3LGzba2qopQJd+bGSOq9F73XxOceEKTSyKhLzin2WidRPbvVeQgeIkI7TExMcOqpL+Wi736rLZ3RX3Jd8HLYPnr91/4oj0PhnfyhkVjwYk9+xkD4q7/6ekcdH4RPNEK4wjP9xd46zeupLv4x2b3Oz6yrzePWtcPc3ByXXnoZMzMzLFq0qOPjGCPGCG1TyANb8GIMhKee2tz1MVxf8FnSeqkTkb7rWuXpBNVwiNkwFir25Gf0hAcffJivfOXr3HTTLaxffxR/9Ef/kzvuuJuvf/1Ctm7dxr777sv09CJ27pwZdFUNB2v3Xcv73ve/2bjxUU477RTe8IbXcemlV/Dtb/8r09OLOPPMN/GiFz1/0NU08kIZqZWaeWBzfkbu3HTTLbzgBacyOzvH3Nwck5MTQLhhfGZmFwCLFk2xa2YrQVBudihjQExMLqFcVoIgYNGiKUTCQaKdO8Oh6sWLp3n3u9/Jued+aJDVHGtynfM7dG+9/svvyONQeKd81Ob8jPHkXe96H9u2ba+8n52t32Q9M7ML1SKIB1oGBKQQ/qvzhD9FPcSbAAI0KAGKeAVUPYQyGunEC3UaRDrxCOMkB1EaFR2UK+WJVwSo0UEBEXXoStFS8Ga6uJ6D0BVRleY6vPAlClqK2twP07REuFgmvA7z89UfJfEPliQ7duzk85//Ku94x1vYb7+1TjswjGHGOj8jV4Ig4LrrbsiUN5ybKkSdXvKDiVQ+H/H9xHsIVzIWa/P5k6kS3Dpa6nDoJhx50joP8avnMmw6Te7jEj+VqbZNsuD7Plde+SPe8pY3tK01hpAxG/a0zs/oCNVwSMxPdBLlchkRYWpq0vm0YCwsPM9jxYrllMvlGjtw2YYx5Chjt8ndVnsabbF9+w7+8A/+lBXLD2HJ4gN4wfNfyUUXXcLzT3oFSxYfwG4rD+WAdesqTpy7JY+Vi0ZMvm1ZKpV531kfZfH0OtbucyRf+PzXeM+ffIjdVh7KksUHcOIJp7Jhw425lmkYeWELXoy2OOUlv8fPfnZD02XwnucxX5qJ5uQqqdRuyo7Tan9ting1w3PxQota0hvaQRC0xbHz1Lk31eely3KsLLr6zfHh7d6+rr58wfcnkcTvZ8/z8DyvxvXc4sXT/Pz6yznooHUY+ZLrgpdD9tLrP59PmFTvtE/YghdjYXHrrb/m5z//Rcv9X0EQ4HuTqAYoQeRlxIs2cgcoVDyPhN5agujrNk5z66jJE6eFnk1iHdHRO9Ol8+SvC9G2dNVzyUHnVducmuuQ/VpVdKknySAI6qJnzM3N8aUv/TWf//y5jczFGBZszs8w3Nxxxz1tDWeKeDVPBuEQpl/zlRl6a/Ez6Ujp6tO8DHkGr0sPP2bRhe/z02Vt81bXqhXz8yV+9avb2tIYRj+wzs+o4bHHnuBHP7qWpUsWc/ILT2BycpIbf3krv/nNXUxMTFg0BaMtJiaKPP3ww7jookuZXjTFC1/0XKamJvnVTbdx++13cvDBB3DMsUfY3O6gUcbOvZl1fkaFL3z+r/nYRz9PcaKAIHi+x/777cNdd91XiYReKBQBW8lpZKNcDvjHb36P73znMgRBBA48aB2/+c1d+J5PoAFPf/phfO/iv2HZsqWDru54Y8Oexjjy85//kk+c+0VmZ2dr5vR+tXlrTT5VRfCjqAoAgicFAg03oYcp4dBYMo+ID5VoDFTmlgItp3QeSnwswZNiNB9Vquji8uL5J08KqMblJXXlmvLSurA8qRzbrfMj3XwHOj9VXiFql+q5iBScunJlo392nYhP0ELnSbHuXMLrUL1+4CNCYsGSRG1cTumk5vq5dKjP3Nw8c3PVSPQ3/vIWktx0422c/b6Pc8HXP41h9Avr/AwALvzb72Tamyci+H6R9EZxzzkX5No4ndYlTVCiRRQuXe3mbS+T6dbncevSG8NdOtf5ZdF1Wl5nOj+DLnkulR8Qma5ffR53PdvbMD83N8dF37mEr13wKRv+HBTa/yjsIrIIuA64QlXPFpHdgE8C9wCHAB9U1UejvOcAy4CVUf6Luy3fOj8DgG3btluQWGNglEpl2xg/aPo/53cu8MvE+/OAH6jqv4jI6cD5wBtF5DjgZFV9uYR+C28TkatVdXM3hdsm9zHmmh9v4M1nnMXpL3sbS5YsZXraYrcZ/UdEOPzwQ3jbm8/hd19+Jl//2j9btI8Fjoi8EbgGuDeRfCpwbfT3NdF7gNPidFWdB24HTuq2DvbkN6Z8+YsXct65X2FmZzjUOXn9BLO75luoek04X2f0E9dG/P6iqtxz18P85rYHALhhw8383d9cxA9++I8sWjQ10LqNFfmN/KwSkaSHkQtU9YL4jYgcDjxNVT8oIs9M5NsD2Bb9vRVYKSKFKP32RL6tUVpXWOc3hmx+aivnfvRLNdEWZnfN4TGBeD7loBSu9vQKCB5lnUc1QMTDlyJKQBCUUBTfK4B6IEFCV0SQTDolINBaXaClis6TAooSRNELPM8HFSAgiDzB+FFUh3pdUAmZlL8OynGkCfHCRTCidTqlXBlOrtUl9s011IWb1dM6T6rRGZI6was7F9CmOj/S1V6ruLx5FMWT0A6QZLu4bKMQrZiv1zW0qWgxzfxc1RvQzp27uO/eB/nWP/0Hbz3zde2at9EJ+cbze7yFh5dXAbtE5M+A5wITIvIe4DFgKbCZcH7vKVUtiUicHrMsytsV1vmNIddffxOTkxOpUEMSfRkX8Pxasyikoyzg4yXnZqJoCZ3q0os0/FSEgXDbxUSdLj1mX69rUM+cdAW/PoKCS5emZzrqzwUkky59rcLy6qNItL7G4LWhq3qgqWXnzl1c+h9XWue3AFHVT8R/i8gUsERVPy8ivwMcDzwInAhcGmW7BPg/Uf4CcDhwdbf1sM5vAbNr1yy+71MsRkveVdm5cxfLli2l7JzctmFHYzgQEVav3o1du2YpFPyKZ6HYhqenp2xlaK4MZLXnawjn7iZE5PXAB4FPicihwEHA2QCqep2IXCUi5xGu9jyr28UuYJ3fguS2W+/iff/rk9z4y9vxPI+XnXoSz37OM/n8Zy/kqSe3smTpNMU6N2VxsFPDGDzFYoEbf3knB+z9Qjzf53df+UKOeOahfPEv/54tW7azfPkS/vSDb+dtb3/NoKu6cOjzak9V/S7w3VTy2xvk/Uze5VtUhwXGpk1Pcvwxv8fWrdVI6p4vaEDNVgYRDSMvkEzzqY964FX2gsUpIc08/rt0GaM61Gygj98HKZ0PlGtUoTPm5rr0sWsdVDfKk01XX6dG9exQh9S0Z7Z6etGVaa7LFtXBda3ysg0BDWps0ZMJfK86hOv7XhQnsJpn0fQUn/jkezjjza9gHMk1qsOBe+r1H3t9HofCe+MXLKqD0X/+8e//g9m5uZq0oFz/A0dV8L1FDu/+cZSF2ugFikbuqaqRAsK0PHU+nuc5IwzEHVnsWUQjbzGhrnFkgl7oJIqC3liXbs+0rhy1gV/TBrW62mM31pWJZkVbRGdw61Ctuw6trxWJtmtH19w20m2QpFyufyqZ2bmLz/zFN8a288uVfBe8jATW+S0wfn37PczummudkUbe/Rt5/K/NA1KX1r1OEjqvQZ5kWiGVNiq69iNbNNa1rpM7OkPs+qzsrEP2a9yprt422o0YAbBx4+Nta4wGWOdnjBqzs3P86L828NRTWznggH1YND1V2b9nGAuZ/dftzf/7jx+zfftOnnvSMey9z+pBV8kYEazzG3Fuvflufu+VZzM/VyIIAkqlMqXyqIYdstWm/ccVEX40EBEe27iV9/7Rpwk0oFwK+OP3vp6z3v+mQVdt9BiAb89BY8v7RpggCHjT732Ap57cyvbtO9m5cxdzc/NoqVA7FEUBjwmqHYvgUaxxSix4+EzUDJl5FCkwWUkL80ziJxwXe/gUmapxjOxXdFKj85rqBJ9CdGxJpBVrjh3mStezEJ1LXF5+Oq9OF+bKovNb6sI2d+uS16pTXXx+NNT5UqDAVE09PSbwE/YStt1UqrxijS48u8ma6+c5dJ3boks3wfxsObT9HbuYnZ3jq1/8Nj+/9maM9tFAc3mNCvbkN8Lc+IvfsG3bzrp0EY+CpDcauz3++02iJUj0n+/w0l+g9vgu7/4uHZl0rc0yexSJwer8DqMl9FUn7jYvpGzDbT8u3WTDPEH0lNmuLTYrL8mumVn++ZuX8ezjj2iazzCs8xthZmfnbKOvYSRQVXbsMKfYHbFAt701wjq/EePhhx7j61/+N27c8GsOOHgf5xJwwxhXFk1PMT21iNeccjbLVy7lzW8/jee/6NhBV2v4sa0OxjBz128e5FUveR+zu+aYny9xy01349jClxvxPqzhY/CRCIzOCDfr9+7azc3M8/3vXcvsbBih5Gc/+RV/9N7/zh+/7/d6VqYxmljnN0Kc9+FvsGP7TMVTS/jU5zPBNCXmAcUj8qRPQJkS8XIFwaPMfGUDtE+RcKtzrU4JCFK6gFJFF87DaHRsTSzECAiijdN+lNa5br6ySdqlCxfzJOsZLrIIzyUur5g4l7g8L4yWUCkvnAuL2yBekBLqwi/PznVhnTrRVdu8sS5sA22qi9suTNNo3tKrpLmvlUR1T7Z5frpCpGvfFqNIIU11PqJS6fgAZnbO8uXzv8Ub3voyVu62LNuNNq7Yk58xrPz82lud0dZFPIo4POmnLm96AUP4hVWvq4uyULdgQuqO5Sqvc106T1InxN5ZWp+f71io4pFeAJJeoDEMurDtmuvCo7t14dOVRkd3RHVw6mprVX+t8tN1aovNdPVu20KKE0V+ueE3vPAlz3J+bmBbHYzhZPvWnczsnGXpsulBV8UwRg4NApYtX8zmp7bbHHkzAs3nNSLYk98Q85vbHuADf/xX3HHrAyCw9/4reaK4hfn56ib24ZyT6yW2ET4LvZ5bG07cG/YLxSLv/O+fYW6uxKLpSd79Z6/ljW9/af+rZwwV1vkNKU89sZUzTv0o27ZW9/E9ePfjBEH9F3/6i67ea344+FUbLSGci2lX5+FFqubREjy8aC4v1oVRFvLR1UcF8JDKHrJYF7SM6lB/LFc0iPBY7UdZSOtix85aV89WdXLpvJrzdem8yrxgMq0QzaM1q2fkSDsHXS9tUSJbjK9V2E4FlFJNnl1bFTR09zc/V+IvP/otFi+e4tX/4wUYVUZpg3oeWOc3pPzrP1/N3Nx8TZoqFGQKXycqk/5Vj/+xl/xaT/qVreoNPPDHXx5ZdfHcTY0H/qjzrU9LeO6vdJp568IOM65TpQ0qeRLn15au0/LcuvhcutHFnflw6WqvVagrR+dJD2zRFQ3CpatGu0BrfzDOzMzy5U9/1zq/JLbVwRgW7r3zEWZ3zTs/c3v8b+0l3+WBPz1s2lQX3Rux35eaPM60VJ0GopMOdcN3fl76mg+FrvZcQl2YJ37izNcWXXbu0jX/atv06OamnxsLH+v8hojZXfP86PKbeOqJbeyz72qmpyfZuXN20NUyjAXH/geu4dLvXMvMjllOeNEz2HvfVYOu0uCxJz9jEPz65gd4+ys/Q2m+TLlcJgiUeXUv3R4YtrfcaINhXXQjIjxyz5N8/H3/QFAOCD6gnPnel/OHfzrGQXHHcKuDdX5DgKryv/7HF9ny1I7UJ4Vo2j/sBL1oc3O8CTzebKyVPOG8R4EJypRqNiR7eJQqG4QlsUG4uiHZp0iZ+ZpN57W66r6rePGD4OPjU65s5ibapCyVY0tUq7Cepcq5ePiUo03SREfPVxfUtF36/Go3XMdtVySgHB2LqA0KKV2xUqdWuhJzifnSatu1ows3tUM5cY3D8ys7bKO2XcLyqjqfAkELncum6nWxA4TGtujh42uBUke2WIhsuJEt1ttGvS26bcPXIqVZpTRbjXn5t1+4jBNe8HSOfPbBWW5ZYwFgnd8Q8OubH2Dr5nTHF/5ydm0Ybua5P57/8JkkHUHB5RG/WOe53xUpwGUmtZubXVEPhkMXIpVupF5XSB3bp1DnNMDd5vVRD+p1i7rSJVfW1m8o7yyqg7vteqebwKczW2wc1aF2xXF3tgHhlMP3vvXT8e78zLG10W9K8+VcozOkFwAYhtEcVWVudm7Q1RgYCuiY7f+3zm9APPLAE1z4+Sv41c/vYd8DVue2b3t4nVEPFmsXoxlTiyYI5uANJ53Hit2X8Po/OJnnnmIxARcy1vkNgPvufJS3vOhT7JqZo1wKuPu2R1AvGnJwOTCJRyOafXcPQpfO06muUT0Wqq6NNo87bEVHyzZ6bFM1i2lysMX5XWV+ePFNzM+Fc5q/+vk9vO19L+WtZ42JJxjb52f0gy9/5N/YuX224qRaVaEsTDAZLuTQcI+Vjx8tBAhvyILGC0nChQfhlH+hslhANZyDkWhZZjnaZ+VrmBYuWNDE3Fd47FDnQUInkY6oPHXpAF89Yq8vZcrRgoKwvFJlYYdU6lmKF3ZEeZrpkm2Qpy6IF1poODcUn59H7bnUtl2Qr07b0EVzWCVKqCo+fkUXn1+tbQSVNk/rkm3XTBfnaUcXz7UlbbhtW4xsKostTlTuj1KkS9sG+FE9SzURRtJt7iHqVTo+gF075/jGZy7jNW87iWUrxsSnrnV+Rq/55U/vckdnQCqLL5IbtNOLCsIvtUKNzsvkgb9Qt0W5QL1v83oP/GkzqeoEGtazmFqIEH5htV7A029d+vyEqheXpGbQuiznV3+t6nXx8dutUzNdcvlJegER9M4WG9l+o7aLF8o0Or8azYTP7b+8n+NOflrTfMZoYp3fAFi+22K2bt7ZNE/8VWLzVIYxGEqlgN1WLx10NfrGuC14sZBGA+CMd7+YYtG1RLyWUej4xmugxBgn9t5vdw5++j6DrkZ/iOf8xiikkXV+A2DliqV4gRetLw5f4TCM1KZp7XugPk88PJpI89RrqfOcuvryfLyWOnc9XToy6OrrSQc6cehc51dXnrbWSRc6Z5sndeSpy3aN665nRptK5hENQ/Kmdb5D52ewxey2n8U2vJrPC/G8eJN6euKxcsnSXLcgGcOFDXsOgL/71BV4JZ+paOJdkPCmJeEMOPqvuplXEKUmLamrScug81SIF1a4y4u+8NLHzlBPz1UnlcSR6+vZ8Fy61qmznk3bTgeo0z7o6FbX3BYD1HmtWpUX21RUQiabal5eY12hiW0AeIHH3Tf/lgfufIz9DtnDcRcvQMZs2NM6vwGw6ZEtQO2ilhi3J31XZILmaVl1nZbXia76rr16dqfrXduZrpFthM+q7evikEjapq798rLkKUz4bHzwqfHo/FQtnp/RGx594Cl+dvntFIo+BzxtDTf/7N5BV8kwjCbMzZZ4auM2Lvry1Rx4xN4c9bwD8TybKVooWOfXB779hR/y9+f9JwiICPOlEiKMmys9w8hEuMNv8BQ9n6/86cXMz5UoFn32OXg1n73knUwvnWwtHkXGbNjTfsb0mHtv28g//MUPmJstMberxOzMPMG8UlC/MsQihAtL/OryDjyEYk1KmKdIUhdO3hdSugkK0ab1xrp4ozQJXSFVXrxZvFbnZdR5LXUTkaf+uA0KjnoWo5omz2Uiccahzq8uYuiBrpjSFSu62muVRVfsQFdw6lzXuFBjUwVnPVvZRtUWY2I7a9c2OtVVr1VzGy6k2s5roCukbLGQOr/0vRa3QWmmzMz2WUpzZWZ2zHH/7Y/yNx//PgsWzek1ItiTX4+56qIbmZ8v1aWHN6rrt0eLaNrSKDJB/fHrYwC4tlc0Ly/M0X4963PUpoazNOJsg3Q9fdJ++3Hq0tum89S52qBTXZbzy6ZzXeMs9WxsG9UlH/Xl5WsbrXVZbAMaX+OkIwmXLabbpdWm9/m5Eld+55f88adf0TSfMRpY59djyqVg7CaSDWOhEpQX6L2sjN33lHV+PWDDf/6G737xxzz56Fb2f8YaipMF5nbVP/2NM+FCdMPoHSLidCPYKX7B48DD1/CeF36FuV0lTnrNM3nlH5zI1OJ692ojyZjN+VnnlzPf+6truPDjVzC7cx6AR+59knKp3EJlGMOBRMtNxusZINsim6AUcM+NjzA/G97P3z7/Kn78bzfz+Sv/iOKEfZWOGnbFcmTXzjn+PtHxQXjDhFPpGm3iFXwET4QAJVCN5vHC56BS5GDPE4lUUI7SfAmn5AOUAMVDKpP05Wi2JpxLCZ+rypWoANXNwO3oCtG+p3JC58d1inZjxeemFV2cRiZdJ/WMv6ha6TRKy67zKq2Q1gWQiLLQvi55Li5dXKd4aYZLF7dBb3TJtgvnvkqVa1Vtz1KLa5zWhW1etY36NuitDfviVfOoNryvUEUFguhJMVnPeHO8h1Q6PoC5XSU23vck11x8Ky947ZGMOuPm2zOXzk9E1gDnAkeq6rOitN2ATwL3AIcAH1TVR6PPzgGWASuBK1T14ij9KOBdwL3AHsDZqloSkSngfODh6FifVNU78qh7njz4m8fwfPekeXjTJRZ7RB2enxr7K0ra1z14kl6wIHUT+IWaQcTw73pP+u3p4kUptXmiOtUcSaLjpz3wZ9O1X0/Xsep1yS/bbnU+rgUneerqF8G0357d6OrbvEi9LU60OLZLR6JO8Qb2ftlwTR6pzVVzX0VuzHwh8pxTLa8Zu3bMceOP7hr9zi/85TZW5PXk91zge8BRibTzgB+o6r+IyOmEndcbReQ44GRVfbmIFIHbRORqYAvwTeDFqrpRRD4LvBn4BvAe4AFV/bSIHBGlPS+nuufGitVLKM3bEKdhjAvFCZ/V+ywfdDWMDshln5+qXgRsSyWfClwb/X1N9B7gtDhdVeeB24GTgAOBRaq60aE5NaG5GThSRJblUfc8Wb12BYcduy/SolVHZT5lVOppGHnSjjNr8YSXnLG+h7XpD0o47JnHa1To5Sb3Pah2iFuBlSJSSKXHn+3RJD19rPRnFUTkHSKyQUQ2bNq0KZeTaJfpRRPOoa8kQv1KR1ea6+J0qkunxbGym+Xppp6uY6VxpXWqc9Uzjes6eKmMrnbJqkvjOpe0zrX/LOt1SJPnteqVLYojLYstNjpWqzxZda7ys9hG0fMoFFqHJxt64mHPPF4jQi8XvDwGLAU2E87vPRXN38XpMcuivI3SafFZBVW9ALgAYP369X1/cHnyt1u57af3URSPghf7lQ9/SRYS8whxWjy3EL+PzqEmLUuedJovXs28hZdRlyyvU1079ey1Ll7mrkRfZl71OmRtg/T1a1eXpZ7NbCPM2Pz8Ytptl0503V8rr20d1NtiGJ+oh/X0GtczXhQTduTClf98A6957wswRotedn6XAscDDwInRu8BLgH+D0D0JHg4EM/5zYjImmjoM6mJj/XjaM7vJlXd2sO6d8STv91KYcJnfraESK3P+PT7rGnt6qRDXb/r2Q9d/GXVbrsM07kkPUY1Or9OyutUNyrt0st6eon2m58tsfHep1gIjJuv4bxWez4feCOwl4h8CPgs8EHgUyJyKHAQcDaAql4nIleJyHmEqz3PUtXN0XHOAD4hIvcTjghdGBXxBeD86NgHA2fmUe88CIKAW398L/fc+DDLVi2hPOAFL7Z53DD6x9TiIqv2XsbFX/wxS3ef5rjTn870sqlBV6sjRmm+Lg8kTw8Iw8T69et1w4YNPS1jduccn3j13/HQbzYxt2ueiakis/Mldg3Ym0t6ztEISc/hjApj5nUqMwGDb5hC0WPRRIHyfEBhwkc84c++9UYOffZ+PS9bRG5Q1VxW2xyzx+7649e+PI9DseT/fjO3evUSi+rQBd/7wo+5/9ZH2bVjjqCs7Noxh84FFBJDJR6h491kQ/uEj9xxmkTv07oJqS6IEMI9SEVJ6aReV5Da8jwJtZI6VnKPYTs6r4UurmdSV5QwX5zmO86v0EDnO3ReTR6p0XkNdMk2D+sZ7hmrOb9UHi/SJfMURGraPL5+vdJNeELRk/rzqzkXt20U07aY0abSuvS18lLXs1ub6kTXqU0VoryN7tHYDgqpPC6dV1LmZkqUSwGzO+fZtX2Oz731WwTlEXuMsgUvRjv8+F9uYn62/inPF8FPzau09mFf1SYpei6P//WkNxtnLi/1PqsuTVInQOiVo5766BOtz891nHpvivXlZVlJibjr2VJH7dxPwzw56vyURwR3NAiHjnx0yWsVP3VlsZfe2WJYv3ZsqtlAV6e2n2Z2Zp57f/VbDjp6nw7Ug2Pchj3tya8LFuqQsWEY3WHfDcOPdX5d8NzXPpPCRO/2+LSx19YwjCGhOFngwCP3HnQ12kY1n9eoYMOeXXDcaYdzyZevGXQ1hg5bcbrwEUbTA5BI77+gj33JoQ19/A4tCgT9uWtF5CBCX9C/ANYCT6jqx6QDf9DdYJ1fF/zHF39MQZTAC1fkxQsBoLpCT4SKF/k4LZ66CQhvxLROo/fxF0xlo2/0rystHq4XmuviKBIunS+466lUPelHaXFMz2Q9XbrkuaTbpVLPxLG61SXz5KGL84ySrpktNmqDLOXV5gkzBYQ21bZtRLowbkg1DZrbcHjPNLbhVjpPpDKi0qie9fejtmVT1198K2/42EtZvHwRhpPdgG+p6vcAROQ2EbkUeDtt+IOOt8h1inV+XXDX9Q9CdNPXua1KvI9vcKdrqyY6x8cN01wuuFy69GKL5Mq4+GZ2nUt68UPBUc9mbdAoTRyVNV13ulbXwXXsLOU581C/SCuTbUQ61WpHlcWGw3vGbcPt6rLbcBi0qzYtVV4irTDh89Dtj3HYc/Z31GQ4UXJd8LJKRJL7zC6IvG+FZalen8rvATsIfTh/Ikq7huo+7xp/0CIS+4Pu6unPOr82UVUef2Azni8s32MJWx/f0VoTxQMzDGPhU5orM7l4gt/euYnd165gYlF67ekwIqjm9h31eNZ9fiLyKuByVf21iDTzB317Qub07dwu1vm1wf2/eoQL/uAiNj8aXp9FKxbhFzzKpeY/mUah47N5OsPIAYFFSyb4zCv+Bq8Q+th9+bufx8ve/VynW7lxRkROBk4mDFkH7fuD7ooRm5UdHDs2z/DZ113IpvufYn5XifldJbY+ug1Sm1nd3uDrZ9g79aSfySN+huGxrJ7003mc9XSU16mu03Zx7RFL63oZ1SG5+TnO06nOdY1bRYNw6eKh7PSx0uRli5l1ySkBqd28HuMa4s10rTLYVCdRJOJ9qy1tQ2HXUzPMz5aY3THH3M55/t8Xf8y1F93kKGGI0GzhivIKaSQipwKnAH8CrBGR46n6cIZ6f9DHR7qkP+iusCe/jFz/77dQTnV0olDMHCkgoYvuoHi+Q5qkdaLzPLdn+6yRCarlVfO0Shukzou8kWQ5vyx5+q3rNApB8tg6IteqlS5r2yWCOtTdD36G+ypOq9Yhy/3YvJ6VCByp37pzM/Nc9sWfcMLrjmKY6dcmdxE5Fvg2sAG4ClgMfIUO/EF3g3V+GXlq41bmZ9w+O7N5jHfpHE9IqbTOdZ17tq8vr3XasOiynF+n7dIrXR7H1kRammG9Vq60rO2SXIDiumda3R9xWpY8Wa9DcpVpmixrA8YFVb0BWNLg47c30Hwm73pY55eRg47dl8nFE8zumBt0VQzDGDEOGHJXZ+F2Dccv7QWMdX4ZecYLD2Z62WRN52eLRIYDuw5UnjzGhfRw6jCQnqaopgunveek/leoHRS0T5vchwVb8JKRR+96nPktM5XJeSGcXK/33F/raT72Bt/a43+9LhmZoFKe1Hvgry9PnVEkkuX51C40cNXTdX6xLn0u9VEBqmnxubTywF8UnNEZCsk2aKBrGX3C1XZe+Epfh37pvCY6aVPnis5Q9MJXXZsn2jN9rSrXoYXOZcOZbKpTndQuDmp2z7juj3TEhuRCFZcNu8prFPEjqZvwwleyTlMTHjdeehvGcGFPfhm58bLbCUoBRccSxVYr2hrmyaBLb8SFzj3iO1dJtqqntK8Twl/BeUV16EaXxtkufoY8I6Jz2kYPdBXPMi5dh7afxRZdUR3aiQaRfDJrqWtgw2lc5U0mhFoK+MXFt/Daj740w9EGx5A9SPcc6/wyEk7QC+M1uGQYRi6MwB4/m/MzKvzqstu58qs/Yeum7ez1tD3xfCEo96/8EbhfDMNogV/02PPA3fjkC75EuRRw1GlP50Xveh5TSycHXbWxxjq/Bvzw6z/l8s/9iPmZeQC2PrZt7II9dootQFn4jOoYSKNFKb0kmA945FcPU54Pv0Cu/tvruPW/7uCsS9/Z05Bo7TJuC16s83Mwv2ueKxIdHwAaLaLw2vfA3yxaQmVjepSWydN8VJ4k5uOSkRA8iSIxtKnzE7qaejp0yXp6DXTxvqd4wDg+FolzjsuH6rGb6byoTV26tHf/tC55bNN1qevSNmJdck45DxtuR5epngk7q7RBm/e/QKXjAyjPldn8yFZuvvx2jj79GQwDoxaLLw+s83Pw+H1PIi7fW1Q7s5q0yv+qxHnixR/JGzGdJ0lWj/iNyotxRprIUZeup0vnasG6RQaOTJ3qsiwEynJs02XXdWobsS75fTuMNlz5IZfoBDuNwJFkbucc9214cGg6P8jVsfVIYFsdHCxdvYTyfD6Te8lfhYZhGADFqQK77bdy0NUYa6zzc7Bk98Uc9vyDGj79tYstXDEMN2N7a4iw/lXPHHQtaggCyeU1Kljn1wCdKzu9wadxeuBPpaXNQRrkSefr1JO+a7iq3577XefcT13ac7+QLRJC5ugM/dalh+iop5NrldUWO7WpLLr0teqmvF7ZVLxvtVU9Xcd2pRVQgvIQraDT6rxft69Rweb8HGx7fDv3XXd/6NEhYbVpj/SeSLR6LLtXftdxmqV14knfE8EfYs/9/dAlr00yXzvRGZLHHiZdJXqAQ9fLNvf6rGtmw9qD82uZx3PrCnGdHNfFdazYEfYv/vUmnv/OEzEGg3V+DnY8sRO/6FOaK2f2Wl/3a69BWrP3neYxXXu6Tq/VsOiCVFqz43Sax3Sd6bRBnnRaabbMlo3b6vIMinBtQn2dFzLW+TlYtW63QVfBMBoijOYeO6PKxPQEBz5n/0FXo4Zx6/xszs9BYbLAc844tmU++wIyjMEwjF/T8ZamLExMF3jaiw7raX2M5tiTnwNV5fZLb6HoKaVojMkj3LtTpnYjrEc4DBUQGn68v6ccb7ylPk88wd9KV3CUJ4k8QujtX6GunkELnSfhH/GcuzTQQevN+JXza6KL00oJnR8dJx7GK3iA1tYzdhBQbqGTqA3S5QU0vg7NrpV2qEtuHu+VLrY7ly2WMthGuzbVyDbSOj/SdWtT7eqa2lSq7SKzB0lsjs/Q5kJ4bG2hK0SPE2lbTNswM7P89tZHWHvkWoaFYMye/Kzzc/DkfU8yszkMX5TFI34vojrE8qLDHjv13O/0dt+px/8sx069F0daVk/6rqgOaVyeojJFu+jyWg2TbsKR1m1Uh3iEYxhtql1ds83qWdo8q85li0kbLs2WuP37tw9P56cydu7NbNjTgfjDFyjTMIyFg4jgFezrd5DYk1/Ept88yo8+dyWP3PQwi1cvZtGyKbbt2j6w+ijuX52jygI6FcPoGq/g8dgtj/DF489neuU0zz7zeI549VHOVaL9IFztOZCiB4Z1fsCmOx/jH17/N6Eja4WZp3aiPe55bMWeMSrEd8JCsFcRkMSc4MCYL/PgdfehgTLz1E5+cO732fzAU5z03hcOrEo25zeG/ORLP6p0fDESKBOSWjAh4a8jl0f85GIWaK2TaBI80PDvRh7qG5VXM3lPfbSEdnWVOiV0fpTWjkf85EKERh74iY7bShdHtkjqPInqNAK6uM270rWyjT7qskb8aKZLdqTt2GKnNuzUpRbYuOrpuo8blRdo8/u/YT3jwoD5mXmu/7ufcdz/PIHJpVMYvcc6P+C3v3rY+bNWHIsvnGnUT5620sXLoltFemhUXpqs9WxHJ8n3KXE7nvtr0sZM57zGnepcttFHXR5tl+4E646Vrmc6LUddIxtODv9lWmAjre//Rmk1xyn6PHHP4+w9oEUw47bPzzo/YNk+y9m2ceugq2EYxhhTni+zdM3ygZQdPvmPV+dny42AE/7webbyysFCmOMxjHYZxJoTEdj/uHUs3XNp/wsfU+wbH5jfPlu3Nyk5pxCTxSO+K6pDVs/96fJcuqye9F1pLeuZHgZqo7w0oxpFwhVhIGudOrlW3ejSeTq5Vv3WeQ3SWtlio3utk3q62txZT8cQb921cnSUmWzKMRRd2j5Tf7B+oeGwZx6vUcGGPYEN3/gpUgooJu6S+Ndfcvw/nebKEy9iyKLzm+TJWl5WnQgUMuiy5HGleX08F9N1p/MHqYvusXZtMSZTeR22XY0uQz0152v1+K8fZcvDm1m+zwoGwcBXwPYZ6/yAmSd3AO7hjixpC0k3jHUy3Xjqku/THceg6pRM0wx52inPK/rMPLFjYJ3fuGHDnsC+xx+ApMc4DMMw+oiWA1YduseASs9nyHOUhj2t8wPWv+WE3FZ32CIRwxgP8v6af8arj6YwVcz5qNmIV3vm8RoVbNgTuOPSXzFZhLlZbR5lgXBcvBzp/Giuq0zCq7tDl4wUANG8oLYXYaCRLulFPs7jU+u5vyBh3jJNPM0nyqtsco/KS+o8wj9a6qj1iN+pLvakn9RVzo/a8yul2g5JePzPqPOBIFGeH6Unr7HfoS62jTod2WyqXV2jyAStrnFdBI5G18qh68SmsugaRfxAaqNWCLX3TKP7Ko7OAFX7zHI/JnVxWjObkjbKu+fymznp/S/B8+2ZpB9Y5wfc9f1bCOYDXLsdkquyhLDBahpNGqxaTI/npzPkpBOiSf5Oozp0WM9eesRvpat0pikmOiwvrfNwbIDutS6Z0KDN29ElRyA6tUVXeVl0PbMpaR59Iu2dpeZYPbLFRrq47ZT6jriRrrSrxBN3PMrqp+3lOGLvGaUhyzywzg/wCq18LxiGYfQWVcUrDu67KNDWeRYS9nwNHP6aY/BdAbhSDOXvoqGslGH0n1G/FaZ3X8JuB60edDXGBnvyA9Y+e3+uKyVnCqB6K6V/DkkizZWnt7rwXSKPRpMM7eqG9Pyy5dEudOnPW9Uzq67T8vLUVfOEn4zKNc5PVz3vftpi+vOqTqL/a82x3Lp91+/HwEIaqQ17jiU3XnA1BQlCp7ZRmhfdQvXRGTThnV0RqXqxF2rNvDLh7dCl86R1HuB7oa6s8XGqeaoe+Ot10kZ59R7xu9f5jrZL6pRExAitLjLyUnmqnjE0odOKLnmt6utZ1bmuVVKXjIhRf37ZdK1sox+6ZBukdWkbLqfaXKHpdUjrcLZ572wqk04SNpW6ZyrRJxJ1b7+8MCV9reIoEo3u/7TtN7rG915+C8efcwqTyxYxCIKRf3ZuD+v8gMdv+y1o/eS2kM1jfBYXRy6dy9Qk/a/UL6KIv/zyKK9n0SAceVy6Tj3iZ2nzrLpOy+vUNnqpy3qt0jd+2Oll1yU74l7WsxNds+gTmuhwBlXPSoeXyuQVfbbc/wR7HDGYqA7jhs35AUvXrhx0FepI/oI3DGPhE8yXWbLX8oGVHw59dv8aFazzA45+5/OR1N6a8RoAyJcRsn/DGA4E9nnOgUyvGkxUByWfDe62yX3EmHtyO0Uf5srVr+146C7p7DXeIJycto43ubar84FSKo9SO5xUcOjS5UV7f2s6nALVTcRxHlJ1SuviY9eXp5QTPwW8aCYjWU8/ypMsz0Nr5hAk0iXr6aeOjUPnoWhKV3/s8MjNdK488V/pY6d16TbIquu0vCw6P8rTG13ra+VHR25mGwAFh21AZ7aYvq860cWOr5Mb6IVwaiHemN6sPFK6LPd/up7x5vja+18oP7kdo39Y5wfcduE1eKUyk4nv4Y69wWfQxRPz7UZQaKe8dj3wu50IhwtXCglhmE/rdAVHmmbQeXV5sunq83Sjo6XOXc/Wuk7Ly6Lze6rTTLpmthGnJW0j7CulA1t0n0tHOs+dJ76vNF6pkrG8XKJIiLL5rkfZ9tCTLF27G4PAFryMIbNbwjhaw+jZfqHrhrFOpuudTjPkybO8TvO0U8+86uQVfGa3zLB0QOtdRmm+Lg9szg/Y+/iDLKqDYRiAu+PqB6rKyoMHFdVh/LDOD/idNzynr6s0xuwHlmGMHIPo/w551dH4kxbVoV/YsCfwwKU3VqI6xMP98eLPchBtTJXqxHXFc7+EeWs8zXvhPFkpkEinVQ/8kWEUJN4ELhVdvBk3ziNoYvJeovLCgvLQ+VHo9XKDesZ5PEKv+cl2USBItUsAlY3F6baTqK1UwjSGXee4xp3qam0jH112W3Tr4nMhyiPaP13BA1GlrJKITBLedxVbJOkkoWrD4fklbFgiXQvbT+s8wgvZTFeI7o/KfRyXJ2EahPeMFy3wacc2GtniQ5fdxDHvfWndyvN+UfUINB5Y5wc89INb0VKAy6es50ir21Rb+Su0bgEmkrPZUZ4C6bT6Z8C4M6s9dnJmPB+dEA7v+Bnq6WqXTFEkRlmXzrOAdJkjjPRUV2/DsS3WRmeozVdI2XnyB2Btec11ZNSl7+NkPeNP8rpWwXyZzXduZOXv7O04opE31vkB3oQ1g2EYg0UDHdiwJ2pRHcaSg17zLLwMUR0WFOM1wmGMEONqmov2WMbSdasGUrbN+Y0pux+5L5SS28JjpEWaILhjJbTSuT9vrQvfuUpUp6JRPcMJCVc9GpXdLD2e3WhX12l5rs9x5KnVudpl+HSOa9XgaK1tqll6q88HpwvNUlNzUK52aWT76XapzdP4OjQ/dtZ7zX3s1u2y+oh9GFRUh3HEOj/gjr++ioIEBFL9Gg9NMLzZ4glvEeqiF3gSDhdUPeKHtOO5PyxPa3SNPM231imeF+7ZCaJMNZ77Ndpk3/D8Wnmob94u7eualxe3Xbs6P/oCTbdn7bWqtl187Ppr1W+dOnTh2cW6eL62tW00b/NOdXGbd6rLZBtRRJNGNlw9lrS8Z+rLa3Sv0VTniabygIhWrhUNzq+RbaTb5bf/dQtzZ72cieXT9B/BFryMIVvucEd1gPohmNDTRe1nWSITJHUxPff4n8oX34jptFbluW6JfupEOtd1EkVi2HU1rrocujRZ26Wfusy20cCGkxuy+xmdodm1iju8rLp0mjfhs/3BJ9htIJ1ff+f8ROTFwKuBxwBV1Y/2r/QQ6/yAJfvtzq7Htg66GrkR/+I0FiZC8wE0YzQJ5spM77Vy0NXoOSIyDfwV8HRVnRWR74rIi1T1v/pZD1vwAhz2P1/Y0d6aYf0Cso7PWOiEU2PDZeld1UZgzxMOYWr3JXlVp200Gvrs9pWB44H7VXU2en8NcGrPTqwB9uQHzG58iqKvNVEdsnqozxq9oFNP+rU6jbzIV3XVeSSpHLtAUFdPobPoDFl0Xobzc+kKopS0eXmF6Djt6nxXPUUJtHk9fcKFFjXtEpVHC126Ddy6VMQPR52yRNIoRsdOziN5KCVHG5QTuoJUN5g30uVpG3nqqGvzILpWVdv3JaCsXtNrnDVyRye6LNEuXPeMB8w/+hSDIp5LzYlVIrIh8f4CVb0g8X4PYFvi/dYora9Y5wfc908/QUplJhI/WtIe6uNFWJ1GL2jXk75L50m4yKDek371WM3K8zusZxadK63Z+XlRT9Bpe3aqa7eeeV6/qo7Wug6P3WkEjlhXjWiQn23kqatpcy/W1dp+gaAaWDV6Quz3Na60ZxNdOm3HfZvY+dATTK/dnRHncVVd3+Tzx4Bk4MJlUVpfGanOr1eTpPPbLKrDIHTSIN+w1XMQukHWSRvk6VV5vdS1cy55lNeNTnyP+a0z9Zn6RB/36F0L7C8ik9HQ54nAV/tVeMzIdH69nCRdfcJhPHTxBrQcNM3XN9MYA2xRjtFrRGqfskeBpQevGVjZ/WoqVd0pIn8IfFFENgG/6vdiFxihzo/Gk6RdN9r+//14Hvre9d0eJicadwvNO4xOu5NGulbH61TX7vEWii4Pel/n2g5jWNq2c1uU1Dx6d/VoHyF7p7Lv7x47Nq4WVfU/gf8cZB1GqaVbTpKKyDuAdwDst99+mQ+88ZINTBaV+bmAIIpoUOfVPYp6oFSjHnieIhp7Z69GQoDaaAleA13V07w2iAYRToqHuiiPpqIz1Hi218qiGJXqMEZrj/gJHdXFM74EmXQFSdcz0gXhYgQPpeCFm3zDdhF80cSkf3is2vIkqpOCxOVJOO+pg9cRXfOsOkErjpTTOpXYLZSmInC4dLStq41e0FpXjKIsZLHFcPFMlCbhgpC8bLE2OkM2W6yPzhBU7sdmtiik7+OqrqY8gSBKc+u0oS5Oi6O6lBO6gqc8/v0bCP74ZXiF/rtaDJ1ijNdYzCh1fi0nSaMVRRcArF+/PvMPvsd/eAuUAorRBHoSp1f3ZB5xe4P3Muni9xonReVV83lobZSFGl0lqc5DPYQ3XZIsHvHr6tlCF3v6qKsncRtU03xij/i155duu7C81LGSaTIcOi83ndZeqyHQtWeL4cKN+Ms+L1vMM6pDM1uMLTnzfeyIhJKHTssBO+76LUt/ZzCh3JtP+iw8RmmfX2WSNHp/InBpHgf2pgbkSX2BkPlXhmEYDdEgwFs0MehqjA0j0/mp6k4gniQ9lxwnSfd59XPGZqw9b6zjM6qM17BZ3kztuYLp/VYPrHxVyeU1KozUN36vJkmXPW0fKJdo/VUeD+xUEdLe5zvTSSVfK52rlunUznSuc8mqy0arejZaHtBq2UB8jHSe0dRlu8b1n0tlCK8/upo6C9RHYmjHptJH79wW29XF917vHDvX34/1CMsO24tBRXVQxm/Yc6Q6v17x4N/+gIKWw8ns6Msg+TWl0UR5uBKu+kUhUElrRxdHXgg0/LUlktKpVDa0h57tBaKFMyL91aXrGTsVDr1BJHWRt32ptkuoC/+uic6g4ftqebFOE3VqpSOzrvZaNdYl6zRculDlaoOqzn0d8tYlzyV9jWNdvW1oZYFNrW3U2mIA0IYNV2yxhS5p5811yUUfSnzFqguIqm1Xe293+r1Rq3vq6puZ37KD4vLFDIJxC2ZrnR+w8+7fAqFBphd7xDd05X3lJqcmrR1d/G/o5cShk1pderHAoHVhnvrzS3vgj724JPMI9Wkuz/256qhNM13nulbXGBrZRtTxa1JXa2c+tWmZbTGDzmnnDXS1Pw7qz8WVlsf3hlcssOuhJwbW+Y0bIzPn10sW7d93t3KGYRg1BPMlpvbebUClh8+ffXJsPRRY5wesfet/6yiqQ6eM2eiCYYwUA5l2E9jthKdRXDmYqA5KFEA5h9eoYMOewNyDj1IswHy5dl4OaqMl+BJugo/TPMKNxCWtdpzOKAteQCmo7s7zEhu8Y6rbl6u4PP6nPc2ndUJ99AKXR/xG55f0iO+RLYqET0CZWp10oIt/N7aqp0dA0ELn9Pifoby8dGHbaV27uNog6FCXzhNGZ2hui+5rFZbXvq79axxvqk9H4NCULrsNtx+5I4suHhIta/Pz8yQg0A5squ5eg9Ijj2L0D+v8gMcu+jFSmmfCr07ex7ijJaTTgpa6yULQ0Kt7+n19eXnrooUHjnrGHvFr09LludolP139+dXryKRrfK2GSdepTbmO3amu0LGus2tcbKLTio6MbZ5ug/zumWIirdH55WWLsw89zq6HNjG1djDbHUZpyDIPrPMDgp27Kn+nhzxcQyBZ0kxnurx1w1inXuiiLmKo6pnsAHtVnvge5e276jP1iVEasswDm/MDlp/wdOjhnN9A5hAMY0QZ19tFRFh00F6DrsbYYJ0fsPoVJyC5xz7p9HiNdK2O144uObfYj/IGqWv3eO3o2jlGv3WNjpGHrh/XuJ3z7n09a3/A5m2LIbu//Nl4xcEMxsV7OG3By5jx1CU/pVBQyvNBuBGW0Du7UvXgLvGGWQm9sxPlQSMv8glv93EejRaMeKINdBLOC9SVV3UWrUjFoCrRErQ9XcELENHQs70KnoRpSZ0nYZfoLC+QygbhuLzwXHLQabjwQACvkU6rnvtjXTXCQUoXp6UiGnSqi/eLNdZpZl1tpID2dDXXuJUu1eZpXc21cpaXRRcg4r7GQcL2Ja1rYhuxzk/YVPX6lRGkEokhvtea2WL9tYrKkw51KVuUqL2q97/DphzfG14iGkRVp2y94jr0nacjA4jqADbnN5Zs/enNSDmgkIrqILi9s9ekSdI7e3VCO5uuNo+7vDCUUhKnZ/smukrkhZRn+8zlDVrnaKv0+6qOpvlGVdfQFp269tu8Y1vsuW3EnXO/ymuhS9hiOFjkjgbhvMatokEEATN3P8z0YdnDsRmdY50f4C2abJ3JMAyjh2g5wJ+eGlj5ozRkmQc25wfs/qqTkFzG2mVoZ+vHzK4NY7QQmNh7FZP7DsbblBIF5c3hNSrYkx+w6IA1SFDC5du+lnhWpEpddAaN/9eeLk6tn+hP61p3ZOljK9WNxO2W16iEznTN6+k6sivVqUsvR3e2b+t6uuvU2uN/v3XZbaqXutbXqlMbTuskkdpuPbOQzRZTdRLAGcmiua5RDabXDS6c0ThinR/wxDcvo0AZ9apGm77ZhPgLNumdPYzOEHuDT0dLaKyrT0t6jG+ti/NpxdGwEk6euzzUJ9O6K689XTyM4jpOfrpoQUYi+kTN+UYLOdLt0jyyRec6V5v3QxfO60pde8a61te4W1vsn03VRxipnqmrDfK1xUYROMIuL0yTzLpk2s5rb6a0ZTuF5QNycTZCsfjywDo/YO7+jUBokOnLX/felUcaeJHPoJM28zTVtfJsn3d5LdI8aZ0nd12qDcTRBlk8/neqc7V5P3RJTyGd2nA/baMrm2rQBkGTNui1LcbRILopTyYKzD/y+EA6P2W0hizzwOb8gIl1trHUMIzBonMlinuvGnQ1xgbr/IDdz3hpOG42cjh+To4dY98ARkTVElrPDPak/G5MUWDJCc8c2JAn2Cb3saR038MUJ5T5XYnhFIJozL/aKRa8MoF6NfMznhdQDryELo5wLdH7cONrvPk41oWbztvXeZ5SLie8yEu4qTbOIxLpAi8xFBaWFyTLqwyXNdcl83SnC7eyN9O52iXc4C0EgVSO7ftB2AZaTYu98tfoVKrtGbVdXbtAbrrktZLEtUq2XSZd2jZSOqdteOEcWNI2PAlAqFz3cKN22ObJ88tuG7X17NSmemWLngfilSPbiNsgHhKt6rzIFpPt6Uf3cbXt6nW+RLaYaAM/ulZJnQDlxPeGS+dJdH5U7+3ygw8ySEao38oF6/yArf9xVRjVwdka5Zp3vtSPjMeeO2qpNaVws24qrV1dtNDB98v1spTOcxzb91uX59K5b4tOdPVtUNXFiwWgQPr86svP0nau8nqpy+Uax2kZrpVb57KNlA0njh3+UpeMttFZPftti3Gba6W64tA1s8Vm5WXTxTZcnQdsbcPlR59g/qFHKa7d01EPI2+s8wOCmdlBV8EwjHHH8whmBhPVQal9yh0HRnGiK3cWH38kZPGnNwS2sXCHJsLGXbjnZ/SbUYumIr4wccDagZWvOb1GBev8gKUvex4EySEIbfA3LdKTl78HOm3W/7ZT3nDqlEGdX+uhtN7rshwvf139T45ObTjL5wPQSfz5YK5xO9Eglrz4uIE5tR5HbNgT2HHZVRSKSjBfjazs+2UkXkCgWlnAQGLxRTjWH0YmiCfd47mMcKGDVib407pw0r+xzvM0nONDKJelsghApBppAqIIEdR6qI9vpsrCnMjTfJCMXtBQV92g20inUVojXTyXUdUpvh+XF3oF8LywXVSjdpHI4398fhovfAj3T1XK8+Izl8rvlWrbJRdMtKvTmv19aZ1GkTR6q6vOH/VC5yUct8c634+iM5QjWxQo+OF9EF8riaKQQKc2NThbLCTaoLI5vmIbiXst1lXapVxpp050yfs/qVNtbIu7rvop+tbX5ORqsU10tFZq5oF1fsCu63+FBGX8yo+u6q+2QiE9mZ003mT+2rQsOt9vrhNppqslS6SAbDqtGw5I6wS3R/zGuvgL27U4QOuOX1ee4Gzz9O6UQiEfHbh1aRaSLsxTTUteqyCg8gjTmU0N3hbj80s6A3DZYqt7NKsuff8ndXE0iDpbVGX+voeYOGQd/UaxTe5jiUwvGnQVnOQeX9cwOsaMseeUy3hLpgddi7HBOj9gye++CGqGGmQoFrcMRR3yYiGdy5gRPvTZBewpAoV99qSw12CiOgCRT9vuX6OCDXsCxb33wNMSQVuT242IOs4aK4hnKRJ5XMdNJUudroG0Lq1eFzvSbU4GnaP8eO6ueXmN2rbV+TnyjLKujiHQiUKTZe7x8Hv9Na5Py1JeL23ReS6JfCLhe60b43Pco3X3savQ9nXu9hQm99k9Xak+IgRj9gPHOj9g50UXU/DKaLFqs3G0BqCy+KLq8T9MT6fFCxHceZLzDZrSKeJVy0JBPJfOVV6jOvVLF9VdQQMh3ogv0ReMqiBe9XzR6OsiVV78cOEuTzvWNTq/jnVSTe9EF9OOLvml35GuQdu1ryP6wtbatMR1Ty646b8tdqZDoRINouE92rrtutPB3C9vJtiyFW/5MozeY50fUH4kEdUh8SWV/KJLpqW/yFy6VnkWkq7SkaQXI3ggiV/JIlG+DsrrVDeM7Vk5pzaP3bEu1zbXep1f/8Xe7NiNyhuULn4KHHQ9pVig/NjjA+v8RmnIMg9szg8o7De4jaULgTG7ZwyjJ+h8CX+NRXLvF9b5AYt/75UwOVGbmA7I1VP6WZZh5I/r6XZkGIa6T0ww9cLn4i0dXFSHccOGPYHiQetY8aGz2H7htynd9yDesqXgldEnN4fzGQAS7tPRpGd7L+GhPn788UKvGS11nhKUqzqRAPFa6Bwe8StzLZ3ooJInq078AFQq5YmneL6i5cRcT7SZulInoiFRldohJl8JylJtg3hjeJ90dW3XSOeFDgka6iTeVJ/QeRA7MuiZLqhtc5Ho/GKcttg7nfihU4ZKWq9t0YtsMZnHq70fndEgUrrY4UKQjBTiRVFIamwjZfuikU15LWzRoZvw8Nauo3T/Q8jiaaZf/mKmX/kyBoltch9TJn7nEHb7iw8BUN70OFvO+QCkNv4C4NiY7mXwUF+raxKdIb1h1+FJ31lep7p0PVvpKvM+Ll2LYzsGSOs3aQ+pzrEpe+A612bygesUvxB3AHGmHtliRZPK02l0Bj+5Cb0NXSfleTB17NOZ/uT/dmgHw5j1fTbs6WRubkSD2xqGMRKUy+iuwURwMELsyc+Bt9caZGoSne1dqCNlOKYaDKNnCCP5OBFuQ+hxIZOTTKw/pseFZEcZv2FPe7xxIJ7H5MknRe8aDQe1spTmusYdX6PjdlJeD3Ta4vNhqWcuulZDge3YRr91jY6RJT1HnTRI71V5ueraPW52m/IWTeIfekiL/H1EqexN7PY1KtiTnwMtlSj98Af4xSBaIBF5tq/sC4onygMQCMqACojiRc6xm+p8BQ+0TOUR0PPDv4MyQLRxPPrlHM+deL4myuuDrhCAB0EpcX6FWBdm9jyNP6pGGPDDdtEyxB7xxXfrwvL6oPPDCAA1usj5cGNdeOzmOtrWBSWJ2jy8LoJDFyTrpJl0zjboVJdou4qug2tcsam4PC8Akf7ZcKe6aG44KGvtvV1pA60sbqlcY6nOG9Z9JyTaTrwAKe2kfMsteEceiTEYrPNzEDy6EcphSCPf4RE/vbjDL0D6V19TnURLw9OtL+C7nsVTIb78fukkylOE9K95P9EG8Y/7Oo/76XZpoEsvpFjoOr9Ybxt1ukT4oaw6Z5061Ul9Hmj/GleFUVocwmeIbT/2blPVtXF+NeU10c3NUbrlZopD1PmN0h69PLDOz4FMT0N53EzBMIy+USggS5cOuhYVFJvzMwBv5W54Bx5Az3buavq3uWEYw0JfNuyLMHHiiX0oyGiEdX4NKKxeGY0/afVVGbZIjIuI1uaBTDrJosuSp1OdtKNL4PpiyJI2jjrX3ZVXeVmOnacua55WOiGxf68TW0zp6ESXvh8dOo/E0HMbukzlgTfhhSNMQ0T6TDp9jQo27OlAd+ygfOMG/GJQ3agbzdPFEQbCtNpoBUDCw0NC54Gi1UH1VJSDeON48jjV8mhcXt66aGKe0LlF/bGcx86QNqS6SiSNums8xrp+Xqsgkvbd9tP3saPtKl51Ejpa67KVp4gXMH/dNUy84L8xLIzbsKd1fg5085NQKECpVPkCiYm93bdMi3VxJwJ1k/BpXeZj90IX11GGvJ690HV7jRegri/XKv6xpSNgG3nXc24O/e0jGIPDOj8HsvtqKDtcj3VC4lfiUKOEK7MHXQ/DGAcmJ/H2P2DQtagQ9u3jdfdb5+dApqYovuhlzF/5/dDVmWEYPaEv3lSGAd+v/qD2fWTxEgrrnzPYOqUYlmFPEfkcsBPYDhwJvEdVN0afnQMsA1YCV6jqxVH6UcC7gHuBPYCzVbXUrBzr/BpQfMVrYeVulK64BN22DW/tvgSP3Atz5cTcXfSKN9DGaRKlEf1dpDbYlRcOpyQ3uVfmBhI68aM5hya6yrxOOZWnlS7arF4xDy9KS+r8xLnE9YyHgRLl4Ufvm+n8+jbIVUf0vpkufa1iXfJaZdG5rnGvdSTy9LK8RtfY651OitTboh/ZaxztIrKNlrafRee6r1rp4rbsxBYnPLx91hE8+igEZQpHHsvEq38fmUiFUTNidqjqhwBE5P3AnwPvFpHjgJNV9eUiUgRuE5GrgS3AN4EXq+pGEfks8GbgG80Ksc6vASLCxPNfxMTzXwTA/H9eTGnjvWFHlsZ3pCVbNrrh2tY1ypOHLp6DSJ9Pp+Vl2GjsPE6eulbHaZTWrzbvVpclzzDUs1NdbIvJzqWT8npgi5UFLB3ZYgDbnmDxZ/9vq4wDZVie/OKOL8IjfAIEOA24NsozLyK3AycBtwKL4qdD4BrgDKzzywedn4egg43v4zWMbhiGi1LTEbiBk1ywmgOrRGRD4v0FqnpBMoOIXA7s6dB+ODGUuQJ4CfCa6LM9gNsTebdGaZuAbY70pljnl5HCEcdQvvL/wXybc4DxkIlhGNkQcv0mHjiej//MYwddi37yuKqub5ZBVU9p9rmILAe+CrxNVZ+Mkh8Dkm5xlkVpjdKbYl/LGZG162DpkkFXIz8W0peLsfAYshET13aGzGgZ//jn51md/NFw2DOPV7eIyCrgK8A5qnqviMRPfpcAx0d5CsDhwNXAPcCMiKyJ8p0IXNqqHHvyy4jefRsyvx0mqhPs8aQ45WgCPZq8F8LJe9UojxelNdPFC1CiCXbxqURLUNckfB66ZJ286laHdD3TiwUq5SV1jcprpPOjj5K6ZJ2y6hLzqd3qJF6MkZMuXiDRC12ja9VvXV2bU1240pVNJds8rQuoLFLJbIslKqMwTXXlRJ1Sukp0hlKqXZroxAMKEPzoUvx172GY0eH5RXwFYd/0jxL6mtsGfFdVrxORq0TkPMLVnmep6mYAETkD+ISI3E84M3thq0Ks88tI+a7bkPm5ysrJGgqJH4UCIoI4FnI11cVar+Zt+G9607nDQ32uuvSm6F6Xt0B1APi90TW6Vv3W1Z0LhKs329Ul20B7cK3SC7u60PnRva2BZtYF996OkQ1VPabJZ59pkH4jcGY75VjnlxFZsgwKRSjND7oqhmGMGosWD7oGTVGGZ7Vnv7A5v4z4Rz3H8XPZhQzP4EEzhmxOxTBGjqz3UHEC/7kv62lV8kBzeo0K1vllRJYsp/jWs2B6KUxOwcQUTE3XR8SMHNvW4GrldB7XhHqnOtdNmfVYrXSeI18WXaf1dNU7S1pWXZY2z0uX9RqPmy5tU9JAl8V+sqRl1bnOL/lWxF3PxUvCUaKpRVAo4D3r+fjHv8iRcbgYlgUv/cKGPdvAO+hwJj70JfSR+0A8ypf8HTxyH+pLdfLeC2f4NbaCaA6w3ht8Ii16Dz3SeSCeF6bFiwVaHbtJeb3SiRedy6jo2rlW/db12qZy1jXME0cYUR0+W1TwxHPo5vHf8QEIAmSPvZHFwxO01qhinV+biO8j+x4EQHlLuP1EROp/FXq1Cc48/dJJIs1vkKeNNNNFacN0jUdc1zBP0l1YO7p+2EaDOuEXQAO8Aw5jlBgLH6sJbNizC2TvdYOugmEYw0ZpHlm996Br0RZK1W1pt69RwTq/LvBe8ArqJ/iGkDH7RWcsQIbxNmtQJzn8WBvqHAFs2LML9P7bYKoAu+Yr8wL4gCdQinocD/AlnAmOvc8n81R00Z2U1pU14TFewvR2dYVo+CjesBvXU6n1wF+I8iQ91Ht91gm1kQLiSBPNdLEVN9PFdSp3qYvPJasu2QbxT82gTV26vDiaQFZdIlpCV7p09IIs17hTXZ1tREOLadtXQluv6KL7qtU9k0UnUZ7KtZKo7bR6LrG9JmxKf3snqgGSaXX48DBKi1XywDq/LtCbrka0DJOp3bIibs/2WSIo9MKTvlCNnJ0mXSdXlJV+67K0Xae6rG3XS12aYa1nt9EZutUlrrEmO882dJ2WV6FZ2wWJTjDJ7AxsfAD2Wuc44JCiNudntEMnUR4Mw1jYiEBQbp3PGCjW+XWBPO3Z4cquEWDMftQZxuDw/NF66sMWvBht4j33d2H5KihOhgmFicj77ZBhPZ8x6gzrghfPr97/fgEKE3ivfCcyjN8DLVDN5zUqjMZjy5AiU9P4b/84+usb0IfuhOW7oddfCvOz0US5gudVPPdTjn4X+VGYh9glgkh14UpSJ9FMeze6KI94ktBR9UxTDsK6+Ylvl3IqT6zzJHxl1cV5nLq4UTLolHCIOalTGp9LUtew7UzXlk6guS06dBXbiG2RyD6z6qq2IRK5DXTaYmT7FZuqtf2aPDU21aGu0gY+FCeQZ54M8/MwvQzvmSciS1diDD/W+XWJ+AXk6cfB049DN95HcMNl0YKX1E/V+FdiEl/qJ8trdNHf3eriRThOneMXap3H/w51rjw1aXE9W+jq2i7StapT5jY3XfZr3G9d9G8cQaHTe8Z5P3aqS+Qpz8PcDvyXvKW+7iPGKA1Z5kFXnZ+Ea3nfDnwceKGq3pL47BzCiLorgSsSoemPAt4F3EsYav5sVS2JyBRwPvAwcAjwSVW9I9KcARxNuKj4blX9Wjf17hnTS6FsE92GMTb4PrJk9J/0FA3dtI0R3T75HQlcB+xMJorIccDJqvpyESkCt4nI1cAW4JvAi1V1o4h8Fngz8A3gPcADqvppETkiSnueiKwFzgaOVlUVketF5EpVvbPLuueOLNsd9lwHj9zFUE20xfuUDMPIGUGe8dxBV8LogK4WvKjqL6MggmlOA66N8swDtwMnAQcCi1R1Y5TvGuDU6O9TE5qbgSNFZBlwCnCDVn+WXAs444OIyDtEZIOIbNi0aVM3p9Y5i6c7izAg1HdQeemkz+VBakiJes/9/dDV5elU5/jl0Kmuk3OJh62TuNZTOHVZynPUsWNdB9chS3ke2WzD9SOvlzY1UQgjOCwALKpDChG5HNjT8dGH46FMB3sQdngxW6O0TYQh6dPpscb1WaP0OlT1AuACgPXr1/f9MuiOLcgjv4EJvz0v8tF76JEn/WRUhzZ0jTzb91s3VFEdkucyaF2hx+X1wBY7to2o2KGzRU/QX/8UOealjDoj1G/lQsvOT1VP6eC4jwFJ53bLorRG6a00B6fS7+qgTr1nx1PhcudyabgiBUify8tZN7RRHRa6bphsI+qoBl3Pujzledj8GMbo0at9fpcAxwOISAE4HLgauAeYEZE1Ub4TgUujvy9NaI4AblLVrcDlwLES/1wL81zWo3p3x4o1tuDFMMaJwiSy98Gt8w05ig17toWIrCRcubkceIeI/JOq/kxVrxORq0TkPMLVnmep6uZIcwbwCRG5n3Dm4sLocF8AzheRDxE+6Z0JoKoPicj5wOdEpAz89TAudgGQiSn06JfAjf8JpblBV6dKNGxjGCOLMBzjcl4Bgsjjt+fDoiVw8PrB1ikPRqzjyoOuOj9VfQo4N3qlP/tMA82NRB1bKn2GsCN1ab5JuEp06JH1p8Ky1eiNV8DMNli1Lzx2V7gJtuIIV0CkuqlXqE6uB1rdWFyJztBEF3uob6Xzvc50deV54abfikf8lC4+l4Wsa7h5vENd8idzOs1VT5cuvlatdFmuca9tsSNdqu0GYRvFIuy+FrZvhXIJDjwaedapSMHlEdsYdmyTe86ICBx2HHLYcQDohv+ATXdC0bE8z88w6uzK0q5OCG/kPMtLW47pOtdljXqQly4Pm+q3Lt0JJunbtVLY/gRyxqccGUcfHYpH6/5hnV+v0fJwDNcYhtE9ujD9oMQP3uOEObbuNQccw6hEfjAMowmeD+uOGnQtjJywzq/HyKr94GkngV8EBMQPX/1kzH7RGQuQQSzYiqI0AFCYhMUr4VmvGEBF+oNFdTByR579SvSg9XD/TeFqsU13hYtgAg2jFYiX8CIfUI0GES0EKJcBSXnSb6XTyPu8QMFLLJBooEvObSR1JMqLjw1QCgCHLjm32KmurOHwUk0bJHSSPJesOj/Vdll0ONo8qXNdq2HR9dimsuriTiuLTTWzxeQinF7ZolLb5l4BjjwNdmwOF7rsfySygEdxgjH7lbxwr+SQIbuvDW8gQC+6OrrBHItQnB74U5cpk07CGz0+RtbyoKprVp5rAU9euoJQMygh4I4KMGhd1mgJw6jr0qba0UH3uvTK5SR52mIyT6EIK9cghz+/Pp8x8ljnNwgml8DczKBrYRhGM4IyLFo26Fr0jVEasswDm/MbBIe/sP4Xey8ZM6M2jFxYuhpZsdega9EXlDCeXx6vUcE6v0FQnGjggT+VliWagEhrnevYHZfn0jnMyKVrlSdPXZY80FkbdKrLtc17qctgU7nqXPVM702UBroWx25UXss8AhMLI1qD4caGPQfBr68M5yC8QtV7RHzzJ7zII1L/3pWnlS72VtGurq3yvO7rOSy6ePynE138lN1PXfy93fNr3G9dwqbi9ulnPbdvQrc+iixzBbVZeFgwW6P37IoiNDm8yNelZcljunx18ZdjJ7phO5eFoou/mPtZT/FhZguMQ+entsnd6Acr1w66BoZhtCIowYq9B10Lo0dY5zcIjnh5ddN7hfTP0hwZs190xgKkh7dHtYzE16FfhIOORyaX9KHgwRMueNFcXqOCDXsOAFmxN3ryu+C2/4SnHoLplTCzCeZ3QSnyBepJuCG4HISbsIXqBt1yEHmfl6qX/FiXjOAQ63yvmhbrClLdJN1MJxJuBtZ+6KKN6LHOk4QH/pSuFJ8L0f4whZLW6uKNzV3pos3bdTpHPQOt3XTeqS59jeON4UnbSOoa2YZL14lNZdUVBMhgi+LQZbEpcrDFOIJDWlechJUHwOZHYGIaDnke7HdM1/f6KDFmU37W+Q0KWb4XHP8mAHTLw3DtN6qdVBLPg/SiM5eH+mY6IfpSc+jSG4uzljdo3YSrDVzve6lLb6SmfuN0p7os17hTXR421QtdM9toGtUhD5sqw5LlyAlvqc9nLEis8xsGxuwXl2EMJeP26JNAR2zIMg+s8xsGlu8VzjGUhyj6u2GME34R9jlq0LUYKOPW99uClyFAxINjXw/+RLQQhvBfyenyjJlRGwuQvBe8iB/5NxXwirD2GNj9gJwLMYYZe/IbEmS3deiLzoFHboa5HUAAD14Hpbn6+Y4gWqARLzyA6gKN5Ib2tC6OOSJSzZfUeV74vt3ykjpXnrg8qC7sSNdzaHVemLeZLj52pzqN0tLXqhNdJtvoVOewqV7o4vZTrdbTS+saXOOstrhkT1jzzHC0ZfVh4Rz8mGPDnsbAkOIi2P/ZAOht/xHuM/K8+ufzOHRMTZrLA3+6gISusgjGsVigk/KcOof/UmekgDHXCfW+XjvVZbKNTnWuOvVAF3eULvt06VwLXlq13c4nkANOqD/WmKJAMGbjnjbsOaxMraAu7IxhGPlQnB50DYwBY53fsLLPUbVParnTy2MbRhf02jS9Iqyzp740mtN/o4J1fkOKTC6Fo8+AqeXhzeoVYHIZpCNJi9R3kq49TWkv+fGwZ7M8kM1zv7O89LFz1GWpZ1bv/s7ysujSe9lSOtexG7V5S12GNnDVKYvOeaxObSqn8hrtSe3UFiemw/vHnwjvpf2eA/s+uz7fmDNuIY1sXG2IkRX7oif+L5jZHHZ6N30T2BWuBFVqv0zTHun9xC8wV57kYoRkvvR717HjtGS+VuU107U6F5euHc/9g9C1ypOHrp1rlYeunWvcaXmQ8DSTky16Ake+EQqLYGo54luoIsM6v6FHREL3Z4DO74wT659OXB7q6w/m1jV73w9dq3MxnTutnTYfFR1U0/OyRfGgNIss37f+eAZAxbfnOGGd3yixdC/Y8sCga2EYo0VQgqVrBl2LIUfHLp6fzfmNEge+MJyz6Fc0CMMYBN2adNI5hFeEvY5GJsYjOoORHXvyGyFkyRr0qDfBA9fA9o0wtRJmfgvBfNVrfU0Eh3gjsx8OBZXL4XxKHBUAgSDpub9QzVPRUT2Wl3C8XSoDGqX5kS4Ijxnvu0rrlDAfUXmxrhyE9Wuo09rN477n1pXKkef+dnQCpVIOOj/aXJ3SlaJ26UpXCNs6Pj8/bvPkNfbCL/1yDrpiIczfUieNbaNG18ymJBqmTNii74XHq5QnUYQKlw0T2n5si8VJWLIvbN8EhUnY51mw+vDM99g4Y8OexlAji1fD014JgG57GO74N5DAvek8TTMP/PHKOafOsZk7U3kZNoE7PfC7dKn3ueom8tE5Vxq62qVTXafXuJe6Lm0j9jrkqkNH9SzDxCTyrHfW5zMaYnN+xoghkKfBxm6jDGOkMRteCIjIh4D3qOqqRNo5wDJgJXCFql4cpR8FvAu4F9gDOFtVS82Ob53fKLNkTbh/KZjv/liq5OZI2zAGhVeEVTbM2Qk6RLv0ROQFhB1cMu044GRVfbmIFIHbRORqYAvwTeDFqrpRRD4LvBn4RrMy7NtuhBHx4ODfrW6CxwOJ/u3ogHnWzjAy4Nqo3pY+is4gke3vfjgs2z+/+o0NYTy/PF7AKhHZkHi9o52aiMiewO8DX0p9dBpwLYCqzgO3AycBBwKLVHVjlO8a4NRW5diT34gjS/dBj3wHPHUnlGZA5+GJm0Nv9fEChhrP9upO87x6T/rxogaXJ/1MUR2a6bwwrRe6pFeYWNdOtIRudOVE2/VLl8zTta7ZNY5soxtd/AsrGWGkps0TthgvBGpmwxMrYfcjQntfvg6ZroyQGYPjcVVd3yyDiFwO7On46MPAK4CzgeWpz/Yg7PBitkZpm4BtjvSmWOe3AJDCJKx+BgD6wBWgpeoXSJIsixNEqmkSv3f8Mu9pNIge6kQcuoztkkU3ClEkOtV5QkfRJ5K6ZMSGNM5r3MKG5zcjezyz/lhGW/R7wYuqnuJKF5H1wDzwTsJhz0Ui8mfAd4HHgKWJ7MuitEbpTbHOb6ExsSIcAtLhGb83jJ7hLxp0DRYMwRDM+anqBmADgIisA85U1U9G7y8B/k/0dwE4HIjn/GZEZE009HkicGmrsqzzW2jsdjg8fmM+nV/sQzGVRCo5na1RnlHUtXPsUdENss1zRQqwx9G9OvqYoagMvvOLEZGDgT8gfPL7EPA5Vb1ORK4SkfMInwrPUtXNUf4zgE+IyP2AD1zYqgzr/BYYMrEEPeB0ePAHML8DUFi0JyxeA0/ckthTVYDybK1T4HhDcuVg0XxMlEchWhGacoXk+2hSB+B5aBDUvEdb6KLyeqJz5RklHdS2gUvnefXtktb5fnjsdsvLco0b1bPmfAWQ2jSXg2q/CFIMbRSF6b1g0Sp44tZIA6w6KnwZCw5VvYtw3u/sVPpnGuS/ETiznTKs81uAyOK90MPOgNIOkAJSmAJA9zwuTCvtgIe+H0WHSHzhiNRuInZ9KXleTYdY1fnNdQ098PdJ186xR1nXy3Zp99jpxSptlefBsgNh92PBS9jwmudAaScUphEL9pwb/Z7zGwbMehYoIgLFWn+G4vkwsQyd2xx1YNQvPGjkSb9VmukWvq7dY4tECzw7qZNCaUedT07xCjCxrF5rdM0wzPn1E9vnN45MrbYFMUbvEen8WUJ8WGwhiIzeYZ3fGCKFRbDyiHDBQDUVhZovK028DKNTmtlU5e+kdyHxwJ+GFYf1qYYGaG7b3EcFG/YcV1ath6lV8FS0IX5qN5h9BEpz1UUMiSgLGm8srkSDSGxITkcmSEZnqNERRYOgNlpCvPerEEWfiKMCdKtTElEksur88KnYFdWhoiuECy4qUSSk2gbJzeqeQwdhVAfIphMJ66kOXTrKgufVnksmnR+m1+k00eYOXRyhohzriPbyafW6J6IzqEOniTbX4kpEJbTFJetgtyMQz+Fw3OgJCgRDtNqzH1jnN6aICCw9IHwBwaarYS4Iw9mkcW4CT0WDyBIVwJXm8tzfaUQDl84VRaKlznNHdWil8zLq0lEdutG1qlNXujaiQVSiM0h9vizRGYKtsN8bEPMva/QJ6/yMEPvSMQaOYxGM0TfGbcGLdX4GALLkYHTnA6FrtHaIf+3b95bRMQKL9g1HI4wBoWPX+dnPfQMAmVoDSw8LV9nhR4thfPAXU4kUIYXo87S4z5U1ho94+LsOL3r5VTvyl9TaVGEpsvtz+lpdw7AnP6OCt/JYdMkhMPNIGCZpet/Qy8au38L8FtSbgG2/AOKFDxoOlyajQQjVIdSk5/44okCQWEgSv483N7siTcRfqC5dI4//o6RLR9KoiZYwJDqR8Jq2usY1bRDbRgGm1iFxhze9X2RTG2F+MxSXwdReNtc3YMJVt+P15Gedn1GDFJeFX0hJFu0dvnY9hHrRSsBm0SCSaZW/G+RxLobIuFCmk6gVw6hr1HbDoktGYmh2jZ3lBVDeiqxIRbhZtFf4MoYEJaDcOtsCwn5uGdnxF2O7/oz2ECgsbZ3NMPqMPfkZ2SmsAH8ZlLZAiyGS+igA4Ttpmgdnnjx07RzbdKk0qW7lS+ZxT/UKtT+QPGTxoc6cxnAxbsOe9uRnZEZEkN1Ogsk1VBYy+Etg6dHhv1Ga+kvA82u9eURzQjVpnpfy+CH1Os+v1YmAX6g/dvpYKR3d6Fzn0pHOy1jPDDq/U12hXieO65Bu80IBTek0jrwQL1yZPhQm90nYxmJk5XORwvJmZmUMAYoSSJDLa1SwJz+jLcSbRFY+Fw3mgTLIJCKCTh8MOouqD1svBwqgTaIAxO+9xFNCJWJEobWu0+gFHUc9aFKnXuj8XuoK1bT0v810vp8qrwDTxyCFFeBNINFKYA1KoCXwJm37gjG0WOdndIR4RaBYfS8CMgWqaDxxPugoBKZrnZbuPNs6lgJzSCqaehhqyL5aRo1xW/BiFmrkioiE84LlrYOuitFzFCnsPuhKGLmgNudnGN3iTR8JRA6Maz9J/O14qkgOsxn9oeGwpKT+Fmqvnw8T+4b79wxjBLEnPyN3pLAb3rKTCHbdFa4MLaxAJg9Ggm0Es/eG80H+EuCJ0It/xeN/uBijEinAtSk73o+W3BzvRx1tkI7OoLVRHeJoCck8UJ82SF20kCWTLhlJw7WpPqlr1nZxebGuUACmIZgAnUWKe8LEOph/BJ17CMRHJg9Aint3byzGUKBAoDbsaRhdI/4y/MXHpFKX4k+EX5jl2Q3h5mmvWC/OsvE9nUckW2QJ17E7KW8YdMn3aS85WXUN02bxpk9EkjEf/YNg6qB6rbEAsGFPw+gTZnqGYQwOe/IzBoJXWEsw9yTktcKs8a5roxO83Wqf+owFj9pqT8PoA95q8PeB8kNRQtRzSRF0nmpvNgnsouo5JF4Uk1ocYx1fA+Lh1HhLwxLQ7VQbLABZDDpDpc1lAm/iiL7X1Bgk4xfSyDo/YyCICP7E09BgPzR4EpFi2CHiQfAkqjsRWQLeCqCMlh8DAlS2gTxWXdwCtZEJnBEGSGyg1+rCjuRet0r0CUfUg3ghSfrYw6yrtIuP6KHAPOKtRLwlqM6i5U2Ah/irCTe9b0GDbYgsAm9325xuLHis8zMGiniLEW9xbaK/O0Jy/1gBKUQLZYIbAHUv2mgZYYDaVZCNdK48WcobUp14KxGZTMgmkcLaVKYViLeivkxjLAjHVOzJzzCGmGlg26ArMWLYbW60QtEx2+pgS+6MkcKT/TCzzYoH7FPxuWkYRhX7FjFGCpHFeHIEsIiq55E9gf2oepUpAgcBuyXyLEU4FJii6q1kDbBvQjcR6VYmdMsi3WRCtxewNvo71h0CrEjolkdpSd3ewD4J3WSUZ3lCtyJKm0jo1kbaqi6s07KEbiVwcEq3D54c2HFbG+NFkNN/o4KNhxgjh8hKPJ5NuE3CQyT8Dad6QJTmRws29o2GcrSybF91rwa6wKGjGqnAqTswpdvHodvboTuoQ93BKd3eDt0+QCnKY79tjayobXUwjFFA4pA6LdP8DHmi+HMLRufwmmMYRg1ddX4i8jlgJ7AdOBJ4j6pujD47h3BcZiVwhapeHKUfBbwLuBfYAzhbVUsiMgWcDzxMOO7zSVW9I9KcARxN+FP4blX9Wjf1NgzDMKoooDo6Q5Z50O2T3w5V/RCAiLwf+HPg3SJyHHCyqr5cRIrAbSJyNbAF+CbwYlXdKCKfBd4MfAN4D/CAqn5aRI6I0p4nImuBs4GjVVVF5HoRuVJV7+yy7oZhGAYwjpvcu5oUiDu+xLG2R3+fBlwb5ZkHbgdOAg4EFsVPh8A1wKnR36cmNDcDR4rIMuAU4AbVSryba4GXueojIu8QkQ0ismHTpk3dnJphGIaxgGn55CcilxMup0vz4cRQ5grgJcBros/2IOzwYrZGaZuo3aQVp8ca12eN0utQ1QuACwDWr19vweEMwzCyoIzdPr+WnZ+qntLscxFZDnwVeJuqPhklPwYsTWRbFqU1Sm+lOTiVflerehuGYRhZsZBGbSEiq4CvAOeo6r0iEj/5XQIcH+UpAIcDVwP3ADMisibKdyJwafT3pQnNEcBNqroVuBw4VqrOBo8HLuum3oZhGMZ40+2ClyuiY/xj1DdtA76rqteJyFUich7has+zVHUzVFZufkJE7ifcXXxhdKwvAOeLyIcIn/TOBFDVh0TkfOBzIlIG/toWuxiGYeSHrfZsE1VNh+pOfvaZBuk3EnVsqfQZwi0QLs03CVeJGoZhGLkzfpvczQWEYRiGMXaYhxfDMAzDhj0NwzCM8WPcOj8b9jQMwzDGDnvyMwzDGHN0DN2bWednGIZh2LCnYRiGYSx07MnPMAxj3FE1356GYRjG+GG+PQ3DMAxjgWNPfoZhGGOPjt2CF+v8DMMwxhxzbG0YhmEYA0REJoD3ATsIw+E9oap/Hn12DmFM15XAFYmA6kcRBka4lzDY+dmqWmpWjnV+hmEYxjAteHk/8CNVvRpARJ4Z/XsccLKqvlxEisBtInI1sIUw6s+LVXWjiHwWeDPwjWaFWOdnGIYx9gzVnN//AB4QkWOA3YEvRemnAdcCqOq8iNwOnATcCixS1Y1RvmuAM7DOzzAMw2hFjp3fKhHZkHh/gapekMwgIpcDezq0HwbWAaqqnxeRFwP/AryAcDjz9kTerVHaJsJA6un0pljnZxiGYeTJ46q6vlkGVT2l0WcishW4Lnr7E+B5IuIDjwFLE1mXRWmN0pti+/wMwzDGHgWCnF5d81/AgdHf+wN3a+h+5hLgeAARKRAuhrkauAeYEZE1keZE4NJWhdiTn2EYxrijQ7XV4Rzgo9EKzqcRzt+hqteJyFUich7has+zVHUzgIicAXxCRO4HfODCVoVY52cYhmEMDar6MPA/G3z2mQbpNwJntlOOdX6GYRhjjjJUWx36gnV+hmEYY89QbXXoC7bgxTAMwxg77MnPMAzDACyen2EYhjFW2LCnYRiGYSx47MnPMAzDIKcN6iODdX6GYRhjj4INexqGYRjGwsae/AzDMAwUHXQV+op1foZhGAbjNudnw56GYRjG2GFPfoZhGAaoDXsahmEYY4WO3ZyfDXsahmEYY4foAn3UFZFNwP0NPl4FPN7H6mRhGOsEw1mvYawTDGe9hrFOMJz1GsY6QeN67a+qq/MoQES+H5WTB4+r6ktzOlbPWLCdXzNEZIOqrh90PZIMY51gOOs1jHWC4azXMNYJhrNew1gnGN56jTo27GkYhmGMHdb5GYZhGGPHuHZ+Fwy6Ag6GsU4wnPUaxjrBcNZrGOsEw1mvYawTDG+9RpqxnPMzDMMwxptxffIzDMMwxhjr/AzDMIyxY0F4eBGRzwE7ge3AkcB7VHVj9Nk5wDJgJXCFql4cpR8FvAu4F9gDOFtVSyIyBZwPPAwcAnxSVe+INGcARwNl4G5V/VqTOnnA24GPAy9U1Vui9HXA94GNUdYbVPV9/ahTs3oNsq1S9fsI8IJE0idU9T/zrl+eiMiLgVcDjwGqqh/Nu4xUeT8DdkVvy6r6IhHZDfgkcA/huX5QVR+N8rfVbm3UYw1wLnCkqj4rSsutHp1cvwZ1+ggDtCkROSiq0y+AtcATqvqxQbfV2KOqI/8Czk38/X7gS9HfxwH/L/q7CNwJrAAEuAVYE332WeDM6O8/A/40+vsI4MfR32uBG6nOk14PHNKkTkcDRwH3Ac9IpK8D3uLI3/M6tajXwNoqVb+PNEjPrX452940cBcwGb3/LvCiHtt7XRsBfwX89+jv04F/6LTd2qjHa6OyNvSiHp1cvwZ1GqhNAc8CXpF4fxtw7KDbatxfC2LYU1U/lHjrET4BApwGXBvlmQduB04CDgQWafR0CFwDnBr9fWpCczNwpIgsA04hfEqLVwhdC7ysSZ1+qao3Nvj4dBE5R0Q+LiKHR2k9r1OLeg2srdKIyJ+LyNki8n4Rme5B/fLkeOB+VZ11lN8rjoja5iMiUneuqTp00m6ZUNWLgG2p5Dzr0fb1a1CngdqUql6vqt9LJHnADgbcVuPOyAx7isjlwJ6Ojz6s1SGBFcBLgNdEn+1BaDgxW6O0TdTeIHF6rHF95kr/ExF5W7M6OdgUfX6riOwJ/ExEjm5Sbrt12iNLWznodVvFmqbXEvgOcJ+q7hCRPwK+BJyZc/22OsrulKbn2iM+pao/FxEfuFpEtqXqsRVYKSIFOmu3bsizHnldv6GxKRF5FXC5qv5aRIaxrcaGken8VPWUZp+LyHLgq8DbVPXJKPkxYGki27IorVF6K83BqfQvqOoX2zyPHcCt0d+PisijhPOUD+VUp7tatVUDet1Wd8Vv2qjflcA5PahfnvSjjBpU9efRv2UR+TFwcqIem6M6PKXhXFAn7dYNedYjlzqq6q2JtwOzKRE5mfBavSd1rM0MSVuNEwti2FNEVgFfAc5R1XtFJH7yu4RwWIroF9XhwNWEE8wz0eQ4wInApdHflyY0RwA3qepW4HLgWBGRKN/xwGUd1PVN0XERkSLh/Nh9g6xTxFC0lYh8JvH2EKqdZp71y5Nrgf1FZNJRfu6IyO+IyJmJpLiNKueaqkMn7dYNedYjl+s3DDYVDU+fAvwJsEZEjmcI22qcWBCb3EXkF4RPsfET3zZVPT367BzCFVMrgcu0dtXUuwkjP+xGddXUIsJVU78lfHo5T2tXMK4nXMF4hzZf7bmScFXW+4B/AP5JVX8mIi8E3km4IORg4Ceq+rf9qFOzeg2yrVL1+wvCRSSPEU7cfzhxzNzqlyci8t8IF1psAua1h6s9RWRvwh96vyD8dV8EziJcEPEpwjY4CPgzrV05mLnd2qjL84E3AS8F/i/hAoxFedWjk+vXoE4fZoA2JSLHAj8CNkRJiwmv4cWDbKtxZ0F0foZhGIbRDgti2NMwDMMw2sE6P8MwDGPssM7PMAzDGDus8zMMwzDGDuv8DMMwjLHDOj/DMAxj7LDOzzAMwxg7/j/FrcHb3uKH4gAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=[7,7])\n", "plt.scatter(x_ep,y_ep,c=opd)\n", "plt.title('OPD')\n", "plt.axis('equal')\n", "plt.colorbar()\n", "plt.show()" ] } ], "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.8.8" } }, "nbformat": 4, "nbformat_minor": 2 }