{ "cells": [ { "cell_type": "markdown", "id": "c12c6b39", "metadata": {}, "source": [ "# Overview of data using View3D\n", "\n", "or the initial overview of what has been measured and for viewing some features in the data, a quick 3D viewer can be initiated. This tutorial shows how to initialize it but not how the interactive parts work; for this see [Viewer3D documentation](../Data/Gui.html), the [dedicated tutorial](AdvancedView3DTutorial.html), or the in-depth explanation in [Interactivity](../../InDepthDocumentation/Interactivity.html). What is done when using the viewer is that all data from all provided files are binned in 3D pixels (voxels) of equal size independent of density of measured points. This does of course not reflect the actual data quality but gives quick insights into the data." ] }, { "cell_type": "code", "execution_count": 1, "id": "5084a1d1", "metadata": { "execution": { "iopub.execute_input": "2023-07-25T11:38:23.091618Z", "iopub.status.busy": "2023-07-25T11:38:23.090616Z", "iopub.status.idle": "2023-07-25T11:38:25.130350Z", "shell.execute_reply": "2023-07-25T11:38:25.129492Z" } }, "outputs": [], "source": [ "# Import of libraries needed\n", "%matplotlib inline\n", "\n", "\n", "from MJOLNIR.Data import DataSet\n", "from MJOLNIR import _tools # Usefull tools useful across MJOLNIR\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "\n" ] }, { "cell_type": "code", "execution_count": 2, "id": "6db44adb", "metadata": { "execution": { "iopub.execute_input": "2023-07-25T11:38:25.132352Z", "iopub.status.busy": "2023-07-25T11:38:25.132352Z", "iopub.status.idle": "2023-07-25T11:38:26.097051Z", "shell.execute_reply": "2023-07-25T11:38:26.096414Z" } }, "outputs": [], "source": [ "# Loading and conversion of data\n", "\n", "files = _tools.fileListGenerator('136,137', folder='C:/Users/lass_j/Documents/CAMEA2018',\n", " year=2018)\n", "\n", "# Define the DataSet object and provide file(s)\n", "ds = DataSet.DataSet(dataFiles=files)\n", "\n", "# Run the converter. If no binning is specified MJOLNIR uses highest possible binning (currently 8)\n", "ds.convertDataFile()\n" ] }, { "cell_type": "markdown", "id": "37a0651c", "metadata": {}, "source": [ "The viewer is an interactive widget which takes all of the data from the data files and bins them into voxel (3D pixels). One can view data from CAMEA as being in a 3D space; either with coordinates ($Q_x$, $Q_y$, $E$) or in the scattering plane ($H$, $K$, $L$, $E$). The binning makes sure that the x-axis of this cube of data is either along projection vector 1 (rlu = True) or $Q_x$ (rlu = False), with the y-axis 90$^\\circ$ to the x-axis and the energy axis out of plane. \n", "\n", "The inputs for the method are\n", "\n", "- dx: bin size along x-axis ($Q_x$ or projection vector 1)\n", "- dy: bin size along y-axis ($Q_y$ or orthogonal to projection vector 1)\n", "- dz: bin size along z-axis (Energy)\n", "\n", "A couple of additional parameters can be set, e.g. if a grid should be plotted." ] }, { "cell_type": "code", "execution_count": 3, "id": "5eedec07", "metadata": { "execution": { "iopub.execute_input": "2023-07-25T11:38:26.099083Z", "iopub.status.busy": "2023-07-25T11:38:26.099083Z", "iopub.status.idle": "2023-07-25T11:38:26.785575Z", "shell.execute_reply": "2023-07-25T11:38:26.785458Z" } }, "outputs": [], "source": [ "Viewer = ds.View3D(0.03,0.03,0.05, grid=True)" ] }, { "cell_type": "code", "execution_count": 4, "id": "18c4ad08", "metadata": { "execution": { "iopub.execute_input": "2023-07-25T11:38:26.788577Z", "iopub.status.busy": "2023-07-25T11:38:26.788577Z", "iopub.status.idle": "2023-07-25T11:38:27.168990Z", "shell.execute_reply": "2023-07-25T11:38:27.168196Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAs4AAAG8CAYAAADO7iN7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAsTAAALEwEAmpwYAABj30lEQVR4nO3de3xcd3kn/s8zo/vd1s26WJIlWfcQSAJxiKH+8WsKoUmbpWUv5VKg927bALttgYWlu79lSylpu6XtUggJhW63LVDaTSAQLtFadi7ghJDo4rt8kxxfZEuyZFuyZp7fHxrZc56R5szIZ3Tm8nm/XvOKj86Zc565ZPTVOZ95vqKqICIiIiKi+AJ+F0BERERElAk4cCYiIiIiSgAHzkRERERECeDAmYiIiIgoARw4ExERERElgANnIiIiIqIEcOBMREREWU1E7hSRwXXc7z0iMhi5PSsiV0WkyvsKKVMI+zgTERFRthKR3wPwLgDzqrrjJvbzlwB+rKqf86w4yjg840xERETZ7AiAt60siMgtIvJU5Czy10Sk0m0HInIHgH4OmokDZyIiIspaqvo1ANeifvR5AP9eVXcB+CaA3xORt4jIsLn9bNR9PgLgv2xc1ZSu8vwugIiIiGgD9QL4KxEBgHwAh1T1WwC+tdrGkUxzt6o+tWEVUtriwJmIiIhyyQEA71bVEyJyN4AGl+3fCOB7qS+LMgEHzkRERJRLfgPAl0QkD4AC+CWX7bsBHE15VZQR2FWDiIiIyGMi8gKA2cjiuKq+1896yBv8ciDRBlurn6iIBETksyLyTOTb3p0+lJe24vVhFZEXonqtPrrBpaUlEckXkS+LyJCI/EBEfsas5/uNKEVEpAjLJyd3RW4cNGcJRjWINlB0P9FVVj8AoEhV7xKRHQAeAvCzq2yXc+I9b9G/oDa6rjT3TgBTqvouEdkM4EUA/ydq/QPg+40oVW4FUCIiT2J5rPURVX3W55rIAxw4E22slX6iX15l3U5EvtWtqs9G+obSsnjPG39Bre4rAL4a+bcAWDLr+X6jrPbm/6dUpy6EUrLv519aGAFwNepHnzM9ni8D+DSAhwFsB/CEiHSrqv3/kDIMB85EG0hVvyYibWusrgAwE7UcEpE8ftC6Pm/8BbUKVZ0DABEpx/IA+qNmE77fKKudvxDCc99uTsm+8xuOXFXVeH9sHgRwWJe/SHZQRKaw3L3jZEoKog3DgTNR+pgFUB61HOAgJiH8BbUGEdkK4OsA/kpV/86s5vuNspwipGG/Dv4+ALcA+E0RacTyH6qn/SqGvMMvBxKlj70A3goAkczpy/6WkzHeh+V8LvgL6gYRqQfwJIDfV9VHVtmE7zei1PkCgCoR2QPgHwC8j3+YZgeecSbyUeRLWw+r6tuwfGbwHhF5GsuZVH4Lew3mefsCgC9GfkEp+AtqxUcAbALwMRH5WORn7wDwGb7fKBcogDD8abmrqosAfsGXg1NKsY8zERERZZ3bbi3UoW9tScm+yxpPPO+ScaYsxTPORERElHUUihBPDpLHmHEmIiIiIkoAzzgTERFRVvIr40zZiwNnIiIiyjoKIMSBM3mMA2fyXE1Njba1tfldRkY5efIkGhoakJfH/yUTpao4deoUn7ckrTxvjY2NCAaDvtby/PPPn1fVWl+LICJKAn/bkOdqa2uxb98+v8vIKA899BBeeukl/M3f/I3fpWSUT33qU9i/fz8eeWS1NsW0lk984hM4duwYPv/5z/tah4gc97UAynqMapDX+OVA8tz09DSOHDnidxkZ5cEHH8Tzzz+P559/3u9SMsoHP/hBPPfcc3jxxRf9LiWj/O7v/i6efvppDA8P+10KUcoogJBqSm6UuzhwJs8VFBTgox/9KMJh36Y6zTh5eXn4sz/7M+zevRvsrZ64vLw8PPTQQ9i9e7ffpWSUgoICfPKTn8TQ0JDfpRClVDhFN8pdHDiT52prazE2NoannnrK71Iyyk/+5E/iTW96EyYmJvwuJaO85S1vwV133YWzZ8/6XUpGue+++3Drrbdienra71KIiDIGB86UEp/5zGfwwgsvYHFx0e9SMkpvby8OHTqEpSXOGJ2M/v5+jI6O8ipHEkQEAwMDGB4e5lUOykoKRShFN8pdHDhTSrzhDW/Av/t3/w779+/3u5SMUlBQgNbWVhw+fNjvUjJKSUkJtmzZgmPHjvldSkapqKhAeXk5r3JQdlIglKIb5S4OnCllmpqaMDs7i9nZWb9LyShtbW04c+YMLl++7HcpGaWzsxPHjx/nVY4k8SoHZSsFM87kPQ6cKWVEBP39/bwUnKRAIIDe3l6Mjo76XUpGycvLw/bt23mVI0m8ykFElDgOnCmlNm3ahOLiYrzyyit+l5JR6urqEAqFMDU15XcpGYVXOdaHVzkoOwlCKbpR7uLAmVKut7cXBw4cQCgU8ruUjNLf34+RkRGerU8Cr3KsD69yUDZSAGFNzY1yFwfOlHJFRUVoamrC0aNH/S4lo5SVlaG6uhonTpzwu5SMwqsc68OrHERE7jhwpg3R3t6OiYkJXL161e9SMkpXVxeOHj2Ka9eu+V1KRuFVjvXhVQ7KNoxqkNc4cKYNEQwG0d3djbGxMb9LySj5+flob2/HwYMH/S4lo/Aqx/rwKgdlEwUHzuQ9Dpxpw2zZsgVXrlzBxYsX/S4lo7S0tGBqagpzc3N+l5JReJVjfXiVg7JJWCUlN8pdHDjThlmZqYyXgpOz8oW3kZERv0vJKCtXOfiFt+TwKgcR0do4cKYNVVFRgYqKCs5UlqTq6moEg0GcPXvW71IyypYtW7CwsMCrHEniVQ7KBoxqUCpw4EwbrqenhzOVrUNfXx/GxsYQDnPeqkRFn63nVY7E8SoHZQOFIIRASm6Uu/jq04Zbmans0KFDfpeSUUpKSlBfX4/x8XG/S8kovMqxPrzKQdmAGWfyWp7fBVBuamtrw9DQEFpbW1FSUuJ3ORmjs7MTQ0NDaG5uRmFh4YYe+y2bf8WxHJqedizndXc6lpcOOKdwDhQUXP+3FBc79zUzE/fYAbN9+MoV53rzHgrf5Ax43wl/5abuny36+/vxgx/8ADU1NQgEeJ6FiIifhOSLlZnKeCk4OXl5edi+fTv279/vdymUA4qLi1FfX49jx475XQpR0phxplTgwJl8U1dXh3A4jPPnz/tdSkZpamrCpUuXMDs763cplAM6Oztx/PhxLC4u+l0KUZIEIQ2k5Ea5i68++aq/vx+jo6P84lYSVtr6DQ8P83mjlFu5ysHJi4iImHEmn63MVHb8+HG0tbX5XU7GqKqqQnFxMU6fPo3GxsZVt7kn8HbHcl7DFsdy6NxU3GPoknMCDJszDpaXO/d32PmlxUBhkdnhjW4gNtMcnX8GADH5bTVTZ9ta4NKhRfLyncsB56XWsDmbap+7eIJlZY7lkGnhlg156aamJhw7dgyzs7OoqKjwuxyihCiAMM8Pksf4jiLfdXV1YXx8nDOVJam3txcHDx5EyAwqibzGqxyUqZhxToyI3CEiX/W7jkzAgTP5bmWmsgMHDvhdSkYpKipCU1MTjhw54ncplAOir3IQZQJVZpwTpar7VPXn/a4jE2Tfq08ZqaWlBRcvXsSlS5f8LiWjtLe3Y3JyEldMezaiVOBVDqLMIyK7ROQZEfmaiOwXkRdE5H4R+Y6InBCRP41sMxzZ/osi8uci8pSIHBaRx0WkzO04PtUdEJH/ISLPicioiIyJyN2Rn39PRD4V2ddPisgpEam/2bqYcaa0ED1T2Y4dO/wuJ22tnb39SwQrKx0/sblevZTc9Mn2/oFNVY7l0FlnNxSbQ7bLEgyueSwNOy//h80fUDYvHXOspfgxH5tpRp7zoy8gznMIUuw8nu1Z7VhnMs221mTy0kB6Z6Kjr3J0dXX5XQ6Rq3AWxirW6bUAXquqPxKRJwB8GMAuABUAJgH80Gx/O4A3AQgDeA7A2wE8umHV3uBW91cBNAK4S1XDIvIhAB9S1ftF5J0AfiQiewH8BYBfUNUzN1sQB86UNjZv3oy8vDycOXMG9fU3/UchEaVAe3s7hoaG0NLSgqKiIvc7EPlkuY8zL6xHjKvqjyL/PgJgRlUXAZwXkVkAm83231LVBQAQkZdXWb9R3OqeAfBRAL8mIh1YHlRfAgBVPS0ivwLgXwB8XFV3e1EQ31GUVvr7+zE2NoZwOOy+MRFtuGAwiJ6eHoyOjvpdChElbsEsu30bPzr/p4Bvp+7d6v5/AXwj8u9/AfBZOGvtB3AGwOu8KogDZ0orxcXFaGhowPj4uPvGROSL+vp6LCws4OLFi36XQhQHvxyYA+4H8Jiq/k8sx00eABAEABF5HYAHAdwBoEpEHvTigIxqUNrp6OjA0NAQmpubUWj6+WYbm3+1PYH1mrM/cXRO2GaGw3PxM8w2i2vZ3shh84XDpUlnNwV7/JhMtckdR+eS7bZu3DLNNlccXrjq3IHJNIcvX3Ysx9Rj1idFnVdLYrLipSWOZdvTOhMy0SvfSXjppZdw9913Q4Q5Uko/7OOcE94P4Esi8hKAEIDdAH5ORCoB/G8Av62qEyLyHgA/EJHdUdGPdeHAmdJOXl4eurq6sH//ftx6661+l0NEq6ioqEBFRQVOnTqFrVu3+l0O0apCyj/qVHUQwEDU8m+Z9TWRf/5FZPk9Zr1jeaMkUfdt5q4rZ5Y7orY9AqDai7r4pxilpcbGRszNzWHGnI0jovTR09ODw4cPY8ll5kYiomzBgTOlpZVLwZypjCh9FRQUoK2tDYcOHfK7FKIYCkEIgZTcKHcxqkFpq6qqCqWlpTh9+jQaGxv9LmddbF7VZnFjcspXF+Kuj2Z7H8esN7lgm2GGWW8zzbZWBM0vC3OWMby4GLee6KyvzfnafLZb3tqymeaYvPVictO5x3ssdt8xPaJNT+hgZaljeWlqKu6x3R57oKDAseyWiU51Brq1tRVDQ0NobW1FSUmJ+x2INlCYX+TLSoODgx2BQOCDAN4ZDofLA4HAJQB/Gw6H/2TXrl0pnU6X7yhKa5ypjCi9BQIB9PX1YWRkxO9SiBxW+jjzjHN2GRwcvDcQCLzU1NT0K6997Wsr3vjGN8prX/vaiqampl8JBAIvDQ4O3pvK4/PVp7RWWFiI5uZmHDmS0j8giegm1NbWIhwO4/z58+4bExGtU+RM81dvvfXWko6Ojvzi4mIEAgEUFxejo6Mj/9Zbby0JBAJfHRwc7HDf2/pw4Expr729HZOTk7jicrmeiPwzMDCA0dFRfieB0oZCENLU3MgfgUDgg01NTfmVlZWrrq+srERTU1N+IBD4QKpqYMaZ0l4gELg+U9ntt9/udzkOb9n0y45lvWonOTJMHtbmkCXP+bdsTD9ihyTjKy6ZaEsKTG/jfOfHRehC/F7HwfJy5/aXLt34t0u3FLH57CTz1jb/bfs+59XWOJaXziV+ptRmmm0e2maQXTPNtge1zTQnmfe24mWgvcw/l5aWoqamBsePH0dbW5tn+yW6GezjnHXe2djYGHcigMbGxvyJiYl3AfiteNutFwfOlBHq6+sxPj6OCxcuYPPmzX6XQ0Sr2L59O/bs2YOmpibk5yc3yQ0ReW94eNjTS0ADAwMberq97UPfcNT/6JtLXSdGKywsRCgUrrD3PfbJn/akdv4pRhlBRDAwMICRkRFeCiZKU/n5+ejo6MCBAwf8LoUIquCU21lmUQULC/Gv7C4sLGAhhXEavvqUMcrLy1FVVYVTp075XQoRrWHr1q24ePEiLkVFc4j8IQin6Eb+eHriGk6cmoi7zfGTE3h6IrkWpMlgVIMySnd3N/bu3YuGhgbk5aX+7evWI9f27HUTk181fW9tv+F4fZxtP+GY/LTJ3sbPS8ceK6b38UL8HtM2rx2+4jxe9GMNX46fj7a54Zi+zOax2JywLsV/rMlkmlfbv2OdeQ3FRBQCEj+3bpfdMs0xz4XJb9uMtWNf5nl9c9kvOteb12U9GeiVyYtGRkawY8eOpO9P5BUFeHY4y3zz6CJ2Nk9iS10tVvuC4MzMDE5OTOKJo/HnFbgZfEdRRuFMZUTpb/PmzcjPz8eZM2f8LoWIssi5K4o/f/4yfvDCj7H/4GFcuXIF4XAYV65cwf6Dh/GDF36MP3/+Ms5dSV2kk2ecKeOszFTW0tKC0tJS9zsQ0Ybr6+vDc889h9raWgQCPEdD/uBkJdnn5fMhfGT3HO5tP4bXn5pAoSgWVPD0xDU8cXQxpYNmgANnykDRM5W97nWv87scIlpFcXExGhoaMD4+jo6OlM1FQLQmhSDMnstZ6dwVxZdGFvClEZcWsCnAgTNlpNraWoyPj+PcuXOora1d935shjlYVuZYds2T5sf/X8j2E47pAeyS9U2GzRgHq6qcx5qbj1/LoksmzOa5zbKsHcde3n/UY43X43k97GOPl/MFADE5YttX2r7u8fLh9j3j1qPajc00B02Oz+7fZqJj+nVreO1jubz/7P8fyWaeOzs7sXv3bjQ3N7u2kCJKBZ5xJq/xHUUZq7+/nzOVEaWxYDCIrq4u7N+/3+9SiIg8wYEzZazS0lLU1tbi+PHjfpdCRGtobGzE3NwcZm7yTDxRshRAWAMpuVHu4qtPGW379u0YHx/HolvMgIh8sTJ50fDwMK8O0QYThFJ0o9zFjDNltOiZym655Zab3l9obs6xHJMfjZMXXV4df2Bgs7ixOeH4QWGJylS79fcNTU/HXS8mByw2j23y227HC1Y7p0IPXXAeP1pMj2gXgVJnr+SYTLJLXtvmve1zE9MH2r5Ocdh+2TZf7ZodN2wtbplp+7rE7C9O3jtef2oAkALna+zW13ytDHRlZSVKS0tx+vRpNDY2xt0HEVE64xlnynhbt27F9PQ0ZyojSmO9vb04ePAgQkn8UUJ0MxjVoFTgq08Zb2WmsuHhYb9LIaI1FBYWorm5GYcPH/a7FMohjGqQ1zhwpqywefNmFBQU4JVXXvG7FCJaQ3t7O06fPo0rLvESIi+oCs84k+eYcaassTJTWV1d3bpnKovp22zyqTYXHLM+if6/gHseNibzLGuf6bA5X5j72swyTH5Vws78ts0V67Ul57LpaR2auuBcb3srR+Vp7fNiH6fNlsdkms329nkLlDhrtz2sY7K9Jrtus72OTc3jTjbD7MY+N/a5kDzn6+iaF496bPZxxxzLvB9tf+3g5k3O9Rcuxj+2EQgE0Nvbi0ceeQS33XZbUvclIkoHHDhT1liZqezo0aPo7Oz0uxwiWkV9fT1qa2vx4IMP+l0K5YAQzw6Tx/iOoqzS2dmJU6dOYWFh46fhJKLE3HvvvSvTcFf7XQtlLwUQhqTkligRqRORkyLSk7pHShuJZ5wpq6zMVDY2NoZXv/rVSd8/UOScFtjGE1xbf5nL6sGC+FNLu13mt3EHjTNFsmtcwd7XPBZbe8wU3SbqESgrde7P1Bq+4owBuMVWosW0eDOPJTaeYDo1LMWPldhly55RcNvecV+XOITb9m5udnryeMd2qzXZaMZaysvL8du//dv4+7//e568oRQSX884i0g+gL8GwFB/FuGHFmWdhoYGzM/Pc6YyojR2xx13AMA5v+sgSqFPA/gsgEm/CyHvcOBMWYczlRGlv4I4E7MQeWG5j7Ok5AagRkT2Rd1+NfrYIvIeAOdU9ds+PHRKIUY1KCtVVlairKwMk5OTaGpq8rscIiLyQSh15wfPq+odcda/D4CKyE8CeDWAL4nIz6gqe6ZmOA6c6ToRCQD4KwC3AlgA8MuqejjR9emmp6cHTz/9NLZs2YKgy1TWa3HLtsZMcz3vzBHbtmkxyy7t7eyU3I62aXa6btOOTs0XJN1aidk8t8002/WBWuf3uvSiMxoTKDbTWC/cqM8t2x2z3iXzHHusBbPeTJ1uW+VtqXMe/5WzceuLV4t9jYMmN29b6yWT/V71+OaxxbyOgbWfd7dj523vcCwvHTriWLbvqZuVbZ9BlLtU9Y0r/xaRQQC/zkFzdmBUg6I9AKBIVe8C8CEADyW5Pq0UFhaipaWFM5URZY4HkEWfQeQvRWpiGpGoBuUoDpwp2k4A3wIAVX0WgL0M5bY+7Wzbto0zlRFljqz7DCJ/hRFIyS0ZqrpLVfen6CHSBmNUg6JVAIi+jhwSkTxVXUpwPQDg1KlT6O7uvr5833334f77709Vza6Wlpbw1FNPocRcUicibw0ODrpu89hjj+Hxxx9fWawxqz35DCICAFUgxLPD5DEOnCnaLIDoxsMB8wvJbT0AoLm5Gfv27UtRievz7LPPoqurC5s3b3b8/BP4S8dyTD7UZJiDDVscy6HTzshazBTZhl0fk2lOgs1f27y0XQ7PzsWvzU7JbQSbnI99eofzS5fXSrY6lqtfNFneH41e/7fNBVu2X3XM1ORx+lmvh8w4eyNrOPFuLMnWYt9TNpvulv+O6bdt+3HHmcbdbYp3KybTXFnpWBaTs9+1a1fc/a1s89BDD0XuL+fNak8+g4iIUoVRDYq2F8BbAUBEdgB4Ocn1aau/vx8jIyNsT0eU3rL2M4j8wYwzeY1nnCna1wHcIyJPAxAA7xWRzQAeVtW3rbbev1KTU15ejk2bNuHkyZNoaWnxuxwiWl3WfgbRxlv+ciDPD5K3OHCm61Q1DODXV1n1Npf1GaG7uxt79uxBQ0MD8vPz3e9ARBsq2z+DaOOFwLPD5C0OnCln5OfnY9u2bTh06BD6+vpW3UZMz1296uwHbDPNavoBS57z7EawssJ5/+npuDXa/QUKnQN8XQojUTb7Giwrc+5r0fR9tn2ae9ody0vlzl7Je//0r+Mev+fh33Astx+/0fPX9pC2uV+bxbWZ45isuH3ekszy2t7KN8M+z6E5Z7Y8UGFehyvxeynHPBaXDjH2sQYKi9bYMnlePk9ERJmI1zAop7S2tuL8+fOYn5/3uxQiIkqhFE+5TTmKA2fKKSKCvr4+DA8P+10KERGl1HLGORU3yl189Snn1NTUIBAI4Ny5c36XQkRERBmEGWfKSf39/fjhD38Y83O1+VGXvsy2H7EuOnsr2x6/Npvrlr0NLzjzr9H3j+kHbPs4m1ywzTDbY8fs77Szxe74O7uRjG+8548dy7/9Z/dd/7dbL2L7WGzPalfmebaPze7P5oiR5/xojH5dY3p1B5K7bGtz826Z5WR6Sq8muj93TN46ySw4UaYJ88uB5DEOnCknlZSUoK6uzu8yiIgoRThzIKUCB86Us7Zv3+53CURElELMI5PX+I6inJWXx78biYiIKHEcORBFkUJnH+fwlbUzxkBshtkKlJoMtMmzuvUrjpftjellHHb2eA4UO/v3hi5dcizH9HU2GWg1edet3zX513evWRoAoCPPuf/w7I18rX1csY8lfq7XLbMM21/b5JBj9i7OcwhiJsiJfp1trYGqSuehz005loNVVc5jm/eMW9/npPPdht2fY98uz3NeU6NjeWli8qZqIdpIyzMHMqpB3uLAmYiIiLISvxxIXmNUg4iIiIgoATzjTERERFlnZeZAIi9x4EwURU02NlBh8qcXZ5zrC505YtsHV69di7vestndYHm58/hROeWYWk1Pabf8tc00x9RiMtJFp5yP/e4P/prz+IvOvGzZcee05hI4eOPY9lgmWw51yXqryXOb18Gut897sNKZSw7PO5+r8Ez81zmazTTbPHVoetq5L9szWuL/YnfrtWx7YsPklm0v8Gj2eQ1u3uRYtpnmeO9HonTErhrkNQ6ciYiIKPsovxxI3uOfYkRERERECeAZZyIiIso6CnbVIO9x4Ew5zfbkxZLpZXzFmRcNFJksrsmTxvQLXnRmSG0eNWz6OtvMqS449xddr80420yzPZYkkX0FAL1gMtSXnbVWmbw3Ck1P6lln/+Bw1PFj+jYvLDiXbR9ml77NbtvH3N/leY/JIZvMdNxa1txyWUzO3WaWbW7evE7J5oyjt3fbNnThYlLHymvdGnd/RH5jVIO8xoEzERERZR121aBUYMaZiIiIiCgBPONMREREWYlnnMlrHDhTzrkn8Pbr/3btgWuysLZXsoacmWibpbX51JgsrUueNaZnb/T24rxgFJO1NfcVm9s17HMRKCt1LC+dO+9YDhaYHLHZnz1+3FrM824zyzH9hstMf+05Z5469njmdTF58JhMtHkuQlF9ne2x7WW7ePcF3Psyw+TkbW3J9k5OZntbm9t9l46fTKoWoo2kYDs68h6jGkRERERECeAZZyIiIspKbEdHXuPAmYiIiLKPMuNM3uPAmXJOdGZU8pz/C7jmQUPOfr6BYmdG2d4/pl+xS3bXrV+x5N+o1+ah7bGC1Zsdy0tnzznXV1Y6azdZXJvvjmGei/CVS2bZef/ox6bXnNnweH2SVz20S6Y5Jktu+2+7bB+ev7zGlu7Htn2ZLZtpDm7e5Kxlbj5ubcmK6UkdvW/bh9xm/k2ted2djuWlA4dvqjaiVGI7OkoFZpyJiIiIiBLAM85ERESUlXjGmbzGgTMRERFlHbajo1TgwJmyXnTfZiA2C+xYZzLGMdlbsxyTabb9gE1PXptxDlZVOfc3Pb1mbYAzG2x7SkOdedXQBee+YvoBm0yzFfNcGDYHbPPewfJyZ3lLN2qP6fFselIHy0wvZLdcsQu3+9se1jaDHb3errOvcUyG2eU1tpnmmNpcen279YWO6RMdj3l/Wjp5xrH8nfBXEt83kQ+UA2fyGDPOREREREQJ4BlnIiIiykrs40xe48CZiIiIso6yjzOlAAfOlNPCl505XZtBttlbDZt8q81Lmwy0xo+MAi79hS1HfabWQFmpc1ub0y11ZqJtxtktS2vZXLDNe9v9xTy3UWKy4CaX65ZBDlY7eyEvnTm75rFWY3tOx+19bHtrx9kWcM+t2wxyTGbZ5XWw71ErWFZ2oxaT9bavkVsW3LXPORFRluPAmYiIiLISvxxIXuPAmYiIiLIQ29GR9zhwpqxj28/Zy9HRcYqYy+RucQUb3XBp32Uvu9sWb0sTk/FrNaQgatrqhQVnrS6X0dVOgW1jJiYuYWu30zPHm1IbWGVK8Ogpt83zFlOLnarcvk5R8QPAPZqRbAwlqRZuLuzz4haZSfbYbu/BePELMS0EYZ4X21KQUQ3KNDzjTF5jOzoiIiIiogTwjDMRERFlHQW7apD3OHAmIiKi7KMxE6oS3TQOnCnrBCsrHcs2Qxqdp7XTVuuis9VYTBs0kz+17eZsVjdm+mVTm5jWZgia9FTILC8527A59mUzyaZlm80N2/Z14VlnFjamRZyt1bDr404Jbqcyt9lxO621eWz2NQ2azLN9nVxb67lMWx29Ptkpre3zYqcqd6vFssez7zm318khzvsJiM0051VXJ75vojTACVDIa8w4ExERERElgGeciYiIKOso2FWDvMeBMxEREWUh9nEm73HgTBkvpm+zzdYa0T2AbVbWcsuLumVjbf40ZpprW6vJnMZM7xyduTb9e21fZZunttuHLlyMu73N4tp+xJbNRNvpzKPZx223tdnymCmx7XTfpldxTO9klym7pdA59XnATmNtM9lx5NXVOpbDc/PO5TjPy/KhnN9mss/rKndIuDb7fo2ZgtvldUl2ingiomzDgTMRERFlJXbVIK9x4ExERERZiRln8hoHzkRERJR1VDlwJu9x4EwZx2aaLV1YcCwHCouc66861zu2dcksxzAZ5pieuiZXbEnAmaUNm+1tPdH9id327ZaBttSsj+nL7JI7tr2Wbc7Y8dyYLLd9nJLn8tEUjn/9NSZvbTPStn+3ec/Y1xVyox63Hp5LZ885j2WeN7f3mH3ePb3SbLPbRkzW3DxPSxOTXlZDRJRxOHAmIiKirMSuGuQ1DpyJiIgoK/HLgeQ1DpyJiIgoKzHjTF7jwJnSXkyfZptZdum1HF64uuY6m8O1PXRtb2Pbh1ld+jbbHLJrH13z2BB0ZlKjj2ePZXsRx+zbZG1hehm7sb2PY9hc8qLzdXG8Dq9/tWNdcP9x57azztqDUb23gVWed/M6umXT3XopW9H7j+mrbPPQ9lg2X21e46Rz9Xb/SWxv/19w69uc7PNERJTtOHAmIiKirKMQnnEmz3HgTERERFmJEWfymltnJSIiIiIiAs84Uxp68+0fdyzbTGi8zPJqYjLRUbnjmH6/tqeu7fls+uDa3sU2B2xzyHb/bplotZnU6Mdis7ZL8Y9tBcpKHcuh6em427uxvZDj9ZkOLDhrnfqZHsfy5h+bLPmBY47lmOfRZJxtbthyywUHy8sdy6FLl24cy2aSbT9rm5tfipP1ToDNptv3mFvG35EPF+d7JvpxrSavZ7tj+YnR/x53e6K0wglQKAU4cCYiIqLsxKwGeYwDZyIiIspKPONMXmPGmYiIiMhDIhIUkUdEZK+I7BGRAb9rIm/wjDNdJyIBAH8F4FYACwB+WVUPm21eADAbWRxX1ffe7HFtn2Y3Nr8aMHlUt6xtIP9Gv2PbDzj2YM6zFYF85/8y9v72WG4545jaSgrNBuE1t9cFZ7bVZmtjsrb2sdpju2TJbZbW9lYOmzy4mIceXXv4+WHHuprJLc5t5+bNfZ3Pg2Vrs1dnbWZZTL9u+zrZ7G/0c+OWUXbLHLuJyWebWu3+8xobHMtLk6ed2yfZF9rh3IX133cd/PoMouzl48yB9y8fX+8WkV0APgHgZ32rhjzDgTNFewBAkareJSI7ADyEqP/RRaQIgKjqLn/KI6Is9wD4GUQeUfgX1VDVfxaRxyOLrQCmfSmEPMeoBkXbCeBbAKCqzwK4w6y/FUCJiDwpIt+P/GIjIvIKP4PIOwpAJTU3oEZE9kXdfjXm8KpLIvI3AD4D4H9t8KOnFOEZ5ywiIvF7cAFQ1XjXbSsARPcBC4lInqqu9L+6DODTAB4GsB3AEyLSHbUeAHDq1Cl0d3dfX77vvvtw//33J/goiChdDQ4O3vQ+HnvsMTz++MqJONSY1Z58BhFtgPOqav+wi6Gqvygivw/gORHpU9V5t/tQeuPAObvMADgNwF6b0sjP6gCU2jtFmQUQHQYNmF9IBwEcVlUFcFBEpgA0ADgZvZPm5mbs27cv4aI/Wf4lx7KUOPvWLp0561iO6cFr8qhht17Kpu9uvG3d8qE2oxysqnIs297Idv8xxzP5bNvDNzQ3d+NYNmNsHlfYJQNtn0d7rJjMtM0BR9WSiLyW5uv/XjpxyrFu6fQrce/r1pfZtXfyFdMPu9iZ57aZ5mBlpXP9jLOvtGNfJSXOWvKdtdj72scitk+zyYrbTLXNottMc8xzFdV73O4rXr9qAFiamnIs79q1Czdr165deOihh5ZLEzlvVnvyGUS0wq+Ms4i8C0Czqv4hlv/gC0dulOEY1cgu31XVdlXdZm7tqroNwPdd7r8XwFsBIHIJ9GWz/n1YzhxCRBqxfHboNIiIvMHPIPKWpujm7p8AvEZEdgP4NoD3q+raZ20oY/CMc3b59yLSEvm3AriiqtfP6KiqW17i6wDuEZGnsXyG+r0ishnAw6r6NgBfAPBFEdkT2f/7VrtEGg7zj2oiWhdPPoOI/BaJZPxrv+sg73HgnF3+wSyXRXLP71bV59zurKphAL++yqq3RdYvAvgFt/1ciROFICJai1efQUTLhBOgkOc4cM4iqnqX/ZmIdAB4FMAbN6qO6elpHD58GJ2dnQlt/62ZR+Kut32eYzLLLpnmeGIyxyYfajPOUmj6LC+anr1L8WuJyTTbHr0mMx1Tb1SWV8Nr93he7ViBMhNvN7XG7M/2Rl6Kn5m26232NzR5BmuxeW2bn47JY5t9w/R5DpSaY5ucsS7E/2UaL9NsqXkPhC9fdtZis+M2K24z0GZ7t7y12/7j5fTd3m9EGY9TbpPHmHHOcqp6BBv80VFQUICPfvSjjGwQEZF/dLmPcypulLs4cM5yIhIEUOm6oYdqa2uxf/9+PPXUUxt5WCIiIqKUYlQji6zSgL0QwM8A+OeNruUv/uIv8Mwzz2Dnzp0otPEGIiKijcCoBnmMA+fs0mCWrwD4I1X97kYXsnPnTrS3t2NsbAyvfvWrb2pf3wl/xbFsM8+2d3JMptM08ozu4Wv7Jsfr8QwgJtMs+c7/hWw21613cmwvZme8JWx6+kbniF2z3OK8oKQme2vzrXY5uu8yEJtRto/dZpxt7jhmffSmSfaEtrnimHx1IP6l1Jhe4B6KyV8b9rHavsz2PRgQ52OJ6QNtlmMy0FH7d6vNsv/vEWUexirIWxw4ZxFV/S+r/VxE/o2q2o4bKdfQ0IDx8XFMT0+jygxuiYiIiDINM8654T/4cVARwcDAAEZGRqB+Td9ERES5y78JUChLceCcG3y7VlVZWYmysjKcPs3JvYiIaINx4EweY1QjN/j6v3lPTw+eeeYZ1NfXI5hEj+W1uOUubQbais4GJ93HVk2LPZM/jenbbLK29oWwedaYTPTC1TX3H6iqcG57xmxrj20fq8lA51VXOZaXTpyKW1u8/sDAKrVH96COk3dejVsvY0uvxZ9MzvY+tq9DcPMm5/EuXFz7WC79rt3Y58lmnm1WPCaf7fI6SFSmH+Y9ELp0KbEiiTKRAmDrOPIYB85ZREROY/mjQnBjjCYANvtWFIDCwkJs3boVhw8fRnd3t5+lEBEREa0boxpZRFUbVLUx6r+NqtoAoMj1zim2bds2nD59mtNxExHRhlFNzY1yFwfOWUREKkXkgyLyXpHl6/AicguAvT6XhkAggN7eXoyOjvpdChER5QpmnMljjGpkl68A2AfgNQC2isgZAB8H8B99rSqivr4e4+PjuHDhAjZvTl165GYz0NFceyXbXsUuvZFj9m9zw0lsr5dNv1/Tz9ct0xyTzTVZ2ZjazPZuz028HLHrvs16t0yzlWy23IqXab5ZedXVjmX7OoWmpx3LNt+dV+3MXy9Nnjbrnftfmppau5amRue2E5NrbkuUkZhxJo9x4JxdylX1IyIiAA4AOAbg1ap61t+ybujv78eLL76InTt3QoQfaERERJQ5GNXILlcBQJebJl8B8DPpNGgGgPLycmzatAknT570uxQicnGJXTcow4mm5ka5iwPn7BL9v/OUql5dc0sfdXd348iRI7h2Lbm2XUS0sYaHhzl5EWWuVOWb+b9ETmNUI7vcLiJPY7kFXV/Uv1VVX+9vaTfk5+dj27ZtOHToEPr6+jb8+NEZaJt3DpSUOJbDly/HXW8lmwuWUpMDNtnaoJmq3JF/XUju76K8WpONPXc+7vaSbz4eYnpYO//uDpY5H0tobm7NfdvnJaZ3cTB+Httye92SZeuJ7rWc7LHc+l/rkrPntGu+2yXvHS/THLOtyTTb7wc8//zzOHPmDLZs2ZLwPomIshkHztnlVX4XkKjW1lYMDQ1hfn4epaWlfpdDRKvo6+vDc889h7q6OgQCvEBJmUb45UDyHD8Js8tfq+rx1W4AICLf8rvAFSKCvr4+jIyM+F0KEa2huLgYDQ0NOHr0qN+lEK0PoxrkMZ5xzi5vEpG1+kkJgE1rrPNFTU0NxsfHce7cOdTW1vpSg700baMbMdEK2/LNXqY3cYaYCIJpk2bbntnj2dZk0S3ebEs1t2mkQxedl/jjxkAQGxmwbHwiNLf+zHpMGz+3aahtlGMx/rGTnbI7Zhrs6Nd9Kf503pZb276Y19E8Vtf7JxHZCZaXO5YTmXK7s7MTu3fvxtatW1FYWJjwsYjSAge55DGecc4iqloQNWOgvTWoqu8zCFr9/f0YHR1FOBx235iINlwwGER3dzfGxsb8LoWIyHccOJOvSkpKUFdXh+PHj/tdChGtoaGhAfPz85hJciIaIt8xqkEe48CZfLd9+3YcO3YMi6bbABGlBxHBwMAA29NRZlEsfzkwFTfKWcw4ZyERqQZQCWBaVS/4XY+bvLw8dHZ24sCBA7jlllt8rcUt8xyTaTYt2QIlJgNq86kuU2q7Tbmt127ka20G2WZl3XK9bplmm60NL8TPEbtlceM9NtcpspNsvRdzbJee4fa50qsLax7fTm1uxWTDTR7b5ort/tweq93eLbds169XZWUlysrKMDk5iaamJk/2SZRqnKyEvMYzzllERF4rIj8A8CSARwF8V0R+JCJp08N5Lc3NzZiensbs7KzfpRDRGnp6enDw4EGEXP7AIyLKVhw4Z5c/BfBzqnq7qv6Eqt4G4GcjP09rK5eCR0ZGeCmYKE0VFhaipaUFhw8f9rsUosQw40we48A5u+Sr6knzs5PIkP/NN23ahMLCQpw5c8bvUohoDdu2bcPp06dxxUSDiIhyATPO2eUbIvJdLEc1ZgCUA3gzgG/6WlUS+vr68Oyzz6K2thZBl8zsRnDt82ynU75sppKurHAsh6ZN/tSll7Jls8COdbaP83z8qaBt3+fo/DSwSi/jONNQA7E9rm3mOq4850eR/Ytew86//dym4LZipk63/bhNptnmsR3bS/zzDTHTrMd5zQAgnOSXYmO2d7l/Ir2akxEIBNDb24vR0VHcfvvtnu6biCjd8YxzFlHV/wrgdwFcAVADYAHA70d+nhGKiorQ2NiI8fFxv0shojXU19fj2rVruHAh7b97TDlONDU3yl0845xlVPVHAH7kdx03o6OjA0NDQ2hubkZRUdrN2UJEWJ686MUXX8TOnTshwvZclKbYOo48xjPOlHaCwSC6uro4UxlRGisvL8emTZtw8qT9WgVRmkjVFwN5xjmn8YwzpaWGhgYcO3YM09PTqDI5YD+5ZZ5tz9zw7JxjOVBS4li2mWabO7Z9n5PNw8ZjM8hufZhtP+JgWZlzfyYnbDlywibTbDPIVkyPabucb/ZnHptb72W33snRvydt1jsm527z17bXt1stGzgRkH0/J6u7uxt79uxBQ0MD8vPz3e9ARJTheMY5B4jIL4rIe/yuIxmcqYwo/eXn52Pbtm04ePCg36UQrY5nnMljHDjnjudFpNHvIpJRUVGB8vJyTExM+F0KEa2htbUVU1NTmJubc9+YaIPxy4HkNQ6cc4Cq/o2qvqyqk37Xkqyenh4cOnQIS0tL7hsT0YYTEfT19WFkZMTvUohi8YwzeYwZZ0pr0TOV9fT0+F1ODLfMs83D6qLp8WvysYEKZ0YaAfO37cyNKclj9mUyyjF9lm0u2GxvexcDJl9teyHbrK/pV2z7Okfnml0zzaa2mF7IJjdsj2XPCNj9hW7i7KhbHjqGqT18xZlhdss8p1Mmei01NTUYHx/H2bNnUVdX53c5REQpw4FzFhGRpwAU2h8DUFV9vQ8leWLbtm3YvXs3WltbUWy/PEdEaaG/vx8//OEPUVNTg4D9g4/ILzw7TB7jwDm7fAjA5wH8KwBZk20IBALXLwXfcccdfpdDRKsoKSlBXV0djh07hvb2dr/LIWIemVKCpwWyiKo+B+DLAF6lqsejb37XdrPq6uoQCoUwNTXldylEtIbt27fj+PHjWEyD+AgRUSrwjHOWUdU/9ruGVOnv78cLL7yAN7zhDWk7U5lr5tlETWy/Yb1seiuXO3slx6Vh531tfrrImeLRa86LErb/sOv9TSZaCpyPLTQz41gOxonZxOSvbS7YZSAmppaYWpMcyNl+3NF5cvu4Y/o2G/Y1jiHxz19IofOxhC5dirt9XmODY1kv3chzf2vmkfi13KS8vDxs374d+/fvx6te9aqUHosoIZw5kDzGM86UMcrKylBdXY0TJ074XQoRraGpqQmzs7OYnZ1135go1dhVgzzGgTNllK6uLhw9ehTXrsU/y0dE/hAR9Pf3c/IiSgvs40xe48CZMkp+fj7a29s5UxlRGtu0aROKiorwyiuv+F0KEZGnmHHOASLyUwAWVXXQ71q80NLSgqGhIczNzaGsLIkMsA/cMs951dXOO+Q5+w3D9PCV6OxuyJlpjsko22xskr2LA2WljuXw3Lxz/wXOXHJMprmy0lnf1YW1D2by2eFF21PayT5voRlnLMBmmm1GOlhV5bz/9LRzOU6O2PZVtiefbI4dJhMdk9cOOs9fxOzf9Ou2zyvCzuduafI00kFfXx+effZZ1NXVIWj6aBNtGJ4dJo/xjHNueA2AMhFp9rsQL6xcCuZMZUTpq6ioCE1NTTh69KjfpVCuSlFMg1GN3MaBc5YSufFVfVX9I1V9XFVP+VmTl6qrqxEMBnH27Fm/SyGiNbS3t2NiYgJXryY52yKRV/jlQPIYB85ZRETaReSfReQUgKMickJEviEiXX7Xlgp9fX0YGxtD2FyqJqL0EAwG0d3djbGxMb9LISLyBDPO2eVhAB+OTIQCABCRHQAeBXC3b1WlSElJCerr6zNqpjLXzHNdrfMOly87FqV68/V/65VJ5zqTjQ2b+7qxvZRjMs2mtzKWnH2gY+4/7zx+cEvdjXWnJhzrbG9kNzaTHDBZd5u3tuxjs32b4/ZKdum77Na3OV6P6FX3t+A8Wxssit/X2e4/1b2b3WzZsgXj4+OYnp5GlcmWE6Uczw6Tx3jGObsURQ+aAUBVn/WrmI3Q2dnJmcqI0piIYGBggO3pyBfMOJPXOHDOLj8WkUdE5F+LyJtF5OdF5BEAL/ldWKpEz1RGROmpoqIC5eXlmJiYcN+YiCiNceCcXX4TwGMA7gTw8wDuAvB45OdZizOVEaW/3t5eHDp0CEsm4kNElEmYcc4iunwd9OuRW86IvhR81113QUTc75QmbOb5zUXviLu9Xpi+8W/TtzlQUuJcLixyLNusbMy+zf6CJjfs1hvZZpzt/pZOrX220fY+tjlhu28xfYHdMs2WrS10KfGZKN2eR9uH2T5PMZlm08PaLUMd0y87mXy2jwoKCtDa2orDhw+jp6fH73IoVzBWQR7jGWfKClVVVSguLuZMZURprK2tDWfOnMHlJL+4SrQu7ONMKcCBM2WN3t5eHDhwAKEkOzQQ0cYIBALo7e3Fo48+imeeecbvcoiIksaBM2WNlZnKjhw54ncpRLSGuro61NXV4cEHH/S7FMoFnACFPMaMM2WV9vZ2DA0NoaWlBUVFRe53SDM2D2tF55jd+irbfeXV1jjXm17GNlccmpuLW4tlc8NWdI5ZCk0vYtOX2ZIC0yPaXOq3eWxbu80dW27Pe/T+9Zp5nk3m2W1fbhnpZKVrpjmee++9F//0T/+EH/7wh9V+10JZjoNc8hjPOFNWWZmpbHR01O9SiGgNZWVlK2ec+TuIiDIKP7Qo62zZsgVXr17FxYsX/S6FiNZw++23A8A5v+ug7CXglwPJexw4U9ZZaU83MjLCmcqI0lR+fr77RkQ3ixln8hgzzpSVKioqUFFRgYmJCTQ3N/tdTsJi+jqX/aJjOTrbm1dX61i3dDb+yTu9YrK4JtPsxq1Pc1IZapdj2x7UNr9tqUsnFTF9otUc3y0jHS/v7da3+WbZ19m+jjbjbN9DRDmLZ4cpBXjGmWKIyJ0iMrjKzwMi8lkReUZEBkWk04fyEtbT08OZyogyULZ8BhFR9uHAmRxE5PcAPAxgtZYUDwAoUtW7AHwIwEMbWFrSVmYqO3TokN+lEFGCsukziNKAT1ENEckXkS+LyJCI/EBEfsbTx0W+4cCZrCMA3rbGup0AvgUAqvosgDs2qqj1amtrw9mzZzlTGVHmyKrPIPKZfxnndwKYUtU3AHgLgL/w6BGRz5hxJgdV/ZqItK2xugLATNRySETyVNWRhTh16hS6u7uvL9933324//77Pa81UUtLSxgcHERJVA/kTBHTr7iy8vq/babZ5oJtv2C3HHCwvNyxbLOzku/8uLAZ59DFGceyBINxj+fY1uSnk+11bHsr29xxaMZZW8C8F5LpWZ1X7Ww9vDQ1lfB918Mtu24NDg6mphCPPPbYY3j88cdXFmvsei8+g4hW+Jhx/gqAr66UAYDv0SzBgTMlYxZA9OgqsNovrObmZuzbt2/jqkrAc889h46ODtTUxPyeTmufwF/6XQJlmF27dvldQly7du3CQw8tJyxE5HySd0/oM4hoA9SISPQvus+p6udWFlR1DgBEpBzLA+iPbnB9lCKMalAy9gJ4KwCIyA4AL/tbTuL6+/sxOjrK9nREmS1jP4PIJ6mLapxX1Tuibp+DISJbATwF4Muq+nepe5C0kXjGmeISkc0AHlbVtwH4OoB7RORpLF96eq+vxSWhrKwM1dXVOH78ONra2vwuJ2G2tdg9gbdf/7eNN8BOwW24tZ+LiWaY/btFL2x0wy5Hxydsy7Z42y4f3Pk3fkwMxR47bqWxERgrXos5r6MZbhEbN9nefi5bPoPIBz72XBaRegBPAvgtVf2eP1VQKnDgTDFU9RiAHZF/X0DkizqqGgbw6/5VdnO6urqwZ88eNDU1cfIFojSWrZ9BlFM+AmATgI+JyMciP7tXVZNroE9phwNnyhn5+flob2/HgQMHMDAw4Hc5RESUYn59OVBVHwTwoD9Hp1RixplySktLCy5evIhLJpZARERZiFNuk8d4xplyioigv78fIyMj2LFjh9/lJC06z2qn49bFa3bzpLhNOx020dtgVZVz++nphNfH5KdNPttmoGOm4Hbh9ljcpsm2y9HTiS+dczaCCJjpvN2y5Mw0E20cTrlNXuMZZ8o5mzdvRn5+Ps6cOeN3KURERJRBOHCmnNTX14exsTGEw2G/SyEiolRhVIM8xoEz5aTi4mI0NDRgfHzc71KIiCgVUjVo5sA5pzHjTDmrs7MTu3fvRnNzMwoLC/0uJ2luvYiT5TbttBQ4c8nqcnybeXbsyyXTbNkcsM0Vi8sU2zH7czmezSFH55pjekwbrllx81jcpjonovWRyI3ISzzjTDkrGAyiq6sL+/fv97sUIiIiygAcOFNOa2xsxNzcHGZczlASEVEGYlSDPMaBM+W0lfZ0w8PDUOWnIRFRNhFNzY1yFzPOlPOqqqpQWlqK06dPo7Gx0e9yEmb7+d4TeHvc7W02N2CytUtTU871JSWOZdsn2uaUvWSPbfPcMb2SzbLNQCMUct4/yUx1XnX1jXVz8/FrMdz6NttMc16T8z34xMn/EXf/RES0cXjGmQhAb28vDh48iJAZYBERUQZjVIM8xoEzEYDCwkI0NzfjyJEjfpdCRERe4cCZPMaBM1FEe3s7JicnccXl0jsREWWAFOWbmXHObcw4E0UEAgH09PRgbGwMt912m9/l3DTJc/Zdjsn1mqyuBIPO7V36NIvJ7trj6ZIzEx2d9bU9oQN5zo8ie2y7b3t/u32g3NlLeensOSQjOtMMANhcdePY5nkM5jtrV7Pe1gqXPs7XttUlUSkREW0knnEmilJfX4+FhQVcvHjR71KIiOhmMapBHuPAmSiKiGBgYIDt6YiIsgCjGuQ1DpyJjPLyclRVVeHUqVN+l0JERDeDZ5zJY8w4E62iu7sbe/fuRUNDA/LyMuN/E9vX+c1lv+hYtpljNa337LIrDTsWbV9nMX2jHf2LTc7Xcut9bB+LlWymOVjmzESHTW/lQFHhjWPbftY2r11b41jWWee+LPs8fXf3f4pfLBER+YZnnIlWUVBQgLa2Nhw8eNDvUoiIaJ0Y1SCvceBMtIbW1lacO3cO8/Pz7hsTeWRpacnvEoiyQ6piGhw45zQOnInWEAgE0NfXh5GREb9LoRxy6NAhv0sgIqI1ZEZ4k8gntbW1GB8fx7lz51BbW+t3OUlx64UcndsFAL3m7OMcsz+TM7Z9od36OCezrT1WwOalbU9qw2akbR475v5i8tklJeb+N04xiXne7HKo2WScn5t0LAf7uxzLT7z8iRt1hcMYGhpCS0sLSktLQUQ3iWeHyWM840zkor+/H6OjowiHw+4bE90EXuUg8o6AGWfyHgfORC5KS0tRW1uL48eP+10K5YCVKxvnz5/3uRKiLMCMM3mMA2eiBHR1deHYsWNYdIkIEHmhv78fIyMjnISHiCjNMONMlIC8vDx0dHTgwIEDuOWWW/wuJyG2r/M9gbc7lkNz8Xshx2SiXXLGbr2V420bk0m225se02612Iy0m9Cl+L2WA1du7C9QYXs+zzm3HR133nnzJuexRtxbHEZf5Whra3PdnohWJ/zjkzzGM85ECdq6dSump6dxyWWQReSFrq4ujI+P8yoH0XqxHR2lAAfORAkSEfT392N4eNjvUigHRF/lIKL14ZcDyWscOBMlYfPmzSgoKMArr7zidymUA3iVg4govTDjTJSkvr4+PPfcc6irq0MgkDl/ewbLyx3LuujMGdscsc0hu51kCVZWOpZDMzNrbhsoLo6/M5daxGSibUY6pg+0y3rLPldSEbWc78x+y9WFuPsKXbjoWLbZ83hWrnKMjIxgx44dCd+PiCJ4dpg8ljm/9YnSRHFxMRoaGnD06FG/S6EcsHnzZuTn5+PMmTN+l0KUcRjVIK9x4Ey0Dp2dnTh58iQWFuKfbSTyQl9fH8bGxjgJDxGRzzhwJlqHYDCI7u5ujI2N+V0K5QBe5SBaJ3bVII8x40y0Tg0NDRgfH8f09DSqqqr8LseV7VVsc8aCoHN9UaFzfaFzOWz2p9cS780cKCt17mv+snPZrQVb0PzNH4p/JjbZvs4x+e/zF9Y8dviys3YrmUxzPJ2dndi9eze2bt2KQvNaENEqGKugFOAZZ6J1EhEMDAxwhjfaELzKQbQOPONMHuPAmegmVFZWoqysDJOTk36XQjmgoaEB8/PzmJ6e9rsUIqKcxIEz0U3q6enBwYMHETIt1Ii8xqscRIkTsKsGeY8ZZ6KbVFhYiJaWFhw+fBjd3d1+l7Mmm7V9c/G7HMsxOWDbMGRuzrGYV13tWF6amnIsBwoKnPfPu/FxE7ro7PEc3FzlrMUlN+y2Plm21mQy0bbns82Sey36KkdTU1NKj0WU8fgHJnmMZ5yJPLBt2zacPn0aV65c8bsUygG8ykGUGJ5xJq9x4EzkgUAggN7eXoyOjvpdCuWA6KscRES0cThwJvJIfX09rl27hgsXLrhvTHSTeJWDyEWqOmrwjHNOY8aZyEMDAwP40Y9+hJ07d0JE/C4nrmR7G0tevvP+Llle24s5EJVxlnznR0942pl5zqutca63fZ5dMs5B01c7ZLpQxGSaba22x3V0Pts8brvsVd9mN9FXOW6//fYNOSZRphFOtkke4xlnIg+VlZVh06ZNOHnypN+lUA7gVQ4ioo3FgTORx7q7u3HkyBFcMzPpEaVCf38/29MRrYVRDfIYB85EHsvPz8e2bdtw6NAhv0uhHFBeXs6rHERrYFcN8hozzkQp0NraiqGhIczPz6O0tNTvctYlr67WsaxXnJloXXSeUQ8UFjl3oM5wYXQu2WaMpcz5HC2dO59UrfbY6tYH2mSaY9Yn8YW7jco0x9Pd3Y09e/agoaEB+fn57ncgygUK9nEmz/GMM1EKiAj6+vowMjLidymUA3iVg4hoY3DgTJQiNTU1EBGcO3fO71IoB7S2tuL8+fOYn5/3uxSitMGoBnmNUQ2iFOrv78cPf/hDvOENb0AgkFl/py6dTW7Ab6eejhEVj4iJSlww7eBM9MKtdZ5dHzPdt4fSIZqxmpWrHMPDw7jzzjv9LocoPXCQSx7LrN/kRBmmpKQEdXV1OHbsmN+lUA6oqalBIBDgVQ4iohThwJkoxbZv347jx49j0eULaURe6O/vx+joKMJhzvxAuU3AqAZ5jwNnohTLy8tDZ2cn9u/f73cplANWrnIcP37c71KI/KWauhvlLGaciTZAc3Mzjh8/jtnZWVRUVPhdDoDYrO49gbc7lmOmvZ6ZdS6bM+i2pZxtXxdPsKzMsRyam4u/vclT22mv3drNJStdc81r2b59O4aGhtDU1ISCFOa9idIdzw6T13jGmWgDiAhneKMNs3KV48CBA36XQkSUVThwJtogmzZtQmFhIc6cOeN3KZQDmpubMTMzg9nZWfeNibIVp9wmj3HgTLSB+vr6sH//foRCIb9LoSzHqxxE/HIgeY8ZZ6INVFRUhMbGRhw9ehTbt2/3u5y49NpS3PU2l7x0+pWUHcuymWbbt/lmM86Zlmley8pVjldeeQUNDQ1+l0O0sRRAmKNc8hbPOBNtsI6ODkxMTODq1cS/PEe0Xn19fThw4ACvchAReYADZ4ohIneKyOAa614QkcHI7dENLi0rBINBdHd3Y2xszO9SKAesXOUYHx/3u5SE8TOIPMOMM3mMUQ1yEJHfA/AuAPOrrCsCIKq6a6PryjZbtmzB+Pg4pqenUVVV5Xc5lOU6OjowNDSE5uZmFBUVud/BR/wMIi8xj0xe48CZrCMA3gbgy6usuxVAiYg8ieX3zkdU9dmNLC5biAgGBgbw0ksv4e6774aI+F1SjND0tGPZ5ohtr+VAcbHZgTMaIFH3t/cNLzhjK5KX79x3UWHcY9tMs+3zLOb+S+fOI5dEX+V4zWte43c5bvgZRERpiwNnclDVr4lI2xqrLwP4NICHAWwH8ISIdKuq45tdp06dQnd39/Xl++67D/fff3+KKs5sV65cwXe/+13k5+e7b0wbZnBw0O8SUmJ+fh4XLlxAMBj0rYbHHnsMjz/++MpijV3vxWcQ0XXsKEMe48CZknEQwGFd7m11UESmADQAOBm9UXNzM/bt2+dHfRlncXERe/fuxc6dO5GX5+//jp/AX/p6/HSya9cuv0tIidnZWd+vcuzatQsPPfQQAEBEkj31n9BnENEKRjXIa/xyICXjfQAeAgARaQRQAeC0rxVluIKCArS2tuLw4cN+l0I5oKKiAuXl5ZiYmPC7lPXiZxAlLlVfDExiMB7vi66UmXjGmeISkc0AHlbVtwH4AoAvisgeLH90vG+1S6Rse5WctrY2DA0NoaWlBSUlJX6Xsya33sjhK1ccy8HKSsdyaGYm4WPp0jXn8mJyZ0dtn2fYZSNb+jYnore3F3v37sWWLVt8v8qRiPV8BhGlg3hfdKXMlf6fmrThVPUYgB2Rf1/A8hd1oKqLAH7B7f5XrlyBqqblF97SUSAQQG9vL0ZHR3HHHXf4XQ5lueirHD09PX6Xs6qb/QwiAgABIP5mnON90ZUyFKMa5LmZmRkcOnTI7zIySl1dHUKhEKampvwuhXJAW1sbzpw5g8uXL/tdClFqhVN0A2pEZF/U7VftoVX1awCu2Z9TZuPAmTxXUFCA//yf/zMjG0nq7+/HyMgIlN8CpxSLvspBlM1ENSU3AOdV9Y6o2+f8fqy0MRjVIM/V1tbiwIED+P73v4977rnH73IyRllZGaqrq3HixAm0trZu+PFtzvfNpe92bmD/EBLzd7eGnZsnkWl2I4XOPsxBsxyTaTZyKcOcqLq6OoyPj2NqagrV1dV+l0NElBF4xplS4jOf+Qx+/OMf4+rVq+4b03VdXV04evQorl3j1T1KPV7loKyWBl01KPtw4EwpsXPnTrzjHe/gpeAk5efno729HQcPHvS7FMoB0Vc5iLKPLk+AkopbohWoHlPVHSl8kLTBOHCmlNmyZQuuXr2Kixcv+l1KRmlpacHU1BTmzLTSRKnAqxxERIljxplSRkQwMDDg+0xlmUZErl9Cv/POO32rw/ZlDhQUONcvOGM4geJiz44dcOlnzUyzd6KvcvT39/tdDpGnOHMgeY1nnCmlKioqUFFRkckzlfmiuroawWAQZ8+e9buUjMSrHMnhVQ7KWj5HNSj7cOBMKdfT04NDhw5haYkTfCWjv78fY2NjCIfD7huTA7/wlpzoqxxEWUMBCafmRrmLA2dKueiZyihxxcXFqK+vx7Fjx/wuJePwKkfyeJWDiMgdM860Idra2jA0NISWlhaUuORX6YbOzk4MDQ2hubkZBSZjnG5sJjqeYFWVYzk0Pe3cl8uMdnmm7/AT5z7rWF5cXMTevXuxZcsW5OXxYy5RfX19+MEPfoCamhoEAjyvQlmAV57IY/xkpA2xMlMZLwUnJy8vD9u3b8fY2JjfpWSUlascnPo9OSUlJdiyZQvGx8f9LoXIG+zjTB7jwJk2TF1dHcLhMKampvwuJaM0NTXh0qVLmJ2d9buUjNLW1oazZ8/issvZa3Lq7OzEiRMnsLi46HcpRDcthVNuU47iwJk2FGcqS95KW7/h4WE+b0lYucrBSXiSs3KVY//+/X6XQkSUdhj+ow21MlPZ8ePH0dbW5nc5GaOqqgrFxcU4ffo0GhsbN+SYXvZlBoBAYdGNhVDIsS5YXu5YVtOBxeanbaZ5LXV1dRgfH8fU1BSqTS6a1tbU1IRjx45hdnYWFRUVfpdDtH482UAe4xln2nBdXV0YHx/nTGVJ6u3txcGDBxEyg06Kb2BggFc5ksSrHJQVFEA4RTfKWRw404ZbmanswIEDfpeSUYqKitDU1IQjR474XUpGKS0tRXV1NU6cOOF3KRll5SrHK6+84ncpRERpgwNn8kVLSwsuXrzImcqS1N7ejsnJSVy9etV9Y7quq6sLR48e5VWOJPX29uLAgQO8ykEZSZCaLwbyy4G5jRln8sXKTGXDw8PYsWOH3+VkjGAwiJ6eHoyOjuK2225L6bFsrjhYWencwKwPmD7TYduVIXjj7/TQpUtJ1fKd8FeS2t6KvsoxMDBwU/vKJdFXObq6uvwuhyh5HOSSx3jGmXyzefNm5Ofn48yZM36XklHq6+uxsLCAixcv+l1KRuFVjvXhVQ7KaKqpuVHO4sCZfNXX14exsTGEw/y2RaJWztbzC2/JERH09fVheHjY71IySvRVDiKiXMeBM/mquLgYDQ0NnKksSRUVFaioqMCpU6f8LiWjVFdXIy8vj1c5ksSrHJSR2FWDUoAZZ/JdR0cHhoaG0NzcjMLCQr/LyRg9PT3Yu3cvGhoakJeX+v+VQzMzjmXXTLMRTmIGv5vNNMfT39+P5557DrW1tQgEeO4gEStXOV566SXcfffdEBG/SyJKCL/IR17jbw3yXV5eHrq6ujhTWZIKCgrQ1taGQ4cO+V1KRuFVjvVZucoxMTHhdylEiWPGmTzGgTOlhcbGRszNzWHGnNWk+FpbW3H27FlcTuJsLi1f5Thx4gQWFhb8LiWj9PT04NChQ1gyMzsSEeUKDpwpLUS3p+MX3hIXCATQ19eHkZERv0vJKLzKsT4FBQVobW3lVQ7KECk628zfUTmNGWdKG1VVVSgtLcXp06fR2NjodzkZo7a2FkePHsX58+dRU1Pj2X6DZWWO5ZBp4+aWabbyqquv/3tpasqxLpWZ5rU0Njbi2LFjmJmZQaXtUU1ramtrw9DQEFpbW1FSUuJ3OURrU3CQS57jGWdKKz09PZypbB0GBgYwOjrKs/VJ4FWO9eFVDiLKZRw4U1opKirC1q1bceTIEb9LySilpaWorq7G8ePH/S4lo1RVVaGkpASnT5/2u5SMUltbi3A4jPPnz/tdClF8bEdHHmNUg9JOe3s7du/eja1bt6K4uNjvcjJGV1cX9uzZg6amJuTn59/0/mRLnfMHh51RjWB5uXO9OWtrox1PnPvsTdeUCr29vXjmmWdQX1+PYDDodzkZY2BgAPv27cMb3/hGtqejtMV2dOQ1nnGmtBMIBNDT04OxsTG/S8ko+fn56OjowIEDB/wuJaPwKsf6lJaWoqamhlc5KL3xy4HkMQ6cKS2tzFR24cIFv0vJKFu3bsXFixdx6dIlv0vJKO3t7ZicnMSVK1f8LiWjdHV1YXx8HNeuXfO7FCKiDcGBM6UlEcHAwABGRkb4xa0krHzhjV/cSg6vcqwPr3JQWlMAYU3NjXIWM86UtsrLy1FVVYVTp05h69atfpeTMTZv3oz8/HycOXMG9fX1696PnjkXd33I5ay2Hy3mbkZ9fT3Gx8dx8eJFbNq0ye9yMsbWrVuxZ88eXLp0CeU2907kK8YqyHs840xprbu7G4cPH+ZMZUnq6+vD2NgYwmF+/TtRbE+3PrzKQUS5hANnSmsFBQVoa2vjTGVJKi4uRkNDA8bHx/0uJaNUVFRcv8pBiYu+ykGUVvjlQPIYB86U9lpbW3H27FnMz8/7XUpG6ezsxIkTJ7CwsOB3KRmFVznWh1c5KC1x4EweY8aZ0t7KTGWjo6N47Wtf63c5GSMYDKKrqwv79+/HrbfemvT93TLMVqZlmtcSfZWjt7fX73IyRvRVjo6ODr/LIbrx5UAiD/GMM2WE2tpaqCpnKktSY2Mj5ubmMDMz43cpGWXlKsfly5f9LiWj8CoHEWU7DpwpY6x8AYlf3ErcSls/fuEtOStXOfiFt+REX+Ug8p8CGk7NjXIWB86UMUpLS1FbW8uZypJUWVmJ0tJSnD592u9SMgqvcqwPr3JQWmHGmTzGjDNllK6uLgwNDaGpqQn5+fl+l5Mxent78cwzz6C+vh7BYHBd+8iWDHMy+vv7sW/fPrzxjW+EiPhdTkaIbuv3+te/ns8b+YcZZ0oBnnGmjJKXl4eOjg48+uijeOaZZ/wuJ2MUFhaiubkZX/rSl/i8JaG0tBSHDh3C17/+db9LyShVVVV4+eWX8dhjj/ldChGRpzhwpoyzdetWlJWV4f3vfz9bXyWhvb0d+fn5+MAHPsDnLQmvf/3r8bGPfQxXr171u5SM8oY3vAEf/vCHce3aNb9LoVzGqAZ5jANnyjgigre85S3o7OzEF7/4Rb/LyRiBQAD33HMPurq6+Lwloa+vDzt37sSnPvUpv0vJKP39/bjrrrvwx3/8x36XQrmMA2fymPCb9uQ1ETkHYCO+wRcEsBnAuQ04Vjbh85Y8Pmfr4/a8tapq7QbWQzmksqBeX1/3b1Ky729NfOZ5Vb0jJTv30PDwsKeDvIGBgQ390kLbh77hWf3HPvnTntTOLweS5/iLkIiI/Mezw+Q9DpyJiIgo+ygAfp+DPMaBMxEREWUnnnEmj/HLgURERERECeAZZyIiIspOPONMHuMZZ8oYIhIQkc+KyDMiMiginWts90Jk/aCIPLrRdaYLt+cr0ecz1yTyvPA9tjYRuVNEBlf5Od9vtMF0eebAVNwoZ/GMM2WSBwAUqepdIrIDwEMAfjZ6AxEpwnKbxV0bX17aeQDxny+39bnqAcR5XvgeW5uI/B6AdwGYX2X1A+D7jTaSAqr8ciB5i2ecKZPsBPAtAFDVZwGs1kPzVgAlIvKkiHw/8gs6V7k9X4k8n7nI7Xnhe2xtRwC8bY11fL8RUcbjGWdKSyLySwA+YH58BsB3opZDIpKnqktRP7sM4NMAHgawHcATItJttskVFQBmopbt8+W2Ple5PS98j61BVb8mIm1rrOb7jTYeYxXkMQ6cKS2p6hcAfCH6ZyLyJwDKo34UWOWX7kEAh3V5SsyDIjIFoAHAyVTWm6ZmEf/5clufq9yeF77H1ofvN9p4/HIgeYxRDcokewG8FQAil8dfXmWb92E5OwkRacTyWa7TG1VgmnF7vhJ5PnOR2/PC99j68P1GRBmPZ5wpk3wdwD0i8jQAAfBeABCRzQAeVtW3Yfks9RdFZA+W5416Xw6f1Yp5vsxzterzSa7PG99jCeL7jXylypkDyXOivIxBREREWaYyWKN3ld6fkn1/+9IXn1fVtP+C6/DwsKeDvIGBAfFyf27aPvQNz+o/9smf9qR2nnEmIiKirKQ840weY8aZiIiIiCgBPONMREREWUjZVYM8x4EzERERZR8F+ziT5zhwJiIiouzEKbfJY8w4ExERERElgANnohwnIneKyGAC231cRH4gIk+LyOtuZlsRWRSRJ0WkTURmRWQwcntWRL4rIpsi272yyn3/QER+3fzs2ci+HheRqyJS5PrAiSirKQANa0pulLsY1SDKYSLyewDeBWDeZbvbAPwEgDsBbAXwNQCvvYltL6jqT4lIG4BRVd0Vdf8/BPBLAD6d7ONR1ftE5Fiy9yOiLKTKqAZ5jmeciXLbEQBvS2C7nQCe1GUnAOSJSK0H2zqIiGB5sH0xke2JiNKRiARE5LMi8kzkalqn3zWRNzhwJsphqvo1ANcS2LQCwEzU8iUAlR5sCwB9kV8sLwE4COAwgL9JoCaL10+JyMHHqMYDAIpU9S4AHwLwUCofJ20cRjWIKBGzAMqjlssBTHuwLRCJaohIMYDHAJxR1aU4218BUGh+Vhb5ORHRDf5FNXYC+BYAqOqzIpL203NTYjhwJqJE7AXwKRH5NIBmAAFVPe/Bttep6hUReQeAF0XkaVX98RqbvgDgQyLyl6q6JCIdAApV9WzSj4qIstYlXPz2d/WrNSnafZGI7Ita/pyqfi5q2V55C4lInstJAcoAHDgTkYOIvAcAVPWLKz9T1edFZAjAM1iOeP37ZLdNhKqeEZH/COCvReT1AKrNL6eHVPV/i8hOAM+LyCwAAfDu5B8pEWUzVX2Lj4e3V94CHDRnBw6ciXKcqh4DsCPqRy8AiLmsqKp/AOAPzI+T2TaRY0NV/xeA/xVZzF/jfh8H8HG3/RMR+WQvgPsB/KOI7ADwss/1kEf45UAisi4AeDQF20bbLCJPruN+cYnI4wC2eL1fIqIkfR3AVRF5GsCfAviAz/WQR3jGmYgcVPVUKrY19ytYz/0S2O99qdgvEVEyVDUM4NddN6SMwzPOREREREQJ4MCZiIiIiCgBHDgTERERESWAA2ciIiIiogRw4ExERERElAAOnImIiIiIEsCBMxERERFRAjhwJiIiIiJKAAfOREREREQJ4MCZiIiIiCgBHDgTERERESWAA2ciIiIiogRw4ExERERElAAOnImIiIiIEsCBMxERERFRAjhwJiIiIiJKAAfOREREREQJ4MCZiIiIiCgBHDgTERERESWAA2ciIiIiogRw4ExERERElAAOnImIiIiIEsCBMxERERFRAjhwJiIiIiJKAAfOREREREQJ4MCZiIiIiCgBHDgTERERESWAA2ciIiIiogRw4ExERERElAAOnImIiIiIEsCBMxERERFRAjhwJiIiIiJKAAfOREREREQJyPO7ACLKbYODgx2BQOCDAN4ZDofLA4HAJQB/Gw6H/2TXrl1H/K6PiIhoRdyBs4jsAvCPAEajfnxOVd9utlsEMKiqP7WeIkTk/wL4L6r6/aif/Q8ALwP4ZQBQ1R1J7FDbfv/x9ZRCRBvolpogfuf2ErQ0NKKluQmFhYVYWFiomJyc/M3JycnffOGFF1BQUOB3mUSUgQYGBsTvGij7JHLG+fuq+m9dtrmw3kFzxOcBvBvA9wFARAoA3A/gIwC+C+Dvb2LfRJSGaosFv3N7CV53262orKy8/vPi4mJ0dHSgpqYGL730EoLBIILBoI+VEhERLfM04ywiu0Tk2yLyf0TkRyLyayLyDyKyX0R+I7LNT4jIHhH5vyLyiIjkA/gqgDeJSElkVz8L4ElVnfeyPiJKH29tL0BLc6Nj0BytsrISDQ0NuHr16gZXRkREtLpEBs5vEpHBqNvvumzfDODnAPwGgI8CeBeAewH8mogIls8uv01VfwLABID3qOpVAP8M4F9F9vFeAH+d9KMhoozx+qZ8tDQ3xd2mqakJi4uLG1QREWUzEckXkS+LyJCI/EBEfmaVbT4gIiNRY57uqHU9IjIjIkUe1fM3IvK+VY7/37zYP6WGV1GNaMOqek1EpgEcUdVFEbkIoAhALYAGAP+4PIZGMYDvRO73eQB/LCKDADap6o+SOCYRZZgCURQWFsbdprCwEKFQaIMqIqIs904AU6r6LhHZDOBFAP/HbHM7gHer6vPRPxSRCgAPAVjwsJ7PA/hvAB6J+tkvAnjAw2OQx1LRVUPjrDsP4BSAn1XVmchfe3MAoKovi0g5gN+B801ERFloUQULCwsoLi5ec5uFhQXmm4nIK1/BcjQUAATA0irb3A7gwyKyBcA3VPUPI1fLP4fl7139i72DiLQB+AcAJwG0Yfl7WQMAXhPZx0dE5BYAfx457hSA96nqHhGpFZFWVT0uIq8F8IqqHvPqAZP3Ehk4vylyFjjavap6JdmDqWpYRB4E8A0RCQCYxfKXAlc8AuCPAbQku28iyixPT1xD26kJdG/vXHObiYkJdtUgIk+o6hwARE7SfRXLcVLr7wH8JZbHJ18XkfsA3IHlAfCPI1fLV9MO4KewfCV9HEATgMsAjmN5wP15LA+WR0XklwD8HoD/BOALWD4T/gkwppoR4g6cVXUQQF2iO4tsPxj5934AuyL/ngbQE/n3kwCeXOP+X8Dym4iIstw3jy5iZ/MkttTVrvoFwZmZGZw+fRoVFRU+VEdE2UhEtgL4OoC/UtW/M+sEwJ+p6kxk+RtYPmv8TgCnIgPeLVgew7zR7Ppo5Er6AoAzqnohso+Vq/C9AP4qMvDOB3Ao8vMvAfieiDyE5THT73j4cOMKhUK4evUqFhcXEQqFEAwGUVBQgKKiIl7pi8OrqMZmEXnyJlvSxRCRTgB/m/QdVeWYl4UQUUoMDg7e++Mf//irTU1N+Y2NjfmRPs6YnJy8NjExcS0cDv/8rbfe+oTfdRJR5hOReiwPen9LVb+3yiYVAIZFpBfAPIA3AXhEVf+/qH0cw/KZZSteTBUADmA5O31CRO7G8ve9oKrnRWQMwMcAfF1VV4uPeG5xcRHz8/NobGxEY2Mjoj57MTk5idLSUl7tW4MnA2dVTcmzq6qHASQ+8QkRZZRdu3Y9MTg4+KqJiYkPTExMvCscDpcFAoE5AF8Oh8N/ypkDichDHwGwCcDHRORjkZ/di+WOXmWq+jkR+QiAp7D8JcDvqeo3PTr2bwD4kojkYXmQ/UtR6z4P4JsAule7o9dCoRDm5+fxqle9ij3010FU3f5IIiIiIqJMMzw8HDPIm5+fR01NDTo6Ota83+HDhzE1NYXS0lLHzzd6Nsa2D33Ds0HqsU/+tCe1ezoBChERERGlr8XFRTQ2Nsbdhj3015aKdnRERERE5LPVzhAPDg6GCwsL4559jfTQDw8MDPia1fDqLLGXeMaZiIiIKEcEAoFLCwvx53FZWFhA5PsmZHDgTERERJQ7/nZycvJavA0i67+8QfVkFA6ciYiIiHJEOBz+k4mJiWszMzOrrp+ZmUGkHeifbnBpGYFdNYiIiIhyyODg4L2BQCBuD/1du3axh/4qOHAmIiIiyjGDg4MdgUDgAwDYQz8JHDgTERERESWAGWciIiKiDCIiPy0iL4nIARH5iohUrGebJI9ZKyJPiMioiAyLyOtv5lgi8pCInBCRFyO3f1hlmztFZJ+IjInI90SkYR11l4jI30X2cUBEHoj8/Cejjv1iZJ2KyO1x98czzkRERESZQURqAYwAuFtVD4nIHwEoV9XfTGabdRz3HwG8qKr/XUReDeAbALYDKF3PsUTkGQD/QVWfXmN9AYAjAP6tqu4Vkd8AcL+qvjXJuj8FoEpVf1VEWgA8C+B1qnrKbPdVAIdU9cNx98eBMxEREVFmEJF3APgFVf3pyHIbgB9jeXCobtsAyAfwRwB+AkAQwI8A/I6qzkYd44sABlX1i5HlPACzAFpV9VzkZ4MA/hxAsVs9qzyGQgAzAL4JoAPAYQAfUNUTUdvcDeBzqtofWS4AcAlAo6pOich/AvBzWE5PHAPwm6o6ucqxDkXq+2HUY3tJVf8kapt3Ang/gB2qurRazSsY1SAiIiLKHFsBnIxaPgWgAkB5gtt8CMASgNtV9VYAkwA+6XLMGgCBlUFz1D6bE6zHagTwfQAfBvBqLJ8F/hcRiZ4p0LFfVV0EcA5Ak4i8G8AtWD5z/GosD8AfXuNYq9XXvLIQGZD/IYD3uw2aAU65TURERJRJ1jrpGUpwm/uwfOb5nsg4tQDAWQAQkecAFAJoAfAmEXk/gL0APhFnf4nU46Cq4wCuRy5E5NMAPgagDcB4go/hdQD2RR5DEEDJGtuvtp/o2n4ewBFV3bNWvdE4cCYiIiLKHCcA3Bm13ATgoqrOJ7KNiAQBPKiqTwCAiJQBKAIAVb0z8rMvIjaqARHZpKoXo/Z5CsuRC7d6HETkVQBuVdXo2QkFQPSMhicANETdJx/LZ74nsDxQ/iNV/Z+RdYUANolII5bPPq94a9R+Xomq78Wobf4NgEfXqtViVIOIiIgoczwJYIeIbI8s/zqAf0lim28D+C0RKRCRAIDPYzmqsKZIhOEbAH4NuD7w7QMwmGA9VhjAn4vItsjyb2A5dxz9hb3nAFRHde94H4BnVHU68hh+Oap7x38F8GVVnVTVV0fdJiO1/Gqk7mYAbwHweGRZALwRwPdc6r2OXw4kIiIiyiAi8lYsD3ZXOk+8G0A7gIcjmd9Vt1HVCyJSDODTAHZh+cztiwB+NfrLgWscsx7LOeJtABTLHTGejHcsl/29E8t56yCWz1z/Epaz198E8FZVnRSR1wH4Cyx37piK7PdYZMD/cQBvj9RyAsAvq+rEKscpA/A/AdwWOdZ/U9W/jayrBTCpqvnxanXsjwNnIiIiIiJ3jGoQERERESWAA2ciIiIiogRw4ExERERElAAOnImIiIiIEsCBMxERERFRAjhwJiIiIiJKAAfOREREREQJ4MCZiIiIiCgB/z/gF0xlaD07rgAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Set plane programatically to a specific plane\n", "Viewer.setPlane(20)\n", "Viewer.set_clim(0,5e-7)\n", "Viewer.ax.get_figure()" ] }, { "cell_type": "markdown", "id": "c537cc69", "metadata": {}, "source": [ "The above figure shows a constant energy cut at an energy transfer of 2.54 meV with the two axis (1,0,0) and (0,1,0) in the plane. Figuratively, one is looking at the cube of data from the top. \n", "\n", "If instead one would like to look at the data as a function of both scattering vector and energy transfer, this is done by pressing 0 or 1. The number signifies the axis (either projection vector 1 or orthogonal to it, respectively) that one steps through with the arrow keys. That is, in this specific example pressing 0 flips the data around so that the scattering vector (-1,-2,0) is along the x-axis and energy along the y-axis.\n", "\n", "Moving the bottom slider to 1.06 then generates the following figure where the dispersion can be seen again." ] }, { "cell_type": "code", "execution_count": 5, "id": "471362a0", "metadata": { "execution": { "iopub.execute_input": "2023-07-25T11:38:27.171000Z", "iopub.status.busy": "2023-07-25T11:38:27.171000Z", "iopub.status.idle": "2023-07-25T11:38:27.428258Z", "shell.execute_reply": "2023-07-25T11:38:27.427389Z" } }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "C:\\Anaconda\\envs\\python39\\lib\\site-packages\\MJOLNIR\\Data\\Viewer3D.py:468: MatplotlibDeprecationWarning: The key_press_event function was deprecated in Matplotlib 3.6 and will be removed two minor releases later. Use callbacks.process('key_press_event', KeyEvent(...)) instead.\n", " self.figure.canvas.key_press_event(str(value))\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtoAAAG8CAYAAADgqwLpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAsTAAALEwEAmpwYAABZFUlEQVR4nO3deZxcZ3Xn/++p6lZr3yzL1mJbWDjCYCMSDJhlSJPgAYwViIMzGQcnxuxMBhNmfoE4MEnmhQOZIJhJwjKExRAgITaQYIGNPUDH2MaAbLyIzUhepZZtLdYutVRV5/dHlUy7dU6pq3Vvl6T+vF+vfrnr3OV5qure0uPbT32vubsAAAAAFKvS7Q4AAAAAxyMG2gAAAEAJGGgDAAAAJWCgDQAAAJSAgTYAAABQAgbaAAAAQAkYaAMAAIxgZs8zs4ExbHepmQ20fm4zs31mNrv4HuJYYORoAwAA/JKZ/YmkSyTtdvdzj2A/H5F0l7t/orDO4ZjCFW0AAIAnWyfpwoMPzOxsM/tO6yr1l81s1uF2YGbnSHoGg+yJjYE2AADAMO7+ZUkHhpX+QdJ/cfd+Sd+Q9Cdm9nIzWzPi51XDtrlC0l+OX69xNOrpdgcAAACOcmdK+qiZSVKvpF+4+/WSro9Wbs3JXubu3xm3HuKoxEAbAACgvZ9L+gN3f8jMXihpwWHWf7Gkb5XfLRztGGgDAAC091ZJnzOzHkku6fWHWX+ZpPtK7xWOeqSOAAAAHAXM7A5JO1oP73f313WzPzhyXNEGAADoMjObrOYF0P5u9wXFIXUEAACg+5ZLmmpmN5jZt81szPndOHowdQQAAEx4L3vJNN+ytV7a/m+/e+jHkvYNK31ieMa2mZ0t6VxJn5R0hqTr1EwuqZXWKZSOqSMAAGDC27y1ru9/c3Fp++9dsG6fu5/TZpV7Ja315hXQe81si5rpJg+X1imUjoE2AACAXHVvdLMDl0k6W9LbzGyhpJmSNnazQzhyDLQBAAC671OSrjKzm9WMELyMaSPHPgbaAABgwnNJDXXve2vuvl/SxV3rAErBQBsAAEBSQ12dOoLjEANtAAAw4blcdZLYUDBytAEAAIAScEUbAABA3Z2jjeMTA20AADDhuaQ6A20U7JgdaM+bN8+XLFnS7W48ye5X/H5YP+0rm8O6T+kN65UtO+IGetq8XXv3xfWh/XG9mswaqlbD8u5nnxrWp/3ksbxPif2L54T1SQ9sijfojV8nf3xbWN93ztKw3rMnT0nqHYz3pZ274voJ8XOoz5oS1it7D4R127Yz7ZPq8ZdyfM6MeP1kbqHVk384askd0Br5l4EGL1gY1k+6LX4e1QcfjXc0uS8s7z9tXljfOz9ef8pjQ/H+JU16LHnvMgfi48OnTQ7r6XuXvG+SpKnx8bHrzBPC+vQ7NsT7mRbvRzuS5zwlfg6SpO3J85gxPa7v2ROWd77w9LC+6L+vDevPnLk+7dJXP/wbYb2avN2bL9ib7ivyK/89+dxKjv3Ggvj9kaTajPjY3L7zyx31aaK6/fbbN7v7id3uB1CWY3agvWTJEq1evbrb3XiSM9/74W53AQBwFDja/n06WpnZg93uw3BMHUHRjtmBNgAAQFFcInUEhWOgDQAAIJGijcIR7wcAAACUgCvaAABgwnM5qSMoHANtAAAAl7KQJmCsGGgXaPGVt4b1JERNlUmT4vUbSUxbb/522emnxAse2hiWfX8SN5fEhk27ezDeTxIn5lncoKRJSdtZtJZv3Ravvi/O+pryk+Q579qd9imLNaxt2RKvXo/fVf/Z42HdTkven+Q5SErjHG1rEv+YHR+74jg2TZ8a17PYP0kL/+nesJ693568rjoQHwOTfhbvp3dNEr3X5j317PzaFR+zPacuDuu2KYmOMwvLjW3b0z5Vkm2m37Iu3lcSpWd7kri+JKJQc2emfbJaEhW4O35tfd7csD793vjY/+H9cTTooycl8YGSpm6Kj8FdC+Jj/NdOiaMC122LY/kOPGV+WK/e/vOwXkliGSVp0oYklhTHHBdztFE85mgDAAAAJeCKNgAAgEx1xX91AsaKgTYAAJjwXFIycxMYM6aOAAAAACXgijYAAIDE1BEUjoE2AACY8FwMtFE8BtrjoPrMM+MFD8aReZUk1q2x+KS0jYfOnx3WT/nfD8dtzI3XzyLLsqg5eedRhJnGljgerHLyiWG9mkSZ+fY4/s5m5RFnqsTPu7InjvWyqXE0XhJmJ2Vxh/U2YVJT43g67U3i5pLIvOy5KYtZrCURcZIau5O4uUo8C81mJBFuSZRjFlFoWRzl5CTmrs2+eqbF71198NGwXs3OleScqCT7l6T6I4/F20yflm4T7yh+/Wz+vHj9LBJSeUSiLYw/b2xvHElZnx+fX76lL6w/cCA+ryXpjM1xG48vjV/zB7bPCet7vh/H+524J45grCxeGNYbG+LIUEmqLMg/l3HsaTgDbRSLOdoAAABACbiiDQAAJjymjqAMDLQBAMCE5zLV+UM/CsZAGwAAQMzRRvH4XzcAAACgBFzRBgAAEx5ztFEGBtrjoH73T8N6z4KTw3ojiaerz0ji3iT1xKlrqj03jhasTY2D6Kau3RrWPYkly+LbfG6bKL09cXSXTY5jwNTbG9dr9bieRA5msX+SZIvi96L+7GVhvfLz9fGOFsSRZY2++DnUpybPTVLvI3F/fV/8+qmevB6nxJFle0+fHdan3BQfr5JUmRvHqKXxhcl7seP5p4X1mXfFx5kn0Y82NY5fbKexKH6PbPH8eIPt8cnlU+LzsZHFMkqq7I7j5uzRLXE9iRbcdVbc16n/b03c8LIlaZ/2n7A4rHvy986hOfE/G9MfjGMnn778wbB+5qw4TlGSbl383LDeSD4iFs6Iz5VfJBGCW54Vfz6deMO2sL79Vc+MG5Y059/j5/es//rhsH7n3/1xui90m6meHfjAGHFEAQAAACXgijYAAJjwXFKD648oGEcUAACAmnO0y/o53pjZOWZ2Tbf7cbTjijYAAJjw3Jmj3Ql3Xy3pNd3ux9GOIwoAAABPMLN+M/uemX3ZzH5mZneY2Qozu9HMHjKzD7fWWdNa/yoz+1sz+46ZrTWzVWYWpyV0v98VM/s/ZvZ9M/uJmf3UzF7Yqn/LzP5Xa18vNbP1ZnbSkfSptCvaZlaV9A+Slqk59ekt7r5mxDp3SDr4dfH7Jb1e0kclLZc0JOkN7r62rD4CAAAc1DgOp3gcgedIeo67/8jMrpP0p5L6Jc2UNCjphyPWf7ak35DUkPR9SRdJ+sy49faXDtfvayQtlPR8d2+Y2bslvdvdV5jZayX9yMxukfT3ki529zwiaRTKnDqyQpLc/YVm1i/pSkmvOrjQzCZLMnfvH1a7UNJkd3++mZ0raeXwbY4W51Uu6mj9ytSpYb228ZGwbtU4eq9nzQNpG4t+kUSKJdF41SSuTyfG8WMZ3xNHevm27ek29T1xXJr1xH21JM6uMmtGWG/sjftUmRGvL0mNhzaE9Z4N8XvhleSPQevj89Es/vCuZrF4knzK5Li+f3/cRk9yOm+M3+vJDzwcr5+0K0mNrXHMXmVacoxvGAzr0x5OjrMsgjE5Jxo7d8XrS3ncYXJsVmYlkZR9caacPbI5rreJkawn711qU9zGtCTu0Gu1uE9rk/daUl92KSN5zaf0JsfZtGlh+cc/PTWsP7QgiYqUtPiO+DyqTYljOCsWx0juOTU+Bk77cvy6aij+rJm96sfx+pKURJye/Plkm7/Ld4XuauZo84f+Ye539x+1fl8nabu775e02cx2SJo7Yv3r3X1IkszsnmD5eDlcv7dLeo+kN5vZUjUH4Tslyd03mtkbJf2bpD9395uOtDOlHVHu/q+S3tR6eJqkbSNWWS5pqpndYGbfbg2sXyTp+tb2t0k6p6z+AQAAIDXy/zwPHGb94Ve5XOranwcO1+/flPT11u//JunjenJfnyHpUUlxoH+HSv0ypLvXzOyzkn5bh06Y3yPpg5I+KekMSddJ+o6a/6dxUN3Metz9kMs069ev17Jlv7yZyAUXXKAVK1YU/AwAACjOwMBAt7vQdddee61WrVp18OG8bvblyfgy5ASxQtK17v6x1uyKd0mqSpKZPVfS5Wpe6P1nM7vc3f/PkTRWeuqIu/+hmb1L0vfN7Onuvru16F5Ja93dJd1rZlvUfKLD/75fiQbZkrR48WKtXr261L5nrtRHutIuAODY1t/f3+0udF1/f79WrlwpSTKzZB7P+CNHe8J4h6TPmdndkuqSbpL0O2Y2S9I/Sfqv7r7BzC6V9AMzu2nYVJSOlfllyEskLXb396t59brR+jnoMklnS3qbmS1Uc5L69ZLOl/Qvrakk95TVPwAAgOHqzpchJcndBySdNezxH41YfvAvEX/fenzpiOVPejxeOuj3r43Y9PLWf5cOW3edpM6+uBYo84r2VyR9xsxuktSr5v9BTDGzL7j7hZI+JekqM7tZzf+RvEzSbZJeYma3qjlf5nUl9g8AAAAoTWkD7dYUkd8NFl3YWr5f0sXB8reU1ScAAICIy0gdQeG4MyQAAICkBl+GPG4NDAwsrVQq75T02kajMaNSqeyU9PlGo/Gh/v7+dWW1y0B7HFSmTOlo/UaSM13ftq38tpMM4EqSFatk/z64MW0jywnPeC1OFKpv2Rqvn+UnJxnDbZcledl+IO6T70/Sj7wRltu9Fo2dO+NdZc8vUZkUZ6yn+0kyv6X82MyOsywbXXf8LN5/ktGcZqP35ZnfljzvbF/1xzbF+8meQ8IbcaazJFWTLHebHmdQN7JjPHkOmXpyLEn5MZgem8nrmmXFT7sv3s/kUw+XFHaobAzUSObVzvx5koOfvK71JCe+Oj2/wV390TinPrt/Ao5e5GgfvwYGBl5RqVSuWbRoUe/ChQt7+/r6NDQ0NHNwcPCNGzZsuHRgYOA1/f3915XRNkcUAAAAjkutK9nXLF++fOrSpUt7p0yZokqloilTpmjp0qW9y5cvn1qpVK4ZGBhYevi9dY6BNgAAmPBcprqX94PuqFQq71y0aFHvrFmzwuWzZs3SokWLeiuVyh+X0n4ZOwUAADjWNFQp7Qdd89qFCxe2nQPYWn5JGY0zRxsAAACSpDVr1uRfNOnQWWedNa6X8pe8++uH9P0zL5umvr6+ttv19fWpXm/MHLn9Ax945RH3n4E2AACY8NzFLdiPQ/vdNDQ0pCltwiGGhoY0VNL0Ho4oAAAAmRol/qA7bt1wQA+t39B2nQcf3qBbN3SehDQaXNEuUDWZaF/bsiWs98w/MaynEWpJrJYkNfbvD+udRpxVZ82MG5iWRFUlUWZZf9qpzo6jzzx5PWTx/ydWsufQhu8bCuuNXbs62k917px4/7t2h/Us1k2SrB5HAmZ9VTWJIkwiB9NjI3u9JfU89fSwXlt7X9xGh1GO2fGURaWlcYrKYyGzWMNOozB9DMd49prXNj4S1ntOmh+3vSc+fxt798XttokorExJIhI9fi9sdvw5l62/e2kcnfnMOZvTPm2ef2pYPzAtHqxs3BWf85YlWM6fF9arB+K+Zq+rlB+bWaTieZWLwvqNjavTNjA+XFzRPh594779etHiQZ08/0RFX4jcvn27Ht4wqOvu6/wzfTQYaAMAAOC4tGmv629v36O36y6dsmihTjtlkVo52nrw4Q16eMOg/vb2Pdq0t7Cp6U/CQBsAAEDcsOZ4dc/muq64aZdecfoDesH6Deoz15Cbbt1wQNfdt7+0QbbEQBsAAEAuS+8yimPfpr2uz/14SJ/7cTL9siQMtAEAAMQVbRSPIwoAAAAoAVe0AQDAhOeSGqSOoGAMtAtU3769o/Vrj20K65W+OG6rMdR5xFS6fhIz5rU43sp6k0Mlib3qWbQwb3xyfIcmf3xbWM+iAtPnnMTiKYl1k/IouOr06fGukti/+tbH0zYilX35KZjGPCbPO4s1bGzdFteTiMe2dsSRZZksVs6SSLksyrGxJz/2O5XF+FUWnBRvkEU/bonf63bnaadRn5odv6e+bUdYz+IU2/XJa/ExaMmd1Hx2fE40piafKcm45cITb0/79HcnnBHWtz07/iyY24gbmX9HElW5K65nx2U1iQOUJCVRi9Us2rLN5xC6zVQn7xoF43/dAAAAgBJwRRsAAEx4TB1BGRhoAwAASEwdQeEYaAMAgAnP3biijcJxRAEAAAAl4Io2AACApDpXtFEwBtoFyqLMsui4TKfrS22i4JIoM88ippIoPX9gfVjPYv9sWpu4wW3J7U8PxNGCqSTGr74njlnM3p92PHk9MtnrnUXp+VB+K9gsqi173o1tSbxkJZ5zmMVIqpr/Q5NFUsribTo9ziyJnaz0xMdZu9cv61MjievThkfiNpLnkB372THQTvbe2VD8OmVxfVkcZcXj11VSGjeXPT9/OH6dLDnGbdLZYb3dn+h7dsd96nk0jhy85Hk/COurel8S1rMYScti/LbnsZa+a3eywMNyFg2K7nNJjS7P0Taz+ZJul3Seu/+sq51BIRhoAwAAyLp6RdvMeiX9X0ljuMkBjlb8jQQAAKD7Pijp45IGu90RFIeBNgAAmPCaOdpW2o+keWa2etjPmw62bWaXStrk7t/s0tNHSZg6AgAAIKle7vXHze5+TrLsMkluZi+V9CxJnzOz33L3+EsROGYw0AYAAOgid3/xwd/NbEDSWxhkHx9KG2ibWVXSP0hapuZfZN7i7muGLe+V9GlJSyT1SXqfu3/NzO6QtKO12v3u/rqy+ggAACBJriemeACFKfOK9gpJcvcXmlm/pCslvWrY8tdK2uLul5jZXEl3mtkNkszd+0vsFwAAwCEaR8FX1xgDHV9KG2i7+7+a2arWw9MkbRuxytWSrmn9bpJqkpZLmtoacPdIusLdbyurj+2cV7mo423S3NkOc7GtL86KzbJfJak6Y0ZYb+yN83azbGBPsnBTWYxxm6zYTrOmK1PjTO4sOzzNy/Y4f7q5UZK5nOU9Z2004uzcLBM7zZluI9sme13TPO6++Hitb9vWcZ8q2bGfvR6V5B+zrK9Jxnr63Nq0nR0HnZ6/2XFZ374jrEuSJZnm6XHW4fGR5r4nx7eUP+9Gkg9tSaZ5Zd4JYX3qz+PPs6kvyjPQh+Ykx+ai+PPsZdN/HNa/Muu8sJ4l6nuWpV7LM/6zz63q3DlhPTtPs39zbmxcnbaNYrlLda5oo2ClztF295qZfVbSb0t6zYhluyTJzGaoOeB+j6Q9asbbfFLSGZKuM7Nl7n7Ip9z69eu1bNmyJx5fcMEFWrFiRVlPBQCAcTcwMNDtLhTu2muv1apVB6/DKblLEHB8KP3LkO7+h2b2LknfN7Onu/sTl0nM7BRJX5X0UXf/opn1SVrr7i7pXjPbImmBpIdH7nfx4sVavXp1af2+Uh8pbd8AAIxGf39/t7tQuP7+fq1cuVKSZGabu9ydJ2GONopW5pchL5G02N3fr+aV6kbr5+DykyTdIOmP3P1brfJlks6W9DYzWyhppqSNZfURAABAOvhlyO7P0cbxpcwr2l+R9Bkzu0nNKXHvkDTFzL7g7hdKukLSHEnvNbP3trZ5laSPmdnNaiaVXBZNGwEAAChaXVzRRrHK/DLkbkm/Gyy6sLX8ckmXB8svLqtPAAAAwHjhhjUAAGDCO3gLdqBIDLTHoOeEOMaqtmVLR/up9E2OFyRRUllclCTVtz7eWdtJxFQapZdFfU1LIuWS/UhSI+lrFlGYRr4l8X6VaXHbNjmOGZOkxo6dcT2J7spi5RpDcfxYpjprVrosi2rL3qNs/frO+LmNRdZ2xrJYNIv/Matv395Zf5LjWJJMSQzi7Nlx20msYXX69LCevd42KQuPk1SPowV75p8Y1muPbcr3FbWdRINmMZ+SpGo8J9WTSMXKSfPj/SQxgZakat6w7ey0Sz1DcTRj5ZH4+d24+8ywvulZcZ+mfafDGYlZVKTafI4n22Sva89pp3TWJ5SAOdooHkcUAAAAUAKuaAMAAEhq8GVIFIyBNgAAmPC4MyTKwEAbAABAYo42CscRBQAAAJSAK9oAAGDCa94ZkqkjKBYD7THoNMYvUzllYVhvDD4S1juN8JPyyLws8i2Ls0tlcXZj6GtR0oi4NslxadxhEt1lSZSZhuJyFuOXRcRJUmXO7LBeG9wY1nsWxcdTz/Rp8X42xsdZFmcnSV6PI/Ns0qRk/STbLZG9D2nM4pQkWk1Ko/Q8i4VM3uvsPWok9bHEcDYe6ywWspK83tl5nUZnjmGbRnL8VZKIwkpyiH9wwffSPv36zF8L65e94tth/WM/eXFYX/LVrWHdk9jJyuL4HPIHHw7rUn4cZH8uzuIfa23awPjhy5AoGlNHAAAAgBJwRRsAAEx43BkSZWCgDQAAIFJHUDwG2gAAAM6XIVE8/tcNAAAAKAFXtAEAwITnInUExWOgXaCeE04I61kcYOOh9XG9TeRbx7IYug5lcWKVE+aG9Sw6Tspj1BpJVGDWdhqrlUXyJbFazWVxG2mUYxLjV5k6NaxnkYPZc2s2HkeQZRqPb4vbmB1HC2Z9zSL8pDYxewfivnrtQNy2FxP757t2p9tkx4f1xMeB15LjL4t+zCL22kRb9ix9Stz2Y5vjDarVuI1t28J6/tzyY8mSNhp749ejemL8Oaeh+KSoJwmMvZZ/Nu0+Of5j68CmM8L6aSfEMX5DJy2I2/5p/Ho0HtoQ1i05VySp0hM/j0YWI5nsKzvOML6YOoKiMdAGAAATHqkjKANztAEAAIAScEUbAABAXNFG8RhoAwCACc9FvB+Kx9QRAAAAoARc0QYAABDxfigeA+1E26gli0/ELAquOn16WPcsnq7DODspj4zS/jherVNZ2402MX7pvpIYvzRmLItpy9bPYgLbRMp58jp1GkXoScRZ2m7D02X1rds621fyHHzHzs7WTyL52jcev7ZpLF+HcYDZOdT2nEjOo+w1z6LxGvs6e0/bRTY2BpPzJelTY38cX5jH+MWvn1XywUMa85gcH43kuKzMmhHWq/Gposfqu/I+JYfg3L44XvJHGxaF9dN2J+f1lPi8zqItG21iJDPpe5R8RmRtv2zy74f1b+77Qsd9wmE4c7RRPAbaAABgwiPeD2VgjjYAAABQAq5oAwAAiCvaKB4DbQAAMOER74cyMNAGAACQ5Ay0UTDmaAMAAAAlKO2KtplVJf2DpGVqfpn3Le6+ZtjyiqSPSlouaUjSGyTdN7Lm7mvL6iMAAMBB5GijaGVOHVkhSe7+QjPrl3SlpFcNW/5qSZPd/flmdq6klZI+G9SGb1O48yoXlbl7SVJ9V54XG2mb4Z1o7NnT0fpZLnGnfc2kud7K+5rlyGY6Xb9d5nKnucRprvPeOOc3y/xul1mdZlBnWclJG0W9p21Z/MexNOs8eZ0y2Xtt7TKrkzaybHTLspX3xkHQWZZ6u+OsKFaJP7orYzmvszz/5HWqTJ8Wt5Hka5/wk/h4fbQeHxuS1Eje1hm98Wu+f3e8QS3uqipJlnWR711+P4T4WB7L/RNQLCdHGyUobaDt7v9qZqtaD0+TtG3EKi+SdH1r3dvM7BxJ64IaAAAAcMwp9cuQ7l4zs89K+m1JrxmxeKak7cMe16OamfW4+yG3kVu/fr2WLVv2xOMLLrhAK1asKKzvAAAcrQYGBrrdhTG79tprtWrVwetwmtfNvozElyFRtNJTR9z9D83sXZK+b2ZPd/eD97LdIWn4PXsrUS0aZEvS4sWLtXr16iPu35X6yBHvAwCA8dTf39/tLoxZf3+/Vq5cKUkys81d7s4wxPuheKWljpjZJWb2p62HeyQ1Wj8H3SLp/Na650q6J6kBAACUzt1K+8HEVOYV7a9I+oyZ3SSpV9I7JE0xsy+4+4WSvirpPDO7VZJJep2ke4MaAAAAcMwp88uQuyX9brDowtbyhqS3BMujGgAAQGlcpI6geNwZMtGz5LR0We2BB+Nt5p8Yr//Ypo7atr6+sF552lPybe4fDOv1bdviehL3VZ0xI6zbrJlh3XfG+7GpeURhZc7seF9z4zbq9/wsrFfnzonX3/p42namXcxe2EaHkXnVWbM6Wl+S6tu3H36lYbK4ufFQSaLxMulxNjk+9tXwsOxjiD7LXqdq1nY1nmGXPoc2x74fiI+zjo/Znvij2xuNsD6WuM3sdbJp8fOrLjo5rD/8m3FfH6rF52877zr5hrD+o02LwvrkX8QRj43k9asm0YW+L44DlCSrJNGWJ8TPr5Lsy3fH70O1TYQlCubNiD+gSAy0AQAAxA1rUDxuwQ4AAACUgCvaAABgwnORo43iMdAGAAAgRxslYOoIAAAAUAKuaAMAAIjUERSPgXYii/Bry4r5k1MWyacfJfVxkMb4zZge1mvrN6T7qvzq08N640c/CevV2bPjHdXjKLPq9LhP7XQc15e0ke0ni3Wz009J2+hJosZUq8XlRx9L9xWpTIlj2hp740i0dttkPHt+994fr7+7HvcpiaCrnrUsb3zNz8NyzxlL4/W3bovryXmdnqc7d+Z96lCnr3d2TlRmx9GZUv7a9pwWv3eNR+K4Ukv6Omtd/Pq9cmoeR/n2X42XLe2Jz7vNG+P4zPm1+HysJOdWY0Ec0Zodr5Lk9fiYzY6nenJ89Jw0P+5T8tl7XuWitE83Nq5Ol6E95mijaAy0AQDAhOfOQBvFY442AAAAUAKuaAMAAIhbsKN4DLQBAADElyFRPAbaAAAAYo42isccbQAAAKAEXNFO9CxamC6rbRiM60m8WqVvcljPYqG8FkfBVSZNSvuU7iuLnkpk0VOp7dvDcs+J89JNalmM36w4osv3xlFfjaE8Hqxsvn9/WM+eQz15nSq/yGMk6x0+v+z4aCR9bRfjl2rEf1fNjpvKT9d21KeOu/Oz+zrf5sGH43qHfcrOa1XaXBHLztPkde30PapMnRrW20U/psfNxkfj9WfFUYG+Iz4Gdi1Om0595AVf6Gj96pQ48lJTks/epK8aw/Fq1Wq8TfLeZet3Gs+J4rmMK9ooHANtAAAASUzRRtGYOgIAAACUgCvaAAAA3LAGJWCgDQAAIDF3BIVjoA0AACCuaKN4zNEGAADoMjOrmtmnzewWM7vZzM7qdp9w5LiiDQAAoK7fGXJFsw/+QjPrl3SlpFd1s0M4chN+oJ1l4WZZ2WNhSZarhobCcpajPZYs18qUKe07N7LtrA2L//iR9bX+eJwb3U5j16646ey5Je9dlq+d7UfK88atpzduI3mdLNlPp31tt02nmemd5mu3O2bS4yPROBDnG3fap+x9aCd7v7PM6qxP2frZsd9pdr2Uv+Y2KX7ennx2ZKOEdhn86XvUG//z0Ni9J6xbb9zX+Xckr8fr0y5pW31aWP/xgW1h3Tf3xTtKXo/0PUqOmeqMGfH6avNeJCw7zjrMTG/3eYaxcXV36oi7/6uZrWo9PE3Stq51BoWZ8ANtAACA5ki71IH2PDNbPezxJ9z9E0/qgnvNzD4r6bclvabMzmB8MNAGAAAo32Z3P+dwK7n7H5rZuyR938ye7u67x6FvKAkDbQAAAHV3jraZXSJpsbu/X9IeSY3WD45hDLQBAACkbudof0XSZ8zsJkm9kt7h7p1N3sdRh4E2AABAl7WmiPxut/uBYpU20DazXkmflrREUp+k97n714Ytv1TSpa2HkyU9S9LJkr4taUerfr+7v66sPgIAADQZN6xB4cq8ov1aSVvc/RIzmyvpTklPDLTd/SpJV0mSmX1EzUH5Pknm7v1Fd+a8F10ZL2gTr1aYLHJr6tSw3ugw5knK48G8EU/vauyJI7rSmLEkorC+9fG43ST6rJ3KrJkdtdFp5FsWqyXl0VqdPg/rS2LGakmEWhLhJ+XRf3lUYDFxc+2Ov04jxbJYuSy2rpLESGavRXX69LTtehIXKSXxiMnaeRRh5zF+qSRCsDHU2WdBpxFx7aTxj8n6NntWWN92eucxdL83PT7nP73jKWF97t1Jr7btCMvZZ0FjV/yds3bnRPY5nkVhdvoetYtmzJxXuSis39i4uuN9TTjcgh0FK3OgfbWka1q/m6RwpGFm50h6hrv/FzN7nqSpZnZDq29XuPttJfYRAABAcm7BjuKVNtB2912SZGYz1BxwvydZ9QpJf9n6fY+kD0r6pKQzJF1nZsvc/ZBB+vr167Vs2bInHl9wwQVasWJFcU8AAIDjyMDAQLe7IEm69tprtWrVwfuyaF43+wKUrdQvQ5rZKZK+Kumj7v7FYPlsScvc/Tut0r2S1rq7S7rXzLZIWiDp4ZHbLl68WKtXrx5ZTl2pWzp/AgAAHCf6+/u73QVJzX6sXLlSkmRmm7vcnSdj6ggKVuaXIU+SdIOkP3L3byWrvVjS8GWXSTpb0tvMbKGkmZI2ltVHAACAX2LqCIqVfa+lCFdImiPpvWY20PpZZGZfGbbOMkn3DXv8KUmzzexmSV+SdFk0bQQAAAA42qVXtM3s1jbbmaSGu78wW8HdL5d0ebDowmHr/M2IbfZLurhNuwAAAOVg6ggK1m7qiEn6vTbL/qn47pTo1ju71nR9+/by20ijzDqTxV5lEXtF6rSNyvRpYd1PWxjWG3f9JN1Xz4KT4wX748i82pYtcRtJbGLPU0+P97P2vrDeThaNV03qNnd23PaDh3z14bA6jQq0yXHcYVHnRLvjvjorjptTEnmpJLqwvm1bWO854YSwnh0b7WTxhZ1KI0OT47LdNp2qDcaz/KoH4ki+sfj8w88L65vPjc/T+d+Iz4n6xvg9yj4HPIn9kyRV4j8MV2fMCOtZVKqS96H2wINxs+3etyRaEKPAQBsFazfQ/qSkx7Lbf5rZa8vpEgAAwDhzScT7oWDt5mifLekeM/uomS0fudDd15XXLQAAAODYlg603f0dks5U85bof2Vmt5jZ682smL8zAgAAHEXcy/vBxNQ2dcTdD7j7Ne7+SkkXqXkTmYfGpWcAAADjyUv8wYR02BxtM5ss6bcl/YGkGZL+pOxOAQAAjDvmaKNg7eL9+iX9oaSXSPpXSf+fu68Zl14BAAAAx7h2V7T/QtInJL3F3YfGpzvjr2fxorBeW7+h/LY7jAdL48rUeVxatq9sP9XZs+P1k+izsciitWobH+mo7Z7e+LBOQt3atpFFdHX63o0lxq86d05Y7zhqcefOjlbvedoZ6TJ/ZFNYz96L8Yi2zKTHcvKednosZ+/1/lc8J91mynd/GtZt6pR4g2Rip58UH3/1NT8P6z0nzU/7lMUdNnbHkYCVU+L4zMrOOALPOkuEbOvbz/haWH/O/oviDZLIy0z2OdAuSi+7/pm9p42dcSRl47H43BpL/GJ2jOPwjCkeKFg60Hb3fkkys6qZvV7SaWp+MXKNu28en+4BAACMA+ZSowSjuQX7/1VzkH2emnO0P1dqjwAAAIDjwGgG2kvd/X9I2uvu10rK5y8AAAAck6z5ZciyfjAhHTZ1RFKPmc2TJDObofZTXQEAAI5NTB1BwUYz0H6PpFskLZB0m6R3lNkhAACArmCgjYIddqDt7v8uaZmZneju8deiAQAAADzJaG5Y82ZJb5Y02aw5x8jdn15yv8bP3n0db1LpmxzWG0PxvrL1s3iwbH3fV1zKYhZ9lrWtemcZXZVJkzrtUsfRWlaJv2KQvU6VKUmEmiQ14ssYNmN6WK8Nbsz31YEsJlDKj4/S7Yjjx6TOI/CsJ45XsyR2rbEnjpQrkg/Fx0fW18rkvrBe3xW/TlNu/lnadraNknp27DfWxMFP2XlX35JHQnb8Xtwf3xy4cnIcIXhgWtp0Ybb9JD6P5j12Z0f7yV7vdp8dnn021pJ6vcPZl7Va3G7ymSVJ6uWy7Jjx0qFgo5k6crmk8yV1GN4LAABwjHDxpUUUbjQD7bslPezuBd52AAAA4OjCDWtQtNEMtL8t6T4zW6fmTbDc3X+j3G4BAAAAx7bRDLTfLOl3JW0rtysAAABdxBVtFGw0A+31kn7o7uRnAwAAAKM0moF2n6S7zGyNWv+v5+4Xl9orAAAA4Bg3moH2+0vvBQAAQJfxZUgULR1om9n17v7y1g1r0uXldW1szqtc1NH6vn9/5410OIsmy9fudP0043oMqtPjfOg057fTPrV5jSozZsQLdu4My0VlK7fN9rY4k7uxOck6T/ZlSQ5vlj89HlnZ6fPuiU//xuPb8n0l77f1xvvKjievHUjbCPefZFx3ei5KkqrVuN6I84qz55CeQ8lxPBa+v7PXqZF8nrXNgU7aqM6aFW9Qjc8VDcVtz3yo/FmH81fHo6PKiXG+tiVZ8dl52un70Nwoft5Z7nZ2nqY53clnliRZcm5n/z7e2Lg63deEQ7wfCtbuivZvmNlgsswkzSmhPwAAAOPPxZchUbh0oO3und/aDwAAAICk0c3RBgAAOP5xRRsFY6ANAAAgvgyJ4nU80DazPncfKqMzAAAAXcNAGwVLv7ZsZl8a9vt/G7boutHs2Mx6zewfzey7ZvYDM/utYJ07zGyg9fMZM6uY2cfN7Hut2lM7ejYAAADAUaLdFe35w35/paSVrd9Hm33zWklb3P0SM5sr6U5JXzu40MwmSzJ37x9Wu1DSZHd/vpmd22rzVaNsr60sqqq+fXvH+0ojtLLos0lxNJklcU5p5FuluNihNMZvHDSyGL+xRC120m6B+8+iuBq7dhfWRmGSqC/VO49dS6Mnq3GsYaaaRDx6LY7Ya+zdG7ebxClKeSykJXFpWURhZUp8Xo8l8i37jEgj3BJZ3KElnxFZ3FtzYbxN9tnYc9L8sK5a/Bxm3bM1b7sgvbvjY7n20Pqwnv17kEZhtovSS46bxr7sD79JvF92TiQxgY0kolCSfChuuzqXsLDD4oo2CjbaqSPDP4lHexheLemaYduP/Bd0uaSpZnZDqx9XSHqRpOslyd1vM7NzRtkWAADAmJkzRxvFazfQ9uT3UXH3XZJkZjPUHHC/Z8QqeyR9UNInJZ2h5pSU70gafhmlbmY97n7IZa7169dr2bJlTzy+4IILtGLFik67CQDAhDYwMDCu7V177bVatWrVwYfzxrVxYJy1G2g/w8y+qObV6OG/P320OzezUyR9VdJH3f2LIxbfK2mtu7uke81si6SqpOF/P6tEg2xJWrx4sVavXn1I/Up9ZLTdAwBgwuvv7x/39laubM5GNbPN49r44XBnSBSs3UD7d4f9/vHk95SZnSTpBkl/5O7fCla5TNLZkt5mZgslzVRz2sj5kv6lNUf7ntG0BQAAcMSYOoKCtbsz5L8f4b6vUPM27e81s/e2ar8v6e/c/UJJn5J0lZndrOahfZmk2yS9xMxuVfPq+euOsA8AAACjwhxtFK20G9a4++WSLg8WXdhavl/SxcHyt5TVJwAAAGC8TJg7Q44lxq9TefRZHA1Vz2L80gby/9VO49KS2LAs+qxT2XMeD1lUWvo+jCEKrjp9eljvNB4xa7vd+9Dp80t5fNwU+t51GBVYTyIeO2W9ccxd+43i8zGNEJwyJV5/DK9ftk119uyw7snxUZkdx9PVHtsU77+vL+9U8t6lUXeN5L3uSyIvp8fHcZGsFh/jPQsXhPXa4MaO9p8dA1L+WZCd81kEYxrrOgbVE04I69kxjmG4oo2CTZiBNgAAQIp4P5SAgTYAAIDEFW0ULr/dFQAAAIAx44o2AACAxBVtFI6BNgAAgJijjeIxdQQAAAAoAVe0x0EW4daz4OSwXtv4SFivzJ2dtpFt06meJBaq0OipuXPCen3r44Wsn2kXpddz0vyw7nviOKzsdfL9+8P6WOLssii4nlMXxxvU4ijHTqPMqrPi6Dgpj8ksO+axZ/6JYd335u1mkXmNZafGG3z/7nj9DiPRsr5Kefxefdu2jtqwofg4S9efEcdUSpJ2dxb1WdsU3zH7wMueE9aHZlc72v9Y7Dsh/qdsaiW+lpQd4+nxPYZYPJucRCoeSF6P5LMjO56yY0mSfGgorGefgedVLgrrNzauTtsAMDoMtAEAACTmaKNwDLQBAADI0UYJmKMNAAAAlIAr2gAAABJTR1A4BtoAAAASA20UjqkjAAAAQAkm/BXtsUQnZXFzGetJXuZqZ7FXRUX4tW0jifGrPuNXwro9tjXeTxIBJnUey2eTJ8f1554d1qsPJ+9dEvUlSbUNgx31qVM9Tz09bnftfR3vy3ftDuvW2xu3/bQz4h0l0YW1h9anbfcsWhgvaDTifRV1zJrF9VMX5JtsSuIi74/jDmvJfrLox9qjj8UbeH5JrNOoyizK0bfHcZE9py8J6405ebxf4/b49ahOj7dJo/Fq8TEwbTB7ZYtTmxIfH9m5ksX4Ze9P23jETPa5n0UzJhGg2b9FtfPOyZv+eXzeVZN/c2xBfIxPNCa+DIniTfiBNgAAgCSmjqBwDLQBAACI90MJmKMNAAAAlICBNgAAgNScOlLWTxtm1mtm/2hm3zWzH5jZbxX91NAdTB0BAACQujlH+7WStrj7JWY2V9Kdkr7Wtd6gMAy0AQAA1NU52ldLuuZgN5QHIeEYw0AbAACgfPPMbPWwx59w909IkrvvkiQzm6HmgPs9XegfSjBhBto9J5wQ1tvlZWc6zYEeD9UZM+IFvclbvP9AWK7v2hXXf3xvMf2RVE/yYitTp4b12mCS87sz7msjeW42Kc6ZHotKX5zt3RjaF9azvOzKlClpGzYlbiM7/rJ9NbK850SWnyy1yfA+Oc6jr+6dHdb3PS/O9u795g/jhmv1sNzuuKw+88x4m7t/mm4TNp28fulnSpsM+Ux27Pvjcd5zdg719E0K6437Hkjbzs5VT86jyrz4effdmbQxe2badlFOuDt+PepnnhbWezbGWeDaHR/fvm1H2nanmdy+byjuU4d57ZP+/e60T558FljyXuw9NXk9JqJyr2hvdvc0AN3MTpH0VUkfdfcvltoTjJsJM9AGAABIjeJLi2Uxs5Mk3SDpj9z9W93pBcpA6ggAAEB3XSFpjqT3mtlA6yf/kyeOGVzRBgAAUPe+DOnul0u6vDuto0wMtAEAACRuwY7CMdAGAAAQt2BH8ZijDQAAAJSgtCvaZtYr6dOSlkjqk/Q+d//a4Zab2R2SDmYp3e/ur4v2P+eOp+i8ykWj7o/X4uz3nvlxLNlYYv/K1i4Kzutx/JnvjePmvBZHd6WRZVu2xH1KYslUrcb1duqNsGw9cSxfGhOYRO9l67fbxqbFr3lhEY9m6aI0xm9SHOGmRmeXYtL3rqfNx4LHbTQeWB/Ws4jCNMYv0UhiBdPXQp3H+HUqPSfaRTYm50UjiXxzxedpJovka6eRfEaksgi87DPo4Q2ddqljNpS0fdc98QYLFyQ7Ss7HRvzZJOUxflmEanZ8+O49cZeSz7/G/v1pnyrJcZYdf5M35PGFEw5XtFGwMqeOHO52oocsN7MbJJm795fYLwAAgCfrYrwfjl9lDrQPdzvRaPlySVNbA+4eSVe4+20l9hEAAAQGBgZK2e+1116rVatWHXw4r5RGxsBaP0CRShtoH+52osnyPZI+KOmTks6QdJ2ZLXP3eN4HAAAoRX9/f2n7XblypSTJzDq/lSpwDCk1deRwtxMdudzM+iStdXeXdK+ZbZG0QNLDZfYTAACAqSMoWplfhmx7O9Fk+WWSzpb0NjNbKGmmpI1l9REAAOAg4v1QtDLj/aLbiS4ys69kyyV9UdJsM7tZ0pckXca0EQAAAByLypyjnd1O9MLDLL+4jP7Ut28vbF/VGTPCehat1RjqLD4rjV1L4u8kySYnMWdJrGH2P+1ZZFkq2f9Y4v3Kfp16Tsy/c9PYHsdbZRF71dmz4/W3bUvbiLSLY0tj4pIYNXn8vLPXKW27zXGWtpGdE9nx0aEsjlKWXyvIov8se12TaDebFMerZdGZquR98gPx87BK0nby3jWSqMo0BrFN5GBj7950WaS2KZ5Sm30uVmbPCutP+fob0jbuf+Unw/r1e+P39MC8+Pn1JfGtaZTevLnx+tvzaND0y3PJcZO93o09cZ+y87cyuS/v08z4vcgiBy05ZrMI3RsbV6dtH/O4oo2CcWdIAAAAiYE2CsdAGwAAwJmjjeJxC3YAAACgBFzRBgAAkJg6gsIx0AYAABBTR1A8BtoAAAASV7RRuAk/0O552hnxgp1xTJYk1TYMltSbpsr0aWHdh/an22Txhe1ivYrQ2B/3qdImdi1TnTsnrGcRe1kcVhbjl8WSSW0itLLn4cV8Gqexdcrju7Lortqjj8U7St6jTtuVpPquXWG90WksZIey96eexNxJUmX69LDuWbxah69Tpuepp6fL6mvvC+tZFGGncZHVWZ0950ItWRTXt8bRmWd8Nj/29cq4/IH7XhHWp65eG9a90SaqMlBbd39Yr/RNTrfJIiwtiTjNovSyYyD7rPF2x2sWX5hFVe6Mz+s0QhXAqE34gTYAAIDE1BEUj4E2AACAi6kjKBzxfgAAAEAJuKINAAAgcUUbhWOgDQAAJjwTc7RRPAbaAAAAEle0UTjmaAMAAAAlmPBXtK/7yV+F9fMqF41zT36p9timwvbV6DA/t+fUxWG99tD6ztod2pe3cfqSuI37HuiojUy7vOxMlsmdavP8ipJlViurF8Smtsle77DtLIe309e7XV52xuv1eMEYMt4j1Vmzwnotycput42SvOee+SeGdd+V5/xHLMnmlyRtTbLwk+zo7Ny2R+LzLnsfhk7I89qv3hW/TgPP+Lew/srZF4T1xrx4P5WNSe57cpzZtDbnRD1+77IM9OqMOAe/0xz3dvdIsBlxnnqWqZ8dH/WS7xlxNLKC7pEAHDThB9oAAADE+6EMDLQBAADElyFRPOZoAwAAACXgijYAAIDE1BEUjoE2AACAmDqC4jF1BAAAACjBMXtF+/Ffu1+rV68ubf83Nq4ubF/djArsVKcxfpmep56et5HEn/WcsTRe/xfrwnrlV58e1hs/+slhehe0veDkuO2Nj3S8r7J12tdNb3tBWD/xo7fG+2kTL1l59llhvXH7mriexPh1up+xqMyZHdZrgxvDes9pp8TrP/hwWLdJvfF+TpyXd2paHKNWe+DBeP0OYw2rc+eEdUtiFiVJWx8Py1mMX8+ihWG9lkTBZefpw78Vx+JJ0kunxu+RFD+Ph343fu9OuX5rWM+OgUw9eY3G4sCvPjWsV276UVjvOeGEsF7bkkQUSqpUk2toSYxfp9Gq2b9pRf672TVc0UbBjtmBNgAAQGGcqSMoHgNtAAAAiSvaKBxztAEAAIAScEUbAABMeCamjqB4DLQBAAAkyRlpo1gMtAEAAMQVbRSvtIG2mfVK+rSkJZL6JL3P3b82bHlF0kclLZc0JOkNku4bWXP3tWX1cbyMR+RR2RGCWWRZbdPmuJ5E+EltYtSSGL+eJafF648hxi9zNMb4Vc/8lbBe++m9He0ni/Ebiyx+r+dpZ4T12s9+0dF+OlWdNStfWK+H5fRYTmL8epbFcWy1n4/hoyk5XwrTiEcJtfUb0k06jo+b3iYqMOpScp6e8fp8m+f+01vD+tyZu8P6wr+Jj/H4CBgf1WeeGS9IYvwy2ftQfUb8+SBJtR/HnxG1884J65PmTA/rO0+P69Ouvi1tG8CTlXlF+7WStrj7JWY2V9Kdkr42bPmrJU129+eb2bmSVkr6bFB7VYl9BAAAaCaOcEUbBStzoH21pGtav5uk2ojlL5J0vSS5+21mdo6kdUENAACgdJbfRwkYk9IG2u6+S5LMbIaaA+73jFhlpqTtwx7Xo5qZ9bj7yEG61q9fr2XLlj3x+IILLtCKFSsK6j0AAOjEwMDAqNa79tprtWrVqoMP29xKFTj2lfplSDM7RdJXJX3U3b84YvEOSTOGPa5EtWiQLUmLFy8u9Rbsx5or9ZFudwEAMIH19/ePer2VK1dKksys5C8udIipIyhYaTesMbOTJN0g6V3u/ulglVsknd9a91xJ9yQ1AACA0pmX94OJqcwr2ldImiPpvWb23lbt9yX9nbtfqOaV7vPM7FY153C/TtK9QQ0AAKBcLnK0Ubgy52hfLunyYNGFreUNSW8Jlkc1HEYUITgwMKArfyOeUpJG7CURZ1mM31hkbVT6Jod137I1rFenx9FTGZs5I13W2LY9rNukSWG9vm1b3Keznxavf8/P4vXbxNPVO4zx61TPGUvDehazKLV5zbdui9c/a1lYr6/5eWf7T9S3x++bJPX0nRgvmJ285skx7o92duz3zE/abaP22KawXpkyJa7P7yyisK1JvWG5OiM5X4YOhOUshnPPmfPD+pSH8vdu5k1xhOCcO+JvqlWSKEw9PBjXzeJ6Jf4jb7vjLFO/+6cdrZ8d++nn1tYd6b6y42byz+MY06GlJ4X1WTc/ENa9w/MUmMi4YQ0AAICY4oHiMdAGAACQ+DIkClfalyEBAACAiYwr2gAAYMIzMXUExWOgDQAA4E7qCArHQBsAAEBc0UbxGGgf56LYv7F4+bw3hfX61scL2b8kNYb2hfXsiwTZ+lksWW1wY8d9qiZxX1kbWYxfpl1smPXEsWtei+PVOtVYH0efZTGLklTftStekNQru3Z33K+ITU1i7ur1dJssMk9ZlN7UOFLOqtX2nRttu8pf2+x48uT5+ZbOzrtKElMpST40FNZt+rR4g/37k0bic2Xq/dvCemNaX9qn6YPx87Z94Y2CZdviqLtG0tfsufkYjtf0fPE4ijDrkx+In1vjsThesjI5f/2sJ/mnvRGPIic9sjNevxb3KesrgEMx0AYAAJBIHUHhGGgDAACIqSMoHgNtAAAAVzq9BhgrcrQBAACAEnBFGwAAQGKONgrHQBsAAEDM0UbxmDoCAAAAlIAr2hiVIvOyM2mOcZZlPW1OWPd9cS5wZUqcxSxJjb17430lmbdZFm6mOmtWvP+kr1KbXPEkt9d649M5zb4u8Es/aU5zh3dZy17vLLe3Mnd2uq/sPc362tizJ6ynx1+Wfb0/zzm3SXE2evY6ZX3qlCXnltT5ud1z6uLOGt+6LSzbxvwcmr4x7m+WhW9z48+CLC/bJsfn0Fg+59LzNHvNs8+UZD89CxfE67fLUq/Gx6wlq9ue+FxRX5zV3RiKs73Pq1yUdqmoezqUjjtDomAMtAEAAMTUERSPqSMAAABe8s8omNnzzGygmCeEowFXtAEAALrMzP5E0iWSdne7LygOV7QBAMCEZ5LMvbSfUVgn6cJynyXGG1e0AQAAJKlR6t7nmdnqYY8/4e6fOPjA3b9sZktK7QHGHQNtAAAAabRXnsdqs7ufU2YDOPow0MaoZNFMv9n/V2G9ctOPOm6j0yizrE/tIqYy1dmzw3p927bO9jN9eljPYt3q27d3tP92shi/nhPnhfXapjiiayyyuMM0miyJ3pPFs9my6DNt2XrYvo1WzxlLw3rtF+vCehrfVslC1CSv18N6duxnEYJKIgezuMh2x3E1icbLou5qD62P99PhOdSz5LS0T97h+zoe8aOZ7HnUHniwkP03Ht8W17NzQvnnUCN5XSuL4vO09uB98f4L+rwEJgIG2gAAAB2kgwCjxUAbAABA3vUb1rj7A5LO7WonUChSRwAAAIAScEUbAABA3BkSxWOgDQAAIHV96giOPwy0AQAAXLJyc7QxAZU+0Daz50n6a3fvH1G/VNKlrYeTJT1L0smSvi1pR6t+v7u/ruw+Yuy+NXBFWB9LxF6nMXRjaSOTxVJlEW5ZHFsWsVdJYt3SiLg2bVTmzo7X3/hIWM8i34rUc9L8sF4b3Biv/7Qzwnrj/ofj9RecHNZ9V5s7FSeRg1kUYSOJ8ctYErFns2el29TWbwjrlSlTwrrvPxDWs2i36pm/Eu//gfh1lTqPxqvOip9fp9Fu7eLvsuOpZ378GdHY+GhYr0yfFtZ9b/z62eS+eP+78+hR31ZMRGf2uqqafJUqi8hUmxjJZBtLjoGsTzY1Pl5FvB9wiFIH2mb2J5IukXTIv4bufpWkq1rrfUTSpyXtk2QjB+UAAAClY+oIClZ26sg6SRe2W8HMzpH0jNZtSJdLmmpmN5jZt82MiBsAADA+vMQfTEilXtF29y+b2ZLDrHaFpL9s/b5H0gclfVLSGZKuM7Nl7l4budH69eu1bNmyJx5fcMEFWrFiRSH9Pl7s2rVLAwMD3e4GAGCCG/5v0bXXXqtVq1YdfBjPB+qSkm/Bjgmoq1+GNLPZkpa5+3dapXslrXV3l3SvmW2RtEDSIRMMFy9erNWrV49bX49FAwMD6u/v70rbV+ojXWkXAHD0Gf5vUX9/v1auXClJMrP4SzjAcaLbqSMvlvStYY8vk3S2pLeZ2UJJMyXF36YCAAAoEle0UbBxvTOkmc01s68MKy2TdN+wx5+SNNvMbpb0JUmXRdNGAAAACuWSGiX+YEIq/Yq2uz8g6dzW71s17MuR7v43I9bdL+nisvuE8t3YuLrjbV7x1P8vrPc04isMtS1bOm4jU509O6x3GlmWyWK1xmRaHgkYqe/c2dH62Wsh5a9H7dHHwnrPooXx+j/7RUd9yqILx0PPCSfECyoWln1nHPHYTqfHR3q8/vTejtvumX9ivCC5spfFbeoFz4rrt97ZcZ/UiEcltXX3h/Xq3Dlh3WvxdZosxi99bm30LFoQL+jws8Mm9Yb1rE/p+yZJSaxh7dS5Yb16d/y6ZtGC3ibuEMCTdXvqCAAAQNeZnC9DonAMtAEAACTmaKNwDLQBAAAkBtoo3Lh+GRIAAACYKLiiDQAAcDB1BCgQA20AAABxZ0gUj4E2jhq1+x4odf/23LPTZfUf3FNq22Py/OVhufa9u8L6jt9/flif+YXvddRsu0jD+m8+O6xXv3V7WK9tGOyo7d0XnRvWJ28+ELf7nbjdInUcI5m8b5Kk5L3r1NbfOjOsz/pcZ++1JNUe29TZ+i89J6z3/L/O7tS7/+XPSZftnlEN69OujqPudr5kWVj3ZHLktKtvC+vbvn5GWJ/9yjyOcs/SODJv0i/WpdtEOo0WbPe+VRYk0X9JVKo/ZVG8+u1rOupT9hl0TGGgjYIxRxsAAAAoAVe0AQAA5FzRRuEYaAMAALgYaKNwTB0BAAAASsAVbQAAAIl4PxSOgTYAAICI90PxGGgDAABIzNFG4Rho46hxY+PqjtY/r3JRIfsZSxuZItsuzD92t/mBgQH19/d3txMTxVXd7kDJvjT2TZ90HHa6H6YTABgjBtoAAACu9KY+wFgx0AYAACBHGyUg3g8AAAAoAVe0AQAAJK5oo3AMtAEAACQG2igcA20AAAC+DIkSMNDGMSuL0isy9q/TNsYjchAAABwbGGgDAADIJSc0HcVioA0AACAxRxuFY6ANAADAHG2UgBzt49i1117b7S4AHIc4KnAcAugGBtrHsVWrVnW7CwDHIY4KHIcYFffyfjAhMXUEAABAYkCMwh2zA+3bb799s5k92O1+HOXmmdnmbnfiaGFmE7LtowDHIY4GHIdHp9O63QGgTMfsQNvdT+x2HwAAwPGCKR4o3jE70AYAACiMS2qQo41iMdAGAACQuKKNwpE6AgAAAJSAK9oAAAASV7RROAbaxzAzq0j6qKTlkoYkvcHd1wbr3SFpR+vh/ZJeP5rtgNEYzXFoZr2SPi1piaQ+Se9z96+NPDbd/XXj1nEcFw53/EXLJd3XbhtMVM6dIVE4BtrHtldLmuzuzzezcyWtlPSq4SuY2WRJ5u79w2oXHm47oAOv1uGPp9dK2uLul5jZXEl3mtkNGnFsAmPwarU//qLlnz3MNpiIXHLny5AoFnO0j20vknS9JLn7bZLOCdZZLmmqmd1gZt9u/aMymu2A0RrN8XS1pPe2fjdJNcXHJtCpwx1/0XI+AwGMC65oHyPM7PWS/nhE+VFJNw57XDezHnevDavtkfRBSZ+UdIak6yR9R9L2w2wHHGKsx6G772ptP0PSNZLeo+DYNLNlHIfo0Ey1/zw7ZPkotsFExdQRFIyB9jHC3T8l6VPDa2b2IUkzhpUqwT8U90pa6+4u6V4z2yKpOortgEMcwXEoMztF0lclfdTdv2hmfTr02Fwg6eHSngCORzvU/vg7ZPkotsFExZchUTCmjhzbbpF0viS1/ux+T7DOZWrOP5SZLVTzSs71o9gOGK3DHodmdpKkGyS9y90/3SpHx+bG8egwjiuHO/6i5aP57ASAI8YV7WPbVyWdZ2a3qjnv9XWS1Pqy2Sfd/UI1rz5eZWY3q3nfq8sk3SbpJSO3A8ZoNMfhFZLmSHqvmR2cq/0qSR8bfmxyVRFjcMjxN+LYi47Pe4MaJjp37gyJwpnzZxIAADDBzarO8+dPW1Ha/r+586rb3f2o/+LtmjVrChsYnnXWWVbUvkZjybu/Xuig9oEPvPKI+88VbQAAAEnOFW0UjDnaAAAAQAm4og0AANC8Y023O4HjDANtAAAAFznaKBwDbQAAAEniFuwoGHO0AQAAgBIw0AaOYWY21cxuMbOnBcvmmdkNZvZdM/uSmU1t1d9oZqvN7DYzu6DduqPsw9+Y2ffM7Idm9sZg+blm9v1WP/+8VauY2cdb2w2Y2VOzddu023bdIp+/mV1qZg+Z2TvN7C/M7N5WvwfM7B4z+7Nh630g6MsDZjZ52OOntbZ9kZndaWb/PJrXGkB5XJI3vLQfTEwMtIFjlJmdI+kmSUuTVf6HpC+6+3+Q9CNJbzazkyW9XdILJb1M0vtbt0I/ZN1R9uElkp7q7s+X9CJJ7zKzOSNW+7iki1vLn2dmvyrp1ZImt7Z7t1p3iEzWzRxu3aKf/xfd/UOt3z/k7v3u3i/pHEmXmdn8Nn0NufvNkt7R6XYASuDenDpS1g8mJAbawLGrT9JvS/pZsvxFkq5v/X6dpJdKeq6kW9x9yN23S1or6ZnJuqPxPTXvNio1LwhVJR04uNDMZkrqc/d13rw71jdb+36iPXe/TdI5bdY9xCjXHY/nL0knSOqVtLeDbQDgCdlf+XDs48uQwDHK3W+RJLP0xlUzJW1v/b5T0qwRtax+sDaaPuyTtM/MeiV9VtIn3H3XiD7sGNHe6UE/6m3WzZ7b4dYt8/m/08z+s6RTJG2Q9AZ339nmvYjwt2TgKNPFKR6vVuuvfGZ2rpp/5XtVtzqD4nBFGzhGmNn7hs0Lro5ikx2SZrR+nyFp24haVj9Yy/qxqtWHv2s9nqPm1eCfuPv72/ShXT8qbdY93HPL1i3l+bd8yN1/XdJFkk6WdO9h1t+r5l8gDpouroADR5/uTR055K98ZT9VjA+uaAPHCHd/T4eb3CLpfElXSXqFpO9K+oGkK1tfzOuTdKakNcm6WT8uOPi7mU2R9C1JK939C8G6O8xsv5ktlXSfmvOi/1LSYkkrJP1L6+rNPW3WjfowmnVLef4j+nF768uP/2xmL2iz6h2SfkfSp1uPXyHph6NpA8D42KnHv/n//Jp5JTYx2cxWD3v8CXf/ROv3Q/7KZ2Y97l4rsT8YBwy0geOImc2V9El3v1DS+yR9tpUEslnSxe6+28z+Vs2BZEXSn7n7PjM7ZN3W/v5Z0jvc/ZGkybeoOWXjjcMSR14n6SmSXuTu/7O1zhfUnL99g7t/38x+KOk8M7tVkrW2UbLuyZL+t7v/XtD2yHULff6j4e6fMrP/JOmtknZL+kMzGz7Hu1/Sn0j6hJm9VVJN0rpW/wEcJdz95V1s/pC/8jHIPj6Yc7tRAAkz+ytJV7r77i72oUfSX7v7f+tiHy6V9DR3f3cJ++6X9JbgfyQATBBm9juSVrj7pa2/8v25u7+iG31Zs2ZNYQPDs846q6MvrhypJe/+eqGD2gc+8Moj7j9XtAG08/FuDrJbTNLfdLkPknSxmT02LOLviJnZiyT9b+XJMQAmhq8q/isfjnEMtAGk3P2ho6APByRlU1fGqw9XqTl/u+j93izpWUXvF8Cxxd0bYjrZcYnUEQAAAKAEDLQBAACAEjDQBgAAAErAQBsAAAAoAQNtAAAAoAQMtAEAAIASMNAGAAAASsBAGwAAACgBA20AAACgBAy0AQAAgBIw0AYAAABKwEAbAAAAKAEDbQAAAKAEDLQBAACAEjDQBgAAAErAQBsAAAAoAQNtAAAAoAQMtAEAAIASMNAGAAAASsBAGwAAACgBA20AAACgBAy0AQAAgBIw0AYAAABKwEAbAAAAKEHHA20z6zezfx5R+4CZXTqWDphZ1cyuMbOXB8ummNmXzey7ZvYNMzuxVV9hZj80s++Z2RvbrQsAAAB0Q1evaJvZUkk3SXpOsspbJd3j7v9B0uckvcfMeiV9WNJ/lPTrkt5kZidF65bdfwAAACBT+EDbzG4ws0nJsgfMbPKw0nRJb5D0nWR3L5J0fev36yS9VNKZkta6++Puvl/SzZJenKwLAAAAdEXPGLf7DTMbGPb4dEn/Q5Lc/T+OXNnMPiPpKZJOlnSDmR1w999097tay7N2Zkra3vp9p6RZI2pZ/WANAAAA6IqxDrS/7e6/d/CBmX2g3cru/rrWeg9I+o/uvm+U7eyQNKP1+wxJ20bUsvrBGgAAANAVR3vqyC2Szm/9/gpJ35X0U0lnmNnc1hSVF0v6XrIuAAAA0BVjvaKdMrMbJF3Qmj/9JO6+pJN9SPqYpM+a2c2S9ku62N0PmNk7JX1Tzf9R+LS7bzCzQ9Yt5AkBAAAAY2Du3u0+AAAA4CiwZs2awgaGZ511VvolvDIseffXCx3UPvCBVx5x/4/2qSMAAADAManwqSPHkoGBgaWVSuWdkl7baDRmVCqVnZI+32g0PtTf37+u2/0DcCjOWwDAsWJMA20z65f0loPJI2b2Gkl/Iel8d3+og/1UJH1U0nJJQ5Le4O5rR6zzRklvllST9D53X2Vm8yR9UdIUSYOSXufue6J1s7YHBgZeUalUrlm0aFHvwoULe6dMnaq9e/bMfGj9hrc9tH7wbSs+eL3u2Vwf7VMBMA7OnlfV2589VacuWKhTFy9SX1+fhoaGZg4ODr5tcHDwbXfccYcmTQpj/AGgrfGe5oCJ4YivaJvZf5b03yX9prs/2uHmr5Y02d2fb2bnSlop6VXD9n2ypLdLOkfSZEk3m9mNamZ2f9HdrzKzd0t6s5n9U7Suuw+NbLR1Reya5cuXT50165dx21OmTNGyM56qk+efqLfrLl1x0y5t2sscduBocOIU09ufPVXP/bXlGnneLl26VPPmzdPdd9+tarWqarXaxZ4CANB0RHO0zewSSX8s6aUHB9lmdrGZvSlZ/9IRmdtP3M3R3W9Tc5A83HMl3eLuQ+6+XdJaSc9UfBfIbN1DVCqVdy5atKh3+D/Ww82aNUunLFqoV5zOlTHgaHH+6ZN06uKFanfeLliwQPv2jTamHwCAch3JQPs/SHqTpLkadmXc3b/o7p8YvqKZ/XrrTpLvlnSxmQ20Bukj7/JYN7PhV9k7uQtktm7ktQsXLuxt9+ROO2WRXrCo7SoAxtELFvXq1MWL2q6zaNEi7d9/SLIoAIyZmT1vxN2wD9ZXmNkPzex7ramrMrOqmf0fM7vFzFab2QXj3mEcVY5k6shGSedJeoOkz5vZK9y9Ea3o7v8uqd/MLpX0NHd/tySZ2a/qyXd5rLh7bdjjw90Fcm9QG7nuIRqNxoy+vr62T66vr099xrQR4GgxyVyjOW/rdb5bAaAYZvYnki6RtHtEvVfShyU9p7XsFjP7mpo3zOt19xea2SJJF41zl3GUOZIr2mvdfZ+7/72aN4j5szHs44m7ObbmaN8zYvkPJP0HM5tsZrMknSlpjeK7QGbrHqJSqewcGjpk6vaTDA0Nacj5XgRwtNjvptGct8zPBlCgdZIuDOpnqjkOerx1g76b1bxT9cskbTCzr0v6B0nXjltPcVQqKt7vMkk/at2VcYGk6SOnj0iSu181ovRVSeeZ2a2STNLrJKl158e17v41M/tbNQfSFUl/5u77zOx9at4F8o2SNqt5x8jd0bpJfz8/ODj4xqVLl6ZzQx58eINu3XBg1C8AgHLduuGAlqzfoGVnPDVdZ8OGDaSOACiMu3/ZzJYEi7LpqvMkPVXNu1u/WNJnWv/FBDWmgba7D0gaGPZ4k6TFY9hPQ9JbgvqHhv3+D2r+X+Hw5Y9Kenmw3SHrRhqNxoc2bNhw6bx588IvRG7fvl0PbxjUdfcx1xM4Wnzjvv160eJBnTz/xPALkdu3b9fGjRs1c+bMLvQOwASTTVfdImmVN2+7/e9m9itd6Ftp6vW69u3bp/3796ter6tarWrSpEmaPHkyf01MTMgb1vT3968bGBh4zV133fVEjnZfva6hvXs1ODh4YMOGDQd6rfGaH/75+dd1u68AfmlgYOAVTzpvmznaT5y3jUbjNcuXL+e8BVC2n0o6w8zmStql5lXrD0o6Wc2prV82s+WSRn1vkaPd/v37tXv3bi1cuFALFy7UsM9fDQ4Oatq0afxFMWDN/+mamFp52n8s6ZJGozG9UqnskvSPjUbjw9xhDjg6cd4CGE+tqSP/7O7nmtnFak2PNbMVat7XoyLp0+7+ETPrk/QxSU9Xc0rsW939jm71fSzWrFlzyMCwXq9rx44deuYzn5n+RfHuu+/WzJkzn3Rle7xvArTk3V8vdFD7wAdeecT9n9ADbQAAAPxSNNDevXu35s2bp6VLl6bbrV27Vlu2bNG0adOeqDHQPoLUETPrN7N/PpLGzWyKmX3ZzL5rZt8wsxODdf7czH5gZrea2XNbtaea2c2t7T7WupV7uC4AAADGbv/+/Vq4cGHbdbiPQazbc7TfKuked/8LM/s9Se+RdPnBhWb2a5J+XdLzJJ0i6ctqZlZ+SNJ73H3AzD4u6VVm9mCyLgAAAEYhugo9MDDQ6Ovra3t1t3Ufg8ZZZ53VtW9FFnEFumhHdAv2iJm908x+K1n2F2Y2PGUkupW6Riy/wZsektTTuur9bEn/PmK7bF0AAACM0WjvP9L6zgyGKXyg7e4fcvevDa+Z2X9q3b70UknvbN2C/TzFt1IfLsupNP/l5PKx3IIdAAAAo/P5wcHBtjcXaS3/x3HqzzGj8IF2xN2/5O79kq6S9CF373f3G/XkHMrolulZTmUjqI36FuwAAAAYndb9Rw5s3749XL59+3a1IlY/PM5dO+qNy0C7jehW6iOXv8zMKmZ2qqSKu29W8y6U/SO2y9YFAADAGPX3969rNBqvueuuu/asW7fuwN69e9VoNLR3716tW7fuwF133bWn0Wi8hojVQ4053q810P2KpPuGlS9W87aja0dOH0n2MVXSZ9W8bft+NW+l/oiZ/S9J17j7D8zsL9QcTFck/bG739y609I/SJqkZmj8G929Hq07picHAACAJ+E+Bp0jRxsAAAAoQbenjgAAAKBEZvZKM7vbzH5uZleb2cyxrNNhmyea2XVm9hMzW2NmLziStsxspZk9ZGZ3tn6+FKzzPDNbbWY/NbNvmdmCMfR7qpl9sbWPn5vZq1v1lw5r+87WMjezZ7fdH1e0AQAAjk+tqOMfS3qhu//CzP5a0gx3f1sn64yh3X+RdKe7/5WZPUvS1yWdIWnaWNoys+9J+m/ufmuyfJKkdZJ+z91vMbO3Slrh7udH67dp539Jmu3ub2p95+82Sc919/Uj1rtG0i/c/U/b7o+BNgAAwPHJzH5fze/AvbL1eImku9QcTPrh1pHUK+mv1bwpYFXSjyS93d13DGvjKkkD7n5V63GPmmlwp7n7plZtQNLfSppyuP4Ez6FPzQjnb0haKmmtmt/Fe2jYOi+U9Al3f0br8SQ1o54XuvsWM/szSb+j5myOByS9zd0Hg7Z+0erfD4c9t7vd/UPD1nmtpHdIOtfda1GfD2LqCAAAwPHrFEkPD3u8Xs17j8wY5TrvllST9Gx3Xy5pUNIHDtPmPDXT3zaN2OfiUfZnpIWSvi3pTyU9S82rzP9mZsPvBPmk/br7fkmbJC0ysz+QdLaaV6afpeaA/ZNJW1H/Fh980BrAv1/SOw43yJa6fwt2AAAAlCe7qFof5ToXqHll+7zWuHaSpMckycy+L6lP0qmSfsPM3qFm3PKVbfY3mv48ibvfr1/GQcvMPijpvZKWSLp/lM/huZJWt55DVdLUZP1oP8P79hpJ60abbMdAGwAA4Pj1kKTnDXu8SNLj7r57NOuYWVXS5e5+nSSZ2XRJkyXJ3Z/Xql2lQ6eOyMzmuPvjw/a5Xs0pIIfrz5OY2TMlLXf34XeeNEnD71b5kJpx0Qe36VXzyvoGNQfWf+3uH2st65M0x8wWqnl1+6Dzh+3nkWH9u3PYOv9J0meyvo7E1BEAAIDj1w2SzjWzM1qP3yLp3zpY55uS/sjMJplZRc37mLy/XYOtKRVfl/Rm6YmB8tMlDYyyPyM1JP2tmT2l9fitas6bHv4Fxe9LOmFYusllkr7n7ttaz+ENw9JN/qekf3T3QXd/1rCfwVZf3tTq92JJL5e0qvXYJL1Y0rcO098n8GVIAACA45iZna/m4PhgMscfSDpd0idbc5bDddx9q5lNkfRBSf1qXhm+U9Kbhn8ZMmnzJDXnQT9FkquZGHJDu7YOs7/XqjlfvKrmlfHXqzl3/BuSznf3QTN7rqS/VzPZZEtrvw+0/gfhzyVd1OrLQ5Le4O4bgnamS/qYpF9rtfU+d/98a9mJkgbdvbddX5+0PwbaAAAAQPGYOgIAAACUgIE2AAAAUAIG2gAAAEAJGGgDAAAAJWCgDQAAAJSAgTYAAABQAgbaAAAAQAkYaAMAAAAl+P8BBGZY3a4u3SIAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Set projection programatically\n", "Viewer.setProjection(0)\n", "Viewer.setPlane(31)\n", "Viewer.ax.get_figure()\n" ] }, { "cell_type": "markdown", "id": "2e73703f", "metadata": {}, "source": [ "Lastly, if one presses the 1 key, the data are plotted with the x-axis being (1,0,0) and the y-axis being the energy" ] }, { "cell_type": "code", "execution_count": 6, "id": "b48d1c04", "metadata": { "execution": { "iopub.execute_input": "2023-07-25T11:38:27.431222Z", "iopub.status.busy": "2023-07-25T11:38:27.431222Z", "iopub.status.idle": "2023-07-25T11:38:27.734128Z", "shell.execute_reply": "2023-07-25T11:38:27.733118Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtoAAAG8CAYAAADgqwLpAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAsTAAALEwEAmpwYAABT3ElEQVR4nO3deZxcV3nn/+9TVa3WLlmWvEgyVhCOWGxEwLFNIKRJcAhGCsTBDD8PDsYhxBAmJuT1A+KYLDN4gASTCROWMZshCRNiswQZDDZLBxtjQHjBYjOWV7VkLMnW0lq61XWf+aNKdrlV57mt6nu7evm8X69+qauee889detW9dGtU99r7i4AAAAAxap0uwMAAADAdMRAGwAAACgBA20AAACgBAy0AQAAgBIw0AYAAABKwEAbAAAAKAEDbQAAgFHM7Ewz6+9gvQvNrL/5c4uZHTSzxcX3EFOBkaMNAADwODN7q6QLJO1z97PG0c4HJN3h7lcW1jlMKZzRBgAAeKLNks49fMPMTjOzbzbPUn/WzBblNWBmp0t6BoPsmY2BNgAAQAt3/6ykQy13fUTSn7h7n6QvS3qrmf2OmW0a9fOylnUulfS3E9drTEa1bncAAABgknuapA+amST1SPq5u39F0lfaLdyck73G3b85YT3EpMRAGwAAIPYzSX/g7g+Y2fMknZiz/Askfb38bmGyY6ANAAAQe4OkT5lZTZJL+sOc5ddIuqf0XmHSI3UEAABgEjCzWyXtad68191f283+YPw4ow0AANBlZjZbjROgfd3uC4pD6ggAAED3rZU018yuN7NvmFnH+d2YPJg6AgAAZrwXv3Ce73ykXlr7P/jh0I8kHWy568rWjG0zO03SWZI+KukUSdepkVwyUlqnUDqmjgAAgBlvxyN1fferK0trv+fEzQfd/fRgkbsk3e2NM6B3mdlONdJNHiytUygdA20AAAC56p51swMXSTpN0hvNbLmkhZK2dbNDGD8G2gAAAN33MUlXmdlNakQIXsS0kamPgTYAAJjxXFKm7n1vzd2HJZ3ftQ6gFAy0AQAAJGXq6tQRTEMMtAEAwIznctVJYkPByNEGAAAASsAZbQAAAHV3jjamJwbaAABgxnNJdQbaKNiUHWgvXbrUV61a1e1uTBpbtmzRypXlBe0XZdmDv5ou9vSkayM5CUezZiVLw08+Ll51YFe6aJbu0nELw3arew4ka/UFs5O12iP7wna1ZzBZ2vZ7T07WPOfVPntn+g/MnlXp/TBvW/yHKaum1/3jt3whWfvsa34rbLe6O71/9z7t2GRtwU92hu1q6FC6Njt9nKkef4lqaOXiZG3Wzv3Jmm1/NGw3W740WTu0sDdZ6713R9ju0Op0u7O2p/urQ/GV9Xx2+kCsz0vv39qjwTYlqZKeCWl708dK9ouH42aPT79/+K7d6RVz5vravLnp4uz087Z98TfDdqeaH/zgBzvcfVm3+wGUZcoOtFetWqWNGzd2uxuTxpo1a6bE/njJ8W/odhcAYMqaCu/zR8PM7u92H1oxdQRFm7IDbQAAgKK4ROoICsdAGwAAQCJFG4Uj3m+aWLduXbe7AAAAgBac0Z4m1q9f3+0uAAAwZbmc1BEUjoE2AACAS3XG2SgYA21MqPqjURxWenacVathu5VjlyRrs7btCdfNtv0i3e7iRcla9UdxTJz1BlFle9LxXT4Yx/vZ3HQs2In/cV+y9rP3HB+2O+uGdORgZTi93tztcZzbwcXp5+5rO5+erHk1Z2bboXTk44Lv3JusZY/uCpuNjiXtT0fM+XAQCyip98DBZK3+UDpirjJ/Xthu5YGHkrVZwTazevy8zdqTft2YBc/NrCCmU5IFUZ0W9Ml6068ZKed1E8TlRRGekuQ7H0kXo0jB4PUvSSNbt6XXrcX7EOVwMUcbxWOONgAAAFACzmgDAADIVFf8CQdwtBhoAwCAGc8lZczRRsGYOgIAAACUgDPaAAAAElNHUDgG2gAAYMZzMdBG8Rhoo1BnV84L65VZ6cgrmzUnWfMgyk2SFMSC+facGL6gTxoKcu1y4ueyE5eliz+/L1mqLJgftuu7g7jCnvRLunZfev9KUmUkPTlx3i/Std4dQ2G7I7PT2/3+7auTtVUL4vi56t50ZFsUIVd50sqw3eh5ze59IN1uFAsoSUHcW2Vh+jm3IGJSUvicayAd/VdZfkLYbLYliJ879pj0itHrSZKyIMYziPAb2TIQNluZExzf9fQ+yoLIxjyVX0nHU9oD6djQ3HZzIh1RnswZaKNYzNEGAAAASsAZbQAAMOMxdQRlYKANAABmPJepzgf9KBgDbQAAADFHG8Xjv24AAABACTijDQAAZjzmaKMMDLQxoSrHLE7Wsr2D6RWf/uS44SBKy5+8Iu7Tjt1BsZos1ZcHEWeSqrvTsWF7znlmer3hdPyZJFWG01F7c+/anqwdd0Y6rk2Sdu88MVlb8tNDydq+lXFs4I5npj84+/sX/0uy9rEPvyRsV57eDz5vbrJme/fF7c5LPx4LIuSy5UvDZrMfbErWaqtOTtZGjovj/Wr3BzFy1fTxO7JsYdzuwXRso+9LH9sWPC+SpJ6edO1Q+jirHXts2OzIznSMpz9nTbrdO+PXW313+v2hsj+I/5wVPM4cFhy/eRGqN2RXd7xdmOrOB/0oFkcUAAAAUALOaAMAgBnPJWWcf0TBOKIAAADUmKNd1s90Y2anm9k13e7HZMcZbQAAMOO5M0f7aLj7Rkmv6HY/JjuOKAAAADzGzPrM7Dtm9lkz+6mZ3Wpm683sBjN7wMz+obnMpubyV5nZ+83sm2Z2t5lda2bzJ2m/K2b2j2b2XTP7sZn9xMye17z/62b2d822XmRmW8zs+PH0qbQz2mZWlfQRSWvUmPp0sbtvGrXMrZL2NG/eK+kPJX1Q0lpJQ5Je5+53l9VHAACAw7JpOMVjHH5V0q+6+21mdp2kv5DUJ2mhpK2Svj9q+edI+k1JmaTvSjpP0icmrLePy+v3NZKWS3quu2dm9nZJb3f39Wb2akm3mdm3Jf2TpPPdPYh2ylfm1JH1kuTuzzOzPkmXS3rZ4aKZzZZk7t7Xct+5kma7+3PN7CxJV7Sug8nPanGkVRgLFkSR+Y/i/295ELtWeTB+jWSDOXFvCbV98XojOx5J1hZuS0e2mcUfNPlwOlLMg2i1h7+bjhSUpOU/SrdbORRFCu4I260OLUvW/sczX5qsrRgKotMkeS19vNjOR9Mrzk1Hp0mS9h1IbzPY95Xtu8JmKyuDmMkD6W1W7ng4bLce9Ml6e5O16k/vD9uNjiUPYvg0MhK2G/UpiuirLl4ctxu999x0e7LkvbM7bjd7cGu63aF0PKIkVWbNSre7PX5NoRyNHG0+6G9xr7vf1vx9s6Td7j4saYeZ7ZG0ZNTyX3H3IUkyszvb1CdKXr93S7pM0h+b2Wo1BuF7Jcndt5nZH0n6D0l/7e7fGm9nSjui3P0Lkl7fvHmypF2jFlkraa6ZXW9m32gOrJ8v6SvN9W+RdHpZ/QMAAEDS6P8tBv/DliS1njFwqWsfD+T1+7ckfan5+39I+rCe2NdnSPqFpDOK6EypX4Z09xEz+6Sk39ORE+b3S3qvpI9KOkXSdZK+qcb/NA6rm1nN3Y84RbJlyxatWfP4RQjWrVun9evXF/wIpo7BwUH19/d3uxsAgC6aCn8HNmzYoGuvvfbwzfhKTxOKL0POEOslbXD3DzVnV7xNUlWSzOwMSZeocaL338zsEnf/x/FsrPTUEXd/jZm9TdJ3zezp7n748/a7JN3tjc8n7zKznWo80AUtq1faDbIlaeXKldq4cWOpfZ9K+vv71dfX1+1u6H/qym53AQBmrMnwdyBPX1+frrjiCkmSmU2aeTLkaM8Yb5b0KTP7oaS6pG9J+n0zWyTp/0r6b+4+YGYXSvqemX2rZSrKUSvzy5AXSFrp7u9S4+x11vw57CJJp0l6o5ktV2OS+lcknSPp35tTSe4sq38AAACt6s6XISXJ3fslndpy+02j6oc/ifin5u0LR9WfcHuiHEW/nz1q1Uua/65uWXazpGPH26cyz2h/TtInzOxbknrU+B/EHDP7V3c/V9LHJF1lZjep8R/JiyTdIumFZnazGvNlXlti/wAAAIDSlDbQbk4ReWWb0rnN+rCk89vULy6rTwAAAO24jNQRFI4rQwIAAEjK+DLktNXf37+6Uqm8RdKrsyxbUKlU9kr6lyzL3tfX17e5rO0y0EahrCc+pOqDg+l1oxztLJ3pK0nV3ixZG8nJpK0EGdyRLMgElyR5uk8aTqck1Q/G+bvRfoq2mZ0S534/9Xd+lqz99LLTwnUjI3PSf7hm9wSZy7v3pGuS6g9vT9aqi9I55Z6Tm57tz3leU+0ObIsXCPLRK7PT2dJZkLEtSZUgBzp6LBYcg5JUmRfkjQdZ2VnOazXKDI/Ud+3qaL3x8pH0fopque3W68laJcp6D3LTMT7kaE9f/f39L6lUKtesWLGiZ/ny5T29vb0aGhpauHXr1j8aGBi4sL+//xV9fX3XlbFtjigAAABMS80z2desXbt27urVq3vmzJmjSqWiOXPmaPXq1T1r166dW6lUrunv71+d39rRY6ANAABmPJep7uX9oDsqlcpbVqxY0bMo8UnnokWLtGLFip5KpfJnpWy/jEYBAACmmkyV0n7QNa9evnx5T7RAs35BGRtnjjYAAAAkSZs2bcr5osXYnXrqqRN6Kn/V2790RN8/8eJ56u1Nfw9Gknp7e1WvZwtHr3/fu1867v4z0AYAADOeu7gE+zQ07KahoSHNCYIPhoaGNFTS9B6OKAAAAJmyEn/QHTcPHNIDWwbCZe5/cEA3D3SeIhThjDYKlRtFNmtWsmZBLYoFlOIYsyj+rLHh9Btg1Kc8UXxaFkT4VRbOD9uN4umyIPqrcte8sN3v3DT6irSPWxTE8I0ctzBsd9sL0p9C9ty2LFk7tieO4autWJ6seXAcZrvj+L7oGFUQrRjGLip+ziO1Y+MrAI/s3JmsVeenj6W811S4D7sQMZf7Og5kQwdL2W5lyeJkbWTbQx1vM+9YQjlcnNGejr58z7Cev3KrTjhumdp9IXL37t16cGCrrrunnPc1BtoAAACYlrYfcL3/B/v1p7pDJ61YrpNPWqFmjrbuf3BADw5s1ft/sF/bDxQ2Nf0JGGgDAACIC9ZMV3fuqOvSbw3qJU++T7+2ZUC95hpy080Dh3TdPcOlDbIlBtoAAABymTLyrqet7Qdcn/rRkD71o86m8XWKgTYAAIA4o43icUQBAAAAJeCMNgAAmPFcUkbqCArGQBtH7ezKeclaFLkmSdkjjyZrUdxYdcGCuFM96UO5HmxTyokcnB9E4tWzsN3q8enoOt+9J71iFn8poxLsiyyIelt0T9zujl9J1wdPTseNrfjP+IpbCzan1933nCBC7rjFYbt2MB05aNX0H8vKvjiCcjxRcJHq4sXJWnic5UTpVebOTbe7KB29WM2LkKvX0+sGF37wQ51n0Vql80GOB/2NVKJ9L0m19H6KIvzCmEjFEYn1vXvjPqEkpjp51ygY/3UDAAAASsAZbQAAMOMxdQRlYKANAAAgMXUEhWOgDQAAZjx344w2CscRBQAAAJSAM9oAAACS6pzRRsEYaKNQ2fYdcf1QFMmWjtHy4TgyzPfvT9YqvbPjdUfSbdd3PpKsWW8ca6cooqsWvPRy4tzqB9Lxc1HUW8++ON6vNpj+A1Ofl44yrO2LY9We/cqfJmvfuuspyVplYHvYbja4L10Ljpe8GLjoOIxXjP9AZ9HxENSi51SSfCh9OeGRga3hup2Kouvy9m8YTxk8p7nt5sTppYwEkZhSTnxirSdZ85yYzk7lRp2iYy4p6/IcbTM7TtIPJJ3t7uk3T0wZDLQBAABkXT2jbWY9kv6PpDjsH1MKn5EAAAB033slfVhSOR9FoSsYaAMAgBmvkaNtpf1IWmpmG1t+Xn9422Z2oaTt7v7VLj18lISpIwAAAJLq5Z5/3OHupydqF0lyM3uRpGdJ+pSZ/a67P1Rmh1A+BtoAAABd5O4vOPy7mfVLuphB9vRQ2kDbzKqSPiJpjRqfyFzs7pta6j2SPi5plaReSe909y+a2a2S9jQXu9fdX1tWHwEAACTJ9dgUD6AwZZ7RXi9J7v48M+uTdLmkl7XUXy1pp7tfYGZLJN1uZtdLMnfvK7FfAAAAR8gmwVfXGANNL6UNtN39C2Z2bfPmyZJ2jVrkaknXNH83SSOS1kqa2xxw1yRd6u63lNVHdKa2ckWylm2PM2nl6TxmmzUn3e7BdFawJFXnpbNuPcjuluI87Ci/O8ovluKMXZud3ma2ZzBsN9qHUX/nDaTztyXpwLHpfbj/Rel846Fj4pznu//+6ek+/VLwR83iM0tZlJ0e5TyPpPefJHkc15xUnZ8+fvN4lJ3ucR6zRY/1QDolrOO88Bx57UZ54nlZ2aHgdTEe0XE2HmEGd5Dtn7ePzq6cl6zdkF2d37EZzF2qc0YbBSt1jra7j5jZJyX9nqRXjKoNSpKZLVBjwH2ZpP1qxNt8VNIpkq4zszXufsRIacuWLVqzZs1jt9etW6f169eX9VAmvcHBQfX393e7GwCASWqy/I3YsGGDrr328Hk4Le1mX4Cylf5lSHd/jZm9TdJ3zezp7v7YqTEzO0nS5yV90N0/bWa9ku52d5d0l5ntlHSipAdHt7ty5Upt3Lix7O5PGf39/err65uQbb1HX5iQ7QAAijNRfyPy9PX16YorrpAkmVl8OeEJxhxtFK3ML0NeIGmlu79LjTPVWfPncP14SddLepO7f71590WSTpP0RjNbLmmhpG1l9REAAEA6/GXI7s/RxvRS5hntz0n6hJl9S1KPpDdLmmNm/+ru50q6VNIxkt5hZu9orvMySR8ys5vUSCq5qN20EQAAgKLVxRltFKvML0Puk/TKNqVzm/VLJF3Spn5+WX0CAAAAJgoXrAEAADPe4UuwA0VioI2j5oPpqLdsKI6QiyKtsiCKrDInjk6rD+ZE4pWguuSYsO4H0vvCg7hC64lflpU5s+OOJQzNjdutHUzHyA0NpuMIa0Nx/Nyjr9+TrPm30/swOyEOI6juS8euRZGOtWVxuyPbO/tuVt4xWJ0/P1209LzQvHajKMO8100kej1GEXPVBQvihoMIykqwbl6cZmVuOmbSgnXz4j8jee93kcr8eclafdeu9DZLihuEJOZoowQcUQAAAEAJOKMNAAAgKePLkCgYA20AADDjcWVIlIGBNgAAgMQcbRSOIwoAAAAoAWe0AQDAjNe4MiRTR1AsBto4eiPpOKwoakySVOvwkMviCLnqokXJmh9Kx4nlqRy3LN2lrdvCdS2I7woF8WdSHK0WRX/17B2ON7sg/dx9ru8Dydobv9ruulOP27cvHUdYXx0cSwfj/mbDQb1aTa+3Jycurzfd38qSxcma743bjWL6xhPDV1mcPvazoE9RfJ8Uv6asN32sjDy8PWw3iuHz4DmtHBfHMmowiL2z9OApitmTpJGdO+PtptoNHqckWXCMRmrHHhvWO+0vGvgyJIrG1BEAAACgBJzRBgAAMx5XhkQZGGgDAACI1BEUj4E2AACA82VIFI//ugEAAAAl4Iw2AACY8VykjqB4DLRRKM+J4YvewqK4Kx/Jibzbl468k2fhurL0Bzsj990frxuoHAxeXlFE4oIFYbvZrl3pdYN4RfvZA2G7e1/1tGTt33edkaztOyH+YOz3n3ZbsvaZjb+aXnHHo2G7UbxfdckxyZofOBi3G8TeZdseCteNVIPntb53b8ftZrt2p2tRBGKO+u50u9X58ztuNzr2I9nDO8K6zQliJMexf6O4x2woOJbyHme1sw+U8+L7ak9a2VG7aGDqCIrGQBsAAMx4pI6gDMzRBgAAAErAGW0AAABxRhvFY6ANAABmPBfxfigeU0cAAACAEnBGGwAAQMT7oXgMtHHU6vv2J2tRRJ8kZUG0WiWI5/KhIL5PkiyIFcyJ96v09qSbrXT+phs91tDgvrCct49TPCfq7ZRPpePcbvjlNclaJU5e1K/OvydZu2bvc9Mr5uz7KGIu2zOYrOVFRUZxj9GxFEUrSp1HzOW1G0X4dXqsSJIF260Ppvdv3jajCFDrSf9J8kNxXJ4P70nWon2YF4EYHS+VuXPT6w0Nhe3Wdz4S1juWE7GKgDNHG8VjoA0AAGY84v1QBuZoAwAAACXgjDYAAIA4o43iMdAGAAAzHvF+KAMDbQAAAEnOQBsFY442AAAAUILSzmibWVXSRyStUePLvBe7+6aWekXSByWtlTQk6XWS7hl9n7vfXVYfAQAADiNHG0Urc+rIekly9+eZWZ+kyyW9rKX+ckmz3f25ZnaWpCskfbLNfa3rYIKcXTkvWbNaOnc6L6M4zD4+cCC/Ywm1Jccka/Vdu+KVg2zk8WQUh9m9Q0HGdk7udyTKKK4uXhSuO7R4TrL2/V/5eLK29utvCNt96qxfJGteTfc3LzfZ5qb7a8G6ucfoooXp4ki63ShbWoqPl/BYqsVv01bPyZhP8Jz1vMPXY167Uof9zXneKnPSx8N4sqWjPPFsf/qaAt0ysmWg212YspwcbZSgtIG2u3/BzK5t3jxZ0q5Rizxf0leay95iZqdL2tzmPgAAAGDKKfXLkO4+YmaflPR7kl4xqrxQUuul6Ort7jOzmrsfcRppy5YtWrPm8avVrVu3TuvXry+s71PN4OCg+vv7u90NAMAkNVn+RmzYsEHXXnv4PJyWdrMvo/FlSBSt9NQRd3+Nmb1N0nfN7Onufvj60nskLWhZtNLuvnaDbElauXKlNm7cWEqfp6L+/n719fUV1t7l+kBhbQEAuq/IvxHj0dfXpyuuuEKSZGY7utydFsT7oXilpY6Y2QVm9hfNm/slZc2fw74t6ZzmsmdJujNxHwAAQOncrbQfzExlntH+nKRPmNm3JPVIerOkOWb2r+5+rqTPSzrbzG6WZJJeK+muNvcBAAAAU06ZX4bcJ+mVbUrnNuuZpIvb1NvdBwAAUBoXqSMoHleGRFu1VSeni0H8nM8PIrYk1X90V7IWxeFFkYJSHA2Yt67N7k3WasvS39PJHnk0bLcyf1661pOOIxzZui1ud+7cdC2KiQv2ryTVe9MzyV7805eG60ae0ZPu70lfSx9LNmd23HAQtRfGJ+bIjYPsUBR7Vz3+uPSK89L7T5L84fQUV1u0IFkbeWBL2G51QXrd+t694bphu4vSMZP13buTtTzh8RJFReYcK1G7fvpT0+vdeFvYbqTy7Gcka9mtP+q43b7feU+y1v+Vt3Xc7rThjYg/oEgMtAEAAMQFa1A8LsEOAAAAlIAz2gAAYMZzkaON4jHQBgAAIEcbJWDqCAAAAFACzmgDAACI1BEUj4E22vJd6Zgt60kfNlHUmJQT4TcniAasxh+++IEgoitnXasFL4Ng3ShSsLFA8I5dGcfHk8FfgjDiLGc/zL31vmTtmjVfStZOHYij71/7wK8na7tel46J691+fNhu7f5fJGuV4FiynJhDP3QoWcv270/WqosXh+1GokhHO+O0cF3/eRCJN464vPFE+IXtBn0aT/RfPSdus1NRu7Xb0sdKOswxn911/zjWTuu5fmMp7U4nzNFG0RhoAwCAGc+dgTaKxxxtAAAAoASc0QYAABCXYEfxGGgDAACIL0OieAy0AQAAxBxtFI852gAAAEAJOKONtqK4vPquIEovR3XJMel2S4rnsmo1rNeHhtLrBhFnUVShJPlIOvor6lOlN4joUxwrGNVqeZ+JWvpMzhzrTdY+8Gv/Gjbbv+dpydrILacma5Xbbg3bzTxL1qw33d+8mLhI9JzXd+2K1507t7Nt/vjesD6eGLluqC5YkKx5cPxGkY1SfOxHr6lsKH4/i/pbVgRidjD9noTyuIwz2igcA20AAABJTNFG0Zg6AgAAAJSAM9oAAABcsAYlYKANAAAgMXcEhWOgDQAAIM5oo3jM0QYAAOgyM6ua2cfN7NtmdpOZpWOZMGVwRhsAAEBdvzLk+kYf/Hlm1ifpckkv62aHMH4MtNGW1ztM57X4Q5Js9570qrWedC0nCzsvCzfiWfqd1YLNZsPDccPBvojazd330T4OsqWzffvjZk9emawd8pFkbVc9nY0uSe8+/o50u+eld8SP/1ecUy5Lv31F+zAv/zw7lH6s0b7Py8nOgmz6qN0oH1qKXzfhsRQcK3nt2qxgm0EuvZSzH4I+ec7rLcrZDjO2c/K5y8rKDrcbPG9MIS6Pq7tTR9z9C2Z2bfPmyZJ2da0zKAwDbQAAAJdU7kB7qZltbLl9pbtf+YQuuI+Y2Scl/Z6kV5TZGUwMBtoAAADl2+Hup+ct5O6vMbO3SfqumT3d3fdNQN9QEgbaAAAA6u4cbTO7QNJKd3+XpP2SsuYPpjAG2gAAAFK3J8F/TtInzOxbknokvdnd4y9pYNJjoA0AANBlzSkir+x2P1Cs0gbaZtYj6eOSVknqlfROd/9iS/1CSRc2b86W9CxJJ0j6hqTD0RT3uvtry+ojAABAg3HBGhSuzDPar5a0090vMLMlkm6X9NhA292vknSVJJnZB9QYlB+UZO7eV2K/MAY+cqij9awSv0l1GqWXF3kXRpH1xIe5B9FfUYRfbpzb/nScns2fl15xX+efFFo1HV2XF5Fow+nn/JKtz0vWPrjiu2G7D9UHk7VT5vwiWftxlo4blKRsKL2fqosXJ2v14HnJEx2HeX+eK3NmJ2vRseK5SZtBlGFvept58X7RazXqb7jNnO1WVp6UrI3cc1/YbF4MYqeq8+cna/XB9LE9nveH3H3YobwoQ6jbU0cwDZU50L5a0jXN301S23BaMztd0jPc/U/M7ExJc83s+mbfLnX3W0rsIwAAgORcgh3FK22g7e6DkmRmC9QYcF+WWPRSSX/b/H2/pPdK+qikUyRdZ2Zr3I+8WsaWLVu0Zs2ax26vW7dO69evL+4BTDGDg4Pq7+/vdjcAAFPQRP792LBhg6699vB1WbR0wjYMdEGpX4Y0s5MkfV7SB939023qiyWtcfdvNu+6S9Ld7u6S7jKznZJOlPTg6HVXrlypjRs3jr57xurv71dfX19h7V2uDxTWFgBgcivy78dYtnXFFVdIksxsx4RteCyYOoKClfllyOMlXS/pTe7+9cRiL5DUWrtI0mmS3mhmyyUtlLStrD4CAAA8jqkjKFalxLYvlXSMpHeYWX/zZ4WZfa5lmTWS7mm5/TFJi83sJkmfkXRRu2kjAAAAwGSXPKNtZjcH65mkzN2TUQTufomkS9qUzm1Z5u9HrTMs6fxguwAAAOVg6ggKFk0dMUmvCmr/t/juYKqrBFFYkqQsHe1V37u34+1WlxyTLh6KPxSpPSkdI+eP7ErWsoNDcZ8WLUoX60HE2bHBY5E0sjU9myqKZazOSkcgSpKCGMTLT/xmsibFMWYnVNPHxMWLBpK1zw/F35GKIvw00vkHYdHzZgvSj2VkS/qxSDn9HUfkYCQbOtjxurVl6f0/sj09rdbmxRFyUcxkXoRfGcYTCxjtI89534me83E9byeekKyNbHuo43ZnDAbaKFg00P6opIdTl/80s1eX0yUAAIAJ5pKI90PBojnap0m608w+aGZrRxfdfXN53QIAAACmtuRA293fLOlpalwS/X+a2bfN7A/NLP6cGAAAYApyL+8HM1OYOuLuh9z9Gnd/qaTz1LiIzAMT0jMAAICJ5CX+YEbKzdE2s9mSfk/SH0haIOmtZXcKAABgwjFHGwWL4v36JL1G0gslfUHS/+/umyakVwAAAMAUF53R/htJV0q62N3jHDNMOWdXzgvrteOPS9Z8aDhZq+/aFbZb6Z2drs2alaxlw+ltSlL9kUfDeiiIFazMTX8lobrs2LjdIGIuikfT7t1xux3Ki0+s3JuOBnzOhjcna/e87MpOuxSqnXxSWB+5/8FStluP9v84npu810ZKGF2p+NivHZs+Rkd27gzbDY/RDvsjSbUnr0oXg22G8YiK9+/B3z0jWZv9xe+F7UY63UfjEcaGSlJvb7JUW7kiWcv7m3BDdnW83WnCmOKBgiUH2u7eJ0lmVjWzP5R0shpfjNzk7hP/7gIAAFAW5lKjBGO5BPv/UWOQfbYac7Q/VWqPAAAAgGlgLAPt1e7+V5IOuPsGSTmfWwEAAEw11vgyZFk/mJFyU0ck1cxsqSSZ2QJJ6WtHAwAATFVMHUHBxjLQvkzStyWdKOkWSW8us0MAAABdwUAbBcsdaLv7f0paY2bL3H37BPQJAAAAmPLGcsGaP5b0x5JmmzXmGLn700vuF7ps5BcPd7Se1XrCuo+kI+QilTlz4u1Wq8lafXAwbjuI8Iuum1vfHsejRY81ijnMhg6G7Ub9jfZTlhPvp2r6KxvLvx7ML3xZ3Gynsofi/9dX589P1mxuej+MPNx5ux7ETOZFUEbxlbL0vh9PdGVehF+nooi5MB5Rkj+cDq2KXhfjuYb13BvuTBej50VSZcGCZK2s/RvJ279RBGUU74cmzmijYGOZOnKJpHMkjSOoGAAAYBJz8aVFFG4sA+0fSnrQ3etldwYAAKBbuGANijaWgfY3JN1jZpslmSR3998st1sAAADA1DaWgfYfS3qlpF3ldgUAAKCLOKONgo1loL1F0vfdnfxsAAAAYIzGMtDulXSHmW1S8/967n5+qb0CAAAAprixDLTfVXovAAAAuowvQ6JoyYG2mX3F3X+necGaZL28rqGbwszfILM6yhmWJOvtDdbtLGN7LNuNZPv3d7RelN0txZni1pP+P67V4yxy1dOzuKJc37yM82qQwf3w6emc5/HYPJLOOLd5cXa6H0jnjdeDrOzw2FZ+7nrHotfNoZGOmw0fTy19nOW93qLjO3q9RTnkUrx/o4z43PzoQLR/87L9sw6zsse1H4LnNC+vPTKyZaDjdWcM4v1QsOiM9m+a2dZEzSQdU0J/AAAAJp6LL0OicMmBtrvHp30AAAAAJI1ljjYAAMD0xxltFIyBNgAAgPgyJIp31ANtM+t196EyOgMAANA1DLRRsGSUgJl9puX3P28pXTeWhs2sx8z+2cxuNLPvmdnvtlnmVjPrb/58wswqZvZhM/tO876nHNWjAQAAACaJ6Iz2cS2/v1TSFc3fx5p982pJO939AjNbIul2SV88XDSz2ZLM3fta7jtX0mx3f66ZndXc5svGuD0chShGS5Lk6f/Wh1FkFsfAhZFiwcVH8+LP8iK6ImEM36x0LS8W0GrpfRFGe/XODtvNhtKxdlEsmGc5p2oq6Zf2cRuDC8P+Qdxs5P0PvzBZ830HwnW9Xk/WwijDnGO0umhRepsH4j5FLIpsG0e7lUULk7WR7TvS6+XEHJZ1nFWCGEnrST9v43ldjOf9IexvFFcaxE/mGU+EXySKHCwt1nKq4Yw2CjbWqSOtf4HHehheLemalvVHj5TWSpprZtc3+3GppOdL+ookufstZnb6GLcFAADQMXPmaKN40UDbE7+PibsPSpKZLVBjwH3ZqEX2S3qvpI9KOkWNKSnflNR6VYK6mdXc/YjTmVu2bNGaNWseu71u3TqtX7/+aLs5bQwODqq/v7/b3QAATENF/n3ZsGGDrr322sM3lxbWMDAJRQPtZ5jZp9U4G936+9PH2riZnSTp85I+6O6fHlW+S9Ld7u6S7jKznZKqkha0LFNpN8iWpJUrV2rjxo1j7cq019/fr76+vjEv/y59orzOAACmlaP5+zKWtq64ojEb1czS85u6gStDomDRQPuVLb9/OPF7kpkdL+l6SW9y96+3WeQiSadJeqOZLZe0UI1pI+dI+vfmHO07x7ItAACAcWPqCAoWXRnyP8fZ9qVqXKb9HWb2juZ9/1XS/3b3cyV9TNJVZnaTGof2RZJukfRCM7tZjbPnrx1nHwAAAMaEOdooWmkXrHH3SyRd0qZ0brM+LOn8NvWLy+oTAAAAMFG4MuQMlRdNF8Z3BVFZeZFhUcRZXoRfKIoGDGLgJMl60i+DKNor10j68VSCmC3PifbqNKKr9uRVYbvZlq3JWllnef5x+feTtRcPrS1lm3lRb2H4Xy19rOS+pqrVsN6pKMIvlNOfsqLgwnZ3707WxiN6HedFK4b1oBY9TklSSZd9i7Zrc4KIROL9GjijjYIx0AYAACDeDyVgoA0AACBxRhuFiy+RBgAAAKAjnNEGAACQOKONwjHQBgAAEHO0UTymjgAAAAAl4Iw22rIo1i6In8tyoulqixYma/XBdExZdfHisN3KrAVhPeJBDF8Y7eXxqY9oX9SOXZJudveesN1ItJ/84TgGrhL0afeqSfh/8l97Vrp28+3JUu24ZXG7xyxKlkZ+dneyVl1yTNhs/ZFH0+suCI7fSrzvO43Eq5x4fFgfuee+ZK127LHp9XbuDNvtNBqwturksD5y3/3JWl6EXxnGE4FYO/GEZG1k20Odb5cIP2DCMdAGAACQmKONwjHQBgAAIEcbJZiEnwcDAAAAUx9ntAEAACSmjqBwDLQBAAAkBtooHFNHAAAAgBJwRnuGiuK5pDiiq3bK6vSKj+zqsEex+q643eqidCRbXvxZ5TmnJmt2VzoyzOv1sN3aU09JF4NowMqiOKrQtz2crEX7qbb8xLDd7Ph0PN3yG/emV/zzsNmO3f3+s8L6Uy+/J1lLBzZK6ukJ2z2wKr0f5uxbkazVT0jHI0pS5eT0/q/f9uNkLS/aMlJbme6v70zHDUpxxFwWPdaceL9ORfF9Us7xPWd2ut3N98btBvtwZMtAuG6nogi/2pNWxus+sKXo7kiSzq6cl6zdkF1dyjYnmokvQ6J4DLQBAAAkpo6gcAy0AQAAiPdDCZijDQAAAJSAgTYAAIDUmDpS1k/AzHrM7J/N7EYz+56Z/W7RDw3dwdQRAAAAqZtztF8taae7X2BmSyTdLumLXesNCsNAGwAAQF2do321pGsOd0M54UmYOhhoAwAAlG+pmW1suX2lu18pSe4+KElmtkCNAfdlXegfSsBAe4aKcrLzZFu2pmsHDoTrWrWarFUXpPOj8zKrVU1/3aAyZ07cp81B7mwl3a4FNUnK7knn/lpvb7LmQ0Nxu8PDYT3p0KGwXBnYnqztef6TO9vmOKzakHNC55h0dnp1f/o4HBlIH7+S1BPUwx7lZCpXgqz3SN7xEMl2PpKu5bxWK8HxYkGfokz7xgLp1039kTjbO5IFueCV45ema3Pnxu0GfarOn5+s1QcHw3Zry9J9Gtm+I10rKSe79pT4NT5ydzq3flop94z2Dnc/PVU0s5MkfV7SB93906X2BBOGgTYAAMAYvrRYFjM7XtL1kt7k7l/vTi9QBlJHAAAAuutSScdIeoeZ9Td/4o9jMSVwRhsAAEDd+zKku18i6ZLubB1lYqANAAAgcQl2FI6BNgAAgLgEO4rHHG0AAACgBKWd0TazHkkfl7RKUq+kd7r7F/PqZnarpD3Nxe5199eW1ceZLC/yTlGcXpb+L38U35dXj2LMKovjyLD6I7vSRc/Cdf1QOrTNKpau5exDD/ZTtndvut1aT9huFEcW7d9scF/c7vx5ydqCn3Qeuxb584eenayNzI2PpezeB9O1oYPJWt6xnxd717HgNVXpnZ2sjac/0bqVWbPilT19/I4nhi+K/4uObR+O4ymj51zb03Gm49m/edGAkWxPHP+X3GZwrEg5+yGye0/+MjMBZ7RRsDKnjuRdTvSIupldL8ncva/EfgEAADxRF+P9MH2VOdDOu5xou/paSXObA+6apEvd/ZZ2jW/ZskVr1qx57Pa6deu0fv364no/xQwODqq/v7/b3QAAzDBH+7dnw4YNuvbaaw/fTF+5Z4JZ8wcoUmkD7bzLiSbq+yW9V9JHJZ0i6TozW+PuR3y2v3LlSm3cuHH03TNWf3+/+vr6xrz8u/Tx8joDAJgxjuZvz+Hlr7jiCkmSmaUvgwlMA6WmjuRdTnR03cx6Jd3t7i7pLjPbKelESenJmAAAAEVg6ggKVuaXIcPLiSbqF0k6TdIbzWy5pIWStpXVRwAAgMOI90PRyoz3a3c50RVm9rlUXdKnJS02s5skfUbSRe2mjQAAAACTXZlztFOXEz03p35+WX1CiyB6TpIsiHrzICausmBB2K7v35/eZhCVle1Ox+FJks1KR+JZTxyX54eC2LCR9P/zssHO4rmkOFrNcmLX6sF2q/Pnp9udkxMLti94bk4o5/tK3/jEmcnaiTf+JFzXg9jGKMLPap2/7UXRdPXdu8N1oxjJ8cQRRq9lDyIFs+HhuN0ORftIkjyI0xtPn6pLjknWojjCvJjDqE9Z8H6WJ3puOn0sUvx4oscysp2p0pKYOoLCcWVIAAAAiYE2CsdAGwAAwJmjjeJxCXYAAACgBJzRBgAAkJg6gsIx0AYAABBTR1A8BtoAAAASZ7RROAba09iLT/+bZK0SxPdJUhZE+MnSU/v9QDqmLI8PDQXFdJSbJPlwEGOWE8FVCWIFw6i9vCiyajW9zWMWJ2u+Px1/JklWS8cVRtF/ecLIwf3BczMOy24PHmtOBGUUVVbpTUcZ1vfGUZGRvAi/iM3uTReDeL8siMOT4uPXxxEbmLfdlPHso+rixel2d+0K1/WD6WM0ir2MYvbK5CPpWNG8CL9I9N6S/eLhjtutHX9cx+sCMxkDbQAAADF1BMVjoA0AAOBi6ggKR7wfAAAAUALOaAMAAEic0UbhGGgDAIAZz8QcbRSPgTYAAIDEGW0UjjnaAAAAQAk4oz2NZbf+KFmrHLcsXjfI362tOjm9Xk5Oa3YgJ3u6Q7Xg8eRle3eaq1x9xi/HC2zbniyNbN2WrNVOPCFsNvrfcV7WcCTKpbZ56VzqqwcXhe2eNz+dq+xVS9byHkuUla1Kut2yVBfF+6HTfOkoW1qK91Pt2GOTtZGdOzvqT5nGc/wOP+/pyVrtho0dtzvVjIwjK7sb7U425pzSRrEYaAMAABDvhxIw0AYAABBfhkTxmKMNAAAAlIAz2gAAABJTR1A4BtoAAABi6giKx9QRAAAAoASc0Z6hfP+Bjtfd/asnJmsLr+s8MqyyNh3Pteu0ODrtmA0/Ttby4vuqS45JFw+NJEu2azBs11ccn6zt+601ydq8q28J2+1U7ZTVYX3k55uTtWx2+q3iRXPTUYWS9PubX5bu056hZK2y+pfCdjWSfm5G7n8wvc2TT4qbDdaN5MX3jbzo9HSfvpaOn7M5QYyhJO0KtjkJI/yiuMLxxPt1GuFnZz4zrPt3f5he94zTkrXqnjhWdOSnP487NsEqz35GWI/iYs+unJes3ZBd3XGfuoIz2igYA20AAABn6giKx0AbAABA4ow2CsccbQAAAKAEnNEGAAAznompIygeA20AAABJckbaKBYDbQAAAHFGG8UrbaBtZj2SPi5plaReSe909y+21CuSPihpraQhSa+TdM/o+9z97rL6OJPVB+NouspzTk3Wovi5es52ayuWJ2v+YDombv4xccRZFK0WxapJkh9Mx8TVbk8ffn4gjkis37k1WVt04MnpFZ+0Mmw3ihwc2fZQsuZzZoXNVhcFEYoHDiVLZ9z4hrBdPTgnWVpdCY7DIL5P6jyG75Ffj/fvkmo1vc177utom1Ic4ReJnlMpfk2NDKSPwTy1NU9Jt/uzzt+Wowi/2soV6W1uGQjbrR23LL3uw9uTtSi+T5Jqxx6bbvd7d6ZrYauxKHK0/sij4bq15en4VdXSx7YPpV/jADpX5hntV0va6e4XmNkSSbdL+mJL/eWSZrv7c83sLElXSPpkm/vSIbwAAABFcJE6gsKVOdC+WtI1zd9NR/4H//mSviJJ7n6LmZ0uaXOb+wAAAEpnWbd7gOmmtIG2uw9KkpktUGPAfdmoRRZKav28v97uPjOrufsRn8Jt2bJFa9Y8fnW9devWaf369QX1fuoZHBxUf39/t7sBAMBj2v1d2rBhg6699trDN5dOZH+AiVbqlyHN7CRJn5f0QXf/9KjyHkkLWm5X2t3XbpAtSStXrtTGjZ3NeZyO+vv71dfX94T7LtcHutMZAACkI/4uHb7viiuukCSZ2Y4J7lKMqSMoWGkXrDGz4yVdL+lt7v7xNot8W9I5zWXPknRn4j4AAIDSmZf3g5mpzDPal0o6RtI7zOwdzfv+q6T/7e7nqnGm+2wzu1mNOdyvlXRXm/sAAADK5SJHG4Urc472JZIuaVM6t1nPJF3cpt7uPkyw7AebkrUo7kq9cYRcFInng/uStVk/eiBud/78ZK26MSeKLHpjDaLeGgmVwaqLFydrI3ffk14vitmT5IfSMVxRLJju3RK2a7N708V9B5Olhd9aEra75CfpdSsHh5M137UnbLe6YEG6GDynS75+b9iuH0j3N4qQ86H0Y5HiCMqw3eE4ds337Q/rHcvS+zA6tqP4PkmqnromWRvZ9LNkrXbiCWG7UQxibdXJ6fXuuz9ud+fOZK0yd26ylu3v/HnJi/CLjGxNx6RGz5seTR+fADrHBWsAAADEFA8Uj4E2AACAxJchUbjSvgwJAAAAzGSc0QYAADOeiakjKB4DbQAAAHdSR1A4BtoAAADijDaKx0B7GqstS1/ZNgui9CQpi2L4RtperLNRy4m0slnp+L/KMYvT/Xl0V9xubzqaLu+xyrNkqXr8cel2d8VxWFEsWy1od+QXD4fthoL9Xzv5pHDVbHsQY1ZLv1XM31oP263tTsfl2aFg3QXpyEYpp79L05GDIwNbw3YjUfSiHxyK1w3iCEce3p5eL4pkU36cXkp0DErSyM83d9RuJXiNS5LflY5XrPTOTtbyXm/hNoP3j9pTT4nXHUjHBtb37u20S13R6bEixbGuUQQiMNMx0AYAAJBIHUHhGGgDAACIqSMoHgNtAAAAV3hFVKAT5GgDAAAAJeCMNgAAgMQcbRSOgTYAAICYo43iMXUEAAAAKAFntKexke07krXq/DijuLrkmHQx+LJIlL8tSYrqQUxulK8rxXm2tVUnx32qpv+/ObI5nfkbZSpLUmX+vHRx3tz0enPmhO3m7uME3zvY0Xp55n//gbC+f+3KZG3uXcFjOZDO35Ykmx1kpwfHft6xFMn2pXPKq8fEx4OCLPJqcDW6vOzj6HiJjpW8vPYoDzsbHu6o1i313ek3l9rc+PXWaVZ2bcXysJ7teCRZq5ywLFkbuf/BeLtl5V0H75PTCleGRMEYaAMAAIipIyjeDPkvKgAAQMBL/hkDMzvTzPqLeUCYDDijDQAA0GVm9lZJF0ja1+2+oDic0QYAADOeSTL30n7GYLOkc8t9lJhonNEGAACQpKzU1pea2caW21e6+5WHb7j7Z81sVak9wIRjoA0AACCN9cxzp3a4++llbgCTDwPtaSyK/bKFC+KVg8ir7IEt6W3OTcfWSZLqwemCiiVLFkSNSZKGgii4nFi7KPKqtvqXkrXswYG4T1n6sdYf3p6s5e3D2rKlyVoU6Wg5sYE+GEwLrNfT6+XE8M39YXo/ZcctTtbsF+l9JEnWm473U3C85B77hw4lS9Fjjfa9FMcKZsHxW128OGw3L/6vU53G9N2QXa3+/n719fW1rZ9dOS+5bu3EE5K1kW0PddQfSaqdfFKydmDN8eG6PR1u97oH/zGsR/vBHtnV0TYlKdufjqAcj5HgPSsv6hSYyRhoAwAAHEU6CDBWDLQBAADkXb9gjbvfJ+msrnYChSJ1BAAAACgBZ7QBAADElSFRPAbaAAAAUtenjmD6YaANAADgkpWbo40ZqPSBtpmdKek97t436v4LJV3YvDlb0rMknSDpG5L2NO+/191fW3Yfp6so/izb+Wi88qO70u0GMXGeFwkWRfjV0odjfffuuNkgEi+LYusk1Y5bli7uP5AseRaf+bCe9OOpHX9cspYF+16SVEl/tSLaDyMDW8Nmw4iu2eljySz9nEoKzxBVBtKRYSODcSyjgnoUbek5+7dy7JJkrZ4T4ReJjofq7PS+96GhjrcZyYuRrASvi5H77k/WDsfWXa4PHHWfogi/vJhDWzg/3e79DyZrPUEtT23F8mQtiu+TpOr8dH/re/d23KfsQPo9qyzR+3PefgCmu1IH2mb2VkkXSDpipOPuV0m6qrncByR9XNJBSTZ6UA4AAFA6po6gYGWnjmyWdG60gJmdLukZzcuQrpU018yuN7NvmBkRNwAAYGJ4iT+YkUo9o+3unzWzVTmLXSrpb5u/75f0XkkflXSKpOvMbI27j4xeacuWLVqzZs1jt9etW6f169cX0u+paHBwUP39/d3uBgAARyN9qdsuKPkS7JiBuvplSDNbLGmNu3+zedddku52d5d0l5ntlHSipCMm0a1cuVIbN26csL5Odu0uefxu/Wt3OgMAwNh0/sUHYArodurICyR9veX2RZJOk/RGM1suaaGkbd3oGAAAmGE4o42CTeiVIc1siZl9ruWuNZLuabn9MUmLzewmSZ+RdFG7aSMAAACFcklZiT+YkUo/o+3u90k6q/n7I2r5cqS7//2oZYclnV92n2aKodOfkqzVvjaOaTdBfFTt2GPDVUd27kzWwni5HNn+/claXozZyMPpiLlIFNEnSerpSdeCyLYollGSPIgcjGLton3U2HA6pm/kgS3JWm1ZPMWyfvLx6drs9FtQz/ChuN1du5K1KOKs9tRTwnZHfvrz9LpRFOSc2XG744iRi0T7fySII8w7HrIgwm88aieflKxF+yh6viVJefUS5EVmRmxx8H6XF23ZoSiOcDyPBUBat6eOAAAAdJ3J+TIkCsdAGwAAQGKONgrHQBsAAEBioI3CTeiXIQEAAICZgjPaAAAAh1NHgAIx0AYAABBXhkTxGGhPY+OK8AtUT12TrI1s+lm4rj//Wcla/abbO+yRZGeclqxl37uz43Zrq38pWRvZfG/H7ZYl+41nJ2v7lsexdovuCiLFfrArWdrdl46RlKSFmx5J1jx4bn76v88K2z3lv92SrEXHwyO/PD9sd2EQ79dpFGSZogi/yajTmMMDv3dmWJ+9YzhZsxtv62ibefzXf6XjbY5sGSi6O7miCL+Dv3tGuO7sL36v6O5MTgy0UTDmaAMAAAAl4Iw2AACAnDPaKBwDbQAAABcDbRSOqSMAAABACTijDQAAIBHvh8Ix0AYAABDxfigeA20AAACJOdooHAPtGeqG7OpS2j27cl5Y/9q3/rKU7Y5Ff3+/+vr6jrg/6vN1P/+7jrd3ytXvTNZW/Zc7krXxPDfRY7llDO2m9tGk9CcltfupuDyl9lEg77Xa6XFYVrvTTVnHUfbQLydrlRPu6rjd6Hkdz3NqZh2vC0wFDLQBAABcUsYZbRSLgTYAAAA52igB8X4AAABACTijDQAAIHFGG4VjoA0AACAx0EbhGGgDAADwZUiUgIH2NNaNGK2pGN1VVp9/ft5l6WKcgNaxqbj/UZ7Vn7k8WdtcUozkDdnV0yYCcSoaT4RfhPcWoDMMtAEAAOSScw12FIuBNgAAgMQcbRSOgTYAAABztFECcrSniQ0bNnS7C5Me+ygf+ygf+ygf+ygf+wiYGRhoTxPXXnttt7sw6bGP8rGP8rGP8rGP8rGPJin38n4wIzF1BAAAQGJAjMJN2YH2D37wgx1mdn+3+zGJLDWzHd3uxCTHPsrHPso3LfaRvSqInxxPu2bSNNlHJWMfNZzc7Q4AZZqyA213X9btPgAAgOmCKR4o3pQdaAMAABTGJWXkaKNYDLQBAAAkzmijcKSOAAAAACXgjDYAAIDEGW0UjjPak5yZVczsw2b2HTPrN7OntFnm1mat38w+MZZ1ppO8x2tmF7bsn1vM7KCZLR6937rV/4lkZmeaWX+b+4/YhzPtODos2Ec9ZvbPZnajmX3PzH63eT/H0RNrM/r96LDgOOL9aNLyxpUhy/rBjMQZ7cnv5ZJmu/tzzewsSVdIetnhopnNlmTu3tdy37nROtPQyxU8Xne/StJVkmRmH5D0cUkHNWq/TXdm9lZJF0ja16b8ch25Dz/Z5r7pfBzl7aNXS9rp7heY2RJJt5vZ9eI4aq3xfqR4H/F+NIm55M6XIVEszmhPfs+X9BVJcvdbJJ0+qr5W0lwzu97MvtH8Q5a3znQzpsdrZqdLeoa7X6n2+2262yzp3ESt3T6caceRFO+jqyW9o/m7SRoRx9FovB81RPtIEu9HwEzBQHvyWyhpd8vtupm1fhKxX9J7Jb1Y0sWS/nUM60w3Y328l0r62+bvR+y3ab6P5O6flXQoUT5iH7a7bybvI3cfdPe9ZrZA0jWSLhPH0Wi8Hyl3Hx02o9+PJi2mjqBgvJAnvz2SFrTcrrj7SMvtuyTd7e4u6S4z2ympmrPOdJO3j2RmiyWtcfdvNu9qt99OlPTgBPR3MjpiH7a7b5ofR7nM7CRJn5f0QXf/tJn1iuOoFe9HY8D70STGlyFRMM5oT37flnSOJDU/TrxzVP0iNeY8ysyWq3H26Cs560w3eftIkl4g6estt9vtt23ldnNSa7cPx7JfZwwzO17S9ZLe5u4fb97NcfREvB+NDe9HwAzBGe3J7/OSzjazm9WYF/ra5hexPuru50r6mKSrzOwmNa5rdZGkWyS9sHWd7nR9wuTtI0laI+melnWO2G8z7SzbqH10xD5U4yzb6PtmlFH76FJJx0h6h5kdnqv9Mkkf4jji/SjC+9EU4c6VIVE4cz4mAQAAM9yi6lJ/7rz1pbX/1b1X/cDdJ/2XgTdt2lTYwPDUU0+1otoai1Vv/1Khg9r73v3ScfefM9oAAACSnDPaKBhztAEAAIAScEYbAACgccWabncC0wwDbQAAABd51ygcA20AAABJ4hLsKBhztAEAAIASMNAGpjgzO9PM+hO19Wb2fTP7jpn9UfO+OWb2WTO70cy+bGbLUssG2wyXNbOnmNlNzW18yMwqzfv/2sy+Z2Y3m9kZ0bKJ7YbLHs1jSy3b0lafmT1sZn9nZhea2QNm1t/8ud3MPtCy3L+16Wu/mT215fZsM7uv+RhuMbNbon0MYGK5JM+8tB/MTAy0gSnMzN4q6aOSZrep9Uj6B0m/Lek3JL2+eXXDN0i6091/XdKnJF0WLNtum2NZ9n2SLmtuwyS9zMye3Vz+TEmvkvSB1LLBQ85b9mge2xHLttneN9z9rc3fP+3ufe7eJ+nZkp5lZkedievud6vx+AFMJu6NqSNl/WBGYqANTG2bJZ2bqD1N0t3u/qi7D0u6SY1LPz9fjctiS9J1kl4ULHs07bZ6jqT/HLWN50u63hsekFRrnkVut2xK3rJH89jaLTtWCyQtlrT7KNYBgLbMrGJmH25+6tZvZk/pdp9QDL4MCUxh7v5ZM1uVKC/UEweCeyUtGnV/u/ta7z+adluZP37Z2dZt7GyzXrtlU/KWPZrH1m7ZyPlm9lxJJ0raI+lyd/+5ma3IWa8Vnx8Dk1gXp3i8XNJsd3+umZ0l6QrFn+5hiuCMNjB97VHjzOthCyTtGnV/u/ta7z+adltlbeqp9dotm5K37NE8tnbLRj7t7r8h6cXN5e/KWf6ApN6W2/Ob9wGYrLo3deSxT9jc/RZJk/5S7RgbzmgD09dPJJ1iZkskDaoxXeK9kk6WdI6k70l6iaQbg2WPpt1Wt5lZn7v3N7fxTUl3S/o7M3uvpJWSKu6+w8zaLZuSt+y3j+KxtdsPudz9XjP7E0lXm9kzgkVvlfT7ku5o3n6JpO+PZRsAJt5ePfrVr/k1S0vcxGwz29hy+0p3v7L5++hP3upmVnP3kRL7gwnAQBuYZszsfEnz3f1KM3uLpK+q8enVx919wMw+JOmTZnaTpGFJ57v7ocSyz5J0obu/+XD7wbJPl/Qmd3+jpD+X9BEzm6XGQPcad6+b2Y2SvtNc70+aTR6xbPNx9De/eNgqtez1ktZJOprHdsSyY93H7v41M/uapL+V9CVJvz3qD+j5kt4j6f+Y2a1qnMl+RFKY5gKge9z9d7q4+dGfvFUYZE8P5lxuFECCmc2TdKm7/2UXtv2/Wgf4Xdh+n6SL3b3whJDmvPp/c/ezim4bwNRjZr8vab27X9ico/3X7v6SbvRl06ZNhQ0MTz31VCuqrbFY9fYvFTqove/dLx13/5mjDSBSU+PMbDdc0aXttvpNM/u7IhtspgkckbsNYEb7vKSDZnazGnGkf9bl/qAgTB0BkOTuXYuvc/cHu7Xt5vb7JR1XQrt3S+JMNoDHuHsm6eJu9wPF44w2AAAAUAIG2gAAAEAJGGgDAAAAJWCgDQAAAJSAgTYAAABQAgbaAAAAQAkYaAMAAAAlYKANAAAAlICBNgAAAFACBtoAAABACRhoAwAAACVgoA0AAACUgIE2AAAAUAIG2gAAAEAJGGgDAAAAJWCgDQAAAJSAgTYAAABQAgbaAAAAQAkYaAMAAAAlYKANAAAAlICBNgAAAFACBtoAAABACRhoAwAAACVgoA0AAACUoOOBtpn1mdm/jbrv3WZ24VG2c6aZ9Sdq683s+2b2HTP7o+Z9c8zss2Z2o5l92cyWpZYFAAAAuqWrZ7TN7K2SPippdptaj6R/kPTbkn5D0uvN7HhJb5B0p7v/uqRPSbosWBYAAADoitIG2mZ2vZnNyllss6RzE7WnSbrb3R9192FJN0l6gaTnS/pKc5nrJL0oWBYAAADoito41//NUdM+nizpryTJ3X87b2V3/6yZrUqUF0ra3XJ7r6RFo+5vd1/r/QAAAEBXjHeg/Q13f9XhG2b27tSCZjZf0rXNmze4++U5be+RtKDl9gJJu0bd3+6+1vsBAACArhjvQHvM3H1QUt9RrPITSaeY2RJJg2pMBXmvpJMlnSPpe5JeIunGYFkAAACgK0obaJvZ9ZLWNedMH81650ua7+5XmtlbJH1VjbnkH3f3ATP7kKRPmtlNkoYlne/uh9otW+gDAgAAAI6CuXu3+wAAAIBJYNOmTYUNDE899VQrqq2xWPX2LxU6qL3v3S8dd/+5YA0AAABQggmboz0Z9ff3v7BSqbxf0qlZlqlSqUjSpizL/rSvr++bXe4eAKBL+vv7V1cqlbdIenWWZQsqlcpeSf+SZdn7+vr6Nne7fwCmhnENtM2sT9LFh5NHzOwVkv5G0jnu/sAY1v8tSe+UdEjSw5L+wN33t9TnSPoXScepEdn3Gnffbmbr1YgRHFFjPvZHUsumtt3f3/+OSqXy31esWKHly5ert7dXQ0ND2rp166kDAwPf+G9XXq8N9xw66n0CAJjaTlta1Z8+Z66edOJyPWnlisN/HxZu3br1jVu3bn3jrbfeqlmz8i4T0d6pp52mTXfeWXCPUYSJnuaAmaGwqSNm9v9J+gtJvzWWQXbTByW93N1fIOnnkl43qn40V4E8YtnURptnsv/72rVrtXr1as2ZM0eVSkVz5szR6tWrtXbtWp331Nl66jHMrAGAmWTZHNOfPmeuznj2Wq055SlH/H145jOfqX379qler3e7qwCmgEJGkmZ2gaQ/k/Qid/9F877zzez1Oav2HV5ejbPrB0fVj+YqkO2WbatSqbx/xYoVWrSo/TVtFi1apBUrVuiCU4+4MjwAYBo758mz9KSVy8O/DyeeeKIOHhz95woAjlTEQPvXJb1e0hK1TEVx90+7+5XRiu6+TZLM7FxJL1TjTHSro7kKZLtlU05dvnx51DUtX75cKxdUw2UAANPLr63o0ZNWrgiXWbFihYaHjyq5FtOAma03s++b2XfM7I/a1J9iZjeZ2Y1m9iEzqzTvv9zMvmtmtzSn3GIGKWKgvU3S2ZL+l6R/OXxgtWNmbzKz/ubPiuZ9fybpzyX9jruPPkVwNFeBbLdsW1mWqbe3N3xQvb298iwLlwEATC+zzMf094GpIzNLMG211fskXdacwmqSXmZmvyLprObPqyT948T1GpNBEQPtu939oLv/kxoXkPnL1ILu/k/u3tf8GTCzv1TjjPiL3H1Hm1W+rcZVIKU2V4E0s1lqTBv5TmLZtiqVioaGhsIHNTQ0JKswRxsAZpJhtzH9fahW+cRzhklNW231HEn/2fz9OjXGNrdJerE3LlpysoKTgJieih5JXiTpj83shXlztJv/E/xrScslXdc8y/2GZu365iD6Q5Ke0bwK5Osl/a27H5J0+CqQ39HjV4E8Ytmgn5u2bt0aPpCtW7dqy17OWADATHLzwCE9sCW+sPDAwEDHqSOYslLTVluZP34VwMfq7j5iZpdLulbSJ8ruKCaXccX7uXu/pP6W29slrRzjur+Q1Padyt1/u/nrsKTz2tQ3SNow6r797ZZtJ8uyPx0YGPjG0qVL237hZffu3RoYGNA/b+LLLgAwk3z5nmE9f+VWnXDcsuTfh23btmnhwoVd6B0mmpm9U42whWdK+m5Lqd0U1SxVd/e/NLN3S7rFzG509ymZxV6v13Xw4EENDw+rXq+rWq1q1qxZmj17Np/yJMzIC9b09fV9s7+//6/uuOOOdjnaGhgYUJZlf/WVt73kf3S7rwCAidXf3/+SO+6445oVK1b0LF++vKfl78OhgYGBQ1mWvWLt2rXXddS4u04tuL8oj7tfJj02R/vHZrZE0qAa00beO2rx28ysr3kS8iWSvmlmvynp9939T9RIVjukJw7Ip4zh4WHt27dPy5cvP2LctHXrVs2bN49Petqwxz/lmHmaedr/KOm0litD3pll2SVcGRIAZq7mlSH/TNIFWZbNr1Qqg5L+Ocuyf+DKkDNTy8XyKmpMW/2AmT1d0pvc/Y1m9suSPqLGp/U/kXQ4meSf1DgjXpX0MXf/yMT3fuw2bdp0xMCwXq9rz549euYzn5n8pOeHP/yhFi5c+IQz2xN9EaBVb/9SoYPa+9790nH3f0YPtAEAAPC4dgPtffv2aenSpVq9enVyvbvvvls7d+7UvHnzHruPgfY4vgxpZn1m9m/j2biZzTGzzzYzJ79sZsvaLPPXZvY9M7vZzM5o3pfKqjxiWQAAAHRueHhYedcfIV++vW7P0T582fS/MbNXqXHZ9EsOF83s2WrkVZ4p6SRJn5X0q3o8q7LfzD6sRlbl/YllAQAAMAbtzkL39/dnvb294dndZr58duqpp3btW5FFnIEuWuFB0Wb2FjP73UTtb8zs4pa78i6b/nxJ13vDA5JqzbPeR2RVBssCAACgQ5VKZe9Y8uWb32VAi8IH2u7+Pnf/Yut9ZvZfzKxf0oWS3tLMzD5b+ZdNT+VWtsuqHEvGJQAAAI7Ov2zduvVQtECz/s8T1J8pY0Iufejun3H3PklXSXpf88qQNyj/sumpy623y6pMLQsAAIAOZVn2voGBgUO7d+9uW29ef+RQlmX/MMFdm/S6fY3xvMumf1vSi82sYmZPklRpXqr9NjPrG7VealkAAAB0qK+vb3OWZa+444479m/evPnQgQMHlGWZDhw4oM2bNx+644479mdZ9gqiL4/Ucbxfc6D7OUn3tNx9vqR1ku4ePX0k0cZcSZ+UdKIaV4E8390fMrO/k3SNu3/PzP5GjcF0RdKfuftN7bIq3b3ebtmOHhwAAACegHz5o0eONgAAAFCCbk8dAQAAQInM7KVm9kMz+5mZXW1mCztZ5ii3uczMrjOzH5vZJjP7tfFsy8yuMLMHzOz25s9n2ixzppltNLOfmNnXzezEDvo918w+3WzjZ2b28ub9L2rZ9u3NmpvZc8L2OKMNAAAwPTWjjn8k6Xnu/nMze4+kBe7+xqNZpoPt/ruk2939f5rZsyR9SdIpkuZ1si0z+46kP3f3mxP1WZI2S3qVu3/bzN4gab27n9Nu+WA7fydpsbu/vvmdv1skneHuW0Ytd42kn7v7X4TtMdAGAACYnszsv6rxHbiXNm+vknSHGoNJz1tGUo+k96hxUcCqpNsk/am772nZxlWS+t39qubtmhppcCe7+/bmff2S3i9pTl5/2jyGXjUinL8sabWku9X4Lt4DLcs8T9KV7v6M5u1ZakQ9L3f3nWb2l5J+X43ZHPdJeqO7b22zrZ83+/f9lsf2Q3d/X8syr5b0ZklnuftIuz4fxtQRAACA6eskSQ+23N6ixrVHFoxxmbdLGpH0HHdfK2mrpHfnbHOpGulv20e1uXKM/RltuaRvSPoLSc9S4yzzf5hZ65Ugn9Cuuw9L2i5phZn9gaTT1Dgz/Sw1BuwfTWyrXf9WHr7RHMC/S9Kb8wbZUvcvwQ4AAIDypE6q1se4zDo1zmyf3RzXzpL0sCSZ2Xcl9Up6kqTfNLM3qxG3fHnQ3lj68wTufq8ej4OWmb1X0jskrZJ07xgfwxmSNjYfQ1XS3MTy7dpp7dsrJG0ea7IdA20AAIDp6wFJZ7bcXiHpUXffN5ZlzKwq6RJ3v06SzGy+pNmS5O5nNu+7SkdOHZGZHePuj7a0uUWNKSB5/XkCM3umpLXu3nrlSZPUerXKB9SIiz68To8aZ9YH1BhYv8fdP9Ss9Uo6xsyWq3F2+7BzWtp5qKV/t7cs818kfSLV19GYOgIAADB9XS/pLDM7pXn7Ykn/cRTLfFXSm8xslplV1LiOybuiDTanVHxJ0h9Ljw2Uny6pf4z9GS2T9H4z+6Xm7TeoMW+69QuK35V0bEu6yUWSvuPuu5qP4XUt6Sb/XdI/u/tWd39Wy8/WZl9e3+z3Skm/I+na5m2T9AJJX8/p72P4MiQAAMA0ZmbnqDE4PpzM8QeSnizpo805y22XcfdHzGyOpPdK6lPjzPDtkl7f+mXIxDaPV2Me9C9JcjUSQ66PtpXT3qvVmC9eVePM+B+qMXf8y5LOcfetZnaGpH9SI9lkZ7Pd+5r/QfhrSec1+/KApNe5+0Cb7cyX9CFJz25u653u/i/N2jJJW929J+rrE9pjoA0AAAAUj6kjAAAAQAkYaAMAAAAlYKANAAAAlICBNgAAAFACBtoAAABACRhoAwAAACVgoA0AAACUgIE2AAAAUIL/B0MqjO9I3BtFAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Viewer.setProjection(1)\n", "Viewer.setPlane(61)\n", "Viewer.ax.get_figure()" ] }, { "cell_type": "markdown", "id": "db9e5dca", "metadata": {}, "source": [ "The reason for the unusual axis (-1,-2,0) being used is that it is 90$^\\circ$ to the (1,0,0) axis and thus is the axis along which the data are binned when doing a regular voxel binning in 3D. Had the data been for a crystal with a 90$^\\circ$ symmetry, i.e. tegragonal, and the scattering plane ($H$,$0$,$L$), the plotting axes would have been (1,0,0) and (0,0,1).\n", "\n", "A small side note is that the alingment of this sample was not perfect and thus the dispersion is not located exactly on top of (1,0,0) as seen most easily in the two cut last cuts. " ] } ], "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.9.7" } }, "nbformat": 4, "nbformat_minor": 5 }