diff --git a/README.md b/README.md index cb71594..2fbee5b 100644 --- a/README.md +++ b/README.md @@ -4,6 +4,51 @@ This package offers the serial interface for communication with an EIT device fr **WIP** Communication with ISX-3 +## Installation + +### Windows: USB Driver Setup + +To communicate with the Sciospec instrument over USB on Windows, you must install the `libusb` driver. + +> [!WARNING] +> Installing the `libusb` driver will prevent the Sciospec instrument from working with the official Sciospec software on the selected USB port (`FS` or `HS`). +> +> After installation, the Sciospec software may no longer recognize the device on that port. +> +> To restore the original behavior, you will need to uninstall the driver manually. +> *(This rollback process has not been fully tested.)* + +#### Install `libusb` using Zadig + +1. Download and install Zadig: + https://zadig.akeo.ie/ + +2. Connect the Sciospec instrument to your computer. + +3. Open Zadig and enable: + Options → List All Devices + +4. Select the correct Sciospec USB device from the list. + +5. Install the driver: + `libusb-win32 (v1.4.0.0)` + +--- + +### All Platforms + +Clone the repository, then create and activate the Conda environment: + +```bash +conda create --file environment.yml +conda activate sciopy +pip install -e . +``` + + + + + ## Contact If you have any ideas or other suggestions, please don't hesitate to contact me. diff --git a/environment.yml b/environment.yml new file mode 100644 index 0000000..a8813b4 --- /dev/null +++ b/environment.yml @@ -0,0 +1,13 @@ +name: sciopy # This will become the name of your environment +dependencies: # The list of packages to include in your environment + - python=3.10 # You can specify versions + - matplotlib=3.10.8 + - numpy=2.2.5 + - Pandas=2.3.3 + - pyserial=3.5 + - setuptools=82.0.1 + - tqdm=4.67.3 + - pip=26.0.1 # Install pip in your Environment + - pip: # Include pip dependencies last + - pyeit + - pyftdi diff --git a/examples/test_EITphantom.ipynb b/examples/test_EITphantom.ipynb new file mode 100644 index 0000000..d6f8779 --- /dev/null +++ b/examples/test_EITphantom.ipynb @@ -0,0 +1,458 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "id": "1e1715cc-c123-45e2-b767-e6412faeebe3", + "metadata": {}, + "source": [ + "# Example code for connecting a Sciospec device" + ] + }, + { + "cell_type": "markdown", + "id": "8a5423a4-f8da-4cdb-8163-1d50d12e6e23", + "metadata": {}, + "source": [ + "## EIT device\n", + "\n", + "**USB-HS**\n", + "\n", + "If you have issues with the connection try [this](https://eblot.github.io/pyftdi/installation.html)." + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "id": "f49a2bed-51db-4627-b5e6-21831a620644", + "metadata": {}, + "outputs": [], + "source": [ + "import matplotlib.pyplot as plt\n", + "import numpy as np\n", + "\n", + "from sciopy import EIT_16_32_64_128, EitMeasurementSetup" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "id": "f8c7590a-9aea-4b22-91ad-b953b540362b", + "metadata": {}, + "outputs": [], + "source": [ + "# create a 'sciospec' class which represents the sciospec EIT device\n", + "n_el = 32\n", + "sciospec = EIT_16_32_64_128(n_el)" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "id": "aebcee3a-f3f0-4709-88d3-525eec9b6bb8", + "metadata": {}, + "outputs": [], + "source": [ + "# connect device via USB-HS port\n", + "sciospec.connect_device_HS()" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "id": "de7aaeb5-afb9-4c0d-a3ad-1a8bcf804d6a", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "No message inside the message buffer\n", + "message buffer:\n", + " ['0x18', '0x1', '0x4', '0x18', '0x18', '0x1', '0x84', '0x18']\n", + "message length:\t 8\n" + ] + } + ], + "source": [ + "# read system message buffer\n", + "sciospec.SystemMessageCallback()\n", + "# should be empty" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "id": "e8609e39-a7bd-4aaf-8f02-56d7c3c527d7", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Command-Acknowledge: Command has been executed successfully\n", + "message buffer:\n", + " ['0xd1', '0x14', '0x1', '0x0', '0x19', '0x2', '0x65', '0x10', '0x3', '0x3', '0x0', '0x21', '0x0', '0x8d', '0x0', '0xba', '0x0', '0xe4', '0x0', '0xba', '0x0', '0xe4', '0xd1', '0x18', '0x1', '0x83', '0x18']\n", + "message length:\t 27\n" + ] + } + ], + "source": [ + "# get device info\n", + "sciospec.GetDeviceInfo()" + ] + }, + { + "cell_type": "code", + "execution_count": 30, + "id": "10d15631-d005-4929-9f4e-38716f0fdfb2", + "metadata": {}, + "outputs": [], + "source": [ + "# create a measurement setup\n", + "setup = EitMeasurementSetup(\n", + " burst_count=1,\n", + " n_el=n_el,\n", + " exc_freq=1_000,\n", + " framerate=3,\n", + " amplitude=0.01,\n", + " inj_skip=0,\n", + " gain=1,\n", + " adc_range=5,\n", + ")" + ] + }, + { + "cell_type": "code", + "execution_count": 31, + "id": "9a80f291-4722-497d-9f16-672b6b50af86", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Command-Acknowledge: Command has been executed successfully\n", + "message buffer:\n", + " ['0x18', '0x1', '0x83', '0x18']\n", + "message length:\t 4\n", + "Command-Acknowledge: Command has been executed successfully\n", + "message buffer:\n", + " ['0x18', '0x1', '0x83', '0x18']\n", + "message length:\t 4\n", + "Command-Acknowledge: Command has been executed successfully\n", + "message buffer:\n", + " ['0x18', '0x1', '0x83', '0x18']\n", + "message length:\t 4\n", + "Command-Acknowledge: Command has been executed successfully\n", + "message buffer:\n", + " ['0x18', '0x1', '0x83', '0x18']\n", + "message length:\t 4\n" + ] + } + ], + "source": [ + "sciospec.SetMeasurementSetup(setup)" + ] + }, + { + "cell_type": "code", + "execution_count": 32, + "id": "4ca492ed-b55d-429c-b1a3-996691918529", + "metadata": {}, + "outputs": [ + { + "name": "stdout", + "output_type": "stream", + "text": [ + "TBD (to be checked)\n", + "Command-Acknowledge: Command has been executed successfully\n", + "message buffer:\n", + " ['0xb1', '0x3', '0x2', '0x0', '0x1', '0xb1', '0x18', '0x1', '0x83', '0x18']\n", + "message length:\t 10\n", + "TBD: Translation\n" + ] + } + ], + "source": [ + "# look inside the docstring of the function and manual\n", + "sciospec.GetMeasurementSetup(2)" + ] + }, + { + "cell_type": "code", + "execution_count": 33, + "id": "d88f4bf1-d618-4024-a533-ce441f141062", + "metadata": {}, + "outputs": [], + "source": [ + "#sciospec.update_BurstCount(4)" + ] + }, + { + "cell_type": "code", + "execution_count": 34, + "id": "7607a215-8bb9-429a-885d-e96707d957b5", + "metadata": {}, + "outputs": [], + "source": [ + "# get data\n", + "data_0 = sciospec.StartStopMeasurement(return_as=\"pot_mat\") # or \"hex\"\n", + "# data.shape = (burst_count, n_el, n_el)\n", + "# check if data shape equals setup\n", + "assert data_0.shape[0] == setup.burst_count" + ] + }, + { + "cell_type": "code", + "execution_count": 37, + "id": "2c4c8ac0-4a41-4d85-836d-b17b0a8b81aa", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAi8AAAGdCAYAAADaPpOnAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAeKNJREFUeJztvXl8VtW1///JAJlIAiGQAcKkMiioCA7grC2Kyu3g12prHTp4r7XWga+1tfbban/X4r3Xerm9tXptbdFqv3pbtN/2aq20CmpBZVSQWYFESAiEkJlMz/n98eR5OE/yJDlr73P2enbOer9eeZGEZ2Wds8/n7LP22mvvk+Y4jgNBEARBEARLSOc+AEEQBEEQBAoSvAiCIAiCYBUSvAiCIAiCYBUSvAiCIAiCYBUSvAiCIAiCYBUSvAiCIAiCYBUSvAiCIAiCYBUSvAiCIAiCYBWZ3AfgN5FIBAcOHEB+fj7S0tK4D0cQBEEQBA84joOmpiaUl5cjPX3g3MqQC14OHDiAiooK7sMQBEEQBEGBqqoqjB8/fsDPDLngJT8/H0D05AsKCpiPRhAEQRAELzQ2NqKioiL+HB+IIRe8xKaKCgoKJHgRBEEQBMvwUvIhBbuCIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC+CIAiCIFiFBC8pzt7DLVi/r55sF4k4WLv3CJrbu8i2nd0R1DW3k+0AoKGtE5GIo2Rrmn11LfivVR+hRaGN3t51GJs/aSDbNbd34c+bq9HW0U22fffjOqzefZhs1x1xsH5fPdq76D6rG9qwrbqRbBfTn0rbNh3rxO7aZrKd4zh4Y3stahqOkW1rG4/hr1sPKml37d4j2H+0Tcnn79ZV4Vgn/bqs/ugwtuw3qz/HceA4Zu/t7oijdKxAVEcqx9vW0Y2qI61KPvcfbUNHV0TJVpVuxf7WcRwcbe3w+WjMIcGLRzq6Ipi35G+Yt+RvaDrWSbLdXduE63/5Dt7bc4Ts96JHVuLqx1ejso52M/32vUpc88QaXPtfa8g+F/3n25jzz3/Fx4doD5Bt1Y047cHX8PVn1pF9rvmoDj/4f1vQ2kF72HV2R/DN327A06v3kn1etvRNLPnzdjz0yjaS3YGjbfjyU+9i0c/eJvv85nMb8I3nNuD+lzaT7Dq6Irj2yXfwpV++S9bf0r/uxNWPr8bdL2wi2QHAvCWvY+F/vEV+OMf19yRdf/OXvI5PPbqK/HD+85YafGXZWpyz5G9knxf+20p8/Zl1+P2GT0h2W/Y34Jon1uDch18n+/zMY3/Ht3//AR5dsZNkV93Qhi/94l1c9Z/m9Oc4Dr781Lu45ok15IBg3d4jOPvHf8Urm6tJdgBw+dI3MeMHr6KRqPnVHx3GrAdewwN//JDs84J/ewPn/+sb5KB9/b56nPvw6/j8438n+/zvtVX47GN/x6Em2qCxoa0TZ//4b0r39oN/2orTf7QCr31YQ7KrbTyGi/7tDSz491Vkn34iwQuB6oZjqG44Bmqc+/Wn1+Hvu+vwBYVAIsbuQ02kz/9+fbQT/vAAfdS8vSbqi9rZPLNmHwDg9e21ZJ9f/MU7eGbNPvzs9d0kuz+9fwAvf1CNHyp0Usc6oyOkdz+uI9mpjLJjrNp5CADw4sb9JDt31qTpGC3A+683PwYAvLKZ1km52VFD01FMf1v20/XX1JOtWbmDpqO3dh0i+4rR1pP9iF0fr6hkRWNU92SI/rbtIMnuk3rz+mvp6Mbfd9dh3b56HCBmtr7y67U42NiO257bQLIDgF09Gbh1e2kDv3/7yw4AwNM9fRKFWABBvS6/X18FQE3z9y7/AJuqjuKRnuP2yv/btB+Hm9vxEvF6AsCyngHfv7y6nWTXGXGwt64Ve4kDar+R4MUA1Js9GWlII33ej+RuWhrNpx9e9xHTtdQHeTKo52k4cx716fqe4bJYoT8/rgu5aX1waoX+XE6pbdTebXYaBfBJC+QbTZ9mYuaZUwvmWycRCV4UIAvGjhKQlIDl4aH9F8xCDSSE1Ma2q2lD8GwrHGepGjwzxHYJSPDiEe4LZQO2jnYEQRgY28ZfHMfL0f+ZLqB2wx1QSvBiAMefoQfRqXlRc0+nqEKNlzg6DB2Xoj/vkB8e+i7JcOuPPoUo+gvMpw9/Q7VpuceZEryoQFQMx1y8H3CI04I+ih0brosvPuU8UxK5LgH6NO+SPniD1LxYBfeFosIyIrQuuawGy1nqZF4sHRHSfcpoOxWdhqWN/NAfPQus7ZLM8cyLTBsJHuBZeZH6qyD8wIbzdHeMDPWSodEffRTqg0/R36CERX9UOPSXKkjwogA1wuaYl/QDlrQpx8PDtntX0tnB+TTv0j79MSDXJThU+1zu5pHgxSPcKTIqtk7h2NDKtk1PsBR4crSRcY/haVthcGwdpFKJa14Kdu2D2l/5s5EV1ae2S5b0MBWWpYLcm9Qx9BpW6C8sS/UZ6xwAkC+MrQGeBUpgatso3O0jwYtHdC6UjJO8w/HwsOKB5cKywxUGwTr9WRA824oNWpCCXWHIYeuS8KGOBM+CLjoasnXDOBt88tCzVFqmjeyD/HYASyvfVdf/myQks0Za6WHRH8En8fNh1B/PFKLoL6lPll19e+z0XWshwYtHuKNMnpUXqb/cIyxz227CsvLCBv1Z6lIL0V+gTlPeZbzmRaaNhCCwbSMrmwjLaFuHsOhPso2piegvQJ+SebEXjgpvlrctM+y5wgFLelgj9c5S1Ez8vK36C830hOhvcJ+iv0Ds/CLQ4OXNN9/EokWLUF5ejrS0NPzhD38Y8PMrV65EWlpan6/t27cHeZie4E6RsYe5hrDp5tXyqWHLsjw8LPqzYGNGX3xq2HK+yVjwHxtqG5ORGeQfb2lpwWmnnYavfOUruPrqqz3b7dixAwUFBfGfx4wZE8ThpTxh7CS4Y0QvcG/PbsynZfrj3PPCqE+7Losyor/U5Php8nbWgQYvCxcuxMKFC8l2Y8eOxciRI/0/IJ8wJVGt19CzVKHb+fAgt60fPi2Yk7NOf/oueVZ7WKA/9tE2S0Ys9fez4QyeuQeaKVnzMnv2bJSVleHSSy/FG2+8MeBn29vb0djYmPAl2EtYVjJQcfr5XrAfG/TnFp0NgzdbYen/yMFzzz4vQRwMgZQKXsrKyvDkk09i+fLlePHFFzFt2jRceumlePPNN/u1WbJkCQoLC+NfFRUVBo84WGx7SHG8JDEkmVoWpGmFGCwjfB3bkKzI4tyMjzvoDnTaiMq0adMwbdq0+M/z5s1DVVUVHnnkEVxwwQVJbe677z4sXrw4/nNjY2PgAYwpwWhtDGXphklUbK0FsaJtQ6g/G4oXRX8ebEV/gfmM2zHnXlIq85KMc845B7t27er3/7OyslBQUJDwFRSmI02Nd6H5Qlg2hrKhbd2dlLHg2fW9DW3ki0/RX1JEf4Z8WqC/VMm8pHzwsnHjRpSVlXEfBgs6nYSkTQP0ad6lde9qkesSoE/zLkV/Xnyad8ma7eaueQl02qi5uRm7d++O/7xnzx5s2rQJRUVFmDBhAu677z7s378fzzzzDABg6dKlmDRpEk455RR0dHTg2WefxfLly7F8+fIgD5OMqVSx248dKy8srZJgKPHnWGFCJYz6YxlNiv768Sn6M4KiU+69zwINXtatW4eLL744/nOsNuWmm27CsmXLUF1djcrKyvj/d3R04J577sH+/fuRk5ODU045BS+//DKuuOKKIA/TM2kwexMn7mpp0HHcZ+ovFfTFp4at4zhq7UQuTHal7U0Fz6I/Mz41bMOiPw64H86msKH+KRmBBi8XXXTRgNHosmXLEn6+9957ce+99wZ5SKHB2iyIBYTlfTY62LrnD9kny3laJoaQEJarkirnmfI1LymJsdVGGrY++Gd5t4fGiNCcT22X7PPFXrBNf/7cl6n/egoO/XEHzzb0C37AcZ5UYj65E1MSvBDgTCPa8LDzh9Svtndjbtm8+wczPt2ER392wRJUMGg+PPpL/anS2GWR4EXoF+vSw5atSBAGxjr9CYFh2x5OLMoNSf8XXyot+7zYB8v22AwpdLpLhikcbY90bN2Mioroz4tPfcKyGRoV0Z8Xn/rQQxCZNrIO09cqUZgWTMD64dK4R73pQNXOQ+dFh8aC54SfRH+B+RT9JfeZ8JPoLzCf5Bfw9tgFcCwUJHhJYXSKsWzNSFjh07xL61bSiP4C9GnepejPi8+Q6Y97KbkELwqwbI/NUW1P/Lyt6WwbNqPiXnos+uvHp+gvMER/Hnwy6E/Xzi8keCFg/N1Gli1PtBWd81TtJHWm8FlWODEQmnfLiP6S+wyh/jjgWKrvBxK8pDLM+2wIyQlLSto2/XEXlRrzad4lD5adaFhW5zkpUvQiwYsCLO82otr6Um1PfZ+IHylMhk3CWDYmY9mBjejRMv2xpO198GmB/rhH2xx3C88bnlP/usRrXsy7TkCCFwKc69q5i6NMwfNiMnXTsGQYwqI/2zCnBMf1HcPgLST6s2GqNJ54kYJdoT+4RztUJIUeHGFZySCkJrbpz7bjVfZp3mU8qOQOJyV4UYBltRHZNiQbJtlWDNJD6k8aif68OWXYDU3058FW9BeY03jmRd+1DhK8UDC+2sidNjXrG7BjwyRffGrYKvcdGhX+5lZ7iP6M+NSwFf0Fh+gvOcdrXmTaSOgHnT5C0qYB+mQ5T959Nsi2or/gfIr+BrcV/QXnUzIv9mJKMNxzvnRx2rlU1Yb0sG2dsejPO6I//32K/rzDHYSoIsELAfPvNnJV+DN0HtxpQVPonKfqiFRnqaqx4DmE+rNhtYcb0V9wSP+XnFTZz0aCl1SGZbSTGsJMZcKStrdNf2G5LqG5RUV/KYkslbYYUw94vTlfNesEM4Zto1ne7WGDz3DELnr688G/HRvj2al5G9DRnx/YcF1kkzoL4Qw0zaVqebFtqoCnFsCu4slU9DNU4HjPEMfgLSz6s6H/i11/7loZCV5SGNsejOwjdcFX9AomVTN/ckGFKLbpLyz9XzzzIsGLfZjbpM78jaS1MZQfr6FnmCqgwrIvlGX1HKK/AH2K/jzY6tuJ/gZ2yl3QLMELAdMXiyVV63Jjw4ZJtkJe7cGw8kL0Z8onx1SB6G8wn6K/5MRfDyCZF6E/ODZp4hjt6MDz9mzz2Fawa9vxKvuU8wzOp46tZcer7JNx/x3ugaYELymM3vyrDylMqs+wrLzwIz1M/DxPx2hefwlvEqbaiv68+yR+XvTnwTY0+lM09BkJXgiYvlYsG0MlpE0NObUY1Q6L44FFJYz64967goroLzifor/URoKXIYos2hha2HY5w5K2t+7CCIERFv3JPi8Ww1GwZgqdkYc/le/mnepsla7hlcEnDW6fYZlOsUJ/frg0iHLNn2SeB0X2ebEQzmtlw/JsP7CtvwjL9vBh0Z9tmGqvxJU/Rlzy+Axh/0dfedZj5/+hkJDgJYVh2edA4961dedX7pvQGxxta5f+hCAR/Q1ua2f/R/YZW20kBbv2wTHaMWWbuEkTyyt2SfB0GIwV/j75p/o0ZRtG/dGnqnj1ZwrR3+Bw6C/WStyDPgleCJiONJ0BfgrMp+tmYKm2Z78laCiv9qD6UXOjhegv9eHQn7HgeYCfAvMZQv2pBrJS8yL0i15ULbUYXhnqo1B1n6I/Ez5Ff/35FP2lpM+ef7kDfQleFDBXPKZh60PalG7LkM7W9mjJZlQup6K//mzDpz9V/1bs86JjK/rz7pP4+eOb1PngXAMJXgiYvlbcFf7cBVk2oPzwsGCpquhv6KKjvyG92ohZfzasNorb+XwcVCR4SWnMV9uHcaWq6lJBk/CswrFLf7a+bdkG/fHAOCeiYhqaaaPUUKAELwqYunQs22Nzv9vDCp9O0u9JPjXeYcKhP1Po6M8PWN4tQ/w8i/5S43nlGeWl0sz6swEp2LWREKiZu5PiviGE5AzlgEkYHO7gOSz641nhRCMUBbtvvvkmFi1ahPLycqSlpeEPf/jDoDarVq3CnDlzkJ2djSlTpuCJJ54I8hBTGr2CNTVrjiI5HTiK5FiwK4MeGv2FBtGfB58axqo+WVZVRX1y95uBBi8tLS047bTT8LOf/czT5/fs2YMrrrgC559/PjZu3Ijvfe97uOOOO7B8+fIgD5OMqcp7noI19X0O/EmhW/DOn36+p6C1z4YF+lNO28s+G4PCrT9TaG1Sp+zTLv1xrHzUtfOLzCD/+MKFC7Fw4ULPn3/iiScwYcIELF26FAAwY8YMrFu3Do888giuvvrqgI7SO8ZXG3EUTLqwoSO3FZ6X8dHgLsyjr4gJx94efqCjP2PBs3uqiqGdOfTH0/+paWFITxtRWbNmDRYsWJDwu8suuwzr1q1DZ2dnUpv29nY0NjYmfA0VuLfHNuWTH4bNFYjY9mC2VX82PDxEf8HZcuvPBmJa4B5oplTwUlNTg5KSkoTflZSUoKurC4cPH05qs2TJEhQWFsa/KioqAj9OloI1hqkCDnRWXij71Fn5w7HnhRmXoj9TPi3Tnyk4XpJom/5Y+r8UCdJSKngB+qYzjxcHJW/h++67Dw0NDfGvqqoqY8cWNFrTRsoFa07S773ZCl7h2EmTSlj0p/MA4Nh51g846myoJAZpDFNVHPqzaGM87k0kA615oVJaWoqampqE39XW1iIzMxOjR49OapOVlYWsrCwTh2ccjoI1lpQ08wOAY28ZKizbs9umP9WpArcdw4Zxor/+fBp3yaI/HWwrpPaTlMq8zJs3DytWrEj43WuvvYa5c+di2LBhTEfVF+4VLkPBT4JP7oeHzgBC8QC0Rr6iv8B82jFt1N8Pwfm0DeUsiPv7ITwlp0PscLklFGjw0tzcjE2bNmHTpk0AokuhN23ahMrKSgDRKZ8bb7wx/vlbb70V+/btw+LFi7Ft2zb86le/wlNPPYV77rknyMP0zFC/4QHN5dk+FMmFoImFAeDYHkAYnNAEz8znZkP/lyr7vAQ6bbRu3TpcfPHF8Z8XL14MALjpppuwbNkyVFdXxwMZAJg8eTJeeeUV3H333XjsscdQXl6On/70pymxTJoDjjfC6y3PVi2S06k5UDZV92neJdN0no6xqpnoLyUR/QVmy60/1cwfd6AVaPBy0UUXDXhhli1b1ud3F154ITZs2BDgUfkBQ/HYUF7V4vreto3xVP3r7SEh+vPVp+t72zYm49GfGTj2ubJNfxz9X8wld8FuStW8pDrGN6njKJjUQHlvBXexPXs8bwYblieGUX9kW0vX2Onpz1DwbNkUYlj6v/g+L8zHIcFLCqOVNeV4t4eynW1p03CsvAij/mzYnl3058HWUv2RbTn0F8+86PvWQYIXBTjeM2Sqw0rwybAxGcf+JzqjUNUm0nu3kZpPKqI/M9imP1Po6E992sgu/bH0f8ctffCujgQvBMxvUmfelmNXyzBCf7cMR82BeVvb9Ger4nX0Zyx4NuMm0Sez/rinYihI5kXoF453e4QRevGieezryP07Ds8+Ve0Sag7UbVUR/fXj07I+zNY24pgq9QMJXhRguXYMqz1Mwf3w0BnuqLrXOU9zl8i8U56gJ0V6Y68kTKeowbAvox5DWX+2yU8Kdu3D/MUKyfJE9ywq9x0hsGJdIBESeB6w5gtZuZeE29D/ScGuMCgcSzg5pgr0OsawDJXsWmFiq/7oHbIveXsLEP0FZcve/ylvUicFu9bBUbBm63yqJz+u7+mFhPr+9Vb+mFnuwbHnhejPDKK/fnz2870pRH/9kCKvB5DghYDpi8UyClB36cvyxLBgx1JV87a26c9W6Yr++rFTdxka/cUzLxK8CP2hJWoGn6p3klaH4UfmhSHbQ4XHp+gvIJcJiP768yn6C8hlAqr6k2kjCzG2YZf7e5ZUraHz1Kg58OfdHkSffqz2oHYYcJJ+HySiPzOI/vrz6fpe9BcYyvqTzItNGN6kjiGHyfI+G8tS734cLsuScJM+Q6I/ntUp+oj+/PcZFv3F4K4zl+AlheF4Pbt1kYQPWPGSRPMurdOfH8tjrXi3kegvUK/KlmHRn75LX5DgRYGhvJEWz6oC5ptXw1Z5sQf5geVK25vSn1sLQ1h/tuHHdIrO9IT0f+Hm+D4vUvNiDcavlWUjLJ59NswjfVpw2LaqJSywrAhkmE4R/Q1OfLUR61FI8JLScAcSZFuGDoOns7Gst1GEexUE2acP+mN5qy/x86I/D7aiP8+o6o97oCnBiwLGplOYU5jmztM9bUS11fevk/5U7bB0pseGsv54fKrrzw9Ef2H3aaf+JPNiEeZnjXQKJhl8Khes+YOpESlPkRzHqhbRH8Wn6M9fRH80n6b1JzUvQr/o7ZekmMLU8mnWDvCnU6XvNsrQkVuWeRP9eUf0579P0Z93VPUnmRcLYdmkyYjHXpG8KZ8JG2CRjdV8JlQJU231/VNJvC6iP1996uhP1afob3Cf/XwfqE/RnwdbNTu/keCFgG6WTGd5LBWWtKlxw15/hvB3WDajckHf4dQ8oj+arejPX0R/NFvT+pPXAwj94tdcqAk723CfJsfeMlT8mM82uZKBJYXO8W4ZH+xEf/341PZIJyz6c8MdhKgiwYsCtnYeqeozIT1MdKn88HBX+Ots2KVop/NWX5ZN6oaw/hJ9GndJxo9pBo6NGXUQ/fnsR8c2XrDry6EoI8ELAeMRqhTJ0f6OP3/Ggx8LnnCWYlthaMLfIfm0S0MsR2tZ8awqtulPCnZDiA2rCjjQW57IkKplyWfr/wnRX3K49WfDu41Ef8Fhq/4k8xIiyEVyoal8T+7flE8r3m3kx5JI0V9yn6K/wf2I/oLzaan+uGtlJHghYDrStC2FybJJU4J/73+JY7TjxobRtm0+RX/eEf357zNs+pPMS4jQKdI0hd6NpFrzonMj6dtZ0ZH78DdEf/3Yif4G9+nD3xD99WNnqf4keLEQY6lahrwpj8/k35uCvFQ14Xu1A+Z4t4zorz+fyb83heivP5+iPxOo60+mjazB/LuNzGNbIavTz/eD2jGkh93Y8W4Z84j+grWNIfrz32fY9CeZFwvhGO2YwrIVnMq4r6EVHbkfD6whrj9VU56NyUR/JnxSCYv+3KgGIbJU2iI436JprsNyV9ubr3w3hc6cb39/h4LWJnVDWAs8Pu3ClykcHf8sWRjRn5/ovQYhNdpFgpeUxq4RlvLIw6eVF8Z2p0yNe1fwCb86Y4oubNMQz+Fa1kiK2KY/mTayGOULb8GqAh1YXkNv2E7XVt2pL3l74y6p2LbDsx/6s+GdP6I/L7bh0p/s8xIibHifDb9PU1NVCbcv1fj4t6oHoPFWX5apKtGfvz5Ff0Z8UuH3aY/+JPMi9It1Iw9ln/44paRfuTMv5I7cMi1w+BT9eUf057/PsOkvFAW7P//5zzF58mRkZ2djzpw5eOutt/r97MqVK5GWltbna/v27SYO1RPKs0ZWrCowf/fq3Uj6qVorNgnj2FVV9OfBVvTnFdGfr2Y9tnz641zAAhgIXl544QXcdddduP/++7Fx40acf/75WLhwISorKwe027FjB6qrq+NfJ510UtCHOii614pjYygqTj/fB+qT4TzdJ6fzwjjVzk6nzkG95Er0l9Sn6M+Dz+Tfk3yK/pL7tFh/nAQevDz66KP42te+hq9//euYMWMGli5dioqKCjz++OMD2o0dOxalpaXxr4yMjKAPNeXgkAhH2lTLUvHG11vhpGwax4q0vXmXoj+iT1VEf/34FP159jmka146Ojqwfv16LFiwIOH3CxYswOrVqwe0nT17NsrKynDppZfijTfeCPIwyRgb7TD0GLZugKXjk/sm9IIvdQ5DXH8sqz1Upwp00vZqLrUQ/XmwtUh/btQ3qePtODOD/OOHDx9Gd3c3SkpKEn5fUlKCmpqapDZlZWV48sknMWfOHLS3t+M3v/kNLr30UqxcuRIXXHBBn8+3t7ejvb09/nNjY6O/J+GC8yFn7DbmSGEmd2/Mj9a7ZRQPWGe+mKdTN89Q1p8OvmxSJ/ob3OcQ1h9PhslfAg1eYvS+URzH6ffmmTZtGqZNmxb/ed68eaiqqsIjjzySNHhZsmQJHnzwQX8PWFBCfat0HZ8axlb5TJUuI2CG8AOD26cOor8h5dKXPpc7Yx3otFFxcTEyMjL6ZFlqa2v7ZGMG4pxzzsGuXbuS/t99992HhoaG+FdVVZXWMXtBvWCN6MeyG0nVVm8UoJqq1ZmhtrMjF/31Y8esP3rmT/QXFKI/7z65Z9sDDV6GDx+OOXPmYMWKFQm/X7FiBebPn+/572zcuBFlZWVJ/y8rKwsFBQUJX0GhO8dHTdX6saqASkJK2lAnmXieRlzqnRlzwaRy8Cz6S+5T9De4S9FfcD4t1R935iXwaaPFixfjhhtuwNy5czFv3jw8+eSTqKysxK233gogmjnZv38/nnnmGQDA0qVLMWnSJJxyyino6OjAs88+i+XLl2P58uVBH2rKwTPyML8Kwq+CNVK1PfOcrx17Xhh3KfrzYqtuGkf0149P0d/gtj3/cu/zEnjwcu2116Kurg4/+tGPUF1djZkzZ+KVV17BxIkTAQDV1dUJe750dHTgnnvuwf79+5GTk4NTTjkFL7/8Mq644oqgD9UzxrbHtqyew48Nk+g+Ve04VhWY7xjdiP6C8KlqJ/oz4ZOK6M+DrZMa00ZGCnZvu+023HbbbUn/b9myZQk/33vvvbj33nsNHBUd7UCTYWMoKgk+jaUwnaTfB+xU2acfox2OFU6iv/58iv4Gte3n75AQ/fXj01L9DeWCXUEPliWJDKMz39KmhJuQpTDP/UNItoQ36VP0R7AT/fnuM2z6497nRYIXJbxf+o8ONce/13mrqik4fOphPlWrCvcKp6GuP57VvAxTBcoPSdFfkD5t0p8fcBfsSvBCQOVabdhXH/9+/Khckm0k4konGtJod7c7hUnDfbwUujTOs1vZZ8QXn5RO3X2emek0Nbmvi2qHlZlBu927GfSXoHmqreJB6uivi0F/CbYW6S+D6tO2/o9Bf+r9n77PjKH+Ysaw09F9vKPJHkZrbretKXR8qtp2avjs7Fa7CTu61H22K9q6z3MYMZDwQwvDMmidjU4bqaLjU1VHevpTC0I6u9Wfxn6cJ4/+iD4Z9OeXFijonGeHov786P+GZ/KGDxK8KGBKJLY9PFRtWXwyd1Ly8OjHp8Z5qgaVWgFTl1oQwhE8c+uP+rDTCSRU4dCfVsCkqD8/gjSqhvxGghcCKuva3R0GVWaJN5KZvGmCT2I+UXU06V8becd905N9dimOdlzHSsygK/vscvmkPjxY9Neloz/9h7op/Wn5VNSC+/60QX+JwYAZ/bWL/jzbSuZliKMT4VqXeVG9kfxK25N8div7VG0jvwImkp374THEMy/KwbNP9yhNf7yZP9FfcvyawqagmrEB1K+pHz4leLEQUifFPG1kss5BtXiMo410bl7lh4dGwKRc5+AKmIZ6nQ2H/ljqHJQDCbv058e0EbUwmUN/HHU2fvgcTry3/UaCFwIql6pdsaAKUC8GdC+J1OkwTC3CY58qIJ5oh2IbsYywXA8s6mqPTsU20tGfTnZAFRb9afQLqvrrUKyPiNqa15/qPerWH3mqNCz686H/k8zLEEe1oArwJ1VrQ5Gmf6uNvLe1zmoP9cI89eWJfsxtk219eDBzjLap+DW1S9IfR5GwTwETyU70Nyi+TVsSLqovq40yMpT/hh9I8KIArUhTo7bChw5Dr0jTDNatMGGZKlBNSTMEzzr6s6zOpoOhzqY9JPpTnxLxJ/NnCu4Aj8rx1UYybWQPCtdKr0hOrcAz4eY1lKrVQaeN/Kk/oXnlGO34sZ+DzkoG1Tovqv5Up1l1do9l0Z/GShrlbI9PeziZ0l+7om1C8Ex8wOpM86uiM1WvvHhAYyWhFOyGBK2Rrw8dI3V5t6rPLobsCaB+43NM5+l1GGoZPI7CULfmqfpTvS6qxZJAOFcbDWX9afV/ij75+j+1Aa6sNgoppkbbqvPi3A8suq3541VNvQPqx+tXJ0WBY4WJXoZJ8SHJtCWBTauNwqI/ntU75gdDAG8tEnXpu99I8EJAabWRRiW5cqqWuQDM1AZYgPq7ZbSq7RVH2xyFoQm1ClT9KW/G588eOrQC2OOfNbWqBRD9DW5nXn++reoz1P+prioF/FmeTV4FK5mXcMAxImQZSbrsMtPNFMn1fhGksdc2cCxP9GHlmSlbnSW5yht9uZfkUqeqFM9TR3/c79AZ2vrjyOQe119Y+j8JXiyEVlthVyW5Lz41Rh7KPolwTFXx1zkQffqwJJds60fArpH5I9kxTVX5U3NF9BkS/akGzwl2VP1x9PM+LJWWdxtZhNK7jTQqyWsajrls1VYVUFdhKBcgJkyPkVz2aiNDqyf82qTO0GhbdeWFzkqGdsXropp6B3TO0z09QXKprD+/BglU/amuiOEoTObQn07/58cKJ+qJdjLoT2vaSGpewoFqhPtJfSsqj7Qa9QnwvBnVt71aGDpyVTudZcuqPk3Zck/JUWHRn6JPx3GsfQGgKVvRnwefPvR/WTJtZB+USLWtU614cdfBZiU7ADim6FPHtq2DwafGeaoeb0dXBKqrcrXaSNFWp41UrwuL/jR8tneqdeQcmm/viijvPyL682CreJ7W9bk+9H9ZmbLDrjWorDZyC4zS6bT2EhfFtsXt07tZH780n10unzSvqj57txEF1ePt22F4t21RPF7HcdDaqXhd2l3nSRSDqm1i29Jo7VDz2dqurr8WRZ86D+bEe9S7076aJ9gqHq99+tPo/zr90DxRf+3m+79WZf0dP8/cLAlehjRuUZuwAxI7crKt6ghL50ZqV7Nt6XWeekGIR58a16Ut4SHp/VjbuyLKSyK1rotFWvCrI6fAoj+de7vdvP50rotqG7Vp9Z2K+lO0A/zri/SCYJrd8Ix0Kdi1ESMiYQhAoraqwZaGz041nzoPSdXj1cv2mA8GtLI9HFpgCCpVz9OvzB8FnWwPSzCqoT/lgZ9GIKHqUzWrpeNT77qoaj5qx511ASR4IUFdbNTVHelVJKee4qWMe1RT70DiCILkUyM9rOpT6yGpeLy9b3rSNAPLCEstnd3eq7aHcl3aNPSnOuWkOvXYx1bRjory9GzvbI+BaQbr9Kc4xQVoaEGj/2tT9KmThWtR1F+sfXKHSfAypNGJxnXSkByjHY60vU4goZs2VUH1wcwx8uXIMDmOo3yuWvpTnqrS0R9HtlE/MKTCoj/Fh3o026ga4DFkezTqIpWzqj33Sm5WppK9n0jwooJHkXAEIDq2vUc7NJ/qRZo8AZNe2lTJVrW2h2MaUGNUp7ySqzuivN2+3jSrPQE7x5QcyzQ0w7R5R7dObZlOkbo9dV6x65k3XDIvVpFGXG/k58iMUmCnmqrtG8mbWUmjuiLLrw5OZxRKslWs7enzYFbMMOmtdiNoQXHlT+/gTlV/xlbYaRXPquqv9xSOuq1nOw39+bWSxkj/55P+KPTO9pjIngA6+ova5UjwMrThSL0DwKtbapTsdG4GjtF2b58mUuFa2TSOJZE+TYmQbBWnS3WmWVVXmHR0aWR7tPTHm+0xpT/1aUB/6jlodjz9n18ruczoL5Z5kWkjK/EqEr+WNVKobmhD4zHzVfrNisswdYKBZsUbMBLRWMmg1ama99msqCOtFRvKPnXO0/zSY9XrolVboXO8qteFQX+qdoB636kTpKn2f371uRT09CeZFyuhrjaqb+1M+Jky2qlv7aA562F/fZtxnwBwtNe5qvqkpHmPtqj5bGjrTGgXShv1Pk+vto7j9NGDsk+CbUPbcVsOLei0Lc3n8eOlPDyOtvVu2+D119zelZDt0bkuFP2ptq+O/hKvi3c71XsFSLymJJ8t5vu/oxr9n+rx6ugv1rYjc4cp+fYTCV4CJNbRqLzAKnbzUm2P9Aha5XXlOjev6sNO5yF5pE9H7u0u7G1HQfV4Wzu6ld8tc0Tjuqja6lwXVVut80wBn171V68Y9ADAEUXbts5u5ZcOclwXrb6IQfNHWtqN++wbyAavv1jbFuUOV/4bfiHBiwJeI9XYTT8qjx6lxmyL8mgiiQk6Ji5SJ6VzI7k6DFrHqH8j6dpRRjuqHbnOA0C1g2vr6Fbe1My6gMknLZiw9es+I/lk0t8xxfdGcehPqy9yZxsN+VQ+T41sT8x2FPG5FAQSvARIPErNywLgXdRd3ZF4uj8mEu8BU2ePT7q4+jzUDUwz9OmMCT45Rtuq01x6ozq169L3WL2jqoXO7giajrlXe3hH9Tx722oF7Ab0pxU8q+rPx4ekEf0pTo91dkcSav5M+AQ0AnaN/k91as2PgF3l+eI3ErwESKyjGU280LF5xbQ0YGQOLWsTz7woiOvXf99LtgGA7ojTp3bAK36OQr3e93qBhOJ59nkAEOa2WaZE/KlhIvlUPN7eDywKfmZevF5RreDZr9F2qutP0VanbkrVp3uwSfbp56BG0U7FpwQvluJ1tFPbGJ0HHT2CdqEPNh4DAIzKHY7MDFqVcMy2eAQtY1PX3I46RVHXNbcrF8DGjpfKsc5u5ZFHdUMvn4TjPdSkNrcd04Iftl71p3qsAHCwdxt5ROs8FY9X5zxVferor49PQtFtTa/7xeu9Vqt4nwHAQcVrWtuk41PNVkvzij4PN6tnbGoU77O2jm7lgElVf8Dx9pXgxTLSiMuN9h1pAQBMLMoF4H20U3WkFQBQ0WMHeH9gVfbYThid5/k4AWBfj11meho+N3ucks/sYXQ5xc41tmOj1/vok/qo3QiFbapjxzuKWDHf2tGFw83RmzdWEO21o6rsdZ5ecRxH2XZfXUuCHWW0HdNufna0fb1qIXY9c4bFfHp2GT/PuK1Hn/vqErVA8RmzzTWov6r6xPP0SkNbZ3xKjuo3dl1yiRqKRJwk96jqdfF+YWJaKCDqL2aXRbw/gWR9Ec0ntW3dPuNt5NWuPrF9KPS+z7zS3N4VH+C6n01cSPASEJ3dERw4Go2qqYFETFwTFQRSpWhb2dPRzJk4iriPMPDcu5U9PmnnCagHW5XJAjyPd36sjSZQ26jHrjBnGAqyaYHPPsXzrG1qR3tXBBnpaSgfmUOyrTzSpuTzWGd3fLQ9cbRaG1HtAJd2ibaq1xMAKnuCNFUtqOgvdq9RfcaCgbH5WcgmPnj2KbbRoWYd/an5TNQfTbuqGgLU71HV80zwqdhfTxit3v+R7+2644M+av8XBBK8KOBFJPvqWtEdcZA9LB1j8rNIf/+j2p6MDVGYR1s74inM46L2puiPD/f16YXO7ghe2rgfAFBRROvcOrsj8Q55AtF2d20zAGAS8XgjEQe7YrbF0U7K62hnz6HjbUTd8+ejHp/UoPKjQ1G78pHZyMygjSZjtlSfsbbNz87EqNiqNY8+P+5pI+rI7GhrR3w+nWobO89JxXTtxvQX073X7ICq/gBgT8+9FrtHvervuM+8uP6o14X6kIzpVk1/an1K7DwLsjPJ+4nE+jHqedY1t8frZSpG0fqi+H1GPM+OrohysB/rw6jBHQB83HO8sTby3P/FdUv3GQRGgpef//znmDx5MrKzszFnzhy89dZbA35+1apVmDNnDrKzszFlyhQ88cQTJg5zUCjPqw2V9QCAU8eNRDrxQbc+Zjt+JGmqKuZzypg8cmS8scd21rjC+O+8dFJbDzTGv//WJSeRfP5u3Sdo74qgMGdY/Cb0+vDYsO8ogGgbUfj4cDOOtnYie1g6ppcWkGw3VUV9nlJ+3M7L4ba0d2HHwSYAwKkVhYN8OpEN+1xaINhFIk78mlJ9xuxOr6DpDziuwdPGx3x6vJ4x7RbnIZ84/bPBdb9Q2F7dhPauCAqyM8kPAVX9HWw8hv1H25CelqgjTz5jbUu8nm0d3dhW3dhjO1LNp47+iG0U83n6hFHx33nVwvF+LObTq/6OAgBOHDuCPP0Yu0epbbu1uhEdXRGMyh0WD9i99n/re3zOJvqsaTiGAw3HkJ4GnEzUX6xtZxLtgiLw4OWFF17AXXfdhfvvvx8bN27E+eefj4ULF6KysjLp5/fs2YMrrrgC559/PjZu3Ijvfe97uOOOO7B8+fKgD9VXVu04BACYM2kU6YWOVUda4yOPMyaMJPl8Y3vU59yJowb5ZCIt7V3xm2HOxCJSlPb69loAwKXTx2I8YcTSHXHw6IodAIAzJ40iZTIOHG3Dqx9G398012XrZY76Lx8eBACcMWEUhhGLod/efRhATxsReGvXYXRHHIwbmYPywmgbee0YV/bo6EziNd1YdRSNx7qQOzwDM8ponc1ft0Wv6ZmTaOe5/2hbfER4BvF44/fLxFEk/dU1t8eDyjMn0Xy+ueu4T4r+jnV24+8fHY77pOhv1c6oz+mlBaS6lUjEiWuBqr81Hx9GV8RBSUEWykdmR4/Vo/5ix0vtUzZU1qPpWBfyhmdgRlk+yTauv4mjSMHzwcZj8UHCHOLxrtwR9TmXqL+G1k5s7NHfnAk0n2/09J1zJhaR9NfS3oU1PfqbO6mIpL+Y5qn6cxwHb+2K+aSdZ1AEHrw8+uij+NrXvoavf/3rmDFjBpYuXYqKigo8/vjjST//xBNPYMKECVi6dClmzJiBr3/96/jqV7+KRx55JOhD9cxgEtlYWY+/9Dxcrzq1zPPf7eyO4CevRR/o551YjNEjjk83DdbZbD3QiOUbPgEALDqt3LOd4zh48s2P0drRjSnFeQkdzWDnua+uBc+s2QsA+IfTyxP+b7ARxAtrq+JTXPddMWMQT8c52tqBL//yXQDAuJE5OMNjh+E4Dv74/oF4+362pyh5sGN1HAfbaxpxw1Pv4sMDjRiWkYZLpo/13L81tHbisTd2AwCuOs27FgDgb9sOYt2+eqSnAQtnHbcd7Loc6+zG0r/uBABcPrMUw9K9pfsdx8GKrQfjHdw/EHQUiTj42eu7AADnTCnC6DzvU6X76lrw4obo1OOVhPslGgDvRGe3g1njCjHJlT0ZTH+1Tcfi2r3C3baDnGd3xMF//G0Xmo51kfQHRAsen3prD4DE8xzsWCMRB798+2PsP9qGguxMXDRtTFx/gz2wjnV24/GVH0V9zionDaTe2F6LtXvp+mvt6MKjK6L6WzirDJkE/b26pRpv9gRM7j7Fi+0Tqz6C40QDSspU/faaxrj+EvrOQewiEQf/+foudHRFML00HyeOHZFwPANRWdeKZ9/Z1+PTu/6OdXbj4T9vR0tHNyYX57kynIOjqj8gOojacbAJwzPTccm0Es8+gyTQV0N2dHRg/fr1+O53v5vw+wULFmD16tVJbdasWYMFCxYk/O6yyy7DU089hc7OTgwbljgd0t7ejvb240u/GhsbERTJouMjLR1YuaMWO2qasOdwC/bVtWJnbRMcB/jUjBKcUl6It3si1sTlxA7W7q3Hhsp67K5txr66FuyqjU5ppKUBd1wanYJJ1tXUNbdjxdaD2F7ThL11Laisa43P9Z41uQjnnlCM3T3zmr3Ze7gFa/ce6fEZPdbYfPhtF5+ItLS0pB1cfUsH/rylBttrGlF5pBUHekbZjgPMKCvAFbPKEjYnc1Pd0IZ3Pq7D/vo2fFLfhm3VjXj/kwYAwL2XT8MJY9w3/XG7jq4I3tx5CFurG3vatgXba5rQ2tGNzPQ0/OQLpyGjnzm5rQca8dauQ9hd24wDDW3Ye7gV+49GC1hPrxiJz54+Lv7wclPbeAy/W/8JtuxvwJ7DLag80prwErMvzK1IWCbofni0d3Vjw76jqDzSgqojbdhd24y3dh1CS0c3CrIz8dVzJ2PNR3V9fHZHHLy9+zA+qDqKvXWtqGlsw/76Nuztqcf44lkTUFKQnVR/h5vbsWrHIew82IRP6tvwydE2fHyoGU3HujA8Ix23XXRC0qWujuNg/b6o/j6qbcG+Iy3Yc7gl/tlr51ZgUnFeUv1VHWnFq1tqsL2mCdUN0TatbjgWnxO//eLj04e9+8V9dS1Yt7ceHx1qRlV9G6qOtGLrgUZ0dEdwWsVInH/SGPzx/QM9bXucjq4I/v7RYWyvbsK+uuh12VHTFF/9sPjTU/sdpW/Z34C/7z6MmsZjaD7WhX1HWrH5kwa0dXZj0uhcLDqtHEv/uquPXVd3BO98fAQf7D+KHTVN2FvXij2HmuN7yiz+9FSk96O/LfsbsHbvEUScaKZwz+EWbKysR31rJ4ryhuOLZ03Aiz2DDTeO42BDZT1W767DztpmVNZFr0vM5+2XnNhvsW51Qxt2HWzG7tpm7O/x+X7VUdS1dCBnWAa+cu6k+JRMMtu3dh7GzoPR89xb1xLPAF83gP66uiNYtfMQPvikIXpNezLHLR3dyMqM6q/P1gSuNnpt60HsPdyC/UfbsPdwS/x6fvGsCZg4Orn+gKgGN1UdxUeHmlFZ14qt1Y3YXhPNutx64Qmu9ky0q2/pwOvba7HjYLS/3nu4BbsPRfuxsyYXYf4Jo+OBTG8+PNCAdXuj/fXeuhbsOtgcX7p+56Un9au/pmOd+Nu2WmyracTenmfErtpmdEccTCvJx8KZZdhW3dTHznEcrNx5CO98VIedB5uwr64VlUda4+8luveyaT39dd9AqzviYN3eI9hW3YiPD7dgb10rPvjkKI4Ooj8guqptw7567DzYFNdRrN/60lkTUJgC7zUCAg5eDh8+jO7ubpSUJEZqJSUlqKmpSWpTU1OT9PNdXV04fPgwysoSR2ZLlizBgw8+6O+Be+T/vleJB//0YdItsC8/pRT/8r9OBdA36Klv6cA/Pbse7+050seueEQW/s9VM3DW5OSp4T9s3I/7XtycdNv3T59cgn+5+tSkHWpHVwT3vbg5np1xMzwzHbdffCKuPmNcn/8DgL98WIO7X9iU9E2k5544Gj+55nQMS/IOpkjPyPjnK3ej91vfM9LTcMM5E/FPF0Q7mt4B0+ZPGnDrs+vjAYebk8aOwL9fezpm9tTnxG9eJ3rT3v/SZjy/tirpeX7l3Em4+1NTk7776e1dh3HLM+v6tG2sOa+YVYbv9ZMlWrH1IL79+/eTbpI1tWQEHr761IQHQCzoOXC0DV9dtjbe8fb2e82cCvyfq05O6vO5d/fh//ufrUn1V16YjYc+Pwsnjs3vE7wcbe3AP/4muf6GZ6bj2rkV+P5Vyc/zN2v24oE/bUV37wuK6HLP7y6cjvNOKsbOg4nn09EVwf0vbcbvN3ySdHR55qRR+M8vnoGM9L7B87q9R3DH/92IA0kegoU5w3D/FTNw8fSxfTbf6uiK4H//7n38qScY6s3JZQX4zy/NThoM7KhpwjeeXR8fFPT2efenTsLne+4Xt/4iEQc/+p+tWLZ6b1Kf40fl4D+um510n4z9R9tw27Pr44G9mxFZmfjGRSfglvOn9Pk/x3Hw0Mvb8Mu39yT1WVqQjUeuOQ0VRbnx4CWmv8PN7fjO7z/Aqp2HEl7WB0T19/kzxuMH/ehv9UeHcd+Lm+NFz24qinLw48/NwpQxI5IGL/++Yif+4299A8bsYem47swJ/d5nbR3duHf5B0mv6fDMdNyzYCounVHSR39AtO/83kubk/ZjF08bg0euOS1pANLa0YW7nt+E17Ye7PN/+VmZWLxgKhbOKku6+dsb22tx1wubku7Jcs6UIvzkC6cn7Yta2rvwtafX4p2P+96j5YXZuOeyaQnZMADxCGbnwSbc9tyGePDppqKof/0BwC/f+hj/vmInWpK00RWzSvHdhdOT2nEQaPASo7cgHMcZcC4z2eeT/R4A7rvvPixevDj+c2NjIyoqKnQOd1Acx8GmqqP43kub4TjA9NJ8nDW5CCeMGYGJo3NxUkk+xg2wrPD7f9iC9/YcQVZmOj51cgmml+RjYnEeJhblYnpZPrIy+3amDoDdtU349u/fR2e3g5PLCnDuiaMxZcwITCzKxZQxI1BamJ3UDgAeXbEzHricNakIM8ryMak4DxNH5+K08SMTpqiOn2d0RBYLXKaX5uPCqWMwuTgP5SNzcFLJCJQVJj9PxwGWrd6Ln/VMmZw2vhBTS/JRNjIHJ4zJw9xJRUnbyEH0xv3H36xDdcMxjMnPwgUnjcEJY/MwsSh6vCeXFfQ74n3u3X14fm0V0tKAS6eXYNa4QlQU5aCsMAenji9EXpJ53pjPu17YhLbObpw2vhCLTivHCWOjbTt+VG6fDsYtxQ2V9fin36xDxAHG5GfhlPICjBuZgyljRmDWuMKeuoi+x+s4Du58fiO21zQhPzsTn5pRghPHjkBZYTbKCqPtNLYgyTV1osWQ3//Dlnjm66xJozBhdB7GjczB+FE5mFaa3yegjD2wvvfS5rj+Lp0xFlNL8jFxdC4mFOViRlkBcocnb6MPDzTgh3/8EBEHOHtyEc49sRjlI3PQ0h6dQjlzchEKe+0I7dbf79ZH9XfmpFGYXlqACUW5qCjKweTiEZhaMiJpGx1p6cAtz6xDfWsnxuRn4Zwpo3HCmDxMKMrFxNG5mDVuZNLO33GAx97YjT+9fwAZ6dGpvopRuRiZOwwTR+diSvEIzBxX0LevQTQ1//Vn1qLqSBtG5g7D+SeNwYyyfEwpzsOEojxMGZPXb/bj2Xf3xQOXsyYXoSh3OMaNysGk4jycNHYE5k4cFV+14/bZHXHigUvOsAxcMn0sTh1fiImjo+fa26d7tdH/fFAdD1yK8objrElFmDA6FxVFuThxzAicMXFk8j7FcXD3C5vitQynV4zEGRNGYXJxLiaOzsO00nyU9KO/TVVHcf0v34XjRN8w/OkZJTipZAQmjo7eoyeNze+TFY3p7//8YQt+0zNl8umTS3DmpFEYNzKqhakl+UnbNqajB/74If70/gGkpUWLiKe6fJ49eXS/00XuvnNaST7OmVIUvb9H5+GEMXkYP6rvap9YkP3wn7fjta0HkZGehgtOKsb0sgJMHp2HCaNzcUp5AfKTLI5wnGgNzu2/3RCf4jnvxGKcMCYPE4vzMKU4L2mReOw8H3ltB975+AhyhmXgs7PHYea4AkwuzsOk0XkoLcjut/9r6+jG156Oajc/OxPnTBmNE8eOwOTReZhUnIczJoxMqj8AeGVzNf755W0A0HNvFWJCUS7KR+bg7MlFOKmEVrsUNIEGL8XFxcjIyOiTZamtre2TXYlRWlqa9POZmZkYPXp0n89nZWUhK4u2FFkVdz/35JvR+dVFp5Xjp9ed7qmwzIGDyrpWvLy5GmlpwO9vnY9Zg8xZuv/sr/++F53dDi6aNga/uunMfgUMJE43NR3rxNM9HepjXzpj0LoCt8/fvluJ1o5unDFhJH536/x+p2l6++yKOHh8VXSu/f4rZuCWC/qOGPvz+cf3D6C64RjGjczBq3edn7RzSEbEAX7x1scAgB9cdTK+cu7kQXwed/ry5mocbm5HRVEO/vvWeUk7+2Q4DvCz16OZpYUzS/HTL85OmoVKxobKo1i7tx7Zw9Lxyh3nD7pE2N1Gv3jrYzhOtC7lPwbRn/t/Kuta8crmGqSlAcu/MT+evfLic9nf98bP8+fXn+HZp1t/P/vSbFx1anlSm944DvD82krUt3ZiWkk+Xvrm/KSBVTKf7V2ReBDx6BdOw2dOT55VjNv20l/VkTaUFWbjlTvOH/QldGlpaYDjJOjv+1fOwNeTZEn62PXw9u7DeP+TBuRnZeKVOwfXgptf9vi845ITsXjBNM92Hx5oxFu7DmN4RjqWffVMzD+heJDjPf59rP+bXpqP33zt7AFrTNzXZe/hlnjg4uV43T7/a9VHeGFdNKP6m6+ejfNOGvh4geMP5qfejvadF08bg6cG6ztd/3W0tQPPvxf1+aubz8SFU8f0b9fr5+fe3YeWjm6cXjESv7t13oD9gttnc3tX3OfjXz4DF00b269d1DYtHmn9yaXdl+84f8CdcN36cxwHP+3JhN1y/mR874oZ5JWGpgm0YHf48OGYM2cOVqxYkfD7FStWYP78+Ult5s2b1+fzr732GubOndun3oWLjq4I/tZTEf9PF0wZ9CK7//e1rdHAbP4JowcNXNw4jhNPW37tvMkD3ny9WbnjENo6uzFlTB6umFXq2Q4AXutZnXPT/EkDBi69Wb+vHoea2lGUNxw3zZ9E8vnqlmgbffmciZ4Cl1j7b6tpRNWRNuQOz8B1Z04g+YytNvhfZ1R4DlwAoOlYV9z225dN8xy4OA7wWk9R98KZZaSHVWe3E9ffP3rQnxu3/gYLXBKP18Fft0W1cPP8SSSfq3b26K84D1f2TnUnwf2n/7w5erxfPW/SgIFLb97bewQNbZ0oLcj2HCwd91kNALj+7Amkt+e69Xf92RNJPl/5IOrz82eMI2nhUFM73v+kAelpwI2E+8xxoiNtAPj0KSWDBi5u3Pr7yRdOIxXHxvR32vhCUqDlOA5+9no0izujrGDQwMWtTsdx4gsobjl/CqnvfH17LTq6I/GsM4WYz6+cO8lzvwAAaz6qQ1tnNyaOziX5dAC8siV6TW+YN5G0hf9Hh6L1hMMz03H7Jf3X76QSga82Wrx4MX75y1/iV7/6FbZt24a7774blZWVuPXWWwFEp31uvPHG+OdvvfVW7Nu3D4sXL8a2bdvwq1/9Ck899RTuueeeoA/VM1sONKK9K4KivOHkvRrW7Y3OOZ93Iu1GqDrSikNN7RiWkUZawhorzASAC04aQxJl47FO7KyNzh1TOjcAWL8vOlc7b8ropGn9/uiOOPF5+fM9jKzcbIgv9x6FHMJW3Y5z/LrMO6Fvdi8ZsbqMjVX1iDjoSe2PGMQqkdh1OfdE2nlu2d+gpD/3eVL1t+dwC+pbO5GVmY7ZhBU2juPEfV4wlaa/to4ubO3Zn+T8k2jHu35vj/5OGE0KuiPu+4X4sFq/V0N/PfeLV58x/cXsppUWoDjJ1O9AHO8X1PV3MmEZvlt/V3gIYt18fLgFTe3RouXHrz+DZPvy5mocaenA8Mx0zCEs83UQXVQB0PuihrZO7DwYrTmh3N9uLcybMpp0v0R69fUknz33y+yKkX2mflOVwGterr32WtTV1eFHP/oRqqurMXPmTLzyyiuYODE6Mqmurk7Y82Xy5Ml45ZVXcPfdd+Oxxx5DeXk5fvrTn+Lqq68O+lAHJdZhfLg/WlBH3cTLcYAtB47bevMZ5cOezeBmlBV42hrcfVhbeo53tsd9Y2Km26sb4TjRJcleRldunx98ouaz6kgrmo51ISszHdNLaXOssTY6zeOmWDGf9a0d8ZeVzSJkIwDggyriefY0kuMcP97TPW46FtPfFqr+XB/5sFpNf7FNzmaUFXgKRnX0F2NbTRO6Iw7G5md52po+wWdcC17bNkrVkTY0HuvC8Mx0z/vjxGxjgRa1bVs7uuK70Xq1jbGVqiGX/o7ryNsD3Q/9bSVulhdv257znD1hZHxnbC84joNn1kSnqaaVJK8nHOBwsfVATLuDt1FakvOsKMrxFFS6fW7vWXlEbaPKumjfGd2Ac/C+0+1z837v55kqGCnYve2223Dbbbcl/b9ly5b1+d2FF16IDRs2BHxU6sS2dJ7s9SbqUUlXxMGBntUzU8aovTjRs88ktlOKaZkB8nkmsaWep/s9Ib0Ly/ojdhPGXkRIPd7Yuz5KCrJII2bgeNueSMy6HGpuR1tnN9LTgAnEd0KpXpeo/qIrP1JZf72vJ/VYAVcbEa+L+11HlFQ/cPx4J5HfxxPtE0bmDktaOD+gz57zpPo80tKBlo5upKXRX6ugqr/o+956+j/Fe3Syx/N0BxJ7e1aM3fUp2g7ggLruYyuwJhP7XED9/WDxV10U5XnuO3v7PEHhXuPCSPAy1FB9mdb++jZEnOgbT8cS33dUGRcmzWdrRzfau6LL9Px4+ZxX9qm+fE7rxX6xFxEq+iQEEbHOsbIu8T01Xom1bVlhDmlaDdDQ39E2dEccZGWmYwzxIan6MkEd/alqCNB/+aHKS1GVHzoKL4c8rj89n2UF2aQar6ituv4iDpTe93b8pYk0n22d3ajv2b6AsqEgEG3b2NYH9H4sdk1p70mKOE78jeWqPqn9tQO9e40LeTEjgViH0dEV3VeD+iLCju6YXa7n6abY59y23ojatfcca352pufNhWKH1tnt9Pj0dp7u/Tlie6UkW4I4kM9Y23q1cxNrI683YO+X21FebxAjth+C5/Ps+Td2nioPLKr+YtfleNvmeC5ajOkvtgeISf21E9vIrb+O7gjS0jDglgUD+aQ8APreLzT9xfbpqVDQfOy6UPUX8zneoP5i7TN+FL3/i+0r5L2Nonax88zP8v6Sx9ihxfaKyh2ekXSLhb4ej5/TJ7GBFFELBxuPobPbQUZ6Wr/bUPRnW6XoMxJx4vtpqQxUuZDgRYOx+X33QEhG7023qKPeBFviiCUGNdOT4FPxeAtzhnmqz0nqk3C8vfvB0YQqezfFGm2UbI8dL6heTwAYM0LNp1fdJvXJoD/V4x2ZM4yc1Yqhc11U9aelBUXbZPu4eEX1urD0fwVZyitoqMW6AFDXEq2hox5vbIfhorzhpEJzILrhIBA9VwpH2zriwaGOBk0jwYsGlGWUbihL2PyyDYvPguxM8nxvjFG53n327laKCLYJdhptNCpPbVWA1nVhOU81W1U7QP14C3OGKeuPEvT01p/XrEJvijS2eg9L//eFufRNT2O77VL6FLedyn2malvX8465guxMcp0XJ/YcaQrg1wOL8tDp7dPrzdB7kEG5eXtnijx3Ur18jiJ0jH18GnhI9rmeisFA9rB0z4W+va+LVsDk8Vz7+NTRn0dbvfPsrQWPx9tb81o+1Wx19KcaDIzIyvRct9JXC8EH7P7qT9GnhhZGerV1mR1pidbKeL5He4zrY0EPqY16bFtjtl59IsFOJ8DjQIIXRXQeWKpBD8CUBWEYbasGWzqjbc+dVC/0rqdawJSVmY4cxSk52/RnU+ZFNQOi41M1A6fjMyvTe//Xx6eO/lK8L4oRDySomZceO8pb2eO28Sknmh5iBc069wsHErwoonMDqopkeGY6chU7DNVMBqDx8GB4qOv5JARMaf4ETDqpd9U5fNUgTSdg0gokLAoqtfoFSnbAde1N+XTDEYxSBox++QTU2jdWQ0L1G1s8oBKQxgqiqdc0VoCtoyMOJHih4OowCjR2IaTsYOh+PhXmDPNepd/rZ8qI0O0iLS1aqU+10/EJAIU5iiNfUtsmOlXdWZJ0ntDw6TpekoZ6/ax6XQp09Keo+ZhfFTuvq5uS2lL01+se9e6z9/SEmv4ofVHfKRHz+lPt/wo8vuMsmU9V/aWlRVfKUe1i5Clm5ynn2vtkvQ5O+vR/GllDDiR4UWSExwc60PdG8rLsLhleb4RkUN4Lk+gzk/QukARbxfMEaO3rbuDcLJ02UrPNU2xbQEMLGm2raku6Jj75zB2eQV51EfepdV0UtaDTRorHq3NdVPsFDi1waD53WIZy/3fjvInK2VG9c1XTrs415UCCF0VU05eAeiel2tEA6oLWOU/VYCAtLZoiVkHngaUc4Gk9PNTaSNVOx1Z1ykjHp47mtfRH3LhN1yegHnhr9UUcPpX7Pwb9adzbKvv2xP0qHm9GehqGK64Y0rnXOJDghYA7hqbcgL2jb9qNf9yWIujePinCdJtSsj29xxgkn267YRnKIxZK2/Z2ofxQ1/BJ0pGqnYYW3F4pDzq/9EfSfK+fKcGA2zaHONpO0K7ieUZtg8/89fE5zLz+bOr/OPQHqPcLucO99519+yL14JADCV4U0RkFqI52dEYBNo18c4g+VTtVN+lp0YJUFXRuetURLIf+qNfFD58so23bsj06U6UM10VdCxz6M68FQP266GSedbJpHEjwoohWh0EY7STama950Xt4qE7D8Nz0lGxP4ujMfJ0Dh/70aq44ghd79EfO9rj1p3g9AZ3pFIb+j0F/HHWGUVuOoFKmjYYsqg+sPqlaUvpd0a7Xz7QO2ZWqpRQm90nVEnymudPDZm76hKkCQx2j2ye1todbf6TpsV4/q09VmdcftbbHj2kGU1mFPvojZHu49ac6DU21dVvTpoT90V/UVnGaS+MZoRP4cCDBiyJaES5HFoRhtG1qqsCPLIhW0a3GSgbV2h6eUahdI0meAmx7MpzUbE+iT7v6P1WfHAsAorYc01ySeQkF1q2kYVhVkGNRelhrtK04nadTQ2LdShqGFXYsqXdDD2bVbE+iTyb9KfZ/OtdFubaMYeuFqC1HgCeZlyGL8qoCtx1xtK06taG1qsWnanvVdDb1AZCWsCKBsj+M2kqG3tBW4ajZAf7pT3UljZb+lKcKeFZsUEhTnGbwS3+q14V+nm5by/o/i/QHUPtO1zSrYv8HSMFuaMhWHG2r2unaqq6k0dnbQ72N1GWpastxXVQzIFo+Nc5TSwuK56p3XTi0YN6n6nXhubct05/WParej5XkZyvZZWn4lOAlJAxXDAZU7QAobz40LCNNeW6b43h1fKoGaVSf7tEOixYYfA5TvJ6Z6Uz6U20jxfMEePSnel1Ef4PD0UZfmDte+XhV9Qfo6Z4Du46WGXe6c3iG2rJGcielaKvTubnPjGLbO1VLOl6XLfl43W2UoZbm1blxabbqQY8f+tNqW0X90R/MLp8M+tO6Rzn0Rzped79Ae0Am6o+hbTn0p6ghHVudNiL1171+Vg0OubDraFMI1ahaRyAcmQzrRmeZaiOWYQznSX14+OGT43pqaZ5hpM6hPxuyjYm2ivcZgxZs0x/HMwLQy9pwYNfRphCUEVainV1pcB3bTNVULdFn7DXyKraqdjojpeN26nPMLPrjCJ5Ff0nRyabFfVqmP5YgjUN/lg02ubDraJlJnE6hCFMnhalmyz1VEP07lFUF6m3U2R1RsvUjAAHUp9bID6wEn/boz1RgGFb9qV8X4rSRok+/9GfsPmPWX5ZGv2BqmisVsOtoUwiOLIjqCIujiCtDcdSh4rMrcnzoqzwK1ZnCUR75Mkwb2ZZ5UbTV0R9VQ77oz7bpWcX7xTr9MfR/OtNGXLYc2HW0KYTUHAzmU+PmtWCFk/oo1GXHMYWoVc9hT8DO0bY6fvUeOgyBhEX9X5j0pzpQTUtTn+biQoIXAqqpSP9WMph6MDNPVWnc+LQN2NTT2d2uQgeSFlzf62iB0klx649axOqP/qgraY5/b4P+urqP6y9LdYWTZfpTLoC1TX86qxBVV2RlpCu/qoQLCV4UCcsKE5apKoa5V+p5dna5gheGjJj6SgaZqhoMncyfKtTzVK2zccOhP+sWLDAU49t2Xbiw74hTBJ4VJvLACgrqeXa4Hh521TmY36WZY0mu1qoW6/Rnz6oWncDQptVGWjV0FkybpwL2HXGKoLoZkLHK9wQ79Ze+qXbk1E7KjxoSKjpTBZ1dag8P1RQvoL4k1y/9qa72MLXpoBsd/anu1ULFL/1lkq6L6I9iaypI02ojty3Dajcu7DtiRjq69acKeFaYmC+eta3qXWfkqzpXTH5guVe1WKQ/jilEvY2+zL/jhb482xn8Q4P5ZJiq4sj2cOiPq2CXY2M8Luw7Ykbcox1Tm1G5YUkPWzS3rVNvRj1ed/Biyie3/kKzJJdBf9SHB4f+3IM3li0JLJo2sm2XXMm8DHESRh6GUrVuaCnBxEpyCn6kIclpUz92rDV0nkBiCt2Uzw5m/alOj+m8W4tDf6YCJpapUg2f3P0fh/6Up6p09EdtI+VN9Vx2knkZ2rhvXtW3fpI3wFJcEqnj0w3LVJVFo20ddB4ept6S64f+OEaSOkW3tk0VqELO9nSFT3/KU1VMmReOQmou7DtiRtw3ryo6DyyOl75xpE1VH1g6NyDHnC99eba+/sg+fdAfx3XJ1AieOfSnE+wr+9Toi1SxTX+q/Z+O/qh1aU3HuuLfK29tIZmXoU234lyBW4rUXQxdNZqkLafdn6RuVe1OuapuuqXlU7F4gO7TZWtogya3G6oWOPTX7YP+qD7dWqD5tFh/6YamjTT0Fwmj/ihTVT7pT0cLU0tGEHwe/17ndQZcSPCS4uivKUi8kU3aKvtk2OnRhs0lVets/KL3i+e822n4ZNECg0/jHumEUX82MSY/C7nDM7kPwxgSvBDw4+bV6RiVb17LHgBh6Ww4AqawaME2/fFogcOn6C8onzrYMHjrjQQvBFRjFx1hOK6IifJ3/BJjqvtMsCN//riFqZvX7YcajIr+Us9ngh3586K/QX2K/rzbkT9vXn9+EmjwUl9fjxtuuAGFhYUoLCzEDTfcgKNHjw5oc/PNNyMtLS3h65xzzgnyMI2iIxL1UYCOT4apAuUMk3mfWhBd+pP507Bl8Sn6CwzRnwefor9UJdAJsi996Uv45JNP8OqrrwIA/vEf/xE33HAD/vSnPw1od/nll+PXv/51/Ofhw4cHeZie8eXm1RCJaipSz6dZO11bdaehcCn6C9hW3WkoXIr+Ara1yacugQUv27Ztw6uvvop33nkHZ599NgDgF7/4BebNm4cdO3Zg2rRp/dpmZWWhtLQ0qEPTQLXaXj095w6YKKY6PhNTyxSf6vhx85A7KcXz1CHxulC9iv4G9qmO6M8Lor+BfarDrT8bCWzaaM2aNSgsLIwHLgBwzjnnoLCwEKtXrx7QduXKlRg7diymTp2KW265BbW1tf1+tr29HY2NjQlfqQy9u1Cb89Xx6Ycxx42sN9oxfyeTteBL5o/oU/QXuF3UVvSX1KfoL3C7qK19kUxgwUtNTQ3Gjh3b5/djx45FTU1Nv3YLFy7Ec889h9dffx0/+clPsHbtWlxyySVob29P+vklS5bEa2oKCwtRUVHh2zn0hn3OVzVtyjHny7CqSgeWFLpl6WEO/aki+vPgU/Q3uK30fykLOXh54IEH+hTU9v5at24dgOQX0HGcAS/stddeiyuvvBIzZ87EokWL8Oc//xk7d+7Eyy+/nPTz9913HxoaGuJfVVVV1FMKHO5OgnojqaZc/bphzVXbu743l7c37lNvtYd5nwlaIFxV7geGjv6MoaE/P1YbkVc4cetP+r+UhVzzcvvtt+O6664b8DOTJk3CBx98gIMHD/b5v0OHDqGkpMSzv7KyMkycOBG7du1K+v9ZWVnIysry/Pe4MXXzCjRsqLb3Z8NChlGdTqduWXbAJp9hwbYpEcsOlw1y8FJcXIzi4uJBPzdv3jw0NDTgvffew1lnnQUAePfdd9HQ0ID58+d79ldXV4eqqiqUlZVRD9V3fHl4cHSMDLbcxWs2wBNIGHcp+vPik2WqQPQXlK1tXZhtxwsEWPMyY8YMXH755bjlllvwzjvv4J133sEtt9yCq666KmGl0fTp0/HSSy8BAJqbm3HPPfdgzZo12Lt3L1auXIlFixahuLgYn/vc54I6VM84DGkQ5bRpvz948JlQJEdI29PcJNpqpJaP/w3i9Jjr88amcBL802xFfwSfRER/gyP6I/gkwq8/+8KXQDepe+655zBr1iwsWLAACxYswKmnnorf/OY3CZ/ZsWMHGhoaAAAZGRnYvHkzPvOZz2Dq1Km46aabMHXqVKxZswb5+flBHqoxqCJxfMj36ExVKY88tArz1G2VfZp3ScafzJ/oL0hbZZ/mXZIR/RF8iv4CJ9BN6oqKivDss88O+Bl3NJ+Tk4O//OUvQR6SFtwvSVT2yXIj6RSvMcBS5xCOVQWiv9R0KvoL0tYu/VmYeJF3G1Hwo/JdZ5MmmtOk3yr8GcpqDx0/+n+HvBlVwvdm7t6EVC3RVvQ3yGdFf4P7FP0R/0w49Gdj7kWCF8OQOww/fOrMbRu6kXpZG/dpwwoTP2oORH+erI37FP3141Pbo+gvaFsuJHgh4MvNy5Kq1Uhhqt68yh7DM+dr27JlZZ8sK2k0bEV/wfm0TH/S/6UuErwYgHsZpM5mVBRTv/bzUC6SI1rqbOSnih/pYbpPDUfKqz38OTljm4T5cLg6+jOFjv7UZ3Ds0p/6Znw293+KThmR4MUwZJHIJnVGsOHmtXafIUPBs1/wFExaIEAfsEF/CbYcPkOiBV0keCFg7WojHVvFG0lu3sGxTQscPsOiP+4sjPjsxzY0fZF95ynBCwU/3rNhaJ8DLZ8+7HNAJfF9Nop/g1pt79OKBGWfhrJwor/BEf15QPQXGOz6sy92keDFNOT+wtrX0EuRnAevpE9zZP5EfxSf6oj++vEp+iP4VEeClyFOWGoOEmyZi41NejXu0TItcPhU3yrdvE89RH+p6DMs+pNpoyGO6lLpxPSwmdFOQtaU6tOdNlUdBTCkMMnv9tBYkaWKTnpY9Bccor/BEf0FB7f+LIxdJHgxDcvGUFSf2h7tK0BkGWuT5//98En7vOiP4JPJVtmn6C+5T22Poj8TSPBCwNZVy7bNhfL45EjVCkHBUXOgQ1j0J5ofnLDoTxcJXgiov55dv5Kc7FNDjP4Uyek8PMxUrHGv9jCFbfpL/Du+/BkzaOjPFDrTE370f1Q4NoxL/DuKdtL/BY4EL4bhSGHqYGMhlwrcDxIv+BJUcqTQiZ9PXKrKUcgajiwIBzZMp7Drz7hHO5HgxTA8nTGPrbJP8y5Z4M7CGPPJoT+t1R7msSF49sVnWPRnGTaepwQvBPzZMInoUzlVq4NrnwPVP0FOYeqvvNArWDNz92qtghD9BYbob3BEf8EdAL/+7EOCF8OQU5h++NTZYVJ1e3Ylq5itok+N4QPLyINj2oj4edEfwU7057dL0R/FTkt/9oUvErwQ8KfmwC6RhCWFzgFP29rVuMqzRrYFEqZgLtIMi/60fBpympDtMePSVyR4IeDPhklUpz6kajVcmtswKfn3qn/D0+d9SNWSSVjtQU3ba7u0Tn+mCI3++vHvhTDqbyj3f/4Z8yDBS4rDsjFUwrs9zKdNVY3t25KbCMNDnV9/ij7VzLSMh7z+/PBJ1gKHT+n/bECCFwoMS1X9wAZh6mwh7rd/8Zk68CxbHrr603olAfe+DUMYjv4vIdtjWb8ASPBCwo9qe45OigrHfiJ+2FLPmePm5X4wDuXz1CGMWwJwrDayYWM8P7Ch/0u0tQ8JXgyj85KwsGDjjaSCqYeHlk+G3ZZZNuPzo85B5+EREtHb8MZk7j43JFLQRoIXw8hUQXL4O4yhO1UQSp9Mtuo+DU0VMGTh3C0aFv1RYSlSZ9GCf0jwQkBdYBpTOKqjbZ1qe3VTZfdun8pFckQzlumUhLltGqK/4AiL/hL8Ez/P8YC1VX829H8JthbmeyR4MYzOhknKPomf597Phmfkm/pOuTfsUvZJ/DyH/hyXU9FfgC5Ff4Pb6run+7QvdpHghQJHRO4HNggz4eERktSyLI8N0KeO7RDWX0KROrnmStGnRrbRD6T/S47WyrMUQIIXAo4PGyYZq7bXSpv6URiqY6uaNqWuNkrsyk2gs8JJ9GeGoay/BP+G9Kfn04yfBJ8h6f96e7UNCV5SHH+yPRzzE+ZtrRttEz/PXQuiiujPg60F+vPHZ+qvNhL92YEELwRYOnJlP/YSmo48JBuw2YZN+qMXCacl/d4LfrwegEoId4qwYvFAKiDBCwE/Xs9uw3w6R5CWsM+Gqk+t1R6mlqry1rlw1FbYAPdSVQ5MTeHo+VScKlWy8g8b+j+/bLmQ4MUw9Mr38I09bHrzrMkOg2flBUOdg7ZHlWJU/dVGelqwR/M62FCwy93jhkULukjwYhhzKUG7bgD+Ik2zdrq2yj7tkoUytgUSpvTHMVXA/fZsG+Do/3RWnqUCErwYILEzTH2RcK8qMDZtpGGrCs/GeAyrqjROjiPbGBb9Jfhn6YrCoT8b+j+/bLmQ4MUwPPPM9owkAfUbSe/FZOFYnmiD/vh3OFX1Kfob1CNLkCb93+C29iHBi2GMpWo1bLmX5Nr0niEblidy5/1smx4byvrjLt6m11z569+TT32XZDj6P44FC34iwYsB9AIJ8xsmcawq8MNWZ7WRKbgDCSqq+tPSPPdLOoew/hL8c/gkZ+Hs1J/WtJEh/dmOBC+GYUlhsjykQ5Kqta0w1AL92bSfUqJPDduQ6I8D6f882FoohUCDl4ceegjz589Hbm4uRo4c6cnGcRw88MADKC8vR05ODi666CJ8+OGHQR6mUUxlQfS2xzaPHyMPHWybFyf50Vhhwp0FMQWP/tyrcBimCox4tE9/3JJn0Z8VeeBEAg1eOjo6cM011+Ab3/iGZ5t//dd/xaOPPoqf/exnWLt2LUpLS/HpT38aTU1NAR5psHCsMEnwT/w8z94e7n02zFcGcXTkHHCsqrIB7qX6qmgVaVqwwsmPl0FywNL/aZw0d3upEGjw8uCDD+Luu+/GrFmzPH3ecRwsXboU999/Pz7/+c9j5syZePrpp9Ha2orf/va3QR6qMWwUiWmU28iy7fZteHj44pP4eY7gOcGW4+GhbGkXNuiPG3lGeCOlal727NmDmpoaLFiwIP67rKwsXHjhhVi9enVSm/b2djQ2NiZ8pTLGXkNv2S3LvT27esCk4d/Ye67U08O+1J/YUGfDMG3Ek5E1P1Wgoz9f/FOLhJn3GTKFvNvIR2pqagAAJSUlCb8vKSmJ/19vlixZgsLCwvhXRUVF4MdJJeGGJXeqfoxCOR5YGraG7aK2DFNVFqxbVn+3jPrJcU/gWKU/hlV9WoREf9b1fxZGL+Tg5YEHHkBaWtqAX+vWrdM6qN4PW8dx+n0A33fffWhoaIh/VVVVafkOGp7liQxOtaZT1GxZlmeru2TBCv35Er2I/lIRK/Tnj1d1Sxb92aYkIJNqcPvtt+O6664b8DOTJk1SOpjS0lIA0QxMWVlZ/Pe1tbV9sjExsrKykJWVpeTPFBwrf7S0yFHh77gL1syQ1s/3pL9hwcjXOv0xY5P+yD4Zpqp0NkOT/i84OPTnJ+Tgpbi4GMXFxUEcCyZPnozS0lKsWLECs2fPBhBdsbRq1Sr8y7/8SyA+TUOOcFl2mGQumOQYhXJ05BbM/4dFfwkPD4v0Z91OwlSDkOjPD1sb+j8/CbTmpbKyEps2bUJlZSW6u7uxadMmbNq0Cc3NzfHPTJ8+HS+99BKAaMd611134cc//jFeeuklbNmyBTfffDNyc3PxpS99KchDNQZPBtOCTsptq5w21VntoejTttVG5l2yIPrzYBsaMUj/N6ithT0DOfNC4Qc/+AGefvrp+M+xbMobb7yBiy66CACwY8cONDQ0xD9z7733oq2tDbfddhvq6+tx9tln47XXXkN+fn6Qh2oMjt1GqbBsDMVSbc9bPGuDT+7iWVOERX8JUwWiv+Q+pf+zgkCDl2XLlmHZsmUDfqZ3NXlaWhoeeOABPPDAA8EdGCM2pGpZ3u2h71ILjnlmDlhGWOSlqgEdx0A+zbv0BfumjahasPPK8LwSg8eWi5RaKh0GbIxwVbAhSPPFp20PDwumqrhfDKrsk8nWJmy7LmHBxueSBC8G0Hu3h+I+B5aJkSdV6y7SVK05IH6eY6pA9DcotuqPim3688MnldBMG5l36SsSvJiGI51I9smwwyTHChO3rbZ3OjzF27SP26o/G7JafujPumyP6C84n5ZlgXWR4MUwNswz++KTJT1sfuWFVoW/ZW1klU/L2lb0N8R8Wta2tmVKAQlejKDTwaimE7W2xw5J2jThfTbKS1V1Pm9+qoBKaPTHkW30QX9UuPWn884pkk/b9CfvNiIjwYthWFKYllW+m3Oa8PQwjh1LVcOhP18eHgz609tPRMPYEGHRny8MdTH0QoIXw4SltoLlTcI6tiFZBWFfF6VGaPSnVedgntD4tE1/vh2FOSR4MYBbGOZWFajbhmWfDYch8aLznhdlnwn+qftsKPq0rDfk3k1E9JccDv2Fpf/jeM+Vn0jwYhhTc746PpneTGbepx9YdtPbkGFikQKDU27J82QkaJ/nbiNj+HCistpICBT7JKKGDfUcibZmCnb9srXKZ0hEb53+LAhkffEZEv3pYGMbSfBiAJ30nHq1vV2wDLZd3xtL27s8sWwSZoX+OPYZMg+3/kzBMT0h+hscG7MtbiR4MU7q73Ng61JBnfNk2WcjJLkXDv3ZsM+G6M+Ux9TXH0f/l2hrHxK8hAELOnJffDLM+epNFajbWoXob3DbkOiPRfOiv0BtuZDgxQA6UwXq75ZRV2NoauRcZ8oxhWMK0Z8Hn8wb44VGf0TbsOiPBYaVZ34iwYthyDcvh0+GyneWzah8SdVq2FpQMCn6C9BnCPVHRfQXnE/bkeDFMCybF1m2YZK6T405X5YUeji0wEFY9Kf3bqNw6E/6Pw+2FnYLErwYgOWV8Bq2YUnbu7G98n4gdPQXmtVuIdSfsb7IMv2Fpf9LS/jetjtWghfj2JZaVvZJTpuax5eUtNZqD/PYkGGyVQs2+HRjQ8Euxwsz/YDlfUpSsCsEiY0Rrgo86WHztrbVHAzlEX6CTwsezH7YaunPgmXL/vg0j23Ts3YdbRQJXgyQkJ4zNPKw7t0eLD5dqz0M+UxMoZvxKvpLTbj1Zwod/Sn7tEx/LNNGlgVYvZHgJcVh2TBJ3yUZ7vSw8vbsQzxVa6v+TL0A0C9Ef8mxVX9U+Ps/VvdKSPASArhHW8Z8stSf2FXzwoHoz4OtVZbq2Dadp+zTvEsLd1vWQ4IXEyRsBkQzVY3HtTZp4h6GGsKP7dmpJHQSDEvPRH/9+OR4n00I9UclLPrjgGM6z08keDEMy3s2ePbHpsGwqipxh1OGtL0FBZOivyBdiv4GQ/QXnEvbkeDFMOFJm9q12oiDsGiBg7Doz7rVbmHxaZn+bCzeleDFADrvluGA5z0v4YC7I6X7D8eVCcdZJmJu1sgu/YWl/+PYPNVPJHgxDD0it7NbNfUCQB04NgBMqHPQd0/GhhEh97J5cz6Nu/RlqkoHjmkjDjj6P9uycLpI8GIYG0WiAs+DmcGnZana0OjPsikcZZ9aNS/mCY3+uA+AiG3HC0jwYgSOd3vowL3yYkj7dI98DfkU/aUmtmYVqNg2PRGW/s/G5dFuJHhJccKyYRfLefrhUyMNbsMo1Fb9cQRpor/BkQ0L+/Hpw9+w7Q3jukjwYhhJZwfo07xL+5aqynUJzqd5l9bpjwPbrouyTx1bC6UgwYsBdHQRlkJCDtg3ozK2MZk6YdEfy5bwor9BCYv+OLBtOq83ErwYhyG1TE7V+lH5zjBVwJGSJn6eO20v+uvHp+jPEKK/pD4Z9OefMQ8SvBiGe58PYz4t6Bi5CUs6m4PwnKdGnYOPx+HZp/R/QXplsORDghcD2FYMJStMhha26Y+DsOiPY58XvfcM+XggqeyTQX/uq2JjFyHBi2F4VkHo/w2yT+LnWeaZWVaYcG8SRoPj3TL2PrAY5o2I8G+SSIO9LkgR2/p5G4u3JXgRhgy2jR4sO1xlbAieffFp2QXlWSodkikcy8Rg2eECkODFCNZV2/vxRyy4G7g3ozLVRKHUnwWEZaoq9XuCRMJxVWDfhelFoMHLQw89hPnz5yM3NxcjR470ZHPzzTcjLS0t4eucc84J8jCNwrORFRFL3y1jQ6rW7ZKnkDAk+gvLlJxt+iN+3lr9UV1yv9tI27t5Ag1eOjo6cM011+Ab3/gGye7yyy9HdXV1/OuVV14J6AjNIxnM4LBhx1C/bJV9mncp55miPqVfCNCneZeh00JmkH/8wQcfBAAsW7aMZJeVlYXS0tIAjogHLWFYWvluw70Qlo3Jwqg/ss+QTI+J/ry4DMfEkTvzJgW7PrFy5UqMHTsWU6dOxS233ILa2tp+P9ve3o7GxsaEr1SGYyMr8sZQ1q5wYjhPjU23bBiph0V/3O8ZMubTbSv6S+4zJPrrZWwdKRe8LFy4EM899xxef/11/OQnP8HatWtxySWXoL29PennlyxZgsLCwvhXRUWF4SOmEZ4Nk2hYWuagBU/NgWWNJAyIbctjw9L/2XaeNvYK5ODlgQce6FNQ2/tr3bp1ygd07bXX4sorr8TMmTOxaNEi/PnPf8bOnTvx8ssvJ/38fffdh4aGhvhXVVWVsu+gsO2BYesoVEiObfrjICz641ntprNJXTguDMdpJrzbyLaRHxRqXm6//XZcd911A35m0qRJqsfTh7KyMkycOBG7du1K+v9ZWVnIysryzV/QkNOJzFXoyj4tuBd8qe0hX091Wz9g0Z9tKXRVn+Tz5N2Z0YZ+wdK9K624txNs/TsMY5CDl+LiYhQXFwdxLEmpq6tDVVUVysrKjPkUzMMRSPhB2JYnqmDDQ9Janzq2NowwfMCGs+QuErZRCoHWvFRWVmLTpk2orKxEd3c3Nm3ahE2bNqG5uTn+menTp+Oll14CADQ3N+Oee+7BmjVrsHfvXqxcuRKLFi1CcXExPve5zwV5qMGiIQxbVyTY8HoAnow0w+sBQqg/sk/mDQuHNJY9GMPS/1l2WfoQ6FLpH/zgB3j66afjP8+ePRsA8MYbb+Ciiy4CAOzYsQMNDQ0AgIyMDGzevBnPPPMMjh49irKyMlx88cV44YUXkJ+fH+Shpiz+VL6b35iMA52VF8o+qW3L/G4ZKrbqj2XDONGf79iqPw50pkptzLwEGrwsW7Zs0D1e3FFuTk4O/vKXvwR5SOzYllq2CsvaVqY2AvTJ4tQul6K/IYZlK890Sbml0kMRnRvW1nS2DYV5HLnaxH02zHQYYdQf2ad5l6I/D9iqPxv6P/f1tzHzIsGLYXTSvOo+aZ/nSGf7AXmqgMEpd52D6K8fn75MG9EQ/Q1OWPTnBxbGH1pI8GKY0KTQiXDMbYcRaaGhhW3XU/q/5HD3fzb2nRK8CH1hGZ7ZWiYnCPbAspJLw2V4egXed07ZF7pI8GIEnY3JwrJhFwcsUwXMDw+eFDrRwNLgmWNKhAp3MEDXQiCHMYjPEO45ZUF/3RsJXoSUgCOQ8AMb062mCUvwLPrrB+46GwvyCtz7vNjQRr2R4MUEOsK0dHRmw7uNeDZDYyCE+iP7DIn+uLHhIRmW/s92JHgxgHu5H30jK/MpCVtvJJ4pOaJPlmAgfPrjeHiI/vzHVv1xoPVuo9SPKfsgwYthwpLO5kBWMgxOWPQn55k6sL+3h9V76iIFu8KgaFXbM4/UVbHhDbu2ti3ZZwj1Z4XPsLStTsG4pW3E0f9Jwa7gOwm7WmrYqmLDhl1+wPMOndTfSlP0Z8inDfpjRvQXHHrF23bpCJDgxTwhSWdTsbXDsA5ppCGFDUuPuUtIbJA8T52NvB5AGGLYWrwoCMLA8KzkUn+Vta0rsmwo3vbLlgsJXggMz4w219SSESQ7nRvQj5uXI6oem59N+vy4kTkBHUn/cLQt+8PDoG0Mnikcok99l1ZsWJjuQ2cwJj9L2acf/oPGj8tSmDOM9PniEbQ2TQZ12iihYNeC69IbCV4I/PH2c/G52ePwyxvPJNm5hTwqdzjJ9s5LpwIArp1bQbIDgJnjCpA9LB1zJo4i2f1g0ckAgNsuOoHs84kvn4GvnjsZi04rJ9ndeuEJ+OJZE/D0V88i+zxxbDSY/MzscSS7c6aMJvuKkZ+VCQA478Rikt3E0bnKPv/XnPEAgNsvPpFkp6O/71w+HQBw8/xJJDs3JxGD/StmlQIAKorUA9p5U2jX5bOnR7UT0xKFWPt+6uQSkt15J0WPcVgG/cFxfo/tTcTrcuakUThjwkhc06MlCv/9T/Nw5qRRePortHt0VN5wfH72OHz29HLyQ/qZr56N/KxM/Md1p5Ps3JxcXkD6/Lcuid5fnzmd1ocBwMOfn4WrTi3D52bT2vebF5+Ay08pxc+vP4Ps89uXTUNJQRb+94KpJLuTywpw1qQizBxXgLnEZ0QqkObYmpfrh8bGRhQWFqKhoQEFBTTRBklzexcAYETPQ88rjuNge00TppbkIyOd1slFIg46IxFkZWaQ7ADgWGc3sofR7Tg41tmNqiOtOKkkn2y7ZX8DSguzyZ3q/qNteH17La6ZM57cTmv3HsHInGHk4+3qjmBbdRNOKS9AOlELze1dcBwH+dm0EaHjONhX14qJo3PJo7P9R9vQ2NaJGWW0+7CzO4K/bavFWZOLUJRHC7bqmttR3XAMM8cVkuwAYHdtM8aPyiFfzyMtHdhR04RzphSR22h3bRPGjMhGYS7turR3deOj2hbMKMu3ctRMIRJxyHoHgIbWTjR3dClldWsajmFsfpaSXw4cxxkSOqA8vyV4EQRBEASBHcrzW6aNBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwCgleBEEQBEGwikzuA/Cb2EuyGxsbmY9EEARBEASvxJ7bsef4QAy54KWpqQkAUFFRwXwkgiAIgiBQaWpqQmFh4YCfSXO8hDgWEYlEcODAAeTn5yMtLc3Xv93Y2IiKigpUVVWhoKDA178tSPsGjbRvsEj7Bo+0cbBwt6/jOGhqakJ5eTnS0weuahlymZf09HSMHz8+UB8FBQVy4wSItG+wSPsGi7Rv8EgbBwtn+w6WcYkhBbuCIAiCIFiFBC+CIAiCIFiFBC8EsrKy8MMf/hBZWVnchzIkkfYNFmnfYJH2DR5p42CxqX2HXMGuIAiCIAhDG8m8CIIgCIJgFRK8CIIgCIJgFRK8CIIgCIJgFRK8CIIgCIJgFRK8eOTnP/85Jk+ejOzsbMyZMwdvvfUW9yFZwZIlS3DmmWciPz8fY8eOxWc/+1ns2LEj4TOO4+CBBx5AeXk5cnJycNFFF+HDDz9M+Ex7ezu+9a1vobi4GHl5efiHf/gHfPLJJyZPxQqWLFmCtLQ03HXXXfHfSfvqsX//fnz5y1/G6NGjkZubi9NPPx3r16+P/7+0rzpdXV34/ve/j8mTJyMnJwdTpkzBj370I0QikfhnpH1pvPnmm1i0aBHKy8uRlpaGP/zhDwn/71d71tfX44YbbkBhYSEKCwtxww034OjRowGfXeKJCIPw/PPPO8OGDXN+8YtfOFu3bnXuvPNOJy8vz9m3bx/3oaU8l112mfPrX//a2bJli7Np0ybnyiuvdCZMmOA0NzfHP/Pwww87+fn5zvLly53Nmzc71157rVNWVuY0NjbGP3Prrbc648aNc1asWOFs2LDBufjii53TTjvN6erq4jitlOS9995zJk2a5Jx66qnOnXfeGf+9tK86R44ccSZOnOjcfPPNzrvvvuvs2bPH+etf/+rs3r07/hlpX3X++Z//2Rk9erTzP//zP86ePXuc3/3ud86IESOcpUuXxj8j7UvjlVdece6//35n+fLlDgDnpZdeSvh/v9rz8ssvd2bOnOmsXr3aWb16tTNz5kznqquuMnWajgQvHjjrrLOcW2+9NeF306dPd7773e8yHZG91NbWOgCcVatWOY7jOJFIxCktLXUefvjh+GeOHTvmFBYWOk888YTjOI5z9OhRZ9iwYc7zzz8f/8z+/fud9PR059VXXzV7AilKU1OTc9JJJzkrVqxwLrzwwnjwIu2rx3e+8x3nvPPO6/f/pX31uPLKK52vfvWrCb/7/Oc/73z5y192HEfaV5fewYtf7bl161YHgPPOO+/EP7NmzRoHgLN9+/aAzyqKTBsNQkdHB9avX48FCxYk/H7BggVYvXo101HZS0NDAwCgqKgIALBnzx7U1NQktG9WVhYuvPDCePuuX78enZ2dCZ8pLy/HzJkz5Rr08M1vfhNXXnklPvWpTyX8XtpXjz/+8Y+YO3currnmGowdOxazZ8/GL37xi/j/S/vqcd555+Fvf/sbdu7cCQB4//338fbbb+OKK64AIO3rN36155o1a1BYWIizzz47/plzzjkHhYWFxtp8yL2Y0W8OHz6M7u5ulJSUJPy+pKQENTU1TEdlJ47jYPHixTjvvPMwc+ZMAIi3YbL23bdvX/wzw4cPx6hRo/p8Rq4B8Pzzz2PDhg1Yu3Ztn/+T9tXj448/xuOPP47Fixfje9/7Ht577z3ccccdyMrKwo033ijtq8l3vvMdNDQ0YPr06cjIyEB3dzceeughfPGLXwQg+vUbv9qzpqYGY8eO7fP3x44da6zNJXjxSFpaWsLPjuP0+Z0wMLfffjs++OADvP32233+T6V95RoAVVVVuPPOO/Haa68hOzu7389J+6oRiUQwd+5c/PjHPwYAzJ49Gx9++CEef/xx3HjjjfHPSfuq8cILL+DZZ5/Fb3/7W5xyyinYtGkT7rrrLpSXl+Omm26Kf07a11/8aM9knzfZ5jJtNAjFxcXIyMjoE03W1tb2iV6F/vnWt76FP/7xj3jjjTcwfvz4+O9LS0sBYMD2LS0tRUdHB+rr6/v9TFhZv349amtrMWfOHGRmZiIzMxOrVq3CT3/6U2RmZsbbR9pXjbKyMpx88skJv5sxYwYqKysBiH51+fa3v43vfve7uO666zBr1izccMMNuPvuu7FkyRIA0r5+41d7lpaW4uDBg33+/qFDh4y1uQQvgzB8+HDMmTMHK1asSPj9ihUrMH/+fKajsgfHcXD77bfjxRdfxOuvv47Jkycn/P/kyZNRWlqa0L4dHR1YtWpVvH3nzJmDYcOGJXymuroaW7ZsCf01uPTSS7F582Zs2rQp/jV37lxcf/312LRpE6ZMmSLtq8G5557bZ2n/zp07MXHiRACiX11aW1uRnp74GMrIyIgvlZb29Re/2nPevHloaGjAe++9F//Mu+++i4aGBnNtbqQs2HJiS6WfeuopZ+vWrc5dd93l5OXlOXv37uU+tJTnG9/4hlNYWOisXLnSqa6ujn+1trbGP/Pwww87hYWFzosvvuhs3rzZ+eIXv5h06d748eOdv/71r86GDRucSy65JLRLIQfDvdrIcaR9dXjvvfeczMxM56GHHnJ27drlPPfcc05ubq7z7LPPxj8j7avOTTfd5IwbNy6+VPrFF190iouLnXvvvTf+GWlfGk1NTc7GjRudjRs3OgCcRx991Nm4cWN8aw+/2vPyyy93Tj31VGfNmjXOmjVrnFmzZslS6VTksccecyZOnOgMHz7cOeOMM+JLfYWBAZD069e//nX8M5FIxPnhD3/olJaWOllZWc4FF1zgbN68OeHvtLW1ObfffrtTVFTk5OTkOFdddZVTWVlp+GzsoHfwIu2rx5/+9Cdn5syZTlZWljN9+nTnySefTPh/aV91GhsbnTvvvNOZMGGCk52d7UyZMsW5//77nfb29vhnpH1pvPHGG0n73JtuuslxHP/as66uzrn++uud/Px8Jz8/37n++uud+vp6Q2fpOGmO4zhmcjyCIAiCIAj6SM2LIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhWIcGLIAiCIAhW8f8DnTJgvjQeNmwAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "for pot in data_0:\n", + " plt.plot(np.real(pot.flatten()))" + ] + }, + { + "cell_type": "markdown", + "id": "79456f12-0f63-4ff8-9130-a1a111e1d4a4", + "metadata": {}, + "source": [ + "Change something on the phantom" + ] + }, + { + "cell_type": "code", + "execution_count": 43, + "id": "273c6b79-39e0-46da-b01c-475c59543a48", + "metadata": {}, + "outputs": [], + "source": [ + "# get data\n", + "data_1 = sciospec.StartStopMeasurement(return_as=\"pot_mat\") # or \"hex\"\n", + "# data.shape = (burst_count, n_el, n_el)\n", + "# check if data shape equals setup\n", + "assert data_1.shape[0] == setup.burst_count" + ] + }, + { + "cell_type": "code", + "execution_count": 44, + "id": "0f0b9cba-e52b-4781-857f-211166252121", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAhYAAAGdCAYAAABO2DpVAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXtRJREFUeJztvXmcHFW5///pnqVnycxkJvsyCQk72YBEIIiAiEhEELevchHx3uvXL15AuKgXcbksisGr16tXFIXrD1GuBpFd1kRIwpJA9n0l2ySZyTJL9yyZnunu+v3RXdVVvdY5dbqqe+bzfr0CPTP19FN16qlznvM8zznl0zRNAyGEEEKIAvxenwAhhBBChg50LAghhBCiDDoWhBBCCFEGHQtCCCGEKIOOBSGEEEKUQceCEEIIIcqgY0EIIYQQZdCxIIQQQogyyt1WGIvFcPjwYdTV1cHn87mtnhBCCCESaJqG7u5uTJw4EX5/9riE647F4cOH0dzc7LZaQgghhCigpaUFkydPzvp31x2Luro6APETq6+vd1s9IYQQQiQIhUJobm42xvFsuO5Y6OmP+vp6OhaEEEJIiZGvjIHFm4QQQghRBh0LQgghhCiDjgUhhBBClEHHghBCCCHKoGNBCCGEEGXQsSCEEEKIMuhYEEIIIUQZdCwIIYQQogw6FoQQQghRBh0LQgghhCiDjgUhhBBClEHHghBCCCHKoGNBipJoTMPaA52IxjQAQCymoaN3AGv2dyCW+F1r8ASCfYPG8ZqmoW8ggl+9sRsH2vsAAB29A3h67UH81+KdWHugE23Bfry7px2r93V4c2GEEDLE8WmaprmpMBQKoaGhAcFgkG83JRa6+wfxjb9sQDgSw7KdxwAAo0dU4gefnIlfL30fmw4F02Smja7Fjz8zG3c+tRFHQ/3oHYgaf/vS/Kn4w4r9WfV95+Nn4Jo5kzC+oUr9xRBCyBDD7vhNx4IUBZqm4bpHVmLlHncjCR8+fQwe/cfzXNVJCCGliN3xu9zFcyIkK6H+SJpTMWlkNc5uHonFW49gIBrDJ8+eiK9/5FSs3d+J9/Z24Mk1BwEAPh9w0qhaXHLaGPx9+xFcdvpYTBtdi+c3HMa+9j589eLp+OqHpqM7HMHqfR346Ws7sa01BCCeKiGEEKIORixIURDsG8Sc+14DAOy6fwH8Ph/K/D4AwImBKI73hDFxZLXxOwAInhjEttYQ5k5tREWZ/XIhTdPw4qZW3PKndZgzuQHP3XKR2oshhJAhCCMWpKQxOxDVlWVobqpJO6ahugIXTB8l/N0+nw81lWUAAFe9akIIGQZwVQgpCjSXh3gffPkPIoQQIgwdC1J0cMgnhJDShY4FGda4W2FECCFDHzoWpChwfYBnWIQQQgoCHQtSdPh8HPUJIaRUoWNBhjVuF40SQshQh44FKQqYCSGEkKEBHQtSdHDQJ4SQ0oWOBSkKXN4A1qTXE7WEEDJkEXIs7rnnHvh8Psu/8ePHF+rcyDDFjdpNFogSQkhhEN7Se8aMGViyZInxc1lZmdITIoQQQkjpIuxYlJeXM0pBlONVRoKpEEIIUYtwjcWuXbswceJETJs2DV/4whewZ8+enMeHw2GEQiHLP0Jy4UaagokQQggpDEKOxfnnn48//OEPePXVV/HII4+gra0NF154Idrb27PKLFy4EA0NDca/5uZmxydNCCGEkOJEyLFYsGABPvOZz2DWrFm4/PLL8eKLLwIAHnvssawyd911F4LBoPGvpaXF2RmTIYlXKQlmQgghRC3CNRZmamtrMWvWLOzatSvrMYFAAIFAwIkaQpTDRSGEEFIYHO1jEQ6HsW3bNkyYMEHV+RBCCCGkhBFyLL75zW9i2bJl2Lt3L95991189rOfRSgUwo033lio8yPDBK/e2eHVxlyEEDJUEUqFHDx4ENdddx2OHz+OMWPG4IILLsDKlSsxderUQp0fGWa4laLwcV0IIYQUBCHHYtGiRYU6D0IIIYQMAfiuEFIcMCNBCCFDAjoWpKhwK0HBVSGEEFIY6FgQQgghRBl0LEhRwHeFEELI0ICOBSkq3HqdOTMhhBBSGOhYEEIIIUQZdCxIUeDdu0KYCyGEEJXQsSBFhWspCuZCCCGkINCxIIQQQogy6FiQosC7d4V4opYQQoYsdCxIUcF3hRBCSGlDx4IQQgghyqBjQYoC71aFEEIIUQkdC1JUuJWi4LtCCCGkMNCxIIQQQogy6FiQosC7d4UwGUIIISqhY0GKC9dWhRBCCCkEdCwIIYQQogw6FqQo8ColwUQIIYSohY4FKSrcSlG49Xp2QggZbtCxIIQQQogy6FiQosCzxRnMhRBCiFLoWJCiwrV3hTATQgghBYGOBSGEEEKUQceCDGuYCSGEELXQsSBFhWvvCnFFCyGEDD/oWJBhDbf0JoQQtdCxIEWB2+M7izcJIaQw0LEgRQUHfEIIKW3oWJBhDRMhhBCiFjoWpCjQXB/iGRohhJBCQMeCFBUc7gkhpLShY0GGNVwUQgghaqFjQYoCrgohhJChAR0LUlTwdeaEEFLa0LEgwxr3i0YJIWRoQ8eCFAVcE0IIIUMDOhakqOCATwghpQ0dCzKs4aoQQghRCx0LUhS4/TIwFokSQkhhoGNBiguO94QQUtLQsSDDGqZCCCFELXQsSFHAVSGEEDI0oGNBigoO+IQQUtrQsSCEEEKIMuhYkKKA7wohhJChAR0LUlRwGSghhJQ2dCxIkeDN8gy3988ghJChDh0LUlS4FbDwsUyUEEIKAh0LQgghhCiDjgUpCrzKSDARQgghaqFjQYoKtxIUrBElhJDC4MixWLhwIXw+H26//XZFp0MIIYSQUkbasVi1ahUefvhhzJ49W+X5kGGKVykJLgohhBC1SDkWPT09uP766/HII4+gsbFR9TmRYQz3sSCEkNJGyrG4+eabcdVVV+Hyyy/Pe2w4HEYoFLL8I4QQQsjQpFxUYNGiRVi7di1WrVpl6/iFCxfi3nvvFT4xMrzwblUIcyGEEKISoYhFS0sLbrvtNjz++OOoqqqyJXPXXXchGAwa/1paWqROlAwPuCqEEEJKG6GIxZo1a3D06FHMnTvX+F00GsXy5cvx4IMPIhwOo6yszCITCAQQCATUnC0hhBBCihohx+IjH/kINm3aZPndP/7jP+KMM87AnXfemeZUEGIXr1ISXBVCCCFqEXIs6urqMHPmTMvvamtrMWrUqLTfEyID3xVCCCGlDXfeJIQQQogyhFeFpLJ06VIFp0GGO3xXCCGEDA0YsSBFhjspCq4KIYSQwkDHghBCCCHKoGNBigLPUiHMhRBCiFLoWJCiwrVVIUyFEEJIQaBjQQghhBBl0LEgRYF37+xgLoQQQlRCx4IUFa69K4QbZBFCSEGgY0EIIYQQZdCxIEUBV4UQQsjQgI4FKSq4KoQQQkobOhaEEEIIUQYdCzKsYSaEEELUQseCFBVurdZgJoQQQgoDHQtCCCGEKIOOBSkKvFsVwmQIIYSohI4FKSq4KoQQQkobOhaEEEIIUQYdC1IUePWuECZCCCFELXQsSFHhXoaCuRBCCCkEdCwIIYQQogw6FqQo4LtCCCFkaEDHghQVPpeWa3BVCCGEFAY6FoQQQghRBh0LUhR4lZHgBlmEEKIWOhZkWMJMCCGEFAY6FoQQQghRBh0LUhR4lZJgIoQQQtRCx4IUFe69K4TJEEIIKQR0LAghhBCiDDoWpCjwLCXBXAghhCiFjgUpKlxLhbijhhBChh10LAghhBCiDDoWpCjw7F0h3qglhJAhCx0LUlT4XEpScFEIIYQUBjoWhBBCCFEGHQtSJHi0QRbfFUIIIUqhY0GKCvdWhTAXQgghhYCOBSGEEEKUQceCFAVcFUIIIUMDOhakqHArQcFVIYQQUhjoWJBhDWs3CSFELXQsSFHA8Z0QQoYGdCxIUcHXmRNCSGlDx4IMazTGSgghRCl0LEhR4HatAwMjhBBSGOhYkKKC4z0hhJQ2dCxIUeDV1tpcFUIIIWqhY0GKC7e29GYuhBBCCgIdC0IIIYQog44FKQq8ykgwE0IIIWqhY0GKCte29HZJDyGEDDfoWBBCCCFEGXQsSFHg2eoM5kIIIUQpQo7FQw89hNmzZ6O+vh719fWYP38+Xn755UKdGxmGuLVag4tCCCGkMAg5FpMnT8YDDzyA1atXY/Xq1bjsssvwyU9+Elu2bCnU+RFCCCGkhCgXOfjqq6+2/Hz//ffjoYcewsqVKzFjxgylJ0aGF169s4PvCiGEELUIORZmotEonnzySfT29mL+/PlZjwuHwwiHw8bPoVBIViUZBri3KoS5EEIIKQTCxZubNm3CiBEjEAgEcNNNN+GZZ57BWWedlfX4hQsXoqGhwfjX3Nzs6IQJIYQQUrwIOxann3461q9fj5UrV+JrX/sabrzxRmzdujXr8XfddReCwaDxr6WlxdEJkyGKRxkJviuEEELUIpwKqaysxCmnnAIAmDdvHlatWoVf/OIX+O1vf5vx+EAggEAg4OwsybDBrdUaXBVCCCGFwfE+FpqmWWooCCGEEDJ8EYpYfOc738GCBQvQ3NyM7u5uLFq0CEuXLsUrr7xSqPMjwwTuj0UIIUMDIcfiyJEjuOGGG9Da2oqGhgbMnj0br7zyCj760Y8W6vzIMMOt1RrMhBBCSGEQcix+97vfFeo8CCGEEDIE4LtCSFHg1eoMjctCCCFEKXQsSFHh2moN5kIIIaQg0LEghBBCiDLoWJCiwLt3hRBCCFEJHQsyLOG7QgghpDDQsSCEEEKIMuhYkKLAu1Uh3uglhJChCh0LUlT4XFoWwneFEEJIYaBjQQghhBBl0LEgRQEzEoQQMjSgY0GKCu6PRQghpQ0dC0IIIYQog44FKQq8fGcH3xdCCCHqoGNBigq3Vmu4tfqEEEKGG3QsCCGEEKIMOhakKPAyGcFMCCGEqIOOBSkqXEuFuKOGEEKGHXQsCCGEEKIMOhakOPAwHcFMCCGEqIOOBSkq3HqdOReFEEJIYaBjQQghhBBl0LEgRYHmYUKCG2QRQog66FiQosK9VSHMhRBCSCGgY0EIIYQQZdCxIEWBl9kIJkIIIUQddCxIUeFagoKZEEIIKQh0LAghhBCiDDoWpCjwNBXCXAghhCiDjgUpLlxaFsINsgghpDDQsSCEEEKIMuhYkKLA09emc10IIYQog44FKSrcylAwE0IIIYWBjgUhhBBClEHHghQFXr6vg6tCCCFEHXQsSFHh2rtCuCyEEEIKAh0LQgghhCiDjgUpCpiNIISQoQEdC1JUcFUIIYSUNnQsCCGEEKIMOhakKOC7QgghZGhAx4IUFW6t1uCiEEIIKQx0LAghhBCiDDoWpEjwcIMsrkkhhBBl0LEgRYV7q0KYCyGEkEJAx4IUBSygJISQoQEdC1JUeFFUSaeGEELUQceCDEu4KoQQQgoDHQtSFDBoQAghQwM6FqSo8KKokk4NIYSog44FIYQQQpRBx4IUBd5u6c2YBSGEqIKOBSkuXMqEsHiTEEIKAx0LQgghhChDyLFYuHAhPvCBD6Curg5jx47Ftddeix07dhTq3MgwwstttZkIIYQQdQg5FsuWLcPNN9+MlStXYvHixYhEIrjiiivQ29tbqPMjwwxu6U0IIaVNucjBr7zyiuXnRx99FGPHjsWaNWtw8cUXKz0xQgghhJQeQo5FKsFgEADQ1NSU9ZhwOIxwOGz8HAqFnKgkQxRvV4V4p5sQQoYa0sWbmqbhjjvuwEUXXYSZM2dmPW7hwoVoaGgw/jU3N8uqJMMAt1ZrcFUIIYQUBmnH4pZbbsHGjRvx5z//Oedxd911F4LBoPGvpaVFViUhhBBCihypVMitt96K559/HsuXL8fkyZNzHhsIBBAIBKROjgwfPM1GMBVCCCHKEHIsNE3DrbfeimeeeQZLly7FtGnTCnVeZJji1moNZkIIIaQwCDkWN998M/70pz/hueeeQ11dHdra2gAADQ0NqK6uLsgJEkIIIaR0EKqxeOihhxAMBnHppZdiwoQJxr8nnniiUOdHhglevq/Dy825CCFkqCGcCiGkkLi3KoTJEEIIKQR8VwghhBBClEHHggx7GIgjhBB10LEgRYVrqRB31BBCyLCDjgUhhBBClEHHghQFnr4rxDvVhBAy5KBjQYoK1zbIYi6EEEIKAh0LQgghhCiDjgUpCrzcpIr7sxBCiDroWJCightkEUJIaUPHwiGxGGe7Q4lwJCodwRiIxBSfTX5k7S8SjSESlTvf4IlBKdloTJNuI9lz1TQNUQdtJGsLjIKR4cyQcSxuW7QO//K/a7C9LSQkF+ofxJcffQ+f+OWb6B+MCsn+9NUdOP37L+PRt/cK67zy58vxyV+9jeCJQSHZFzYcxmX/uRR/XLlfSC4SjeGrf1iNbz25AX0DESHZ1uAJLHx5G97adVxIDgB+sWQX5i/8OzYdDOY8rhhWhazZ34lz7luMBb94U3gwum3ROpz576/g6bUHheRaOvpw4cK/45oH30I4ImZ///HKdpz+/Zfxe0H7O9YdxoUPvI4P3L8Ex3vCQrKPLN+DOfe+hn/760YhuRMDUVzykzdw7g8WY/fRHiHZp9YcxJn//gq+/ZSYzlhMw2ceegez7nkVaw90Csm+vv0IZtz9Kv7vH1YLyQHAVx5bjRl3v4rlO48Jya090Ilz7nsN33xyg7Bj8uTqFnzlsVXYfCj3c5ZK8MQgblu0Dvc8v0VK5yd++SZe29ImJNfdP4j/89sVuOQnbwjb31/XHMR59y/Bw8vfF5IbiMTwD4+sxFf/sFq4z93eFsKtf16Hv208LCQHAN9/djM++9A72N/eKyS38WAXrnnwLXzxf97Foa4TQrL//PtV+Kffr0JX34CQnEqGjGPx3PrDeGlTGx57Z5+Q3NIdx7B0xzFsPhTCPsGb/+AbuzEY1fCz13YKyW1sCWJ7Wzc2tHQJdwT3vrAFe4714vvPbhaSe/9YL17begRPrjmIpTvEOrz7XtiK3y7bg6/+UayTDfUP4r+W7ERrsB8r97QLyXrBmv0d6BuIYntbt/BD+dz6w4jGNDz69j4hudX7O3A42I+NB4PYe1zM/n699H0MRjX85NUdQnI72rpxtDuMzr5BbG/tFpL90cvbAABPrzskJLevvRcHO0+gJxwRtvk/rNiHwaiGRatahOQ6+gaw9kAX+gaiWLtfzLF4cWMbwpEYlmw7KiQHAEu2HUHfQBTPbxAbiFbv60Bn3yD+uuag8OD3rb9uxJJtR/EfgrawfOcxPLf+MH7/zj68f0zM/r71143YfCiEe1/YKiS36WAQ7+3twP72PmH7u++FLTjaHcaPXtouJLfneA/eeb8dr209gjd3ifV/P355O17YcBi3/GmdkFxH7wD+uHI/Vu/vFLajFze1YuPBIN7afRxLth4Rkn19x1G8vv0oBqPezdaGjGMxd2ojACAsGGo1h1hlZ819gpGOmElRTFDp8R45L9Q8AxcNR+9r7wMA9A2IXWd40P3UgBPMt0L2kRS9nzFTE8naX7/g/TSfo2jRrOw5WttW7EsGJDtIy3UKfoWKVIa4zuTniGT6pl0wAjBo6v8iMbnntS3UL3S8+dpEbSHULxZt1TFfWkTQntpCYm1q6DEpFU1Zmo8XtQXdjrwsIxsyjsWVM8YDEH+YzffMrXC8WY17OuWdGWmdAgOYp6mQhG5v7stw0Sk/yDtQmlG/8NcInLATh8SL+2Lu/yT9CmG8tj+3+j8n9mdxxCXtz++hZzFkHAu9DUUfbJHBTxXWGaM7WA3VJZ1Z9OfCzdUaqapkH2YneG1/7jmZZv2uqLS0phOdIk3kRI+TSJIsXtif5rH9FXP/l0lWBLP9ebnubQg5FvFmFL0hXnjPsHSyHjxYrmj0RqcTLLNqt3R60OHBa1tw6UKdpELMiIg6ilh4Mfh5rdMdlZ7otDrwTmTlIhZMhShAb0PRG6gp6nyEdHrQs3sRChTp2L10PPS28WZWM/zsz5vBxKVUiLSWVJ2KviifnmFof95M5pykQgTkTJ+93KtnyDgWfulUSPKzawOuuWDPtfCj+QdXVFqjQTaVuvko5NLlnvOV/OxaKsRkf55cZxE6trm/R06nsB4HRY2yqHK+ZHV6Mch70v+5FDGLMWKhFiMV4iCX5UEk2rWCKc/ztyWQC7Hmm93SmfzsRe2Be+FvD2bGipwZkedF1QDili3EPLA/LwrmvagrsjiKoqvFpFMhyc+ssVCAUbwpOCRYZzUeDLiuaPSm0xIK4Xm6LEQ/B9Ov3LIFDD/78+JOOxrwS8wREtJJ+3NHv3CKXl5Wh6tCFCAdsfDYY/dk6WcRd5RuPgupOUgvZvJezxi9sD8vanzcWxXiIBXi9TPqjsqUSc7w6P+cFG/KyjEVooBk8aZojYW3yRBvjNyL2ZArKh3hzQqN4Wd/3gyapZUK8SZFNVwGefcdW+HiTfP3CHgl1lQIIxaOSe5jISbnzcOc7QwKqNN1jamFibnx0u/QjP+7P6v22v7cK6R0X6eqtpWtyneiyLU6qMzqxb7DQcH8UO7/nLStbCTJfCwjFgrwy+5j4Xko2iWdHlRMycxMvFwV4kVY2HxfPCmec0elJ3l1VfVTIk6ms1RI5s+FxBObH5b9n3wkXdb+6FgoQG9DJ9W3br0i2Zu140m8GcCKPxdinWF4MWMcHvYnPnuTbZfM+oW/R1KnuB730xKWOhTJjkF0vwSv7a8UHHjZ1KH5WBZvKsCvoHjTvdmbWadb3rP5s/vBwLwqPfQ7jHNTVOwnpRtezRg9GMA8sD8nOjWBJeHOIiOm75H+FjFU2F8p7B3kpP+Tv6fyDpTshnKWnTfFVCplyDgWeisKF2+aPssauajHrqpaXQSv92iwq9PLd4VY70XpDLjCM0YPbMGy2aKgUlmbUBX9dq140wPn30kkSRZv3pXkbcpR+K3HkrZgKd5kxMI5RipEUE5T4LI7ez+A+x2IFw9WKSwL8XwVyzCxP0fPqJCcvM5s35P3WCd6PFmt4zz94mQAKwX7k70+J/YnnQoxffazxsI5pZQKcS8UnMSb2ZD9mYmXNRiZ3hUytFMh7g9gnmxEpygVIls85+j9EEKS8niu04vnzKVUiBNH0dp3ytkfIxYK8EmnQrwIRSc/F+MLwVQh04G4uyokxwZZQzoVYvrsSY2PmKx0KsRcGyHcscuJarKCSA3Xe+D8e5IKGbr9n+WdUI7qUOTkvGTIORaiqMjDOjEaT6ryS2jQdBP5SmwHqYgs+oW+Y4jbn2y43skbVWU3cXKmM/Pn/HJOojHmz8XvZErrtHwWUyqdCvHQ/rxMgwBDyLEo1VSIe/UODsJysh27QNt66XfoumWjV6q2ix4u9ufW9tqWWZ+gUhXFc84K9uzLKSsYlf8aMZ2mz8Ol/3Pb/rxMgwBDyLHQKYlUiPmzBy67k4dZVtZ2KsTTHbJMH12bMQ4/+3OSCpFNSwjPGLN8FtIpnH6RdWy9tT8nOody/+el/XkcsBg6joWKl5C5VpXvySxBy/jZDpaO3aVCNi+QPUNlV+bBlNGLvLqjPSVkUyGOZqkCcopSIUJycmLKvoj9nw2dbkVJjFQIIxZK8MsWb3peSOSSTkshkZis7L4HIk6bl26HrlvFNrqimEOkntifB8WbThCZMTqpH5BNaTjT6b79sf8TQ7b/c2Z/9uWMY1ljoQa9wl/UZiwdgcAOexbdjkLRcjpFURX+FnlA5IJBLm6QlfKzfPGc/DmYRWXtT3xLb7OsnE5RnNif7ODnrKhW1pmWn47L3hcn95D9nxiuOXyS0Q792li8qQijIR2kQlwLy3miUz4sZw0F2pfzJJfqAMsgL1mkJazTA1sYNsWbWfTbwTzIyr62erjsY+HWUkonOOn/5NNimb/Dnk7z94jbn5evTAeGkGPh5T4WonhTMJX87EmnlS8V4qHfoT/0sp2skxAy97HIj6xj62gwkWwXy6ApGAGQtQV1M2MPUiGupV8yfxZFNpIkqlN2FYvhWDBioQrnqRDp4qkSCEVDkQMlNpMXDye7+UCk6pIPP8qfgwpbKDX7c7SPhcTsTU5n5s955ZxcZxb9eeVUpUIkv6cUUiFO7M+6Kkm2/3NifwI6WbypFj0VIl6hbDFzZeeTW6X4gOtcZ+bP9mRVhALzHet9qkT2YXYUdfLa/jyImLm23E9SLvV4+eJNMZ2apC04uoMK7M9Jn+va9u5F3v/llBWQ022Oy00VoXuVwg+z6bNrHrvkAOYE2Q1X0r+nsDNGLx4I/TylV7840O3JjNFyX9xBc9DLSi93jskNCAmJDJ9sSDnK5Sc/i60EkL+LXkQPvEjLet3/OXEypaLETIWoQfbtpir2kHdU5e7Sk+Vkfb3sBkVe7M8vQs53hbgUilbxroahbn9qHD63QtGqdMrJiaKiqFbcsS0t+1PR/7lmC4n/MxWiCL9+JS5V3zqh1Pbnlw4FutQ5qkKTvC/OXltu+jyE7U+VTk2gIFJdIaWITvmBWnoVlRPHVoH9Oanx8aR43YP+z237Y/GmIvTZZ0mkQsyfXQs/yoXWUpHfec5mKsTDB8KL4s3hUjynbLmzW8VzOX7KrTOzflGdwyYV4oH9udb/eWh/rLFQhN73iN5AFXstONvwx30chQIlB9wiCEikkbYqRHL25uTaVCz3c2J/nkTMhGUlHb4s+u3plJOVD357Y39eLH0vtf7PIlsi9sdUiCKM4k3R3eM8sGxvijflZ4wqtvTOp7MYHA/ZjsBRLYnH9ueaTk+WO5vl5PXIRg+cFOzJ2p8w7P+kvycf6uzPvrB+LFMhipAt3rSGH+UsTjwUrcbIhXQqepjFtvQWDwV6sWOcfkkqiudEW1a2AzHfw1KzP0c1AQLHqtqISXrvDAepEOlcvqgzI6nTTKnZn7Ot3u2jzP6kbIERCyXooR9Ro5FdkuZkLwAvCvZUpSXEZlKZPxcL6Y+e80HekS3IOjOiNm+5L046PJGZlLxOTVLWyaBpcfgEoqBevBDMSTrN3P95YvNDuP/Tsv4gpkfGmeG7QhRh1Fi45LGr2uTH2fay4iEyQOzdB3E9cjqFQrtexOdznILI2TgZwKRnqWJqrLJedLJOUnGQO2EvBj9ZJyh+fObPIjrZ/2XTKd//ydbhOallko9kxv/PVIgikqkQt26gmtCqs3yfiM7Mn+3pMc+k5HTalfTigdDvh7wDlf5ddill+5N2vkT1yDp8DhwoJ4O1rJx1kC8d+/MqFeJW/+e1/cl4/yzeVITPSIWIycmHnMT0ZNepJnSa/+DM+gVFXZulukFqhyjb4Tkr2DPrF9EppiebTi+W3rk2e4P8fZEt9nNSJCjvwDuIBimJ0slHg9yyP0/6Pyf3RTItaxRviqlTzhByLOL/d6tCWpVD4NaM0cmyMqdhudTPGY8VOqPCIB3+tnyH+/bnZOdDJxR7uN5JHYp1wihn88J9keT5KtvwTNKBElajyP5Kqf9zkhaTiZKI9gmqGTqOReL/oiarwmhEcdLhmZEtZBNGSSjQnqQnqZDEqSmJzAjqVuO0yc8YVaVU8h/rAAVRHWH9ko6QluWzPVm5CYfX9ufVu5KKvf/z0v5KrsZi+fLluPrqqzFx4kT4fD48++yzBTgtcfx+d1MhqnKDqkLaIsc6WwYnHpZL1V8spD57muR9cTQbHyb2p2zGKHDCToo3ne4jIKXTtPpE9tl2a/WBqv7EC/tzq//z0v5KzrHo7e3FnDlz8OCDDxbifKQxIhbCRiPpsTtJhSgqXhLBSR2ArIMg0t8Vg+Mh2xE46WTl2zZ5sFfFc/KOraAescMNHL0fQlK/s+t0Hj0QXsorGyVx0P95MeGQnTTEj1fQ/7lsf17sB2SmXFRgwYIFWLBgQSHOxRFG8aagnGznoywVougBzavT8lmw85FuI/EH0pMNslL+n/o5v7zpOiV1A96kQpy9Z0Jy9uaS/cGBAyVfB+XAFmSjBw50ytZnsP+zI6fG/sRON36w1/tYCDsWooTDYYTDYePnUChUED2yxZuy909V8ZJ7qRAHRm7puCQHv6Ioz0wh7V0hchELVZs/eWF/Tu6K7OzNre2NnThQ1gmHO7bghU6v7c+L/q/U7E8sFRL//5Av3ly4cCEaGhqMf83NzQXRk0yFiErKDSbmI528XdKJ6y8TIgMkHixT7ld27Xi+k3XrnQE5z8H82YOZSanZX6FWCOWUFZBzUjwnHb1ycJ2yERZnq8zkTtiRQ2r5wf3+r+TsT8gW4v8f8stN77rrLgSDQeNfS0tLQfT4JfexUFEw5WxXP3nkN1RSM9vIh9S6fE9WhWiJ/5t+JyDvKJcvG2a12K3obEjRjFFSp1tvY3VUS2LWKdBIsntRpKgUTIVk/g47SE8aFEUahnL/p8r+RKK9xVK8WfBUSCAQQCAQKLQa05begp2WbMRC0UzPWUhRRKfps+AbYGUHIifOlxukrQoxfRa7Lw4GTemImRq7cf/lSDKpOBU6RZ2vzPoLq1PSgXLitLH/syHrfv/ndHv3IZ8KcQu96M/JLME1j13R94hcrNVfl4+wiCgVmSMUg9shG1J2MgNT8X4IUVTZn2zqxi37cxQ9ULBawhMHSkylfJSE9pcXL+2v5Io3e3p6sHv3buPnvXv3Yv369WhqasKUKVOUnpwIsi8hk52ZOPG0nRic7DkoK94U6vDE5bx4HrSU/wNyBVOicoCT9IuqGaP013hUPCcmm+k7bB0vmWpyEg2SjpJk0S+qU7p+ShDZyVyu78l/rPv9Xzb99o6Xk03WWJTYctPVq1fjwx/+sPHzHXfcAQC48cYb8fvf/17ZiYkivSpEQfGmMA4GItlzcPIwK3nrophKV0gLF8p2Go5mqd7an2vL/bTMn23JZvke13SKyJmfFcGQu/yqJC3jZ1s6Je3Pi9fDp32PkM7kZy/6Pyc6Za6z5GosLr300qLMl/sl97GQLIp2Fn7M+oPg98gWEjkIBYpt42y/0/LSpHTd8pEZ02dJ3a7qzPqDg+/Je6z79ueskFLSFiCv06pf7lgn/Z+QWBHYn3Qhb8nZn/h1ssZCEbLFmzFJ79lJaM1cZe7kVdnSOXkHvY9QwZTEihs3H4i0gIVkJEnViodSsz/Z4jnh6IF5Ri4QBXBSPCdvC/I65fuizJ9FdbqWWlWgU1Svpa90cL7yETMn9mdfztjHQkibeoaOYwG55aYWexOamaT8LFtIJOscQNR7Nn8W7fDkZK3PcvFFuVJREj0YTvYnnQoR7GQtnwt/nal6hCIzjgbc5GcRB0rdezDk5ET1qtrt2K23zpqPlu7/XLY/r1MhQ8ax0KtgRYcvVVX5boQuU43arRCtvMduX85bt0NL/Fe+jWTlvLY/sdlQiv3ZF/UkLSEbAUjV48XqAy/sT2wAy/2zXdmh3P95aX9+pkLUIL+lt9zNT9UjGzp38urfop9VS+h083HI/XbTwt/PuE5v7c/Jq89F7SipX/54IZ2SEYC4TufOtFtbR6uyP1mdTmS96P+Gg/0xYqEMj1MhYmrldKbOEiRz8uJL72QHInmP3Qtk87BFkQoRUyunM+VY1zZxknaEHMwYs3xPXjnJmWaqTtpf/mPla5KGvv2xxkIRftmdNyUNzkkhkeyqmvRZgohO02dBvbKdj4in7+VKo6Rq5zN5t+wvdZtp+RSV/CxV5EpVbXsuU8gmoxOSsk7eD2HRWeT2x/4vP17aH1eFKMIn+a4QWW/fyezNSSGl9RzkdLpVIS/zMLv5PKQ+fPIDmLlzFjsHJ7PNbOeQ/1iTTnmV8sWbTmZvsjoFL1R+G2d5W5DV6YX9sf8Tk3NmC+IOH1MhijDebiooJ+vJOhkEVOyXICwrqTN+vHkm745OL5CdmTiaDZWQ/SkrnnPJ/hwVz2XRXyi5dFkBuWFifyXd/zmyP3E5Fm8qQm9I8aWUcp6hs+Kl5GcnBXvSnY/ogyWp0zrrKl7PQj8zFevrS83+ZGdgmX4uhM5UWdlB3q1txGUHoXRZsSEs03fYQb6QUk1f5MTmS6n/cy+SlIhYiKlTzpBxLGTfFeLUM5SStTxYDnS6VbypSXZcEiFaT1eFZPmcD1kbciKryv5krzP+c2FtIZMeIVlZu03TI+d8OdsvpvBygKpJgzupuFLu/xzplCneZCpEDYZjIRpykp0NpXrPYloz6hfWKTkQic+k5GSd6PQCTbKXtbaPqP3JzXC9sD9HxXNZPtuSlbQ/J7Uk0jPGLN8hKuue/WXRn1eO/V8+vLQ/Fm8qwmekQsTk5PO3qT8LeOySb050otOTUKB5mWoeQSc5W6fous2nIFs8J25/5s/u2IKKN3eKyjrb9tz0PQIN7GRPCVlZZ4WUsn2RIp0O7M+NXVidRUmctJHpezyxP/HrZCpEEUZDujRopioS8yoljTz1YRYQVfU+C/ktve3h7btC5NpIk7nQDIcPZftz5thKRnWynYCgrGv7WEgPYJk/25LN8j1iki45X8rsT9QWvLY/AbnEwUyFKCL5dlOXBs3UQ0U6dkuITM7rzngOuXRm+WxLVrLjctLheY1YB6Jl/GwH+WVlaSchJeuseK6wTmZST+bP+fA8FSKqU3oAk7c/+ZfgpZ2ElKzshlMZzyGnbObPtmRL0P64KkQRyS29xeRkb76q4iVHxZuSg4JbxUsi74fwcmdOXbcme18k5VJlvbA/BxNG6foMt+zPi9U6zna4Neu3L+fE/lQUjMZ/Lnz/5yT9Usr2J6OTEQtFGPtYOAk5CQ3yKd6zkE6z0QjIpc4YBXSaD3YULhUTNem0J+nu85CyQRZkH2bzZ9pfFqVSOlP1yNufvE5ZPc7SErS/NLlhan9CbZs42OdxlcXQcSyMVIgYSvLqgrIWIxdKheQ+B5sqPUmFlAKy56vuzYkCtpDyUiNp+5OcacZ/dmfGqCJ1I9yxS9q8qre4Csmpsj8HqRD2f+koKxgVSjkyYqEU6X0sLEZuHydL76TfXJdWMCUyEDkP5wnrFBhMvHRCdN3yA1j6d4nqBuRnfannkA/ZwU/VpkgiOHsnRfKzW2/9dBQ98Nr+JG0o9RxEZIdy/6fK/oTuS2LCweWmijAXq4gVlSmKWAgVIZm/R9LtznAO9nUKyKVdZ4F1uvquEOvPKl6VnennnLIlbH/SK4Q8sD9RpM9XUi4uS/uzLZjhHOzrFJArUfvzM2KhBnM7CnnPCgqm4r+Qk3VUvCQZupTtBOLfIzt7sy3mGVqWz3nl0gZc+7Kq7E9stimn00ko2puty+V0OpF1lgoxfRaRy/E9+bDuaSKgNN9J5DpUkf3JvxBR3v7k926xLZYmy30sPMQ8+5T32O3rK7VQtMVQYzkOzCHnRGe+h9lLv8M4TUuHJ2dD8e9z3/5kZ7hO7E8EWScz7VDJAcytFJWsDaUeL7sSRVgWZp22xdSlQuyLObI/2f7PidPmLBUidLhJMP4/pkIUYW5I2dmmo+IlF0KBTmapVv3uXGeu78mGm9XM6e8KkRzk09pITrbU7E/2ldee2J/o8ZKehaO3aFr0C8gNU/srqf5P9HiHdShMhSjC7KC5sdW1o/cmyBbmpEYPhB4QRZEZkWiHuWDKvphnWG1BJPyY+vMwsT/ZiJmT65S0P/dWopi/w51ZqqPVEorsT+jcldm8yH3x1v7cL+RlxEIJ5maUHThldxIUllU0S5DNT7o2GxIQdJTfdYixQZb5d5IDtbBsCduf7MzYiZdZKPvLJSsdARDVaQnXy00aAPlUiBP7k69DGR725yh6JSKXOJjLTRUhu4WpbC7VSVjOSedj/R7JDsTJw1zgAczNByLnu0IEvseJT+S1/Qnlm51EZjywP1V7Skg7iqI6s3wWkRPFa/sbPv2fGNL7KzEVohbZVIiqgin5sLDcucZ/FtEpK+dEp1zn7BXmU3SreM5r+3PrXQ1e2J9sODlVj/x21aI61dhCsduftWB0ePR/TuxP5jnjzpuKMDekmJGbPkvKxX8W6aAldaYdK+k9Owh5CukUkPL0XSGa9f+pn+3KS8kqkIv/TPvLdaR4KFoy8uDAmR4u9qcqLUv7yyQYP5qpEEWoKN50ax8B6fBjnnOwKy2k08HMxNrJ2hN083lI9epFHKFscqKy0i/Bi8nfF1mbV1Uw6oVO8bd+ZvshN9b7Ka/TLfvzPBUiojPPOdiVdqv/U2V/MhvR8e2mirDsYyEgNxAxJ5lFOpDUDk/EmREfcDPrtC0qvRFT+ooHOVnv4hECKLsvcr2sW8VzVj3yOl0ZTBQNYCK1JGmyLjhQqbJu2V9p67Qt6kn/Zz7daNSBkyl0nYmDGbFQg2VLb5udSHtPGFtbQ8bPTgrZZD1vZ1ESkUFB0lLTrlNOZz6VxVCDITlJVRaKLjX7c2P1QeqFOdm4TATZOgAt6w9iskWfCikC+yul/u9wsB8H2vukZGXsjztvKqLcVAY7ELXXQ29v67b8PCjQsw+m6Ej92a6srBwARASeyojJYx4Q8J5T2yQqqdOulLerQpKfRe5L6oMvImuOmInYX6qNy9qRNWKXR2fKsbKyInLhaNTys2zb2u0TgPisb9BkuxGR5yVqvZ8ig59Zp8h1RmJObMHcLzjoiwTayNxvObvOwvd/Tp6zVDtf19JpSy7V/mSe7cpyb4f2IeNY+Hw+VFeUAQD6B6N5jo5zYiCa8+ecsoMOZE3H2j1XxzpNsmERnSk6RBwLy/kWQ0giD+ZZQmpb5yK1AxG5L+b7L2QLDmy3z3SsyHWm2Z+krMh19g+ktK2sToH2CafeTwGd5rbVtPTvyoWsLfSlXJuYHUWMz+FB++fqyBZMOt26zj7ZPjftOZNvI7t6U20m9bpzoR+rj4VeMWQcCwCoqYw3pl0j73PwcKQaiWyHJ/ZAqtEpYqipOkRmJpZO1raU++g+j3lCJDJQp0aOZDs8sfsSsf4s6djKygHy1+nImRE63+QN9cKBAuy3kaZp6DMNuEK24MSxkOyL0vo/yT5FSGeKjtT+27ZOJxNIB46Q3Xsqa0NmWToWCqlKNKbdG5hqqLKzPsDBw+ygA5F1Spw4MyLpF7MeuwELb98VkkSkjSJRNbNq2UFTVKd0lGRQrqNM1SM2A7M6UGLnm5SNxDTbIWUnzkxan2LzfAeiMUu9gZNBXtYpSW1ru3KZzsG+TgdOmwvOjJNBXnbymXofBqP2bVdvk+pKOhbK0BvTbkeg3+iqingzOAmRyT7M/U7Cjx48WKl5zlyYH3zZF1a5iTkfLtK2qble2QHXUVrMlSiJGvtzy5mRfUZVpZlEZFNTPl6kSE8I9EWys3Egpf9z0raSg7yTyZyM82Wk6O3aQuJca0zOgX2nhI6FckRrLPSbMKo2YPnZDqoeZiezBNnIw0AkZrtWIi1iIZsKybsqxDvHI1Olt1DEIiYXPRiMxixOiaxDIiorn4pLSb9IOjNCMzAP0hKqUo4iss4cKDX2JxsBENGZeqyTtpV1ZgaisbQoY3adqSlHceerqbbS8nNenQkdDdUVxtbcdu8NUyEFoFqwxkI/Tr/xTsKPstEO1/LNijo8kYiFJRViU8bdVSEpG2RJhqJTIxainYCMTmcpAlnHQt7+ZFMEzgoT5c7XyX2RrT1wNlDL2UKaczpof5LjRY2ZE51pdmSzqFZFynHUCLHxRX+2qyvLDAfBdn/N4k31GDfBtmcYv4GNCcfCSYWyXW80FtMsVb/9gzHbbzJUmu+TDAuL1FiYdeaLSBRDpsSyKkSkeDPFsbDbCaQOtmKDplzleDSmpS39tB29GlRTSwIIOF8KV265WTyX7edspDqKbtReOVrxkGq7Av2fOf0r1P85eF7SI3z2nCgVttBYk3AsJKIO1ZXlAMRtl6kQheiOhd1qYf0mjBIMVZlls/1sVw6wvyTNSV5TNveb2payqRC7eLETbaZ3hYhFLOQGXEehXUUz42y/yyib0hm78byoTIXID7iFf85UXqcXDpTtgvmIuv7PiZOZWtOSVc6Bw5dMtYtFLMw1FtWVfiFZpkIKgO6l2Z8NxY3LnAqxm+sP9g1av8uBY1HoTlbTNPmO3eVVIV5iWRUiUryZWmNh12nLMNO0a3+6LeiOmP26ovSZmqhTIqrTPEvVZe0ORKqKBAGBGotBBw6UdPRAfi+UtEiSpGPRNxCxbX+qoiS6XjvI1hVl6v/spn1S5XoFauL08xVNtScXFZShpiIesRCN8DFioRC9Me177HEj0W+83Q1t+gYieHFTKwBgZE1F/LsEb3xVhR+BckFvNCHbUF0hJBeOxIyB3ThfSUcoKrkqZMeRbhzt7s96rNd+x6GuE9h9tMf4Wah4U3JVSDJUGr8nIhsq6Z1xU41YtE2fqVVXlCX3fREc5JMRPpsDgmmWKhodTI0qyqQIdJsXLZ5rFHxWgORgYjyjgtcp+mwDyedMv07RItX6qvjgFdPs776Zer72nbb4cYHyZP8naguiOjP1f7K28Oau42mTymzoOhoFbd68yVVV4hntZfGmd9QF4g9IT9jeje/uj3eM4+qrjJmU/rtc7Dft+X7FWePiOm3Ixc8tflxtZTlGJM63u9/u+caPG1sXsH2uZp1AciCy2+GZZQH7qZBYTENPysDz01d32JB0PxeiAfjxy9stv+sNR2znflMrzHvDYjOwxtpKIfszHzdxZHXiZ3s2FEocV1dVjrrEgBKybX+pOu2da+hE/Lhyvw+jR8RtN9Wu8umcMLJKSC4SjRmd8fh6XdbederHjasX02k+33H1YteZ1BmXE7G/7oQOo1+wrTMhl7hOQLwf08/Xri30huP3pDaQtD+7UYCelLa1a7f6ufp8yTay+4zqz9WFp4w2frflcDCv3GA0ZjgIkwSfF/06RwTKDafP/hgRl62rqrB1fKEYUo5FQ8Kr7LLpUerHjaqtRH3iRgRPDNiWO3lMLc6d0piQE9M5sqbCOF+7HrCuY+qoWiGd+nH1VeVGVMfuzKSrz9oedlMh3f2RtPRHe0/+tnUTcz1H6oMb0+x30IOJNtH357drf3oRZU1lmZD9mXVMHVUT/1nQFkbWVGBkdaXld3l1ptif7ecscU0jayqMQrZUu8oqq19nk5jOkKkT19vI7nOW2rbBvkFbKQJN09LaSLRfmJK4ThH7CybaUpcNCeocVVtpTMpkz9e+XMIWqitQXy3YX6e0reh11leZ7U/seTlrQj3OmTLSch65MJ/bFP0ZtWvzxjNaiZEC5xuLacn2raFjoQy9o7TbyeodXkNNhRFes/OABI2OstK4gfY7dpNstZisftxJCUMVfbBG1lQa4Ue7YUS9PfRQoN31310ZBshMhVs6XtdgZHK07A5EepuMSczGM117JvS0R2WZP2lHwk5mjeXnfBi2UF0p7IjrA9hJkjobqiuM67TvFMd1il9nXK4uUI5Rxn2RG8AGojFbqYnegaixwka2jcbVB4wN+2w/3yn9guhAbZ7kiPadyesUGzQbaiqM/k/0np4k7Ewn+nmT/YnawsiaCmO/Izvtq8vVVZVjtICc+biRNRVoqC63fb49AxFj91Z9PPOKIeVYGM6B6A00G5wdozHJNVTLzcDisoIRC30mNVqfvdnTaX6wAuXxiIXdXL4xq0l0znYjFpna0c5LjrxZFaJlPDe7DoKeHho9Qmw2pN+DQHlZ0sm0bbv6gJuYMQpGD8wdu93rTI9YCNq82REXfEZPSug8MRi15RRnGsDsO1Dx4yY0VKGizGdbVm+PynK/kUYRdRQbqpORJDuy4UjUCLlPFXQsgkb0IHlf5PsiMbl4n1sppjPRRlMkI7aWKJ3dvtPkiCedkvyywQxOm13bDZoifCLnq59rVYXfeL2FVwwpx0LkxkeiMSMfNbKm0hjkRbzRhhrzDMxe2NIqq0dYxDp282zITh7W7AEHErMhu2841R8QfdC0W2OR6aHPFbHwArMTkzFiYbPj0lMhuvNlV05PhVSW+9FQYz/aNmiqH9AHXNszMEln2nycbn+9A1Fbr0A3O7YNkumXSY3Vxi6EdmbyQZPN69dpNwKQjNIJ9gsZBk2pFJXEAObzAZOb5NJiohHbaEwzUk1Tm8QiM8nrNLWtzf7PiNI1yfV/DdWm9LNgZGakOcJiwxaCJp11gXIh280U4bNzT83RSK8ZUo6FSCdgzsHWV5ULdQSWcLLxQA7YysM+kCgStMrm16lpmunBig8mdvOwZkPVUyGyEYsnVrfYuk69Hcv9ydE713tRMm2r7SbmiEVlmf1aiVhMwwsbDgOQiVgkK+STs2qxwaS5qdqQs9PJhiwDmP1oW3wwics2N9UYTpmdQUGFM2MZ5O08oyeSs3ERObOseIrUPDPWByG7UZ0M5yswgNVXVRiF2aIpUmuUJP/5WuoHTM6MvX7B7GTav05N09JqXzTNXkGkxbFI6Hzkzb22ztcSSRJwSsz25/f7DL2dIg6q8IS3OOorgCHmWOhbp3b0DuTdTVBf+jiypgLlZX40JW7G8Z5wXj3HusOGPr0YaDCq5R3kzQ/t1FE1hmxHr73BRL+msaY8rB3vudNU0COSCtE0DR0JWb1+AAAWbz2SV7Yj0Y76UisA6LPhBHmQCYEG67p2/cFst2ELr2xpMz7r4e+OXnuDvO7MBCrKjCXPx20UuOod24hAuZH3jWnJ+5wL3XYbayuNJZxHu/NfZ3tPGJoG+H3x5dn6IHbMhuxR0/NiyNlo20g0ZlzTqBGVhi0dt6FTP6+m2mQBXLuN5wxIFhk31iSf7/be/Dr1vqOxphKNtXEb6rDpWOh9gLnA1U6/oF9To0muvTdsa9Ds6E1ep27zIjrrAuUYnVhlMRCJ2VpO2dGbHKh1m7dT1B3qjxhp2LF1VcaKuuM27ov5OqcnUjcAsPd4b065WExLytYm7c+OzZvtD4CxGsqWbMKOmmoqDTk7z2iqTi8ZUo7F2LoqlPt9iMS0nHsmAMChzhMAgMmN8RmfvoTucNeJvHoOdcWXm04aWY3qyjLjocwne7Az+ffrzpuC8Q1xo2kL5j5Xs+zoEZWoqigzHiw7y9kOJc5rQkO1KWKRvxNo7x3AQCQWD7Mm2gkANh3Kv9zqcOKa9KVhgP1O1i3Mr2jvMHVuenTmsI378r5p74vJiXD9QDRmy0HV0y+VZX5jSdohG/anL5UbEShHZbnfWEJnR1Y/ZtLIakxK3NNDnfnlDibkxtdXoaLMLybbGX9eJjfWCD1nbaF+RGMaKsv8GDMiINRG+nlNaqzGxAb7OgejMRwJxe/7xJFVmCBwvoe74nKTRlZjfEJnW7DflpN5qEvXWW0srT0czK+zNZh8tscmnrX+wZitmbHRL4yswoQGXWd+m9flJo6sxohAubGipFXgfCdadNpp2/gxjTUVqK4sS8raui/J67xixnjj9wfz2O7xnjAGoxr8PmBcnZj9HTbdTwDJ56WrL6sMEC+q1x2EiSOrTM9Zbrn4Mcln22uGlGNR5vdhfMLg8nV4B1NugkiHd9DUaZm/w67OOc0jUVnuNzqfVoGHWddVG9D3kLfhWJicqGSNRf6IhS43ti6AgKkYyE6o9ZDhCCUdi/7BWNacvJerQsKDMcsGNHpKw04HYq7NGIjEjD0TDtqQTUYskgP1QRsdiLEXSsIGRAZ5i2MhOVDr8rZlszgz+WbVus4JI6vg9/sc6zzWHc5bPNcW7EdMi9e9jK5NDib6QJFbZ2LC0ViNcXUBlPl9GIxqeaMz/YNRwxG13Bc799N0X6oqyjAm4WTamyAl+gVBJ/Ow4VhUGbrj32ejjTqTTomQ3aban4RTPGlkNcr8Pnz49DGJ883TX5uc6fIyvzG5smO7smOEPhZUVfjRVFtpyIX6I3n37TBsvpGOhXL0G5HPG23p0KMO8XxdsmPPLTcYjRkRBl2X3Q7PmLkljte97rZQ/gfy9kXrLedZWykesZg0stpIhfzi77vyvrr6sEnOvMw0ZCOvqT+UZscCyB+d8WJVSOpMS08R5LMFwJo66+wbtG1H/YNR/NeSnQDiNRZ27RawbjIE2Lf5aExL2m5jNSY3xm2/LdSfd8Bt6UxG6cz/P9CR2xHSNA0H2pMDrm7zJwajecPuLSnOv0gbtXQkO9nGmgpjJ8J8Trz5Ov1+n+37CVhnjOVl/qSTmUdWf85qEtFPMQfKOjOeaLONQv2Dlg3PjKiOjehBViczj05N0yx9yuRE33uo60TeqI5+Xvp56jq//fQm2w7q5BSnJJ8Tn36d8fPtDkfy1jyk9vX6d+R7Xsw25PP5UBsoN5b6H2jPLZv6jHrJkHMsThk7AgCwtTWU8zg9nH/GhDoAMHJvrcH+nHULO9q6EYlpqK8qNzqOaQnZnUe6c+rccjh+TqeOi5+jHl3p6hvMuRNm8MSgsY7+Q6fGvW09FaIPMtnoCUeMDmrKqBpLDcO7ezpyyu48Eg/znzSq1pI/zdeBxGKakSKY2FBl+duuo7nbyE10J2ZHyn3TI0m7j3TnrdXR2wiIR3J0+9uWx/7e3HXc+Oz3+XDy2BHw+eKz6nxpvGQqJD5YnjYubsP5dgTs7BtAJKbB54vXzIweEa+z0LT8z8uWQ/G/nzY+ruv08brO3HItHScQ6o+gssyP6WNqUVVRZhT8bWvNbQubE8+ormv6mHjbbm/LrbN/MIrdx+L35fRxdfD5fDjJ5jO6PXFOen8wbZQ9OSB5Pfp56kWGu/PYvG5DU0fVwufzGct5dx/tyTvg7kqcl75S5ySbOncldI6pC6A2UG60z8HOE3l35d1h6IzLTDF09mSVAeIObO9AFGV+n5HyqSz3IxyJYX+eAVdvf/08LzltjPG3fTkG3J5wxHBK9HY9daw929V16v17dWWZ4Zzkel7CkSjeT9if3h+cOb4eALA5j86trUGLHGB+1rI/35qmGdejH+8lUo7Fr3/9a0ybNg1VVVWYO3cu3nzzTdXnJc2cySMBAOsPdGU9JhyJGo7F2c3x40fWVBodwYaD2WXX7O+M62keCX9ixcOsyQ0AgI0Hc9/4tQc6LTrrAuWordRnUtkH642J85k0shrXnTcFAFCbGFR+s+z9rHIAcNMf10DT4tGRsXVVlsp43fjz6Z01ucEyo91yOJhzwN1zvAc94QiqKvzGErjkd2ZuIy/XhGxssZ7Th04djZrKMvQORHN2lpFozPKw33jhSYb9rUvc62zsPZ783p7+CEYEynFqojNZl8N2ozENtz+xHgBQlYg+zUnY09oDXTlnb+aixPIyP3w+nyG7Zl/u8914KH5O+vXpcpsOBnPW66zaF3dez5xQZ0TLZieel1zPGZBsQ/150eW2t3bnjLBsaOlCNKZhTF3AiJDMnqQ/o7l1rm/psuiclZA70NGHzhwRliOhfrSF+uH3ATMn1Vtk89Uk6e0wJ3F9p44bgUC5H939kZwDbiQaw+aE/c2ebD3fXH1R/O9xnXq7jK0LYExdANGYZgxumdA0LSmbprMrj8749546dgSqKspQUebHWRPibfWNv6y3JavrumLGeONzrmdt86EgNC0+wdGjp7rtrk/YSTY26DoT12mWzaVz6+EQBqMaGmsqDEdEt93dR3tyRup0+9P1mD+vzvGM7j3ei66+QVSW+XFGwonxEmHH4oknnsDtt9+O7373u1i3bh0+9KEPYcGCBThw4EAhzk+Y86c3AQBW7+/A0Swphle3HEHfQBTj66twypikZzhvalz279syr3rQNA3PrDsEALjItHf8zImJDq8t+wx3w8Eg9rf3obLcj3OnxrcB9/l8GJfo+J5NfG8mXt4cX3lwwfRRxu/0MOamQ0Hsb89c3dzVN4D39sY79nknxa/N3Dku33ksq871LV34+/ajAIBzpjRaZjG9A9GcD9bSHfHvnTN5ZNpGLW/tPp5JxMDnwbqQFXvaAQCfOmcS7r1mBi48eRTmJu7Rkiy2AADv7e1A/2AMI2sqsOOHV6K5qQbnTUvY377OnBXgehsByXcezE/c379tbM0q9+7eduOznrs/d0q8ZudAR1/OQez+l7YBgLEUEgAuPjVux0+vO5R1dtzdP2ikFvTZ0CljRmB8fRVODEaxZOvRrDpfSrys72LTDPP8RBu9srktqyO051gPNhwMwudL2v3kRColEtMs7ZfKc4nlvx88eRR8ibDUuVNHAgCW5bD5/sEo3tgRv5YPJM6xoabCmD0u35Vd9tXE6qCZkxpQk0hT6s7JW7uO53T43kg8Z/rrASrK/MagmesZ1e2vobrCiLDo206v3NOe0+HT20E/3ufz4ZzE+S7fmf0Z3d7WjSOhMCrL/IYDpV/n6v2d+NvGw1ll3048+7pOIGkX61q6sqZJj/eEjdm4WfaihO2+tKktkxgA03Umnmcgvj13fVU5uvoG8c77ma+1JxzBu4l+Ye6UpKz5Gc12T1/dEu8zzpvWZNjfqBEBzJgYb69sff2JgSiWJexaf0YA4OJElPq1rUeyvuNE7zPmndRovFrAS4TP4Gc/+xn++Z//GV/5yldw5pln4uc//zmam5vx0EMPFeL8hJk6qhZzpzYipgFffnSVJa3RFuzHX1a14L4XtgAAPv+BZiPqAABXzY5XDD+2Yr9lKdLxnjBe3dKGb/xlA9a3dCFQ7senzplk/L05MSsfiMTwo5e2GQYXi2nY396LFzYcxh2JWeZVsyYY74UAgD3H4nr++/Xd2GOKIAxEYth5pBtPrDqAv64+CAD49LlJnea6jEff3md81jQNx3vCWL2vA995ZpNRXPgfn5kNwLoyY+nOY9iaITR3sLMPX3t8DYB4Zz5nckPalsa/Xb4nbSDSNA0vbWrFD1+MD2AfnzUBV84Yj8vOGIt/+uA0+HzxiI/u7Jixu8W4SswuTJnfh+9ddSZuvPAk+Hw+XD1nIgDg/3trL9bs70i71mhMw2+W7wEALJg53piNTx8zAmc3j0QkpuFbf92ALYfTZ/Qr97TjnfeTDoLuWHx2bjMA4G8bD+NXb+zGziPdFtnO3gH8+JXki9z0jriuqgILZsZt95tPbsDSHUcta+1jMQ3PrDtoDFJ7TLZ99ZyJqKksw7bWEL6+aB1W7mk37kU0pmHlnnbMufc143i9/sTv9+Fz8yYDAO59YQte3dJmyA1EYlizvwPffWaT4Zx+zFSNv2DWBFSW+bHpUBC/e2uv0bZ6PcbTaw/isv9cBgD48OljjWW8Pp8P1yTuy8+X7LQsB+7qG8CqfR345pMb8Kd345Oc//OBZuPvHz1rPMr9Pmw+FDKcB51oTMOx7jC+/+xmdPdHMGlkNc47Kdmxf3zWBADAw8v3pKUJYok2+uXruwEA156dfEYvPm0Mqir82Nfel6ZTv96n1hzE9rZuVJb5ccWMccbfrkzczz+s2JcxNdE/GMWDb+xOnN94ox87u7kR4+oDCPVH8JNXdmR8rtbs7zRsYUHi2sw6f/H3XYazYyYSjeGXr+8CAFx6+hjDgZo+ZgROS6R3b/nTOqx4vz3teWnvCeO59XGnw2wL/3r5qZjSVANNA+5+fnNGZ/x3b+1FNKZh9uQGoy4IAD6d6IOXbDuC/3lzT1oK8Vh3GH9Z1RJvo5nJ66ws9+PahOz3n92M17cfQf9g1OIoPLR0N8KRGKaPqcWZE5KphatmTUB1RRm2t3Xj+89txsHOPsu1rjvQicdX7gdgtQUARrT5P1/bgZ8t3on97b2Gzr6BCH7w4lb0DkQxpanGcDKBuGM9dVQNgicGccV/LceGli5D50Akhpc3teLhRF/02bmT09rPC3yanQXPCQYGBlBTU4Mnn3wSn/rUp4zf33bbbVi/fj2WLVuWJhMOhxEOJ40lFAqhubkZwWAQ9fWFCdks3XEUX350FYD4Bk3nTmlEa+iEMesCgBkT6/HkTfONhwOIe6kz737V+Pm8k5pwrCectt75B9fOxA0XTLX87qRvv2h8PmtCPfz+eNjLvCnU2LoAnrvlg5jQkCyuufOvG/HE6rjx11SW4ZwpI9Ha1Y/9HX2W6McVZ43Db2+Ya3jAp3/vZcteFHOnNiIciWL/8T7Lfhp+H/DHfz4fH0xEWFbuaccXHl5p/D1Q7sf8k0cZu+l19Q2gpaPP2HP+qa/Nx9ypTdjWGsKCX7yJyY3VONx1AjENmD6mFqeMGYHBxC6mBzr6jPXWExuqsPiOS4wCQwD4xl824Km1B1Hu9+GcKSPRVFuJ/sEYTgxE8V4iZP7lC0/CPdfMgBuc+4PFRljyH86fgh99apbxt4FIDNc8+Ba2t8XzrFUVfoweETA22WkL9WPPsV5UlPnwyu0X42RT5GvN/g584eGVGDTtUtpQXYFRIypRW1mOba0hy9boH5sxDr+9YR4A4HvPbsLjK63Rv1G1lRhTF8CeY72WVSinjh2BxXdcAiBeAPfJB9+y7IMxubEaU0fVYM+xXkvRos8H7F14lfHzk6tb8K2/bjR+riz34/RxdWjp7EsrUtv3QFKuJxzBtb9620gXVVX4cerYuoRDZC0M3nX/AlSUJecxP1+yEz9fEh+kxtUHcPKYEdh1tCdtYHn6Xy60dLLHe8L46M+WobNvELWVZThnSiP2Hu9NK3T89DmT8LPPn2353d3PbcZjK/ajzO/DBdOb0D8YX1raFuw37ofPB/zmi3Mtg9/xnjA+/NOl6O6PYEJDFU4fX4dj3WG0Bfste2OcOnYEXrj1IkuU7od/24r/eWsvyvw+zJzUgDEjAugJD+JYdxjvH0v2KzddcjK+veAM4+fgiUFc9tOlaO8dQLnfh/OmNWFEoBzhSAzhSBSbD4XQE44gUO7Hq7dfbNQeAPGB+Ad/2wogXjw9Z3ID6qrK4ff7sOdYLzYe7EJMizskv75+riHXPxjFx36+3Hh78xnj6zC5MV5EOBCJYVtrCEe7w/D7gGf+5YOWcP2mg0Fc/eBbxs+Bcj80xGs+6qsq8P6xHnT2DeKM8XV48esfQplpQrd6Xwc+//BKS383rj6AqU21GIjGjPTAQ9efa3GEAODfn9uMP6zYb/ldU20l6qrKcTQUxonBKE4eU4tXbr/YYn/He8K4+pdvWZ6LQLkfExqq4Pf5DOf7v687x3Bmdf64Yh++/9wW4+eqCj9OHhOvkdqcqEWaN7URf/l/8y0T10g0hi8/usoStR1bF8C00bXYfChorEx7+Ia5lmWxAPDO+8fxD4+8a/w8qrYSp4wdga2HQ0Z/f/60Jvzp/15gaVvVhEIhNDQ05B2/hRyLw4cPY9KkSXj77bdx4YUXGr//0Y9+hMceeww7dqS/Fvuee+7Bvffem/b7QjoWg9EY5i98PW0vAb8POHNCPa6aPQFfvvAki1Oh84WHV2BlhqLG08aNwNypTbh6zgRcePLotL//+JXteGhper1DZbkfp40bgQ+eMhpfuWi6sRxM53hPGE+sasHLm1sNo9QZESjHyWNH4JJTR+NfPnyKpcP6z9d2GLOkMr/P8lD6fPHq6eljavHVi6cbBZ86+la4/+/xNRmjB0DcSBd+epZRiAbEZ4X1VRV4dv0h/PtzWzKuSAmU+3Ht2ZPwjStOs7yOGYiH1L/+53V4I0sYe+qoGvzqH87FzEQYuNB8/rcr8O7eDlx0ymg8/KW5afZwrDuMH720DS9uas24TLa2sgw/+dwcY0ZrZu2BTvzy77uwel9nxo3TrjhrHK6eMxGPr9yPn35ujhH1isU0/HXtQfxlVQu2tYYsy2CBuEP86XMn46VNrbj3mhmWtjrcdQK/emM3lu08lrYqoK6qHPOmNmJ/Rx+++/Ez8ZEzx1n+vr6lC4++vRfvvN9uGdzrqsoxdVQNNh8K4RsfPQ23fuRUi1zwxCB+vXQ3nl13CEdCSbnGmgrUVVXgQEcfLj9zHP7nxnkWuVhMw/+8tQe/fH23ZffEijIfRgTK0dk3iClNNVj6zUstnTMQr/H51pMb0wro6gLl6A5HUFNZhrfuvCxto6D+wSi+8ZcNeHFT9lTTL75wNj6ZMtME4g75bYvWWa7RzCfPnojvf+KstFVQJwai+Nr/rsmauqmuKMOX5k/Ftz52OsrLrAHk1fs6cNPja7JumjahoQo//sxsS5oJiEdgfvraDvz5vQNZVy9ccdY4/PT/zLFET4G4DS18eTte2tSaMa3bWFOBe66ZkbGNgn2DWPjyNjyz7lDGDfgmN1bj0S9/AKeOSy8uXL7zGH7y6o6MqTyfD7j9I6fhtstPTftbNKbh8ZX78diKfUb018zJY2rx8JfmWRx/naOhfjz4xu40xwSI7y3zrx89DV+79OS0vwHAkq1HcPsT69P6QJ8P+OScibjnmhnGxmxmwpEorv7lW5aib53mpmrceeUZ+MTsiWl/A4DbF63Ds+vTU01j6gL4zLmTcetlp1gmcoWgoI7FO++8g/nz5xu/v//++/HHP/4R27dvT5PxImIBxKuiO3oH0VBdgU2HghhbF8A5U0bmfU99Z++AMeveeLALY+urcPbkkcYe87noTyyhW/F+O+qrK3DymFpMaapJ6zAyEYtpeG9fB/a392LiyGqcOrYO4+oDRoQilYFIDKv2dWDu1Ea0Bvuxck87RtVWYtroWjQ31dh6CY1eULqjrQeBcj/qqytQX1WOaaNr05yCVLr6BvD27nZ09g0gUO5HTWU5JjXGnZnUziqV7W0h7GjrRqg/gqry+AtzJjdWY9akBlttpYpg3yBaQyfyFjtFojEc6jqB9t4BBPsGETwxCL/fh0tOG5P3LYKaFt+9r6N3AMd74v8/eWytrQIrfRvj1mA/jnT3Y3x9Fc6cYO+Z6eobwI62buzv6MO4+iqcP63Jtk3sOd6LzYeCmDSyGmc3j0R5mR+apmW1RV1ue1s3dh/twZkT6nHymPgKh7ZgPxqq45saZSIcieK9vR04GgqjuakGsyc3oKqiDLuPdmNEoMJYOZVKLKZhXUsXdrR146TRNZg5qQH1VRXY0dYNvw8ZBy+d9S1d2NEWQkN1BUaPCGBSYzV6+iMI9UeM2ppM9A1E8PbudhzvCWNcfQDj66vRH4mis3cAl50xNmv7aJqGNfs7cajrBLr7I6ivroAP8RTYtWdPyjkY6CnRba0hDERjqPD70R+JYs7kkZg1qSHN6UqVfW9vBw509KFvIIJwJIYpiTbWV0hkoy3Yj82HgkZfWO73YXJT3B4yTcjMRKLxVR7LdhxDU238jcoVZX586LTRRsowG939g3hxYysqy/0IlJfhxGAUF0xvsqRAsvHG9qMYjMYwdVQtQv3xnYrnTm20RCoysaGlC6H+QZw0qhaHu06gs28Q805qTHMSUznQ3ocDHX2Y3FiNXUd70DcQwfyTR2FsXe6+81h3GHuO9WBO80hsaOnCgY4+nDquDnMmN+R8xvoHo9hzrBenjB2BNfs7caCjFzMmNuDMCfUFjVKYKYhjIZMKkT0xQgghhBQPdsdvoelhZWUl5s6di8WLF1t+v3jxYktqhBBCCCHDE+GEzB133IEbbrgB8+bNw/z58/Hwww/jwIEDuOmmmwpxfoQQQggpIYQdi89//vNob2/Hfffdh9bWVsycORMvvfQSpk6dml+YEEIIIUMaoRoLFbDGghBCCCk9ClJjQQghhBCSCzoWhBBCCFEGHQtCCCGEKIOOBSGEEEKUQceCEEIIIcqgY0EIIYQQZdCxIIQQQogy6FgQQgghRBl0LAghhBCijMK+vD0D+kafoVDIbdWEEEIIkUQft/Nt2O26Y9Hd3Q0AaG5udls1IYQQQhzS3d2NhoaGrH93/V0hsVgMhw8fRl1dHXw+n7LvDYVCaG5uRktLC99BUgDYvoWF7Vt42MaFhe1bWIqhfTVNQ3d3NyZOnAi/P3slhesRC7/fj8mTJxfs++vr62nUBYTtW1jYvoWHbVxY2L6Fxev2zRWp0GHxJiGEEEKUQceCEEIIIcoYMo5FIBDA3XffjUAg4PWpDEnYvoWF7Vt42MaFhe1bWEqpfV0v3iSEEELI0GXIRCwIIYQQ4j10LAghhBCiDDoWhBBCCFEGHQtCCCGEKGPIOBa//vWvMW3aNFRVVWHu3Ll48803vT6lomfhwoX4wAc+gLq6OowdOxbXXnstduzYYTlG0zTcc889mDhxIqqrq3HppZdiy5YtlmPC4TBuvfVWjB49GrW1tbjmmmtw8OBBNy+lJFi4cCF8Ph9uv/1243dsX2ccOnQIX/ziFzFq1CjU1NTg7LPPxpo1a4y/s33liUQi+N73vodp06ahuroa06dPx3333YdYLGYcw/YVY/ny5bj66qsxceJE+Hw+PPvss5a/q2rPzs5O3HDDDWhoaEBDQwNuuOEGdHV1FfjqrBdS8ixatEirqKjQHnnkEW3r1q3abbfdptXW1mr79+/3+tSKmo997GPao48+qm3evFlbv369dtVVV2lTpkzRenp6jGMeeOABra6uTnvqqae0TZs2aZ///Oe1CRMmaKFQyDjmpptu0iZNmqQtXrxYW7t2rfbhD39YmzNnjhaJRLy4rKLkvffe00466SRt9uzZ2m233Wb8nu0rT0dHhzZ16lTty1/+svbuu+9qe/fu1ZYsWaLt3r3bOIbtK88Pf/hDbdSoUdrf/vY3be/evdqTTz6pjRgxQvv5z39uHMP2FeOll17Svvvd72pPPfWUBkB75plnLH9X1Z5XXnmlNnPmTO2dd97R3nnnHW3mzJnaJz7xCbcuUxsSjsV5552n3XTTTZbfnXHGGdq3v/1tj86oNDl69KgGQFu2bJmmaZoWi8W08ePHaw888IBxTH9/v9bQ0KD95je/0TRN07q6urSKigpt0aJFxjGHDh3S/H6/9sorr7h7AUVKd3e3duqpp2qLFy/WLrnkEsOxYPs6484779QuuuiirH9n+zrjqquu0v7pn/7J8rtPf/rT2he/+EVN09i+Tkl1LFS159atWzUA2sqVK41jVqxYoQHQtm/fXuCrilPyqZCBgQGsWbMGV1xxheX3V1xxBd555x2Pzqo0CQaDAICmpiYAwN69e9HW1mZp20AggEsuucRo2zVr1mBwcNByzMSJEzFz5ky2f4Kbb74ZV111FS6//HLL79m+znj++ecxb948fO5zn8PYsWNxzjnn4JFHHjH+zvZ1xkUXXYS///3v2LlzJwBgw4YNeOutt/Dxj38cANtXNarac8WKFWhoaMD5559vHHPBBRegoaHBtTZ3/SVkqjl+/Dii0SjGjRtn+f24cePQ1tbm0VmVHpqm4Y477sBFF12EmTNnAoDRfpnadv/+/cYxlZWVaGxsTDuG7Q8sWrQIa9euxapVq9L+xvZ1xp49e/DQQw/hjjvuwHe+8x289957+PrXv45AIIAvfelLbF+H3HnnnQgGgzjjjDNQVlaGaDSK+++/H9dddx0A2q9qVLVnW1sbxo4dm/b9Y8eOda3NS96x0El9BbumaUpfyz7UueWWW7Bx40a89dZbaX+TaVu2P9DS0oLbbrsNr732GqqqqrIex/aVIxaLYd68efjRj34EADjnnHOwZcsWPPTQQ/jSl75kHMf2leOJJ57A448/jj/96U+YMWMG1q9fj9tvvx0TJ07EjTfeaBzH9lWLivbMdLybbV7yqZDRo0ejrKwszRM7evRomudHMnPrrbfi+eefxxtvvGF5pf348eMBIGfbjh8/HgMDA+js7Mx6zHBlzZo1OHr0KObOnYvy8nKUl5dj2bJl+O///m+Ul5cb7cP2lWPChAk466yzLL8788wzceDAAQC0X6d861vfwre//W184QtfwKxZs3DDDTfgX//1X7Fw4UIAbF/VqGrP8ePH48iRI2nff+zYMdfavOQdi8rKSsydOxeLFy+2/H7x4sW48MILPTqr0kDTNNxyyy14+umn8frrr2PatGmWv0+bNg3jx4+3tO3AwACWLVtmtO3cuXNRUVFhOaa1tRWbN28e9u3/kY98BJs2bcL69euNf/PmzcP111+P9evXY/r06WxfB3zwgx9MWx69c+dOTJ06FQDt1yl9fX3w+61DRFlZmbHclO2rFlXtOX/+fASDQbz33nvGMe+++y6CwaB7be5KiWiB0Zeb/u53v9O2bt2q3X777Vptba22b98+r0+tqPna176mNTQ0aEuXLtVaW1uNf319fcYxDzzwgNbQ0KA9/fTT2qZNm7Trrrsu4/KnyZMna0uWLNHWrl2rXXbZZcN2OVk+zKtCNI3t64T33ntPKy8v1+6//35t165d2v/+7/9qNTU12uOPP24cw/aV58Ybb9QmTZpkLDd9+umntdGjR2v/9m//ZhzD9hWju7tbW7dunbZu3ToNgPazn/1MW7dunbE1gqr2vPLKK7XZs2drK1as0FasWKHNmjWLy01l+NWvfqVNnTpVq6ys1M4991xjySTJDoCM/x599FHjmFgspt19993a+PHjtUAgoF188cXapk2bLN9z4sQJ7ZZbbtGampq06upq7ROf+IR24MABl6+mNEh1LNi+znjhhRe0mTNnaoFAQDvjjDO0hx9+2PJ3tq88oVBIu+2227QpU6ZoVVVV2vTp07Xvfve7WjgcNo5h+4rxxhtvZOxzb7zxRk3T1LVne3u7dv3112t1dXVaXV2ddv3112udnZ0uXaWm8bXphBBCCFFGyddYEEIIIaR4oGNBCCGEEGXQsSCEEEKIMuhYEEIIIUQZdCwIIYQQogw6FoQQQghRBh0LQgghhCiDjgUhhBBClEHHghBCCCHKoGNBCCGEEGXQsSCEEEKIMuhYEEIIIUQZ/z88mcsNfTLDngAAAABJRU5ErkJggg==", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "for pot in data_1:\n", + " plt.plot(np.abs(pot.flatten()))" + ] + }, + { + "cell_type": "code", + "execution_count": 40, + "id": "ea137107-de11-42f1-acd7-31a34744cc37", + "metadata": {}, + "outputs": [ + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAGdCAYAAAC7EMwUAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAG21JREFUeJzt3X9slfXd//HX4ddl0dMTG+z5MY5No+imIMnEQTuUwkJj74yAbIk/cpsSNyMKJE01ODS5afYHRYyNJp1scwuTTFf+mIiJKHTBlpmOpRCIDRiDsc7utmcVgueUiqcUPvcffjlfj0Dpac/h3XP6fCRXYq/r03M+1z4dz1w951z1OeecAAAwMMl6AgCAiYsIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM1OsJ/Bd58+f1+effy6/3y+fz2c9HQBAhpxz6u/vVyQS0aRJw1/rjLsIff7554pGo9bTAACMUU9Pj2bOnDnsmJxF6OWXX9bzzz+v3t5e3X777XrxxRd19913X/H7/H6/JGmh/ktTNDVX0wMA5MiQzup97U79ez6cnERox44dqqur08svv6wf//jH+t3vfqeamhodO3ZMN95447Dfe+FXcFM0VVN8RAgA8s7/uyPpSF5SyckbE5qamvSLX/xCv/zlL/WDH/xAL774oqLRqLZu3ZqLpwMA5KmsR2hwcFCHDh1SdXV12v7q6mp1dHRcND6ZTCqRSKRtAICJIesROnHihM6dO6dgMJi2PxgMKhaLXTS+sbFRgUAgtfGmBACYOHL2OaHv/i7QOXfJ3w9u2LBB8Xg8tfX09ORqSgCAcSbrb0yYMWOGJk+efNFVT19f30VXR5LkeZ48z8v2NAAAeSDrV0LTpk3TnXfeqdbW1rT9ra2tqqyszPbTAQDyWE7eol1fX6+HH35Y8+bNU0VFhX7/+9/rs88+0+rVq3PxdACAPJWTCN1///06efKkfv3rX6u3t1ezZ8/W7t27VVZWlounAwDkKZ9zzllP4tsSiYQCgYCqtJwPqwJAHhpyZ9WmXYrH4youLh52LHfRBgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzObl3HArcCP5ufJrxdWcoAOMIV0IAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmplhPAPnHN3lyRuPd0FCOZgIg33ElBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAz3jkPmMrx3nLh3HIDL4EoIAGAm6xFqaGiQz+dL20KhULafBgBQAHLy67jbb79df/vb31JfT8701zcAgAkhJxGaMmUKVz8AgCvKyWtCx48fVyQSUXl5uR544AF98sknlx2bTCaVSCTSNgDAxJD1CM2fP1/bt2/Xnj179MorrygWi6myslInT5685PjGxkYFAoHUFo1Gsz0lAMA45XPOuVw+wcDAgG666SatX79e9fX1Fx1PJpNKJpOprxOJhKLRqKq0XFN8U3M5NYySz/MyGu++tb4ACt+QO6s27VI8HldxcfGwY3P+OaFrr71Wc+bM0fHjxy953PM8eRn+owYAKAw5/5xQMpnUhx9+qHA4nOunAgDkmaxH6KmnnlJ7e7u6u7v1z3/+Uz//+c+VSCRUW1ub7acCAOS5rP867t///rcefPBBnThxQjfccIMWLFigAwcOqKysLNtPBSOTMvz16TleEwJwGVmPUEtLS7YfEgBQoLh3HADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYyfmfckDh8U0vyuwb+Gu5AC6DKyEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMDNub9vzP12HdJ3fvpHn5Bvx2LNuckaPfdaN/H/+xPlrMnrsL4aKRzy26dhPMnrsTB2t2DPisSfODWT02P3n3YjHfu0y+3k6n8HaZ/Jzks/Ou4lxnrk0yTfyn1lJOpvBz+1kZfbYuXK6/7wWzRnZWPt/5QEAExYRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAz4/becZPlRnwfpEzvxZSJazK4F9Mk39mMHntyBvMuUyKjx/Zf88WIx/6iIpbRY9/W8d85G79+9t6MHrtk8ukRjy2e9HVGjz3VN5TB2HMZPXam9/jK5T3Bcvn/n1zKZN75fM+78XI/uFzhSggAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAICZcXvvuI1z5mmKb6r1NHAJ5aGR35dOkoZi/xnx2L8okul0AIwzQ+6spP8d0ViuhAAAZjKO0P79+7Vs2TJFIhH5fD69+eabacedc2poaFAkElFRUZGqqqp09OjRbM0XAFBAMo7QwMCA5s6dq+bm5kse37Jli5qamtTc3KzOzk6FQiEtXbpU/f39Y54sAKCwZPyaUE1NjWpqai55zDmnF198Uc8++6xWrlwpSXr11VcVDAb1+uuv67HHHhvbbAEABSWrrwl1d3crFoupuro6tc/zPC1atEgdHR2X/J5kMqlEIpG2AQAmhqxGKBb75i90BoPBtP3BYDB17LsaGxsVCARSWzQazeaUAADjWE7eHefzpf8pXefcRfsu2LBhg+LxeGrr6enJxZQAAONQVj8nFAqFJH1zRRQOh1P7+/r6Lro6usDzPHmel81pAADyRFavhMrLyxUKhdTa2praNzg4qPb2dlVWVmbzqQAABSDjK6HTp0/r448/Tn3d3d2tI0eOqKSkRDfeeKPq6uq0adMmzZo1S7NmzdKmTZs0ffp0PfTQQ1mdOAAg/2UcoYMHD2rx4sWpr+vr6yVJtbW1+tOf/qT169frzJkzeuKJJ3Tq1CnNnz9fe/fuld/vz96sYcolB62nAKBA+JxzznoS35ZIJBQIBFSl5dw7bpyafP31GY0/d+pUjmYCYDwacmfVpl2Kx+MqLi4ediz3jgMAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM1n9Uw6YGNyZM9ZTAFAguBICAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADPctgcZc0ND1lMAUCC4EgIAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGe8chY+68s54CgALBlRAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmOG2PcicO289AwAFgishAIAZIgQAMJNxhPbv369ly5YpEonI5/PpzTffTDu+atUq+Xy+tG3BggXZmi8AoIBkHKGBgQHNnTtXzc3Nlx1z7733qre3N7Xt3r17TJMEABSmjN+YUFNTo5qammHHeJ6nUCg06kkBACaGnLwm1NbWptLSUt1yyy169NFH1dfXd9mxyWRSiUQibQMATAxZj1BNTY1ee+017du3Ty+88II6Ozu1ZMkSJZPJS45vbGxUIBBIbdFoNNtTAgCMUz7n3Kj/VrPP59POnTu1YsWKy47p7e1VWVmZWlpatHLlyouOJ5PJtEAlEglFo1FVabmm+KaOdmrIJZ8vs/Gj/xEDkIeG3Fm1aZfi8biKi4uHHZvzD6uGw2GVlZXp+PHjlzzueZ48z8v1NAAA41DOPyd08uRJ9fT0KBwO5/qpAAB5JuMrodOnT+vjjz9Ofd3d3a0jR46opKREJSUlamho0M9+9jOFw2F9+umneuaZZzRjxgzdd999WZ04ACD/ZRyhgwcPavHixamv6+vrJUm1tbXaunWrurq6tH37dn355ZcKh8NavHixduzYIb/fn71Zwxav8QDIkowjVFVVpeHey7Bnz54xTQgAMHFw7zgAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzGUWosbFRd911l/x+v0pLS7VixQp99NFHaWOcc2poaFAkElFRUZGqqqp09OjRrE4aAFAYMopQe3u71qxZowMHDqi1tVVDQ0Oqrq7WwMBAasyWLVvU1NSk5uZmdXZ2KhQKaenSperv78/65AEA+c3nnHOj/eYvvvhCpaWlam9v1z333CPnnCKRiOrq6vT0009LkpLJpILBoJ577jk99thjV3zMRCKhQCCgKi3XFN/U0U4NAGBkyJ1Vm3YpHo+ruLh42LFjek0oHo9LkkpKSiRJ3d3disViqq6uTo3xPE+LFi1SR0fHJR8jmUwqkUikbQCAiWHUEXLOqb6+XgsXLtTs2bMlSbFYTJIUDAbTxgaDwdSx72psbFQgEEht0Wh0tFMCAOSZUUdo7dq1+uCDD/SXv/zlomM+ny/ta+fcRfsu2LBhg+LxeGrr6ekZ7ZQAAHlmymi+ad26dXrrrbe0f/9+zZw5M7U/FApJ+uaKKBwOp/b39fVddHV0ged58jxvNNMAAOS5jK6EnHNau3at3njjDe3bt0/l5eVpx8vLyxUKhdTa2praNzg4qPb2dlVWVmZnxgCAgpHRldCaNWv0+uuva9euXfL7/anXeQKBgIqKiuTz+VRXV6dNmzZp1qxZmjVrljZt2qTp06froYceyskJAADyV0YR2rp1qySpqqoqbf+2bdu0atUqSdL69et15swZPfHEEzp16pTmz5+vvXv3yu/3Z2XCAIDCMabPCeUCnxMCgPx21T4nBADAWBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAMxlFqLGxUXfddZf8fr9KS0u1YsUKffTRR2ljVq1aJZ/Pl7YtWLAgq5MGABSGjCLU3t6uNWvW6MCBA2ptbdXQ0JCqq6s1MDCQNu7ee+9Vb29vatu9e3dWJw0AKAxTMhn87rvvpn29bds2lZaW6tChQ7rnnntS+z3PUygUys4MAQAFa0yvCcXjcUlSSUlJ2v62tjaVlpbqlltu0aOPPqq+vr7LPkYymVQikUjbAAATw6gj5JxTfX29Fi5cqNmzZ6f219TU6LXXXtO+ffv0wgsvqLOzU0uWLFEymbzk4zQ2NioQCKS2aDQ62ikBAPKMzznnRvONa9as0dtvv633339fM2fOvOy43t5elZWVqaWlRStXrrzoeDKZTAtUIpFQNBpVlZZrim/qaKYGADA05M6qTbsUj8dVXFw87NiMXhO6YN26dXrrrbe0f//+YQMkSeFwWGVlZTp+/Pglj3ueJ8/zRjMNAECeyyhCzjmtW7dOO3fuVFtbm8rLy6/4PSdPnlRPT4/C4fCoJwkAKEwZvSa0Zs0a/fnPf9brr78uv9+vWCymWCymM2fOSJJOnz6tp556Sv/4xz/06aefqq2tTcuWLdOMGTN033335eQEAAD5K6Mroa1bt0qSqqqq0vZv27ZNq1at0uTJk9XV1aXt27fryy+/VDgc1uLFi7Vjxw75/f6sTRoAUBgy/nXccIqKirRnz54xTQgAMHFw7zgAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzBAhAIAZIgQAMEOEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADBDhAAAZogQAMAMEQIAmCFCAAAzRAgAYCajCG3dulV33HGHiouLVVxcrIqKCr3zzjup4845NTQ0KBKJqKioSFVVVTp69GjWJw0AKAwZRWjmzJnavHmzDh48qIMHD2rJkiVavnx5KjRbtmxRU1OTmpub1dnZqVAopKVLl6q/vz8nkwcA5Defc86N5QFKSkr0/PPP65FHHlEkElFdXZ2efvppSVIymVQwGNRzzz2nxx57bESPl0gkFAgEVKXlmuKbOpapAQAMDLmzatMuxeNxFRcXDzt21K8JnTt3Ti0tLRoYGFBFRYW6u7sVi8VUXV2dGuN5nhYtWqSOjo7LPk4ymVQikUjbAAATQ8YR6urq0nXXXSfP87R69Wrt3LlTt912m2KxmCQpGAymjQ8Gg6ljl9LY2KhAIJDaotFoplMCAOSpjCN066236siRIzpw4IAef/xx1dbW6tixY6njPp8vbbxz7qJ937ZhwwbF4/HU1tPTk+mUAAB5akqm3zBt2jTdfPPNkqR58+aps7NTL730Uup1oFgspnA4nBrf19d30dXRt3meJ8/zMp0GAKAAjPlzQs45JZNJlZeXKxQKqbW1NXVscHBQ7e3tqqysHOvTAAAKUEZXQs8884xqamoUjUbV39+vlpYWtbW16d1335XP51NdXZ02bdqkWbNmadasWdq0aZOmT5+uhx56KFfzBwDksYwi9J///EcPP/ywent7FQgEdMcdd+jdd9/V0qVLJUnr16/XmTNn9MQTT+jUqVOaP3++9u7dK7/fn5PJAwDy25g/J5RtfE4IAPLbVfmcEAAAY0WEAABmiBAAwAwRAgCYIUIAADNECABghggBAMwQIQCAGSIEADCT8V20c+3CDRyGdFYaV/dyAACMxJDOSvr//54PZ9xFqL+/X5L0vnYbzwQAMBb9/f0KBALDjhl39447f/68Pv/8c/n9/rQ/hpdIJBSNRtXT03PFexHlM86zcEyEc5Q4z0KTjfN0zqm/v1+RSESTJg3/qs+4uxKaNGmSZs6cednjxcXFBf0DcAHnWTgmwjlKnGehGet5XukK6ALemAAAMEOEAABm8iZCnudp48aN8jzPeio5xXkWjolwjhLnWWiu9nmOuzcmAAAmjry5EgIAFB4iBAAwQ4QAAGaIEADATN5E6OWXX1Z5ebmuueYa3Xnnnfr73/9uPaWsamhokM/nS9tCoZD1tMZk//79WrZsmSKRiHw+n95888204845NTQ0KBKJqKioSFVVVTp69KjNZMfgSue5atWqi9Z2wYIFNpMdpcbGRt11113y+/0qLS3VihUr9NFHH6WNKYT1HMl5FsJ6bt26VXfccUfqA6kVFRV65513Usev5lrmRYR27Nihuro6Pfvsszp8+LDuvvtu1dTU6LPPPrOeWlbdfvvt6u3tTW1dXV3WUxqTgYEBzZ07V83NzZc8vmXLFjU1Nam5uVmdnZ0KhUJaunRp6v6B+eJK5ylJ9957b9ra7t6dX/dGbG9v15o1a3TgwAG1trZqaGhI1dXVGhgYSI0phPUcyXlK+b+eM2fO1ObNm3Xw4EEdPHhQS5Ys0fLly1Ohuapr6fLAj370I7d69eq0fd///vfdr371K6MZZd/GjRvd3LlzraeRM5Lczp07U1+fP3/ehUIht3nz5tS+r7/+2gUCAffb3/7WYIbZ8d3zdM652tpat3z5cpP55EpfX5+T5Nrb251zhbue3z1P5wpzPZ1z7vrrr3d/+MMfrvpajvsrocHBQR06dEjV1dVp+6urq9XR0WE0q9w4fvy4IpGIysvL9cADD+iTTz6xnlLOdHd3KxaLpa2r53latGhRwa2rJLW1tam0tFS33HKLHn30UfX19VlPaUzi8bgkqaSkRFLhrud3z/OCQlrPc+fOqaWlRQMDA6qoqLjqaznuI3TixAmdO3dOwWAwbX8wGFQsFjOaVfbNnz9f27dv1549e/TKK68oFoupsrJSJ0+etJ5aTlxYu0JfV0mqqanRa6+9pn379umFF15QZ2enlixZomQyaT21UXHOqb6+XgsXLtTs2bMlFeZ6Xuo8pcJZz66uLl133XXyPE+rV6/Wzp07ddttt131tRx3d9G+nG//WQfpmx+Q7+7LZzU1Nan/njNnjioqKnTTTTfp1VdfVX19veHMcqvQ11WS7r///tR/z549W/PmzVNZWZnefvttrVy50nBmo7N27Vp98MEHev/99y86VkjrebnzLJT1vPXWW3XkyBF9+eWX+utf/6ra2lq1t7enjl+ttRz3V0IzZszQ5MmTLypwX1/fRaUuJNdee63mzJmj48ePW08lJy6882+iraskhcNhlZWV5eXarlu3Tm+99Zbee++9tD+5UmjrebnzvJR8Xc9p06bp5ptv1rx589TY2Ki5c+fqpZdeuuprOe4jNG3aNN15551qbW1N29/a2qrKykqjWeVeMpnUhx9+qHA4bD2VnCgvL1coFEpb18HBQbW3txf0ukrSyZMn1dPTk1dr65zT2rVr9cYbb2jfvn0qLy9PO14o63ml87yUfFzPS3HOKZlMXv21zPpbHXKgpaXFTZ061f3xj390x44dc3V1de7aa691n376qfXUsubJJ590bW1t7pNPPnEHDhxwP/3pT53f78/rc+zv73eHDx92hw8fdpJcU1OTO3z4sPvXv/7lnHNu8+bNLhAIuDfeeMN1dXW5Bx980IXDYZdIJIxnnpnhzrO/v989+eSTrqOjw3V3d7v33nvPVVRUuO9973t5dZ6PP/64CwQCrq2tzfX29qa2r776KjWmENbzSudZKOu5YcMGt3//ftfd3e0++OAD98wzz7hJkya5vXv3Oueu7lrmRYScc+43v/mNKysrc9OmTXM//OEP094yWQjuv/9+Fw6H3dSpU10kEnErV650R48etZ7WmLz33ntO0kVbbW2tc+6bt/Vu3LjRhUIh53meu+eee1xXV5ftpEdhuPP86quvXHV1tbvhhhvc1KlT3Y033uhqa2vdZ599Zj3tjFzq/CS5bdu2pcYUwnpe6TwLZT0feeSR1L+nN9xwg/vJT36SCpBzV3ct+VMOAAAz4/41IQBA4SJCAAAzRAgAYIYIAQDMECEAgBkiBAAwQ4QAAGaIEADADBECAJghQgAAM0QIAGCGCAEAzPwfcqKkDt4laM0AAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "data = data_1-data_0\n", + "for pot in data:\n", + " plt.imshow(np.abs(pot))\n", + " plt.show()" + ] + }, + { + "cell_type": "code", + "execution_count": 15, + "id": "50c63c37-9e0d-4d0a-b7d4-729b78f67d62", + "metadata": {}, + "outputs": [], + "source": [ + "d0 = np.diff(data_0.flatten(), append=data_0.flatten()[-1])" + ] + }, + { + "cell_type": "code", + "execution_count": 16, + "id": "769f249b-be6f-47c0-828c-c39bbad90b3d", + "metadata": {}, + "outputs": [], + "source": [ + "d1 = np.diff(data_1.flatten(), append=data_1.flatten()[-1])" + ] + }, + { + "cell_type": "code", + "execution_count": 17, + "id": "c057f22c-949c-4f86-9fb0-bf88a9028ef2", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "(0.0, 35.0)" + ] + }, + "execution_count": 17, + "metadata": {}, + "output_type": "execute_result" + }, + { + "data": { + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGdCAYAAAAxCSikAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAATC1JREFUeJzt3XuYHOV9J/pvVfX9rrnPSCNpJIEQEkggYRAx941sccKamONgnz2BOAlZ2WAvlolj2ZvFcZwVGxMWc8zFxGDM8tgmOcIsWUiCvIDkmKtAgCwECHQZSTOj0Yw0fb/Xu39UVUsjzUjTM11dVd3fz/PMo5me7qp3ehr6O/W+v98rCSEEiIiIiGxMtnoARERERGfCwEJERES2x8BCREREtsfAQkRERLbHwEJERES2x8BCREREtsfAQkRERLbHwEJERES257J6AFOhqioGBgYQDochSZLVwyEiIqIpEEIgmUyip6cHsjyzaySOCCwDAwPo7e21ehhEREQ0DQcOHMCcOXNmdAxHBJZwOAxA+4EjkYjFoyEiIqKpSCQS6O3trbyPz4QjAosxDRSJRBhYiIiIHKYWyzm46JaIiIhsj4GFiIiIbI+BhYiIiGyPgYWIiIhsj4GFiIiIbI+BhYiIiGyPgYWIiIhsj4GFiIiIbI+BhYiIiGyPgYWIiIhsj4GFiIiIbI+BhYiIiGyPgYWIiKgJCFXF/g/ehlpWrR7KtDCwEBERNYHX//H7mPfzK/DGP/43q4cyLQwsREREzWD4Pe3fI7usHcc0OSqwqOWy1UMgIiJyJKWYAgC4CgmLRzI9jgosqWTc6iEQERE5kruY0P9NWjyS6XFUYMkkRq0eAhERkSN5StoVFl+JV1hMl0kctXoIREREjuQv64FF/9dpHBVY8qkxq4dARETkSH41DQAICQYW0+VTx6weAhERkSOFRLryr1Cd14vFUYGllB2zeghERESOUyoWEJRyAACPVEY247yFt44KLCLDKyxERETVSifGv38m484rYnFUYEHOmSubiYiIrJSKjy9ayTCwmEvKsQ8LERFRtbLJ8YEl68CqW0cFFtmhzW6IiIislEuNDyiFFAOLqZzanY+IiMhKxfT4GYqiA6tuHRVYPCUGFiIiomoV0+MDSjnLwGIqp3bnIyIispKaHX+FRWSdtybUUYElIDJWD4GIiMhxxElFK1JuzJqBzICjAotT2wkTERFZScprbUFKQnvbVwrOaxPiqMDil4rI53iVhYiIqBqyHliG5TYAgKvAKSHTJcec1+yGiIjISi69yvaYuwuAM4tYHBdYMgkGFiIiomoYASUTmA0A8DOwmI+BhYiIqDo+PaCUIr0AgICatnI40+K4wFJIOq92nIiIyEo+PaC4WuYBcGYRi/MCS9p57YSJiIisFBRaYAl0LAAAhKQsSsWClUOqmuMCSzkzZvUQiIiIHEOoKkJ6H7NZPQsrt5+8g7PdOS6wCAYWIiKiKctmknBLZQBAtLUTaeEDAKQdtibUcYEFDuzOR0REZJV0Qlv7WRIyAsEIUlIQAJCJM7CYSs47r9kNERGRVYxgkpICkGQZGTkMAMg5rIjFcYHFVXReO2EiIiKrZFNaMEnrV1aySggAUEhxDYupPAwsREREU5bXA0tW1oJK3h0BAJQzvMJiKm/JebXjREREVimmxwAAef3KSsmtTQmpDCzmCqjOaydMRERkFaMdSMGlBZayNwoAEDlnrQl1XGAxmt8QERHRmanZMQBAyaNNBQk9sDitiMVxgSUs0hBq2ephEBEROUNeW/tZ9mhTQZI/BgBwFZy1JtRxgUWRBFJJZ6VCIiIiq8h6YDGurCiBGADA7bAiFkcFlrxwAQDS8RGLR0JEROQMSkFb+yn5tCkhd7AFAOAtOWtNqKMCi9GdL+2w/Q+IiIisYlxJkfUrK57QLABAoMzAYhqj6U0+6ax2wkRERFbxlrViFVdACyr+iHaFJeCwIhZHBZa80fQm5azacSIiIqv49CspnqC2hsUfaQUAhEUKQlUtG1e1HBVYcnoNeSnNKSEiIqKpCKjalRRfWLuyEo61AQA8Uhm5rHOusjgqsBT17nxGExwiIiI6vZA+9eMP62tXghGUhPb2nxxzThGLowKLUUMOvQkOERERTa5ULCAg5QEAQX0qSJJlJCVtxiKTcM6aUEcFFuHVSrIkh3XnIyIiskLqhKraYGRW5fPKzs0MLCbxxQAAssO68xEREVkhrQeSjPDC7fFWbs8o2oyFk4pYHBVYZL+2wtnNwEJERHRG2aQWSIw+ZgZj5+ZiyjlFLI4KLC696Y3TuvMRERFZIZfUAklGHh9YCm5tiYWTilgcFVg8wRgAwO+w7nxERERWKKS1NZ85/YqKobJzs4OKWBwVWLxGO2E1ZfFIiIiI7K+U0aaE8icFFtUbA+CsIhZHBRZ/WCvJCjmsnTAREZEV1KwWSIr6FFCFvhGi0qiBZePGjbjooosQDofR0dGB66+/Hh988MEZH7dlyxasXLkSPp8PCxYswEMPPTStwQaj+hUWKY9CPjetYxARETULoQeWSh8znazvK+RyUBFLVYFly5YtuPXWW/Hqq69i8+bNKJVKWLNmDdLpya947N27F9deey0uu+wybN++Hd/61rfw1a9+FZs2bap6sCF9wybAWd35iIiIrCDltUCiesdfYTGKWDwl5wQWVzV3/pd/+ZdxX//kJz9BR0cH3nzzTVx++eUTPuahhx7C3Llzce+99wIAlixZgm3btuHuu+/GDTfcUNVgFZcLSeFHWMoikxhFa+ecqh5PRETUTBTjCspJgcUT0i4A+MrOWRM6ozUs8bh2qamlpWXS+7zyyitYs2bNuNs+9alPYdu2bSgWixM+Jp/PI5FIjPswpIx2wnHndOcjIiKygquoVdXK/ti42736RojBZggsQgisX78en/zkJ7Fs2bJJ7zc0NITOzs5xt3V2dqJUKmFkZOJpnY0bNyIajVY+ent7K9/L6CudndSdj4iIyAoePbAogei42/36vkJB0QSB5bbbbsO7776Ln//852e8ryRJ474WQkx4u2HDhg2Ix+OVjwMHDlS+l9PbCRfSzunOR0REZAWvfgXFHZw17vagviY0LGVRLpXqPq7pqGoNi+ErX/kKnnnmGWzduhVz5px+HUlXVxeGhobG3TY8PAyXy4XW1tYJH+P1euH1eif8XsEVBgpAOT02naETERE1Db/et8yrN141hGNtlc9T8VFEW8fPhNhRVVdYhBC47bbb8NRTT+GFF15AX1/fGR+zevVqbN68edxtzz//PFatWgW3213daAEUPdplLTXLKSEiIqLTCep9y/zh8WtN3R4vMkK7MJCKO6PqtqrAcuutt+KJJ57Az372M4TDYQwNDWFoaAjZbLZynw0bNuCmm26qfL1u3Trs378f69evx65du/Doo4/ikUcewR133DGtAVdKs3LOaXZDRERUb0JVERIZAIA/cmpxTNIoYkk4Y4lFVYHlwQcfRDwex5VXXonu7u7Kx5NPPlm5z+DgIPr7+ytf9/X14bnnnsNLL72EFStW4K//+q9x3333VV3SbBC+GABndecjIiKqt0w6AZekAgBC0VMDS0bWAkvOIYGlqjUsxmLZ03nsscdOue2KK67AW2+9Vc2pJiX5tCkhJ3XnIyIiqrd04hiCAIpCgT8QPuX7WSUMqEDRIUUsjtpLCABc+kpnd5GBhYiIaDKZhNavLCUFIMmnvt0X3FqIKaWdsSbUsYHFSd35iIiI6i2X1IJIWgpO+H1jQ0Q1O1avIc2I4wKLT28n7C8nLR4JERGRfRkNVnPyxIGl7NECi2BgMYdPX+lslGoRERHRqYqZMQBAznXq+hXgeBGL7JAiFscFlmBUa3YTFmkIVbV4NERERPZU1gNL0RWa8PuS31lFLI4LLKGo1h3XJalIp5yRComIiOpNzWrvkcZalZMpAaOIxRlLLBwXWHz+IApCq8ZOjTmjOx8REVHd5bQrJ6pn4ikhVyAGAPCWeIXFFJIsI6mveHZKdz4iIqJ6kwvaFRbhi074fW9Ym7Ew9huyO8cFFgBI6935snqNOREREY3nKmhTPZJv4ikhn76/UJCBxTxZPbAUUrzCQkRENBFjbYrsj034/YC+JtQpRSyODCx5l5YWiw7pzkdERFRvnpJ25cStN1w9mVHE4pFKyGXt3yrEkYGlqLcTdkp3PiIionoz1qa4gxOvYQmGoigJLQak4vZfYuHIwFJyWHc+IiKievOr2lUTo0P8yU4sYknH7V9168jAonq1tCjl2IeFiIhoIiG9I7w/PHFgAYC0pK0JzTmg6taRgUXSFxApBQYWIiKikxULeQSkPAAgFJ08sGQVPbA4oIjFkYHFWPHslO58RERE9ZSKHw8gwcjEi24BIK8HlpIDilgcGVhc+opnb9EZ3fmIiIjqKZPUAkta+OByeya9X0Fv28/AYhJ3ULu85Ss7o9kNERFRPRmd4NNS4LT3qxSxOGBNqCMDiy9idOfjlBAREdHJ8voVlow88U7NhnKliGXM7CHNmCMDS0APLMYKaCIiIjqumBkDAOSU0wcWSd9nSM7bf4mFIwNLKNoGAAhIeRQLeYtHQ0REZC+l9BgAIO86Q2Dxa2tC3Q6ounVoYGmtfJ4cs3+zGyIionoqZ7UAUtI7w0/GFYwBADwl+y+xcGRgUVwuJIUfAJB2QDthIiKiejIW0Z4psBwvYmFgMY3RnS+btH+zGyIionqS9DUpRmf4yXhD2pRQwAFVt84NLPpCojwDCxER0TiKsYjWd/rAEtCXWIQEA4tpjJXPBQe0EyYiIqonl94JXvafIbBEtMASlrIol0qmj2smHBtYCi69O59eukVEREQaYxGtEjh9YAnH2iqfp2y+JtSxgcXozqdm7N9OmIiIqJ68Za1PmScw+T5CAODx+pARXgDj9x+yI8cGlrJHT405+ze7ISIiqqeAqq1J8YRiZ7xvSgoCADIJe7cJcWxgEZXufGPWDoSIiMhmgvoiWn+45Yz3Ndr3272IxbGBRfLHAACuAq+wEBERGYSqIiQyAI5vZXM6WUXr1WL3IhbHBhYlEAMAeIoMLERERIZ0Kg5FEgDGd4afTN6lBRa7F7E4NrAY3fm8Dmh2Q0REVC/phHalpCgU+PzBM96/6JAiFscGFm/Y6M5n/3bCRERE9ZLRA0tKCkKSz/w2r+qBxWjnb1eODSy+sHaZKyjSFo+EiIjIPvJJ7UpJWjrz1RUAUH0xAICcGzNpRLXh2MAS1OflwiIFoaoWj4aIiMge8mntCktWnlpgkfSqW6Vg7xkLxwaWUEwLLC5JRSbNhbdEREQAUExrUzt5V2hK95f15nJ2L2JxbGDx+0MoCAUAkByzd7MbIiKieinr1T4FvfrnTNxBLbB4S7zCYgpJlivd+bIJe+9/QEREVC+qvni25J5aYPGEtKpbv8rAYpq0pF3uyibs3eyGiIiobvTAonojU7q7T++GG1DtXcTi6MCSVbTAUkjZu3aciIioXuS8thZFeE+/U7PBKGKJ2LyIxdGBJefS0mMxzSssREREAOAqalM7km9qV1iMwOKRSsjnMqaNa6YcHViK+vxc2ebthImIiOrFrQcWYwubMwmFYygLCYC9i1gcHVhKHu1yl8iOWTsQIiIim/CWtC1rXFMMLJIsI6kXsWTi9i1icXRgMRYUSXl7txMmIiKqF5+qBRaPXq48FalKEQsDiykkvZ2wwsBCREQE4Hi1jy889cCSlbXAkrfxmlBnBxZ/DMDxBUZERETNLqTvsefXy5WnIqc3mSumxswYUk04OrC4gtovw2vzdsJERET1UMjn4JcKAI5X/0zF8SIW+7YJcXRg8YRiAAB/OWXtQIiIiGwgdcKi2VBk6lNCRhGLauMiFkcHFl9YS48Bm7cTJiIiqodMQrtCkhJ+KC7XlB+nevQilpx914Q6OrAY83PGfB0REVEzyya1RbNpKVDV44S+JlS2cRGLowNLKNYGAAhKORQLeYtHQ0REZK2cvlVNRq/6mSpZDyxuG68JdXZgiR5fAZ2K27cUi4iIqB6KaS2w5JRgVY9T9J4tHhtX3To6sLjcHqSEHwCQitu3nTAREVE9lNJjAICCXqY8VZ6AFlh8ZQYW06T0dsLZBK+wEBFRc1P1RbNGmfJUefU1oX7VvlW3jg8sxjxdPmXfdsJERER1kdUCS9kztZ2aDf6IUcTCwGIaoztfIWXfZjdERER1kdcWzRp77U1VMKoVsYREFmq5XPNh1YLjA0teDyzlNAMLERE1N6WgBRapysAS1qtuZUkgadMiFscHlpJb+6WUs/atHSciIqoHY289Y6+9qfJ4fcgILwAgHbfnEgvHB5ayV2snLOXGrB0IERGRxTwlLbC4AtGqH2sUsWQS9qy6rTqwbN26Fddddx16enogSRKefvrp097/pZdegiRJp3y8//770x3zOMKnBxYbd+cjIiKqB19Z6/zuDk59HyGDUcSSS9pzSmjqGw3o0uk0li9fji9+8Yu44YYbpvy4Dz74AJHI8Tm19vb2ak89IeOyl4uBhYiImpxRlmxsDlyNrBICVKCQapDAsnbtWqxdu7bqE3V0dCAWi1X9uDNRjHbCJfs2uyEiIqqHoL63nl/fHLgaeVcEKB5vPmc3dVvDcsEFF6C7uxvXXHMNXnzxxdPeN5/PI5FIjPuYjDuk1Y77GFiIiKiJqeUyQiIDAAiGq58SKunN5tSsPatuTQ8s3d3dePjhh7Fp0yY89dRTWLx4Ma655hps3bp10sds3LgR0Wi08tHb2zvpfT0h7ZfiL9u32Q0REZHZ0qk4FEkAAEKx6q+wGEUswqZVt1VPCVVr8eLFWLx4ceXr1atX48CBA7j77rtx+eWXT/iYDRs2YP369ZWvE4nEpKHFH9F+KUEbd+cjIiIyWzpxFGEABaHA6wtU/XihBxbZpmtCLSlrvuSSS7B79+5Jv+/1ehGJRMZ9TCagtxMOizSEqtZ8rERERE6QTWj9U1JSEJI8jbd3fU2o0XzObiwJLNu3b0d3d3dNjmV053NLZWQzXMdCRETNKZcaAwBk9H4q1ZIDMQCAu2DPKyxVTwmlUil89NFHla/37t2Lt99+Gy0tLZg7dy42bNiAQ4cO4fHHHwcA3HvvvZg/fz6WLl2KQqGAJ554Aps2bcKmTZtq8gP4A2EUhQK3VEZybASBUPXNcoiIiJwur++pl1VC03q8O6jNWHhL9lxiUXVg2bZtG6666qrK18Zak5tvvhmPPfYYBgcH0d/fX/l+oVDAHXfcgUOHDsHv92Pp0qV49tlnce2119Zg+IAky0hKQbQggUziKICFNTkuERGRkxjlyLlpBhaP3mzOr9pztqLqwHLllVdCCDHp9x977LFxX3/jG9/AN77xjaoHVo20FEKLSFTm74iIiJpNOTsGACjqmwJXy2cUsaj2vMLi+L2EgOOXvwope9aOExERmU1ktcWyJc/0Aksgol1hCevN5+ymIQJLTk+TRZu2EyYiIjKdvgmw6pm8svZ0QjFtyxyvVEQua7/Q0hCBpejWfjnG5TAiIqJmIxW0tSfGpsDVCoaiKAsJAJA6Zr8dmxsisJT0NCkYWIiIqEm59P4pkm96V1hkRUFSL4lO23BNaEMEFlXvzifl7Fk7TkREZDZjE2BjU+DpSOuBxY5FLA0RWOCLAQAUm7YTJiIiMpvRP8UVmH4/soysrQnN27CIpSECi6ynSVfRnu2EiYiIzObXy5E9oZZpHyPv0qpu7VjE0hCBxRWMAQC8JXs2uyEiIjJbQNUqe3yh2LSPUbBxEUtDBBajnbCPgYWIiJpUSO+f4tcbwE1HSQ8saoZTQqbwhbXAErBpdz4iIiIz5XMZ+KQiACAYnX5gsXMRS0MEFiNNhmzanY+IiMhMqfjxNSehcGz6B9KLWGQbFrE0RGAJ6WkyJGVRKhYsHg0REVF9ZZLaFE5S+KG4qt4msEIKxAAArqL9llg0RmCJHb/8dWLKJCIiagbZhPbeZ/RRmS5FDyxeG1bdNkRgcbs9SAsfACA1Zr92wkRERGbKp8YAAFk5NKPjePQiFm/ZfmtCGyKwAEBS0n5J2SSvsBARUXMpprX3vpwysyssnrC2Y3NA5ZSQaYxUmUvar50wERGRmUqZMQBA3hWe0XECNi5iaZzAomiBpaBfFiMiImoWalar6im5ZzYldDywZKCWyzMeVy01TGDJ681uSmlOCRERUXMRet+Usmd6OzUbwnoRiywJJBP2ah7XMIGl5NYugwkbthMmIiIyk5zXqnpU78wCi9cXQFZ4AABpmxWxNExgKevd+URuzNqBEBER1Zlc0AKL5Jv+Ts0Go4glk7DXmtCGCSxCDyxGyiQiImoWrqJWhiz7YzM+VkbWKo3sVsTSMIFF0n9JroL92gkTERGZyatv/qsEZn6FJatoSyyKaa5hMYXRnc9tw3bCREREZjIavXmCs2Z8rIJeGl1MMbCYwq135/OVOCVERETNJVAJLLEZH6uoV92qNitiaZjA4glpqdJvw3bCREREZgoiAwDwR1pmfCyjNNpuVbcNE1iMX1JQMLAQEVHzUMtlhIQWWAI1CCzCFwMAyHl7rQltoMByvJ2wUFWLR0NERFQfqVQcsiQAAKFo68wP6NcW7ioFey2xaJjAEo61AQA8Uhm5DK+yEBFRc8jEtfLjgnDB55/Z5ocAIAe0JRbuIgOLKQLBCEpC+3GScXvVjhMREZklk9C2pElKMw8rAODSA4vXZlW3DRNYJFk+3p0vbq92wkRERGbJ6eXHmRoFFm9IWwfjV+01W9EwgQUA0vovK5vkBohERNQcCnpgySoz26nZ4AtrV1iCKq+wmCajd+crMLAQEVGTKGXGAAD5GgWWQPR4EYudNFRgMX5ZBZu1EyYiIjJLWQ8sRXdtAkswqhWx+KQicln7hJaGCiwFvTuf8csjIiJqdGpO65dS0t8DZyoUjqEsJABAykZFLA0VWEo27c5HRERkFimnlR+rnnBNjicrClJSAACQHrNPEUtDBRbVGwMASDbrzkdERGQWWW/wJnwz36nZkNKrbu1UxNJQgQU+7QqLwsBCRERNwqUHFqmGgSUra4Elz8BiDqM7n8tm7YSJiIjM4ipp/VKUQKxmx8y5tOmloo2KWBoqsLj0X5anxMBCRETNwVfS+qW4axhYCi6jiIWBxRQevTufr2yv7nxERERm8evveZ7QrJod0yhiUW1UxNJQgcUb1gJLkIGFiIiahF9kAADeGgYW1auth5EYWMzhj2jd+YKCgYWIiJpDWO9IG4jULrDAFwNwvALJDhoqsAQj2hWWsJRFuVSyeDRERETmymXT8EpFAEBA71BbC5Jfu8JipyKWhgos4djxX5aduvMRERGZIRXXyo5VISEcjtXsuEpQu1rjKdlnA8SGCixujxcZ4QUApOL26c5HRERkBqOxW0ryQ1aUmh3XrQcWHwOLeZJ6d75Mwj7NboiIiMyQTWplxxkEa3pco4gloNpnTWjDBZaM3p0vx8BCREQNLp/SA4tc28DiD2tFLCEbFbE0XGDJKkZ3PgYWIiJqbIXUGAAgp4Rqetxg1AgsGajlck2PPV0NF1gKbi2wlGzUTpiIiMgMZb1PivHeVyvhmBZYZEkgmbDH+2nDBZai237d+YiIiMxgvNeVXLW9wuL1BZAVHgBA2iZVtw0XWMp6O2HBwEJERA1OZOMAjr/31VJK0tbFZGxSddtwgUUY3fnycWsHQkREZDI5rzV2U721Dyxpo4glZY81oQ0XWOzYnY+IiMgMiv5eZ7z31VLOKGJhYDGHEtCa3biL9ml2Q0REZAZXSSs7ln21Dyx5l1HEMlbzY09HwwUWVyAGAPCWeIWFiIgam0f/49z4Y72WinrlUTnDKiFTePVmN34bdecjIiIyg6+svdd5QrGaH7vs0a7aiJw91oQ2XGDx6e2EgwwsRETU4AJqGgDgDdX+CouqTzPJubGaH3s6Gi6wBPTufGGRhlBVi0dDRERkngC0wOIP1z6wSP4YgOMLe63WcIElpAcWj1RCLpu2eDRERETmUMtlhEQWAODXZxdqSdEDi7vIwGKKYCiKktB+rJRNuvMRERHVWio5BlkSAIBQ1ITAEtSO6S3ZY4lF1YFl69atuO6669DT0wNJkvD000+f8TFbtmzBypUr4fP5sGDBAjz00EPTGeuUSLKMpN6dL22T7nxERES1ZrzH5YUbPn9td2sGji/k9Zft0Sak6sCSTqexfPly/PCHP5zS/ffu3Ytrr70Wl112GbZv345vfetb+OpXv4pNmzZVPdgpj1HSu/Ml7NHshoiIqNYy+qaExh/ptebXq24DNilicVX7gLVr12Lt2rVTvv9DDz2EuXPn4t577wUALFmyBNu2bcPdd9+NG264odrTT0lWCQEl+7QTJiIiqrV8SgssWSlgyvH9EaOIxR6BxfQ1LK+88grWrFkz7rZPfepT2LZtG4rF4oSPyefzSCQS4z6qYbQTLqXt0eyGiIio1gpGYFFqu1OzIRRrAwD4pKItilhMDyxDQ0Po7Owcd1tnZydKpRJGRiZeY7Jx40ZEo9HKR29vb1XnNLrzMbAQEVGjKmbGAAB5swJLOAZVSADsUcRSlyohSZLGfS2EmPB2w4YNGxCPxysfBw4cqOp8RZt15yMiIqq1sh5YjD/Sa01WFCT16aZ03PolFlWvYalWV1cXhoaGxt02PDwMl8uF1tbWCR/j9Xrh9XqnfU5jm23JJt35iIiIas34o7zkiZh2jrQUQlSkkU1YX3Vr+hWW1atXY/PmzeNue/7557Fq1Sq43W5zTmq0E87bo9kNERFRrUl6YFFNDCwZWZtuMhb4WqnqwJJKpfD222/j7bffBqCVLb/99tvo7+8HoE3n3HTTTZX7r1u3Dvv378f69euxa9cuPProo3jkkUdwxx131OYnmIDs11oUu4ucEiIiosYkF7T+KMJnXmDJubTppqINqm6rnhLatm0brrrqqsrX69evBwDcfPPNeOyxxzA4OFgJLwDQ19eH5557Dl/72tdw//33o6enB/fdd59pJc0AoARjAI5vu01ERNRoXHrLfFmfVTBDwRUB8sfXy1ip6sBy5ZVXVhbNTuSxxx475bYrrrgCb731VrWnmjaP3k7YZ5PufERERLXmLmr9UZRAzLRzlDxhIA2I7Jhp55iqhttLCDi+zXagbI9mN0RERLXm1d/jXCYGlrJXv3pjg6rbhgwsAX3H5iCsb3RDRERkBp8eWDz6H+nmnEQLLEp+zLxzTFFjBha9nXAEGZRLJYtHQ0REVHtBof1R7jMxsBhFLErB+iUWDRlYwno7YQBIJawvxSIiIqq1sB5YApEW085hrI/xlKxvE9KQgcXj9SEjtMZzqTHrm90QERHVUi6bhkfSZhDMDCxuo4ilxCsspknp221nk9bvf0BERFRLxt4+qpAQCsdMO483rBexqNYXsTRsYDG68+USDCxERNRYMgmtkVtKCkBWFNPO49fXhAYFA4tpsorWna/AHZuJiKjBZPXAkkbA1PME9arbkMhALZdNPdeZNGxgyevthEsMLERE1GAK6TEAQEYJmXqekB5YFEkglRwz9Vxn0rCBpahvBqXaoDsfERFRLRmzB3mTA4vPH0ROaBsVp+PWFrE0bGAxdq+0QzthIiKiWipntM6zxmyCmZKSFooycWvXhDZuYPHFAABy3vracSIioloyZg9KbnOvsAAnFLEkrV1i0bCBRaq0E7Z+/wMiIqJaEjntj3FjNsFMWX3aKZ86avq5TqdhA4sc0GrH3UVeYSEiosYi63+Mq17zA0uliCXDKyymcAe1wGKH7nxERES1pBS0P8aN2QQzFd16EQsDizk8Ib2dsA268xEREdWSq6i9t8l+8wNL2Zh2ylq7xKJhA4svrAWWIAMLERE1GK8+e+DSNyc0k+rVQpGUGzP9XKfTsIHF6M5n7GZJRETUKHxl7b3NWP5gJskfA3B8GsoqDR9YvFIR2TTXsRARUePw67MH3pD5gUUxilgKY6af63QaNrCEwjGkhB8AMNz/ocWjISIiqp2gPnvgj7SYfi5/Rx8AoCV/yPRznU7DBhZJljHg7gUAjO7fYfFoiIiIaqNcKiEsZQEAgbD5V1i6FiwHAPSog8jnMqafbzING1gAIB7UUmF+cJfFIyEiIqqNEzchDEXNv8LS2tWLpPBDkQQG975n+vkm09CBpdRyFgDAfewji0dCRERUG2l9T5+ccMPrC5h+Pm3GYi4A4Og+62YsGjqw+LqXAABi6b0Wj4SIiKg2svqePikpWLdzJowZi6H363bOkzV0YGmdvwwA0FM6CLVctng0REREM5dLanv6ZOT6BZZS69kArJ2xaOjA0t13LgpCQUDK4/BBTgsREZHzFfRNCLOy+Ts1G3zd5wCwdsaioQOL2+3BgNIDADiyl5VCRETkfKWM1iI/76pfYGmzwYxFQwcWADjqnw8AyAywUoiIiJyvnB0DABT1XZTroXv+ksqMxfChj+t23hM1fGDJxxYBAKQRNo8jIiLnEzmtRX7ZU7/A4jphxmJ4jzUzFg0fWNyd2rxbOLnH4pEQERHNnJTTpoRUYxflOjnq1yqFMgPW9GJp+MASnbsUANBd7Ld4JERERDMnG5sQ+mJ1PW9lxmJ0d13Pa2j4wDJ70fkAgFlIYOzIoMWjISIimhlXQdvQV/LX9wqLu3MxAOtmLBo+sARCUQyhHQAwuOddi0dDREQ0M/7CCADAFWqr63lj87RKoa7C/rqe19DwgQUAhn3zAADJAzstHgkREdHMtBcHAACR2Yvret6ehecBAFqQQHx0qK7nBpoksGQiCwAA6pEPLB4JERHR9KWTY2jDGACgY96Sup573IzFR/WfsWiKwCK1aynUH7emdpyIiKgWhvZpe/kcQxjRWfWdEgKAYZ+2CWLiYP0rhZoisITnaJVC7Xlr5t2IiIhqIX5ICyzDrh5Lzp+JLAQAqMP13wSxKQJL50KtUqhLPYJsOmnxaIiIiKanMKzNFCQDvZacX2rTZywS9Z+xaIrA0tLWjTGEIEsCAx+zUoiIiJxJHtM2HyxG51ty/tCccwEAbbn6z1g0RWCRZBmDbm3e7Vg/K4WIiMiZgimtCaqrbaEl5+/SZyy61WHkMqm6nrspAgsAJENapVBpqP7zbkRERLXQWtBKmsM9Z1ty/pb2HsQR1GYs9tT3AkDTBBa1VfvlesasaSlMREQ0E7lsGh1CaxrXXueSZoM2Y6H1Nju2v76bIDZNYPH3aL/cWRlWChERkfMc3v8BZEkgJfxoaeu2bByJkLYJYuFwfXubNU1gae/T5t1mlw+hVCxYPBoiIqLqjB3UAsKQqweSbN3bt9p6FgDAc6y+MxZNE1i65p6FrPDAI5UwuJ8db4mIyFmywx8BAOL+OZaOozJjkd5b1/M2TWCRFQUDLu2XPLq3vvNuREREMyUd1XZJLkTmWTqOtvnankI95UMol0p1O2/TBBYAGAvMBwDkBndZOxAiIqIq+fWSZqXVmpJmQ9fcxcgLN3xSEUP99ZsWaqrAUpilzbspR1kpREREztKSPwQACHZbU9JsUFwuDCizAQAj++o3Y9FUgcXTpc27RVJ7LB4JERHR1BULeXSqwwCAtrmLLR4NcFSfscgO1G8TxKYKLLPma5sg9hT7IVTV4tEQERFNzfDBj+CWysgJN9q751s9nOMzFqMf1u2cTRVYZi9YhrKQEJayGB06YPVwiIiIpuRov9alfUjphqwoFo8GcHdpV3lCqfpVCjVVYPH6AhiUuwAAQ3vesXg0REREU5M5rJU0H/NZW9JsmDVXqxTqruOMRVMFFgAY8c8HAKQP1W/ejYiIaCaEXtKcD1tb0myYvXAZVCEhhhSOjQzW5ZxNF1iyUb0c7AibxxERkTP4ktq2MlLrAotHovEFQhiSOwAAQx+/W5dzNl1gUTq0ebdgkpVCRETkDLHcQQCAv/Msi0dy3BGfdrUnebA+MxZNF1givVqlUEe+3+KREBERnZlaLqO7PAQAaO09x+LRHJeNaDMWok4zFk0XWLoXLgcAdOAoEmOjFo+GiIjo9IYH9sIrFVEUCjp7re1yeyJZn7EIJD6uz/nqchYbic5qwwhiAOo370ZERDRdo/36Ls1yJ1xuj8WjOS4851wAQHtuf13O13SBBQAOe7R5t/iBnRaPhIiI6PTSQ1pztqPe2RaPZDxjxqIbR5BJxU0/X1MGlpQ+71Y6zEohIiKyt/KIViSSC821eCTjxdq6cAwRAMDAx781/XzTCiwPPPAA+vr64PP5sHLlSvz617+e9L4vvfQSJEk65eP999+f9qBnrE3bOMoX/8i6MRAREU2BN7kPACBa7FHSfKJBtxaixvptGFiefPJJ3H777fj2t7+N7du347LLLsPatWvR33/6qpsPPvgAg4ODlY+zzrKuNCvQo827tWX3WTYGIiKiqYhmtZJmX+cii0dyqmRYC1HFw+ZfhKg6sNxzzz34kz/5E/zpn/4plixZgnvvvRe9vb148MEHT/u4jo4OdHV1VT4UC/dC6Fx4PgCgWx1CIZe1bBxERESnI1QVXaUBAMCsOfYpaTaIVu3ig3fM/BmLqgJLoVDAm2++iTVr1oy7fc2aNXj55ZdP+9gLLrgA3d3duOaaa/Diiy+e9r75fB6JRGLcRy21d89DSvjhklQM7OXCWyIisqejRw4hKOWgCgld8862ejinMGYsWjL7TD9XVYFlZGQE5XIZnZ2d427v7OzE0NDQhI/p7u7Gww8/jE2bNuGpp57C4sWLcc0112Dr1q2Tnmfjxo2IRqOVj97e3mqGeUaSLGPArR3z6L4dNT02ERFRrRzZr021HJba4PUFLB7Nqdr6tE0Qe8oDKBULpp7LNZ0HSZI07mshxCm3GRYvXozFixdXvl69ejUOHDiAu+++G5dffvmEj9mwYQPWr19f+TqRSNQ8tMSDfUD8QxSGLFz8S0REdBrJAa2kedQ7G90Wj2UiXb2LkBFeBKQ8Duz/AL2LzjPtXFVdYWlra4OiKKdcTRkeHj7lqsvpXHLJJdi9e/ek3/d6vYhEIuM+aq3Uos27uY6xUoiIiOypPKJ1kU0H7VXSbJAVBQOuOQCA0b3mzlhUFVg8Hg9WrlyJzZs3j7t98+bNuPTSS6d8nO3bt6O729qs6OteAgCIpfdaOg4iIqLJuBP7AADqrD5rB3IaY4H5AIDc4C5Tz1P1lND69evxh3/4h1i1ahVWr16Nhx9+GP39/Vi3bh0AbTrn0KFDePzxxwEA9957L+bPn4+lS5eiUCjgiSeewKZNm7Bp06ba/iRVau07D3gZmF06ALVchmxh1RIREdFEwpkDAABvh/1Kmg3FlrOA5P+GcnTymZNaqDqw3HjjjRgdHcV3v/tdDA4OYtmyZXjuuecwb57W7n5wcHBcT5ZCoYA77rgDhw4dgt/vx9KlS/Hss8/i2muvrd1PMQ0985egIBT4pQIGD36E7nmLz/wgIiKiOurQS5pjs+37HuXpPAfYD0RS5s5YSEIIYeoZaiCRSCAajSIej9d0Pcu+7y7DfPUA3r3iEZx/1f9ds+MSERHNVPzoEUTv066sZO7oRyAUtXhEE9v73hvo+4d/hwQCCP+XQ5Dk46tNavn+3ZR7CRmO+rU5wczAexaPhIiIaLzh/dqakCOYZduwAgA9C5aiLCREkMHo8EHTztPUgSUf05KrPPqhxSMhIiIaLz6gbdA74u6xeCSn5/UFMCh3AQCGPn7HtPM0dWBxd2pzguHkHotHQkRENF5pWGu7kbJpSfOJRvzzAQDpg+bNWDR1YInNWwYA6CoesHgkRERE4ynx/QCAUnS+tQOZglxE30l6xLwZi6YOLD0LtY58s5DAsSODFo+GiIjouFBaq7h127ik2SB3ahszBhIfm3cO047sAIFQFINoBwAMmjjvRkREVK324iEAQLTHfpseniwyR9sEsSPff4Z7Tl9TBxYAOOLT+sekDnLXZiIisodMagxtGAMAdMxbYu1gpqB70QoAQCdGkUocM+UcTR9YMpGFAAD1CCuFiIjIHob2aRVCxxBGtKXd4tGcWXRWG0YQAwAMfPSuKedo+sAitWuVQoG4efNuRERE1Ygfeh8AcMRlxz2aJ3bYo1UzxQ+YM2PR9IElrM+7tef2WTsQIiIiXf6wVtKcCNi/pNmQCmuVQqXD75ty/KYPLF0LzwcAdOMIsqmExaMhIiICpLF9AIBiZL6l46iGaNMWB/tMmrFo+sDS0jEbxxAGABz6eIfFoyEiIgKCekmzq22hxSOZuuBsbcaiJbvPlOM3fWABgCG3dsltrP+3Fo+EiIgIaMtrJc3hnrMsHsnUtfdpvc16yoMoFvI1Pz4DC4BkSNsEsWjSvBsREdFU5XMZdIgRAECbA0qaDZ2zFyAtfHBLZQzs3VXz4zOwAFD1eTfv2EcWj4SIiJrd0P4PIUsCKeFHa7u9Nz48kSTLGHD3AgBG99V+iQUDCwB/jz7vltln7UCIiKjpjR3UrvYfdnVDkp31Nh0PzAcA5Id4hcUU7X1apVBP+RBKxYLFoyEiomaWPbwbADDmm2PxSKpXatFmLFxHd9f82AwsALp6FyErPPBIZQzuq30qJCIimirp6F4AQMFBJc0Gb7fWjDWa3lfzYzOwAJAVBQMuLcmO7uOeQkREZB1/SitpVtoWWDyS6rXM0yuFigcgVLWmx2Zg0R0LaJVC2UFeYSEiIuu05A8CAAJdzilpNvQsWIqiUBCSshge2FvTYzOw6IqztBeGa5SbIBIRkTVKxQI61WEAQNtc55Q0G9weLwaVLgDA8J7a9jZjYNF5us4BAETTeyweCRERNavDBz6GWyojJ9zo6Jlv9XCmZdQ3HwCQPlTbJRYMLLqW+csAAN0mzLsRERFNxdGD2rKEIaULsqJYPJrpycW07QSk0dpWCjGw6HoWLEVZSAhLWYwM9Vs9HCIiakKZIe1N/pgDS5oNrg5txiKUqO0miAwsOq8vgAG5GwBw+ON3LB4NERE1IzGqLVTNhedbO5AZiM5dCgDoLOyv6XEZWE4w6p8HAEgfes/ikRARUTPyJbU3ebmlz+KRTF/3Qq0ZaxvGED82WrPjMrCcIBvVt/EeYaUQERHVXyx3AADgd2BJsyEcbcEwWgAAh/fWbk8hBpYTKPq8WzDJSiEiIqovtVxGV3kIANDae47Fo5mZYe9cAEDi0Ps1OyYDywkivfq8W762825ERERnMjK4Hz6piKJQ0Nm7yOrhzEg6rM1YqEdqVynEwHKCnkXLAQDtOIbEWO3m3YiIiM5kpF8raT4sd8Dl9lg8mhlq1zZB9Cb21eyQDCwniMRacQSzAACDH7FSiIiI6ic9qK2fHPU6t6TZEJx9LgCgNVu7GQsGlpMc9mqVQvED3ASRiIjqpzSqrZ/MhuZaPJKZ61qgzVh0i8M1OyYDy0lSYW13zPJw7RYKERERnYk3oV+NcHBJs6G1qxdJ4YciiZodk4HlJFKbNu/mj9e2Qx8REdHpRLNaSbOvw7klzQZJljHgru2VIgaWk5gx70ZERHQ6QlXRWRoEAMzqXWzxaGojEaztlSIGlpN0LtTm3XrUQeRzGYtHQ0REzeDYkQGEpCxUIaFzbmMEllLr2TU9HgPLSdq65lbm3Xb860+sHg4RETWB4f3auslhqRU+f8Di0dSGr7u2ze8YWE4iyTJ2tVwDAFi1/Vt49e+/inKpZPGoiIiokSX1kuYRj/NLmg1zz78CGVG7fjIMLBO48Ms/wSudXwAAXHLop/jt3Z9G/NgRi0dFRESNqjyiFXqkG6Ck2dDaOQcD1/60ZsdjYJmAy+3B6i89hG0r/xY54cby3BtI3ncZ9u3aZvXQiIjIRoSq4rf/9gzef+NXMzqOK74PAKDG5s98UDayaMVlNTsWA8tprLruP+LgZ/8nBtGOOWIQHb+4Ftv/tXZpkYiInOvo8CFsv+czWParP8Si//U5fPzuy9M+ViTTDwDwdjh7DyEzMbCcwaLlvwPvl7dip2c5AlIeF7zyVbz649uhlstWD42IiCzy1r/+D+CBS3BhaisAwCWpUJ/5Ksql4rSO166XNEdnO3uXZjMxsExBS0cPFv/5r/Bqx+cBAJcc/Al23P1pxI+NWDwyIiKqp/jRI9h2zw248JXb0IIE9slz8c7lf48EgjirtBtvPPlfqz/msSOYhSQAoGt+Y5Q0m4GBZYpcbg8u+fKP8MYFd2nrWrKvI3HfZdi/602rh0ZERHXwzgv/gMJ9F2FV4lcoCwmv9NyE7m+8huVX/wF2nffnAIDzP7wfA3ur29rlyD5tl+YRxBAMx2o97IbBwFKliz7zJRz4/acxhDb0igG0/eJabH/+CauHRUREJknGj+L1H/w/WL71FrTjGA5IPdh93Sas/rP/D16f1jPlouv/E3Z6zkNAymPkyVshVHXKx48PaCXNR9yzTRl/o2BgmYazVnwS7i9vxU7P+QhKOVzw8q149cfrua6FiKjB/PbX/xPp//4JfOLYs1CFhFc7bkTbHa/jnFXXjLufrMiIfO5+5IUb5+e24c3/9aMpn6N4RCtpTgYap6TZDAws09TaMRtn3/ErvNr+BwCASw4+gnfvvhaJsVGLR0ZERDOVScXx2g+/iGX/+yZ04QgGpE68/+mf45IvPwx/MDzhY3rPWo63+m4BACx8629w7MjAlM6ljO0FAJRj82oz+AbFwDIDbo8Xl9z693hjxX9FXrixIvsq4j/4JN7e/DNebSEicqhdr/0rjv7dJ3DxyFMAgNdar0d0/es4d/XaMz525Re+gz3yfMxCEh8/8dUpnS+Y1nZpdrezpPl0GFhq4KLrb0X/9U/hMFrRKwaw4jdfwoG/OR9vbLqXGygSETlELpPCqw+uw+LnbsQcMYQhtGHH1Y/j4q/8dMqLYT1eL4r/171QhYRV8c3YsWXTGR/TXjwEAIj0sELodBhYauSsCy6H59bf4JXum5AUfsxTD+KiHXcicde5ePWn/5kl0ERENhU/egSv/fxvMPL9Vbjk8M8hSwKvx65F4PbXcd7ln6n6eItXXoXXOz4HAGh76ZvIphKT3jeTiqMdxwAAnfPZg+V0JCGEsHoQZ5JIJBCNRhGPxxGJRKwezhkl40ex85/uw4KPfooOHAUApIUPO7p+H33X3YHOObzsR0RkJaGq2PX688i8/GMsi78En6Q1fBtBDIcu/29YfvXnZ3T8VHIMqb9biS6M4LWuL+DidQ9NeL+9O19D3z+uwRhCiH3n0IzOaUe1fP9mYDFRIZ/D2//8CNrf/RH61P0AgKJQ8HbsGrStuQN9Sy+2eIRERM1l7Mgg3n/+YXR//A+Ypx6s3L5Hno8ji7+AJZ+6BZFYa03O9fYL/4AVW29BWUjY+/vPYNGKy0+5z/bn/wcuePk2fOg6G2f/5zdqcl47qeX7t6tGY6IJeLw+fOL6WyH+/Zfwzkub4Hr1PiwtvIuL4s8D//g83vmni+C67Hacu/paSDJn54iIzCBUFe+98hyyrz6K8xNbcIlUAgCkhRc7W/4dYpf9Gc5acTkW1Pj/wyuu/gO8sf3nuCj5AqR/+k8onfsqXB7vuPvkD+8GACT8LGk+EwaWOpBkGcuv/hxw9efw4VtbkHrh77A8uRXLc28Am/8Ddr+4CIkLv4SlV34evkDI6uESETWE0cMHsfv5hzF7zz9iqdBLjCVgt7IIR8/5As5d88f4RLTF1DHM/3/vw9iDl2BheQ9ee/J7uPgP/3rc9+VjWklzMcqS5jNhYKmzsy+8ArjwChz8eCcO/vP3sfzIszir9BHw+tdReO0v8J53CRKdFyN8zlVYeMGVDDBERFXIZdPY/dq/oLjtp1iW/DdcImktJlLCj52ta9By+S04a8VldRtPe2cvXlu+ARe/820s/+hBDHz8efQsXFr5fiCt7dKstC2s25icimtYLHZ0+BA++Kf/joUHNlUW6BoKwoWPTgwwF14Fnz9o0UiJiOwnPjqEfdtfRPbj3yA28iYWFD6ER5/yAYAPXWdjbMl/wNI1f2TZPj1CVfHbu67GeYXt2OldgXP/4sXKMoDBv1qEbnEEu679/7HkE79ryfjMxEW3DUioKg7u2YnBtzdD7v8N5iXerJS6GfLCjY+85yDReQki51yFhRdeyQBDRE1DqCoG93+IgR0vQN33CjrGtmO+euCU+40gho/arkLbFf8Ri85bbcFIT3Xg451oe/xK+KUC3lzxPay8/ivI5zJwb+yBLAkcWbcD7V2Nt46FgaUJGAFmYPvzUPp/g3nJtyYNMMnYuZBaFyLQfQ7a55+LjtkLICuKRSMnIqqNcqmEfe+9jpH3tsB16DX0Jt855Uo0AOyX5+BwdAWkeavRc/5V6Jm/xJaFDC8//l9w6Z4fII4gyl96Dan4KOb+7AqkhQ+BOwdtOeaZYmBpQkJVcfDjHRh4+1dQ+n+D+ck30YaxCe+bE24MKj045p+LfLQPrrZFCM0+Bx3zl6Klvach/6MgIucRqorR4UM4OrAHqcN7UBjdDyl+AJ70ACL5IXSVBhCUcuMeUxQK9rgX4VjrhfAu/CTmrbgKLR3O2OW4WCxg/10XY1F5D96KXA3Xis/j/K1/ho+VBVj4l9utHp4pLA8sDzzwAL7//e9jcHAQS5cuxb333ovLLpt8EdOWLVuwfv167Ny5Ez09PfjGN76BdevWTfl8DCynEqqKAx+9i6EdL0Id/gC+5D605PrRVR6CR5p8H6Ok8GPIPQfxwDwUgz2QQm1QQu3wRjsRmNWJcGs3Ym3dlS3TiYimo1QsIH50GKljw0gdHUJmeC9Kx/qhJA7CnxlArDCEdvVIpWHbZFLCjz3+pUh3rkLk7MuxYMXlk24+6AQfbP81Fj797+GSVOzwrcJ5uW14K3g5Lvzzf7J6aKawtA/Lk08+idtvvx0PPPAAfud3fgc/+tGPsHbtWrz33nuYO/fU+be9e/fi2muvxS233IInnngCv/nNb/DlL38Z7e3tuOGGG2Y0+GYmyTLmnr0Cc89eMe72UrGAQwd2Y2T/e8gOfQhp9CMEkvvQWjiILvUIwlIW4dJuILEbmLxbNJLCj7gcRUqZhaxnForeFpT9rZBC7VACLVD8Ebj9EXgCYfiCUXhDUQTCsxAMRaG4WHxG1AjKpRIy6QRy6QRy6TiyyTHk4kdQSI2ilBqFyByFnD0KpRCHpzAGfymOYDmBkEgiggxaAZy2BZsEqELCiDQLR12dSPu7UQjNgRybA1/7PES7F6H3rBU4v4H+n7L4gsvw8qufx6WHf4bzctsAAPkIS5qnouorLBdffDEuvPBCPPjgg5XblixZguuvvx4bN2485f5/8Rd/gWeeeQa7du2q3LZu3Tq88847eOWVV6Z0Tl5hqY1cNo2hfbtw7MD7yA99AKSG4ModhTc/ikBpDJHyGGIiAfdprtBMRUZ4kZH8yEl+5OQACkoARSWAkisAVfZCVTwQihdQPIDLB+HyQHL5ILl9kFxeyG4fJJcPiscHxa19rbjckF1uyIoLsuKGorggu1yQFQ9kRYbi8kBWXHDp93O5XFBcbkiSDEmSAACSJGlfa1/oX0uA/n0iOxJCQFVVlEtFlEtFlEpFqOUSSqUiSsU8irksivksSoUsSoUcyoUsyoUs1GIO5WIeaiEHUcpBFLV/USpAKmYgFdOQSxm49A93OQuvmoVXZOETWfhFDn6pMOPxJxBEXIoi7ulAxt+DcmQOlFlzEeiYj1ndi9A+uw8er68Gz5RzpFMJjP3dKswWhwEAr533HVx8w9csHpU5LLvCUigU8Oabb+Kb3/zmuNvXrFmDl19+ecLHvPLKK1izZs242z71qU/hkUceQbFYhNvtPuUx+Xwe+Xy+8nUicZpLATRlPn8Q85eswvwlqya9j1BVxMdGkRgdQOroELJjwygmDkNNHYGcGYErNwp3MQFPKQOvmoFPzcCPLAIiW5mKCkh5BJAHxBhQhvbhAKrQgosAUIaMMhSokFGWFJQhQ8UJ/0rav6qk6P/KUCUFAhIEJED/V0jHP4ckQfvrQNZvR+W+ACCLMiSokEQZslAhowxJqJCFdlYJKhT988r39SMa5z1+blTOrZ4wHui3i8rPoaAsuaBKLu1nkVxQZe1rIbsg9K+F5IJQ3ICkL+YWKiC0sUIIQB8nhABwwudChYQyJOPvIun4GCf8XDr1e2JcoJzgc2n87QLjA6jxHGnfPPHvs4k/18YqACH034eqP05onxvPtPEcaK8e/XelVh4jVX5HKiQhIKNc+b72+xv/+YmvMpcwPlfhlspQANR9Gf0JT2NJyMhIPmQQQFqJIOuKoOCOouidBdUXgxRogRJshTvcBl+4FYFYO8KzOhGZ1Y6Iy4UIgN56j9/GgqEIdl9xF2a/9EXt607uLzcVVQWWkZERlMtldHZ2jru9s7MTQ0NDEz5maGhowvuXSiWMjIygu7v7lMds3LgRf/VXf1XN0KhGJFlGtKUd0ZZ2AMuremw+l0EmGUc2FUcuHUchHUchk0Axm4SaS6CcTwGlPEQpD5TykEo5SOUCpHIeklqAUs5DVguQ1QJcagGKWoBLFOAWhcobtCLKUPQooZz09UyvDMnS8Tct5UxJS5z0rx1NNjY7j5k0U7zoVxAKinCjIHlQgBtFyYOS5EZR9qAkeVCWjQ8vVMULVfZAKB4Ilx/CE4LkDULyhKD4QlB8Ybh8YbgDIXgDEXgDEfiDEfhDEXi9fkRkGby+XTsrrvwsfr33DfhGfovzL2q8/itmmNbEoHTSJXQhxCm3nen+E91u2LBhA9avX1/5OpFIoLeX+dzuvL4AvL4AZrWfGkLrQagqVFVFqVRAuVREsVg8/he1UCuvOyEEhFBP+Fz/ixoAVAEBAVUto1wqQS0VoaolqOUS1HIZarlY+VeoZailEoT+faEWIVTtL3Pj+BBl7fjGbUI7PoQAVO2v9cp/D4oLkqxAVlyApIz7WpIVSLILiuICFAWy7Krcrv1s6gnn0X4eoZb1f/WfH9o5jZ9flEsolwsQpSJEqQBVLUGUi9rXZf1DLQHlIlAuAWoRUEsAJEBWAEnWrrhIsnZlxLhNViBJMiBr35ckWbtde0LGPRfG18d/P0K/enPC55Vf8MRXSsblrxPvD4HjV2G084/7P47x/5/K/4dO+K4xZmNKcaKvZRmAPs2ofw5FgSS5IMkSJNn4vUnabYqi31f7XJbkcb9LxeXSPtenOhWXAkVxQ6pMdboqU52KrN3fIyvwAGA3Jme67IunLqOgyVUVWNra2qAoyilXU4aHh0+5imLo6uqa8P4ulwutrRMvx/J6vfB6vRN+j2gykixDkWUu+iUiakBVNeTweDxYuXIlNm/ePO72zZs349JLL53wMatXrz7l/s8//zxWrVo14foVIiIiopNV3UFs/fr1+PGPf4xHH30Uu3btwte+9jX09/dX+qps2LABN910U+X+69atw/79+7F+/Xrs2rULjz76KB555BHccccdtfspiIiIqKFVfe38xhtvxOjoKL773e9icHAQy5Ytw3PPPYd587Q68sHBQfT391fu39fXh+eeew5f+9rXcP/996Onpwf33Xcfe7AQERHRlLE1PxEREZmilu/f3FSGiIiIbI+BhYiIiGyPgYWIiIhsj4GFiIiIbI+BhYiIiGyPgYWIiIhsj4GFiIiIbI+BhYiIiGyPgYWIiIhszxHb2hrNeBOJhMUjISIioqky3rdr0VTfEYFldHQUANDb22vxSIiIiKhao6OjiEajMzqGIwJLS0sLAKC/v3/GP3CjSSQS6O3txYEDB7jP0kn43EyMz8vk+NxMjs/N5PjcTC4ej2Pu3LmV9/GZcERgkWVtqU00GuWLYRKRSITPzST43EyMz8vk+NxMjs/N5PjcTM54H5/RMWowDiIiIiJTMbAQERGR7TkisHi9Xtx5553wer1WD8V2+NxMjs/NxPi8TI7PzeT43EyOz83kavncSKIWtUZEREREJnLEFRYiIiJqbgwsREREZHsMLERERGR7DCxERERke7YPLA888AD6+vrg8/mwcuVK/PrXv7Z6SJb7zne+A0mSxn10dXVZPSxLbN26Fddddx16enogSRKefvrpcd8XQuA73/kOenp64Pf7ceWVV2Lnzp3WDLbOzvTc/NEf/dEpr6NLLrnEmsHW0caNG3HRRRchHA6jo6MD119/PT744INx92nW181Unptmfd08+OCDOP/88yvN4VavXo1//ud/rny/WV8zwJmfm1q9ZmwdWJ588kncfvvt+Pa3v43t27fjsssuw9q1a9Hf32/10Cy3dOlSDA4OVj527Nhh9ZAskU6nsXz5cvzwhz+c8Pt/+7d/i3vuuQc//OEP8cYbb6Crqwu/+7u/i2QyWeeR1t+ZnhsA+PSnPz3udfTcc8/VcYTW2LJlC2699Va8+uqr2Lx5M0qlEtasWYN0Ol25T7O+bqby3ADN+bqZM2cO7rrrLmzbtg3btm3D1Vdfjc985jOVUNKsrxngzM8NUKPXjLCxT3ziE2LdunXjbjvnnHPEN7/5TYtGZA933nmnWL58udXDsB0A4pe//GXla1VVRVdXl7jrrrsqt+VyORGNRsVDDz1kwQitc/JzI4QQN998s/jMZz5jyXjsZHh4WAAQW7ZsEULwdXOik58bIfi6OdGsWbPEj3/8Y75mJmA8N0LU7jVj2ysshUIBb775JtasWTPu9jVr1uDll1+2aFT2sXv3bvT09KCvrw+f//znsWfPHquHZDt79+7F0NDQuNeQ1+vFFVdcwdeQ7qWXXkJHRwfOPvts3HLLLRgeHrZ6SHUXj8cBHN9kla+b405+bgzN/ropl8v4xS9+gXQ6jdWrV/M1c4KTnxtDLV4ztt38cGRkBOVyGZ2dneNu7+zsxNDQkEWjsoeLL74Yjz/+OM4++2wcPnwY3/ve93DppZdi586daG1ttXp4tmG8TiZ6De3fv9+KIdnK2rVr8bnPfQ7z5s3D3r178Zd/+Ze4+uqr8eabbzZNx04hBNavX49PfvKTWLZsGQC+bgwTPTdAc79uduzYgdWrVyOXyyEUCuGXv/wlzj333EooaebXzGTPDVC714xtA4tBkqRxXwshTrmt2axdu7by+XnnnYfVq1dj4cKF+OlPf4r169dbODJ74mtoYjfeeGPl82XLlmHVqlWYN28enn32WXz2s5+1cGT1c9ttt+Hdd9/Fv/3bv53yvWZ/3Uz23DTz62bx4sV4++23MTY2hk2bNuHmm2/Gli1bKt9v5tfMZM/NueeeW7PXjG2nhNra2qAoyilXU4aHh09Jsc0uGAzivPPOw+7du60eiq0YlVN8DU1Nd3c35s2b1zSvo6985St45pln8OKLL2LOnDmV2/m6mfy5mUgzvW48Hg8WLVqEVatWYePGjVi+fDl+8IMf8DWDyZ+biUz3NWPbwOLxeLBy5Ups3rx53O2bN2/GpZdeatGo7Cmfz2PXrl3o7u62eii20tfXh66urnGvoUKhgC1btvA1NIHR0VEcOHCg4V9HQgjcdttteOqpp/DCCy+gr69v3Peb+XVzpudmIs3yupmIEAL5fL6pXzOTMZ6biUz7NTPjZbsm+sUvfiHcbrd45JFHxHvvvSduv/12EQwGxb59+6wemqW+/vWvi5deekns2bNHvPrqq+L3fu/3RDgcbsrnJZlMiu3bt4vt27cLAOKee+4R27dvF/v37xdCCHHXXXeJaDQqnnrqKbFjxw7xhS98QXR3d4tEImHxyM13uucmmUyKr3/96+Lll18We/fuFS+++KJYvXq1mD17dsM/N1/60pdENBoVL730khgcHKx8ZDKZyn2a9XVzpuemmV83GzZsEFu3bhV79+4V7777rvjWt74lZFkWzz//vBCieV8zQpz+uanla8bWgUUIIe6//34xb9484fF4xIUXXjiuvK5Z3XjjjaK7u1u43W7R09MjPvvZz4qdO3daPSxLvPjiiwLAKR8333yzEEIrUb3zzjtFV1eX8Hq94vLLLxc7duywdtB1crrnJpPJiDVr1oj29nbhdrvF3Llzxc033yz6+/utHrbpJnpOAIif/OQnlfs06+vmTM9NM79u/viP/7jyXtTe3i6uueaaSlgRonlfM0Kc/rmp5WtGEkKI6q7JEBEREdWXbdewEBERERkYWIiIiMj2GFiIiIjI9hhYiIiIyPYYWIiIiMj2GFiIiIjI9hhYiIiIyPYYWIiIiMj2GFiIiIjI9hhYiIiIyPYYWIiIiMj2GFiIiIjI9v4P/OmayPFCWBMAAAAASUVORK5CYII=", + "text/plain": [ + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "plt.plot(np.abs(d0))\n", + "plt.plot(np.abs(d1))\n", + "plt.xlim([0,35])" + ] + }, + { + "cell_type": "code", + "execution_count": 18, + "id": "a526b930-1c20-4f6f-a500-602451eca9d7", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([1.99374064e-04, 4.57524450e-04, 4.32687240e-04, 1.24907308e-04,\n", + " 6.33593869e-05, 4.56193763e-05, 4.75973902e-05, 5.54942255e-05,\n", + " 5.74022984e-05, 5.52140546e-05, 6.66276446e-05, 6.13030226e-05,\n", + " 6.02718391e-05, 6.69965106e-05, 6.11999090e-05, 6.76837959e-05,\n", + " 6.11215231e-05, 6.26887322e-05, 6.25784182e-05, 6.09588752e-05,\n", + " 6.45932617e-05, 6.03679218e-05, 6.35512850e-05, 6.24106123e-05,\n", + " 6.16357695e-05, 7.33561331e-05, 6.06342353e-05, 7.11609493e-05,\n", + " 7.34195697e-05, 7.50249677e-05, 9.92650223e-05, 1.28357489e-04])" + ] + }, + "execution_count": 18, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "abs(data[0,1,:])" + ] + }, + { + "cell_type": "code", + "execution_count": 19, + "id": "c6c247dd-a41c-4ebd-9c65-867c1b995dcf", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([1.02363030e-04, 1.46313849e-04, 3.11154612e-04, 2.78583974e-04,\n", + " 8.09806018e-05, 3.68793542e-05, 2.63866548e-05, 2.65659120e-05,\n", + " 3.36026736e-05, 3.41811469e-05, 2.98064237e-05, 3.58745262e-05,\n", + " 3.01800618e-05, 3.89231851e-05, 3.18832099e-05, 3.74617245e-05,\n", + " 3.47063111e-05, 4.39307486e-05, 3.74252928e-05, 3.30476321e-05,\n", + " 4.23797958e-05, 4.32881682e-05, 3.41088402e-05, 3.81523427e-05,\n", + " 3.94577691e-05, 3.87632164e-05, 3.99456672e-05, 4.02552912e-05,\n", + " 4.69936571e-05, 4.14156458e-05, 4.99788884e-05, 5.93829145e-05])" + ] + }, + "execution_count": 19, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "abs(data[0,2,:])" + ] + }, + { + "cell_type": "code", + "execution_count": 20, + "id": "2cf6296e-107b-48da-8727-7f0569eb0a40", + "metadata": {}, + "outputs": [ + { + "data": { + "text/plain": [ + "array([0.00095033, 0.00089374, 0.00028051, 0.00016192, 0.00013652,\n", + " 0.00013013, 0.00013085, 0.00013495, 0.0001451 , 0.00015001,\n", + " 0.00014344, 0.00014393, 0.00014325, 0.00014962, 0.0001415 ,\n", + " 0.00014056])" + ] + }, + "execution_count": 20, + "metadata": {}, + "output_type": "execute_result" + } + ], + "source": [ + "np.abs(data.flatten()[:16])" + ] + }, + { + "cell_type": "code", + "execution_count": null, + "id": "188053fa-48bf-40ec-a65f-63860fbf955d", + "metadata": {}, + "outputs": [], + "source": [] + } + ], + "metadata": { + "kernelspec": { + "display_name": "Python 3 (ipykernel)", + "language": "python", + "name": "python3" + }, + "language_info": { + "codemirror_mode": { + "name": "ipython", + "version": 3 + }, + "file_extension": ".py", + "mimetype": "text/x-python", + "name": "python", + "nbconvert_exporter": "python", + "pygments_lexer": "ipython3", + "version": "3.10.20" + } + }, + "nbformat": 4, + "nbformat_minor": 5 +} diff --git a/requirements.txt b/requirements.txt deleted file mode 100644 index 3af6f5a..0000000 --- a/requirements.txt +++ /dev/null @@ -1,8 +0,0 @@ -matplotlib==3.7.2 -numpy==1.25.0 -pandas==2.0.3 -pyeit==1.2.4 -pyftdi==0.55.0 -pyserial==3.5 -setuptools==78.1.1 -tqdm==4.66.3