{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Simulating an Add-Drop Filter" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "## Imports" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "slideshow": { "slide_type": "skip" } }, "outputs": [], "source": [ "%matplotlib inline\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import photontorch as pt" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "## Simple Add Drop filter\n", "" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "## Simulation & Design Parameters" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "slideshow": { "slide_type": "-" } }, "outputs": [], "source": [ "dt = 1e-14 #[s]\n", "total_time = 2000*dt #[s]\n", "time = np.arange(0, total_time, dt)\n", "c = 299792458 #[m/s]\n", "ring_length = 50e-6 #[m]\n", "transmission = 0.7 #[]\n", "wavelengths = 1e-6*np.linspace(1.50, 1.6, 1000) #[m]" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "-" } }, "source": [ "## Circuit" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In the [all-pass notebook](01_all_pass_filter.ipynb), we saw how to create a network by subclassing pt.Network. Although this is the preferred way of creating a network, sometimes you want to get rid of the boilerplate of creating a network. This can be done by creating the network using a context manager:" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "slideshow": { "slide_type": "-" } }, "outputs": [], "source": [ "with pt.Network() as nw:\n", " nw.term_in = pt.Source()\n", " nw.term_pass = nw.term_add = nw.term_drop = pt.Detector()\n", " nw.dc1 = nw.dc2 = pt.DirectionalCoupler(1-transmission)\n", " nw.wg1 = nw.wg2 = pt.Waveguide(0.5*ring_length, loss=0, neff=2.86)\n", " nw.link('term_in:0', '0:dc1:2', '0:wg1:1', '1:dc2:3', '0:term_drop')\n", " nw.link('term_pass:0', '1:dc1:3', '0:wg2:1', '0:dc2:2', '0:term_add')" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Simulate Time Domain" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "slideshow": { "slide_type": "-" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAawAAAEGCAYAAAA5T6EkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deXzU1dn38c+VBAi7bLKEHRMlICggi7fWKn3UWG6X2iq2ta2tRWtttS5FX61KXfrYFp6WblK07vamLthaq6i1ttBbqYJlERJoRBCEyBIJ+5Lkev6YCR2TmclkMr+ZDPm+X695MfM7Z37nmiGZK+f8zpxj7o6IiEhLl5PpAERERBKhhCUiIllBCUtERLKCEpaIiGQFJSwREckKeZkOoKl69uzpgwcPznQYItJKLF26dLu798p0HJKFCWvw4MEsWbIk02GISCthZhsyHYOEaEhQRESyghKWiIhkBSUsERHJCkpYIiKSFZSwREQkKyhhiYhIVlDCEhGRrJB138NKp3U71/Hi+heJtQVL+7z2fGH4F8jPy09zZCIirY8SVhyPlz7OU2ufwrAGZU4oiRX3KGZSv0npDk1EpNVRwoqjxms4tsOxvPq5VxuUrdi2gi+88AWqa6szEJmISOuja1hJyrHQW1fX0xIRkWApYcXh7lGHAwHMQsdrvTadIYmItFpKWEnKtVwgNGwoIiLBU8KKw/EjPan66npesWYQiohIailhJanuGpaGBEVE0kMJK45417DqhgSVsERE0kMJK454MwA16UJEJL2UsBoRq4dVNySoSRciIumhhNWIWJMu9D0sEZH0UsJKkiZdiIiklxJWHPGmrGvShYhIemktwTjiTroIX9taX7Wetz98O2qdoV2Hckz+MYHEJiLS2ihhNSLWpIv2bdqTYzk8tOohHlr1UNQ6k/pOYu7Zc4MMT0Sk1VDCiiNeD6tL2y489d9PsWP/jqjls9+ezZ7De4IKTUSk1QksYZnZg8AUYKu7j4xSbsBs4DxgH/AVd48+tpZBsWYJAhR1K4Ju0cseXf0oOw/sDCgqEZHWJ8hJFw8D58YpLwEKw7dpwH0BxpKUeCtdNMYwTXkXEUmhwBKWuy8EKuNUuQB41EMWA8eYWd+g4kk3M9MMQhGRFMrktPYCYGPE403hYw2Y2TQzW2JmS7Zt25aW4CD+au2NSbZnJiIi0WVy0kW0T/SoY2juPheYCzBu3Lj0jbM1o6XmDAlW11azvmp9zOfn5+UzoPOA5IMTEclCmUxYm4DIT93+wOYMxRJT0j0lS36vrPuW38fcFfGnwz9y7iOM6T0mqfOLiGSjTCas54BrzWweMAGocvctGYyngeZMmmhOD2vngZ10atOJH5z6gwZlW/ZuYeaSmew4EH06vYjI0SrIae3/A3wS6Glmm4A7gDYA7j4HeIHQlPZyQtParwgqlkzIsZykE1aN15Cfl8/Zg89uUPZe1XvMXDKTwzWHmxuiiEhWCSxhuftljZQ78M2g2k+F5k66SHZI0PEji+vW1yanDQCHa5WwRKR10UoXcSSbcCA0rT3Z59fU1jSasA7WHIw5bd6wpBOtiEhLpYTViOZMT092SNDxI6vB15eflw/AXYvv4q7Fd0WtUzKkhB9/4sdNbressox5ZfNilrfNbcs3Rn+DbvkxlvcQEQmQElYcmZp0UeM1MRNl13ZdufPUO6nYVxG1fMF7CyjfWZ5Uu8/++1nm/3s+vdr3alBW7dVUHqhkVK9RTBk6Janzx1JTW8OG3Rti9kjb5LRhQOcB6jWKtHJKWI1IemmmZgwJ1notuTnRe1gAFxVeFLNs3c51lFWWJdXuvup99O7Ym1c++0qDsqqDVZw27zS+/4/vc8f/3hH1+ZMHTubHZzTs2bk7c1bM4YPdH0R93t82/Y2qg1VxYyvoVBBauzGK7vndGdRlUMN2cf619V8tYsWRlvBF8pYQQ6ZDOK3faVx6wqWZDUKSpoQVkJxmLCJS67VJf7jk5eRRXVsdtay6tpqr/3I1W/ZE//bA9v3b6d2xd9Syru26ctvE2/hgT/Sk88bmN3j1/Ve58uUrG5TtOriL0spSAPp2bLj6Voe8DvTv1J+vjPxKgzJ358F3HsTd2byn4df0Pjr4EVv3bY0aU51+HftldF+y5lwLPZq0hLU1dx7UgtTZTAmrEUkPQ1nyv6C1XhvzGlZj2uS0iTmDcOfBnfxzyz8Z2WMkA7sMjFrn1H6nxjz3JcdfErNsYt+JzFk+h0M1hxqU5eflM6HvBG6feHvMduMpGVISt/xA9YGYvajcnFza5bZrcpsi0vIoYcXRrFmCJL/4ba3XJp0o83LyYias/dX7AZh6wlQuOO6CpM4fy6R+k5jUb1JKz5mouokoInJ0U8KKo1mTLpp7DSvJHlbb3LZUHqjk1P9p2FOqS6Ad2nRI6twiIpmkhNWIZn1xOE7Cq9hbwUcHPopatuvQrpjfw2rMxYUXx227XW47JvSdkNS5RUQySQkrjuYOCcay9/BeSuaXxJwcATC299ik2i3sVsj08dOTeq6ISEumhNWIIKa17zu8j+raai49/tKYkxxO6H5CUu2KiBytlLDiCOqLwzVeA0Bxj2LOGnhW0m2IiLQmmdxxuMVzvFk9rFizBOuGApOdWCEi0hopYQUkkR5WvNUsRETk45Sw4vFmzBI0I9aIYk1tKGHlmUZkRUQSpYTViKSHBOP0sKo9PCSoHpaISMKUsOJo7tpnMYcEwz0sXcMSEUmcxqTiaE7CyrEcKg9UUvJMw3XwDtWG1tvLy9HbLyKSKH1iBuSC4y7gQPWBmEmvfV57RvcaneaoRESylxJWHO6e9KSL0b1GKyGJiKSQrmGJiEhWUMKKozlfHBYRkdRSwmqEEpaISMughBVHS9jSW0REQpSw4mnGShciIpJaSlgiIpIVAk1YZnauma0xs3IzuyVKeVcz+5OZLTezVWZ2RZDxNJUmXYiItByBJSwzywV+BZQAxcBlZlZcr9o3gdXuPhr4JDDLzNoGFZOIiGSvIHtY44Fyd1/n7oeAecAF9eo40NlCF4o6AZVA7H3j08zdUQdLRKRlCDJhFQAbIx5vCh+L9EtgOLAZWAlc5x5j18MM0JCgiEjLEWTCivZJX3+e+DnAMqAfcBLwSzPr0uBEZtPMbImZLdm2bVvqIxURkRYvyLUENwEDIh73J9STinQFcK+7O1BuZu8BJwBvRlZy97nAXIBx48al9MtRr3/wOjsO7Ihatm3fNjq37ZzK5kREJElBJqy3gEIzGwJ8AEwFPl+vzvvAZGCRmfUGjgfWBRjTx+zYv4Or/nJV3DpnDjgzTdGIiEg8gSUsd682s2uBl4Bc4EF3X2VmV4fL5wB3AQ+b2UpCQ4jT3X17UDHVd7DmIADXjbmOcwadE7VO74690xWOiIjEEej2Iu7+AvBCvWNzIu5vBs4OMoZ4asPzO3q278mALgMaqS0iIpnUqle6CF0601b1IiLZoFUnrFpCPSytFygi0vK16oRV4zUA5LTut0FEJCu06k/quiHBHGvVb4OISFZo1Z/UdZMulLBERFq+mLMEzey5BJ5f6e5fSV046aWEJSKSPeJNax8OXBmn3Aitxp616nYU1qQLEZGWL17C+p67/z3ek83sBymOJ62O9LBa98ioiEhWiPlJ7e5PNvbkROq0ZBoSFBHJHkl9UpvZtFQHkglKWCIi2SPZT+qj4qKPEpaISPZIai1Bd/9NqgMJyrs732Xrvq0xy0CTLkREskGjCcvMbo923N3vTH04qffIqkd4tvzZuHU6t9GeVyIiLV0iPay9EffzgSlAaTDhpN7XTvwaFx53YczyDm06cHy349MYkYiIJKPRhOXusyIfm9lMIJEvFbcIg7oMYlCXQZkOQ0SOIkuXLj02Ly/vAWAkrXzFoBSqBd6prq6+cuzYsVGv4yRzDasDMLRZYYmIZLG8vLwH+vTpM7xXr14f5eTkeKbjORrU1tbatm3biisqKh4Azo9WJ5FrWCuBuv+QXKAXkBXXr0REAjJSySq1cnJyvFevXlUVFRUjY9VJpIc1JeJ+NfChu1c3OzoRkeyVo2SVeuH3NOYQayLXsDakNCIREZEkJLvSxfOpDkRERBKzffv23HvvvbdXpuOob/z48ccvXLiwQ/3jP//5z3t86UtfGtjc8yc7u+XrzW1YRESSs2PHjtzf/va3xyZav7a2lpqamiBDSotkV7rYkupARESy0c1PLx+wtmJ3g15FcxT16bzvJ58dvTFW+Y033th/48aN7U444YTiM844Y9exxx57+Nlnn+1+6NAh+/SnP73zpz/96eY1a9a0LSkpKTz11FN3L126tNPMmTPfv/baaweNHz9+z9tvv91p+PDh+7761a9uv/POOwt27NiR9/DDD68788wz90Vr77XXXutwww03DDxw4EBOfn5+7cMPP/ze6NGjD+7Zs8emTp06ZO3atfmFhYUHDhw4cGTZoNmzZ/f46U9/2rdXr16Hhw0bdqBt27bNvubXaA/LzArN7GkzW21m6+puzW1YRESSM2vWrE0DBgw4WFZWtvqcc87ZVV5enr9ixYrS0tLS1cuWLevw4osvdgJYv359/hVXXLGjtLR09XHHHXdo48aN+TfeeOPWsrKyVe+++27+E0880WPJkiVl99xzz6Z77rmnb6z2Ro8efeDNN98sKy0tXX3HHXd88N3vfrc/wMyZM49t37597dq1a1fffvvtW1avXt0RYMOGDW3uvffefq+//nrZokWL1q5du7Z9Kl53Ij2sh4A7gJ8CZwJXcJQsfisi0lzxekLpsGDBgi4LFy7sUlxcXAywb9++nLKysvyhQ4ce6tu376HJkycfWa2ooKDg4Pjx4/cDFBUV7T/rrLN25eTkMGbMmH133313v1htVFZW5l566aVD1q9fn29mfvjwYQP4xz/+0enb3/72VoAJEybsLyoq2gewcOHCjhMnTtzdr1+/aoDPfOYzlWvXrs1v7mtNJGG1d/dXzczCMwZnmNkiQklMREQyyN25/vrrt9x8883bI4+vWbOmbYcOHWojj0UOy+Xk5JCfn+8Aubm51NTUxOyITJ8+veCMM87Y/corr7y7Zs2atmedddaR9exiLR4exKLiiUy6OGBmOcC/zexaM7sISPhin4iIpFbXrl1r9u7dmwNQUlKy67HHHutZVVWVA/Dee++1+eCDD5KanxDLrl27cvv3738I4De/+U3PuuOnnXbanscff7w7wFtvvZW/du3aDgCf+MQn9i5evLhzRUVF7sGDB+3ZZ5/tloo4EklY1xNajunbwFjgi8CXU9G4iIg0XZ8+fWrGjh27p7CwcMRLL73U5XOf+1zlKaecckJRUVHxRRddNGznzp25qWxv+vTpFTNmzOg/ZsyYEyJnG950001b9+7dm1tUVFT8wx/+sM+JJ564F2DQoEGHp0+fvnnixInDTzvttKJRo0ZFnczRVOYe3Je1zexcYDahJZ0ecPd7o9T5JPAzoA2w3d3PiHfOcePG+ZIlSwKIVkSkITNb6u7jIo8tX758/ejRo7fHeo4kb/ny5T1Hjx49OFpZzB6Wmc1o7MTx6phZLvAroAQoBi4zs+J6dY4Bfg2c7+4jgM811qaIiLRO8cY5rzSzXXHKDZgKzIhRPh4od/d1AGY2D7gAWB1R5/PAfHd/H8Ddo28NLCIigZs9e3aP++67r3fksVNOOWXPY4899n6mYooUL2HdDzS2Fe/9ccoKgMjpnpuACfXqFAFtzOxv4bZmu/uj9U9kZtOAaQADBzZ7dQ8REYniuuuu23HdddftyHQcscRMWO7+g2aeO9qcxvoXzPIITeSYDLQH3jCzxe6+tl4sc4G5ELqG1cy4REQkC6V06mM9m4ABEY/7A5uj1Nnu7nuBvWa2EBgNrEVERCRCkFs7vwUUmtkQM2tL6HrXc/Xq/BE43czyzKwDoSHD0gBjEhGRLJXIjsPd3b2yqSd292ozuxZ4idC09gfdfZWZXR0un+PupWa2AFgB1BKa+v5OU9sSEZGjXyI9rH+a2VNmdp41ca0Nd3/B3YvcfZi73xM+Nsfd50TU+Ym7F7v7SHf/WRPjFxFpdTK9H1ZBQcGJW7ZsCfKSUlSJJKwiQhMeLgfKzeyHZlYUbFgiIhJLS9wPKx1tNJohPbQUxivAK2Z2JvA4cI2ZLQducfc3Ao1QRKQl+8M3B7B1dUr3w+LY4n1c+KsWsx9WRUVF7sUXXzy0srKyzcknn7y3boWk+m388Y9/LJ81a9axf/3rX7uamd98881bvv71r3/0/PPPd54xY0a/bt26Va9bty5/woQJux977LH3c3ObtoJUIvth9TCz68xsCXAT8C2gJ3Aj8LsmtSYiIs2W7v2wbrnlln6TJk3aU1pauvr888/fuWXLlrZ1ZZFtvP766x1WrlzZvrS0dNWrr7669vbbb++/YcOGNgArV67sOHv27I1r1qxZtX79+naPPvpokxfETWQM8g3gMeBCd98UcXyJmc2J8RwRkdYhTk8oHdKxH9bixYs7z58/vxxg6tSpVVddddWRsb/INhYtWtT5kksuqczLy2PAgAHVEyZM2POPf/yjQ9euXWtPPPHEvcXFxYcALrnkkspFixZ1uuKKKz5qymtNJGF9392fjDxgZp9z96fc/UdNaUxERFIrHfth1dWPJrKNeIup15+zl8x+WYlMurglyrFbm9ySiIikRLr3w5o4ceLuBx98sAfAk08+2WXXrl1RLz6dccYZu59++unu1dXVbN68Oe/NN9/sdPrpp++F0JBgWVlZ25qaGp5++unup59++u6mxhHzRZlZCXAeUGBmP48o6gJUN7UhERFJjcj9sM4666yquv2wINTjeeKJJ97Ly8tL2TJ299577+aLL754aHFx8fBJkybt6du376Fo9S6//PKdr7/+eqfhw4ePMDP/wQ9+sGngwIHVK1as4KSTTtpz44039i8rK2s/YcKE3ZdffvnOpsYRcz8sMxsNnATcCdweUbQbeM3dmzT2mCraD0tE0kn7YTXf888/33nWrFm9X3vttfLG6sbbDyve4rfLgeVm9oS7q0clIiIZFW9I8El3vwT4l5lFdsOM0NezRgUenYiIpE1Q+2FNmTJl95QpU5p8zaq+eBfmrqtrq7mNiIhIy9fS98OKOUvQ3beE724HNrr7BqAdoe0/6m8TIiIiEqhEprUvBPLNrAB4FbgCeDjIoEREROpLJGGZu+8DPgP8wt0vAoqDDUtEROTjEkpYZjYJ+ALw5/CxtC8rLyIiIZneXiRTEklY1xFa2eLZ8AaMQ4HXgg1LRERiaYnbi6RDownL3Re6+/l16wa6+zp3/3bwoYmISDSR24tcddVV/W+77bbeI0eOHF5UVFT8ne98px+E1hIcOnToiC9+8YsDR4wYUbxgwYJOQ4YMGXHppZcOKiwsHHH++ecP+cMf/tB5zJgxJwwaNGjka6+9FnOLlBtuuKHfhRdeOGTixIlFgwYNGjlr1qyeAFVVVTmTJk0qKi4uHl5UVFT8+OOPHwOwa9eunE9+8pPHHX/88cWFhYUj7r///m4A11xzTcGwYcNGFBUVFU+bNq1/U193o0N74c0abwIGR9Z397Oa2piIyNHmtv+9bUD5R+Up3Q/ruG7H7bvrv+6KuQr8rFmzNk2ZMqV9WVnZ6vnz53d56qmnuq1YsaLU3fnUpz513Isvvthp6NChh9avX59///33r3/88cffX7NmTduNGzfm//73v183duzYDaNGjRpet73I7373u2Puueeevmeeeea7sdosLS1tv3Tp0tLdu3fnnnzyycUXX3xxVUFBweE///nP5d27d6/dsmVL3oQJE074/Oc/v3P+/Pld+vTpc/hvf/tbOYR6hB9++GHuCy+80G3dunXv5OTksH379qZthkViQ4JPAf8Cvg/cHHETEZEMi9xeZMSIEcXvvvtufllZWT58fOsP+M/2Irm5uQ22F9m0aVO7eO2UlJTs7NSpk/ft27d60qRJuxYtWtSxtrbWrr/++v5FRUXFZ555ZtHWrVvbbtq0KW/MmDH7Fy1a1OUb3/hGwYIFCzr16NGjpnv37jXt2rWrnTp16qBHHnnkmE6dOtXGay+aRCZPVLv7fU09sYhIaxCvJ5QO6dpeJNr2IL/5zW+679ixI2/lypWl7dq184KCghP379+fM2rUqINvv/326meeeabr9773vYK//OUvu2bOnLll2bJlpc8991yXefPmdbvvvvuOXbx48dqmvNZEelh/MrNrzKyvmXWvuzWlERERSZ10by8C8OKLLx6zb98+q6ioyF28eHHn0047bW9VVVVuz549D7dr187/9Kc/dd68eXNbgPXr17fp3Llz7TXXXFN5/fXXf7hs2bIOVVVVOZWVlbmXXnpp1Zw5czaWlpY2eRg1kRf15fC/kcOADgxtamMiItJ86d5eBODkk0/eO3ny5MLNmze3vemmm7YMHjz48JVXXllZUlJy3MiRI4ePGDFi35AhQw4ALF26tP2tt97aPycnh7y8PP/1r3+9YefOnblTpkw57uDBgwZw9913N7lnGnN7kZZK24uISDppe5HQLMFOnTrV3HnnnR8G3Va87UUaHRI0sw5m9n0zmxt+XGhmWhBXRETSKpEhwYeApcCp4cebCM0cfD6ooEREJP2C2l4kVRJJWMPc/VIzuwzA3fdb/ekiIiKS9bJ2e5EIh8ysPaGJFpjZMOBgoFGJiLRstbW1tfrDPcXC72nM72clkrBmAAuAAWb2BKEtRqYn0riZnWtma8ys3MxuiVPvFDOrMbPPJnJeEZEMe2fbtm1dlbRSp7a21rZt29YVeCdWnUaHBN39ZTNbCkwEDLjO3RudHWNmucCvgP9D6LrXW2b2nLuvjlLvR8BLjZ1TRKQlqK6uvrKiouKBioqKkST2h780rhZ4p7q6+spYFRJZS/BVd5/Mf7YWiTwWz3ig3N3XhZ8zD7gAWF2v3reAZ4BTGotFRKQlGDt27Fbg/EzH0drETFhmlg90AHqaWTdCvSuALkC/BM5dAER+MWwTMKFeGwXARcBZxElYZjYNmAYwcODABJoWEZGjTbwe1lXA9YSS01L+k7B2ERrqa0y0sd3631L+GTDd3WviTTx097nAXAh9cTiBtkVE5CgTM2G5+2xgtpl9y91/kcS5NwEDIh73BzbXqzMOmBdOVj2B88ys2t3/kER7IiJyFEtk0sUvzOxUGu6H9WgjT30LKDSzIcAHwFTg8/XOPaTuvpk9DDyvZCUiItEkMuniMWAYsAyo22PZgbgJy92rzexaQrP/coEH3X2VmV0dLp/TnMBFRKR1SWSli3FAsSexSq67vwC8UO9Y1ETl7l9p6vlFRKT1SOT7A+8AfYIOREREJJ5Eelg9gdVm9iYRSzK5u76DICIiaZNIwpoRdBAiIiKNSWSW4N/TEYiIiEg88Va62E3DL/pC6AvB7u5dAotKRESknnhfHO6czkBERETi0SrDIiKSFZSwREQkKyhhiYhIVlDCEhGRrKCEJSIiWUEJS0REsoISloiIZAUlLBERyQpKWCIikhWUsEREJCsoYYmISFZQwhIRkayghCUiIllBCUtERLKCEpaIiGQFJSwREckKSlgiIpIVlLBERCQrKGGJiEhWUMISEZGsoIQlIiJZIdCEZWbnmtkaMys3s1uilH/BzFaEb6+b2egg4xERkewVWMIys1zgV0AJUAxcZmbF9aq9B5zh7qOAu4C5QcUjIiLZLcge1nig3N3XufshYB5wQWQFd3/d3T8KP1wM9A8wHhERyWJBJqwCYGPE403hY7F8DXgxWoGZTTOzJWa2ZNu2bSkMUUREskWQCcuiHPOoFc3OJJSwpkcrd/e57j7O3cf16tUrhSGKiEi2yAvw3JuAARGP+wOb61cys1HAA0CJu+8IMB4REcliQfaw3gIKzWyImbUFpgLPRVYws4HAfOByd18bYCwiIpLlAuthuXu1mV0LvATkAg+6+yozuzpcPge4HegB/NrMAKrdfVxQMYmISPYy96iXlVqscePG+ZIlSzIdhoi0Ema2VH9Itwxa6UJERLKCEpaIiGQFJSwREckKSlgiIpIVlLBERCQrKGGJiEhWCHKli+zz0QZYcCtUH4hebjnwiZth4IT0xiUiIkpYH7PhdVjzZ+hzIuS2a1j+wRLoM1IJS0QkA5SwItUeDv172TzoGmWnk7t7Q21NemMSERFA17A+riacsHJi5HHLBa9NXzwiInKEElak2urQvzltopdbjhKWiEiGKGFFqktYuTF6WDk5GhIUEckQJaxIR4YEY/WwNCQoIpIprW/SxZNfgrUvRy+rm3SRG29IUD0sEZFMaH0Ja/O/oNsgKDw7enmPYbETVk6uhgRFRDKk9SUsd+g3Bs6+q+nP1ZCgiEjGtMKEVRuaPJGMoGYJ7tkGr90N1Qdjtzt+GvQ7KfVti4hkidaZsCzJhNWcWYLby+HDd6KXvf8GLH0YuvSPnkx3boR2nZWwRKRVU8JqCsuFqo1Q9ufo5V0KYieVp78CFStjn7tNR/jWUmiT37Bs5vFwaG+TwwVg7Uvw3sLoZe6w58PYaycCTLgahpze9HZ3boStpR8/ZlavksV4aNBtMLTpEKcBbzwGT6BOtp0nrbEk0lS64knRa2p/DHTu03g9aZGUsJqiQ3fY8L+hWzS5beHWDyCvbcOyg7uh6FyYfEf053bsGT1ZAbTtAO/9HZ65Mnp5flc47YYoCYHQYr4frYc27WOcuyN07BW97MN3oOz5UDKFeuePuB/t+MGq6OcUyaSxV8B//yzTUUiSlLCa4ovPwM73o5eteBLe+CVU74+esKoPhZJS7+Kmt3v8ebDmBfhgacOyynWhf996IPbzJ14D5/7fpre79uVQoqzvY3/JeozjQO8RcGxxw3rR6kaWH9wNOzfQoAdWX7QE3bBSAlWy7TwtKZZUnicV52ikTrfBzY9DMkYJqynadwvdotn4ZujfWBMnag6FemDJOOee0C2a6oOw6tk4w3oW6tklo+js0E1EpAVQwkqVvPB2JDET1uHkE1Zj7Y6emvrzioi0MK0wYXlACSt8/WnuGdFXez9YFfsLySIi0qhWmLAC6mENOys0oy7ebsWjP5/6dgN0F4YAAAj1SURBVEVEWolWmrBScAG4vo49oeRHqT+viIgArXG19qB6WCIiEqhAP7nN7FwzW2Nm5WZ2S5RyM7Ofh8tXmNmYIOMBlLBERLJUYJ/cZpYL/AooAYqBy8ys/peQSoDC8G0acF9Q8RzhtaEVK0REJKsEeQ1rPFDu7usAzGwecAGwOqLOBcCj7u7AYjM7xsz6uvuWVAXxxgPfoe/mV448HlxbzeP/fJ9Hlkf5QqyIHNU+Paov13+qKNNhSJKCTFgFwMaIx5uACQnUKQA+lrDMbBqhHhgDBw5sUhA5XfpSuXPIkcfbOY73u02mML9Tk84jItnv2M4xlj+TrBBkwoo2Fa/+ejyJ1MHd5wJzAcaNG9ekVTsnXPJd4LsfOzauKScQEZEWIcjZB5uAARGP+wObk6gjIiISaMJ6Cyg0syFm1haYCjxXr85zwJfCswUnAlWpvH4lIiJHj8CGBN292syuBV4CcoEH3X2VmV0dLp8DvACcB5QD+4ArgopHRESyW6ArXbj7C4SSUuSxORH3HfhmkDGIiMjRQd+gFRGRrKCEJSIiWUEJS0REsoISloiIZAULzXvIHma2DdjQxKf1BLYHEE62tN8SYsh0+y0hhtbefkuIIZn2B7l7ryCCkabJuoSVDDNb4u4ZW+Ai0+23hBgy3X5LiKG1t98SYsh0+9I8GhIUEZGsoIQlIiJZobUkrLmtvH3IfAyZbh8yH0Nrbx8yH0Om25dmaBXXsEREJPu1lh6WiIhkOSUsERHJCkd9wjKzc81sjZmVm9ktKTrnADN7zcxKzWyVmV0XPj7DzD4ws2Xh23kRz7k1HMMaMzsn4vhYM1sZLvu5mUXb1DJWHOvDz11mZkvCx7qb2Stm9u/wv92CiMHMjo94ncvMbJeZXR/0e2BmD5rZVjN7J+JYyl6zmbUzs9+Hj//TzAYn0P5PzKzMzFaY2bNmdkz4+GAz2x/xXswJqP2UveeNtR8nht9HtL/ezJYF+B7E+v1L28+BZIi7H7U3QtuavAsMBdoCy4HiFJy3LzAmfL8zsBYoBmYAN0WpXxxuux0wJBxTbrjsTWASod2XXwRKmhDHeqBnvWM/Bm4J378F+FGQMUS8zxXAoKDfA+ATwBjgnSBeM3ANMCd8fyrw+wTaPxvIC9//UUT7gyPr1TtPKttP2XveWPuxYqhXPgu4PcD3INbvX9p+DnTLzO1o72GNB8rdfZ27HwLmARc096TuvsXd3w7f3w2UAgVxnnIBMM/dD7r7e4T2/xpvZn2BLu7+hod+Mx4FLmxmeBcAj4TvPxJxviBjmAy86+7xViBJSfvuvhCojHLuVL3myHM9DUyO7PFFa9/dX3b36vDDxYR2zo4p1e3HkfLX31gM4bqXAP8TL7Bmvgexfv/S9nMgmXG0J6wCYGPE403ETyxNFh4qOBn4Z/jQteGhoQcjhiRixVEQvp9sfA68bGZLzWxa+FhvD+/aHP732IBjgNBfoJEfUOl8DyC1r/nIc8JJqAro0YRYvkroL/U6Q8zsX2b2dzM7PaKNVLefqve8ua//dOBDd/93xLHA3oN6v38t6edAAnC0J6xofxGlbB6/mXUCngGud/ddwH3AMOAkYAuhoZF4cTQ3vv9y9zFACfBNM/tEvHCDiMHM2gLnA0+FD6X7PYgbXhJtJh2PmX0PqAaeCB/aAgx095OBG4DfmVmXANpP5Xve3P+Py/j4Hy+BvQdRfv9iycT7IAE42hPWJmBAxOP+wOZUnNjM2hD6ZXnC3ecDuPuH7l7j7rXA/YSGJOPFsYmPDx81KT533xz+dyvwbLi9D8NDHXXDLluDjIFQsnzb3T8Mx5LW9yAsla/5yHPMLA/oSgJDcGb2ZWAK8IXw8BLhIagd4ftLCV07KUp1+yl+z5N6/RH1PwP8PiK2QN6DaL9/tICfAwnW0Z6w3gIKzWxIuCcwFXiuuScNj2X/Fih19/8XcbxvRLWLgLpZVM8BU8Mzj4YAhcCb4WGL3WY2MXzOLwF/TDCGjmbWue4+oQv/74Tb+nK42pcjzpfyGMI+9hd1Ot+DCKl8zZHn+izw17oEFIuZnQtMB853930Rx3uZWW74/tBw++sCaD+V73mT24/wKaDM3Y8MswXxHsT6/SPDPweSBs2dtdHSb8B5hGYRvQt8L0XnPI3Q8MAKYFn4dh7wGLAyfPw5oG/Ec74XjmENEbPggHGEPmDeBX5JePWRBGIYSmjm03JgVd1rIzTO/irw7/C/3QOMoQOwA+gacSzQ94BQctwCHCb0V/DXUvmagXxCw5vlhGaQDU2g/XJC1zvqfhbqZpddHP6/WQ68Dfx3QO2n7D1vrP1YMYSPPwxcXa9uEO9BrN+/tP0c6JaZm5ZmEhGRrHC0DwmKiMhRQglLRESyghKWiIhkBSUsERHJCkpYIiKSFZSwREQkKyhhSYtiZj3sP1tRVNh/ts3YY2a/DqC9h83sPTO7Ovy4l4W2k/hXxLp39Z/zk3BsN6U6HhGJLS/TAYhE8tAyPidBaJ8nYI+7zwy42Zvd/enw/cmEVmv4cqzK7n6zme0NOCYRqUc9LMkKZvZJM3s+fH+GmT1iZi9baLPAz5jZjy20Ed+C8DpzdZvz/d1Cq9m/VG8Jo2htnERoT6Xzwr26juEe2Dvhc38n+FcqIrEoYUm2GgZ8mtC+RY8Dr7n7icB+4NPhpPUL4LPuPhZ4ELgn3gndfRlwO6HN+k4CTgAK3H1k+NwPBfZqRKRRGhKUbPWiux82s5WEdjxeED6+ktAut8cDI4FXQuuakkto/bumWAcMNbNfAH8GXk5B3CKSJCUsyVYHAdy91swO+38Wxawl9HNtwCp3n5RsA+7+kZmNBs4BvkloJ92vNi9sEUmWhgTlaLUG6GVmkyC0f5KZjWjKCcysJ5Dj7s8AtwFjUh+miCRKPSw5Krn7ITP7LPBzM+tK6Gf9Z4S2ukhUAfCQmdX9YXdrisMUkSbQ9iLSqpnZw8DzEdPaE33eDNIz5V5EwjQkKK1dFXBX3ReHE2FmPwG+COi7WCJppB6WiIhkBfWwREQkKyhhiYhIVlDCEhGRrKCEJSIiWeH/Ay1zKoDJ0okUAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "with pt.Environment(wl=np.mean(wavelengths), t=time):\n", " detected = nw(source=1)\n", " nw.plot(detected);" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Simulate Frequency Domain" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "slideshow": { "slide_type": "-" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAawAAAEGCAYAAAA5T6EkAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOx9eZwdVZn2c2q59/ae7qSzh+wdshAkIYQgyCZCEAFF2Ub4xI8Bx0FBBVFnQGTEYb6BYQCVTSEYUEdZlDXAIBogBEgweyehEzrppLN0p5Ne71LL+f44dWq5W99ObtXp3K7n98svfe+tvvX2qar3Oc/7vuc9hFKKECFChAgRYrBDEm1AiBAhQoQIUQhCwgoRIkSIEEcFQsIKESJEiBBHBULCChEiRIgQRwVCwgoRIkSIEEcFFNEGDBQjRoygkyZNEm1GiBAhhghWr17dTimtF21HiKOQsCZNmoRVq1aJNiNEiBBDBISQHaJtCMEQhgRDhAgRIsRRgZCwQoQIESLEUYGQsEKECBEixFGBkLBChAgRIsRRgZCwQoQIESLEUQHfCIsQ8jghZD8hZEOOzwkh5AFCSBMhZB0hZJ5ftoQIESJEiKMffiqsJQDOy/P5YgDTrX/XAXjIR1tChAgRIsRRDt8Ii1K6HEBHnkMuAvAbyrASwDBCyBi/7BGGHe8B8UOirQhRCHauBPauF22FF21bgAPbRFuRiS3LgEMtoq0IMcQgMoc1DoD7jt9lvZcBQsh1hJBVhJBVbW1tgRhXFCS7gSfOA575umhLAEMDHjoVWPM70ZYwtH8MPHH+4HLGj58LPHyqaCu8+MVJwIODLFpuGsDvLgN+fY5oSzKx+kngv2azZy9EyUEkYZEs72XdTZJS+iil9ERK6Yn19UdRh5TO3ez//Y1i7QCAvg5g33rgT98QbQnDjnfZvxUPiLYkxEDRZwVOuveItSMbXvw20LUL6GoVbUkIHyCSsHYBmOB6PR5Aad1lXRZhRSvF2gEAqR7RFnhh6qIt8MI0RFuQH4NpZ/CevaIt6B+hwipJiCSsFwBcbVULngygk1I6CKdsRwAtzv4fDM4w1SvaAi/0pPVDNqEtAIlO0RZkwjSdn+3xGgRIdIm2IA+s+ykkrJKEb81vCSG/A3AGgBGEkF0AfgxABQBK6cMAXgFwPoAmAH0ArvHLFmHQE9b/g8DZaH2iLfCCEygZJITFrxXA1MxgsMt9zfQ4oMbE2eKGHhdtQW4Qwq5fSFglCd8Ii1J6RT+fUwD/7Nf5Dxe6qWNn107s7N4J3dRRrpRjcs1kjK4YDTJQJ8aJ6jCdX9JIoulgE/b17YNmahgeG45JNZMwomzEwL/sCEOClFLs6t6Flu4WdGldKFfKMaZiDKbUTIEsyYdvj546bHtae1vR3NmMbq0bMTmGkeUjMb12OlRJHfgXugnL0AAlMuCv0AwN2zq3YU/PHmimhpgSwzFVx2BC1YTDGyP3REdLAGUD/wpKKfb07kFzVzN6tV4QENSX12NqzVRURg4zVF2ECRilFG3xNnzS+Qm6Ukyx1UZrMXXYVNTGao/gm61nbbCFwEMUBUfd9iJ+oaW7BY9veBxv7HgDncnM8NDE6om4aOpFuHLmlahQKwr7Uu4EB0hYH+z5AEsbl+Kd3e9Az5LrmVE7A19u+DIunnYxYkqBs+7U4SmsPT178HTj03jlk1fQFs+s0KxUK3HWMWfhqzO/ipnDZxb+xdyeASq/vb17sXTTUrz6yatZ7YnJMZwy9hRcMfMKnDzm5MK/2O2E9UTBhEUpxft738fvGn+H9/a8h3gW9VEXq8MZE87AVTOvwrTaaYXbZLhtGpiqae1pxVONT+HNHW+itTczNSwRCXNGzMHF0y7GF6Z8ofD7CHBC3YeBvb178dvG3+L1Ha9jd8/urMdMqZmCxZMX47IZlw2cvKgVfh8MUY0QRceQJyyTmliycQke/OhByJKMcyaeg1PGnoKJ1RMRlaPoSnVh68GteKvlLTzw9wfwVONT+PdT/x2njDul/y+3FVZhs+uuVBfuWHEH3tjxBkaUjcAVx16BE0aegHGV46BICtrj7Wg80Ig3dryBu96/C0s3LcXPTvsZjq8/vv8vNwamZCilWLppKX6+5ufQDA1nTDgDi8YuwrRh01AVqUJcj2NH1w58uPdDvNb8Gl7c9iIum3EZbllwCyJyAc6eO+MCnZ9JTTyx4Qk8tPYhGNTA6eNPxyljT7HtSegJ7O7ZjdX7VuP1Ha/jLy1/wZkTzsQdp9yBulhd/yfwEFZhzm5/337c+d6d+Nuuv6EuVocLp16IE0ediHGV4xBTYujVevFJ5ydYuWclln2yDM9//Dy+3PBl3HzizShXywdmk5bIfZz7V0wdj61/DI+ufRQgwKljT8U1c67B1GFTUROtgUlN7Ovdh/Xt6/FWy1u487078ev1v8bti27HKWMLuKfT7SoQuqnj0XWP4lfrfwVKKRaOXYirZl2FacOmYVh0GADgQPwAGjsa8W7ru/jFml/gyY1P4oYTbsCVx1458OiGoQ3YxhBHASilR9W/+fPn02LBNE1654o76Zwlc+h33voO3de7L+/x6/avoxf/6WI698m59OVtL/d/guX3UvrjakrvP6HfQzviHfSi5y+in3ryU/TRtY/ShJ7Ia/eK3Svouc+cS+cvnU//1vK3/m1Z8ztmy4+r+z3UNE16x4o76Jwlc+g//+8/093du/Me35nspP/+/r/TOUvm0KtfuZr2pnr7t+e565ktv/liv4fqhk5vXX4rnbNkDr3pLzf1a09ST9Jfr/81nb90Pj33mXPpzq6d/duzY6UzPgf7P37boW307D+cTU9ceiJdsmFJ3utFKaUH4wfpf3zwH/S4JcfRy1+8nB6MH+zfpn2Njk0tq/o9vE/ro9e/fj2ds2QOvXX5rXRPz568x/P76AvPf4Eet+Q4+vSmp/u3iVJKP3is4HuJUkp7U7302teupXOWzKHf/9v3+71+lFLadLCJXv/G9fbvaIZWmG0/rmF2vftgYccXAACr6CDwfeE/OrQJ6/H1j9M5S+bQ/1r1X9Q0zYJ+pzfVS7/26tfo8U8eT/++7+/5D/7Lz9jD88C8vIeljBS98uUr6fyl8+nK1pWFmk874h300hcvpQueWkAbDzTmP3j1bwp2Mr9c80s6Z8kcet+q+woeF0opfWX7K3Tuk3PpN//3m9QwjfwH//EaZsuSC/r93rvfv5vOWTKHPrL2kQHZs27/OnrKb0+hFzx3AT2UOJT/4G1/dcan7eO8hx6IH6Dn/PEc+pnff6b/cU/DX3b8hc77zTz61Ze/SlN6Kv/BrWscmz55J++hhmnQG968gR635Dj6zJZnBmRTb6qXfuvNb9E5S+bQF7e92P8vrPi5Y5eh5z1UN3R63evX0blPzqXPbn12QHaZpkkfWfsInbNkDv3R2z8q7NrfMYzZtfzeAZ0rH0LCGjz/hmy39i0dW/DA3x/AORPPwU3zbio45FCuluPBsx7E6IrR+MHbP0BfvhwMz2H1U9b+xIYnsK5tHX766Z9i4ZiFhf4JqI3V4udn/RxVahV++PYPoeULg7hDgjT3mp4N7Rvw0JqH8IUpX8CN824cUChm8eTFuHXBrVi+azl+t7mfjhrcHiP/eqx3dr+Dpxqfwj/M/AdcN/e6AdlzXP1xePCsB7Grexf+34f/rzB7AG8BRhoopbjt3dvQkejALz/7Sxxbd2zB9gDAmcecibtOvQtr2tbggb/3s2jaHXoz84e4lm5air+2/BXfX/B9XNJwyYBsKlfLce/p92LB6AW4Y8Ud2Nm1M/8vuMO4ecYKAB5d/yhWtK7AbSffhi9N/9KA7CKE4Lq51+Ebx38DL2x7Ac99/Fwhv8X+C0OCJYkhS1j3fXQfKtQK3H7y7QOOj1dGKvFvn/437O7Zjd9u/m3uA3nBRJ6HZ3/ffjyy9hF8buLncN7kfL2Cs6O+vB63L7odTYea8D9b/qd/W/LYQynFT1f+FCPKRuBHC3808LwBgCuOvQKfHvdpPPj3B3EokaeHIrchjyPWTR0/e/9nmFozFd+Z/50B2wIA80bNw9eP+zpe2PYCPtjzQe4D3Y43jxNevms5lu9ajm+d8C3MHj77sGw6b/J5uGT6JVi6aSk+PvhxHptchJXnHmrra8Mv1vwCZ4w/A/8w8x8OyyZVVnH3aXdDlVTcufJO0DyTGo8teexq6W7Bo+sexeLJi3HJ9IGRqBvfmPsNLBy9EPesugcH4gdyH0ipU3QxwJxtiKMDQ5KwNrZvxLu738XX53wdw2LDDus7FoxegDPGn4EnNjyBRC4HV4BTfmLDEzCocdgOGQBOn3A6FoxegMc3PI6kkSMh7nYsOex5e/fb2HhgI74979uHXfJMCMHN829Gn9aH32z6Te4DbYWVe2xe3v4yWrpbcOO8GxGVo4dlDwBcN/c6jCofhV+u/WX/9gA5u3BQSvHfH/03JtdMxpUzrzxsewDgpnk3oUwpw8NrH85jU2GE9fDah6GZGm5ZcMthTTI4RpaPxA0n3ID397yPD/d+mPtAszDCevCjB6EQBTefePMR2SVLMn508o+Q0BN4bP1jeexyRTJCwipJDEnC+uPWP6JMKcNXGr5yRN9z9eyr0ZXqwmvNr2U/wFZY2R+ehJ7An5v+jHMnnYvxVeOPyJZr51yLtngb3mp5K4ct/TuZpzY9hZHlI/H5yZ8/Ilum1U7DmRPOxLMfP5s7TMnfz6P2lmxcghm1M3DGhDOOyJ6oHMXXZn8Nq/etxsYDG3PY078CXblnJZoONeH/zvm/h7fWy4VhsWG4dMaleGPHG9jVvSv7Qe41ajkmGZ3JTvx5259x0dSLcEz1MUdkEwB8ueHLqC+rx+MbHs99UAGTn9aeVixrXoYrZl6BkeUjj9iuKTVTcP6U8/Hcx89lXXaSYUsYEixJDDnCShkpLGtehs9N/ByqIlVH9F0njjoRk6on4c/b/pz9AP4A5cjTvNXyFrq1bnxx+hePyA4AWDhmIUZXjMafmv6U/QC3DVkURHu8He/vfR8XT7sYqnxkzhgAvjLjK+hIdOQmUB7uyuHwNndsRtOhJlw649Ijmp1zXDjtQkSkCF5oeiH7Ae4xyWHT803PY1h0GBZPXnzE9gAsfAoAL25/sX+bcjjgF7e9iKSRtL/rSBGVo7ik4RKsaF2Bvb05egZ6lEx2u57Z+gwA4PIZlxfFLgC4etbViOtxvLT9pewHeEKVocIqRQw5wlq1dxV6tV58btLnjvi7CCE4Z+I5+GjfR9lnfWb+ePqbO99EfVk9Thp90hHbIksyzp98Pla2rkR3KktbGrcNWZzMa82vwaQmzp98/hHbAgCLxixCbbQWb+58M/sB/YQEX/nkFSiSgnMnnVsUe6oj1TjzmDOxrHkZTGpmHuCZnWcSekJP4K8tf8VnJ362sHVmBWB0xWgsGL0AL29/OXvOyG1TjjDl6ztex7F1x2JG3Yyi2AQAF029CBQ0NzH0YxelFMual2HR2EUYWzm2aHbNqJuBhtoGvPrJqznschN8SFiliCFHWH/b9TfE5FhRSAJg+SODGnhn9zuZH/IHKMuMXTd1vNf6Hj497tOQSHEuw2njToNBDby/5/0stuQP47y7+11Mqp6EqcOmFsUWWZJx2vjTcnbrcPJ72R3xu63v4sRRJ6ImWlMUewDgzAlnoiPRgcYDWbZ76UdhvdfKulh8buKRT3Tc+NzEz2FH1w7s7M5SmWfkd8Dt8Xas2b8GZx9zdlFtGl81HjPrZmL5ruXZD+in6KLpUBNauluKbhfAKlHXtq3Nrv4KLAYJcfRiyBHWB3s/wPxR8wfWiiYPjhtxHKoiVdmT1PyhoWZGafvGAxvRlerCp8d9uih2AMDxI49HhVqRnTzzPMyaqWH1vtUDKqkvBKePPx1dqS5saN+QxZ7cCqutrw0fH/y48M4LBWLR2EUgIDnGJ3/47YO9HyAqRzF/1Pyi2sT/xhWtKzI/7Cck+F7re6CgR5zjy4bTxp+GtW1rc0QO8pP733b9DQCbIBQbnxn/GQDsb8+0KwwJljqGFGH1pHqw7dC2wloZFQiJSJhbPxfr2tdlfpjH4azdvxYAy4MVC6qkYt7IeVjbtja/LWmqZmP7RvTpfUVTnRzzRrGdcvPak8WxfLCXlZ8PqBdgAaiL1eHYumPt7/fakz/M9eHeD/Gp+k8VLRzIMaF6AsZXju/fAWexafW+1aiKVKGhtqGoNgHAqeNOhUnN7BOxfoh09b7VmFozFfXlxd9sdfqw6aiL1WHlnpWZH4Y5rJLHkCKsjQc2goJibv3con7v8SOOR9PBJvSkd4jOE1PfcGADRpWPOrzO63kwZ8QcbO/cnrmgOY/C4grohJEnFNWWEWUjMK5yXH7CyuKINx7YiJgcw/Ta6UW1BwDm1s/FpgObMvNYeQid95OcP7q46opj3qh52VVoAcQwb+S8ooWU3Zg1fBYUSck+EctzLxmmgTX71+CEUcW9lzgIIVg4ZuFhEWmIox9DirDWt68HwJx6MTG3fi4oKDYd2OT9oB9VU2w7AGD28NkwqYnGjrQ8jZG7RHpzx2YMjw33ZUY8t35uPwor07FsbN+IGXUzoEjF7808e/hs9Gg92NG1w/tBnpDglo4toKA4bsRxRbcHYOTQFm/Dvt59uW1Ku2adyU40dzXjUyM/5YtNUTmKY2uPzUGkufOhTYea0KP1YN7Ieb7YBbAwfFu8DW19ad36Q4VV8hhShLXt0DaMrhhd1EQ+AFsJbOvc5v0gh8LqTnVjZ/dOzBo+q6h2AMDsEaz7wsb2tPVGHlu85Ln14NYBtxgqFLPqZmF/3/7MXEiOghTDNNDY0XjYXST6Ax+fDEecJy+zpWMLALatiy82WX9rxhqxPJWLWw9uBQDfrhvAWlttaN8AI721WJ6y9i0H2Vj5df0A2M9N5gQxf9HFg39/EK83v+6bXSH8x5AirObOZkyqnlT0760vq0eFWoFPOj/xfpAjdMJn91NrilOR58aIshGojdbik648trgebM3Q0HSoqahl0W5MGTYFAHKPjal7ehvu6tmFuB73zRFPrZkKVVIzWyLlIYctB7egLlZX9PAtx4y6GZCIlKmK84SUuf3ThxU/bMoxs24m4no8c98qI3durelgEyJSpCiLmPPZRUAyCcujkr0dXwzTwJMbn8yu9kMcNRgyhEUpRXOXP4RFCMGUminY3rnd+0GOVjGcsCZWTyy6LQBwTPUxmQ1Mc8w+W3paoJkapg0bwMaCA8DkmskAmLr12pN9ls7HZlLNJF/skSUZx1QdkyUkmDvMtfXgVjTUNhRlAXM2lCllGFsxFs2dzWk2cQdMMmz6+NDHqI5UF6WLRC7wa5cx2TA1gHf6SFMyWw9txZRhU3wJ53KUq+WYWD3RVnMeuwBmWxrBt3S3IGkkfbvPQwSDIUNY7fF29Gg9vjnCyTWT8cmhLA+2/bMz+9vZtRMEBBOqJ/hiy8TqiWjuava+mSMfwtsCTajyx5axFWMRlaNZyDy7PZxoj6nyb4Y+sXpiJmGZBsCLF1xOmFKKHV07MKVmim/25LbJskMtyyCGjw9+jOm1030jUSAPYRkas8lto4Wmg02+qj6OrOPFxyhSnjFe/Fj+N4U4OjFkCIs78MnV/tywk2smY398v7c6L0fVUnNXM8ZUjDmihq75MLF6Ivb3pdni6XTh2MUJ60h7GeaCLMmYVD0pC2FpAF8Ll6awKtXKwnYJPkxMrJmInd07vbkZUwMU7oSd8elMdqJX68W4ynG+2QM4DtjT8cLUAUlhiiEt9NbS3eJLtMCNmmgN6mJ1meFl03AIy3Xt+rQ+7OvbZ4eB/cQx1cegpbvFW+1pE3xFhsLihOVnqDKE/xgyhLWndw8AFLVVjBujK0YDAPb2uVbgmzqc/XmcB2hn105fHxz+3S3dLS5bNIBvy+5WWD27UKaUYXhsuG/2jKsc5+1MQCkbGzWTIHZ2s7HxUzlMqp4EzdTQ2tvqvGnqgBrLsGdXj7+EzjGxeiL69D4cSLi2zzA0Rliy4iGGuB5HR6LDdxIF2ERs+6Esk40shNXaw8ZzfKW/YwWwa5g0ktjft995k0/EIpmEtbN7J6rUKtRGa323LYR/GDKExW9sv2L+o8pHAYC3NNnQXSThOMG9fXt9dTYTKll4jzsQx5ZMJ7O7ezfGVY7zlSBGV4xOIyxrVszHJo3MJ1b5k9vj4OFPD6EbGiBHWFjQyAyZBkFYALx5LFNn6kpSPZOMPT3+Tr7cmFA1IXvRhVph2egiLGsCMKZyjO928UmZJyxo5g4J7uzaiQnVE3y9z0P4jyFDWPt696FKrUI5d5JFxuhyprD29bkIy60iLKesmzoOxA/4suaJg393W9y1TsWtsAyvwvLbGY+qGIUercdZWM3PnxYSpJRiX98+W636BX6tPLPzHORgKyyfVQMnnwyFLsmMSN2TDItAglBYoytGoz3eDi09H6tmhnODtItPajw9GLktahbC6vZ/IhTCfwwZwtrft9/XiqqRFey7PQrLzFQ1B+IHQEFRX+YfYdXF6iARKS1ckj1Rvrd3r+3A/UIGmXO1maY+u1Jd0EzNVzIHXITuXnhq6iz0JqueHN/unt2ojdb6NtHh4PemxyZDY/akhQQ5MQShsEaVjwIFRXtfu/OmmT1y0NrTiqgc9TW8zFFfXg+JSJnPG8Duc9c9blIT+3r3BTJeIfxFSFhFQlSOojZam6mwuIqwEvxc9fhJWIqkYHhseBaF5SVPzdDQlerC8DJ/HYyd3+NhQduxeGfpnGD9JqyYEkN1pNp7rXi+SFI8zq6trw2jKkb5ag8AVKgVqFAr0lSfllX1tfa0IiJFfFsX5ga/dhljlS283LMbYyrGBBJ2y3qPexSWQ6QdiQ7oVPf1+Q8RDELCKiJGVYzKEhL0JvL5DNpvW0aUjciisLyzYp7gF0dYXsXHx8ZPMucYWT4yU2FJqkVYjrNrj7cHohgA9nd7Ccuwii68qm9/fL+tMPwGz816cpBmbrU+psL//BVHfXn6eLmWAaRNOgD/n7kQ/mNIEJZu6mhPtPt+w44sH5lWdOEq3eZOmSssn1VEVoecHp7khOWzQ64vrwcByR0StPJ7fGxGlvnvWEaWj8ySw5ItcnCc3YHEAd8J3W1ThmKQM1XfgXhwNnF16Z2IGc4SABeR+p2bTcfIspGZIVQgY91aUM9cCP8xJAjrUPIQTGr6HkKpi9XhYPKg84YnJMge7P19+0FAfF1nBLCH0xsS1F3rjJx8GgDfx0WVVNREa9CR6HBsAXKGS0eU+x/qGlk+EvvjaQpUVj1rniiljBwCUljZSVS1ii6cSsqORIfv9w9HlVqFcqXcq7AMzRXOZXZRSgO1C8hxjwPWRIg691VfcBOhEP5iSBBWV7ILADAsOszX89RGa3EocchZ/OlZYMkepvZ4O+pidb62rgHYw9mR6HCqu7KUtXPCCmK2XhurdQjLPRN2vd7ftx9VahXKOLH6iPqyehyIH3AWD9vk4BQ4dGvd0EwtMDVTX16Ptr421/2jZw0JdiQ6AiNRQgjqy+vRHncXXWTm1vr0PiSMRKCENbLcusf5/ZTrvrImJkHk/EL4iyFBWJ0p1im82F3a01Ebq0XKTKFPtzpMuLs5WLO/Q8lDqI35v3iRn8Puku7JpzEnHVRIEGBkfjBx0HP+9LFpj7cHRg7Dy4bDoAa6Ul2ODXZXieAJHWAkmjJTjk1ZQoImNXEwcTBQYvBcO4CRp6x6wqcdcTYZCZqwADhkauewvAvk2/vYJFGV1cBsC+EPhgRhcYXlN2FxBec45sxuDl2pLlRHqn21A3D+Vg9hyVGwRqpW0UX8ACrUCsQ4cfiIulidd1yAjPxeV7LLdxXMwc9zKHmIvcHJwaVmbMIKSM1k2JSlEKQr2QWDGoGRKAAMiw3LDHVLsid8yic/QRIWV0x2WJCr0CwKK1RXpQFfCYsQch4hZAshpIkQ8oMsn9cQQl4khKwlhGwkhFzjhx22wor4r7AAy+Hw9kOK1S/Q5ZSDJCzH+bnLth2HHJQzro3VOk7PngnzPIiLzKP+jw2Qixy8aqY9wWbuQZFD5oRHc4UELSWTCF7JeCYbtl3e8Cm3K1AizbiGaQvSeVQjEUxUI4T/8I2wCCEygF8AWAxgFoArCCHpOxb+M4BNlNLjAZwB4F5CSKTYtnCV4bcz9DgcO+zFFRZ7HZRT5rY4CsvIcMhdWjDkCTDC4sUvjsISoz4BphoA5syYDZqjZiwnzJ10UOSQcc146E1SMsK4gYcEkwdZbo1PxGSv8hNhV1aVzAmevwabsPo9WQ0RDPxUWCcBaKKUbqeUpgD8HsBFacdQAFWErTSsBNABQEeR0ZnsBAFBVaSq2F/tgUfV5FgcKzYkaD3MlvPrTnWjMlLpuy0Acy4mNdGj9ThkrmaGBAMjrIzZueEaH3bteCupwEk0j+oTQlixWuimjm6t27mvJW9FpYgcVu5Jh+q8Bruv/E4HhAgGfhLWOACu7qLYZb3nxs8BzATQCmA9gBspde8XwEAIuY4QsooQsqqtrS39437RmexEdbTa94WWnBCZU+bhCUdF6KaOXq03EAeYc/YpyR6H7DeJc/DzdKe6Xb0EHfVpmAa6tW5xIUFPpwuH0KNyFBG56KK/MJvcIUFXDst9bBCwQ90J10RMkj0hwYPJg6hUKwMbK4CV3MtEdl1DVzGI9ZpSis5UZ2CTjhD+wk8Pnq0/C017fS6ANQDGAvgUgJ8TQjLuLErpo5TSEymlJ9bXD3zxX1AhAZuwUj2ZCsvUmbOG/6FJAChXyqFIipNPo66QoOVkhBFWFvVpj01AjqVcKYcqqWnkoHpDpqkuVKrBKFAAqFQroRAl0wFLip3n4+MUlDIGnGvimWykhQS7U92B3UschBDURGuy52kBwNQR1+PQTT1UWCUCPwlrFwD3NrbjwZSUG9cAeI4yNAH4BMCxxTYkqJBAVI5ClVQWOuEVS64EMC9XDsIpE0JQHam2yNMKwaWFcbq17sAcssfpZclhBTk2ABufqkiVTQCeNkhcgWrBETq3qTpa7Q3jSrKHGHq0HihEQUz2v7KTg49BV6orS0jQWrMWYHjZDQ9h2c2CnZAgv+PgkQIAACAASURBVK9CwioN+ElYHwKYTgiZbBVSXA7ghbRjdgI4GwAIIaMAzACQtlvckYOHBINAVaTKq7DsKkHdDucE5ZQr1AorPOkK41ghL83UENfjgTmZrE7PlcMS4Vgq1Uo2PoB3zZNLzQQdSvLYlEX1cWIIcl+nrOFc2dpY0poM9Wg9qFKDVVgAq/zlz5Un5wcAhuYUXIUhwZKAb4RFKdUB3ADgNQCNAP5AKd1ICPkGIeQb1mH/BuAUQsh6AG8CuJVS2p79Gw8f3Vp3YA9TpVrpJSzXOhpbRQREnpVqJXq1XiefZuewNPSmepktAT3IWUOCrn50QZM5wMJq9h5dnjVPTsg0aNXgscmtGFwKq4JvnhgQsqrjLOFlEQqrMlLJIhqAKw/JFZZuE1aosEoDvvYHopS+AuCVtPcedv3cCuBzftoAAHEtHthDXhWpsqqp3CTBHuygw14ZCsvl/PhDHlRIkJ+nO9UNEJm96VafAY8Nt8mrZrzr1LpSXYHsnutGlVrlssmddwy+UIbDS1j8vs4MCU6umRyoXQAbL3vnaNO1fxjAFJbJus6ECqs0MCQ6XfTpfYH0pwNcM2SeN7IbqhpOwjxAkmAKi+ewHIccdPLeQ1jpPd9cIcGg1Ce3yQkJZrYb6tF6Ai26AFyTDCCNRF0hwYBtKlPKIBOZXSMjffLjCgkGTKSApbBSboXlLWsPqstNiGBQ8oRFKUVcjwdGWPYM2ZM3YmG4uB4HAN93r+WoiFQw8uQEwW0xdDvsFFSoVJZkVKgV1iw9s5cgd9KBVuVlhAQVe3IBCMphpduUVo3Xq/UGHnrzFKh4QoIyYGiglLKQYMBECmQbL8VT1s673IQKqzRQ8oSVMlMwqBEYSZSr5UzVGGkhQVO3m+IGpvZsheV2Mt6QYJCz4gqlgpF2+gaOho4+jY1NEH0NOex8I+BSM2xykTJSSBrJwMnBLtoBmJpxVy5SKqy4oSpSZRXMuO9rFhJMGAnoVBeSw6pSq5AyU0gZKZfC4mXtGnpSPZCJHNgzF8JflDxhxTWmaoK6YcuUMq9Tlpy8UVyPQ5VU37cW4cisElSFhQQBRuZ9Wp/j9OQIeDNeroKD2EWXozLCQoLUMABqenN8qeAJHXDClJRSF4nyEJchrHy8Qq1gEy73OixrrHo1VsAjgkj5WNj5Nc/CYQ19eh/KlfJAqypD+IeSJyyuasqVgBSWUp5GWE7SPKEnAp3pVaqV0EwNKYu03WqPO5kgwzhlShm7Hu4iECs/E2TYlqNSrQQFRR/fzoOX/RtOiFIEYVFQZ5xkldkFgBopFhIUEHorV8rZ5C998mOImfxw8LHo0XosRerNYcX1OMrUUF2VCkqesHjeKDCFpZYhaSRh8B1i7RyWHrhT5mHQHp6U5u10LFvcxwQBh7BcRSAu9Rk0YfFJjENYjgLlhFWhBFtCXhFh57Nzj65mrvFUNwxqCCluKFPLvOFl1z5dPIQpgki93WW0jBxWn9YX2GQ1hP8YMoQVWA7LejjiloJxb7seOGFZtiR0ty0OYREQRKTger/ZIUFPEQibpYsgLD7zjvOckUvxJfSE55jAbLLGgN0/1JOT6bY6OohQMuVKuTck6Jps2EskBIUqAThkmpbD6tP7Ap2UhfAXJU9YPJkfZA4LcBGWqyw5aKfMCxgS1hi414Ql9ARiSizQ2L7t9DJyamxsgp4J29eKE5Zr7ylOWEG2QPLaZKlirmQA9CbFqD7AlcPyrMNSPBWnIhSWPV56PEvzWy1UWCWGkieswEOCNmExh2dvdGcawROW5WzjnhwWK9tOGsnAnbFTdJGW3zNZ6CZwhWVfKxdhSQoAapO8sDBlkodxHcJK6D1CbOJ2MXWcGRIUEV522wXAIdO0ThehwiotlDxhBV10wZ1JHw/D2WEvDXEtWMLi4ayE7iYsJ58WZAk5kKMgxdqSXkhIkF+rdDUMRyEHPUaZJOoohkRKjE2AM9mgWZZr2OFTAUTqVViZnS5ChVVaKHnCClph2TksPbMyT5jC0tNCgqYTEgwS5Wp6SFAWNjaAO/yWluMDkOCEJSokmIVEEwGv43OjXCmHTnVohjtyYIVPDTHhUyCNsLJs4Bhkl5sQ/qPkCYvnsIIKC9iJfE5YAosu+Lk8CstKlCeMRPAhQaUcuqlDM1LMFkK8+T1RBQ62GnapGQELmYFsRTtuwmLXMSpHA7UJcJ4f/jyl39cAEFUE22Wkd7pgUY0wJFg6KHnCEpbD0vlM1BuGE0FYcT4rtqvgjMDXhLnt6dPjTiWXq6w96NCNQw6c0GV7ds7JIfBCEDVdYblVH7NJSEjQnSsCPPd10khCkRSoXNkECFVSIRM5q8Kihm4vHA5RGih5wurT+yATObCHyXbKhouwBK014o7NIU/eS5CFcYJ2fHwGnjBSTthGUkB1MQuHHTXsVg1skW5cY/dNUF1JbJvskKDLJq76BIYE+Vj1ugnLVVEpIhwIsD6HdncZnsOSJIBISJnJQNuyhfAfJU9YfOYeVPk2D9ckjSR7w8pBmFbPtSDDXk5I0K32VDtRHnRoiTu1pJmyiQGSgpSpwaCGwCpB7oQdcojrfYGX/QPuvGNaZSeciYeIkKB97XS3WlcBaggp4HHD00HFngip6Au4LVsI/zEkCEtEoUOSd7qwwnAJqxxYyDost9pzVXYFrrBkt8JSbJviJhuroB0L7+sYTx8fsJCgCNUgSzKicjStUIaROycLEeTg3Et8IuYOVfYJU1gAy2PFNa6wnFBz3GS2hiHB0kGw8Q4BOG/SeTi+/vjAzheRWecIm7AsJxinwROW45DTnIwZfF9DwHF6STPlJMZllREWEVcW7RTIOGpGBKG7berLFhI0BoHCsidiik0OiUGgsDw5LACQFPRZBB+GBEsHJU9Yi8YuCvR8dkgwbcfhuPU68LCXXMYUDQBPL0EjeAXhKCwtTWElAFmMY4nJMVf41slhJQTk1GyblJhL9bmUjJUrCrKjPYdz7dLCyxA7VgAn+F6n475lX58RKqxSQ8mHBIOGLLFEfcr0KqyEpbCCnoky5+fNp4kOCabMtJAgFUPm3CZPvpHnsASU/XPE5Bjb3wnwqL64kRRSOg7kDwmKzmHF5JiXSAFAVm3CCnNYpYOcCosQ8kIBv99BKf1a8cwpDUTlKBKmd9+gJDUBEESlYB1ORI5Ac9siKTBMHZqpicth8d19LZsSSd22NWgwwuLj485hJRGLBt8bD2Dj4C3asXJYAtppcTghQfd9zcYqqSdQJ3AL+qgcxcF4u2MXAEiqHeUICat0kC8kOBPAtXk+JwB+UVxzSgNROYqUJyQoQzN1AGrgTjkiR5BM23E4aVW+lcnBPshcHXhyWJKClLXdiIjcTESOuNSws4YnYSRRq9QHbg9gKSzdZZPs2CTK+dpLEswUQCRr0beTWxOpsKJK1Kv8AEBWbMJS+b0W4qhHPsL6F0rp3/L9MiHkJ0W2pyRgkwSR7Qc7RQ2IIKyoHEUqbR1W3CKsoMNL9iyd6p6QoGaFS4Pc6oQjKkeR0Ph+WK5CAiMpzAlH5AgSybQtT2CFBAWQOpC+JMGZbAAQ0jXFjagcRYoTVhaFJWrMQhQfOXNYlNI/9PfLhRwzFBGTY0ilOeUUmIoIuhtARIowWwA775CQiG1nkPCEBN0Ki7KxERUSTKWFTAEgbgpUMxkK3VIyZkoYiXoKZvh1stesiSN3wJp0uEOoACCrSFmRhZCwSgeHVXRBCLmu2IaUEiJyJItTNu3PgrYlmdYdPcFDgoLK2hmZO2Xt3LGICN2w8ckkh7jAfFFUjnpt4rkiI/i8I4cssW4xSdO11olXVApUowBXWK51j4D1zInLjYbwB4dbJRjs8v+jDGyGrDvdHGSHsAJXWHLEJgS7YpGHBAOeeToKy0hTnxaZCwgJxpRYxvgATEmIqshjhOVdeA4whRV03tGNmBxjxUSubhIUFmEJDgk6SzcyJ0KhwiodHBZhUUofKbYhpYSoHEWSGp5YvyYo7BWR3SFBlsNKCSIsRVKgEIWNjewiLIEhQY/CclW+paguzAlnK5QBgLgpjkQBq7jBHTmQVWgATFChlXhRJQqd6izobt9XKrvPECqsUkK/C4cJIbdne59SemfxzSkNROUoutMLC0ABBK8iHPJU7AIQzSIsUSG4hOYlc64+RcyE2fh4c3wmAJ0aQrqPA2mqz6OwNKFKJibHkNQSnvvazocKXocFAElCUO6pEtQhSRIUUvL9EYYMClFYva5/BoDFACb5aNNRj4gcQcKz6l61VU3QJMGKLtJCcNwWAQ45psTYmjTXLJ2HBEXYE5Hcaka2CJ29FFUOze6fLDksqoslBiVmqWNnspGwum6IDLvZ7dAI8T5z1EBEigTewDiEf+h36kEpvdf9mhByD4BCFhUPWUTlqEUSvCO5E4YTExL0EpZIhaVKKjSkESiltq1Bw75WgL29iAZxhG7bZK1Nc3eUSJqG0PBWRI4wNeqqEkxaXDBYFJY7h5VMGYj4pEhXr149UlGUXwGYg7BjULFgAtig6/q18+fP35/tgMPRyuUAphyRWSWOqBxFAqY3AcwJS0hI0HSV+yqw5u5CihwYgXp7vqVgQiJS4HtPcXuS1AAFQKxFuqImFxxROQodJnQAiqs1kyYwTAkwAk9R0xPO1QSqdQ6vwvJOhPxSfoqi/Gr06NEz6+vrD0qSRH05yRCDaZqkra1t1t69e38F4MJsx/Q7MyCErCeErLP+bQSwBcD9hRhACDmPELKFENJECPlBjmPOIISsIYRsJITkXah8tCAiR6C5SUJS7DDToFJYgnaI1UA962U0UCHkCTBlYIJCB+wqQdFO2O65yBWDpIAC0GAKV1gpT0hQ3ETMDXsXgDSFlYKvinROfX19V0hWxYMkSbS+vr4TTLVmRSFT2gtcP+sA9lHKs9S5QQiRwVo3nQNgF4APCSEvUEo3uY4ZBuCXAM6jlO4khIwswJ5BD9spy968kULkwDttR+QIUjBBJYUFutw5LEFFFylQzyw9CbFqBgBSkgxVkiwn7NgqAm7FUC4pgCRBt+4boUpGiqDLEzlQnPCpwPZH9g4JaTmspI8KC4AUklXxYY1pTidZSA5rx2Ge+yQATZTS7QBACPk9gIsAbHIdcyWA5yilO61zZY1bHm1ghOVVWClChBUVmAB0WYUKeKsERSksamaQeUSQI7bJQVZRATAnLFhheXIylgPWLDtFKhlVVlm+0R3OHQQhQZuwJOLtdAEzXINVYjjcThcvFXDYOAAtrte7rPfcaABQSwj5KyFkNSHk6hznu44QsooQsqqtre1wTA4UqszDXt4HW4RTdhSEdakl2QlPCnB+TkjQnQcBIqIq8qwxSMmZeRlR5OAoLMku3ElZBC9SyUQkptadiZironIwEJZHYSlIgpbsGqz29nb57rvvFtOdOQ9OOumkGcuXL8/YgOyBBx4YfvXVVx9zpN9/uPGpfyzgmGy1pOkSWgEwH8DnAZwL4DZCSEPGL1H6KKX0RErpifX1g+4aZUCVVOgAqF0lyJygKqioAHAcHlsTJrBKUFaRArxl7QIVlmorGMcJpwSHufh5ddm5X7h9oosbmMLiVYKK8AIVwDVeQFoOq3QJ68CBA/Kvf/3rglMopmnCMAw/TQoEh+VBKaV7CjhsF4AJrtfjAbRmOaadUtoLoJcQshzA8QC2Ho5dgwXcqeiSzMJwMnOCEYGElbTbRIkNCUakCHoIdZX8q8LCpdweANAGSVGK+7wpF2GlpMFBWCw3O7iqBPm5tfR1WKCoDoCwbnlm7YSte7uLuq1xw+iqvv/88vEtuT7/3ve+N76lpSV67LHHzjr99NO7Ro4cqT3//PN1qVSKfP7znz903333tW7ZsiWyePHi6aecckr36tWrK++5556dN9xww8STTjqp56OPPqqcOXNm39e//vX2O++8c9yBAweUJUuWbD/zzDP7sp3vrbfeKv/ud797TCKRkGKxmLlkyZJPjj/++GRPTw+5/PLLJ2/dujU2ffr0RCKRsIXK/fffP/y+++4bU19fr02dOjURiUSOOOdXSJXgdELIM4SQTYSQ7fxfAd/9IYDphJDJhJAIgMuRuX7rzwBOI4QohJByAAsBNA70jxhssB8g2wmyMFxEwIp7x5bB4WRY0QVcIUEZSUKEjA2Q/VqJLrrIuGZwKSyBIUHVIgFvSHAQKCxO8MSdw2LFPEFvmBoU7r333l0TJkxIbt68edO5557b1dTUFFu3bl1jY2PjpjVr1pS/+uqrlQDQ3Nwcu+aaaw40NjZumjZtWqqlpSX2ve99b//mzZs3btu2Lfb0008PX7Vq1ea77rpr11133TUm1/mOP/74xAcffLC5sbFx049//OPd3//+98cDwD333DOyrKzM3Lp166bbb799z6ZNmyoAYMeOHerdd989dsWKFZvffvvtrVu3bi1K765CvMQTAH4M4D4AZwK4BgU0v6WU6oSQGwC8BkAG8DildCMh5BvW5w9TShsJIcsArANbNPYrSumGw/tTBg/sMFOGihBHWHpaDkshspAOACyHBW8/OkIQFRx+s0OChNjkJUz1Wc5fk2X7PW6TyKKLiBxhZO65duzHQamwSDAdOPIpoSCwbNmy6uXLl1fPmjVrFgD09fVJmzdvjk2ZMiU1ZsyY1Nlnn93Ljx03blzypJNOigNAQ0ND/KyzzuqSJAnz5s3r++lPfzo21zk6Ojrkyy67bHJzc3OMEEI1jV35d955p/Lb3/72fgBYuHBhvKGhoQ8Ali9fXnHyySd3jx07VgeAL33pSx1bt2494lXchXjQMkrpm4QQYlUM3kEIeRuMxPKCUvoKgFfS3ns47fV/AvjPAdg86OHMkJ0cFsvTCCCsXORJ5Dy/5R8iXOGlFaREBSksvlg5JTnjwUNx4hWWKyQoywCo8LJ275IE2c73DQaFpRF4l0uQodH4llKKm266ac8tt9zS7n5/y5YtkfLyctP9njssJ0kSYrEYBQBZlmEYRs4Z7K233jru9NNP737jjTe2bdmyJXLWWWfN4J/lmvj6MSEupOgiQQiRAHxMCLmBEPJFACWxXsov2CEKW9UwJx0RQBK88aebPDUQqAGvB+NQiRVycxc5CAqXAlkmF66fheewstkkskpQjkAnBKZr8sNDgiK6lHDYkzK4Ol1YeeOoQIL3EzU1NUZvb68EAIsXL+5aunTpiM7OTgkAPvnkE3X37t1FvSBdXV3y+PHjUwDwyCOPjODvn3rqqT1PPfVUHQB8+OGHsa1bt5YDwGc+85nelStXVu3du1dOJpPk+eefry2GHYV4rZvA2jF9G6yi76sA/k8xTl6q4A+v5t4Pi0AIYWUqLEVYPg1ghOVVWLKlPgUpvqzhN/azqPCbTaKDMCQIuCsqB0eVoF04k0W5i4hqBIHRo0cb8+fP75k+ffrs1157rforX/lKx4IFC45taGiY9cUvfnHqoUOHivpA3XrrrXvvuOOO8fPmzTvWXW1488037+/t7ZUbGhpm/exnPxt93HHH9QLAxIkTtVtvvbX15JNPnnnqqac2zJ07N2sxx0BRyMLhD60fe8DyVyH6gU0SxFFYKUJQIYKweA6Ly3OrFFmswvIu8BSlPgF3OMkVEhSsZmxicNnE7yXRRRcAkCIyooAnhyV6QTNghQQ9jXnFFfMEgRdffPET9+vbbrsto/HCxx9/vJH/PGPGjJT79bPPPtuc67N0fPazn+1tbm626wvuv//+VgCorKykL730UtYivBtvvPHAjTfeeKDAP6cg5PRahJA7+vvlQo4ZisiZwxJIWFr6mjBBhBXhhOWapScF5tSc8XHGY7CEBD1hSlms6gNca/pcBTypQdAyyg57u4ouKFGQkgiiJaqwhiryXc1rCSFdeT4nYKXqdxTVohKAE6Kw3pBUaBCUw+LhSZfa0wUqmgiRoRMCSnhvQ9UquhBLWJ6iC8EKKxuJpqzxEZkrcu5rN7lLkEEgCwrpAoAsyZBB2HWzIgkpa2NJUff50Yr7779/+EMPPTTK/d6CBQt6li5dulOUTW7ku/sfA1DVz+8/VkRbSgaZeSPZCsOJKSMHYCfHYZWVq4K28OFjoEkSIoDdQFVYSDA9fGv/bIrLYWWzySIvsbmitGIisHtc1L3khgoJuruq0hq7kLAGBj/CeMVETsKilP4kSENKCRkkYcX6IwIe7GzrsESRJ+A4kBSxCEviBSmCQpR2pwtnPJhDNgdVSJArBrHrnZxrx6ERWdi95IZKiGe8+P2uDALbQhQP4qdGJYhMVSOu0CFbAYhGiLBZsWKd1x0uTQm0x7lWXoUlA8LCXLlscn8mAnzClUoj98GhsNIIyxovhYaEVUoQf6eVIGyHAy9hiVBYTg7LrfZEKizL6VmnN4gESsjhNbUsAuzGqZ7wG0Gk/2YuvoHb5CYGTlhCF+haY+Im0pQkCR0rDgXEYxevii3NVVhDF4X0EqwLwpBSQkbSXFJggEAR8GBnsyVFHOcTNDgx6datZ/BQlyB7PH3oLKSIJORacShEAaGuSQYGh8LiY6Kn2SXq2rkRAfEUqfAJSJjBKi0UMuV/nxDyR0LI+URE87mjELxEW+NvSAp0AsgZu6sEYAvPYXGnQiQrJCiIsKwh0C2i4naJciwZCtT6WWRDH0IIVAy+ogt7suEeq0GisFSk5dZ4DkuQPX5D9H5Y48aNO27Pnj2BD28hhNUA4FEAVwFoIoT8LNueVSEcqGl5GpNIMAmxnXWgtqQrLMJCJxFBsX3FIm2D8DGyWvsIsQaQiAQlTc2kCIEqOPfBHLDzmjtjoS2QrGunp+XWVAETsXSo1Fs4o1v/i3jmgsBg3A8riHMU0umCAngDwBuEkDMBPAXgm4SQtQB+QCl9z1cLj0Ko1nPDFZZhJYNFhJlyKQhRgSXuQLg9/PaWqTjPooJmKizBjk5NDwlKBAqlkARVUwIudZxG7oNFYbkVqRGkwvrTP0/A/k1F3Q8LI2f14eJfDJr9sPbu3StfcsklUzo6OtQTTjihl1rPa/o5/vznPzfde++9I//yl7/UEELoLbfcsucf//EfD7700ktVd9xxx9ja2lp9+/btsYULF3YvXbp0pywPLLZSSA5rOCHkRkLIKgA3A/gWgBEAvgfgtwM62xCBEzphw6tbM1ARTjmjYtH6WRVEEHYexHrNnZ/ImTAjB+e1yPHhiCBLrkgwiSqUKyzvvaQIHisgk+B5qHkw2OYHgt4P6wc/+MHYRYsW9TQ2Nm668MILD+3Zs8eOTbvPsWLFivL169eXNTY2bnzzzTe33n777eN37NihAsD69esr7r///pYtW7ZsbG5ujv7mN78ZcEPcQiYg7wFYCuBiSuku1/urCCEP5/idIQ3Z5A82e23YYa/gHx4ZgERphpMRFcZxnB57zUODIqu5VErTii4Gg8KiTpUprDCl4NAbv0aa670UISgbBJygUppWdBFgqDmPEgoCQeyHtXLlyqrnnnuuCQAuv/zyzuuvv96O/bnP8fbbb1ddeumlHYqiYMKECfrChQt73nnnnfKamhrzuOOO6501a1YKAC699NKOt99+u/Kaa645OJC/tZDr+a+U0j+43yCEfIVS+kdK6X8M5GRDBXxQ+RUVGk83NMv5ORCqsNJm6SLVJ4dK00OCEK6wGIk6rzUAEVOsTfzapav1GtPM9SuBQQVFj+s1f+ZEX8cgEMR+WPz4bHCfg+YZ7/SavcOp4SskIP6DLO/9cMBnGkJQrAfYccq8Ek7Aw2NqUOANeaUg7kFWwMbGSCuRFlnNpVBqOziArZ9TqFgnLFNqK3MA1lIECgh0wHYOy/WeyHvJDTbpcF7zn0s1JBj0flgnn3xy9+OPPz4cAP7whz9Ud3V1ZU0+nX766d3PPPNMna7raG1tVT744IPK0047rRdgIcHNmzdHDMPAM888U3faaad1D9SOnH8UIWQxgPMBjCOEPOD6qBreezZEGmTKhke3CwuY8xPy8NgKyzUrBqAKmhUr1ml5SNBRWIIIglIoaeRgDAKFpVDT64Bh2WRogCKmtJ2TuO6ySydAZBAorIhpeqIITlSjNAnLvR/WWWed1cn3wwKY4nn66ac/UZTixXTuvvvu1ksuuWTKrFmzZi5atKhnzJgxqWzHXXXVVYdWrFhROXPmzNmEEPqTn/xk1zHHHKOvW7cOn/rUp3q+973vjd+8eXPZwoULu6+66qpDA7UjHwu3AlgF4EIAq13vdwP4zkBPNJQgmyaIa9auWwQmxCkbKaiU2sRgUhOUQJiC4AqLj42d3xNFWIYGBV6FpQOICVZYTPV5y7QjlAJGUiBhUdsWjhQoVMFjBTDbPIRV4goLCHY/rNGjRxvvvvvux663WgBgzJgxnnNIkoRHHnlkF4Bd6d9RVlZmvvzyy1n3zioU+ZrfrgWwlhDyNKU0VFQDgZGCDMcZGybLZolRWCkWEuRraEx2KRVRCovnQbjzo/x9Ieaw8aHIICxR4wMAMHSLsJxBYcQAprAEIZvCEqnW3VBNw77HAbfCEm9biOIhX0jwD5TSSwH8nRDidicEbHnWXN+tO1phheHSFZYQJ2hozPnxEJxFWLKwkKDl9OxOF6bn/cBhpJjCcoe5AMjU3wWQeWGkIIPaeT4A0AkLXUJPCjNLtYo+vMRAoZgCx8qCSk2vXdaEKCSsgcGv/bAuuOCC7gsuuGDAOat05AsJ3sjPdaQnGXKwFRZ76ZCEAKFqaJCpU7FoUK72xIhm7ty4etAte2RRIt4aH4/C4uRgmkCOyih/bUpBpUxV2TbBWhZhZE0dBALV1GxbOHRQKCLJ3UKGwrLm2IrgysqjDYN9P6ycTyOldI/1YzuAFkrpDgBRAMeD5bdC5IKR8lSe2SFBIQqLKQi7xN4mTzFORrGcnmHderb69LmlS04YyYzwmwErdCmKHOz7x6tkZNEhQYNPNhzooMLUuhuqmV1hDYb8WojioZDp43IAMULIOABvdD3y8QAAIABJREFUArgGwBI/jTrqkTZrd/JGYhSWO0djKyzRIUHeU9Amc3EKSwFgpDk7hYIVOAixieUd9XSbQMXZBICYKchpSwB0SsVdOxdUU0fKRU520cUgsC1E8VAIYRFKaR+ALwF4kFL6RQCz/DXrKIeRtFSN5ZR52EvEw2NqkEHtXJFNnoagkKDBw0rp9ohSWNnVDCMHQWomTaEDXPVBaEgQuqVGiXusTDYJEamyTNMqazfthav8ygm7r0L4goIIixCyCMA/AHjZeq9Uu/YXB3blWVplngiSsGwZPCFBdn6N50PskKDA8BucEBLAnLAsssDB0Cwl41IMMK2iC4GExStOqWuJBMTn1mBq9ro5fj/ZuVGRdvkI0duLiEIhhHUjWGeL5ymlGwkhUwC85a9ZRznS1vbY67BMATN2XnFGB4Hag0PanDgdAhWlZrTMHBalbEYmytnpSYtEXYQl2iYA0JPWIvS0iZjI8KllF+/TqXEFbwqeCPmMwbi9SBDol7AopcsppRfyvoGU0u2U0m/7b9pRDCPFcliWw7HzNEIUFs9heW1RBRGEbDkUTpyiQ5TQk0zNeMjBFFx0YZEo9SosmYrNYdmRA5pGWCLDp4B9jwMu5c5tE6lIfYR7e5Hrr79+/G233TZqzpw5MxsaGmZ95zvfGQuwXoJTpkyZ/dWvfvWY2bNnz1q2bFnl5MmTZ1922WUTp0+fPvvCCy+c/Kc//alq3rx5x06cOHHOW2+9lXOLlO9+97tjL7744sknn3xyw8SJE+fce++9IwCgs7NTWrRoUcOsWbNmNjQ0zHrqqaeGAUBXV5d0xhlnTJsxY8as6dOnz37sscdqAeCb3/zmuKlTp85uaGiYdd11140f6N/db2jP2qzxZgCT3MdTSs8a6MmGDKwcREYpuYiH2nJ+CR7btx5oWRBBqGkKS+jYAK6QoEMOBkyx+SIjBRXwhAQNaloKSyAxWErGCbvxDi4Quj6MP2+AS7nT4Ajrtndvm9B0sKmo+2FNq53W92+f/recXeDvvffeXRdccEHZ5s2bNz333HPVf/zjH2vXrVvXSCnFZz/72Wmvvvpq5ZQpU1LNzc2xxx57rPmpp57auWXLlkhLS0vsf/7nf7bPnz9/x9y5c2fy7UV++9vfDrvrrrvGnHnmmdtynbOxsbFs9erVjd3d3fIJJ5ww65JLLukcN26c9vLLLzfV1dWZe/bsURYuXHjslVdeeei5556rHj16tPbXv/61CWCKcN++ffIrr7xSu3379g2SJKG9vX3AG40XEhL8I4C/A/hXALe4/oXIBUODDMfhCCUJviaMKyzhIUHmQNJDguIIS7MmFw45aNRkjYpFzc4NVo1nuPNqlCss0Tks4hCWHRIcBMrP+tF9XykUICLtCgju7UVmz549a9u2bbHNmzfHAO/WH4CzvYgsyxnbi+zatSua7zyLFy8+VFlZSceMGaMvWrSo6+23364wTZPcdNNN4xsaGmadeeaZDfv374/s2rVLmTdvXvztt9+u/qd/+qdxy5Ytqxw+fLhRV1dnRKNR8/LLL5/45JNPDqusrBxwpU4hxRM6pfShgX7xkIbVv8/ICAkKcDY8gZ9hi8Y6fx9Gi/8jATGtkv/0HJbwogs2PpRS6NQQrrAUCmiuBbk6Ndh+VCJDXDy3xtUxv5cA4SHBjKILU2d26QnfT59PCQWBoLYXybY9yCOPPFJ34MABZf369Y3RaJSOGzfuuHg8Ls2dOzf50UcfbXr22Wdr/uVf/mXc//7v/3bdc889e9asWdP4wgsvVP/+97+vfeihh0auXLly60D+1kIU1ouEkG8SQsYQQur4v0K+nBByHiFkCyGkiRCSbZsSftwCQohBCPlywZYPZliFDtwJOk5ZTNGFAtjkaYdKADFORk95nJ7wXAMvIbfIwbTGSRa55om3i3IRljEoFFaSKaz0aye4ZRQfL7dNjLCI2PHyEUFvLwIAr7766rC+vj6yd+9eeeXKlVWnnnpqb2dnpzxixAgtGo3SF198saq1tTUCAM3NzWpVVZX5zW9+s+Omm27at2bNmvLOzk6po6NDvuyyyzoffvjhlsbGxgGHUQv5o/6P9b87DEgBTMn3S4QQGcAvAJwD1rn3Q0LIC5TSTVmO+w8ArxVq9KAHL3Sg3jCcmBxWypPAt8mTO7+gO39bYSUeJrXHRihhOdeKE7rQRrP8mrmWHmi26hNZjZeCCimzsAEQrLBSdtEFV32aqTHCCkBhiUDQ24sAwAknnNB79tlnT29tbY3cfPPNeyZNmqRde+21HYsXL542Z86cmbNnz+6bPHlyAgBWr15d9sMf/nC8JElQFIX+8pe/3HHo0CH5ggsumJZMJgkA/PSnPx2wMu2XsCilkwf+pwEATgLQRCndDgCEkN8DuAjAprTjvgXgWQALDvM8gw9GCjIIUumxfhFOWU9kzWEJC3kZKciEZKkSFJgvgkMO3OGxHJYgctCtvKO7EIQag6AaLwnVrbDoIMphpRUVGdSAQohY5eczgtxeBACmT5+e+N3vfrfD/d6YMWP0NWvWbE4/dsaMGalLLrkk3ddj/fr1jfnO0R/6DQkSQsoJIf9KCHnUej2dEFJIQ9xxsPZMsbDLes/93eMAfBHAw/3YcB0hZBUhZFVbW1sBpxYMS0Vw52dXU4lwyrxs2/QShCxqsaeehEJIZpWgKIWlJy2F5V3QLLY1E1vvZMK0Q5S6aVi9BEXmsFJQiJSZfxQdqrRya0B6Dqt0FdZQRSEhwSfANnA8xXq9C6xy8KV+fi9bAi9dov43gFsppUZ6Qs/zS5Q+CuBRADjxxBMHf/tlQ2NOmaYlp4UorCRUILMARKDCcudBnApKUWrGaqOVHr6lAtWMRaIAc7wROeIoLKG5omRWwmLFDSIJK5FZ1m7qUIhU0grLD/i1vUixUAhhTaWUXkYIuQIAKKVxko9dHOwCMMH1ejwyu7yfCOD31teNAHA+IUSnlP6pgO8fvDBSUCBlqghTC37LCj0BmUgZa2eEhZeMFNQ0p0cAyMLCbwnWuooarELQVqAQ5+y0uKeIQJVURliCu7XDYAorkX5fi1ZYWjwPYYUKayAY7NuLFEJYKUJIGSx1RAiZCqCQJ/lDANMJIZMB7AZwOYAr3Qe482OEkCUAXjrqyQpgOQiSWU1lh06kWHC2GCnIRM5iCwQqLMVTGi1DYK7BaugKsLGxO4GIdMJ6kl0fsAmGM8mQhBddKEquogvRrZmsHz1VgrKf95VpmiaRJGnwR3yOIpimSQDkXJ9VyFT/DgDLAEwghDwNtsXIrf39EqVUB3ADWPVfI4A/WL0Iv0EI+UYB5z16oSegEDmzsAAI/sFOs0V4UQHPg7hyDSqRAGoAIhry6gnIEnN3OtW9CksYYSVY/gVsfGybJEl80QWRbcJy8n2Ci0H0uLMOy0y7r/y7xze0tbXVWA42RBFgmiZpa2urAbAh1zGFVAm+TghZDeBksLzUjZTS9n5+jf/uKwBeSXsva4EFpfRrhXznUQE9kUPVCOieoKegEMUmKk5cwsq2jSRkUu0JlyrE6tCiJ4FIUTvc9A89AYUTlqmnVb4JJCy3TXzCQ3xVDAXYlYTiuq/trikAoMXF2uXK+QGARjUokn/jpev6tXv37v3V3r1756CwiX+I/mEC2KDr+rW5Diikl+CblNKz4Wwt4n4vRDboCaiSnEESohSWLGUhT6FVgo49mqlB5s+7kQQggLBIJjkImVy4bFJdhGUXysgRQBdIDFocilSbY+GwwFxRWs6P/68QxTe75s+fvx/Ahb58eYicyElYhJAYmPcYQQiphVP1Vw1gbAC2Hb3Qk5CJ4gl7ySBsAIOeIetJyJICnSattkPuHJaA2boWz5ilK7wIRVAVJVczBjVcakYSqrBkiS3oNkzDUX2SKljJMOWXXv2qUgBan1i7eLd26uTXFElmoWZDB+RwC79SQL6reD2Am8DIaTUcwuoC62ARIhf0BBTVKSzQqQ6ZcBURcBjOSNohN5OarrJ2AQrCNABT8zg93dQh2yFBAbN0PcGUC1JMzXA1LKniVIOWgCIrAExoVINqqgDASEwkYWlxKJIK3WQ2ePJ9mkCFpSfYBANehVXmvq/kSlHWhSgichIWpfR+APcTQr5FKX0wQJuOfmgJyFHFfrAN05WnCTwkmGTOF2kKAgg+vGQ5W1VSkHJVCTpjI1JhpaCZmjM+ckQcYekJKFIEQMKbw5JVsUpGi7PJhuHkiZhdMcF2JaDIrNG4JyRo3ffQk0A0JKxSQCFFFw8SQk5B5n5Yv/HRrqMbegKKNMypzKOGS0UEXXSRgBxRAKtDumftTNCzdYsAZMkpAmG5BlfRRdDQE5DlCIA+Fn6zySEGpAQ5YT3JyMkiLPuaKVFxSsY0WZWgHIGmWe2PuFpXBCs/PRdhKfbnIUoDhRRdLAUwFcAawN6TkAIICSsXLFWja+4EsECFFY0xwkov2w56Vmw5NUVSnLU81OVYhOTUElAUNhN3VwnKSkScarAVFtJINArEe4TZBACqHMks4FHKBiVhqVxhDYE9sYYKCslEnghgFqU0XCBXKPQ4ZKs7Ae+eoEiCVISehCxXAqbX+clCFBb72xVJzZHDEtMcWI2UASYjT6ciLyrOCesJKCojLPckQ5EFKhmLsJQshKUqokOCcahKFIAT0tWpDkWxKk7D9kwlg0LWD2wAMNpvQ0oKdkjHcTiyqDyNa7buCS8BAgiLKyw1e+hGxExYT0K2rpVHzShlYhWW4igGR/XFBBaCsLFQpIg9EXNClTHx67AU1j0mew4rDAmWCgpRWCMAbCKEfABXSyZKabgGIRf0BGTJcYIGNexuCoE/PEbKJk9edCETGUQSkMDXXLN0vdO2yXYsIvIzesIqGmBrwuwSciUGJAU5YV5EoDObnFxRVByJWtdG5kRKdadxsSqQ3AFAj1uE1e3pus/ve6FkGqKoKISw7vDbiJKCabIGr9aDbVADhmnYC0EDf3i0uFVU4MzWZSIDannwBMEVVlpYyXEsvcHaA1hq2DU+PGSqxoCeg8HbY9tk3T/uMK5IJWNfO5fyc6tRUQUqgJWHZArL3edQka2enSJtC1FUFFIl+LcgDPEN+zcDvW3A5NOCOZ9dCedN5PPXSAXslLU+5ujgdE1QJAUQMSu2FZbqrRK0CEOIY9H6bNXgKftXywWqmV5GmPASg8rLxykFCtowoZg2WYTlClXaIUG1AugV2OBb64VcORpE9zac5iQmZCIUwhfk63TRjcz9qwC2gJhSSqt9s6qYWPEgsP0t4LsZm1/6Azs57Z2J2k45yBmyobOQoPXgGpSHJ2VAFTBbd83S3dudRCU+NmIIS1XLmC2eEnKBaibVx1QL4OnWLlt2Qk+y6xckbMKy7iVPvk9w0UWqF4iUQzGV7IQV9CQxhG/It3C4KkhDfEOkPNgb1qpI4ol87gRthRXkg80T5Yp3tq4QxQoJilFYspW45zYpUetWC9qx2IRebtviKKwKMU7YNAE9DjVSkWmTRWLQ+oInLHsiZt1L7iUSarnYPFGqD4hUQEmlE5Y1XiFhlQxKv8FWpCJgwkqbiVphJkVSACUWrC2WE7FDgpaTYSFBRZjCkhUnh2WYnMyJAAJl10JWLcLyOGFBhQR8kqE6hGUXXXCFJYIc0hQWJ1KJSJAigsjdtq0XUCtYFw53WTsfr5SgtWshio7Sb4uvVgCmFtwan1SmqnHyRgGrGutcsjXTdCoWedFF0AUgfPFpzJPDUmU1+IkF4FwrNZvCKgdMPfjej/yauQjLCQnydUUCqimta8PHildUMrUueOFwqg+IlEO1lkuY1IRJTbY+DCRUWCWE0icsvr9SUIlXLc0JWuW/siRbTllESNAiLFdZO9QyAb0EOYHGnI7fvG2VGnDo1mWPEmF95jIKCVzHBAZODDwk6FrMrHLCEqFmLLtk63niOSx7ImakWIg1aOgpNiGNVECxdkiwFamkApHKkLBKCEOAsCzHExRRWOEHrmrsHBYniSArlngYJ+IoCFZir7LwZNCz4lQvQGTIShS6qYNSytbLSAqbWIgiB4sIMqoEAQEqNJNEnfVOgmwC7Ps6PVTJ1DovBhERqrSeJ1dI0N4JWVKsSWIYEiwVlD5hBT1T5mEmiygN0/CGBINUWHxW7CJP7zosAQQRqbQ7W5jUtBYOK+w6BV3WzskhmkVhWYQR/Bhxm6oybRKtsIhkqzye71Ml1SEsEQu/+T0TKbd7VDr7hyliQs0hfMPQKLoAgptl2WE4V16Ex/qDTk6nKSw7JCjJgCwg75DqYaGbtA0TZWKFS4NeL2MrrCxqxiasoBWW1ybPwmF78iUohxWpsruScCK1c1iA0CIVPhHy5CFtwgoXDpcKSl9h8RxWYCFBbw7CUyUYdJ7Gzqc5zs9RNAKKLqz1MryvoqcgJRKw+gScnJpbYZkGCIg4wkqlqT5381vbJgGKIdUNRCrYZAcOudsFPIAYJcMnomp5BmGxiVBlGBIsIQwBwrIe8qAeJrvKy5U34pV5QedpstliumbFvGtCUEj1ehSWt8xeYJVghK2B587Ok5cRVQgSdWyyQ4Ixa61+UoADTr92puvaRbld3QLs8oYEPZ1BJEHVpyF8Q+kTlhpwlaBdTeVUeTmVeUFXCVrrntQ0W7hDpmawZdtWDosrLMM0XOHS8uCVAw+/Wfkiroa9eRkxhSCSq+jCVlixWnZMsitYm7hdViUeYFUvcrVuE6kIuyzyjlTlCQmGhFUqKH3CCrxKsBeQI54WNp6wV5BO2Zrx8pm5TRA8JAgEGy5Jy2F5FY2AkGAqsyLPrugMWpnbNrHrQaJZHHC0BgABEoIIy7IJcHJrMpEdhZXoDN4uPhaxarus3UtYYVl7KWEIEVaAIUErng6khU6CzhsluwFJ8awJM0zLIfN2SEGGcbS+rEUX9kw4aDVjKQKljCkXnpdhqqGGHRO0E+YOOFoF1doE1OnWrrLrJkrJZJlsqJIqVmElresTrbYXDodl7aWLoUNYgYUEmVO2Cwuoq5Q8UsG6FFgLG31HshuIVtl7cdmzYkkW42SssBIfG83QQEG9oZsgc2qJLiBSBUlWQUCyh7mCJqwkswmSDJnI3hwWsfJFInJFySyERdNyWCKUn1thcUXKy9p5qDlUWCWD0iesoNuzWNVU7tCJpzIPCE5JWITlLnIwqAGVqGKcTKoHUJ1Ks4RhNVTlzXhBg1egFjG5w28ykdl98//bO/Mwuaoy/3/eruotezr7AmQFkpCwhSAghiCgILIJCCIMKgPMDDq4IOK4wKC/RxFFERgEBlDBgYAIiKwisoUlhEAWwpIFspN0d5JO0kl6O78/3nur7q26tXR33VtV3ef7PPVU1b2nbn3PXd7ved/znnNiVcXxsDycWjtak4tuiui+YoTe9jRBzcBEH1arSfKiqi9IrEgeVhM483QG9mFV99eZMOwijj0CPV+wRKIdi7G7CaoH+NJ/kw92xOm/e7Yrl5QkB/WwnJBXVEamo8MRiIEJAd3TrjPbJwwLRBu+2bMt8b/xiri/v1HE8WYiNsJ7mhKcYhLTkKDryUBxOBkDu7b6rp1vaiYR5VwsD6t6AIj4BB6c+8oJ97Jra/TcLAqOni9YEO2gVKclWinJZekT0yFF7dU4xi8tyUFi0ffRtOzQrERPK323u9ilxKBmkJaL0rA4gg74OuwT4lAzsDghweoMXh84HlbEwtC6S72UmoG+hlhiaiaXV7E8rAAv2f2eEKzdVrB6AnqHYEU5YNcJ6fg8rIRX4xjlqB6elJCgK55F6XdwDX9N0uPzeVi1rmBFuCz97gBj56bZO1wjF4cgTh1F9rAS126gNrwggNfA4nhY3kZH0DismiLcVxahoXcIVk2ED5PTQk4kFnj6IJJGOVrB8s0skTZ2JiIPwmP0XCOX6MOqiEcv5uAPv1XE/Gn2DtfoPaztmcOUUCQRTb92rri7kQT1sIqQDLK7KREtCExrT4QErWD1BIQqWCLyWRF5T0SWi8j3AvafJyKLnNc8ETkwFCK1g6K7YXdvS2QsAf4U2yJ5WGnT6UgM4tWaWBCVQfYYPZfPnjavh1WEvoaA1nlC0KGIfVgpYUqv11dkD8uX/eoV9+oB0TV+vNi1JSlY2UKCtg+rRyA0wRKRGHAzcCIwFThXRKamFFsFzDbGzACuBW4LhUzt4GgEq61F09arMyQWRO5hOX1YEhASBDUykYcEk31Y7rnxeZ9RelhO4wJUHBLjwqRIfVjG6H3qnAtfH5a3r6i9JdoJcBPXbnDw+EKXVzFCgrsaoU8dkDxfrcZpJEqRQs0WoSFMD2sWsNwYs9IY0wLcB5zqLWCMmWeMce+kV4GxoTCJSrDclq/Hw3K9CF9iQRRGuaVZxbO2DhHxjelJduAPjK61njMk6CSBRCXmnvPjckg3whELVstOFaM+QwB/mNLXyIBovSz3fvXc14HeaNThU2OguTFxvoLT2gdoyr0VrB6BMAVrDLDG832tsy0TvgY8EbRDRC4WkTdE5I3Nmzd3nokrWGEPSnUf2OoMiQXxKk0AicIo72rUd6f1GZNYQBZcMTysQQlvoaW9BXDOTcxZHTYqDytxftKNna8Pq7U5uvkWXU4pIuoThmL0yTQ3JHgFZpwC9B2q1y7KVYdbdmj2onu+JECwRNTLslmCPQJhCpYEbAtUDBGZgwrWlUH7jTG3GWNmGmNmDhs2rPNMageDaQ+/U9gVotqkUfaFvUC9rCgeHtegOQ9zrCLmH8QM0baKXT7VA9LS2pMhuEHReVjNGQTdpIQEIbpz5ApDn6QB9q0ZBioMADvro+Hk8pIKqB3sCy8npmby8nLrEBUvSAsJ+gQLoouwWISOMAVrLbCX5/tYYH1qIRGZAdwBnGqMCeduj6pVmniAhiaNsjfsBU4CSARG2TXITt3jFQHGL8qQV3O9conFg/v3INqWsMdrcDmkhd8c7ysyI9ycwcNyx/EB9HUabDs3RcMJVBxr66CiIi28nDhXxeDVHOAlp2YJgtMQsoLVExCmYM0HJovIeBGpAs4BHvUWEJG9gYeA840x74fGJHLBSvYb+TLhwPGwIhCJlJCgN1ySCHn1HQY7uxBi7Qp21kMfbYW7gumKeYJPlC3h1PPjpJD7BMs1wjsiMsJu3VPClK2m1eNhucIQpYdVn/SgvLw6gnhFdD9BYAgVPA0h11PuUxet52cRGkITLGNMG3AZ8BSwDJhrjFkqIpeKyKVOsR8BQ4BbROQtEXkjFDKRC5bTaS6xpFF2H+yovIhESFDrHquI0drhmWwWoN9w5eKIaqhobkgYvcC0dlCPLyrBSmmduyFBXxZlv+H6vuPjiDl5PKxUTrV1gEQrDDsbEo0NLy+/uDvnKkohdRsS/YYleIEn1Jy4jiNge0TX0CJUxMM8uDHmceDxlG23ej5fBFwUJgcgutHuzQ06EadnXIgvsQBUQNa/FS4PUCMDvtZnWp+Rt1U8MJwEzSSfehgyMcEFSBfzfsPho3nh8nDR3ABIQtArKypp6Wjx92EljHBE4rBjo2a01Xr61VK94lhcRTZKwWquh+FTEl+9vJKC5fatRchr+wZ97z8KIBE2TQvD9x+locr2Nj1/FmWL3jHTRcLDagz3f5rr1ZiI5pvEKmLpYa9+I/Th6egIl8uOjSrUlbqQZExi6X1G/UY4ZSMIeXnCSqnjsBL9M/1H6TWKwuNrWq+CHdP/DgwJ1g5WAYkqJLh9I/QfCRUVCU5u6C1xjiDaUC6oh+mKdyZeNQOhojK6cwV6vmoGJVaHzuhh9R+p81hGec4sQkHvECy38zzscMXOhuR/oYY4EfYST2uvoy38mPr2jYmWJzgeVmrLs19EHkRHuzNeJjgkmPCw+o/U9yhCcNs3JP+PZMg0sYAjqHD0Gx5dIkHTev818w5mrvB4Bn2HRhd627ND+1wHJkekxCvi6Y0fEUdIIwwJBpwvgF1tupRI8r5yyrgemUXZoncIVrxKH6amtCTFwmK7/wEK9GpcIxn2w+O21h3EJZ7e8owqqWDHxzqsYMBo3//valfDUul4OfRzz83GcPmAnn+HD2QTh2GwI6KWeeo1C8pcTHCKyutz7tMBScGqrKhMv5dAhTSq/j4IPF+gnnu8Iq7rhwH0H5Esb1HW6B2CBSokYYtE03qfEfSFBFO9iLAfniAPK60z2vWwQjZ+29bpu2P0EoOq07zPCAWrKd3DChSHfsOjM8LbN6RdM3fyZF9IsP9ILRvF6szb1uq79772NMR8vAbuBU3rwufkIuAeBw0J+s+X9bB6CnqPYA0YrUYqLLS1aKvX0xKNS0DSRSLsFaJR7ujQ46e0PtOSLiprdfBw2B6Ea8ScsFLGcVgJwxKyYLXt0T61FGOXNnUVqNcXhYC27tKMzRQRdRfdTFwzgEF76wwczSH3yUIyKuH1RoMaPy6vraujEdIM9zio5+47X32HA2I9rB6A3iNY/UdpyC4s7NgImNwPdr8IwhPNDdpPlqsPC/SBb1obHhdICtYAv2C5fQ2+gboV8XDFHJIek+f8VFZUps/eAGqEd2wMf7LZROjNH6ZsM220trcmw6YuJ4CtH4bLCZKC1d9/XyeuXaqQtuyIZmjCzk16j6c8b6Aelu8ej8X1Wm9bk3oUizJD7xGsAaPVkIeVgZZoifo7p12RSBjBeLWmLYcZnkik+/pbn66R8RnkweNgy0fhcQENCcZrk2PCUgYOJ/hUVKhHsy3ksFKTPx3a5eSG33zGbvA4fd+6OlxObugtpZGRlrkIMGifaDiBNmb6DE1km4KeK7f/Mc3DiopX40p9rxuf2JRRsECHVDSsCJ+XRajoPYIVdhw74UX4Y/2BoZOBY5IGKgy4BiMlsyvNowFHsD4MN4zTtFa5uOn+qX1YXj5145PGKCxsWaXvrhiRYaJZb5ktH4bLyTWmQyalcUoPCToznkUhDFs+Sv6fA2/2a7DnFwGvhuX6XjenVnfgAAAZ0ElEQVTRxwugua3Z3yiDaO4ri9DRewTLNd5bQwoLbAsQrIoY7aYdSPFq6iaE29prdI6d8jAHiufgcbpURZhhnK1rfAOTRUSzFoNClHUTkvzDQsNynczVI1gZGxeDXW8mZC+0YbkuqJniobvjnXycagbq+KMohKH+Axi6r29TYH8oeAQr5HMF+vxUVGqih8vLk9buE1LQZ6G5PvolUCwKit4jWG7L1W2ZFRoNy7UPxl0wDr/h8xmcIZO0xR7WshUNyzWMk8LFOJPlR+pBGBNo9GIVSYHwifmQiRq6DXOC4IblGlaLVyU2BY4tAu1zjNdE42HVTUgMGnZ5pE2B5GLwuPBDXHt2qHc8dLJvc7wiHhwSrB2k9119eNOCJtC4Qs9BLP0Z29W2K8DDmuD8znpZ5YzeI1gDxmo/Sv0H4Ry//v30lqin9en3sCbquKSwWsgNK32hpdT/DxasVeFwaVoPLdvTBUtiwQLqeoVhelkNy9OMsJt0ASleg4gau7DuGxeNK5JG1YF3TbU0wRo+BTYtC5dTg1Pnoful8Qr0Rl1eH78TLi9w7vGJvk0JwWoNECy3bH1IDVaLSNB7BKuiQo14Q4SCldHDch6esLy9xhXpD7PHCPsMct0EDY9tDqlVXP+evg/zG73c5yYkwTImUNATU2elfAZgxDT4eEk4fECHRDSuSuPknpcO05FugIdP1ezFMFPbXZEOuK9bO1p9HH28Ni0Ld+qx1t16X3nmN/RyaeloSQ8JDpmkGaibIhBTi9DQewQLYOikcFrKOxs0jDUsvSXqwh/2cgxTGKGTXVs1sSSD8Uv9TGUtDJkMGxcXngskhXBosGDFJEaFeG7DuokQqwqPT+NKaN2ZVUDTxGHEAZpUE5Y4bF6mK+eOmuHnJBmuGaiIAny8NBxOABsX6bVI8fwy3ksAI6bq+d0WYv/apqWa0j7qoIy80q5hvBqmnurLnLUoP/QywdpPO4Rbdhb2uJud0ExAP42LtClsBowJZ9b2Dc4xR/sfZm+LM+1hHjk9PIH4eLGm8fcb7tvsinmawYtXqUBsCGlG+/UL9X30IYF8AjmNPEDfw/KyEpwO9m3OLgyuYIXo+a17E0bO8PX1pXIJFHcId0WC9cH3eFZeAGfeCYdfEh4vi9DRuwRrzCE6a3OhH6a1zjJenTE4ow+G9W8WlgeokYH01me21vrI6doiDsODWDMfxh6WSGlP5ZDGBdQQrX87nFT79Qs1iSJDOCmQ00jH8wnLCK9fqJl/g8f7Nmc1wP1GaL/smtfC4dTeprzGzkzblbFvFhyBqwmPF+hzUzMoOR7N5ZVLsCzKHr1MsA7V97XzC3vctfM1bOJZlRVyPNijD9bwVKHTyde/qYbPWQQwwSWXeEJSeAuFXVu0r2Gvw9J2ZfSwXD57toXTx7dugQp0Sh+HLylFUjj1G66hyo9eLjwfUFEffUiaqGf00EHL7nOErh8WhrBvWqrTP7nPjAdZ76V4FYyZGe66ZqtehH2OSjtflZIlimDRI9C7BKvvUDXmhRQsY/R4YwOMcjaD45Zf83rhuHR0wEevBHLJamT2mqV9FR++UDguoIYYYOysjHwCDcs+R+n7yn8Wls+e7Xqtxn0ybVfWkCDobz6ap0ulFBLbN6o4TDgmbVfGRBkX+xyp00yFkaq94h/6Pu7odF7Z7iVQId24KJyhCY2rNKw/4ZisvNKSLix6BHqXYIEaww9f1JBHIbBpmRqNfY5M2+U+QIL4J1QF2OtwqOwDy/9eGB4AG9/WwZGTPp2RS+pnQBMvxs6Clc8XjgvAB0/rUIK9MgtWoMGrm6DhnhXPFZbPqhe0s35i9vOTliUIMP5TOsB6XYHDuG4dJx6blVOgAZ5wjL6//1RhOQF88Ix6ogNGpe3y3UtBQjr5BA29f/BM4Xm5z8vEOVl5WQ+rZ6L3Cda+n9HR7mteLczx3n9C3yd/Jm2XN+wlKeELKmu09frBM4UL6SQe5nTjlzXkBWoANi4q3EwgxsD7T+pxnRVhg/gEchFR0V31PLQ0F4YPKJ+qftpYSIFXpAKN3aRP68wKyx4pHB+Adx/T+RPdZAUPcnoydRP0d8v+WlhOOzbD6ldVeAKQ0xsdM1NnSH+3wLwAljwEw/ZPy4JN5WIFq2ei9wnWxDka/lr2WGGOt+wxTXAIaIm6D03Gh2e/E3XAbiEy4oyBxQ+qp5SSkQd5GL9pp+v70r90nwtof9i2NbDfSYG7Kz1L0wdi2hk68/d7jxeGT+suWPowTPl8WtYb5BF+qx2sDYGlDxcuLNjcqN7R9DN9M1wkeOTyZACmnAKrXynsBMaL5+rA9ulnB+7OeS9VVGgK+XtPFjaRZ8tHsHqenq/UBiBWsHoDep9gVfdXoVh0f/eXjNi4WJMcDjwncHfWxAJQkYjXwMJ7uscDlMfmd+GgLwXuzmlkhkzUDvaF9xRm0OeCu6GyL0w7LXB3VYWKRsa+hn2O0nniFtzdfS4A7zyiIb0M1yrn+QH97bY1GuosBBbeo+OvMnHKltnp4uDz1Hi/cWdhOHW0wxt36b0wfP/AIhlnTfFi5legfQ+8+YfC8AJ47Xc6+PfAcwN32z6sno/eJ1gAh14Iuxph8QPdO86rt6rgZBIsJ8yUsbVXOwimngZv/Ql21nePy8u/0cUYDzgjcLeXQ0Y+sy7RrL7uGuRt6/TczjhbGwhZ+GT0HCoqYNbF2t+4upsp0h0d8NKvYdgUGPepwCJ5CdaUUzSV/KUbuh/Gbd0Fr9ysfWMjp3ed08Cx6jXO/9/CLMS59C86G8yRX89YJGf4FHSc2IRjYN6NsLup+7y2b9TGy7TTfRMpe2E9rJ6P3ilYE+ZoGO+F69RwdAWb3oW3/wQzv5pY5ykVWRMLXBz9LeXwwvVd4wGaqv3OI2rgawZm5ZKVzwFnaLLDMz/S6YK6iud+qp3uR38rYxG3BZzVsBz2NZ1M9e9Xd8/rW3CnDu6efUVg6C2VR8bzE4vD7O/qGKMlf+46H1AB3bERZl+ZsUh1rDqQXxrm/EBT0J/7afc47W7Saz98Gkw5NWOxvO4lgOOu1hlgXriue7wAnv4htLfAMVdl5pVtGIlFj0DvFCwROP6/dfLZZ/+7879vb4WHL1VxOPrbGYtlTd12MWw/Fb3XbtWU9M6iZSc8cpmu93XUNzIWy88gV8KJ16mX1VUj896T8Na9cMR/JJebyMInq8Gr6gvH/Vj7LV69pWt8GlbAM1dra39asPcJUBWrCvychoO/rI2dJ77b9cmLP3oFXrweDjgzMMU+iEfW8zRsXzj8UlhwF7z7t65xMgYe/45OVnzKjRmFHfL01kHH0x1yAcy7CZY/2zVeAAvv1X61o7+VNkemFyLJbFwrWD0TvVOwACbM1hDYq7fAm3/M/3cdHfDoN3QWgM//Jm2wsBduCzmrsQE4/hpdd+n+L3dursO2PfDni7Tv6pSbMnpXqRx8c/elYr/PwkHnwQu/UEPRGax9Ax78qs52kKUlDHkKFsDB58P+J8MzP+x8Qsi2tXDPFzTJ4vO/Ceyod5G3YFXE4Av/q42We8/WUFVnsOFtuP889WRP/lXWonkLFqiwjz5Y74dVnRxP19EBT/9A+3XnfD9wdgsvvJ5f1nMF8Nmf66wi958PH77UOV6gSS5//YZm1GbxRl24KwDYPqyeid4rWAAnXKtjch79OvzjJyoA2bD9Y7jvXA0FzvkvzYTKAu+Dnb1gfzjPCTHdcZymKefqI2lcBX84TbPoTrwOJh+XtXhOg+fFyTfA+NnwyL+rIcs192JHO8y/A+7+nAr4l+b6llQPQs4MShcicMZtOhj6ga/As9fmTnU3Bt59HG4/VkNS597nW6wxCN5rlfO6DZ0E59yrHtZtc+D9PPr82ts0meGuk3Rs2nkPZG1gQDIxJS9O8Wo49371av94Ojz/i/ySirZ8CH86C165CQ77V/jUFTl/4uWS876q6gPn/0UXNv39KfB8nmH4XVvgyavggX/RWUDO+ZM2FnLAOM+N9bB6JjphxXog4tVqeP72HfUo3vo/mHkhTDhW10uqrNUH5+OlGmp5614deHriL+Dwi3Me3jU47jpLWTF0EnztaZh7gXpaY2dpxt/en9DWeKxSByhveBveeVT7UOI12tqffmbuqmZKbggsXK0G9YnvwrzfwqK5GgqbeKwuH1E9QGfkblihSREL79GZ58fP1glGs3idLnKmtXtR1RfOf1hDVi9eDwv/qF6gy6dmoPbhbF2tKd5v3wfr3tAkiwvuSps3MPAvYp0QB9Bkia8+AQ9drAZ/zEyYfhaMO0rFMV6r64DVL4eVzymnxhXqKZz+u+QK2IXk1H8EfPVJ+Nu34bmfaCNixtkw6Tg9B7WDtXHRtE4HQL/3OCx7FGLVcNL1cNhFWb3QTnHx8RoJF/0dHvum9rO9frves5OO0+SMWmcasZ2b9P5e/qwO0djTpJxO+GnOBpCLwDXWLHoM7FWtrIXTboYZZ2nr7x8/0Vcq4jXqUc2+Mmsc3Qv3we4weSYMDJkIF/9Ts6FeuxUeuzy4XFU/7ff65OXacs0DnQ6RxKs1jDbjHM2Ke/FX8OIvg8uOmalCNe2MvAweJFvAeRu/qj5w2i0aInzxl/Dyr+GlDCG1ofuqAT70wrQ5AzMeviLPkKAXow6ES16ABb+H+bfDk5lCVqKDlY+/Bvb7XNb+IR8nD4+aeH4Gm9rBei0O+RfN0Hv1Fn3PVHbWJdrfmIeAuui0YIFmxJ51lwrQvN+qmGbql4zXwpST4aj/zJhBmQvWw+qZsILlYsIx+mpar1l3DSt0jExVf/W29jocqvt16pCuwclbsEAN7Kx/1Qd787va4ty+QftM+g5TLmMPU0HpBGpieRq8VOxzhL6aG3Xewy0fwu6tKpoDx6gn2Alj56JPvA8AtfH0WTBy83lQ+ax7Uz27Pdu1Bd5/tPa/DB6Xt3C66FRI0It4tXrbh1+s52btG9p31t6qjaG6CTqmqf+ITvGBLnhYXkyYra/d23T+xIYV6rEgmqAzYqr2NeYRZktFlwTLxbij9LW7SfuB69/X+8mgnvnwKdoQCJgdJR/EJEa7abeC1UNhBSsVA0bn7bXkgvtgu2GKTkFEH948wln5oE9ln24eoE4TMgoEV6i6bPz61Gm/XY6+u3zhS7qoyNPDSsXgcTn7yjoDL48uNzhqBmrobVJhzhNAdScbS4GoGZAU1QLCfdb6VvYt6HEtSgO9O+kiZLjGuBRae65HUyrotoAWGF7hTJv3sUjwep8FEYkCoVseVshwky76VXYuGmJRHghVsETksyLynogsF5HvBewXEbnR2b9IRA4JOk65wjUypdAB3OnQW8jolvcZAkrRwHlnlCglkShl78W9n0qtQWRRGIQmWCISA24GTgSmAueKyNSUYicCk53XxcD/hMWnGBhWOwyAXW1dnE2jgCi1B9gdC1YqglXq43ayjp2LGENrc2eBFhul2ACx6D7CbPrPApYbY1YCiMh9wKnAO54ypwJ/MOrHvyoig0RklDFmQ6FIXPPXpbyzvgBzmXUB7eyEGti6vZIv/q4Ls1gUEB20Qg0MaTuh6FwAdkk7UlXN/MXj+eKi4vMBwOkmKoXzk0AJcjK0KS8jJcULoC5+PI3xZ7hy7iripM+tOHX0AH78+WlFYGZRCIQpWGMA7+JKa4HUhYiCyowBfIIlIhejHhh77515up9SQ4y+jG69kL7t+xWbChVUsv/um6mgNDyJWjOO/ffciFAa/UUAY1ouIkZpeaJ7t1yOmOKHlL0Q4oxoPYs+HfsWm0oaRradxbC2z5XcdbQoDMJ8EoIsUWr8J58yGGNuA24DmDlzZqdiSMVvTR1R5P+3yB+leK1KkROULi+LnowwA+Nrgb0838cC67tQxsLCwsLCIlTBmg9MFpHxIlIFnAM8mlLmUeACJ1vwE8C2QvZfWVhYWFj0HIQWEjTGtInIZcBTQAy40xizVEQudfbfCjwOnAQsB5qBr4TFx8LCwsKivBFqb64x5nFUlLzbbvV8NsB/hMnBwsLCwqJnoHQGd1hYWFhYWGSBFSwLCwsLi7KAFSwLCwsLi7KAFSwLCwsLi7KAmFxLsZcYRGQz8FEnfzYUqA+BTrFg61O66El1AVsfgH2MMcPCIGPROZSdYHUFIvKGMWZmsXkUCrY+pYueVBew9bEoLdiQoIWFhYVFWcAKloWFhYVFWaC3CNZtxSZQYNj6lC56Ul3A1seihNAr+rAsLCwsLMofvcXDsrCwsLAoc1jBsrCwsLAoC5StYInInSKySUSWeLZdLSLrROQt53WSZ99VIrJcRN4Tkc94th8qIoudfTeKSORL4HamLiJyvIgscDgvEJFjS6kuna2PZ//eIrJDRL7j2VaW9RGRGSLyiogsdfjXlEp9OnmvVYrI7x3Oy0TkKs9vil6XTPVxtn/dedaXish1nu0lawcs8oAxpixfwKeAQ4Alnm1XA98JKDsVeBuoBsYDK4CYs+91dPlUAZ4ATizxuhwMjHY+HwCs8+wrel06Wx/P/j8DD3jLlGN90BUQFgEHOt+HlPG99iXgPudzH+BDYFyp1CVLfeYAfweqne/DnfeStgP2lftVth6WMeYFoDHP4qeiD94eY8wqdP2tWSIyChhgjHnF6F37B+C0cBhnRmfqYoxZaIxxV2VeCtSISHWp1MXh2Jlrg4icBqxE6+NuK9f6nAAsMsa87fy2wRjTXir16WRdDNBXROJALdACNJVKXSBjff4N+JkxZo9TZpOzvaTtgEVulK1gZcFlIrLICRUMdraNAdZ4yqx1to1xPqduLxUE1cWLLwALnQez1OsCAfURkb7AlcA1KWXLsj7AvoARkadE5E0R+a6zvdTrE1SXB4GdwAZgNXC9MaaR0q/LvsDRIvKaiDwvIoc528vVDlg46GmC9T/AROAg9CH7pbM9KB5tsmwvBWSqCwAiMg34OXCJuyngGKVSF8hcn2uAG4wxO1LKl2t94sAngfOc99NF5NOUdn0y1WUW0A6MRkNo3xaRCZR2XUCvwWDgE8AVwFynT6oc7YCFB6GuOBw1jDEfu59F5HbgMefrWmAvT9GxwHpn+9iA7UVHlrogImOBvwAXGGNWOJtLti6QtT6HA2c6HeODgA4R2Y32aZVjfdYCzxtj6p19j6N9LPdQovXJUpcvAU8aY1qBTSLyMjATeJESrYuDtcBDTnjvdRHpQCe9LTs7YOFHj/KwnFi0i9MBN3PoUeAcp69nPDAZeN0YswHYLiKfcFpgFwCPREo6AzLVRUQGAX8DrjLGvOwWKOW6QOb6GGOONsaMM8aMA34N/D9jzE3lWh/gKWCGiPRx+n5mA++Ucn2y1GU1cKwo+qIey7ulXBcHDwPHAojIvkAVOkN72dkBixQUO+ujqy/g/9DwRSvaQvoa8EdgMZql9SgwylP+v9CsoPfwZAChLcYlzr6bcGb/KNW6AD9A+xXe8ryGl0pdunJtPL+7Gn+WYFnWB/gymkCyBLiulOrTyXutH5q5uRR4B7iilOqSpT5VqEe7BHgTONZTvmTtgH3lftmpmSwsLCwsygI9KiRoYWFhYdFzYQXLwsLCwqIsYAXLwsLCwqIsYAXLwsLCwqIsYAXLwsLCwqIsYAXLossQkRtE5HLP96dE5A7P91+KyLcK+H93i8iZhTqe57jf93welzrzdxYuq0Tk0gJx+KaIrBaRmwpxPAuLnggrWBbdwTzgSAARqUBnE5jm2X8k8HLA70oN389dJBBXGGNuLQQBY8wNwI8KcSwLi54KK1gW3cHLOIKFCtUSdMaAwSJSDUwBForIj0RkvogsEZHbnJkTpojI6+6BHM9mkfP5UGfS0gWO1zYq9Y8zlRGRf4rIz0XkdRF5X0SOdrb3EZG5zgSv9zsTo84UkZ8BtaLrQN3rHD4mIreLrqX0tIjU5joRjsd1o4jME5GVricoIsc4POc6fH4mIuc5/BaLyMSunnwLi94GK1gWXYbRZU7aRGRvVLheAV5D1xWaiS6z0QLcZIw5zBhzALpMxcnGmGVAlTOZKsAX0UlKK4HfAmcaYw4F7gR+6v3fPMrEjTGzgMuBHzvb/h3YYoyZAVwLHOrU4XvALmPMQcaY85yyk4GbjTHTgK3orPj5YBQ64e3JwM882w8E/hOYDpwP7OvwuwP4ep7HtrDo9ehRk99aFAWul3Uk8Ct0WYYjgW1oyBBgjugyG32AOnSqn78Cc4GzUeP+Ree1H7ow5TM6rRsxdOodL3KVech5XwCMcz5/EvgNgDFmievNZcAqY8xbAcfIhYeNMR3AOyIywrN9vtH56hCRFcDTzvbF6GKDFhYWecAKlkV34fZjTUdDgmuAbwNNwJ2iy8PfAsw0xqwRkauBGue39wMPiMhDgDHGfCAi04Glxpgjsvyn5Cizx3lvJ3mPd2bJ8z2ez+2oV9jZ30mG7R2e7x3YZ9DCIm/YkKBFd/EyGgJrNMa0G13gbxAaFnyFpDjVi0g/IJHlZ3RplHbgh6h4gU5KOkxEjgAN/4mu/eVFPmVS8RLqzSEiU1GBddHqhBktLCxKGFawLLqLxWh24Ksp27YZY+qNMVuB251tDwPzU35/Pzq7+VwAp8/rTODnIvI2Ohv9kd4f5FMmALegIrcIXeF4ERq2BLgNWORJurCwsChB2NnaLXoFRCQGVBpjdjuZec+iyQ8tXTze3cBjxpgHC8jxQjR0elmhjmlh0ZNg4+cWvQV9gOec0J8A/9ZVsXKwDbhWRIYWYiyWiHwTuBRdadnCwiIA1sOysLCwsCgL2D4sCwsLC4uygBUsCwsLC4uygBUsCwsLC4uygBUsCwsLC4uygBUsCwsLC4uywP8H4cGslYSffo0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "with pt.Environment(wl=wavelengths, freqdomain=True):\n", " detected = nw(source=1)\n", " nw.plot(detected)" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Add Drop Filter With Grating Couplers\n", "" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "## Simulation & Design Parameters" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "slideshow": { "slide_type": "-" } }, "outputs": [], "source": [ "dt = 1e-14 #[s]\n", "total_time = 1000*dt #[s]\n", "time = np.arange(0, total_time, dt)\n", "\n", "distance_x = 100.0e-6 #[m]\n", "distance_y = 30.0e-6 #[m]\n", "center_wavelength=1.55e-6 #[m]\n", "bandwidth=0.06e-6 #[m]\n", "peak_transmission=0.60**0.5\n", "reflection=0.05**0.5\n", "\n", "c = 299792458 #[m/s]\n", "wg_length = 89.22950569e-6 #[m]\n", "ring_length = 50e-6 #[m]\n", "transmission = 0.7 #[]\n", "wavelengths = 1e-6*np.linspace(1.50, 1.6, 1000) #[m]" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "skip" } }, "source": [ "## Circuit" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "with pt.Network() as nw:\n", " # components\n", " nw.src = pt.Source()\n", " nw.through = nw.add = nw.drop = pt.Detector()\n", " nw.dc1 = nw.dc2 = pt.DirectionalCoupler(1 - transmission)\n", " nw.wg1 = nw.wg2 = pt.Waveguide(0.5 * ring_length, loss=0, neff=2.86)\n", " nw.wg_in = nw.wg_through = nw.wg_add = nw.wg_drop = pt.Waveguide(\n", " length=wg_length, loss=0.0, neff=2.86, trainable=True\n", " )\n", " nw.gc_in = nw.gc_through = nw.gc_add = nw.gc_drop = pt.GratingCoupler(\n", " R=reflection,\n", " R_in=0,\n", " Tmax=peak_transmission,\n", " bandwidth=bandwidth,\n", " wl0=center_wavelength,\n", " )\n", " \n", " # links\n", " nw.link('src:0', '0:gc_in:1', '0:wg_in:1', '0:dc1:2', '0:wg2:1',\n", " '1:dc2:3', '1:wg_drop:0', '1:gc_drop:0', '0:drop')\n", " \n", " nw.link('through:0', '0:gc_through:1', '0:wg_through:1', '1:dc1:3', '0:wg1:1',\n", " '0:dc2:2', '1:wg_add:0', '1:gc_add:0', '0:add')" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Simulate Time Domain" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEGCAYAAAAjc0GqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nO3deZgU5bn38e89M6wyIMvIMoAgAQFZHUQ07kZFYyRqPJITozFxD8YkmlfNOUdN9I1bEqPGSFCJMckbjSZGDpq4xaMxxwUQUBGIuIVBRzaRfZievt8/ugeaYZbqmq6p6Z7f57rm6u6qeqruYul7nqWex9wdERGR1lYUdwAiItI+KQGJiEgslIBERCQWSkAiIhILJSAREYlFSdwBZKtPnz4+ZMiQuMMQkXZiwYIFa929LO44ClHeJaAhQ4Ywf/78uMMQkXbCzD6IO4ZCpSY4ERGJhRKQiIjEQglIRERioQQkIiKxUAISEZFYKAGJiEgslIBERCQWefcckMRn4eqF/GPVP0KVLbIipn1mGuXdynMclYjkKyUgCeyuhXfxStUrGJZ1Wcep9VounXhpBJGJSD5SApLAEp7goH4HMfuE2VmXPei3B1FTWxNBVCKSr9QHJIG5e6jaD0BxUTEJT+Q4IhHJZ0pAEpgTPgGVFJWQSCoBicguSkASmLsTMv9QbMVKQCKyGyUgCaylNaBar81xRCKSz5SAJLCWJKAORR1UAxKR3SgBSXCeep4njGIrpiapUXAisosSkATmhO8DKikqoTapJjgR2UUJSAJLelKj4EQkZ5SAJLCW9AEVm54DEpHdKQFJYO6OWfhBCGqCE5FMSkCSFTXBiUiuaC44CaxFTXBFxWyo3sBrH7+WddkiK2J079F0LO4Y6toSjLtzx8I7+HjLx6HKb9yxkaE9hoYqe/Sgozmw74Ghykr+UgKSwFrSBNe9Y3fmVc3jnL+eE6r8jAkzuHD8haHKSjCVmyq594176dW5F11KumRVdtXmVQD874f/S0lRdl8r1bXVPPn+kzz1paeyKlcnkUxkfU1pG/S3JoElCT8K7tpDruXM/c8MVXbGszPYXLM5VFkJ7uG3HwbgtqNuy7o2cu5fz2X+x/O5+3N3c3D/g7Mqe8PLN/DQ8od4ofIFhu89PKuyj7/3OHcvupsJ+0yga0nXZo//bPlnmT5yelbXkOhEmoDMbCpwO1AM3OvuNzVwzFHAz4AOwFp3PzLKmCS8ltSAenbuySEDDglVtqSohKQnQ5WV4Ko2VwEwtmxs1mWvPvhqnvngGQ7cJ/tmtG+M+QYPLX+Ibz77zazL1tm0YxObdmxq9riNOzaGvobkXmQJyMyKgbuA44BKYJ6ZzXH3tzKO2Rv4BTDV3f9lZvtEFY/kRtgaUEsUW3HeJaDaZC23LbiN9dvXhyp/+ojTqehbkXW5JeuW8MwHz/CVUV+hQ1GHrMqu3raaCWUTsi4HMKLnCEb0HJF1OYD+3fpzz/H38NHmj0KVn7jPRIb0GBKqrMQryhrQZGCFu78LYGYPAtOAtzKO+XfgT+7+LwB3Xx1hPNJCLakBtYSZ5V0C+nDzh/z6rV+H6k+pGwQQJgHNWjyLv638G/e+cW/WZQGO2/e4UOVaakr/KbFcV+IVZQIqB1ZmfK4E6jcOjwA6mNn/AKXA7e7+QP0TmdkFwAUAgwcPjiRYaZ7jsVy3yIryLgFV11YD8P2Dv88JQ07Iquypj53KtsS2UNetK3flQVeGKn9Y+WGhyomEEWUCauhX5frfYCVABXAs0AV4ycxedvd/7lbIfRYwC2DSpEnxfAsKjoeejLQl4kpAO2p38MHGD0KVXbp+KQCdijtlXbZTcSee+dcz3PjKjVmXXf7Jco4aeBRnjT4r67IirS3KBFQJDMr4PBD4sIFj1rr7FmCLmb0AjAf+ibQ5LVmSuyWKrIgkrZ+Abnz1Rh755yMtOkeYZ5fWbV8HwKMrHg3VHzO5/+Ssy4jEIcoENA8YbmZDgVXAdFJ9PpkeA35uZiVAR1JNdLdFGJO0QEseRG2JIopSq7G2snXb1jFgrwFcPunyrMte/nyqTMei7BNQXTPaj4/8MUcMPCLr8iL5IrIE5O4JM5sBPElqGPZsd19iZhel989096Vm9lfgdSBJaqj2m1HFJC3TkiW5W6KoqCiW1VQTyQR7d96b44ccn33h51MvYZrgqhOp/qOenXpmf12RPBLpc0Du/gTwRL1tM+t9vhW4Nco4JDfirAHF0QeUiyfswzTBTR06lT+v+DP9u/Vv0bVF2jpNRiqBtbdh2AlPUGItS0BhakDXHnItz/3bc/Tp0qdF1xZp65SAJLC4akBxPYiaSCZCDQIAGNJ9CACdSzpnXbakqETJR9oFzQUnWYlrGHYcgxASyUSg+cUa8tOjfsqiNYvo27VvjqMSKRxKQBJYHEkAUk1wcQxCqEnWhO4DGt5zOMN7Zjexpkh7oyY4Cawls2G3RLEVxzILg6b5F4mWEpAE1u4GISgBiURKCUgCi3MQQr41wYlI85SAJDgnnhoQFtsghLCj4ESkefr1TgKL7UHUFkxGuj2xnfc+fS9U2eraatWARCKk/10SmBNPH1BLngO6/uXrmfPOnNDX7tahW+iyItI0JSAJLOnxjIJrySCEddvXsW/3ffluxXezvy5GRb/sF4UTkWCUgCQr+bYcw47aHfTp0odjBh+T46hEpKU0CEECi2sYdkv6gKprq0PNxyYi0VMCksDycUnuHbU7Qs1ILSLRUwKSwPJxQTrVgETaLiUgCczdY5uMNOyDqDtqdygBibRRSkASWFzDsFsyG3Z1bbWa4ETaKI2Ck8Dc43sQden6pZz4xxOzLrtu2zrVgETaKCUgCSyuQQjTR06nR6ceocpO3GciX9jvCzmOSERyQQlIAotrGPYRA4/giIFHtPp1RSRakfYBmdlUM1tuZivM7KoG9h9lZp+a2aL0zzVRxiMtE9coOBEpTJHVgMysGLgLOA6oBOaZ2Rx3f6veoX9395OjikNyK45RcCJSmKL8NpkMrHD3d919B/AgMC3C60nE4hqEICKFKcoEVA6szPhcmd5W3yFmttjM/mJmB0QYj7RQ0pMo/4hIrkQ5CKGhr6r6w6heA/Z1981mdhLwZ2D4HicyuwC4AGDw4MG5jlMCUh+QiORSlDWgSmBQxueBwIeZB7j7RnffnH7/BNDBzPrUP5G7z3L3Se4+qaysLMKQpSlKQCKSS1EmoHnAcDMbamYdgenAbiuDmVk/S4/rNbPJ6XjWRRiTtERMS3KLSGGKrAnO3RNmNgN4EigGZrv7EjO7KL1/JvAl4GIzSwDbgOkeds4ViZxqQCKSS5E+iJpuVnui3raZGe9/Dvw8yhgkd+KaC05ECpMe6pDANAxbRHJJCUgCUw1IRHJJc8G1Mx9s/IDX17weurxqQCKSK0pA7cwPX/ohr1a9Grp8r869chiNiLRnSkDtzNaarVT0reD6Q6/PumxRURED9hoQQVQi0h4pAbUztV5Ltw7dGNR9UPMHi4hESIMQ2pmEJygp0u8dIhI/JaB2JpFMUGzFcYchIqIE1N7UJmspLlICEpH4KQG1M4lkgg5FHeIOQ0RECai9Sbia4ESkbWi0N9rM5jS2L8N6d/9a7sKRqNUmazUIQUTahKa+iUYB5zWx34C7chuORE01IBFpK5pKQP/h7s83VdjMfpDjeCRiqgGJSFvRaB+Qu/+hucJBjpG2JZHUc0Ai0jaEGoRgZhfkOhBpHWqCE5G2IuwoOE2JnKfUBCcibUWobyJ3/2WuA5Hg1m5by9ufvJ11OXfHcT2IKiJtQrMJyMyuaWi7u/8w9+FIEFe9cBWvVL0Sunz3jt1zGI2ISDhBakBbMt53Bk4GlkYTjgSxpWYL4/qM4/JJl2ddtsiKOKDPARFEJSKSnWYTkLv/JPOzmf0YCPKQqkTEcbp36s6BfQ+MOxQRkdDC9AF1BfYLcqCZTQVuB4qBe939pkaOOwh4GTjT3R8JEZOISKtZsGDBPiUlJfcCY9CUZo1JAm8mEonzKioqVjd0QJA+oDcAT38sBsqAZvt/zKyY1EwJxwGVwDwzm+PubzVw3M3Ak82dU1IcxzQQUSQ2JSUl9/br129UWVnZJ0VFRd58ifYnmUzamjVrRldVVd0LnNLQMUFqQCdnvE8AH7t7IkC5ycAKd38XwMweBKYBb9U77lLgj8BBAc4pItIWjFHyaVpRUZGXlZV9WlVVNaaxY4L0AX0Q8vrlwMqMz5XAwZkHmFk5cCpwDE0koPSDrxcADB48OGQ4hcPdMVMNSCRGRUo+zUv/GTXaRBl2JoS5QQ5rYFv9v7CfAVe6e21TJ3L3We4+yd0nlZWVBQ2zoKkJTkTyXdjOs/MDHFMJDMr4PBD4sN4xk4AHzex94EvAL8zsiyFjEhGRtDvuuKP32Wef3WCTUdeuXSe2djwNCTsTwkcBDpsHDDezocAqYDrw7/XOM7TuvZndD8x19z+Hiak90SAEESkEQUbBDQduBEaTehAVAHdvcii2uyfMbAap0W3FwGx3X2JmF6X3z2xJ4CIi7dnnPve5YR999FHH6urqoosuuujjK664Yu3tt9/e+7bbbutfVlZWM2zYsO0dO3Z0gGXLlnWcPn36folEwo499thP4469TpAa0K+Aa4HbgKOBcwk4Gam7PwE8UW9bg4lHK6sG5+6aDlakjfjeI4sH/bNqU9dcnnNEv9Ktt35p/Mqmjvnd7373ft++fWs3b95sEydOHH366ad/etNNNw1YsGDB0l69etUeeuih+48ZM2YrwCWXXDL4vPPOWzNjxox1N954Y5vpSA/SB9TF3Z8FzN0/cPfrSI1aExGRmNx88819999//9EVFRWjqqqqOtxzzz29p0yZsmnAgAGJzp07+2mnnba+7tjXXnut2/nnn78e4MILL1wXX9S7C1ID2m5mRcDb6Sa1VcA+0YYlTVEfkEjb0VxNJQpz584tff7550vnz5+/rLS0NDl58uT9R40atX3ZsmWdGyvTFoeNB6kBfZvU9DvfAiqAs4BzogxKREQat2HDhuIePXrUlpaWJhcuXNh58eLFe23durXo5ZdfLq2qqiqurq62Rx99tGfd8QceeODme+65pxfAPffc0zu+yHfXbAJy93nuvtndK939XHc/3d1fbo3gpGGqAYm0b6effvqniUTCRowYMfr73//+gPHjx28pLy+vufLKKz+cMmXKqMMOO2zEuHHjttYd/4tf/OJfs2bN2mfMmDGjPv300zazIFijTXBmdl26v6dRQY4REZHc6tKli7/wwgsNrkp52WWX7dHHM3LkyB2LFi1aVvf5Rz/6UVWU8QXVVB/QeWa2sYn9RurZnutyGpE0S1PxiEghaCoB3QOUNlP+nhzGIiIi7UijCcjdf9CagUh21AckIvlOCymJiEgslIDykPqARKQQNJuAzKxXawQiIiLtS5Aa0Ctm9rCZnWT6tbtN8D2WVRKR9uy73/3ugGuuuaZv3HFkK0gCGgHMAr4KrDCzH5nZiGjDkuZoEIKINKWmpibuEJoVZCYEd/en3f3LwHmkpuF51cyeN7NDIo9Q9qAakIhceeWV/YYMGTLm0EMPHfH22293Apg8efL+M2bMKD/ooIP2v+GGG/o+9thjpaNGjRo9YsSI0WecccaQbdu2GUB5efnYiy++uHzs2LGjxo4dO+rNN9/sFMc9BFkPqDep+d++CnwMXArMASYADwNDGy8tUdAgBJE25M/fHMTqt3K6HAP7jN7KF+9qdJLTv//9710fffTRXm+88cZbNTU1TJgwYfTEiRO3QmqeuHnz5i3funWr7bfffmOfeuqp5ePGjas+9dRTh9x6661l11xzzWqA7t27177xxhtLf/7zn/e+9NJLBz333HMrcnoPAQRpgnsJ6A580d0/7+5/cveEu88HtKiciEgre+6557qddNJJG0pLS5O9evVKHn/88Rvq9n35y19eD7B48eLOAwcOrB43blw1wNe+9rV1L7744s7JBc4555z1AOeff/76hQsXdmvte4BgyzH8p7v/IXODmZ3h7g+7+80RxSXNUB+QSBvRRE0lSo21gpSWliYhvXBlE4qKdtU/zCyWdv0gNaCrGth2da4DERGRYI455pjNjz/++N6bN2+2Tz75pOjpp5/eu/4xEyZM2L5q1aqOdf07DzzwQO/DDz98U93+Bx54oBfAfffd13PixIlbWi/6XZqaDftE4CSg3MzuyNjVHUhEHZg0TssxiLRvhx122NZTTz11/ZgxYw4oLy+vnjx58ub6x3Tt2tVnzpz5/hlnnDGstraW8ePHb73iiivW1O2vrq62cePGjUwmk/bggw++27p3kGKNVdPMbDypgQY/BK7J2LUJeM7dP4k+vD1NmjTJ58+fH8el24yTHz2Z0b1Gc8uRt8QdikjBM7MF7j4pc9vixYvfHz9+/Nq4Ymqp8vLysfPnz1/av3//yCsTixcv7jN+/PghDe1rajLSxcBiM/udu6vG04a4O6oAiUi+a7QPyMzqBh4sNLPXM37eMLPXg5zczKaa2XIzW2Fme/Qlmdm09DkXmdl8Mzss5H2IiEhAq1ateqM1aj/NaWoU3GXp15PDnNjMioG7gOOASmCemc1x97cyDnsWmOPubmbjgD8AI8Ncrz1RH5CIFIJGa0Du/lH67Vpgpbt/AHQCxgMfBjj3ZGCFu7/r7juAB4Fp9a6x2Xd1Qu0FesRfRKS9CDIM+wWgs5mVk6qxnAvcH6BcOZA5Pr4yvW03ZnaqmS0DHge+3tCJzOyCdBPd/DVr1jR0SLuimRBEpBAESUDm7luB04A73f1UYHSQcg1s26OG4+6PuvtI4IvA9Q2dyN1nufskd59UVlYW4NIiItLWBUpA6UlHv0KqlgLBZlCoBAZlfB5IE0137v4CMMzM+gQ4d7umPiCR9m3t2rXFN910UxnA3LlzS48++ujPtOb1c3XNIAnoMlIzHzzq7kvMbD/guQDl5gHDzWyomXUEppOaxHQnM/tM3RpDZnYg0BFYl80NiIi0N+vWrSu+77779smmTCIR+6C3PQRZjuEFdz+lbt639KCCbwUolwBmAE8CS4E/pBPYRWZ2Ufqw04E3zWwRqRFzZ2YMSpAmqAYk0n5dfvnlA1euXNlp5MiRo6+66qqBW7ZsKZ46dep+Q4cOPeCUU04ZmkwmgdQDp1dccUX/ioqK/WfPnt3zl7/8Za8RI0aMHj58+AEXX3zxzj75rl27Tqx7/6tf/arn6aefPgRgyZIlncaPHz9yzJgxo7797W8PyDyusWtmI8hyDCOAK4Ahmce7+zHNlXX3J4An6m2bmfH+ZkATmoagQQgibcN//eO/Bq34ZEVOl2P4TM/PbL3+s9c3OsnpT37yk8qTTz65y7Jly96aO3du6Ze//OVhixYtenfIkCE1FRUVI59++uluJ5xwwmaAzp07JxcsWLD8/fff73DIIYeMXLBgwdKysrLE4YcfPuI3v/nN3l/96lc3NHadGTNmDLrkkktWX3jhhetvueWW3Trgly5d2qWxawYVpAnuYWAh8J/A9zJ+JCaqJIpIprFjx24ZNmxYTXFxMQcccMDWd955p2PdvrPPPvsTgBdffHGvKVOmbBowYECiQ4cOnHnmmeuff/75JpdhWLhwYbevf/3r6wHOO++83bpHmrpmUEEGEyTc/e5sTyzR0YqoIm1HUzWV1tKpU6edXwrFxcUkEomdTSRBlmfIbFGpWzW1JdcMKkgN6L/N7BIz629mvep+sr2QiIjkRo8ePWq3bNkS5Pt7pyOOOGLLK6+8UvrRRx+VJBIJHn744V5HHXXUZoDevXvXvPbaa51ra2t57LHHetaVmTBhwub777+/J8Ds2bNz/r0fpAZ0Tvo1s9nNgf1yHYwEo2HYIu1bv379aisqKjYPHz78gE6dOiXLyspqmiuz77771lxzzTWrjjzyyBHubscee+ynZ5111gaAH/zgB6umTZv2mf79+9eMHDlyW11yu/POO1d+5StfGXrHHXf0O/744zd069atNpf30ehyDG2VlmOA4x45joP7HcwNh90QdygiBa8Ql2MIatOmTUV77bVXsqioiFmzZvV86KGHej377LPvZHOOUMsx1DGzrsB3gcHufoGZDQf2d/e52QQhuaOpeESkNfzjH//oetlllw12d7p37157//33v5/L8wdpgvsVsAA4NP25ktTIOCUgEZECNnXq1M3Lly9/q/kjwwnSiTXM3W8BagDcfRtaDi1W6gMSiV0ymUzqP2Ez0n9GjT6hGiQB7TCzLqQnEjWzYUB1bsITEclLb65Zs6aHklDjksmkrVmzpgfwZmPHBGmCuw74KzDIzH4HfJbUkgwSF9dMCCJxSiQS51VVVd1bVVU1hmC/yLdHSeDNRCJxXmMHNJuA3P0pM1sATCHV9HaZuxf86A8RkcZUVFSsBk6JO45812zmNrNn3X2duz/u7nPdfa2ZPdsawUnD1AckIoWg0RqQmXUGugJ9zKwnuwYedAcGtEJsIiJSwJpqgrsQ+DapZLOAXQloI6mlEyQmmgtORApBownI3W8HbjezS939zlaMSURE2oEggxDuNLND2XM9oAcijEuaoJkQRKQQBJmK5zfAMGARUDcRnQNKQDHSIAQRyXdBngOaBIzWUtlth/qARKQQBHmA6k2gX9SBSHZUAxKRfBekBtQHeMvMXiVjCh5310NYIiISWtCpeKQN0SAEESkEQUbBPR/25GY2FbgdKAbudfeb6u3/CnBl+uNm4GJ3Xxz2eiIikj+amglhEzTY222Au3v3pk5sZsWkHlg9jtQaQvPMbI67Z64t8R5wpLt/YmYnArOAg7O8h3ZHgxBEpBA09SBqaQvPPRlY4e7vApjZg8A0YGcCcvf/zTj+ZWBgC68pIiJ5IsppxMuBlRmfK9PbGvMN4C8N7TCzC8xsvpnNX7NmTQ5DzE+ajFRECkGUCaihb8gG247M7GhSCejKhva7+yx3n+Tuk8rKynIYooiIxCXIKLiwKoFBGZ8HAh/WP8jMxgH3Aie6+7oI4ykYGgUnIoUgyhrQPGC4mQ01s47AdGBO5gFmNhj4E/BVd/9nhLGIiEgbE1kNyN0TZjYDeJLUMOzZ7r7EzC5K758JXAP0Bn6R/o0+4e6TooqpUKgPSEQKQZRNcLj7E8AT9bbNzHh/HtDoeuEiIlK4omyCk6g46gMSkbynBCQiIrFQAspD6gMSkUKgBCQiIrFQAspDmgtORAqBElAe0oOoIlIIlIBERCQWSkB5SIMQRKQQKAGJiEgslIDylGpAIpLvlIBERCQWSkB5SKPgRKQQKAGJiEgslIDykEbBiUghUAISEZFYKAHlIXdHFSARyXdKQCIiEgsloDykPiARKQRKQCIiEgsloDykGpCIFAIloDylB1FFJN9FmoDMbKqZLTezFWZ2VQP7R5rZS2ZWbWZXRBlLQdF6dCJSAEqiOrGZFQN3AccBlcA8M5vj7m9lHLYe+BbwxajiKERqghORQhBlDWgysMLd33X3HcCDwLTMA9x9tbvPA2oijENERNqgKBNQObAy43NlelvWzOwCM5tvZvPXrFmTk+DymasNTkQKQJQJqKE2olDfnO4+y90nufuksrKyFoYlIiJtQZQJqBIYlPF5IPBhhNdrN7Qcg4gUgigT0DxguJkNNbOOwHRgToTXExGRPBLZKDh3T5jZDOBJoBiY7e5LzOyi9P6ZZtYPmA90B5Jm9m1gtLtvjCquQqBRcCJSCCJLQADu/gTwRL1tMzPeV5FqmhMRkXYm0gQk0YmkD2jrevjvy2DHlnDl+4yAE2/KbUwiUrCUgGSXDxfC0jlQNgo67pVd2U0fwTvPwvHXQ3GHaOITkYKiBJRn3FMj2SPpA6rZlno9bRb0H5dd2Rdvg2eug9oaJSARCUSTkcoudQmoQ9fsyxalk05Sk1qISDBKQHmmbhaEaGpAW1OvHTpnX7a4Y+q1NpG7eESkoCkByS4tqQEVp1tza3fkLh4RKWhKQHmmrg8okseAEnUJqEv2ZdUEJyJZ0iCEPNVoE9zaFTD7eNixNfuTJmvAiqCkJU1wSkAiEowSUJ5pdibs9e/C1nUw7kzo1jf7C5TtD2GeMdrZBKcE1KYlqmHtP8FDzAvcqRR6Dc19TNJuKQHlmWYHISTTgwCmXAwDJrZSVKgJrrW4w8LfwuaP6zbs9rL7NmePna8/lPolJawjr0w9cAwZv6hk/Ftsalt5BfQIMfHJpo9h04fp+8m8X9/9HoO879YP+nwm+xgkEkpAhaYuARW18rM4dc/+5FMNKFENT18L2zeEK1/aH469JlyNMay1b8OcGeHLl3SBcdNh1BeyK7dlNcz9Djx/c/hrQ/a/FCUTUPVGy66ZqeJc+MLPcnc+aREloHxTNwahsS+9nQmolf9q40pA7vDu/8D2T7Mvu/ZteOXuVCLJ9uHZHVtSTZ0VX4Oe+2Z/7cUPph7ezVZd3943nob+E1Lvs62JhE2YI0+GbXXJun4tq5ltK1+F5X+BMEuC9dkfhh8PnbsDtvu9mbHzHo2M99bw++4Dsr++REYJqNAka1OvrZ2A6mpc7/8dtq7Nvnz3chgwIftyq5fCb76Yfbk6JZ3hkpegS8/syi3/C/x+OmxZGy4BLZubmr5ov6OyLzv8uFRNorVnnOi2T+onjL4HwKRzcxuP5D0loDwTuA+oqLiVIkrbq0/q9W/Xhytf3BGuroSSTtmV27ou9TrtLhhwYPbX7dor++QDsFf6i/iZa6G0X/bll/437PtZ+LcHsi8rUiCUgApN3SCA1q4B9T0ALn0NdmzOvuySP8OLP001o2X7G3bdzN37jIa+o7O/dlhlI6B8EmxclfrJxtb1qde6pC3STikB5ZmdNaC21gcE0HtYuHJrlqdet28MkYDSCa9jt3DXDqtTKZz/bLiy7qkaUHmIGptIAVECKjR1fUD5NCN1p9LUa3WIgQR1NaBsl4+IkxmMPiXuKERipwSUZ7y5Bwjj6gNqic49Uq+/PX3XjApB1Y0Ky6cEJCKAElDhibMJLqwBB8Khl0L1pnDl994Xuuyd25hEJHJ59C0lkM0ouDz6q+3QGY6/Ie4oRKSV5dG3lASSjwkogE+27CCRzP4hxiZ5TKEAAAlpSURBVG07allcuYFEMhnqusP3KWVIn+yb95LufLRhO7UhYgbo16MzvfbKsjkyB1Zv2k51ze5/Vu67fvFJva/b7hnv4e2PN/H26hCjIIFhZd34/Lj+IaOWfBXpt5SZTQVuB4qBe939pnr7Lb3/JGAr8DV3fy3KmPLdziW5GxsFV7cgnOVRH1AzHl1YyXceWhx3GK3uW8cOpyjEpAU9u3akR5fsB6G89M46Hpq/MvsL5sj/fbwzJcVFOJ5KdBm5uy7Z1SXDukTo6Tnhdn1OJ8q6Y+uSZHr/GZMGcu0XDmjdG5NGRZaAzKwYuAs4DqgE5pnZHHd/K+OwE4Hh6Z+DgbvTrxJWMpFaUqGocJZ6WrJqI51KivjPk8M95zO4V1f27ZX9InubtidY8MF6amrD1WL27tqB0s7ZJ4JX31vPr196nzuefTvUdVvqptPGUpzOfHW/6Bi7z+RT1wSc+XtQx+IiDh9RRpcO2f3ys25LNbf+dTk7apPp69jOWXUyr2PsurbtnJHHGty3K+bM8kbFviEeOpbIWLOjqsKe2OwQ4Dp3PyH9+WoAd78x45hfAv/j7r9Pf14OHOXuHzV23kmTJvn8+fMDx3HV7FN4vbYFs/+2MW5Q2cE49xPntI177u/pG+ji2/h89z+1fnARqdq4nbJunfjbFUfFHUqrCfv/0h1WbdgWqrkSoHe3jnQPkTQLmZktcPdJccdRiKJsgisHMuvzlexZu2nomHJgtwRkZhcAFwAMHjw4qyD27lxGv51T1xcAh/Ltxgjvwfque/YRrAcqOw1jeK9WfjAzQsP7duPo/UPOQZanGm1ibbYcDApR2xOJQ5QJqKH/QfV/LQtyDO4+C5gFqRpQNkFc9e/3ZXN4QTgQ0GOOItLWRdlRUAkMyvg8EPgwxDEiIlKAokxA84DhZjbUzDoC04E59Y6ZA5xtKVOAT5vq/xERkcIRWROcuyfMbAbwJKlh2LPdfYmZXZTePxN4gtQQ7BWkhmFrwRARkXYi0ueA3P0JUkkmc9vMjPcOfDPKGEREpG0qnIdFREQkrygBiYhILJSAREQkFkpAIiISi8im4omKma0BPsiyWB9gbQThxEX307YV0v0U0r1AuPvZ193Logimvcu7BBSGmc0vpLmcdD9tWyHdTyHdCxTe/eQ7NcGJiEgslIBERCQW7SUBzYo7gBzT/bRthXQ/hXQvUHj3k9faRR+QiIi0Pe2lBiQiIm2MEpCIiMSi4BOQmU01s+VmtsLMroo7noaY2SAze87MlprZEjO7LL29l5k9bWZvp197ZpS5On1Py83shIztFWb2RnrfHRZ2ac2W31OxmS00s7n5fi/pWPY2s0fMbFn67+mQfL0nM/tO+t/Zm2b2ezPrnE/3YmazzWy1mb2ZsS1n8ZtZJzN7KL39FTMb0hr31S65e8H+kFoG4h1gP6AjsBgYHXdcDcTZHzgw/b4U+CcwGrgFuCq9/Srg5vT70el76QQMTd9jcXrfq8AhpFab/QtwYkz39F3g/wFz05/z9l7SsfwaOC/9viOwdz7eE6kl798DuqQ//wH4Wj7dC3AEqYV/38zYlrP4gUuAmen304GH4vp3V+g/hV4DmgyscPd33X0H8CAwLeaY9uDuH7n7a+n3m4ClpL4oppH64iP9+sX0+2nAg+5e7e7vkVpPabKZ9Qe6u/tLnvrf80BGmVZjZgOBzwP3ZmzOy3sBMLPupL707gNw9x3uvoH8vacSoIuZlQBdSa1CnDf34u4vAOvrbc5l/JnnegQ4Ns7adyEr9ARUDqzM+FyZ3tZmpav7E4FXgL6eXiE2/bpP+rDG7qs8/b7+9tb2M+D/AMmMbfl6L5CqQa8BfpVuVrzXzPYiD+/J3VcBPwb+BXxEahXip8jDe6knl/HvLOPuCeBToHdkkbdjhZ6AGvqtpc2OOzezbsAfgW+7+8amDm1gmzexvdWY2cnAandfELRIA9vaxL1kKCHV5HO3u08EtpBq5mlMm72ndN/INFLNUQOAvczsrKaKNLCtTdxLQGHiz5d7y3uFnoAqgUEZnweSam5oc8ysA6nk8zt3/1N688fppgLSr6vT2xu7r8r0+/rbW9NngVPM7H1STZ7HmNlvyc97qVMJVLr7K+nPj5BKSPl4T58D3nP3Ne5eA/wJOJT8vJdMuYx/Z5l0M2UP9mzykxwo9AQ0DxhuZkPNrCOpDsU5Mce0h3T78n3AUnf/acauOcA56ffnAI9lbJ+eHq0zFBgOvJpuethkZlPS5zw7o0yrcPer3X2guw8h9ef9N3c/Kx/vpY67VwErzWz/9KZjgbfIz3v6FzDFzLqmYziWVJ9jPt5LplzGn3muL5H6N6waUBTiHgUR9Q9wEqlRZe8A/xF3PI3EeBipKv7rwKL0z0mk2p2fBd5Ov/bKKPMf6XtaTsboI2AS8GZ6389Jz3YR030dxa5RcPl+LxOA+em/oz8DPfP1noAfAMvScfyG1AixvLkX4Pek+q9qSNVWvpHL+IHOwMOkBiy8CuwX17+7Qv/RVDwiIhKLQm+CExGRNkoJSEREYqEEJCIisVACEhGRWCgBiYhILJSAREQkFkpA0qaYWW8zW5T+qTKzVen3m83sFxFc734ze8/MLkp/LktPwb/QzA5vpMyt6diuyHU8Iu1JSdwBiGRy93WkHvrEzK4DNrv7jyO+7Pfc/ZH0+2OBZe5+TmMHu/v3zGxLxDGJFDzVgCQvmNlRtmtxu+vM7Ndm9pSZvW9mp5nZLenFxf6anlevbsGx581sgZk9WTdXWBPXmEBqXZmT0rWuvdI1pDfT5/5O9Hcq0n4oAUm+GkZqzaFpwG+B59x9LLAN+Hw6Cd0JfMndK4DZwP9t6oTuvgi4htQCZBOAkUC5u49Jn/tXkd2NSDukJjjJV39x9xoze4PUyrd/TW9/AxgC7A+MAZ5OryVWTGr+sGy8C+xnZncCjwNP5SBuEUlTApJ8VQ3g7kkzq/FdkxomSf27NmCJux8S9gLu/omZjQdOAL4J/Bvw9ZaFLSJ11AQnhWo5UGZmh0BqvSUzOyCbE5hZH6DI3f8I/BepNYBEJEdUA5KC5O47zOxLwB1m1oPUv/WfAUuyOE05qWW4635RuzrHYYq0a1qOQdo1M7uf1JpFjzR3bL1y19E6Q8RFCpaa4KS9+xS4vu5B1CDM7FbgLEDPAom0gGpAIiISC9WAREQkFkpAIiISCyUgERGJhRKQiIjE4v8Dv6GvEx2te8UAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "with pt.Environment(wl=np.mean(wavelengths), t=time):\n", " detected = nw(source=1)\n", " nw.plot(detected);" ] }, { "cell_type": "markdown", "metadata": { "slideshow": { "slide_type": "slide" } }, "source": [ "## Simulate Frequency Domain" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaAAAAEGCAYAAAAjc0GqAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+WH4yJAAAgAElEQVR4nOy9e5xdZXn3/b33ec/sPcdMMsnkMElIyJEEAxR5QAH7KCoKijyKx+qDihaLb6WPvD6VVnuCWttK0RfRgpW2WmmhHC3YFlFRlIQEEkJCjiSZZDLnmT37vNa63z/utdZee447k32YSe7v5zOfzF577bXumST7t6/r+t3XJaSUaDQajUZTbXy1XoBGo9Fozk60AGk0Go2mJmgB0mg0Gk1N0AKk0Wg0mpqgBUij0Wg0NSFQ6wWcKvPmzZOdnZ21XoZGozlL2LZtW5+Usq3W6zgTmXMC1NnZydatW2u9DI1Gc5YghHi91ms4U9EpOI1Go9HUBC1AGo1Go6kJWoA0Go1GUxO0AGk0Go2mJmgB0mg0Gk1N0AKk0Wg0mpqgBUij0Wg0NUELkOaMQkrJYwceoy/dV+ulaDSaadACpDmj2HpyK1/6xZe468W7ar0UjUYzDVqANHOanlQPOTPnPj40fAiAA0MHarUkjUZTIlqANHOW7mQ3b3nwLXxrx7eKjgEgarQojUZTMlqANHOWV/pfAeChfQ+5xwazgwCM5kZrsiaNRlM6WoA0cxYn2vH7/O6xZC4JwEhu5NQv+NT/hT9uLMvaNBrN9GgB0sxZ+tP9APhE4Z/xaF5FPsPZYaSUAOwf3M/n/vtzHB05OvUFf3W3+tN+nUajqSwVFSAhxFVCiL1CiP1CiNsmOedyIcQOIcQrQohnK7kezZlFMq+inUQuMe5Y3sqTNtIAPHbwMX569Kc89fpTpV3YyJZ3oRqNZkIqNg9ICOEHvgn8T+AY8IIQ4lEp5W7POU3At4CrpJRHhBDzK7UezZmHIzZpI03ezBP0B90ICFQari5Yx8nUSQBOjJ4o7cL5FAQjZV+vRqMpppIR0EXAfinlQSllDvghcM2Ycz4IPCSlPAIgpeyp4Ho0ZxiOAAEM54bdY9FAVB3LqmNDmSEABjIDpV1YGxg0mqpQSQHqALxJ92P2MS+rgWYhxE+FENuEEB+d6EJCiE8JIbYKIbb29vZWaLma2c6BoQNFIlIU7WRH3GML6xeqY7YRwUnROSI1Lbnk9OdoNJrTppICNNFOjLHV3QCwBXgn8Dbgy0KI1eNeJOW9UsoLpJQXtLXp0exnI8PZYa595Fr+4Nk/cI8l80mE/c9sOKdMB6O5Udrr293noSBETkQ0LblUGVeu0Wgmo5ICdAxY4nm8GDg+wTn/IaVMSin7gJ8Bmyq4Js0cxels8Jvu37jHkvkki2KLABUBZcwMpjRZULcAKERIpyxAllGuZWs0mimopAC9AKwSQiwXQoSADwCPjjnnEeAyIURACFEH/BbwagXXpJmjnEiONxCM5keLxMaJeJwIaDQ3ipTSTc+VvDdIC5BGUxUq5oKTUhpCiJuBpwA/cJ+U8hUhxE328/dIKV8VQvwH8DJgAd+VUu6q1Jo0cxdv7ce0TPw+P8l80hWgZD7pdj/wilLaSGNIg3gwTiKfIGfmCPlDU9/Mylfmh9BoNEVUTIAApJRPAk+OOXbPmMdfA75WyXVo5j5j7dWN4UYlQPVKbBK5hBsBtURaCIgAyXzSjXoWxxfz6sCrjORGmBedN/XNdASk0VQF3QlBMydIepxpg9lBd5Npa6QVv1DRkCNSsVCMWChGIpcoEiAosQ5kmWVevUajmYiKRkAaTbnwRkDD2WHqA/UA1IfqqQ/WM5ofLQhQMEZ9sF5FQHb9Z1H9onHXmRRTp+A0mmqgIyDNnMArHIOZQTfdFgvGiAVjJPPJomPxUJzR/KgbATluuZK6ZOsUnEZTFbQAaeYEo/lRWiItAAxlh1yxqQ/WUx+qZzQ36oqLGxXlRt1NqI4AJfKJCa5OcQNSLUAaTVXQAqSZE4zmRt00WiKXKEq3TRQBOcdKjoC8aTddA9JoqoIWIM2cIJlPMr9uPgJR5HirDxbXgIK+ICF/aJwJwbs3aEI8Y71LsWEPZgZ536Pv4+Xel0/vB9NozmK0AGnmBKP5UeKhuCssTgRUH6wvioDqg8qc4EZA2RHiwTixYEyJ12QpuCIBmj4F99OjP2Xv4F7+8oW/PO2fTaM5W9ECpJkTJHNJYqEY8WC8qOuB43hzDAcNoQYA91gil6Ah3IBP+IgFYyVGQNML0PGk6ioV9odP7wfTaM5itABpZj2WtBjNj1IfrCceijOSGylKwXmjHUeA4qE4eStPb7rXPRYLxSa3YXsFyJxegHpSanJIb1p3Z9doZooWIM2sJ22kkUjiwTjxUNxNwYV8IYL+IPWhetJGmoHMAA3hQgQEcDRxlOZIM2AL0GQRkHFqEVBfug84hQanGo1mHFqANLMer73aEREntQYqDQeqYWljqLHoWNdoF03hJgA3fTchp5iCcxucZkeQcuyUEY1GUwpagDSzHm+9pyHUQCKXIJFLEA/FAVyBGcoOjRMloCgCcvYFjcPMFr6fQIC+u/O7PLTvIfex464zpOG2BbKkxX8f+W9yXjHTaDSTogVIMyt5+vDTfPVXXwUKm0edGtBYAfI2F/XWexyaw7YABaeqAXn3ARULUNpI840Xv8Ef/fKP3GjHO9rB+f75489zyzO38Dfb/uaUf16N5mxEC5BmVvKFZ7/Ag689yInRE24j0ngo7orIcHZ4QgFqDKsUXGuk1T3mRkBTueCMySOgo4nCZHlnLMRIdoTFMdXg1BGg3QO7AXh95PVT+2E1mrMULUCaWc3hkcNFe37ioTgSSXeym4aginYmioDa6gqj250UXSwUI5FPTFyzmaIGdHy0MMi3a7SLjJEhZ+XoiHcAhXpQf7ofUOPBNRrN9GgB0sw6nJoPKJuzt+2OIzD9mX43AnKiHoCFsYXuuQ5L4moyfDwUx7AMst56j0ORDbu4E4IjLKAiIHfEw5gIaCg7BGhnnEZTKlqANLOOk6mT7vc9qZ6CCy5YXyQ2jgD5ROGfsSMKQggu7bgUgOWNy4GCKPVn+vnR3h8VmwWKIqDiXnCOsIBqweNEPM6MIS1AGs3M0POANLOOk8mCAPWmeoucbc4EVCgIEMDvbv5ddvXtoiPW4R678013YlgGkUBEvd42Jvz11r/m6defxi/8XLf6OnXyFPuAvCk1bwTk3MsRSNeanRvBklaRMGo0mvFoAdLMOvozKuXlEz56Uj0IIYgFY/h9ftrr2t3zvLWfmzbdNO46TrrOwakFbe/ZDsC+oX2FJ6doRjqcHWZ+dD4juREVAdkCtLBepfscl54TAVnSIplPFgmkRqMZj/6Ippl1OBHFisYV9KX7iizXzkwggPl180/pum1RZUxw2ud0J7sLT7p1ITEuAhrKqP1FzZHmogioJdJCNBB11zuUHSIaiAI6DafRlIIWIM2sw9ksurxxOX3pvqImo0II97zOhs5Tuq43YgKV3nNxjAeh+glrQE3hJiVA2QE31dYQanBt4aZlksglWBpfChTvE9JoNBOjBUgz60jkEoT9YRbWL6Q/069GKnjSWV/6rS9xxZIrXMdbqXijp4AvUCwSTgouWDe+BpQddgVoKDPEUHYIgSgaDzGSG0EiWdqw1H2NRqOZmooKkBDiKiHEXiHEfiHEbRM8f7kQYlgIscP+ur2S69HMDUZyI8SCMeZF55E20pxMnSwSoBvW3MBdV951ytcVQrCycSUAly++vFiAnI2oobpxNuzh3DCN4Uaaw80MZYfclj9+n1/1l8uNuvUfx/Ltbfmje8VpNBNTMQESQviBbwJvB9YBNwgh1k1w6s+llJvtr69Waj2auYNT83FSZl2jXWUr6D/wjgd46N0PsaxhWXEjUUd0xkRAUkqGskM0hhtpCjcxmBlkODtcaHoaKnRmgII129m7dGL0BJf/6HL+4/B/lGX9Gs2ZRCUjoIuA/VLKg1LKHPBD4JoK3k9zhjCaH6Uh1EBrtNBOZ6yjbabEQ3FWNa+iIdxQ1EgUMwu+IPgCRTWglJHCsAyawk20RFpIGSl6Uj2F7gpBlYJzBcjeh+REQL88/ksGMgP84NUflGX9Gs2ZRCUFqAM46nl8zD42ljcKIV4SQvxYCLG+guvRzDYGDsJw17jDiVyCWChWZBrw1m/KgSNobhrOzEMgbAtQIQXnCEtjuJGmiBKdwyOH3Q2x8ZAa8eCk4BbFFqnR37YAOX3hjBJGPGg0ZxuVFCAxwbGxyfAXgWVSyk3A3wH/PuGFhPiUEGKrEGJrb6+eQHnGcNf58Dfjs7KJXMKtATmcquV6OhwBcs0CRhb8TgRUEAtHWJwaEKjNqN4IyFsDao40u+PAoTC4ztnbpNFoClRSgI4BSzyPFwPHvSdIKUeklKP2908CQSFEsVdWPXevlPICKeUFbW1tY5/WnGE4NSDnTR6gvb59ilecOk4EU4iAcuAPKRHypOAcgXJScGNfHwvFyJgZ+jP9+ISPWDBWNHdoIKu6ZxdZvjUaDVBZAXoBWCWEWC6ECAEfAB71niCEaBf2xg4hxEX2evRHxbOcRC5BQ6gBn/C53QbWtKwp6z0cU4PrVjNz4A+Dz18UATlteBpDjUVtgBbULSi6ztGRozSHm10RcjenZlRklLNyZIxMWX8GjWauU7FWPFJKQwhxM/AU4Afuk1K+IoS4yX7+HuB9wGeEEAaQBj4gtWf17MCYoCM1kDNz5Kyc27ft/qvupyfVU9SEtBzUB+sBT+dtM1dIweVS7nnOplPHBecwtuv264nX3blDTl0IihuZJnIJty+dRqOpcC84O6325Jhj93i+vxu4u5Jr0MxSshOPxnYiEueNvSPWUdRgtFw4ApTK22JjZG0TQrA4AvKYEIL+oHv8nKZzgMKwuwNDB7hwwYXu2p3az0BmgPnR+fSkexjJjRTNKdJoznZ0JwRNbZhGgCrdyLMuUAdA0nAioLzHhFCoATn93UL+EADfuOIbfGjth1jRuAIomCMsabk1ImdvUMbIkDbSbncE7+ZUjUaju2FrakUuOeFhJ3VVaQGKBqL4hM+Tgst6akDFNmxv+u/KpVdy5dIr3cdOg1OAlqgSoLHdEZY2LGXrya26P5xGMwYdAWlqg1eAPBHH2BRcpRBCUB+o9whQXrngxtiwh3OFrgcT4aTgANem7Yz+HsgoB9yyhmWAblCq0YxFC5CmNtguMaCo91q1IiCAumBdsQkh4NiwCwI0kh2Z0gBRNI3VbsPjjP52Bustiy9zr6XRaApoAdLUhqIIqCBA1aoBgTIiuAJkZD0RUPE+oOkceI5JwrGKx4O2NTuhGoEsaVDb4XQEpNEUo2tAmtrgFSBzvAA5NuxKUh+sL7jg3BScv2g9/Zl+toS3THmdv7n8b9jVv4uVTarTtrP2I4kjgNoz5B1cp9FoFFqANLXB+2bsSXmN5kcRiIrXgECl4JyUnzIhhIps2Kl8iqHsEItii6a8ztrWtaxtXes+dqK3wyOHCYiAmhtkD67TaDQFdApOUxuKakA599vR3Cj1wfqi2kqliAVjxSYEtxmpEqCuUdUo9VT3ITnieWTkCM0RuztCKFaUgutL9+k5QZqzHi1AmtowSQpuJDdSlfQbjEnBTdCM9Pioal04XQQ0FicCOpE84e4NiofibgpuR88OrvjRFTxy4JFy/BgazZxFC5CmItz289v44s++OPkJRSaEQgquP9NPa6R1gheUn7pAnWcjaq5QAxoTAZ2qAHk7d7sCFCy05/nV8V8B8J+v/+dprV+jmetoAdKUnVQ+xRMHn+DJQ08WIoyxTGLD7k/3F41hqCSOC05KOaYbdkGAwv7wKQtiQ6jBTcO5m1NDcddgcWjkEACD2cFy/SgazZxEC5Cm7Dj2YygMZBvHJDbsvnRfVQXIsAxyVs4TAdkpOCk5PnpcDZgTE422mhwhhNvzzakfeUc09KR6AOge7S7jT6PRzD20AGnKzsnUSff73vQkc3CKakAq4jAtk4HMQNEo7kridsTOJpToOCYEAMuka7Rrxo1Qc7axYnnjckCl4BwBcmYDDWWHtBFBc1ajBUhTdrzD1yYdxJb12rBVBDSYHcSSVlUjIICkMzLBMSEAWAbdyW53HtGp8n8u/D+saFzBZR2XASoFl7NyZM2sK8o5K0faSJ/eD6HRzGG0AGlOm5HcCK/2v+o+9tY2etI9E78oNwpB1ZHasWEfSxwDmPGb/qnijmTI2GO5/YUIyDJzDOeGi2YAnQpXLr2SR659pGhyKkBPsoe0kXYjI10H0pzNaAHSnDZf+eVX+F+P/y9XQIazw4T9YVoiLZNHQLkkROw3d9uEsHdgLwDnNp9b8TWD2ogKMJp1BCjkClAyO4wlrbINwnOs2Y4BwZkn5B1Yp9GcbWgB0pw2T7/+NAA7+3YChf5pbdG2qQUoaguQ7TrbO7iXeChOe317xdcMnhScs0HUk4IbsUWpIdRQlns5/eEODh0EYHXzaqAwslujORvRAqQ5bcL+MKBaz4D6VN8UbqIl0uKOJBhHLglRe5SBEwEN7uXc5nNP2XU2UxyrdMppkeMPgt8RICUM5RKghrC6zr6hfQCsal4F6BSc5uxGC5DmtEjlU2TNLKD28EAhAmqNttKf6R//IilVDcgRINuEcGjokNvQsxoUIiBbgHyBQgSUUcLgCMfp4hgrdvfvBmBVkxIgZ+S3RnM2MmkzUiHEoyW8fkBK+TvlW45mruF9A/UK0IqmFZNHQEYWpOmpARlkjAyJfIIFdQuqsWygUANy+8H5AiAtoDC7p1wRkDM5df/QfuqD9W53BT0jSHM2M1U37LXAjVM8L4Bvlnc5mrmGt8GmE+0MZYdoCDXQEmkhbaRJ5VPumz1Q6ILg1oDy9KX7AKpmwQbVigcg6aTgfH6QQaDwc5VLgCKBiNsNYWH9QgK+APFgnOGcjoA0Zy9TCdD/lVI+O9WLhRBfKfN6NHMMJwJqjbTSn+5HSslwVtmXnQ2lA5mBYgEyMupPp+momXMFqFqbUAECvgARf4Sk0y7IF1TpQQqTWcvZGLUt2kYil3Cjn4Zwg07Bac5qJq0BSSl/NN2LSzlHc2bjRAormlYwkBkgZaQwpOGaEIDxaTin+Wgwqv408276rpoRENj94AxHgAo1oJypRNIxWJQDpz2Ps8+pIaQFSHN2MyMTghDiUyWed5UQYq8QYr8Q4rYpzrtQCGEKId43k/VoaofzBrq8YTmj+VG3DU9juNFt4umIi4vpCJAdFVlGTVJw4AiQ3Y3AF1BOOAqtdIK+YNnu5RgPzms7D1C/Iz2mW3M2M1MX3LQ+WSGEH1UjejuwDrhBCLFukvPuBJ6a4Vo0NcSpYaxoWgHAgaEDAK4LDmDPwB4++MQH+a8j/6Ve5DQf9URAbsqrCpNQvaiO2E4E5PdEQDmCvmBZLeE3bbqJP7v0z3h759sB9TvSEZDmbGZGAiSl/HYJp10E7JdSHpRS5oAfAtdMcN7ngH8DJunZopnNjGRHCPgCLI4tBmDfoNrn0hptpTXSil/4+cdX/5GdfTt5YPcD6kXO+IWQEwHlXSt3JBCp6vrrgnUkTU8E5PMDkDOzhPyhst6rMdzIu1e+m6AdZTWGdASkObuZyoQAgBDi9omOSym/Os1LO4CjnsfHgN8ac+0O4D3AlcCFU6zhU8CnAJYuXTrdkjUVZEfPDu55+R5uv/h2FsUWMZwbpiHU4KbOnHY6rZFWgv4gSxuWcmhYtZ9x32ydCCjgREAGGTND0BesyihuL/XBenpH7e7dvgCgTAh5OwKqJE4EJKWcPNIaOgondsDad1V0LRpNLSjlf3vS82WiUmqdJbxuov9RY3vP/y3wRSmlOdWFpJT3SikvkFJe0NbWVsKtNZXiB3t+wHNdz/HM0WcAFQE1hhtZUK/27+wZ2AMUJoEub1juvnbQ3tyJZf91+0MgfGDmyBpZIv7qRj8A9YF6Unb0pVrx2DUgK0/IV94IaCwNoQZMaZIyJhnaB/C3G+BfPuy68zSaM4lpIyAp5de9j4UQfwWUskn1GLDE83gxcHzMORcAP7Q//c0D3iGEMKSU/17C9TU1wDEUOKm24dwwjaFGWiItBESA48njRANR13Z9/vzz+e+j/82KxhW8PvI6lrTwOSk4f0C94Vt5sqZJOFA+x1mphANhMpYyHKh9QOqNPmfm3VRZpXAanQ5nh92uDJNi5iFQWUHUaKrNtAI0AXXAihLOewFYJYRYDnQBHwA+6D1BSul+PBZCfA94XIvP7MbpXeZM9RzJjtBW14ZP+JhXN4/uZLcb/QB8dP1HuXzJ5Tx3/Dnu+M0dDGeHaXZScL6gijpMg4xpldXyXCphf5iM6QhQwBWgvJUvew1oLE6bn+HssLs3aFKMjBagWcS2bdvmBwKB7wIb0C3NJsMCdhmGceOWLVsmrPGXUgPaSSF15gfagOnqP0gpDSHEzSh3mx+4T0r5ihDiJvv5e0r7GTSzCWd8gCtAuRF3tMD8uvl0J7uL2un4hI/Oxk5eHVDzggYyAzQ7Nmyn+7SVJ2sYNUnBRQNRstJej6/w36EaKbjGkB0BldINwchWdC2aUyMQCHy3vb19bVtb26DP59P50QmwLEv09vau6+7u/i7w7onOKSUCutrzvQGclNL5Hzs1UsongSfHHJtQeHRPudmP0+UACmO3R7Ij7if5pfGlvNz7MksbxhtFvHuCVroRUEDVgcw8GTNbmxScP0zGyiMB4RGgvJUnFIpW9N7eFNy06Mmps40NWnymxufzyba2tuHu7u4Nk54z3UWklK97vrpKFR/NmUfaSJM1s9QH6xnKDpHKp0jkE+4n+fPnnw/Ahe3jDY3um21uuNAJwdn4aduwaxEBObbvrBDF+4CsfMVdcE6fudIESEdAswyfFp/psX9Hk+rMTDshPD7jFWnmLM4bpTOxdP/QfqBQy7h+9fU8+K4HedeK8ZZh5802kUsU9gE5KTjTIGtka1IDckRPCVCwsA/IMipeA2qJtCAQ9KZ7eWT/I7zS/8rkJzv98zSaM4iZFs8+WdZVaOYEjgHBmebpWK6bw2qujxCCNS1rJtzT4jT1TOQSngjIMSHkyJiZmrngADJC2BtRHRt25QUo6A/SXt/O1u6t/OFzf8gHHv9A8Qle67WOgDSnyF133dX60Y9+dMKNk3V1dedXez0TMdNOCCfKvRDN7McxIKxuKRagUjpY1wfrEYgxEZDXhl2jFFxRBFRoRpqXRsVTcABL4kvYenKr+7goHee48wDyugakOfOYVoCEEKuEEP8qhNgthDjofFVjcZrZxVBGCZCTgnu1XznbvLbryfAJH7FgTPV8K7JhB9yBdLWyYQOkx9WAjIq74AAWxxcXPe5OdhceeNNuzu9Mo7H57d/+7ZXr169fe84556z/q7/6q3kA3/jGN1o7Ozs3XHjhhef+8pe/dBsr7tmzJ7R58+Y1GzZsWHvLLbdM4/mvHqW44O4H/gj4G+AK4OOU0IxUc+bhREAdsQ6igei4rgfT4QxkI2x3vPaNiYCq3AcOPCYEnx0B2enDvDQrvhEV4A3z38BD+x7iDfPfwIs9L3IieYJzW5TAF6XdTC1As5U/+NeXlrzWnaib/szSWd0eT33tfZuOTnXOP/3TPx1esGCBOTo6Ks4///x111133fAdd9yxaNu2ba+2tLSYl1xyybkbNmxIAXz2s59deuONN/befPPN/X/xF38xa9rJlJKCi0op/wsQthPuj1G92zRnGY4ANYYbWVC3AEMa+ISPpnBTSa+PhWJ2Cs6zD8gfBFMJUO1NCJ55QNKseA0I4F0r38W3/+e3ueOyOwA4kfRkt71pN286TqMB7rzzzgXnnnvuui1btqzt7u4Ofuc732m9+OKLE4sWLTIikYh873vf6w7ievHFF2Of/OQnBwA+/elP909+1epSSgSUEUL4gH32xtIuYH5ll6WZLfSn+5VbSwh31HbAF2B+3XwOjxymKdyE33aOTYcbAUWK9wFJK1+7FJzXhOAPgqUioJw0q5KC8wkflyy6BEtaBHyBYgEqioC0AM1WpotUKsHjjz8ef/bZZ+Nbt27dE4/HrYsuuujctWvXZvbs2TNpGmE22sZLiYA+j2q/83vAFuDDwMcquSjN7GBX3y4u/9Hl/Nu+fwNUDciJdtrr2wG1+bRU4sH4hDbsvJlDImuTgrMjoIwQIPweE0J1IiAHn/DREm5x62xAcQ1Ip+A0HoaGhvyNjY1mPB63tm/fHnnppZfqU6mU7/nnn493d3f7s9msePjhh5ud89/whjeMfuc732kB+M53vlO9uffTUMpG1BeklKNSymNSyo9LKa+TUj5fjcVpastzXc8B8NiBxwCVgmuKKAFa16pmC07U9WAy4qG4bUKwu2HbG1GdXmw1ScHZopfx+cHn86TgrKq44Lw0RZpcqzugIyDNpFx33XXDhmGI1atXr/vSl760aNOmTcmOjo78F7/4xeMXX3zx2ksvvXT1eeed57ZZ/9a3vnXk3nvvnb9hw4a1w8PDpaUsqsCkKTghxB/b9Z5JKeUczdzl8MhhoDDHZyg7RFudql9es/IajiaO8uG1Hy75erFQTF3L24rHFyRr7wuqhQ07YAuO4aQRhcDwBbGQVY2AAJrCTcU2bEPXgDQTE41G5c9+9rN9Ez13yy23jKvxrFmzJrdjx449zuM///M/7x57Ti2YqgZ0oxBiqnGNAtXh+o/LuiLNrGEgo2qYxxLHsKTFUHaIVc2rACUmt1102yldLx6Kk8wnsYwcPsdx5g+QsQWpFhtRA8JOuXnqWLmgWkfVI6BwE68NvlY4oF1wmjOcqQToO0B8mtd/p4xr0cwynAFyGTNDT6pHpeBKdLxNRDwYx5IWKTNNzHlz9wXdbtS1qAE5VmvDM4k1b0c+1Y6AmiPNrtNQLcpbA9IRkObMY1IBklJ+pZoL0cw++jP9tEXb6E33cnDoIGkjTXOkefoXTkI8pD7PjBoZYs4eG3+QrF0TmhUpOCBnR2LVcMF5aQw3MpIbwbRM5SzUNSDNGY4epKSZECklg5lB1s9bD8BLfS8B0Bad+R42px/ciJlxm37iD2H297MAACAASURBVJKxI6BamBCcFJzh80RAASWOVY+Aws1Y0lJOQdAuOM0ZjxYgzYQk80nyVp71rUqAXu59GYB50XkzvqYbAZkZt+lnzVNwvvEpuJwdnVWjE4IXx2HoOuH0RlTNGU4pveBK67OiOaNw6j8L6xfSEmkpjwAFlQAlzKzaAwR2BKRScDWJgCZKwTk1oCqn4Jz6muuE0yYEzRlOKRHQr4UQDwoh3iEm6rOvOSPpzygnZ3OkmY5Yh2vFLkcElLCyhfHXvgBZWbsakBACP2B4/mnn/WpttUjBQUH83RRcKKYjIM2U/P7v//6i22+/fUGt13GqlCJAq4F7gY8A+4UQfy6EWF3ZZWlqjfMm2BJpYWH9QgD8wn96LjhHgMzcmAjIAmpjwwYIIIpqQE4KruoRkJ2Cc51wRhaED4J1WoA0p0w+P/uj5lI6IUgp5U+klDcAN6La8PxGCPGsEOKNFV+hpiY4dYjmSLPb7aC9vr3kvm8T4QqQzHkioCBZbAGqQQoOlADli2pAam1VrwGFxwpQGgIRCIR1Ck4zji9+8YvtnZ2dGy655JLV+/btCwNcdNFF5958880dF1544bl/+qd/uuCRRx6Jr127dt3q1avXXX/99Z3pdFoAdHR0bPzMZz7TsXHjxrUbN25cu2vXrpr855u2GakQohXV/+0jwEngc8CjwGbgQWB5JReoqQ3OJtTmcDMXtl/Id3d+l/PmnXda1wz5Q4T9YUYto2BC8AdVJ2pqk4IDOwLypOByvgBY1d+IWheoI+gLFkwIRlaJjz01VjNL+fffXULP7rKOY2D+uhTXfnPSJqc///nP6x5++OGWnTt37s7n82zevHnd+eefnwLVJ+6FF17Ym0qlxIoVKzY+/fTTe88777zse97zns6vfe1rbbfffnsPQENDg7lz585X77777tbPfe5zS5555pn9Zf0ZSqCUFNyvgAbgWinlO6WUD0kpDSnlVuCeyi5PUysGM4NEA1HqgnW8ceEb+eZbvsmX3/jl075uLBhjRObVIDpQKTifevOvXQqOMQKkorxq14CEEMXteIyMioD8IS1AmiKeeeaZ2Dve8Y6heDxutbS0WG9961vdHcw33HDDAMBLL70UWbx4cfa8887LAvzO7/xO/y9+8Qu3ucDHPvaxAYBPfvKTA9u3b4+NvUc1KGUcwx9KKX/kPSCEuF5K+aCU8s6pXiiEuAr4BuAHviulvGPM89cAfwJYgAF8Xkr5i1P5ATSVYSAz4BbFhRC8afGbynLdeCjOaKYXfFF1wKciIL/wVz3icBgbAeX9fjCqXwMCtRnV7YhtZG0BCuoU3GxmikilkkzmCYvH4xaovXxT4fPUPYUQNRnVUEoENFHDr/93uhcJIfzAN4G3A+uAG4QQ68ac9l/AJinlZuATwHdLWI+mCgxmBk+r68FkxENxEtIsNiEIQbjK0YaXgATDM+Q3L2rjggNVB3JrQPm0joA0E3LllVeOPvHEE02jo6NicHDQ95Of/GScO2jz5s2Zrq6ukFPf+f73v9962WWXJZznv//977cA/P3f/33z+eefn6ze6gtM1Q377cA7gA4hxF2epxpQ0cp0XATsl1IetK/3Q+AaYLdzgpRy1HN+PTDrBiadTbzQ/QJBX5DN8zczkBmgNVr+sSHxUJxRzEInBF+ArBBEamRAAAgyNgWnPpfVSoAODR9SD9wakBYgTTGXXnpp6j3vec/Ahg0b1nd0dGQvuuii0bHn1NXVyXvuuefw9ddfv9I0TTZt2pS69dZbe53ns9msOO+889ZYliV++MMfHqzuT6CYKgV3HNgKvBvY5jmeAP6fEq7dAXhD02PAb409SQjxHuAvUFNW31nCdTUVIG2k+cRTnwDgxY+8yGB20O18XU5iwRjHsYpMCBkhCNco/QYT1YDU97VICTaGGz0uuIxOwWkm5c477+y+8847i8YqfPWrXz3pfXzNNdckrrnmmt1MwK233tr79a9//cREz1WLSVNwUsqXpJT/AKyUUv6D5+shKeXgZK/zMFGCclyEI6V8WEq5BrgWVQ8afyEhPiWE2CqE2Nrb2zvRKZrTZM+AOyqE/YP76Uv3VSQCigaiat+P39OKRwjCNai3OASkLArpc6K2EdBwbljl740sBCdOwW3v2c4Hn/ggB4dq8sFVoykLkwqQEMIxHmwXQrzs+dophHi5hGsfA5Z4Hi9GRVUTIqX8GbBSCDFuq72U8l4p5QVSygva2mbeDFMzOccSx9zvf3bsZxiWQXtde9nvEwlEyAhZ2Adk27AjtawBAXmvCcH+vhYC1BhuxLAMUkZqShfcfTvvY2ffTh587cGqr1Ez9+nq6tq5cOHCUkopFWWqFNwt9p9Xz/DaLwCrhBDLgS7U8LoPek8QQpwDHJBSSiHEG4AQMG6an6bynEypyD0gAjx77FlAbTwtN9FAlCx4IqAAGZ8g7CvFkFkZxkdAdgquBr16vZtR642MqgEBGMUCtHtAZVV29e1yj/3q+K/oSfVwzTnXVGexGs1pMtU8ICc32AekpZSW3YJnDfDj6S4spTSEEDcDT6Fs2PdJKV8RQtxkP38PcB3wUSFEHkgD75fTeQc1FeFk8iTxUJyVjSvZ0bsDqIwAhf1h0gKk8KscrT+kIqBa1oCkJO9JGOeEICAlPjNfEMoq0RhuBJQAdTgRkLSKIqBUPkVPqgeAI4kjABiWwad+8ikALl54MQvq51xbMM1ZSCkfO38GXCaEaEbZprcC7wc+NN0LpZRPAk+OOXaP5/s7gSn3EmmqQ0+qhwV1C1jdvLqiAuSMXMj6/ETANiH4iItaRkCQ8jzOAUEpwcwC5d3gPh1uR+zMcGEfkGUUmRAOjxwG4Lx55/Fy38skcgmOjxay2zv7dmoB0swJSskxCCllCngv8HdSyveg9vVoziB6Uj3Mr5vPhnkb3GMtkfJP4ogG1AbU7Bgbdi1TcEEkhscfkxcQknJc2qsaFPWDK3LBFdZyePgwgLs5+MjIEQ4MHXCfd6IijWa2U5IA2U1HPwQ8YR+r3buFpiI4AnTFkitYEl/CzZtvrsh9nH5vaUeAAmGVgmPmTU5Pl4BlFdWA8tgCZGYne0nFcFJww7lhyGcm3Ad0aOQQAsGlHZcCSnAOjRzCJ3xE/JGiaEhzZtLX1+e/44472gAef/zx+BVXXHFONe9frnuWIkC3oDofPGzXcFYAz5zujTWzB8My6Mv0Mb9uPk2RJp54zxN8etOnK3IvJwXn9H8jGFX7gGo4aiogLQxvDQhJSFKTCMitAWWGxrjgPCm44cMsii1iRdMKQEVAh4YP0RHrYFnDMrqT3RNeW3Pm0N/f7//7v//7+afyGsOoueltHKWMY/iZlPLdTt83KeVBKeXvVX5pmmoxkBnAkhbzo+rfcyXnDjoRUNbpQxWsI+sTRGSNBciTgsshPTWgKq/FFyAWjDGcGQDkhN2wD48cprOxk2ggyvy6+RxJHOHg8EFWNK5gYWwhx5M6AjrT+cIXvrD46NGj4TVr1qy77bbbFieTSf9VV121Yvny5evf/e53L7csNeKko6Nj46233rpwy5Yt5953333N3/72t1tWr169btWqVes/85nPdDjXq6urO9/5/v7772++7rrrOgFeeeWV8KZNm9Zs2LBh7ec///lF3vMmu+epUMo4htXArUCn93wp5ZWnfDfNrKQv3QfAvLqZTzstFScCSjsiF6xTEVDF7zw5QcvCkF4BsuwaUPUFCJyGpAP24qJKfKw8SIkhTQ4OHeS32lVTkaXxpRwaPsTh4cNcuuhSMmaGbd3bpri6ptx8+bkvL9k/uL+sbpVzms9J/cn/+JNJm5x+/etfP3b11VdH9+zZs/vxxx+P33DDDSt37NhxsLOzM79ly5Y1P/nJT2Jve9vbRgEikYi1bdu2vYcPHw6+8Y1vXLNt27ZX29rajMsuu2z1Aw880PSRj3xkaLL73HzzzUs++9nP9nz6058e+Mu//MuiTZivvvpqdLJ7lkopKbgHge3AHwJ/4PnSnCE4AtQWrfwm34jd8cAeNo0ZCGMIQbiG7nu/NIv3AUmLELJm7W+awk0MOWO5nRoQgJnjSOIIOSvH6hY1lHhpw1J29u0kb+U5t+Vc5tfNJ5FPkMqnJrm65kxk48aNyZUrV+b9fj/r169PHThwwN1F/dGPfnQQ4Be/+EX9xRdfnFi0aJERDAZ5//vfP/Dss89OOYZh+/btsU984hMDADfeeGPRHs2p7lkqpZgJDCnl/3eqF9bMHdwIKFqFCMh2u2XsAMhxw4WtGgmQlPgtC9PrgpNWzUwIYLfjSdktpwIRFQUB5NO8NvgaAKuaVJ++zoZO93VrW9Zi2ePNe1I9dDYWntNUjqkilWoRDofdf8B+vx/DKFQ1SxnP4E27O1NTT+eepVJKBPSYEOKzQoiFQogW5+tUb6SZvfTab3bVEKCoUILjCJDTdSAsTz1/XBbMPH4oEqAcFsEamRDATsHl7KF0gYj6AjAy7Bvch1/4XQPCRe0Xua9b1rCMBXVq/4/T2eJo4ijbe7ZXb/GaqtDY2Ggmk8lTatXxpje9KfnrX/86fuLEiYBhGDz44IMtl19++ShAa2tr/sUXX4yYpskjjzzizmHZvHnz6Pe+971mgPvuu6/s7/ulREAfs//0pt0ksKLci9HUhr50Hw2hhqr0PgvZn3mytgBlLPUmH7Zq5NAxcwTkGAGyTGI1jIDqg/WkjLR6EIhA0E4F5lO8NvgayxqWEbbHV6yft57bLrqNVU2r8Pv87gZUp1PCJ5/+JF2jXfzkfT+pyMZiTW1ob283t2zZMrpq1ar14XDYamtrmzZfvGzZsvztt9/e9eY3v3m1lFK85S1vGf7whz88BPCVr3yl65prrjln4cKF+TVr1qQdcfu7v/u7ox/60IeW33XXXe1vfetbh2KxmFnOn2NaAZJSLi/nDTWzj750X1XqP6Da3kBh/EHOdneFzbL+uy4dK48fieGJwPKYNTUhRANRMo74BSLKgACQVxHQxnkbi87/0NpCU5L5dcrJeDJ1kv50P12jXQBsO7mNd67Q007OJB577LFDEx3//ve/7+5E7urq2ul97qabbhq46aabBsa+5uMf//jgxz/+8XFTDjo7O/M7duzY4/P5uPfee5s3btyYBLj66qsTV199tXe43Yx2P5figqsDfh9YKqX8lBBiFXCulPLxmdxQM/voS/dVJf0GyvIMhQmkGVPZEWoXAeXxj4mA8pZh14Bqk4KLBCKkzRwSEIGwasUDjGYG6Brt4r2r3jvpa6OBKPFQnJPJk+wd3OsePzKiuyNoTp3nnnuu7pZbblkqpaShocH83ve+d7ic1y8lBXc/aiDdJfbjYyhnnBagM4STqZNsWbClKvcKThoB1WjgmpkjgMTy2rAdAaphBCSRqkOE04wUODqiat3LG6dOSiyoW0BPqoe9A0qAYsGY2z9OozkVrrrqqtG9e/dOONCuHJRSxFoppfxLVIcSpJRpJh42p5mDZIwMJ5InWNawrCr3c1Jwjtxk7VRT2KiVAKkIyMByXUI5aSgTQo0iIKdfXloINZDOdsE5dR3HaDAZC+oWcDKlIqAFdQvY1LaJ10der+yizz4sy7L0++A02L+jSR1GpQhQTggRxZ5mKoRYCdTmo6Gm7DhvTF47byUJWKrW4xg2s3aUEa7Rm71ywSnhcSzMOStf0wioLqD2NKZ9osgF15NRdnmnzjMZ8+vm053s5tX+V1nTsoaFsYWcSNZ08vKZyK7e3t5GLUKTY1mW6O3tbQR2TXZOKSm4Pwb+A1gihPgn4H8AHy/LCjU1x+miXK09IwEmjoBC+VoJkHLBAZjSxI+fvJm3N6LWRoDcfnlC2APp1Htcb2YAgZh2VPrKppU8vP9h+jP9XNV5FX6fn4HMABkj415bc3oYhnFjd3f3d7u7uzdQ2gf5sxEL2GUYxo2TnVCKC+5pIcQ24GLU/4RbpJR95Vujppb8+NCPmRedxzlN1Wmm67dMfFIWIiD7TT5So2jDccGBasoa8ocKEVCN6lKFFJxPRT9Cvb/1ZAdpibQQnGZ43/rW9e7361rXqc7aqFpftVKtZzpbtmzpAd5d63XMdaZVbiHEf0kp+6WUT0gpH5dS9gkh/qsai9NUlpPJk/ys62dce861BKo1j8eOOBzPWyECykz+moquR9WAQEVApmViSpMgoqYmBICUz46A7MeD+STNkeapXgoo0XHYsmAL7XVq/4/ukq2ZbUz6riOEiKDGQc6zp6E6uc4GYFEV1qYpN7kUpPqhaQkA/3nkP7GkxbXnXFu9NZh5AkjydsHfjYDy6eqtoWg9OTcCMi2TnL0xNij8NbVhg21CCETBp+Q6Y2api0ydfgOoC9bxz+/4Z/w+P7FQjIX1CwEtQJrZx1Qfez8NfB4lNtsoCNAI8M0Kr0tTCf79Jtj9CHzxMESbebn3ZRbULahuWsbME5ASwzbGuBFQrlYClHdrQIY08FkqKRAS/ppHQBnhU6MYnP55Vq7kGs7GtsJmVac7gjYiaGYbk6bgpJTfsLsg3CqlXCGlXG5/bZJS3l3FNWrKwdHfKPEBOPkKoDagLopVOZg1cwQYn4KL5GrUvdnjgjMt092XFBKBmpkQ3BpQIARCgM8H/jBpKz8jE0HIH6Il0uJGQA/ve5j//dT/ZiQ3UtZ1azSnSikmhL8TQlzC+HlA36/gujTl5sB/F77PqKL0QGaApfGl1V2HlScopTsAzo2ArLwq+vunLrCXHTNXVAOyTBWZhXyBmjUjdYf2BTy/i2CUjJV3nztV2uvb6U4pAbrjN3eQMlI8ffhp3rf6fae9Xo1mppTSiucBYCWwA3AadklAC9BcYsjTMT6t5k8NZAbY1LapuuuwU15O77WsmSUk/Cq/m09XX4CsvGsNNy0Ty04NBn21i4Acl1ve2xw2GCVtGW50dKosrF/I4eHDdCe7SRkq2ny1/9XTXqtGczqUYn26AFgnpxomoZn9ZIchvhASJyAzhCUthrPDtESqPFnDzBGUkrxTAzKyhIX9zzCfhkhDlddTcMEZ0sCw+66FfMGaRUBOV/K83/PfMxAhI80Z7+Npr2/n+RPPu6MZwv4wB4YPnPZaNZrToZQNVLuAGfVxF0JcJYTYK4TYL4S4bYLnPySEeNn++qUQosofx88iMiPQYNd7sglGsiOY0qy+ABlZAkh3BHbWzBJ2LOC1mOJp5vDL8S640CyIgHJea3ywjgzmjCOg9rp2kvkkzx57lrpAHVd1XsXBoYPlWK5GM2NKEaB5wG4hxFNCiEedr+leJITwo9xybwfWATcIIdaNOe0Q8GYp5XnAnwD3ntryNSWTHYFosxrvnE8zkFEd2UvZV1JW8ml7H5Ddd83MEbbHdFMLK7ZtigBVA8rbm09VBFQbAXL2ZOX9fveYFQyTQc44AnIG2D1x8Ak2z9/MOU3nMJgd1EYETU0ptRXPTLgI2C+lPAgghPghcA3gdlaVUv7Sc/7zwOIZ3kszHdlRaO5UjS1rLUAUUnAZM0PYX0MBymeKUnCOCy7oD0KNGqQKIQgCeU8ElLWFZ6YmBG+t7w3z30BHvAOArkQXDa1VTntqNDaluOCeneG1OwDvrPRjwG9Ncf7/Bn480RNCiE8BnwJYurTKrq0zBSOjNjUG6yCfYjCrZk+1lrCxsazkUwSlwJDKz5Izcx4BSlZ3LfZ6vDbsvOVEQGEwa7Aem6AU5HyFBEUmEAaLGUdAjeFGLmy/kB09O7hq+VWk7HTnsdFjrG1dW5Y1azSnylSdEBLARMYDAUgp5XQfmybqEjuhkUEIcQVKgC6d6Hkp5b3Y6bkLLrhAmyFmgpFRbV2CdZBPM5hRAlSTCEj43FSXioDC7nNVJ5/Gb2eiTWl6IqAQGOMGRFaNIJAXhRRcJhiCLDOuAQHcfeXdDGWHWBRbRCKnhll2JbpOd6kazYyZVICklPHTvPYxYInn8WLg+NiThBDnAd8F3i6l7D/Ne57VPHbgMXb37+bWC27F7/MXP5nPqPSbLUD9GfWrbg5XW4BSBISPrFRuMxUBRdznqk4+TcAWQMMyCiYEf7hmJgSAkIS8JwJK21HiTFNwoFr01AXVqId4KE5DqIFjo8dOb6EazWlQyQ6ULwCrhBDLgS7gA8AHvScIIZYCDwEfkVK+VsG1nPFkjAxf+sWXALhs8WVcsuiS4hPcCCiqUnCZQeKhuKp1VJN8ioDwu3bnvJmnrqY1oBT+gLq/NwIK+cM1MyEABJHkRSGJkLEt2eUcp7A4vlgLkKamVGyOhZTSAG4GngJeBX4kpXxFCHGTEOIm+7TbgVbgW0KIHUKIrZVaz5mOM9cH4OfHfl78pGWClbdrQFE3BVd1CzZAaoCAP+gKkCENgm4KrhYRUAq/ff+iGlCgxgIkJXnhqQHZHxROJwU3lsWxxUUpuHytxqJrzloqOkhJSvmklHK1lHKllPLP7GP3SCnvsb+/UUrZLKXcbH9dUMn1nMk4mwo7Yh08f+L54icNe9SBWwNSEVDV028AyT4C/nBBgCyDQKCGEVAuScAeeV1UAwpEapqCC1pWcQRkp1TLKUAd8Q66RruwpMVdL97FJT+4hD0De8p2fY1mOvQkvzOE/UP7CfqCXHvOtewf2u+aDABV/wG7BqRScP2Z/uobEABSfQQDETfSMCzDjUBqIkDpIfyhuLuWQgouUrNOCKAioJzHxpO29wRFphlGdyosji0mb+U5ljjGfbvuI2NmeHjfw2W7vkYzHVqAzhAODB2gs7GTixdeDMD9r9zPH/3yjxjODk8QAdUoBSclJHsJBOuKIyB/0N4gW4MUXGYYf1gZOk3p6YQQjNYuApKSoGXiTYilbUNCRE5kLp0Zi+Nq292/vvavmLYtflf/rrJdX6OZDi1Ac5Q9A3v46I8/ysu9LwNKgM5pPIf1reuJBWPcv+t+Htr3EA++9qBHgFQEZOVTDGWHqi9AmSGwDIKBKKal3vAMaRAQAbc2VXUywwTCKgIyrUInhGAgCtIC05jq1ZXBMghJSd6jNRk7HReZcHfDzFjTsgZQH1bqg/Vcv/p69g3uw7IbxWo0lUYL0BzlH175B7b3bOf+XfeTyqfoGu1iRdMKgv4gn3/D51nRqFqv7O7fPSYCijJiZjClWf0UXFJZv/3BOgzpiYB8ARWZ5aq88VNKSPXhjzSptUiDvJXHL/z4HbdZLaKgfJogkpxn25wjQNEyikNLpIW1LWoT6psXv5mN8zaSNtIcTRyd5pUaTXnQAjRH2dWnUiXPHX+OV/rVgLlVTasAeP+a9/PItY/wpsVv4vDI4TE1oDoGpPqUX/UIKNkLQCBUX5yC89UoAsoMKxdcXE0MNS2zsJ6AXZeqhRPOyBKUkC8SIPVnxCrvPuybNt3EprZN3LTpJrdf3KHhQ2W9h0YzGZXcB6SpEHkzz9HEUda1rmN3/27++dV/BmBN65qi8zobOvn1iV9j5VPqk4YdAQ3ajqrqR0BKgPyhWJEABX1BtzZVVUaUBTkQU83eTals2AFfQNWkAMwaGBGMjErBeSagpJEIKQlZ5U0JXrn0Sq5ceiWAqhcCh4cPF28h12gqhI6A5iCvj7yOKU2uX309QV+Q/zzyn8SDcRbVF4/XXtawjKyZpSfVow4EIhCsY9Cv/tqrHgGl+tQyQjG36F0cAVXZhNCrLMf+1nPcteStvBLEmkdAhYatgOqELSXCSadWgMZwIy2RFhU1azRVQAvQHOTgsJrjsq51HVsWbAFgy4ItCFFcoG6vV5/sTxYJUJQB29Jb9X1ASUeA4rMjBXdyNwg//nnnAioCctfjWMNrEgGllQBJrwBZRKWs+O+os6FTp+A0VUOn4OYgjgB1NnRy8/k3kzWzfHrTp8edt6BO1TZOZlTqqyBANYqAkr0QaSQQCGNKEymlcsE5JoR0lZt/9u6B1pX4nY2oRTUgOwVXyxqQV4CkSdSSBUNJhVjWsIxnj820Ab5Gc2roCGgOcnD4IIvqF1EXrGNT2ya+//bvs2HehnHnORFQd9ru8WqbEAZ9fuKBuur3gUv2Qd08/HaX56ztMFM27LqCWaJaDByClpXuADhTmgVbuBsB1UKAMsoFZ6cpAdLSJCKtiqcpOxs7GcgM6EF1mqqgBWgOMJgZ5GM//hj377ofUC6l5U3Lp31dQ6iBiD/CSXv4HKF6ZULw+2gOxiq55IlJ9kJ9m/uGn7E/zbsRUDVTcFLC4GFo7nQF0bCMCSKgWpoQCoaDjGUQkbLiIr28Qf272j+4v6L30WhAC9Cc4NEDj/Jiz4v89ba/JpVPcXj4sLvPZyqEELTXt3MyN6QOBOuUDdvvpzlQV+FVT0CyD+rnFQTI9ApQtLoD6VL96n7Ny9zRFc5I7qA/WOMIyDYhWJ4IyMopC7ZRWZE+r+08ALb3bMeSlt6UqqkoWoDmAFtPFpqEP7D7ATJmhnObzy3ptQvqFtCdTwDC7QU34vPR4AtXaLVTkCoWIDcFVwsTwtAR9WfTUjcCMi1PCq6WLrh82q4BGUjbip2x8nYEVNnfUWu0lc6GTu59+V6u+NEVfPLpT2oR0lQMLUBzgNcGXuOty95KU7iJu3fcDcDm+ZtLeu2C+gWcNJIq/SYEBOtI+AQNZWxqWRJSKpNBtMV9w3dTcE4NyMiAVaU3u7Sdlqybh0/48AkfhjQK+5JcAapyXcq+Z8gWHsctmDFztguu8ut518p3kTJSDGQG+E33b3il75WK31NzdqIFaJaTyCU4njzO2ta1XLHkCkC515bGl5b0+gV1C+i1spj2JEyCURI+H3FRZQNkPgWWoVxwE6XgQvb6KpxicknbaclIIwB+4S92wTm/rxqNiAjaAuQ0R02bGSKSqvx+PrHhE3zzLd/kifc8gV/4+emxn1b8npqzE23DnuXsH1LF4FVNq3hb59s4MHyAD5z7gXF7fiajvb4dE0lfOMoCwAqElQDhn/a1ZSVju6oiDQUXnOFNwXnegoEEJAAAIABJREFU8EP1VViP2vVPtMldQ9E+IGc91e5PB5AbJWi34cmbeQiqaDEqfFURxIAvwJsWvwmAc1vOdRveFnHo59B2LsTmV3w9mjMXHQHNcl4bUJPKVzevZkl8Cf/0jn/iXSvfVfLr3b1AYbXXJSXAEoLGMnZVLgnnDT/cMLkJAarXDSHjREBKgPz2mHC3E4ITkdViREQuSdD+r+nMTcqYGSIiUPWIzGn3JD1tgeh9Df7hanj6D6u6Fs2ZhxagWc6+oX3Eg3F3T8+p4u4Fsm3FI3bhP17mppbTknUioKZJbNi2AOWq9IafHrI35qqu136ff0wEVF/d9XjJjhK0e9HlrBxSSjKGLUBVrkmta13HSG6EY6PHCgcH7PHvr/+yqmvRnHloAZplGJbBh574EO9//P1kjAz7BvdxTvM5JafcxuJGQAGV9krkRwFoyFW783QhBTdxBFTliCMz5NZ/oFADKjQjtRuSVtMa7pBLEvQrYcybeQzLwJQmUV+g6hHZ+tb1AG7HdcBtKouZn+AVGk3paAGaZewd2MvLfS+zu383jx54lH2D+9wxCzOhMdxIWMJJvxIwZ4d7LFPlne6ZQtF/rAvO7YYN1UsxZYbd9BsoJ15RBAT2jKJapOBGCdkuvJyVI22q30nEF6x6t4hVTasI+oJqrpSD3dOPGX4o0mgctAlhlrGjdwcATeEm/uT5PwHgwvYLZ3w9IQRtFvQKez+J/aYfSVW575qTggs3EMieLFqLsmHbtvBqCVB6yDUggErBOZ0Qgo5FPVRfsxpQKBAB8uTNfOHvzBeqnkvQJugPsqp51cQCVIvfjeaMQkdAs4ztPdtpr2/n97f8vnvsko5LTuuacctkVKi/amfzZyTVf1rXPGW8KTgxyUZUqFkKzid8xfOAoDZTWsEWIPX7yFm5ggD5QzWxhY8zIjgpuMywTsNpTouKCpAQ4iohxF4hxH4hxG0TPL9GCPErIURWCHFrJdcyF5BSsr1nO+e3nc/VK6/ms5s/y99e/rc0hBpmflHTIG7kGbX/pp03/XBqsLpvHplhsGs9zht82v40XyxANUrB+QLF+4BAOeFqEgGNEnQEyMy5v6doIFL9hq2oOlAil+BYwjYiOAIE1e9grjmjqJgACSH8wDeBtwPrgBuEEOvGnDYA/B7wV5Vax1yiO9lNT6qHTfM3EfQF+cymz/CWZW85vYtmhohZFgl7uJkbAVkWDB893SWXTnYEwg0ghNt7zVmL3+f31ICqFHGMTcEJv+qE4LTiAeWEq0kNKEnY/n1kzawrQBF/lfvl2axrVf9td/WrMfAkugtPVjuS1pxRVDICugjYL6U8KKXMAT8ErvGeIKXskVK+AOg4HpV+Azh//vnlu+hojxIgZz+Jnc4JS6nGEVSLzLCb8pqwF1yoiiYEyypaD9g27AkjoBqk4DLDhEKqW3nezLtuwWgoBtnRqi9nVfMqIv4IL/W+pA6MHIc2e/y7FiDNaVBJAeoAvB+xj9nHThkhxKeEEFuFEFt7e3unf8EcImtm2dq9FcMyeLHnRaKBKKubV5fvBiNdxC3JqC1AbgQkJQxWU4BGIKJSiU6EMa4XHFQn5ZUdAeSkLjh3TlKovvpv+JYJqX5CdfOAMTWgUAyyiequB+VS3DBvg/qAlB2F7DC0b1RPagHSnAaVFKCJPJoz2v0opbxXSnmBlPKCtra201zW7OLLv/gyH3/q4/zp83/KL7p+wUXtFxU+gZeD4aPELIukmcGSlvtpOuwPq3k41cJJwcG4fUBBXxD8QVUjqkYENKYNDxRScHkrX0jB1bVW/w02NQBIQvWqxU3WzLoCFA03KBecaUxxgcpw/vzz2Tuwl9SgmsarBUhTDiopQMeAJZ7Hi4HjFbzfnGM4O8xTrz8FwL/t+ze6RrvchqPlu0kXcQkSSTKfJGtkCfgC+JuW1SwFN2EvOKjevpsxbXjAtmGbBpa0Cjbs+vmqa3Y1zRojqtAfii8Eik0IkbCdMsxVPw23cd5GTGnyWvc2dWCBPYFXC5DmNKikAL0ArBJCLBdChIAPAI9W8H5zjl+f+DWWtLjvbffxlqVv4bKOy3jnineW9yYjXcRDcUB11s6aWRX9NHfC0OvlvddUZEbG1YCcDZauAIXjhf1ClWRMJ2xQolhUkwKI2dF2Nd9k+1Tz2VCr2nyctwo1oFoKUGdjJwBH7N6ENHdCKA5JLUCamVOxjahSSkMIcTPwFOAH7pNSviKEuMl+/h4hRDuwFWgALCHE54F1UsqzYiD98yeepz5Yz+b5m09rs+mUjHRRF2kGhknlUwUBii6AEzsqc8+JyAy7KTjXBTc2Aoq2VMfWO0EKLuALkMglitdTbwvQaA/EZ9aL75R57ccQihFqVXXAnFkYCR5xIrYa1IEWxxbjF35eT9iD/BoWqd9J4oR7jpSSnX07WdG4glioBiPfNXOOinZCkFI+CTw55tg9nu+7Uam5s5LnTzzPhQsuLKR8KkGim2jLQsgPkzEzZM0sEX9EpZeSfcoR5qvwfmTLhFzCjTicn7eoFxxAXbNdA6kwE6XghN9NdRWl4KB430sl2fMk7Po3uOwLhMLqDTxrZt0NoJFoizqvBk64oD/IotgijqROQrRZ7dtqWKQccTZf2/o1Htj9ABvnbeSBtz/gftDQaCZDd0KoEUdGjnA0cZSLF11c2RuNnCBif5JPG2kyRoZwIKzmuEizMBm0koyJOMZNRC2KgKqwnolScL4JUnBO1DN8jKrw8r9ArB0u/xIBXwCf8KkakJkm7A/jsyNIctWPgEA1tu3JJ6DBNrM2dLgC9HLvyzyw+wE6Yh3s7NvJk4fU586hzBAHhg7UZL2a2Y8WoBrxzNFnANzBXxUhm4Bcgki96oidMbwRkJNeOlm5+zs4abVoMzDeBee6zqJVjICEX9WcbAIi4AqQGwE1LVNpw8d+D374IRXJVQop4eAzsOq3VSduIOQLqRqQkSESiIAdFdUiBQfQGm1lwMyoyAegYSEkTmAaOf7s13/G/Oh8/vVd/8rKxpV875Xv0ZPq4bpHr+PaR67lJ6//pCZr1sxutADVgO5kNz/Y8wM2ztvIkviS6V8wU0ZUfj4aU46qtJEmY2ZUDciZZDnaU7n7O4wd/jZZDajOrgHJCs8qSvUri7Wnm7Pf5y9uDQQqNfnG31Xf73kcjjxfuTWlB1WkOH+9eyjkD7mdEKKBqIoQoToiPQGtkVb6McF26NGwCKTJv+/+R3b37+YLF3yBWCjGx9Z/jNcGX+NdD7+LRD5BQ6iBu7ffXTzUTqNBC1DVOTpylPc++l4GMgN8/g2fr+zN7AJx1P7EmjEzZI2snYJTUVFV6htuBGSPvxZjIiBvCk6ahZRdpUj2KQHyMKELDuDy2+CLh9X3x35TuTU5jsSmpe6hkD9EzlQbUSP+iP2hQRS3wqkireFGEj5B1klN2qm4f9n/EGtb1vL2/7+9Mw+PqsgW+K96S2ffQ0JCSAg7grIqCKKyiIw7OOrwdNRR4SmM6+io89QZZ1xwQ8VxQHFEZ1BRERFQQFEWQVYh7IQtQDYIIXs6vdX7o26STkiAQJoEqN/35evuutX3VnV36tyz1DmpVwNwTdo1DE0eis1s49XBr/J438fZU7SH9YfWN8u4NS0XLYDOMO9sfAev9DLzmpn0S+jn34sZAsgeprSs5jPB1daAfCuimoQJk5Gpu9rn4hNZ5RfKCyA4plaT2WTGK1W+vGOCQgIjlcA+sst/Yyo0ost8BZCPCS7QEqg26wbH+v/zaYBoI2apINDwRYW15pDZzLaSTEamjqwummg1WXnjijdYestSLku6jGFthxFsDWZWxqxmGbem5aIFkD/Z9QO8cwkUqN3jbq+bZQeXMSR5SPW+Cr9iOIjt4WpRq2WCs4erip9n0gQXWFsASWSN/wcgMlU9bv4SFjztv1II5fnK3OeD7zjqzUQR3R6O+NGZXmhkraqjAVV6KqnwVCgfEBihz82kAXmUgD4SYKRNCksky6i02z6y/TH9qwRSkDWIESkjWJS5iFJnKXN2z+HBxQ/q4ASNLkjnV9b9Gw5vg1VT4OqXST+cTrGz2L+BB76U5EBAGPZgZW5yuB04PU61mAmh7urPhACqMsEZGpBJmBAIJYB8F/vYTupx6SvqMTQeBkxo+vFU+YB88A0Zrl8ApakwaX9RuF8FPPjsTaoSQA63gyCLseiHJjSbBhRllIIosKpqrQRFk2tTgjE+6Pj7pG7qcBNfZnzJ2EVjSc9PByC3PJdPfvNJjQasOe/Q37w/KTX8KxtmQGUpy7KWYRZm+rfuf2auX5IDoQlYTVYsJkt1GLbNbFPH62wk9Bulh9XiarVXN1Ut8rUW+4AQ6HIdGLVw2LOk6cfi9SiBWI8PqIp692VFpirNyV97cAr319J+AOwWOw63oyYKDiCyrdKo/RmR1wCRlUojLTAZwRtCkBusIhvjg48vgLrHdGdEygjS89MZ3nY4z/Z/lq1HtrI8a7lfx6xp2WgB5E8K9kBMJ5VeZvOXLD24lJ5xPU+vwFxjKM5RobJAoDmw9kZUUALoTPiASrJrIqcM6hVAAKOmwSNbocctcGhb04+lOBukt2YvS53x1DsmUAs/+C99UT0CKNASWG02rRZASX1VKp5d3/tnHMchukyZUo96aori5dmDCZbihJkPhBBMvGwi34/+ntcuf43r064nLjCOj7d+7Ncxa1o2WgD5i8pSKDsEF94CQTHkHljOzqM7z5z5DQwNSEXAVd9NV/mAQG16PBMaUHF2zd4RgyqfSy0fEIDFpvwzcV1UYs6qAIamoqoERVRqrebqz4QGNKCIFOP9+5p2PKDCzo8jgKrDsAE6jYTwZJjxW9j0RdOP5TgEluQSIKGgsiZdUq7FSrz35MKrhRC0MvakWc1WbutyG7/k/MLOozspdBQyf898jlTo3HLnE1oAnSYOt4Npm6Yxc8dM3F6fNPnVUU1tITqNXwozABjQesCZGZjXo5zVRmRZoCWQMlcZbq9bhWGDOuYo8n8JhOKcYwRQlc+lwdITsV3U4+Edp399rxfmPgI/vlAdEFId8GBwQgEUmaIej/pBA6o4qrIb1COAqkxw1QIoIATGLVOa9S/vNv1YjoMoySVKmClw1OxDyjVJ4p2Vp7R36+aONxNoCeSpZU8xas4onlj2BGPmj6G8Ocqga5oFLYBOk+dWPsek9ZN4/pfnuePbO5j862RVuKtKAEWmQFQ7VjkPE2WPokNkhzMzsLJ8tacmrEYDKqxU2kSNCc4wi/kzqsrjVma+kzXBVRFrFOXL33n6Y9i3FNZOgyUvw7rpYLJCeO0UhNUmrobGFBQFthD/aEBFx0bAAQRZgqr9dtXfGahAhW43QPb6M5O8tYribCJN9loCKM/rpJXbdUoZw8MDwqv3CIUFhPFo70fJKs1i+pbpTTlqTQtGC6BGkl2azTe7v6HQUcjK7JXM2zOPsT3G8srgV8gqzWJK+hTGLRpH0RFj4YxIhqh2rDZ76RfX68xF/BRnqUdj4bdb7BRVqg2e1Xf7ocZmVH8IoLIjsO5DyFqrBGFM7SqvVU7/BgVQRFsVJn4k4/THslulPcIeoRbtpL5QJ1GmrwZU75iE8F8Ji6rw7iotyyDQEkipqxSn11mjtVaRMkj5sg6safrx1IfXAyU5RNlCqwWQy+PiiKeCeLe75vfWSEZ3HM2qMauYdd0s7rzgToYmD+XDLR/WEnKacxcdht0I8ivyuWXuLRRWFhIREIFE0ia0Dff2uJcAcwBDk4eyo2AHt867le9yV3KLJRCCYykPT+KQxUIn+xms5lptAlSbUAPNgRwuV1F5NXtKqjQgP9QJnPcwbP265nVir1qHT6gBmcwQlQb5TSCADqyGxD4wYDx8/xwMevSYLr4aRoPZySPa1pjwmpKqDa5RabWaq3xAQE0YdhWtLwIEZP8KHYc3/ZjqUpwN0kNUQCS7HUrryivPQyKJd3vU8YQLT+nUvp/3hJ4TWHxgMRPXTOSG9jewaN8iusV048b2N1bvK9KcO2gB1AheWPUC5a5yXhj4Al9mfEmJs4R/DPxH9d2zxWSha3RX2oW3Y0HBdm5J7AVCkB2oIoQSPWcwF5avDwq1mFWZ4Krv9iNTQZiaxs/ii7sSdi4w9hnlQeteaiOnDycUQAAxHeDQ1tMbi8eltJ4+d0O3G9VfPVSHph9vTJEpsHux0gaastRA/k4IbwO22kLG1yxY7QOqIiBUaZXZZyi9Tb76jUSHJnEkezdSSvLKVQRlK4/nlDWgurSLaMc93e9havpU5u2Zpxp3qmKKv+/2+ya5hqbloAXQcVh6cClvrX+L1PBU2kW0Y1HmIib0nMC1addybdq19b5HCMGV8f35d+FuipL7Ew5kW9TH3LryDDpXCzNVtgNjY6PdYq8p7Vx1t28LgugOkJOuFuqKozVJSk+HrHXgdsA1HygBGJFcK/En1Cz4xxVAsZ1g+zyVEcFRpDb09v49RLU7+bHkbVZjSWq44F9RhQsTJyGAEnqAuwIOb1ch4uunw9WvQFznkx9PFUtfUfujhj8PuZtrNuH64Ct0jhFAoLTK3YtVAIC/tQPjJiUuuiPOg4sorCwkt0yZbuM9slZdoNNl/EXj6R3Xm1JXKQMTB/L08qeZtG4SF8VdxIWxp6ZlaVom2gfUABlHM3hsyWMUOYtYnrWcf274J73ienFXt7tO+N7LCcQjBD+HqsU/y6V8L4m52/065loU7K3l1Pa9m67lT0i4UC1i/x4Jr3eBfU2wMXDfz4CA5P4QfwHYj933VCUELcKC2+OtP1NyUj/lP8paB988BD9PgrkPH//aLgf89DLs+Um9PrjWOFf9Aqiw3Mmw15cwaeHe6rZnZm+n598W8uv+Og7+qnNs+Qq+eRD2LlUmvcaStQ4W/x1WT4E5E1S2jJRBx3Q7oQBq3VNpmE24+DdI7iYIiiHBCKLJKcupFkCtAuOaNDpQCMGAxAEMTxlOkDWIv136N1oFt+KxJY+x5MASFu9frH1E5whaAzIod5Xz8daPKXeX07tVb15Z8wrB1mBmjJyBzWxjT9EeLoi54KSql3bPyyDK4+WnskxGogIXAjARvXsJHN5ZE+HlL7xeZZrpfE11k6+Pw9fhzgWjYNPMmkzPKyZDysDTu37mcmjV7Zh8awBHy5z8a8luXG5lwiqu8NL12QW0iwnms7H9CQ/0+Xzb9FMmwvmPq0XaHqEES36GMs9VHIU176u9Ma2MMgZr3oOfXgBzAExYCwdWKVNg+LGFd90eL49/kc6hkkrMXhdBxl7K7zYfBq+d1xft5OM/XFzzhqh2yuy19BV1/i7XKQ2t9DCENMK/lz5TBVh0uU4VoQPoMOyYbicWQIZfLXs9hCcee7wpOLQdtn+j5tnxKhKM0h45ZTnklecRag0lOLG98rP5iVBbKK8Nfo17Ft7D+MXjAQizhTHpikn+K2WvOSNoDQglfMZ9P47JGyYzfct0HvjhAfLK85h42URig2IJDwinZ1zPkyudLSWmPUsYbI1mefYKXB4XWaVZJIQkIuxh8MVd6i7dn+TvVItzm5rFM9Bas4DVCunteBXc9incuxgufQgyFqoItsbg9dbcAbsqYP+qeoXYoWIHv52ykilL97DvsBOA/UcchAda2XWolFcX1PFFBUYon83hbWrfy9ilSiBVLdrfPqE0if/erK7r9cDqqWqjJsCiZ5UvKm3IMSYqKSWPfb6RhVvzGNUrCbw1363NZOG+y9qxLCOfAwU+ZlMhYMSLKljgurdg8BNKQ9v+jTp+ZDd8/UBtLdJVofYfbZmtXnvcsHmW+tyvext63ApDnq0RoL7TP5EAir8ATBYViOAPSvLg3yPUZywEXPoQCcGGACrNIacsR20sTRkERfuVKddPdIvpxnejvmP6iOlMGz6N2MBY7v/+flblrKrdseKo/5LYapqc81IASSn5KuMrnlz2JJPWTeK+Rfex8fBGXh38KstuXcbkKycz54Y5p3Z3VbAHivYzKKE/Ja4SthVsI7s0m8SwZLjhX8onseb9pp8UwIZP4B+t4YPhquJn2hXVhwLNNQtYsDW45j1CQKerIbG30oZ8F1RnmTK9nGiT4dyH4M0e8OOLSkNxV0CH2pFZlW4Pt09bTVZhBX3aRlLqUAKhvBLuujSF0b2T+GztAQ6V1BHOV78CQ5+DMZ+rdDjtLoeNn6lxpc9UZrHiLFj9HuyYr4Ivrvo7XHwfbJml0iBd9LtjhrxgSx6zN2Tz8NCOPH9DN6SsMQakRIdx54AUhIDP1x6o/cb2Q+GP6+HCW5XQiO6gBIrLAZ/+Dn79D3z2PzVF45ZMVPuPPv+9Mk3u/kFlyOhxi/LB3TQFBj1S6xKF5U7u/nANCzbVZIHwuK2syzxa21RpDVQZIzJXqu8oY1HThmWv/UBlohi7DB7LgPgLiAiIwG62k1OWw8GSg6qg4gWjVC2nGbf4J3+fQXhAOL1a9aJfQj8+GPEBSaFJjP9hPN/u/VZ1cBSp7PPvXtpsVWM1jeO8E0D5Ffk88tMjPLPiGVZkr2D6lunsK97HS4Ne4qqUqwi1hTK4zeATJldskN2LAejUUQUp7C7cTXZpNq1DWqtw2eT+ykzk9TbVlBSVpbDwL8oU1Ko7/Oa1WmYnXx9QkDWovjNAfHdlZkqfqRbQ94fCvwYq30sVGz+Fz26vCY/OXKGc8QBLXoKvxqoUP3U0oA9/3seOvBLevq0nt/dvizQ0DimttI8NYezgNNweLx8s31d7TMHRMPDhmlxsF41Rd9vThqsgizGfQ9qVsOxVVcIhMhU6/QaufAYGPaYEWGpt/0qF08Pzc7fSOT6UB65II8hmIS6kJpdZanQorSMCuaxDLP9ZtZ9DxQ1orEIoQbRvGfxnlApOGP53tRAueVndjKycrEyhkakwexwsfh6C46DjiHpPWVbp5o+fbmDx9kN8vqamWOCfv9zOqHdXMGVpnTDwbjfC/hUw9XL472h189EUfjy3U2Vz7zBMBV8YfkMhBAkhCbUFUFAU3P4V2ILVGA43webhExBlj+L94e/TJboLjy99XOWUWzIRSnNVuqUfnvf7GDSnz3nhAyp2FjN141QWZi4ktywXs8nMI70f4c5ud+KWbszC3HQbRDMWQkQyiYmXYDPZ2Jy/maOVR5UAAuh3L3xxt0om2ZT7N1ZPVdmab/tE+U4Mih0uFm3JY3+pq7rtzvfTCbbt4Y7+bbn+Ih/fgRDQ9x5Y8BS8008tpLGdlQmm/VBlkpn9v8YGyNVw5zyV4iY8WaWHmX2/usO/9s3qBQsgr9jB24t3MbRLHEO6tGJLdlG1AMJrIzUmmNSYYK7unsDHK/cxIC2arq3DiAmps/kSoNtNqjz2rsVw4xRVLG74P5SpqOgg3DoDzMbPesj/1ftRvfPjLrIKK5g5tj8Ws/rekyPDqUp92i5OCaMnR3bmhnd+ZshrSxjVO4l7BqWSFFlHeF88VmU7z1wOfe9V5SOO7FYa2davwWKHka8ozWz6derx+n+q4nI+eL2SyT/u4v1leyh2uLn/8jTe/Tm/+vieQy7C7BZeW7iDyzvF0jneCOzoey/sXKg060GPweYv1Pdw/0olEDwupcn6lHmoFylhz49KqCf2Vn7B0jzoN/aYrgnBCaQfTsfhcdSUlG99Edw1H97uAwufVjcGfiY6MJppV03j0Z8e5Y21rzHoQBYpPf8HrMGweiqebjeyOSiI9XnrySrNIj44nivaXEFaRNqJT645I5zzAmj2rtm8vvZ1CisLuTL5SkZ1GMWwlGG0C1ehvFZxEn6dk6U4WxWhGzAes9lCanhqdbr5xBBjoe98rXKKr3hL3VnuXQYx7VVEU2M4uk8tcG0uVtFuyydB+2G1hM+K3flMmPErR8qcWMJzCDRk4N7DLlJj3Dz46QbyS538YaBPXrS+98DBNcqvcP07ysT1Tj/4+EalZcV1g2snwcc3weTe6j1jvlAL3G0zlGZnUou61ysprHDx1KxNOD1e/u+argCkxYYgpBIu0htAmyi1qP9peCfW7C3gjg+UQ3tk93hevKlH7cAEkwlu/rD2Z9GqK4xfB67yGk3JByklK3YfIauwgvzSSt5dspubeiXSL7UmSKJtVDjbitXz1BhlouwcH8bXDwzknz/t4r+rMpmzMZsP7+pLjySfxTwgVPmmjmRAwkWqbdhfVRj8oe1w03sqHVJYaxi/RqVISupda3xuj5c/z9rEF+sOMrxrK8Zdnkav5Ejmb82gyhsXGhDE9w8NZuSbyxjz3irGXNKWG3smkhoTBn9YUBOKnXYlfDgSvntSmSsX/kX9Lgc+pLRCk0nVNdr2jerbfTR43TDvEVj/kbpYp5Gwf6X6TbYfcsznmRCcwIrsFQAkh/mkDwqJg8sehUXPKFNcu8HHvLepsZqsPNP+Vq7J/IGXWrfh3REvI4DyHfO554dxbDIrS0OoLZQSZwkHSw7y3IDn/D4uzcnhVwEkhBgBvAmYgfellC/VOS6M4yOBcuBOKWWT7qw7WJRJangqT/R7gq7RXZvy1Aq3U5mgstarSDKzFXqrUO12Ee2q7dNJIYY5zGKDSx9UWsZrPns/+t4LQ55RC/+uRUrzuGCUuoPet1wJm8gU6HKtSm8z50GVwLIKa5BykKMW3Hmbcnjks40kRwfx9m09uffLzdVdR/VO5m/XdWP8jF95fu5WMvJK6JcaRWSQjY7xoSTWXeB/95mKRAuKVs730Hi4Y7ZKtdN+aO0ILpOJzVlFvLFoJ0szDuMyNt8+d21X2karhd1uNRNmC6cCCLbZsFtVRFxKTDCLH7uc1XuPsHbfUaYu3cOW7OU8cHl7nB4vh0oq6RIfyuBOsQTZLDhcHg4eLScqOIAoIwqt2OFifeZRKpweYkMDKHG4eefHXazNrAmpHpAWzV+vq+30T4uJBkMAtYup8ZF1ig/lzVt78tDQjtw+bRW3Tf2Fewa1QwiocHmuosI0AAAQ90lEQVTolRzJ5Z1iCfC9gbCHK5OUQZXfRkS2hci25JdW8sO2PDLySvFIyZbsYlbvLeChoR14cEiH6h3/A9sl8bWK4OeytETiQu38955LeH7uVt5enMHkxRmM7p3Ew8M6khBu+PhSLlVay+op6ncZ1xWS+sDyNyBno8qDl7FARfFtnKE0Z6S66bjkfqU1rZqiBOaoaSAERRUu9hwu5Yt1B9l4sJA2qTX+xM17A8k8uJ9reiQQarca135Phalf8bSqhuv1QOeRx+S6O22khNI8Yr66n/stXl4JcfF97kqGtR3GO10vY1PuMp6O6MXwoa8QFRxHfkU+To+zacegOS1EvfsvmuLEQpiBncAw4CCwBrhNSrnVp89IYAJKAF0MvCmlvLie01XTp08fuXbt2pMeh3PuI1i3zUHE91D/iKmXqQXC6wGkMuEERinnu8cNHqfx51KPXlf1c+l24i3OwXs0E4/Xi9fjxrb9ayyFe3AHx+Oyx7AxbRxznRexPacER8gCMr2zADDtf4EQazApMUF0iA3hN5ZVRLnyKI/rRaushcRtmVY9ZilMCOnFExCOtNixlOXhNVkxeWvMaMXRF7Lt4heJPLIBe8EONseOZGV5EttyitmRW0JJpZseSeF8dHc/IoJsjPnPDNI9SkBN7r+YwR1jcbq9vDB/GzNW7cfpqfFJXdgmgpToIMqdHkodbqJCbCRFBhIVZKPM6aGs0o3b4yU80Ep+mZNtOcXsOlRafRNe4nATZrcwuncbWkfY6ZkcSe+2kbW+l99M/wf7+ZQozxUsufuter+7dZkF/PGTDWQV1s7WbbeaiAkJIKfIgccoBdA63I7NYiKzoPyYmInY0AAeHNKBwR1jEQISIwKPSeuydOchHlip7vaXjF5HVLCNuuQVO5gw41dW71MBBlazwOWRhNktJEcH4XB5Ka5w4ZWS9nEhRAcHcLCwgoy8EkxCkBQZiNsr2X1YfVZ2qwmr2USwzcKEIe0Zc3Ft7W3hllweXauE+98v/K6WuTSv2MF7S/fw0cpMhIDO8aFYzCbsVhNto4K4zJRObJDgaOvBCJOF5F0fkbrlHQByu93H/k53EbXrK9pumYwXweYOD7Az4RoqXR4qXR6cHknmkTJ+PVBI5hEVCWgzm0iNCWZ3xXICEz8BoGTbi4DAYhJ0ax1G5/gwBgftYfjGh7D4lG2QCGTalZiS+6vAkLJ8pT1GtlU3VuYAJaycZRAco6wEjiJlVi3NU5knLHak9CKzNyAO/IIoOqjObbFT8tuZ3L19MocrDvFs/2d5+KeHuckczbMZ65C2UESbvsoX1/V6df5GIIRYJ6Xs06g3aU4Kfwqg/sBzUsqrjNdPAkgpX/TpMwX4SUr5ifF6B3C5lLLBIjWNFUAfTXuTwL3f082USSexH7No2vlu8qbwuvtmfvQaubmAYJuZbq3D2V6yAhn3IV5nNH3ME4kMsrI3v4ydeaVUuGpXtOwjtjPIvIkd3jb84O1FL1MG15t+JlA4WerpwVzvJcSLAoaYfuWIDGW+9xJcdRTY0AALnRNC6RwfRs/kCK69sDVWw8fx8Zp0Jm4dg9cdyq+/X47NUuPzqnB6yC12UFDmZPXeAhZuzeVIqZMgm5ngAAtHSivJLnRUC6lgmxmzSVDscBMRZKVTq1A6xYdiMZnwSklSZCA392lT23RWh8fnf8K3h1+gT/A4/j36gQb7ebxqIQy0mYkKtrE+s5CFW3MpKHOSHBVEu9hgDpdUsjmrGI+UdIgLoV9KFBFBNg6XVuKVkv7toqu1rIYoKncx8HO1ryb9jvTj5h0rd7qrP9efd+Xz7aZcDpdWYreaCLNbkRJ25JVQ7HCREG6nQ1woAAePKkHaIymcIV3i6JoQdtzrFDtcXPqZGtOym9cREXSsUDxQUM6Upbs5UFCB2+ulrNLD3vwyiipcx/RVSKp+pyciLjSAnskRxg1JMP3bRRMWaOXu/85lvfcppNfMg2mz6JsaxYItuWzYX8iOvBIKypzYqaStyOOoDMUmXNxsXsrN5p9oLQpwSgtHCCeEckJF40uB5MkI1ng7sUcm4JYW5nkvZrdMxBSQQ1DKZITJA65onJkPMNizhUGmTVxq3ko7kc3m6BFcMOGzRl1PCyD/4U8BNBoYIaW8x3h9O3CxlHK8T5+5wEtSyuXG6x+AJ6SUa+uc6z7gPoDk5OTemZknv+t6wZZctmQXI6UkwHmU2NLtWD1leLDglWB3FWN3F+LFjFtYcGPBZTx6hBVhsSHMVoTZhsliwxUQSXlIMharDZtZYDWru1irxUSAxUTb6CDax4ZgMZvYe6SIRxZOZGjyMO4fMLh6sXF7vGQcKqXc6cHt8VLmdFNc4abY4UIAFrMJi0md22JWd5cWk3puNY5ZzAIhBF6vxGI2ERsaQOtwe4MLmtvj5cFvpjM0rRc3dm98OhOvV1Lh8hBoNWMySjJ7vBKz6eQWs7oUlDl446cV3D+wf435qJn5y6KPiAmK5KFL60+z1BxMXDIbi8nCI4OuOXFnAyklBWVOcosdSAleKXF7JaUON063lwCrCZvZRIDVbDyq367NYiLAYlbPzabq77kuZZVuxs55mUtaX8z4S4cec+2cIgfbcoopcbipcHlwur04XB6cLg9uVyUOrxmPBI/Xi91dTLgjCwseKi2huM2BhLiPEuouoNIaQmlAPOW2WMxCYpEuLLhxWiMwmdTvXwgwC4HJeL6pYBU7ilfTNeh6WgXHExJgUdd3ugkr2UnnxCiGDGpcUUgtgPyHPwXQzcBVdQRQPynlBJ8+84AX6wigx6WU6xo6b2M1II1GozkdtADyH/7cB3QQaOPzOgmom7TqZPpoNBqN5hzEnwJoDdBBCJEqhLABtwJz6vSZA9whFJcARcfz/2g0Go3m3MFvYdhSSrcQYjywABWG/YGUcosQYpxx/F/AfFQE3C5UGPaJU01rNBqN5pzAr/uApJTzUULGt+1fPs8l0HAIlEaj0WjOWc67XHAajUajaRloAaTRaDSaZkELII1Go9E0C1oAaTQajaZZ8NtGVH8hhDgMNLYAfQyQf8JeZw96Pi2Xc2kuoOcD0FZK2Yia65qT5awTQKeCEGLtubSTWc+n5XIuzQX0fDT+RZvgNBqNRtMsaAGk0Wg0mmbhfBFAU5t7AE2Mnk/L5VyaC+j5aPzIeeED0mg0Gk3L43zRgDQajUbTwtACSKPRaDTNwlkrgIQQHwghDgkhNvu0PSeEyBJCbDD+Rvoce1IIsUsIsUMIcZVPe28hxCbj2FvieDWSW8BchBDDhBDrjDGvE0Jc2ZLm0tj5+BxPFkKUCiEe82k7K+cjhOghhFgphNhijN/eUubTyN+aVQgx3RjzNiHEkz7vafa5NDQfo32C8b++RQgx0ae9xa4D5yVSyrPyD7gM6AVs9ml7Dnisnr5dgY1AAJAK7AbMxrHVQH9AAN8CV7fwufQEWhvPLwCyfI41+1waOx+f418Cn/v2ORvng8ownw5caLyOPot/a78DPjWeBwH7gJSWMpfjzOcK4HsgwHgdZzy26HXgfPw7azUgKeVSoOAku1+P+keqlFLuRdUf6ieESADCpJQrpfoVfgTc4J8RN0xj5iKl/FVKWVU1dgtgF0IEtJS5GGNszHeDEOIGYA9qPlVtZ+t8hgPpUsqNxnuPSCk9LWU+jZyLBIKFEBYgEHACxS1lLtDgfP4XeElKWWn0OWS0t+h14HzkrBVAx2G8ECLdUM0jjbZE4IBPn4NGW6LxvG57S6G+ufgyCvjV+Edr6XOBeuYjhAgGngD+WqfvWTkfoCMghRALhBDrhRCPG+0tfT71zeULoAzIAfYDr0opC2j5c+kIDBJCrBJCLBFC9DXaz9Z14JzlXBNA7wJpwEWof5rXjPb67LnyOO0tgYbmAoAQohvwMjC2qqmec7SUuUDD8/kr8IaUsrRO/7N1PhZgIDDGeLxRCDGElj2fhubSD/AArVEmq0eFEO1o2XMB9R1EApcAfwJmGj6ds3EdOKfxa0XUM42UMq/quRDiPWCu8fIg0ManaxKQbbQn1dPe7BxnLgghkoCvgDuklLuN5hY7FzjufC4GRhuO4gjAK4RwoHxCZ+N8DgJLpJT5xrH5KB/Ff2ih8znOXH4HfCeldAGHhBA/A32AZbTQuRgcBGYZ5rTVQggvKgnpWbcOnOucUxqQYcut4kagKjJmDnCr4StJBToAq6WUOUCJEOIS4w7pDuDrMzroBmhoLkKICGAe8KSU8ueqDi15LtDwfKSUg6SUKVLKFGAS8IKUcvLZOh9gAdBDCBFk+E4GA1tb8nyOM5f9wJVCEYzSKLa35LkYzAauBBBCdARsqAzYZ906cM7T3FEQp/oHfIIyF7hQdzB/AD4GNqGikOYACT79n0ZFvezAJ8IFdUe32Tg2GSM7REudC/AXlF1+g89fXEuZy6l8Nz7ve47aUXBn5XyA/0EFVGwGJrak+TTytxaCikzcAmwF/tSS5nKc+dhQGudmYD1wpU//FrsOnI9/OhWPRqPRaJqFc8oEp9FoNJqzBy2ANBqNRtMsaAGk0Wg0mmZBCyCNRqPRNAtaAGk0Go2mWdACSHPKCCHeEEI85PN6gRDifZ/XrwkhHmnC630ohBjdVOfzOe9TPs9T6mZWPs5Y9gohxjXRGB4WQuwXQkxuivNpNGcDWgBpTocVwAAAIYQJtdu8m8/xAcDP9byvpfHUibvUy5+klP9qigFIKd8AnmmKc2k0ZwtaAGlOh58xBBBK8GxG7SiPFEIEAF2AX4UQzwgh1gghNgshpho767sIIVZXncjQPNKN572NJJLrDK0qoe6FG+ojhPhJCPGyEGK1EGKnEGKQ0R4khJhpJNz8zEhU2UcI8RIQKFQdnP8apzcLId4TqpbMQiFE4Ik+CEMjeksIsUIIsadKUxNCXG6Mc6YxnpeEEGOM8W0SQqSd6oev0ZztaAGkOWWkKgvhFkIkowTRSmAVqq5KH1RZAicwWUrZV0p5ASqt/zVSym2AzUhuCXALKmmkFXgbGC2l7A18APzD97on0ccipewHPAQ8a7TdDxyVUvYAngd6G3P4M1AhpbxISjnG6NsBeEdK2Q0oRGUdPxkSUAlIrwFe8mm/EHgQ6A7cDnQ0xvc+MOEkz63RnHOcU8lINc1ClRY0AHgdlcZ+AFCEMtEBXCFUWYIgIAqV2uUbYCbwW9RifYvx1wlVaG+RSsuFGZVqxZcT9ZllPK4DUoznA4E3AaSUm6u0rQbYK6XcUM85TsRsKaUX2CqEaOXTvkaqfGMIIXYDC432TajiaRrNeYkWQJrTpcoP1B1lgjsAPAoUAx8IVY76n0AfKeUBIcRzgN1472fA50KIWYCUUmYIIboDW6SU/Y9zTXGCPpXGo4ea33hjSixX+jz3oLS2xr5PNNDu9XntRf8Pas5jtAlOc7r8jDI5FUgpPVIVLItAmeFWUiNs8oUQIUB1FJtUpSQ8wP+hhBGoJJGxQoj+oMxtQtU+8uVk+tRlOUrbQgjRFSUwq3AZZj2NRnMG0QJIc7psQkW//VKnrUhKmS+lLATeM9pmA2vqvP8zVPbomQCGz2g08LIQYiMq2/cA3zecTJ96+CdKaKWjKrCmo8yEAFOBdJ8gBI1GcwbQ2bA15wVCCDNglVI6jMizH1DBAM5TPN+HwFwp5RdNOMY7UabK8U11To2mJaPtz5rzhSDgR8PUJoD/PVXhY1AEPC+EiGmKvUBCiIeBcahKsBrNeYHWgDQajUbTLGgfkEaj0WiaBS2ANBqNRtMsaAGk0Wg0mmZBCyCNRqPRNAtaAGk0Go2mWfh/ZbuG+imLcdsAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "with pt.Environment(wl=wavelengths, freqdomain=True):\n", " detected = nw(source=1)\n", " nw.plot(detected)" ] } ], "metadata": { "kernelspec": { "display_name": "ptdev", "language": "python", "name": "ptdev" }, "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.7.7" } }, "nbformat": 4, "nbformat_minor": 4 }