{ "cells": [ { "cell_type": "markdown", "id": "ffde92de-43bc-40e2-9d63-6d8992580cc0", "metadata": {}, "source": [ "# Time dependence - amplitude phase jump" ] }, { "cell_type": "code", "execution_count": 1, "id": "d7e6742d-2fa0-45fe-b5af-05ad2d529b28", "metadata": {}, "outputs": [], "source": [ "import nexus as nx\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "code", "execution_count": 2, "id": "e1deab56-29e2-41ce-aeb9-45dec1e67a44", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlcAAAHgCAYAAACIHEjUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/GU6VOAAAACXBIWXMAAA9hAAAPYQGoP6dpAAB4TElEQVR4nO3deVhU1f8H8PcM+77vIqAoiqAgCuJekkiuqbmWpmZqaipZSplb31+Yptlimi3aormlZmom4q64Ae6CG4vKviPIOvf3BzE5gQpyYWDm/XoeHp1zzz33c+cC8+Gcc8+VCIIggIiIiIhEIVV2AERERESqhMkVERERkYiYXBERERGJiMkVERERkYiYXBERERGJiMkVERERkYiYXBERERGJiMkVERERkYiYXBERERGJiMkVkRIcPXoUEokER48elZe98cYbcHZ2bpDjx8fHQyKRYOPGjQ1yvErVnTc1Dc7OznjjjTfkrxv6Wvbu3Ru9e/eu9+Mo62eDVAuTK1Ib33zzDSQSCfz8/JQdSo0UFhZi8eLFjT4ReeONNyCRSJ759fgHM6mmzZs3Y/Xq1coOo0aaUqzU9GgqOwCihrJp0yY4Ozvj3LlzuH37NlxdXZUdkoLvvvsOMplM/rqwsBBLliwBgAb5i/15TZkyBQEBAfLXcXFxWLhwId566y306NFDXt6yZUv4+fnh0aNH0NbWVkaoJKKePXtWuZabN2/G1atXMXv2bOUFVkNPitXJyQmPHj2ClpaWcgIjlcDkitRCXFwcTp8+jZ07d2LKlCnYtGkTFi1apOywFDTVX+b+/v7w9/eXv75w4QIWLlwIf39/vPbaa1Xq6+rqNmR4VE+kUqlKXkuJRKKS50UNi8OCpBY2bdoEMzMz9O/fH8OHD8emTZuq1Kmca/HZZ59hzZo1aNGiBfT19dG3b1/cu3cPgiDg448/RrNmzaCnp4fBgwcjKytLoQ1nZ2cMGDAABw8ehJeXF3R1deHu7o6dO3c+M8bH51zFx8fDysoKALBkyRL50NrixYsBPHn+SXXztnJycvDGG2/AxMQEpqamGD9+PHJycqqNISYmBsOHD4e5uTl0dXXRqVMn7Nmz55mx11R183R69+4NDw8PXL58Gb169YK+vj5cXV2xY8cOAMCxY8fg5+cHPT09uLm54dChQ1XaffDgASZOnAgbGxvo6OigXbt2+PHHH6vUS0xMRExMTI3j3LZtG5YsWQIHBwcYGRlh+PDhyM3NRXFxMWbPng1ra2sYGhpiwoQJKC4ulu/fq1cvdOjQodq23dzcEBgY+NTj//HHH+jfvz/s7e2ho6ODli1b4uOPP0Z5eblCvbq+d4sXL4ZEIkFMTAxGjBgBY2NjWFhYYNasWSgqKqrRe1R5LXv37o19+/YhISFB/v1a+b24ceNGSCQSxMfHP7WNSuvXr0fLli2hp6cHX19fnDhxotoYiouLsWjRIri6ukJHRweOjo54//33Fa5FdZ4Wa3Vzrt544w0YGhoiMTERAwYMgKGhIRwcHLBmzRoAwJUrV/Diiy/CwMAATk5O2Lx5c5Vj5uTkYPbs2XB0dISOjg5cXV3x6aefKvRWk+pgzxWphU2bNmHo0KHQ1tbG6NGjsXbtWpw/fx6dO3eutm5JSQlmzpyJrKwsLF++HCNGjMCLL76Io0ePYt68ebh9+za++uorzJ07t8qH+K1btzBy5EhMnToV48ePx4YNG/Dqq6/iwIEDeOmll2oUr5WVFdauXYtp06bhlVdewdChQwEA7du3r9V5C4KAwYMH4+TJk5g6dSratm2LXbt2Yfz48VXqXrt2Dd26dYODgwPmz58PAwMDbNu2DUOGDMHvv/+OV155pVbHro3s7GwMGDAAo0aNwquvvoq1a9di1KhR2LRpE2bPno2pU6dizJgxWLFiBYYPH4579+7ByMgIAJCamoouXbpAIpFgxowZsLKywl9//YVJkyYhLy9PYdhn3LhxOHbsGARBqFFcoaGh0NPTw/z58+XXXEtLC1KpFNnZ2Vi8eDHOnDmDjRs3wsXFBQsXLgQAvP7665g8eTKuXr0KDw8PeXvnz5/HzZs3sWDBgqced+PGjTA0NERwcDAMDQ1x+PBhLFy4EHl5eVixYoVo712lESNGwNnZGaGhoThz5gy+/PJLZGdn4+eff67R+wQAH374IXJzc3H//n18/vnnAABDQ8Ma71/phx9+wJQpU9C1a1fMnj0bd+/exaBBg2Bubg5HR0d5PZlMhkGDBuHkyZN466230LZtW1y5cgWff/45bt68id27d4saa3l5OYKCgtCzZ08sX74cmzZtwowZM2BgYIAPP/wQY8eOxdChQ7Fu3TqMGzcO/v7+cHFxAVAxxN+rVy88ePAAU6ZMQfPmzXH69GmEhIQgOTmZc79UkUCk4i5cuCAAEMLCwgRBEASZTCY0a9ZMmDVrlkK9uLg4AYBgZWUl5OTkyMtDQkIEAEKHDh2E0tJSefno0aMFbW1toaioSF7m5OQkABB+//13eVlubq5gZ2cneHt7y8uOHDkiABCOHDkiLxs/frzg5OQkf52eni4AEBYtWlTlnHr16iX06tWrSvl/29i9e7cAQFi+fLm8rKysTOjRo4cAQNiwYYO8vE+fPoKnp6fC+chkMqFr165Cq1atqhzrSc6fP1+l7UrVnXevXr0EAMLmzZvlZTExMQIAQSqVCmfOnJGX//3331XanjRpkmBnZydkZGQoHGvUqFGCiYmJUFhYWOVYz1IZp4eHh1BSUiIvHz16tCCRSISgoCCF+v7+/grve05OjqCrqyvMmzdPod4777wjGBgYCA8fPnzq8R+PudKUKVMEfX19hetT1/du0aJFAgBh0KBBCsd6++23BQDCpUuX5GVOTk7C+PHj5a+ru5b9+/dXeB8qbdiwQQAgxMXFKZT/t42SkhLB2tpa8PLyEoqLi+X11q9fLwBQ+J7/5ZdfBKlUKpw4cUKhzXXr1gkAhFOnTlWJ43FPirXy98Dj79P48eMFAMInn3wiL8vOzhb09PQEiUQibNmyRV5e+f4//nP78ccfCwYGBsLNmzcVjjV//nxBQ0NDSExMfGqs1PRwWJBU3qZNm2BjY4MXXngBQMWcipEjR2LLli1VhlkA4NVXX4WJiYn8deXdha+99ho0NTUVyktKSvDgwQOF/e3t7RV6eYyNjTFu3DhER0cjJSVF1HN7lv3790NTUxPTpk2Tl2loaGDmzJkK9bKysnD48GGMGDEC+fn5yMjIQEZGBjIzMxEYGIhbt25VOU8xGRoaYtSoUfLXbm5uMDU1Rdu2bRXu7qz8/927dwFU9Mz9/vvvGDhwIARBkMedkZGBwMBA5ObmIioqSr7/0aNHa9xrBVT0dD0+F87Pzw+CIGDixIkK9fz8/HDv3j2UlZUBAExMTDB48GD89ttv8uOVl5dj69atGDJkCAwMDJ56XD09Pfn/K69Hjx49UFhYWGVY83nfu8dNnz5d4XXl98f+/fufGqfYLly4gLS0NEydOlVhonzlsPbjtm/fjrZt26JNmzYK1/3FF18EABw5ckT0+N588035/01NTeHm5gYDAwOMGDFCXl75/j/+Pm/fvh09evSAmZmZQqwBAQEoLy/H8ePHRY+VlIvDgqTSysvLsWXLFrzwwguIi4uTl/v5+WHlypUIDw9H3759FfZp3ry5wuvKX+qPD0k8Xp6dna1Q7urqColEolDWunVrABXzOWxtbetwRrWTkJAAOzu7KkMebm5uCq9v374NQRDw0Ucf4aOPPqq2rbS0NDg4ONRLnM2aNavynpmYmDzzPU9PT0dOTg7Wr1+P9evXV9t2Wlrac8dVm+8FmUyG3NxcWFhYAKhIzLZu3YoTJ06gZ8+eOHToEFJTU/H6668/87jXrl3DggULcPjwYeTl5Slsy83NVXj9vO/d41q1aqXwumXLlpBKpVXmSNW3hISEauPR0tJCixYtFMpu3bqFGzduyOcm/lddrnt1dHV1qxzLxMTkie//4+/zrVu3cPny5QaLlZSPyRWptMOHDyM5ORlbtmzBli1bqmzftGlTleRKQ0Oj2raeVF6bnhCxSCSSao9bXU9cTVROqp07d+4TJ1vX59IVz/ueV8b92muvVTuPDKj9PDUx4gKAwMBA2NjY4Ndff0XPnj3x66+/wtbWVmHZiurk5OSgV69eMDY2xtKlS9GyZUvo6uoiKioK8+bNqzIBuj6+X/+bLNTVk9p73u9XoOLae3p6YtWqVdVu/29yWVd1eZ9lMhleeuklvP/++9XWrfzji1QHkytSaZs2bYK1tbX8rp7H7dy5E7t27cK6desUhmHqqrIX6PEPlJs3bwJArVZgf9oHnJmZWbXDO5V/+VdycnJCeHg4Hj58qNB7FRsbq1CvsldAS0vrmR/+jYmVlRWMjIxQXl7e6OLW0NDAmDFjsHHjRnz66afYvXs3Jk+e/MQP40pHjx5FZmYmdu7ciZ49e8rLH+95FdutW7fkk6+Biu9hmUxW6ycGPOl71szMDACq3KVa3fdrZTyVw3sAUFpairi4OIU7MFu2bIlLly6hT58+z5UMip1APk3Lli3x8OHDRvc9SvWHc65IZT169Ag7d+7EgAEDMHz48CpfM2bMQH5+vqhLDQBAUlISdu3aJX+dl5eHn3/+GV5eXrUaEtTX1wdQ9QMJqPhlHRMTg/T0dHnZpUuXcOrUKYV6L7/8MsrKyrB27Vp5WXl5Ob766iuFetbW1ujduze+/fZbJCcnVzne48dpTDQ0NDBs2DD8/vvvuHr1apXt/427pksxiOX1119HdnY2pkyZgocPH1a77td/VSZfj/d8lJSU4Jtvvqm3OP/7x0fl90dQUFCt2jEwMKgybAlUfL8CUJhbVF5eXmUot1OnTrCyssK6detQUlIiL9+4cWOVn4MRI0bgwYMH+O6776oc79GjRygoKHiuWOvDiBEjEBERgb///rvKtpycHPlcPVId7LkilbVnzx7k5+dj0KBB1W7v0qULrKyssGnTJowcOVK047Zu3RqTJk3C+fPnYWNjgx9//BGpqanYsGFDrdrR09ODu7s7tm7ditatW8Pc3BweHh7w8PDAxIkTsWrVKgQGBmLSpElIS0vDunXr0K5dO4U5OgMHDkS3bt0wf/58xMfHy9fcqu5DZc2aNejevTs8PT0xefJktGjRAqmpqYiIiMD9+/dx6dKlOr839WHZsmU4cuQI/Pz8MHnyZLi7uyMrKwtRUVE4dOiQwlpktV2Koa68vb3h4eEhn3zdsWPHZ+7TtWtXmJmZYfz48XjnnXcgkUjwyy+/1GvMcXFxGDRoEPr164eIiAj8+uuvGDNmzBPX6noSHx8fbN26FcHBwejcuTMMDQ0xcOBAtGvXDl26dEFISAiysrJgbm6OLVu2VEkqtLS08L///Q9TpkzBiy++iJEjRyIuLg4bNmyoMufq9ddfx7Zt2zB16lQcOXIE3bp1Q3l5OWJiYrBt2zb8/fff6NSpU61jrQ/vvfce9uzZgwEDBuCNN96Aj48PCgoKcOXKFezYsQPx8fGwtLSsl2OTcrDnilTWpk2boKur+8S1paRSKfr3748DBw4gMzNTtOO2atUKW7duxf79+zF//nyUlpZi69atz1w4sjrff/89HBwcMGfOHIwePVq+OGTbtm3x888/Izc3F8HBwdizZw9++eWXKh/eUqkUe/bswdixY/Hrr7/iww8/hIODA3766acqx3J3d8eFCxfQv39/bNy4EdOnT8e6desglUrl6zc1RjY2Njh37hwmTJiAnTt3YsaMGfjiiy+QlZWFTz/9VNnhYdy4cQBQo4nsAGBhYYG9e/fCzs4OCxYswGeffYaXXnoJy5cvr7cYt27dCh0dHcyfPx/79u3DjBkz8MMPP9S6nbfffhtjxozBhg0bMGbMGIW7Ujdt2oSuXbti2bJl+OSTT/DCCy9g2bJlVdp466238M033yApKQnvvfceTpw4gT179lSZQyWVSrF7924sW7YMV65cwdy5c7FkyRKcP38es2bNeuY8pqfFKjZ9fX0cO3YM7733Ho4ePYpZs2Zh2bJluHXrFpYsWVLlTkhq+iSCMmbjEqkoZ2dneHh4YO/evcoOhRqJL774AnPmzEF8fHyVuw+VbfHixViyZAnS09PZc0IkIvZcERHVE0EQ8MMPP6BXr16NLrEiovrDOVdERCIrKCjAnj17cOTIEVy5cgV//PGHskMiogbE5IqISGTp6ekYM2YMTE1N8cEHHzzxpgoiUk2cc0VEREQkIs65IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiUnsbN26ERCJBfHx8k2i3sR6XiCowuSIiaoJOnz6NxYsXIycnR9mhENF/MLkiIqonr7/+Oh49egQnJyfR2z59+jSWLFlSbXJVn8clomfjswWJiERWUFAAAwMDaGhoQENDo8GPr6zjElEF9lwRkdKsX78eHTt2hL6+PiQSicJXixYtqt1n8eLFkEgkiImJwYgRI2BsbAwLCwvMmjULRUVFVepHR0cjKCgIxsbGMDQ0RJ8+fXDmzJlnxpaQkIC3334bbm5u0NPTg4WFBV599dUq85gq47l+/TrGjBkDMzMzdO/eHYDi3Kf4+Pgq5/j4V22Ou3jxYrz33nsAABcXF3kblXWqm3NV0/eh8nxu376NN954A6ampjAxMcGECRNQWFj4zPeNiNhzRURKMmfOHKxevRp9+/bFhAkTcP/+fXz++ecoLS3FgAED4OPj89T9R4wYAWdnZ4SGhuLMmTP48ssvkZ2djZ9//lle59q1a+jRoweMjY3x/vvvQ0tLC99++y169+6NY8eOwc/P74ntnz9/HqdPn8aoUaPQrFkzxMfHY+3atejduzeuX78OfX19hfqvvvoqWrVqhU8++QSCIFRpz8rKCr/88otCWWlpKebMmQNtbe1aHXfo0KG4efMmfvvtN3z++eewtLSUH6M6z/M+jBgxAi4uLggNDUVUVBS+//57WFtb49NPP33ie0ZE/xCIiBrY8ePHBQDCtGnTFMqXLFkiABDOnTv3xH0XLVokABAGDRqkUP72228LAIRLly7Jy4YMGSJoa2sLd+7ckZclJSUJRkZGQs+ePeVlGzZsEAAIcXFx8rLCwsIqx46IiBAACD///HOVeEaPHl2lfnXt/jdmDQ0N4fDhw7U+7ooVK57Y9n+PW9P34fHzmThxokL5K6+8IlhYWFR7HkSkiMOCRNTgPv/8c5ibm2PFihUK5b169QIA3Lx585ltTJ8+XeH1zJkzAQD79+8HAJSXl+PgwYMYMmSIwhCjnZ0dxowZg5MnTyIvL++J7evp6cn/X1paiszMTLi6usLU1BRRUVFV6k+dOvWZMT/u559/xjfffIPly5fjhRdeeO7jPsvzvg//PZ8ePXogMzPzqe8ZEVVgckVEDaqsrAxhYWEICgqCgYGBwraSkhIAgLGx8TPbadWqlcLrli1bQiqVyucZpaeno7CwEG5ublX2bdu2LWQyGe7du/fE9h89eoSFCxfC0dEROjo6sLS0hJWVFXJycpCbm1ulvouLyzNjrnTx4kVMnToVo0ePRnBwcJ2O+yzP+z40b95c4bWZmRkAIDs7u9YxEKkbzrkiogYVHx+Phw8fwsPDo8q2yMhIABUf+rX1+KRwMcycORMbNmzA7Nmz4e/vDxMTE0gkEowaNQoymaxK/cd7nJ4mOzsbw4YNQ+vWrfH999/X+bj15Ul3GwrVzCcjIkVMroioQeXn5wOAwiRuoOJDe/v27WjXrh1cXV2f2c6tW7cUeotu374NmUwGZ2dnABWTu/X19REbG1tl35iYGEilUjg6Oj6x/R07dmD8+PFYuXKlvKyoqKhOi3bKZDKMHTsWOTk5OHToUJVJ8bU5bk2Tybq+D0RUexwWJKIGVTncdOjQIYXy1atXIyoqCiEhITVqZ82aNQqvv/rqKwBAUFAQgIqel759++KPP/5QWJIgNTUVmzdvRvfu3Z86/KihoVGll+arr75CeXl5jeKrzpIlS/D333/jt99+e+IwYk2PWzmk+qxkr67vAxHVHnuuiKhBWVhYYMiQIdi9ezfGjh2Lbt264eTJk/jtt9/w5ptvYuzYsTVqJy4uDoMGDUK/fv0QERGBX3/9FWPGjEGHDh3kdf73v/8hLCwM3bt3x9tvvw1NTU18++23KC4uxvLly5/a/oABA/DLL7/AxMQE7u7uiIiIwKFDh2BhYfFc533lyhV8/PHH6NmzJ9LS0vDrr78qbH/ttddqddzKpSo+/PBDjBo1ClpaWhg4cGCVeWx1fR+IqPaYXBFRg9uwYQNMTU3xxx9/4Pfff4e7uzt++OEHTJw4scZtbN26FQsXLsT8+fOhqamJGTNmVLn7sF27djhx4gRCQkIQGhoKmUwGPz8//Prrr09d4woAvvjiC2hoaGDTpk0oKipCt27dcOjQIQQGBj7XOWdmZkIQBBw7dgzHjh2rsr0yuarpcTt37oyPP/4Y69atw4EDByCTyRAXF1dtclWX94GIak8icHYiETUhixcvxpIlS5Ceni5fPJOIqDHhnCsiIiIiETG5IiIiIhIRkysiIiIiEXHOFREREZGI2HNFREREJCIuxdDAZDIZkpKSYGRkJPrjOoiIiKh+CIKA/Px82NvbQyp9et8Uk6sGlpSUxEdNEBERNVH37t1Ds2bNnlqHyVUDMzIyAlBxcfjICSIioqYhLy8Pjo6O8s/xpxIaqa+//lpwcnISdHR0BF9fX+Hs2bNPrb9t2zbBzc1N0NHRETw8PIR9+/bJt5WUlAjvv/++4OHhIejr6wt2dnbC66+/Ljx48EChjczMTGHMmDGCkZGRYGJiIkycOFHIz89XqHPp0iWhe/fugo6OjtCsWTPh008/rdV55ebmCgCE3NzcWu1HREREylObz+9GOaF969atCA4OxqJFixAVFYUOHTogMDAQaWlp1dY/ffo0Ro8ejUmTJiE6OhpDhgzBkCFDcPXqVQBAYWEhoqKi8NFHHyEqKgo7d+5EbGwsBg0apNDO2LFjce3aNYSFhWHv3r04fvw43nrrLfn2vLw89O3bF05OToiMjMSKFSuwePFirF+/vv7eDCIiImpSGuVSDH5+fujcuTO+/vprABWTwB0dHTFz5kzMnz+/Sv2RI0eioKAAe/fulZd16dIFXl5eWLduXbXHOH/+PHx9fZGQkIDmzZvjxo0bcHd3x/nz59GpUycAwIEDB/Dyyy/j/v37sLe3x9q1a/Hhhx8iJSUF2traAID58+dj9+7diImJqfY4xcXFKC4ulr+u7FbMzc3lsCAREVETkZeXBxMTkxp9fje6nquSkhJERkYiICBAXiaVShEQEICIiIhq94mIiFCoDwCBgYFPrA8Aubm5kEgkMDU1lbdhamoqT6wAICAgAFKpFGfPnpXX6dmzpzyxqjxObGwssrOzqz1OaGgoTExM5F+czE5ERKTaGl1ylZGRgfLyctjY2CiU29jYICUlpdp9UlJSalW/qKgI8+bNw+jRo+XZZ0pKCqytrRXqaWpqwtzcXN7Ok45Tua06ISEhyM3NlX/du3ev2npiKCuX4UZyHhphZyQREZHaaHTJVX0rLS3FiBEjIAgC1q5dW+/H09HRgbGxscJXfXlvx2UEfXEC3x6/W2/HICIioqdrdMmVpaUlNDQ0kJqaqlCempoKW1vbavextbWtUf3KxCohIQFhYWEKiY6trW2VCfNlZWXIysqSt/Ok41RuU7Zd0Q8AAGsO31ZyJEREROqr0SVX2tra8PHxQXh4uLxMJpMhPDwc/v7+1e7j7++vUB8AwsLCFOpXJla3bt3CoUOHYGFhUaWNnJwcREZGyssOHz4MmUwGPz8/eZ3jx4+jtLRU4Thubm4wMzN7/pMWWTmHBYmIiJSm0SVXABAcHIzvvvsOP/30E27cuIFp06ahoKAAEyZMAACMGzcOISEh8vqzZs3CgQMHsHLlSsTExGDx4sW4cOECZsyYAaAisRo+fDguXLiATZs2oby8HCkpKUhJSUFJSQkAoG3btujXrx8mT56Mc+fO4dSpU5gxYwZGjRoFe3t7AMCYMWOgra2NSZMm4dq1a9i6dSu++OILBAcHN/A79HTlMiZXREREytIoV2gfOXIk0tPTsXDhQqSkpMDLywsHDhyQTx5PTExUeK5P165dsXnzZixYsAAffPABWrVqhd27d8PDwwMA8ODBA+zZswcA4OXlpXCsI0eOoHfv3gCATZs2YcaMGejTpw+kUimGDRuGL7/8Ul7XxMQEBw8exPTp0+Hj4wNLS0ssXLhQYS2sxoAdV0RERMrTKNe5UmW1WSejtpzn7wMAaEoluP3Jy6K2TUREpM6a9DpXVHecc0VERKQ8TK5UEHMrIiIi5WFyRURERCQiJldEREREImJyRURERCQiJldEREREImJyRURERCQiJldEREREImJypaJO3c5QdghERERqicmVivrq8C1lh0BERKSWmFypKD67mYiISDmYXKmoc3FZKCuXKTsMIiIitcPkSoVtu3Bf2SEQERGpHSZXKuyDXVdQyt4rIiKiBsXkSsUN+vqUskMgIiJSK0yuVNyN5DyM+/EcyjnDnYiIqEEwuVIDx2+mo+UH+xGbkq/sUIiIiFQekys1Erj6OMb9eA4FxWXKDoWIiEhlMblSIVLJs+scv5mOdov+xif7b0DGoUIiIiLRMblSIVJJDbKrf6w/fhctPtiPVQdjUVxWXo9RERERqRcmVyqkNslVpS8P34bbggOY//tlZDwsroeoiIiI1AuTK1VS+9xKbsv5e+j0v0MI+uIE7mUVihcTERGRmmm0ydWaNWvg7OwMXV1d+Pn54dy5c0+tv337drRp0wa6urrw9PTE/v37Fbbv3LkTffv2hYWFBSQSCS5evKiwPT4+HhKJpNqv7du3y+tVt33Lli2inXddCILiHKr97/RA/LL+8q87n7yMv2f3xBtdnZ/Yxo3kPPRYfgTO8/fhwNXkKm0SERHR0zXK5Grr1q0IDg7GokWLEBUVhQ4dOiAwMBBpaWnV1j99+jRGjx6NSZMmITo6GkOGDMGQIUNw9epVeZ2CggJ0794dn376abVtODo6Ijk5WeFryZIlMDQ0RFBQkELdDRs2KNQbMmSIaOdeF4+vZdXC0gDu9sYK2zWkErjZGmHxoHaIX9YfUR+9hLd6tnhie1N/jYJLyH6sCruJRyWcl0VERFQTEqERdk34+fmhc+fO+PrrrwEAMpkMjo6OmDlzJubPn1+l/siRI1FQUIC9e/fKy7p06QIvLy+sW7dOoW58fDxcXFwQHR0NLy+vp8bh7e2Njh074ocffpCXSSQS7Nq167kTqry8PJiYmCA3NxfGxsbP3qEWnOfvk/8/5uN+0NXSqPG+KblFWH4gBjujHzyxjpejKVaP9IKzpUGd4iQiImpqavP53eh6rkpKShAZGYmAgAB5mVQqRUBAACIiIqrdJyIiQqE+AAQGBj6xfk1ERkbi4sWLmDRpUpVt06dPh6WlJXx9ffHjjz8+deisuLgYeXl5Cl8NoTaJFQDYmuhi1UgvxC/rj70zu6O1jWGVOhfv5aD3Z0fhPH8ftp5P5HMLiYiIqtHokquMjAyUl5fDxsZGodzGxgYpKSnV7pOSklKr+jXxww8/oG3btujatatC+dKlS7Ft2zaEhYVh2LBhePvtt/HVV189sZ3Q0FCYmJjIvxwdHZ87pobi4WCCg3N6IS70ZfxviEe1deb9fgWtPvwLU3+JRFpeUQNHSERE1HhpKjuAxujRo0fYvHkzPvrooyrbHi/z9vZGQUEBVqxYgXfeeafatkJCQhAcHCx/nZeX1yQSLKBiCPS1Lk54rYsT0vKLsOTP69h3OVmhzoFrKThwrSKJ3TChM3q3toLkOZaEICIiUhWNrufK0tISGhoaSE1NVShPTU2Fra1ttfvY2trWqv6z7NixA4WFhRg3btwz6/r5+eH+/fsoLq5+jSgdHR0YGxsrfNUXY936y5WtjXSxZkxH3PnkZWx9qwsMtKsOO07YcB4uIfsRsvMKMrlmFhERqalGl1xpa2vDx8cH4eHh8jKZTIbw8HD4+/tXu4+/v79CfQAICwt7Yv1n+eGHHzBo0CBYWVk9s+7FixdhZmYGHR2d5zqWmDo4mgIAhno71NsxNKQS+LWwwLWl/XBlcV/MDmhVpc5v5xLh879D6PhxGM7FZdVbLERERI1RoxwWDA4Oxvjx49GpUyf4+vpi9erVKCgowIQJEwAA48aNg4ODA0JDQwEAs2bNQq9evbBy5Ur0798fW7ZswYULF7B+/Xp5m1lZWUhMTERSUhIAIDY2FkBFr9fjPVy3b9/G8ePHq6yTBQB//vknUlNT0aVLF+jq6iIsLAyffPIJ5s6dW2/vRW1Uzqvv2frZSaEYjHS1MDugNWYHtMbVB7kI2XkFVx7kyrdnFZRgxLcVNxW82d0FM/u0gomeVoPERkREpCyNMrkaOXIk0tPTsXDhQqSkpMDLywsHDhyQT1pPTEyEVPpvp1vXrl2xefNmLFiwAB988AFatWqF3bt3w8Pj38nYe/bskSdnADBq1CgAwKJFi7B48WJ5+Y8//ohmzZqhb9++VeLS0tLCmjVrMGfOHAiCAFdXV6xatQqTJ08W+y14LgIqsitlTHnycDDBnzO7o6i0HH9cfIB5v19R2P79yTh8fzIOALD1rS7o5GwOjZo8aZqIiKiJaZTrXKmy+lznasx3Z3D6Tia+GOWFwV71NzRYUxkPi/FV+C38FJFQ7fZXvB2weGA7mOizN4uIiBo3paxzVVpainv37iE2NhZZWZxnQ4CloQ6WDPZA/LL+2PpWF9ga6yps3xX9AB2WHoTz/H3YHf2Aj9ohIiKVUKfkKj8/H2vXrkWvXr1gbGwMZ2dntG3bFlZWVnBycsLkyZNx/vx5sWKlZ6jMTRrjUgh+LSxw5oM+uPV/QXi/n1uV7bO3XoRLyH6M/f4M4jMKlBAhERGROJ47uVq1ahWcnZ2xYcMGBAQEYPfu3bh48SJu3ryJiIgILFq0CGVlZejbty/69euHW7duiRk3VUM+50rJcTyNloYUb/d2Rfyy/jg1/0X4OpsrbD91O1O+CvyGU3F8piERETU5tZpzpaGhgfLyig+70aNHY8GCBWjXrt1T9ykuLsaGDRugra2NiRMn1i1aFVCfc65GfhuBs3FZ+HqMNwa0txe17fpULhMQfiMVb/0SWe325ub6WD3KCx2bmzVwZERERBVq8/ldq+RKKpVCJuPz5OqiPpOrEd9G4FwTTK4el19UipUHb2Lj6fhqt4/2dcR7gW1gbqDdsIEREZFaq7cJ7Y1xLg9VJWnUA4NPZ6SrhcWD2iF+WX/8PbsnXK0VHyD927l76PhxGFw/2I99l5P58GgiImp0GuU6V/Sc5BPalRuGWNxsjXAouBeKy8rx15UUzN56Ub6tTCZg+uYoAEDH5qZY8WoHtLQyfEJLREREDUeU5OrRo0cQBAH6+voAgISEBOzatQvu7u7VLsZJ9aMpTGh/HjqaGhji7YAh3g5Izy/GyoOx2HL+nnx7VGIO+qw8BgCY0M0Zc/u6wUCHfzcQEZFyiLLO1eDBg/Hzzz8DAHJycuDn54eVK1di8ODBWLt2rRiHoBoQVKznqjpWRjpYNqw94pf1x663u8LTwURh+4ZT8Wi36G90+SQcey8nQSbj2llERNSwREmuoqKi0KNHDwDAjh07YGNjg4SEBPz888/48ssvxTgE1YC6pRHezc3w58zuiP1fP3z2ageF5xam5BVhxuZotPigYu2sO+kPuUgpERE1CFHGTgoLC2FkZAQAOHjwIIYOHQqpVIouXbogIaH6R59QfVLhrqtq6GhqYLhPMwz3aYbMh8X4+shtbDgVL99+6nYm+qw8Bk2pBK91ccK7fVvDSJeP3CEiovohSs+Vq6srdu/ejXv37uHvv/+Wz7NKS0sTfbkBerLKnhlVHhZ8FgtDHSwaWHG34d6Z3dHbzUq+rUwmYOPpeHguPoheK45g89lElHPYkIiIRCZKcrVw4ULMnTsXzs7O8PPzg7+/P4CKXixvb28xDkE1UJkmqHFupcDDwQQbJ/ji9v8F4ZuxHdHG1ki+LSGzEB/suoKWH+zHmO/OIDIhm8OGREQkiudeof2/UlJSkJycjA4dOkAqrcjZzp07B2NjY7Rp00acaFVAfS4iOmTNKVy8l4PvxnXCS+42oratKvKLSvHrmUT8cPIuMh6WKGzTkEowrKMDZgW0hoOpnpIiJCKixqg2n9+i3a9ua2sLW1tbhTJfX1+xmqcaYL/LsxnpamFa75aY1rsl7mUV4pujd7Dtwj2UywSUywRsu3Af2y7ch52JLkb7NseEbs6cn0VERLUiWnIVHh6O8PBwpKWlVXlEzo8//ijWYagGOCxYM47m+ggd6olPXvFAVGI21h69i0M3UgEAyblFWBV2E6vCbqJDMxNM6OaClz3toK0pykg6ERGpMFGSqyVLlmDp0qXo1KkT7Ozs+JgcZeGE9ucikUjg42SO78ebo6RMhsMxqfjuRBwiE7IBAJfu52L21ouYvfUiXva0xbRervBsZvKMVomISF2JklytW7cOGzduxOuvvy5Gc/Sc5BPamVw9N21NKfp52KGfhx0KS8qw+WwiNp9LxN30AgDA/isp2H8lBZaGOhjiZY/XujjB2dJAyVETEVFjIkpyVVJSgq5du4rRFNWBfIV2DgyKQl9bE2/2aIE3e7RAen4xvj95F/uvJONe1iNkPCzG9yfj8P3JOLSwMsArXg4Y5+8ME33OzyIiUneiTCB58803sXnzZjGaojoQ5E9uVm4cqsjKSAchQW1x4v0XcSakD97s7gJLQ20AwN30AqwMu4kOSw9i+NrT2Hw2EUWl1d9VS0REqk+UnquioiKsX78ehw4dQvv27aGlpfjX+6pVq8Q4DD0Dl2lqGLYmulgwwB0LBrjj6oNcbDgVjz8vJ6GkTIYLCdm4kJCND3ZdQUBbGwzxtke/drbQ1OBEeCIidSFKcnX58mV4eXkBAK5evaqwjZPbGx7f8Ybj4WCClSM6YMXw9jgck4Zd0Q+w70oyAODQjVQcupEKTakE/dvbYWRnR/i3sODPBBGRihMluTpy5IgYzVAdyedc8cO7wUmlEgS42yDA3Qary2XYdzkZu6If4NjNdJTJBPxxMQl/XEyCkY4mBnvbY4yvE9zt+WgoIiJVVKexioULFyIyMlKsWBSsWbMGzs7O0NXVhZ+fH86dO/fU+tu3b0ebNm2gq6sLT09P7N+/X2H7zp070bdvX1hYVPQcXLx4sUobvXv3hkQiUfiaOnWqQp3ExET0798f+vr6sLa2xnvvvYeysrI6n68Y+PibxkFLQ4oh3g74aaIvbizth6WD28G7uSkAIL+4DL+eScTLX55A19BwLPsrBrdS85UbMBERiapOydX9+/cRFBSEZs2aYdq0afjrr79QUlLy7B2fYevWrQgODsaiRYsQFRWFDh06IDAwEGlpadXWP336NEaPHo1JkyYhOjoaQ4YMwZAhQxSGKAsKCtC9e3d8+umnTz325MmTkZycLP9avny5fFt5eTn69++PkpISnD59Gj/99BM2btyIhQsX1vmcxcAHNzc+etoaGOfvjF1vd8P5DwMwr18btPunxyoptwjrjt3BS58fx4srj2JV2E0kZhYqOWIiIqqrOj9bUCaT4dSpU/jzzz/xxx9/IDk5GS+99BIGDx6MAQMGwNzcvNZB+fn5oXPnzvj666/lx3B0dMTMmTMxf/78KvVHjhyJgoIC7N27V17WpUsXeHl5Yd26dQp14+Pj4eLigujoaPk8sUq9e/eGl5cXVq9eXW1cf/31FwYMGICkpCTY2FQ8u2/dunWYN28e0tPToa2tXWWf4uJiFBcXy1/n5eXB0dGxXp4t2G/1ccSk5OOXSb7o0cpK1LZJXMm5j7Dt/H3sufQAd/5ZQ6uSu50xBnawx6udmsHSUEdJERIR0eNq82zBOt/CJJVK0aNHDyxfvhyxsbE4e/Ys/Pz88O2338Le3h49e/bEZ599hgcPHtSovZKSEkRGRiIgIEDhGAEBAYiIiKh2n4iICIX6ABAYGPjE+k+zadMmWFpawsPDAyEhISgs/LcnISIiAp6envLEqvI4eXl5uHbtWrXthYaGwsTERP7l6OhY65hqi+tcNX52JnqYFdAK4e/2xrH3emNKzxawN9EFAFxPzsOnB2LQ6X+HMGztaWw8FYfcR6VKjpiIiGpKtGcLVmrbti3atm2L999/H+np6dizZw/27NkDAJg7d+4z98/IyEB5eblCAgMANjY2iImJqXaflJSUauunpKTUKvYxY8bAyckJ9vb2uHz5MubNm4fY2Fjs3Lnzqcep3FadkJAQBAcHy19X9lzVh38ntNdL81RPnCwMEPJyW4S83BaX7+dgZ1TFHYfp+cWITMhGZEI2Fv95HV1amKN/e3sMam/PxUqJiBqxOiVXb731FlatWgVDQ8Nqt1tZWWHSpEmYNGlSXQ7TYN566y35/z09PWFnZ4c+ffrgzp07aNmy5XO1qaOjAx2dhhnaqVxElLlV09W+mSnaNzPFooHuuJCQjT0Xk7DvSjKyCkpw5m4WztzNwke7r6K7qyUGtLfDwA72MNAR/W8kIiKqgzoNC/7www8Kw2ZisLS0hIaGBlJTUxXKU1NTYWtrW+0+tra2tapfU35+fgCA27dvP/U4lduUTeDtgipDIpGgs7M5Ph7igcgFAdg2xR+jOjvC8J9E6uTtDMzfeQXtFv2NiRvPY0fkfTwq4arwRESNQZ2Sq1rMha8xbW1t+Pj4IDw8XF4mk8kQHh4Of3//avfx9/dXqA8AYWFhT6xfU5XLNdjZ2cmPc+XKFYW7FsPCwmBsbAx3d/c6HUsMXKBdNUkkEvi6mGPZsPa4srgvNr3ph+E+zWD0T6J1OCYNc7dfgvuiA5i08Tx+j7yPwpLGsTwIEZE6apTjCcHBwRg/fjw6deoEX19frF69GgUFBZgwYQIAYNy4cXBwcEBoaCgAYNasWejVqxdWrlyJ/v37Y8uWLbhw4QLWr18vbzMrKwuJiYlISkoCAMTGxgKo6HGytbXFnTt3sHnzZrz88suwsLDA5cuXMWfOHPTs2RPt27cHAPTt2xfu7u54/fXXsXz5cqSkpGDBggWYPn16gw391QQntKsuiUSCbq6W6OZqiU+Htcfxm+nYdyUZB66m4GFxGcJj0hAek4a5Oy7hRTdr9G9vh34ettDXbpQ/6kREKqlOv3ErF9oU28iRI5Geno6FCxciJSUFXl5eOHDggHzyeGJiIqTSfzvdunbtis2bN2PBggX44IMP0KpVK+zevRseHh7yOnv27JEnZwAwatQoAMCiRYuwePFiaGtr49ChQ/JEztHREcOGDcOCBQvk+2hoaGDv3r2YNm0a/P39YWBggPHjx2Pp0qWivwfPg+tcqRcNqQQvtLHGC22sKxKtW+n482IS/rqagkel5f8mWtsv4cU2Nujf3haB7ZhoERHVtzqtcyWVShEUFPTMXpvKu+2odutk1NaLK4/ibnoBtr7VBX4tLERtm5qOcpmA4zfT8eelJBy8noqHxf8OEUolwAtu1gjytENgOxsY6fKuQyKimqjN53ed/oQdP358XXYnsfHZggTFHq2ychlO3MrA/ivJOHAtBflF/w4dvrcD6NnKCv097RDoYQsTPSZaRERiqPMK7VQ79dlz9cJnRxGXUYBtU/zh61L7lfFJtZXLBJy6nYG9lyuGDvOL/u3RkkiA7q6WCPKwQ//2dky0iIj+o956rurj7kASHzuuqDoaUgl6trZCz9ZW+HRYe0TcycTeK8n4+2oKMgtKcOJWBk7cysAHu66gSwtzBHlUTIa3MdZVduhERE1KrXqurly5gj179kBPTw99+vRBhw4d6jM2lVSfPVe9VxxBfGYhdkz1Rydn9lxRzchkAs7FZ+HA1RT8dTUZqXnFCtu9m5vi5X96tOxN9ZQUJRGRctXm87tWyVWlkpIShIeHIzo6GpqamujWrRu6dev23AGrk/pMrnqtOIKEzEL8Ps0fPk5Mrqj2BEHApfu52HspCQeupeB+9iOF7R4OxujXzhb9POzgal39kxmIiFRRvSdXjxMEASdPnsSpU6cgCAI6duyIPn36QFOTt3tXpz6Tq57LjyAxqxC/T+sKHyczUdsm9XQjOQ/7Lifjr6vJuJNeoLCtpZUB+nnYIsjDDh4OJkqKkIioYTRocvVf0dHRCA8PR0lJCdq0aYOgoCDo6XEooVJ9Jlc9lh/GvaxHTK6oXtxOe4i/r6XgwNUUXHmQq7DN3kQXfdvZ4mVPO3RyMoNUyol/RKRaGmwphup4e3vD29sbQEWiFRQUBDMzM4SGhqJNmzZiH44eI8iXYlBuHKSaXK0N4WrtiukvuCI1r6hiHa1rqTgXn4Wk3CJsPB2PjafjYW6gjZfa2iDQwwbdXa2grVmnp2wRETU5oidXeXl52LNnD3bs2IG7d+8iKCgIw4YNY2LVgJhbUX2zMdbFmz1a4M0eLZDxsBhh11Px97UUHL+ZjqyCEmy9cA9bL9yDvrYGXnCzRt92NnixjTUXLSUitSBaclVYWIjhw4cjOTkZgwcPxtKlS+XP5KOGIXARUVICS0MdjPZtjtG+zfGopBxhN1Jx8FoKjsSkoaCkHPuuJGPflWRoSiXwb2mBvu1s0a+dLayMGs/zOImIxCRacqWvr48+ffrg3XffFatJek5MrUhZ9LQ1MKiDPQZ1sEdJmQynbmfg4PUU/H0tFVmPraX10e6r6OBoir7uNghsZwNXayNlh05EJBrRJ7TT09XnhPauoeFIyi3CH9O7oYOjqahtE9WFIAiISszB39dScPBaCuIzCxW2u1ga4KV/Eq2Ozc3Y+0pEjY5SJ7ST8lRmyfxcosZGIpHAx8kMPk5m+ODltriVmo+D11MRdj0VF+/lIC6jAOuP38X643dhpq+Fl9xt8JK7Lbq7WkJPW0PZ4RMR1cpzJ1eJiYlo3rx5jes/ePAADg4Oz3s4qgUJBwapkWtlY4RWNkaY/oIr0vOL8fe1FIRdT8XpOxnILizFtgv3se3CfehoStGjlSUC2tqgT1sbztMioibhue+R7ty5M6ZMmYLz588/sU5ubi6+++47eHh44Pfff3/eQ1ENcSkGaoqsjHTwWhcn/DTRFxcX9sU3YzviFW8HmOhpobhMhkM30jB/5xV0/r9DeOWbU1hz5DZup+UrO2wioid67p6r69ev4//+7//w0ksvQVdXFz4+PrC3t4euri6ys7Nx/fp1XLt2DR07dsTy5cvx8ssvixk3VUMAp89R02ago4mXPe3wsqcdyspluJCQjbDrqTh0IxUJmYWITsxBdGIOVvwdC0dzPQS0tUFAWxv4uphDS4PraRFR41DnCe2PHj3Cvn37cPLkSSQkJODRo0ewtLSEt7c3AgMD4eHhIVasKqE+J7T7/t8hpOUXY+/M7nwcCamcO+kVK8SH30hDVGI2Hv/NZaSjiV5uVghoa4MX3Kxhos/1tIhIXEp9/A09XX0mV53/7xDS84ux753uaGfP5IpUV+bDYhyOScOhG6k4fjMDj0rL5dukEqCTszn6tLHGS+42aGHFB0wTUd3xbkE1xwntpOosDHXwaidHvNrJEcVl5Yi4k4nwG2kIv5GKpNwinIvLwrm4LIT+FQMnC3282MYafdpUDB/ycTxEVN+YXKkQTmgndaSjqYHebtbo7WaNj4d4IDYlH2HXUxAek4aL93KQkFmIDafiseFUPAy0NdCztRVebFNRn3cfElF9YHKlUiqyKyZXpM7cbI3gZmuEGS+2QsbDYhyJScPhmDScuJWBh8Vl+OtqCv66mgIA6OBoihfdrNGnrTXc7YwhlfKHh4jqjsmVCuHsOSJFlo8NH5aUyXA+PgvhN9JwJDYNcRkFuHQvB5fu5eDzQzdhZaSDF9wqerV6traCvjZ/PRLR8xFl8sGjR49QWPjv4ywSEhKwevVqHDx48LnbXLNmDZydnaGrqws/Pz+cO3fuqfW3b9+ONm3aQFdXF56enti/f7/C9p07d6Jv376wsLCARCLBxYsXFbZnZWVh5syZcHNzg56eHpo3b4533nkHubm5CvUkEkmVry1btjz3eYpJvkI751wRVaGtKUU3V0ssHOiOI3N74+jc3ljQvy26uVpAQypBen4xtl24j6m/RsFj0d8Y890ZfH/iLu6kP1R26ETUxIjyp9ngwYMxdOhQTJ06FTk5OfDz84OWlhYyMjKwatUqTJs2rVbtbd26FcHBwVi3bh38/PywevVqBAYGIjY2FtbW1lXqnz59GqNHj0ZoaCgGDBiAzZs3Y8iQIYiKipIvBVFQUIDu3btjxIgRmDx5cpU2kpKSkJSUhM8++wzu7u5ISEjA1KlTkZSUhB07dijU3bBhA/r16yd/bWpqWqvzqy+VN35yWJDo2ZwtDfBmjxZ4s0cLFBSX4cStdByOScPR2HSk5Rfj9J1MnL6Tif/tuwFHcz28+M+8ri4tLPhIHiJ6KlGWYrC0tMSxY8fQrl07fP/99/jqq68QHR2N33//HQsXLsSNGzdq1Z6fnx86d+6Mr7/+GgAgk8ng6OiImTNnYv78+VXqjxw5EgUFBdi7d6+8rEuXLvDy8sK6desU6sbHx8PFxQXR0dHw8vJ6ahzbt2/Ha6+9hoKCAmhqVuShEokEu3btwpAhQ2p1TpXqcykG76UHkV1YirA5PdHKxkjUtonUhUwm4HpyHo7EVAwfRt/LURhy19GUoksLC7zgZoXebtZwtjRQXrBE1GAafCmGwsJCGBlVfJgfPHgQQ4cOhVQqRZcuXZCQkFCrtkpKShAZGYmQkBB5mVQqRUBAACIiIqrdJyIiAsHBwQplgYGB2L17d+1O5D8q38DKxKrS9OnT8eabb6JFixaYOnUqJkyYAMkTuouKi4tRXFwsf52Xl1enmJ6GD24mqjupVAIPBxN4OJhgZp9WyC4owfFb6TgSk4ZjN9ORXViKYzfTcexmOvDndbhYGqBXayv0drNClxYW0NVirxaRuhMluXJ1dcXu3bvxyiuv4O+//8acOXMAAGlpabXuncnIyEB5eTlsbGwUym1sbBATE1PtPikpKdXWT0lJqdWx/xvHxx9/jLfeekuhfOnSpXjxxRehr6+PgwcP4u2338bDhw/xzjvvVNtOaGgolixZ8txx1AYntBOJz8xAG4O9HDDYywHlMgGX7ufgaGw6jsWm4fKDXMRlFCAuowAbT8dDV0sKPxcL9K7s1bLQf+IfXkSkukRJrhYuXIgxY8Zgzpw5ePHFF+Hv7w+gohfL29tbjEM0qLy8PPTv3x/u7u5YvHixwraPPvpI/n9vb28UFBRgxYoVT0yuQkJCFHrV8vLy4OjoWC9x/zvCy1/mRPVBQypBx+Zm6NjcDMEvtUbGw2Icv5mOo7HpOH4rHTmP9Wot+fM6mpvro1drK/RqbYWurha8A5FITYjykz58+HB0794dycnJ6NChg7y8T58+eOWVV2rVlqWlJTQ0NJCamqpQnpqaCltb22r3sbW1rVX9p8nPz0e/fv1gZGSEXbt2QUvr6c8o8/Pzw8cff4zi4mLo6FRdkFBHR6fa8vrAYUGihmVpqIOhHZthaMdmkMkEXLyfg2Ox6Th6Mx2X7uUgMasQv5xJwC9nEqCtIUVnFzP0bGWFXm5WcLMxYq8WkYoS7c8oW1tb2NraQhAECIIAiUQCX1/fWrejra0NHx8fhIeHyyeNy2QyhIeHY8aMGdXu4+/vj/DwcMyePVteFhYWJu9Bq6m8vDwEBgZCR0cHe/bsga6u7jP3uXjxIszMzBosgaoJ/romanjSx3q15rzUGjmFJTh2Mx3Hb2bg2M00ZDwswanbmTh1OxOhf8XAxlgHPVtZoWdrK/RoZQlTfW1lnwIRiUS05OqHH37A559/jlu3bgEAWrVqhdmzZ+PNN9+sdVvBwcEYP348OnXqBF9fX6xevRoFBQWYMGECAGDcuHFwcHBAaGgoAGDWrFno1asXVq5cif79+2PLli24cOEC1q9fL28zKysLiYmJSEpKAgDExsYC+DcpzMvLQ9++fVFYWIhff/0VeXl58snnVlZW0NDQwJ9//onU1FR06dIFurq6CAsLwyeffIK5c+c+/xsnJvnjb5heESmbqf6/c7UEoeIOxGM303EsNh3n47OQmleM7ZH3sT3yPiQSoH0zU/RqZYmera3g5WgKTQ0+A5GoqRJtztWqVaswc+ZMeW9RREQE5syZg8TERCxdurRW7Y0cORLp6elYuHAhUlJS4OXlhQMHDsgnrScmJkIq/fcXT9euXbF582YsWLAAH3zwAVq1aoXdu3fL17gCgD179siTMwAYNWoUAGDRokVYvHgxoqKicPbsWQAVE/QfFxcXB2dnZ2hpaWHNmjWYM2cOBEGAq6srVq1aVe26WcrA+exEjZNEIkE7exO0szfB271dUVBchlO3Myp6tm6l417WI/lq8V8evg0jXU10a2mJHq0t0bOVFRzN9ZV9CkRUC6Ksc2VlZYUvv/wSo0ePVij/7bffMHPmTGRkZNT1ECqjPte5arfwAApKynF0bm+uvUPURAiCgITMQhy/lY7jN9Nx6nYmHpWWK9RxttBHj1YVw4f+LS1gpPv0uaBEJL4GX+eqtLQUnTp1qlLu4+ODsrIyMQ5BNcAJ7URNj0QigbOlAZwtDTDO3xll5TJcSMjG8ZvpOHErA1ce5CI+sxDxmRUT4yvuWDRFd1cr9GhtCU8HE2hxCJGoURGl52rmzJnQ0tLCqlWrFMrnzp2LR48eYc2aNXU9hMqoz54r94UHUFhSjuPvvYDmFhxGIFIFlYuYnryVgRO3MpCSV6Sw3VhXE/4tLdCjlRW6u1rCiWtrEdWLBum5enztJolEgu+//x4HDx5Ely5dAABnz55FYmIixo0b97yHoFoS5BPalRsHEYnn8UVMAeB2Wj6O38zAiVvpiLibibyiMvx9LRV/X6tYjsbRXA/dXS3R3dUKXVtawMyAdyESNbTn7rl64YUXanYAiQSHDx9+nkOopPrsuWrz0V8oKpXhxPsvcAIskRooLZfhQnw2Tt6u6Nm68iAXssd+o0slQDt7E3RztUR3V0t0cjbj43mInlNtPr9FGRakmqvP5MptwV8oLmNyRaSucgpLEHEnEyduZ+DU7QwkZBYqbNfWlKKzs5k82fJ0MOEQIlENNfiE9krXr19HYmIiSkpK5GUSiQQDBw4U8zD0BJzQTqTeTPW1EeRphyBPOwDAvaxCnLydgZO3MhBxNxNZBf8uZLocsTDW1UQ3V0t0bWmBbq6WcLE0YLJFJAJRkqu7d+/ilVdewZUrVyCRSOTPuKv8IS0vL3/a7iQy/nIkIgBwNNfHaN/mGO3bHOUyATeS83DqdgZO3s7A+fgs5BWV4a+rKfjrasVD7u1MdNG15b/Jlq3Js59SQURViZJczZo1Cy4uLggPD4eLiwvOnTuHzMxMvPvuu/jss8/EOATVROWEduVGQUSNkIZUAg8HE3g4mGBKr5YoLitHVEIOTt+pGEKMvpeD5Nwi/B51H79H3QcAtLAyqEi0WlqiSwtOjieqKVGSq4iICBw+fBiWlpaQSqWQSqXo3r07QkND8c477yA6OlqMw9AzCKjsMVRyIETU6OloasC/pQX8W1rg3b5ueFhchjN3MnHqTgZO385EbGo+7qYX4G56AX49kwiJBHC3M4Z/i4peLV8XcxjoiDqzhEhliPKTUV5eDiMjIwCApaUlkpKS4ObmBicnJ/kz/Kj+8dYEInpehjqaCHC3QYB7xWPGMh4WI+JOJk7fycDpO5lIyCzEtaQ8XEvKw/cn4yCVAB0cTdGtZcWq8T5OvBORqJIoyZWHhwcuXboEFxcX+Pn5Yfny5dDW1sb69evRokULMQ5BNSCf0M6BQSKqI0tDHQzsYI+BHewBAPezC3H6TiYi/vlKyStCdGIOohNz8PWR29DWlMLb0RRdWliga0sLeDU3hY4mky1ST6IkVwsWLEBBQQEAYOnSpRgwYAB69OgBCwsLbN26VYxDUA38eyOBkgMhIpXTzEwfIzrpY0QnRwiCgLiMgopk625FspVVUIKzcVk4G5eFL8JvQUdTCh8nM/i3qBh67OBoysf0kNqot3WusrKyYGZmxjvX/qM+17lqEbIPMgE490EfWBvzLh8iahiCIOBm6kOcvpOBiDuZOBuXhdxHpQp1dLWk6Oxsji4tLNClhQU8HUygrclki5oOpa1z9Thzc/P6apqeQJ4lM58logYkkUjgZmsEN1sjTOjmIl/24czdTJy5W5Fs5ReV4cQ/z0cEAANtDXR0Mvsn2TJHh2am0GTPFqkI3uqhQjihnYgag8eXfXizRwuUywRcT/o32ToXl4X8YsVkq7Jny8/FHH7/9Gxxgjw1VUyuVBAntBNRY6IhlcCzmQk8m5lgcs+KZOtaUu4/yVYWzsdX7dnS09KAd3NT+LqYw8/FAh2dOEGemg4mVyri8alznOZGRI2ZhlSC9s1M0b6ZKd7q2RJl5TLEpOTLk61zcZnIKyrD6TuZOH0nE8AtaGlI4O1oBr8W5vB1MUfH5mZcZ4saLX5nqiDmVkTUlGhqSBWGEWUyATdS8nD2bhbOxmXifHw2sgpKcC4+C+fiswAAWhoStLM3gZ9LRbLVydkcJnpaSj4TogpMrlTE4/OteIcmETVlUmlF4tTO3gQTu7tAJhNwK+0hzsVVTI4/F5eFtPxiXLyXg4v3cvDt8bsAgLZ2xvB1NoNfCwt0cjaDtRHvmiblYHKlIjiXnYhUlVT6792Ir/s7AwDiMgoqkq27FWtrPch5hBvJebiRnIefIhIAVDwbsbOTOTq7mMPX2RyO5nr845MaBJMrFaEw50qJcRARNQQXSwO4WBpgZOfmEAQBSblFOB9XMWx4Pi4Lt9Mfyp+NuPXCPQCAjbEOOjtXDCN2djZHG1sjJltUL5hcqYjHe674u4KI1IlEIoGDqR4cvB0wxNsBAJBVUILz/yRa5+KzcPVBLlLzirH3cjL2Xk4GABjpaqKTkxk6u5ijk5M52jfj8g8kDiZXKohLMRCRujM30EZgO1sEtrMFABSWlCEyIRvn47JwPj4bF+/lIL+oDEdi03EkNh0AoK0pRXsHE/g4m8HX2Rw+TmYw1ddW5mlQE9Vol8Nds2YNnJ2doaurCz8/P5w7d+6p9bdv3442bdpAV1cXnp6e2L9/v8L2nTt3om/fvrCwsIBEIsHFixertFFUVITp06fDwsIChoaGGDZsGFJTUxXqJCYmon///tDX14e1tTXee+89lJWV1fl860phAVHmVkRECvS1NdGjlRWC+7rht7e64NKivtg9vRs+fLktXnK3gbmBNkrKZLiQkI1vj93FpJ8uwGtpGF5adQwhO6/g98j7SMgsQD09MY5UTKPsudq6dSuCg4Oxbt06+Pn5YfXq1QgMDERsbCysra2r1D99+jRGjx6N0NBQDBgwAJs3b8aQIUMQFRUFDw8PAEBBQQG6d++OESNGYPLkydUed86cOdi3bx+2b98OExMTzJgxA0OHDsWpU6cAAOXl5ejfvz9sbW1x+vRpJCcnY9y4cdDS0sInn3xSf29IDQic0k5EVGPamlJ4OZrCy9EUk3tWLP9wN6MAF+IrerYuJGQhIbMQt9Ie4lbaQ/x2LhEAYGWkg05OZvBxMkMnZ3N42BvzsT1URb09uLku/Pz80LlzZ3z99dcAAJlMBkdHR8ycORPz58+vUn/kyJEoKCjA3r175WVdunSBl5cX1q1bp1A3Pj4eLi4uiI6OhpeXl7w8NzcXVlZW2Lx5M4YPHw4AiImJQdu2bREREYEuXbrgr7/+woABA5CUlAQbGxsAwLp16zBv3jykp6dDW7tq93FxcTGKi4vlr/Py8uDo6Cj6g5uLSsvR5qMDAIDLi/vCWJfrvRAR1UVafhEi47NxISEbF+KzcOVBLmT/+cTU1pTC29EUnZwrEi5vRzOYGXAoURU1igc3P6+SkhJERkYiJCREXiaVShEQEICIiIhq94mIiEBwcLBCWWBgIHbv3l3j40ZGRqK0tBQBAQHysjZt2qB58+by5CoiIgKenp7yxKryONOmTcO1a9fg7e1dpd3Q0FAsWbKkxnGIgaOCRER1Z22kiyBPOwR52gGo+CM2KjFbnnBFJ2Yjr6gMZ+MqloOo1NLKAD7/9G75OJmhpZUh70pUM40uucrIyEB5eblCAgMANjY2iImJqXaflJSUauunpKTU+LgpKSnQ1taGqanpE9t50nEqt1UnJCREIfGr7LkSGxcRJSKqX7paGuja0hJdW1oCAMplAm6nPURkQsUwYmRCNhIyC3EnvQB30guw7cJ9AICpvhY6Nv+nZ6t5xWN/DPnoHpXGq1vPdHR0oKOj06DHZGpFRFT/NB5b3HSMX3MAQMbDYlyIz67o4UrIxtUHucgpLMXhmDQcjkkDAGhKJWhtYwQfJzN0dDKFp4MpWloZ8A9jFdLokitLS0toaGhUuUsvNTUVtra21e5ja2tbq/pPaqOkpAQ5OTkKvVePt2Nra1vlrsXK49bmWPWBE9qJiJTP0lAH/Txs0c+j4jOhqLQc15LyEP1PshWVmI3UvGJcT87D9eQ8/HKmYjV5KyMddGhmCl8XM3g4mKBjczOuudWENbrkSltbGz4+PggPD8eQIUMAVExoDw8Px4wZM6rdx9/fH+Hh4Zg9e7a8LCwsDP7+/jU+ro+PD7S0tBAeHo5hw4YBAGJjY5GYmChvx9/fH//3f/+HtLQ0+V2LYWFhMDY2hru7+3OcrXgUhwWVFwcREf1LV0tDPvfqzR4VZfezCxGZkI3oxBxcvp+Da0l5SM8vxqEbqTh0o+IPdi0NCVytjeDnUrG4qY+TGRzN9CGV8hd8U9DokisACA4Oxvjx49GpUyf4+vpi9erVKCgowIQJEwAA48aNg4ODA0JDQwEAs2bNQq9evbBy5Ur0798fW7ZswYULF7B+/Xp5m1lZWUhMTERSUhKAisQJqOhxsrW1hYmJCSZNmoTg4GCYm5vD2NgYM2fOhL+/P7p06QIA6Nu3L9zd3fH6669j+fLlSElJwYIFCzB9+vQGH/r7L8VlrvjDR0TUWDUz00czM30M9qpYTb64rBzRiTm4+iAXZ+OycPl+DlLziuXPSqxkYaCN9s1M0MnZHJ4OFQmXAeduNUqN8qqMHDkS6enpWLhwIVJSUuDl5YUDBw7IJ48nJiZCKv13XZGuXbti8+bNWLBgAT744AO0atUKu3fvlq9xBQB79uyRJ2cAMGrUKADAokWLsHjxYgDA559/DqlUimHDhqG4uBiBgYH45ptv5PtoaGhg7969mDZtGvz9/WFgYIDx48dj6dKl9fl21IjCswWZWxERNRk6mhro0sICXVpY4M0eLSAIAu5nP0L0vRxcTMxBZEIWbiTnI7OgRGFFeamk4hmLvi7m8HQwha+LOVwsDaDB3i2la5TrXKmy2qyTURv5RaXwXHwQABDzcT+O1RMRqZCSMhku3a/o3Tofn4XL93NxP/tRlXo6mlL5XYnudibwa2EOS0Pljqyoiia9zhU9H2bIRESqS1tTis7O5ujsbI4J3VwAAGl5RRW9W/dyEJmQjSv3c/GotByn72Ti9J1M+b52JrrwdDCBr4s53O2M0d6RS0HUN767KoIT2omI1Iu1sa7Cw6krHuHzEBF3s3A9KRcX7+XiRnIeknOLkJxbhIPXKybLa0glaGNrhLZ2xvB1NodnMxO0tDKEtiYf4yMWJleq4vHkihPaiYjUjlRacYehq7WRvCz3USliU/JxISEL0Yk5uJ6Uhwc5j3AtKQ/XkvKwI7JioVMdTal8KLGDY8WkeXsTXa699ZyYXKmIx9e54s8CEREBgImeFnxdzOHrYi4ve5DzCFcf5OLy/Rycj8tGTEoe8orKcOZuFs7c/fcxPka6mujkZIb2zUwrerlczGHO5ybWCJMrFcTcioiInsTBVA8Opnry4URBEHAnveIxPjeSK3q5YpLzkV9UpnB3IvDv/K1OzmZwszWGl6MpTPS0lHUqjRbvFmxg9XW3oEwmIKuwBEDFWijsyiUioudVUibDzdR8nIvLwo3kPFy6n4ObqQ+r1JNIgBaWBmhjZ4zOTmZo52CCNrZGMNJVvYSrNp/fTK4aWH0lV0RERPUp91Epbqbm40J8Ni7fz8GN5DzEZxZWqSeVAG62xmhjawQvR1N4NzeFq7Uh9LWb9mAZk6tGjMkVERGpirT8IsQk5+PKg1xciM9CbEo+knKLqtSTSoC2dsbwdDBBWztj+DiZoZWNIXQ0m86ajEyuGjEmV0REpMpS84pw8V4OYv6Zv3UjOQ8ZD0uq1Kvs4fJubgo3GyN4N6+YOK+l0TiXhGBy1YgxuSIiInWTkluEyIRsxKbk4dL9XEQnZiOvqKxKPU2pBC2tDNHRqSLRamdvDA8Hk0bRw8XkqhFjckVEROpOEASk5hUjMiEb15IqFjuNTHhywuViaYD2zUzRwdEErW2M0L6ZSYPP4WJy1YgxuSIiIqpKEAQk5xbhyoNcRCfm4FZqPi7ey0FmQdUhRQ2pBM4W+mhjZ4yOzc3gZmMEDwdjmOrX3zpcTK4aMSZXRERENSMIAh7kPMLN1HxEJeQgJqViZfnkaibNA4CzhT5a2Rihu6slRnRyhJ62eMOJTK4aMSZXREREdZOSW4Sbqfm4fD8Hl+7n4lZqvsKyEIY6moj66CVRn5dYm8/vpr3oBBEREakdWxNd2JroomdrK3lZVkEJYpLzcPlBLh6VlCv1QdRMroiIiKjJMzfQRldXS3R1tVR2KGici0kQERERNVFMroiIiIhExOSKiIiISERMroiIiIhExAntDaxy5Yu8vDwlR0JEREQ1Vfm5XZMVrJhcNbD8/HwAgKOjo5IjISIiotrKz8+HiYnJU+twEdEGJpPJkJSUBCMjI0gkElHbzsvLg6OjI+7du6eSC5Sq+vkBqn+OPL+mT9XPUdXPD1D9c6yv8xMEAfn5+bC3t4dU+vRZVey5amBSqRTNmjWr12MYGxur5A9MJVU/P0D1z5Hn1/Sp+jmq+vkBqn+O9XF+z+qxqsQJ7UREREQiYnJFREREJCImVypER0cHixYtgo6OjrJDqReqfn6A6p8jz6/pU/VzVPXzA1T/HBvD+XFCOxEREZGI2HNFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQi0lR2AOpGJpMhKSkJRkZGkEgkyg6HiIiIakAQBOTn58Pe3h5S6dP7pphcNbCkpCQ4OjoqOwwiIiJ6Dvfu3UOzZs2eWofJVQMzMjICUHFxjI2NlRwNERER1UReXh4cHR3ln+NPw+SqgVUOBRobGzO5IiIiamJqMqWHE9qJiIiIRMTkioiIiEhETK6IiIiIRMQ5V0REYsm8Dpz9BCgt/M8GoYYNPG0uR2Npo6bt1DWOuuxfX8dWVtxP2+95Y1XFY/2zj0QCaBkCQT89I676w+SKiEgs0V8DNzYpOwoi0jVX6uGZXBERiaWsqOJf1yGAc786Nibg2T1I6tJGXfZX1r512b8h96vvYz3W0yQIFb1KgOK+j999JzzeMyU8eb9q9xH+/b+GTg1iqz9MroiIxGbXBegwRdlREJGScEI7EZFoKv+C5qOtiNQZkysiItHUdMI4EakyJldERGKpnO/Bh7ITqTUmV0REouGwIBExuSIiEh97rojUGpMrIiKxCJxzRURMroiIRMRhQSJickVEJB5OaCciMLkiIqoHTK6I1BmTKyIi0bDnioiYXBERiYcT2okITK5qLTQ0FJ07d4aRkRGsra0xZMgQxMbGKjssImoUOKGdiJhc1dqxY8cwffp0nDlzBmFhYSgtLUXfvn1RUFCg7NCIqNFgckWkzjSVHUBTc+DAAYXXGzduhLW1NSIjI9GzZ08lRUVEjQLvFiQiMLmqs9zcXACAubl5tduLi4tRXFwsf52Xl9cgcRGRMnDOFRFxWLBOZDIZZs+ejW7dusHDw6PaOqGhoTAxMZF/OTo6NnCURNRwOOeKiJhc1cn06dNx9epVbNmy5Yl1QkJCkJubK/+6d+9eA0ZIRErBYUEitcZhwec0Y8YM7N27F8ePH0ezZs2eWE9HRwc6OjoNGBkRKY3AnisiYnJVa4IgYObMmdi1axeOHj0KFxcXZYdERI0G51wREZOrWps+fTo2b96MP/74A0ZGRkhJSQEAmJiYQE9PT8nREZFS8W5BIgLnXNXa2rVrkZubi969e8POzk7+tXXrVmWHRkSNBpMrInXGnqtaEvh4CyJ6IvZcERF7roiIxMM/vogITK6IiETEuwWJiMkVEVE9YHJFpM6YXBERiYV3CxIRmFwREYmIc66IiMkVEZGIOOeKiJhcERGJj8OCRGqNyRURkVj4bEEiApMrIiIRcc4VETG5IiISD+8WJCIwuSIiEhGHBYmIyRURkfjYc0Wk1phcERGJhc8WJCIwuSIiEhGHBYkI0FR2AA2ltLQUKSkpKCwshJWVFczNzZUdEhGpGk5oJyKoeM9Vfn4+1q5di169esHY2BjOzs5o27YtrKys4OTkhMmTJ+P8+fPKDpOIVA6TKyJ1prLJ1apVq+Ds7IwNGzYgICAAu3fvxsWLF3Hz5k1ERERg0aJFKCsrQ9++fdGvXz/cunVL2SETUZPHYUEiUuFhwfPnz+P48eNo165dtdt9fX0xceJErFu3Dhs2bMCJEyfQqlWrBo6SiFQLJ7QTkQonV7/99luN6uno6GDq1Kn1HA0RqQXOuSIiqPCwIBGR8jC5IlJnKttzNXTo0BrX3blzZz1GQkTqgz1XRKTCPVcmJibyL2NjY4SHh+PChQvy7ZGRkQgPD4eJiYkSoyQilcJFRIkIKtxztWHDBvn/582bhxEjRmDdunXQ0NAAAJSXl+Ptt9+GsbGxskIkIpXDuwWJSIV7rh73448/Yu7cufLECgA0NDQQHByMH3/8UYmREZFK4rAgkVpTi+SqrKwMMTExVcpjYmIgk8mUEBERqSSBPVdEpMLDgo+bMGECJk2ahDt37sDX1xcAcPbsWSxbtgwTJkxQcnREpDo454qI1CS5+uyzz2Bra4uVK1ciOTkZAGBnZ4f33nsP7777rpKjIyKVwXWuiAhqklxJpVK8//77eP/995GXlwcAnMhORPWIyRWROlOL5OpxTKqIqP5wzhURqVFytWPHDmzbtg2JiYkoKSlR2BYVFaWkqIhIpXCdKyKCmtwt+OWXX2LChAmwsbFBdHQ0fH19YWFhgbt37yIoKEjZ4RGRyuCcKyJSk+Tqm2++wfr16/HVV19BW1sb77//PsLCwvDOO+8gNzdX2eERkcphckWkztQiuUpMTETXrl0BAHp6esjPzwcAvP766/jtt9+UGRoRqRT2XBGRmiRXtra2yMrKAgA0b94cZ86cAQDExcVB4BwJIhILf58QEdQkuXrxxRexZ88eABULis6ZMwcvvfQSRo4ciVdeeUXJ0RGR6uDdgkSkJncLrl+/Xv6Ym+nTp8PCwgKnT5/GoEGDMGXKFCVHR0Qqh8OCRGpN5ZOrsrIyfPLJJ5g4cSKaNWsGABg1ahRGjRql5MiISOXw2YJEBDUYFtTU1MTy5ctRVlam7FCISOVxzhURqUFyBQB9+vTBsWPHRGnr+PHjGDhwIOzt7SGRSLB7925R2iUiFcBnCxIR1GBYEACCgoIwf/58XLlyBT4+PjAwMFDYPmjQoBq3VVBQgA4dOmDixIkYOnSo2KESUZPGYUEiUpPk6u233wYArFq1qso2iUSC8vLyGrcVFBTEVd2J6OnYc0Wk1tQiuaq8U1AZiouLUVxcLH+dl5entFiIqJ5xnSsigprMuVKm0NBQmJiYyL8cHR2VHRIR1RsOCxKRCidXW7ZsqXHde/fu4dSpU/USR0hICHJzc+Vf9+7dq5fjEFFjwOSKiFQ4uVq7di3atm2L5cuX48aNG1W25+bmYv/+/RgzZgw6duyIzMzMeolDR0cHxsbGCl9EpOI454pIransnKtjx45hz549+OqrrxASEgIDAwPY2NhAV1cX2dnZSElJgaWlJd544w1cvXoVNjY2yg6ZiJo6LiJKRFDh5AqoWGJh0KBByMjIwMmTJ5GQkIBHjx7B0tIS3t7e8Pb2hlRau867hw8f4vbt2/LXcXFxuHjxIszNzdG8eXOxT4GImhROaCciFU+uKllaWmLIkCGitHXhwgW88MIL8tfBwcEAgPHjx2Pjxo2iHIOImiguIkpEUJPkSky9e/eGwNutieipmFwRqTOVndBORNTw2HNFREyuiIjEw15tIgKTKyIiEfFuQSJSs+SqpKQEsbGxKCsrU3YoRKTKOCxIpNbUIrkqLCzEpEmToK+vj3bt2iExMREAMHPmTCxbtkzJ0RGRyuA6V0QENUmuQkJCcOnSJRw9ehS6urry8oCAAGzdulWJkRGRauGcKyJSk6UYdu/eja1bt6JLly6QPNZd365dO9y5c0eJkRGRamHPFRGpSc9Veno6rK2tq5QXFBQoJFtERKLg7xUitaYWyVWnTp2wb98++evKhOr777+Hv7+/ssIiIlXDOVdEBDUZFvzkk08QFBSE69evo6ysDF988QWuX7+O06dP49ixY8oOj4hUBudcEZGa9Fx1794dFy9eRFlZGTw9PXHw4EFYW1sjIiICPj4+yg6PiFQFny1IRFCTnisAaNmyJb777jtlh0FEaoHJFZE6U4ueq6ioKFy5ckX++o8//sCQIUPwwQcfoKSkRImREZFqYc8VEalJcjVlyhTcvHkTAHD37l2MHDkS+vr62L59O95//30lR0dEKoPPFiQiqElydfPmTXh5eQEAtm/fjl69emHz5s3YuHEjfv/9d+UGR0QqhHcLEpGaJFeCIEAmkwEADh06hJdffhkA4OjoiIyMDGWGRkSqiMOCRGpNLZKrTp064X//+x9++eUXHDt2DP379wcAxMXFwcbGRsnREZHK4DpXRAQ1Sa5Wr16NqKgozJgxAx9++CFcXV0BADt27EDXrl2VHB0RqQ7OuSIiNVmKoX379gp3C1ZasWIFNDQ0lBAREakkrnNFRFCT5OpJdHV1lR0CEakkJldE6kwtkqvy8nJ8/vnn2LZtGxITE6usbZWVlaWkyIhItXDOFRGpyZyrJUuWYNWqVRg5ciRyc3MRHByMoUOHQiqVYvHixcoOj4hUBudcEZGaJFebNm3Cd999h3fffReampoYPXo0vv/+eyxcuBBnzpxRdnhEpCo454qIoCbJVUpKCjw9PQEAhoaGyM3NBQAMGDAA+/btU2ZoRKRSOCxIRGqSXDVr1gzJyckAKh7gfPDgQQDA+fPnoaOjo8zQiEgVseeKSK2pRXL1yiuvIDw8HAAwc+ZMfPTRR2jVqhXGjRuHiRMnKjk6IlIZfLYgEUFN7hZctmyZ/P8jR45E8+bNERERgVatWmHgwIFKjIyIVAuHBYlITZKr//L394e/v7+ywyAiVcMJ7UQENUqubt26hSNHjiAtLU3+EOdKCxcuVFJURKSamFwRqTO1SK6+++47TJs2DZaWlrC1tYXksb8qJRIJkysiEgl7rohITZKr//3vf/i///s/zJs3T9mhEJEq44R2IoKa3C2YnZ2NV199VdlhEJHK44R2IlKT5OrVV1+Vr21FRFT/mFwRqTOVHRb88ssv5f93dXXFRx99hDNnzsDT0xNaWloKdd95552GDo+IVBLnXBERIBEE1Zwk4OLiUqN6EokEd+/eredo/pWXlwcTExPk5ubC2Ni4wY5LRA3gazOgOAeYEAOYuyk7GiISUW0+v1W25youLk7ZIRCR2uGcKyJSkzlXjxMEASraWUdEjQWHBYnUmtokVz/88AM8PDygq6sLXV1deHh44Pvvv1d2WESkSgT2XBGRmiRXCxcuxKxZszBw4EBs374d27dvx8CBAzFnzpznXkB0zZo1cHZ2hq6uLvz8/HDu3DmRoyaipoe94kSkwnOuHrd27Vp89913GD16tLxs0KBBaN++PWbOnImlS5fWqr2tW7ciODgY69atg5+fH1avXo3AwEDExsbC2tpa7PCJqKngswWJCGrSc1VaWopOnTpVKffx8UFZWVmt21u1ahUmT56MCRMmwN3dHevWrYO+vj5+/PHHKnWLi4uRl5en8EVEqo7JFZE6U4vk6vXXX8fatWurlK9fvx5jx46tVVslJSWIjIxEQECAvEwqlSIgIAARERFV6oeGhsLExET+5ejoWPsTIKImgj1XRKQmw4JAxYT2gwcPokuXLgCAs2fPIjExEePGjUNwcLC83qpVq57aTkZGBsrLy2FjY6NQbmNjg5iYmCr1Q0JCFNrPy8tjgkWkqngnMhFBTZKrq1evomPHjgCAO3fuAAAsLS1haWmJq1evyutJ6uGvTR0dHejo6IjeLhE1RrxbkIjUJLk6cuSIaG1ZWlpCQ0MDqampCuWpqamwtbUV7ThE1IRxWJBIranFnCsxaWtrw8fHB+Hh4fIymUyG8PBw+Pv7KzEyIlI+9lwRkZr0XIktODgY48ePR6dOneDr64vVq1ejoKAAEyZMUHZoRKRMnHNFRGBy9VxGjhyJ9PR0LFy4ECkpKfDy8sKBAweqTHInInXDnisiYnL13GbMmIEZM2YoOwwiaow454pIrXHOFRGRWPhsQSICkysiIhFxzhURMbkiIhIPny1IRGByRUQkIg4LEhGTKyIi8bHnikitMbkiIhIL17kiIjC5IiISEYcFiYjJFRGR+DgsSKTWmFwREYmOyRWROmNyRUQkBoX5VkyuiNQZkysiIlFwMjsRVWByRUQkhsd7rjjnikitMbkiIhIdkysidcbkiohIFOy5IqIKTK6IiMTABUSJ6B9MroiIRMG7BYmoApMrIiKxcViQSK0xuSIiEgPXuSKifzC5IiISBedcEVEFJldERGLgOldE9A8mV0REomNyRaTOmFwREYmCPVdEVIHJFRGRKDjniogqMLkiIhID7xYkon8wuSIiEh2TKyJ1xuSKiEgUnHNFRBWYXBERiYHPFiSifzC5IiISBedcEVEFJldERGLjsCCRWtNUdgAkkqJs4FsHZUdBpMbYc0VEFZhcqZKyR8qOgIisOgBS/molUmf8DaAqtI2BN+Oq3ybGEIUg1P9QR1M/RlNruym0KWZ7DdWWoQOHBYnUHJMrVSHVAEyclR0FERGR2uOEdiIiIiIRMbkiIiIiEhGTKyIiIiIRMbkiIiIiEhEntDcw4Z9HZOTl5Sk5EiIiIqqpys9toQaPumJy1cDy8/MBAI6OjkqOhIiIiGorPz8fJiYmT60jEWqSgpFoZDIZkpKSYGRkBInIa+Hk5eXB0dER9+7dg7GxsahtNwaqfn6A6p8jz6/pU/VzVPXzA1T/HOvr/ARBQH5+Puzt7SGVPn1WFXuuGphUKkWzZs3q9RjGxsYq+QNTSdXPD1D9c+T5NX2qfo6qfn6A6p9jfZzfs3qsKnFCOxEREZGImFwRERERiYjJlQrR0dHBokWLoKOjo+xQ6oWqnx+g+ufI82v6VP0cVf38ANU/x8ZwfpzQTkRERCQi9lwRERERiYjJFREREZGImFwRERERiYjJFREREZGImFypiDVr1sDZ2Rm6urrw8/PDuXPnlB3ScwsNDUXnzp1hZGQEa2trDBkyBLGxsQp1evfuDYlEovA1depUJUVcO4sXL64Se5s2beTbi4qKMH36dFhYWMDQ0BDDhg1DamqqEiOuHWdn5yrnJ5FIMH36dABN89odP34cAwcOhL29PSQSCXbv3q2wXRAELFy4EHZ2dtDT00NAQABu3bqlUCcrKwtjx46FsbExTE1NMWnSJDx8+LABz+LJnnZ+paWlmDdvHjw9PWFgYAB7e3uMGzcOSUlJCm1Ud92XLVvWwGfyZM+6hm+88UaV+Pv166dQp6leQwDV/kxKJBKsWLFCXqcxX8OafC7U5HdnYmIi+vfvD319fVhbW+O9995DWVmZ6PEyuVIBW7duRXBwMBYtWoSoqCh06NABgYGBSEtLU3Zoz+XYsWOYPn06zpw5g7CwMJSWlqJv374oKChQqDd58mQkJyfLv5YvX66kiGuvXbt2CrGfPHlSvm3OnDn4888/sX37dhw7dgxJSUkYOnSoEqOtnfPnzyucW1hYGADg1VdflddpateuoKAAHTp0wJo1a6rdvnz5cnz55ZdYt24dzp49CwMDAwQGBqKoqEheZ+zYsbh27RrCwsKwd+9eHD9+HG+99VZDncJTPe38CgsLERUVhY8++ghRUVHYuXMnYmNjMWjQoCp1ly5dqnBdZ86c2RDh18izriEA9OvXTyH+3377TWF7U72GABTOKzk5GT/++CMkEgmGDRumUK+xXsOafC4863dneXk5+vfvj5KSEpw+fRo//fQTNm7ciIULF4ofsEBNnq+vrzB9+nT56/LycsHe3l4IDQ1VYlTiSUtLEwAIx44dk5f16tVLmDVrlvKCqoNFixYJHTp0qHZbTk6OoKWlJWzfvl1eduPGDQGAEBER0UARimvWrFlCy5YtBZlMJghC0752giAIAIRdu3bJX8tkMsHW1lZYsWKFvCwnJ0fQ0dERfvvtN0EQBOH69esCAOH8+fPyOn/99ZcgkUiEBw8eNFjsNfHf86vOuXPnBABCQkKCvMzJyUn4/PPP6zc4kVR3juPHjxcGDx78xH1U7RoOHjxYePHFFxXKmtI1/O/nQk1+d+7fv1+QSqVCSkqKvM7atWsFY2Njobi4WNT42HPVxJWUlCAyMhIBAQHyMqlUioCAAERERCgxMvHk5uYCAMzNzRXKN23aBEtLS3h4eCAkJASFhYXKCO+53Lp1C/b29mjRogXGjh2LxMREAEBkZCRKS0sVrmebNm3QvHnzJnk9S0pK8Ouvv2LixIkKDypvytfuv+Li4pCSkqJwzUxMTODn5ye/ZhERETA1NUWnTp3kdQICAiCVSnH27NkGj7mucnNzIZFIYGpqqlC+bNkyWFhYwNvbGytWrKiX4Zb6dPToUVhbW8PNzQ3Tpk1DZmamfJsqXcPU1FTs27cPkyZNqrKtqVzD/34u1OR3Z0REBDw9PWFjYyOvExgYiLy8PFy7dk3U+Pjg5iYuIyMD5eXlCt8sAGBjY4OYmBglRSUemUyG2bNno1u3bvDw8JCXjxkzBk5OTrC3t8fly5cxb948xMbGYufOnUqMtmb8/PywceNGuLm5ITk5GUuWLEGPHj1w9epVpKSkQFtbu8qHlo2NDVJSUpQTcB3s3r0bOTk5eOONN+RlTfnaVafyulT3M1i5LSUlBdbW1grbNTU1YW5u3uSua1FREebNm4fRo0crPBT3nXfeQceOHWFubo7Tp08jJCQEycnJWLVqlRKjrbl+/fph6NChcHFxwZ07d/DBBx8gKCgIERER0NDQUKlr+NNPP8HIyKjKdIOmcg2r+1yoye/OlJSUan9OK7eJickVNWrTp0/H1atXFeYkAVCY5+Dp6Qk7Ozv06dMHd+7cQcuWLRs6zFoJCgqS/799+/bw8/ODk5MTtm3bBj09PSVGJr4ffvgBQUFBsLe3l5c15Wun7kpLSzFixAgIgoC1a9cqbAsODpb/v3379tDW1saUKVMQGhraJB6zMmrUKPn/PT090b59e7Rs2RJHjx5Fnz59lBiZ+H788UeMHTsWurq6CuVN5Ro+6XOhMeGwYBNnaWkJDQ2NKndEpKamwtbWVklRiWPGjBnYu3cvjhw5gmbNmj21rp+fHwDg9u3bDRGaqExNTdG6dWvcvn0btra2KCkpQU5OjkKdpng9ExIScOjQIbz55ptPrdeUrx0A+XV52s+gra1tlRtMysrKkJWV1WSua2VilZCQgLCwMIVeq+r4+fmhrKwM8fHxDROgyFq0aAFLS0v596UqXEMAOHHiBGJjY5/5cwk0zmv4pM+FmvzutLW1rfbntHKbmJhcNXHa2trw8fFBeHi4vEwmkyE8PBz+/v5KjOz5CYKAGTNmYNeuXTh8+DBcXFyeuc/FixcBAHZ2dvUcnfgePnyIO3fuwM7ODj4+PtDS0lK4nrGxsUhMTGxy13PDhg2wtrZG//79n1qvKV87AHBxcYGtra3CNcvLy8PZs2fl18zf3x85OTmIjIyU1zl8+DBkMpk8uWzMKhOrW7du4dChQ7CwsHjmPhcvXoRUKq0ylNZU3L9/H5mZmfLvy6Z+DSv98MMP8PHxQYcOHZ5ZtzFdw2d9LtTkd6e/vz+uXLmikCRX/qHg7u4uesDUxG3ZskXQ0dERNm7cKFy/fl146623BFNTU4U7IpqSadOmCSYmJsLRo0eF5ORk+VdhYaEgCIJw+/ZtYenSpcKFCxeEuLg44Y8//hBatGgh9OzZU8mR18y7774rHD16VIiLixNOnTolBAQECJaWlkJaWpogCIIwdepUoXnz5sLhw4eFCxcuCP7+/oK/v7+So66d8vJyoXnz5sK8efMUypvqtcvPzxeio6OF6OhoAYCwatUqITo6Wn633LJlywRTU1Phjz/+EC5fviwMHjxYcHFxER49eiRvo1+/foK3t7dw9uxZ4eTJk0KrVq2E0aNHK+uUFDzt/EpKSoRBgwYJzZo1Ey5evKjwM1l5h9Xp06eFzz//XLh48aJw584d4ddffxWsrKyEcePGKfnM/vW0c8zPzxfmzp0rRERECHFxccKhQ4eEjh07Cq1atRKKiorkbTTVa1gpNzdX0NfXF9auXVtl/8Z+DZ/1uSAIz/7dWVZWJnh4eAh9+/YVLl68KBw4cECwsrISQkJCRI+XyZWK+Oqrr4TmzZsL2tragq+vr3DmzBllh/TcAFT7tWHDBkEQBCExMVHo2bOnYG5uLujo6Aiurq7Ce++9J+Tm5io38BoaOXKkYGdnJ2hrawsODg7CyJEjhdu3b8u3P3r0SHj77bcFMzMzQV9fX3jllVeE5ORkJUZce3///bcAQIiNjVUob6rX7siRI9V+T44fP14QhIrlGD766CPBxsZG0NHREfr06VPl3DMzM4XRo0cLhoaGgrGxsTBhwgQhPz9fCWdT1dPOLy4u7ok/k0eOHBEEQRAiIyMFPz8/wcTERNDV1RXatm0rfPLJJwqJibI97RwLCwuFvn37ClZWVoKWlpbg5OQkTJ48ucofqE31Glb69ttvBT09PSEnJ6fK/o39Gj7rc0EQava7Mz4+XggKChL09PQES0tL4d133xVKS0tFj1fyT9BEREREJALOuSIiIiISEZMrIiIiIhExuSIiIiISEZMrIiIiIhExuSIiIiISEZMrIiIiIhExuSIiIiISEZMrIiIiIhExuSIitXH06FFIJJIqD3dtKOHh4Wjbti3Ky8tFaW/+/PmYOXOmKG0RkXi4QjsRqaTevXvDy8sLq1evlpeVlJQgKysLNjY2kEgkDR6Tj48PgoODMXbsWFHay8jIQIsWLXDx4kW0aNFClDaJqO7Yc0VEakNbWxu2trZKSaxOnjyJO3fuYNiwYaK1aWlpicDAQKxdu1a0Nomo7phcEZHKeeONN3Ds2DF88cUXkEgkkEgkiI+PrzIsuHHjRpiammLv3r1wc3ODvr4+hg8fjsLCQvz0009wdnaGmZkZ3nnnHYWhvOLiYsydOxcODg4wMDCAn58fjh49+tSYtmzZgpdeegm6urryssWLF8PLywu//PILnJ2dYWJiglGjRiE/P19eZ8eOHfD09ISenh4sLCwQEBCAgoIC+faBAwdiy5Yt4rxxRCQKJldEpHK++OIL+Pv7Y/LkyUhOTkZycjIcHR2rrVtYWIgvv/wSW7ZswYEDB3D06FG88sor2L9/P/bv349ffvkF3377LXbs2CHfZ8aMGYiIiMCWLVtw+fJlvPrqq+jXrx9u3br1xJhOnDiBTp06VSm/c+cOdu/ejb1792Lv3r04duwYli1bBgBITk7G6NGjMXHiRNy4cQNHjx7F0KFD8fhsDl9fX9y/fx/x8fHP+W4Rkdg0lR0AEZHYTExMoK2tDX19fdja2j61bmlpKdauXYuWLVsCAIYPH45ffvkFqampMDQ0hLu7O1544QUcOXIEI0eORGJiIjZs2IDExETY29sDAObOnYsDBw5gw4YN+OSTT6o9TkJCgrz+42QyGTZu3AgjIyMAwOuvv47w8HD83//9H5KTk1FWVoahQ4fCyckJAODp6amwf2WbCQkJcHZ2rvmbRET1hskVEak1fX19eWIFADY2NnB2doahoaFCWVpaGgDgypUrKC8vR+vWrRXaKS4uhoWFxROP8+jRI4UhwUrOzs7yxAoA7Ozs5Mfq0KED+vTpA09PTwQGBqJv374YPnw4zMzM5PX19PQAVPTAEVHjwOSKiNSalpaWwmuJRFJtmUwmAwA8fPgQGhoaiIyMhIaGhkK9xxOy/7K0tER2dnaNjl95LA0NDYSFheH06dM4ePAgvvrqK3z44Yc4e/YsXFxcAABZWVkAACsrq5qcLhE1AM65IiKVpK2tLdp6Uo/z9vZGeXk50tLS4OrqqvD1tCFIb29vXL9+vdbHk0gk6NatG5YsWYLo6Ghoa2tj165d8u1Xr16FlpYW2rVr91znQ0TiY88VEakkZ2dnnD17FvHx8TA0NIS5ubko7bZu3Rpjx47FuHHjsHLlSnh7eyM9PR3h4eFo3749+vfvX+1+gYGB+Omnn2p1rLNnzyI8PBx9+/aFtbU1zp49i/T0dLRt21Ze58SJE+jRo4d8eJCIlI89V0SkkubOnQsNDQ24u7vDysoKiYmJorW9YcMGjBs3Du+++y7c3NwwZMgQnD9/Hs2bN3/iPmPHjsW1a9cQGxtb4+MYGxvj+PHjePnll9G6dWssWLAAK1euRFBQkLzOli1bMHny5DqdDxGJiyu0ExE1kPfeew95eXn49ttvRWnvr7/+wrvvvovLly9DU5MDEUSNBXuuiIgayIcffggnJyf5hPW6KigowIYNG5hYETUy7LkiIiIiEhF7roiIiIhExOSKiIiISERMroiIiIhExOSKiIiISERMroiIiIhExOSKiIiISERMroiIiIhExOSKiIiISERMroiIiIhE9P+tgy8e1CjnAQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "class PhaseJump(nx.FunctionTime):\n", " def __init__(self, t0):\n", " super().__init__(\"phase jump\")\n", " self.t0 = t0\n", "\n", " def Function(self, time):\n", " phase = 0\n", " if time > self.t0:\n", " phase = np.pi\n", " return phase\n", "\n", "func = PhaseJump(100)\n", "\n", "detuning = func.CreateDetuning(400, 4000)\n", "\n", "mat = nx.Material.Template(nx.lib.material.Fe)\n", "\n", "site = nx.Hyperfine(magnetic_field = 0,\n", " isotropic = True)\n", "\n", "mat.hyperfine_sites = [site]\n", "\n", "lay = nx.Layer(id = \"iron\",\n", " thickness = 1000, # in nanometer\n", " material = mat)\n", "\n", "sample = nx.Sample(id = \"my sample\",\n", " layers = [lay],\n", " geometry = \"f\",\n", " function_time = func)\n", "\n", "beam = nx.Beam()\n", "beam.LinearSigma()\n", "\n", "exp = nx.Experiment(id = \"my experiment\",\n", " beam = beam,\n", " objects = [sample],\n", " isotope = nx.lib.moessbauer.Fe57)\n", "\n", "amp_time = nx.AmplitudeTime(experiment = exp,\n", " time_length = 200,\n", " time_step = 0.2,\n", " id = \"my amplitude time\",\n", " fft_window = \"none\")\n", "\n", "time_axis, amp = amp_time.Calculate()\n", "\n", "amp_sigma = amp[:,0]\n", "amp_pi = amp[:,1]\n", "\n", "amp_time.Plot(sigma=True, pi=False, polar=True, unwrap=False)" ] } ], "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.11.9" }, "widgets": { "application/vnd.jupyter.widget-state+json": { "state": {}, "version_major": 2, "version_minor": 0 } } }, "nbformat": 4, "nbformat_minor": 5 }