{ "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": "iVBORw0KGgoAAAANSUhEUgAAAlcAAAHgCAYAAACIHEjUAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/SrBM8AAAACXBIWXMAAA9hAAAPYQGoP6dpAAB4qUlEQVR4nO3deVxUVf8H8M8M+77vsom4ICiIiriXBJJraq65Z2lqJllKmVtP4k8fl+oxyTKt1DRNzdRMxV1xxV3BjUVll2UQZJ37+4OYnEAFvTAw83m/Xrx07j33nO+dC8yXc849VyIIggAiIiIiEoVU1QEQERERqRMmV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0REREQiYnJFREREJCImV0QqcPjwYUgkEhw+fFixbcyYMXBzc6uT9hMSEiCRSLBu3bo6aa9CVedNDYObmxvGjBmjeF3X17J79+7o3r17rbejqp8NUi9MrkhjfPPNN5BIJAgICFB1KNVSUFCAefPm1ftEZMyYMZBIJM/9evKDmdTTxo0bsWLFClWHUS0NKVZqeLRVHQBRXdmwYQPc3Nxw5swZ3L59G02aNFF1SEq+++47yOVyxeuCggLMnz8fAOrkL/YX9e677yIoKEjxOj4+HnPmzME777yDLl26KLZ7eHggICAAjx8/hq6uripCJRF17dq10rXcuHEjrl69ig8++EB1gVXT02J1dXXF48ePoaOjo5rASC0wuSKNEB8fj5MnT2Lbtm149913sWHDBsydO1fVYSlpqL/MAwMDERgYqHh97tw5zJkzB4GBgXjrrbcqldfX16/L8KiWSKVStbyWEolELc+L6haHBUkjbNiwARYWFujVqxcGDRqEDRs2VCpTMdfiv//9L1auXInGjRvD0NAQwcHBuHfvHgRBwOeff45GjRrBwMAA/fr1Q1ZWllIdbm5u6N27N/bt2wdfX1/o6+vDy8sL27Zte26MT865SkhIgI2NDQBg/vz5iqG1efPmAXj6/JOq5m3l5ORgzJgxMDMzg7m5OUaPHo2cnJwqY4iNjcWgQYNgaWkJfX19tG3bFjt37nxu7NVV1Tyd7t27w9vbG5cvX0a3bt1gaGiIJk2aYOvWrQCAI0eOICAgAAYGBmjWrBkOHDhQqd4HDx5g3LhxsLOzg56eHlq2bIkffvihUrmkpCTExsZWO85ff/0V8+fPh5OTE0xMTDBo0CDk5uaiqKgIH3zwAWxtbWFsbIyxY8eiqKhIcXy3bt3QunXrKutu1qwZQkJCntn+77//jl69esHR0RF6enrw8PDA559/jrKyMqVyL/vezZs3DxKJBLGxsRg8eDBMTU1hZWWFadOmobCwsFrvUcW17N69O3bv3o3ExETF92vF9+K6desgkUiQkJDwzDoqrF69Gh4eHjAwMED79u1x7NixKmMoKirC3Llz0aRJE+jp6cHZ2Rkff/yx0rWoyrNirWrO1ZgxY2BsbIykpCT07t0bxsbGcHJywsqVKwEAV65cwauvvgojIyO4urpi48aNldrMycnBBx98AGdnZ+jp6aFJkyb4v//7P6XealIf7LkijbBhwwYMGDAAurq6GDZsGFatWoWzZ8+iXbt2VZYtLi7G1KlTkZWVhcWLF2Pw4MF49dVXcfjwYcycORO3b9/G119/jRkzZlT6EL916xaGDBmCiRMnYvTo0Vi7di3efPNN7N27F6+99lq14rWxscGqVaswadIkvPHGGxgwYAAAoFWrVjU6b0EQ0K9fPxw/fhwTJ05EixYtsH37dowePbpS2WvXrqFTp05wcnLCrFmzYGRkhF9//RX9+/fHb7/9hjfeeKNGbddEdnY2evfujaFDh+LNN9/EqlWrMHToUGzYsAEffPABJk6ciOHDh2PJkiUYNGgQ7t27BxMTEwBAWloaOnToAIlEgilTpsDGxgZ//vknxo8fD5lMpjTsM2rUKBw5cgSCIFQrroiICBgYGGDWrFmKa66jowOpVIrs7GzMmzcPp06dwrp16+Du7o45c+YAAEaOHIkJEybg6tWr8Pb2VtR39uxZ3Lx5E7Nnz35mu+vWrYOxsTHCwsJgbGyMgwcPYs6cOZDJZFiyZIlo712FwYMHw83NDRERETh16hS++uorZGdn46effqrW+wQAn376KXJzc3H//n0sX74cAGBsbFzt4yusWbMG7777Ljp27IgPPvgAd+/eRd++fWFpaQlnZ2dFOblcjr59++L48eN455130KJFC1y5cgXLly/HzZs3sWPHDlFjLSsrQ2hoKLp27YrFixdjw4YNmDJlCoyMjPDpp59ixIgRGDBgACIjIzFq1CgEBgbC3d0dQPkQf7du3fDgwQO8++67cHFxwcmTJxEeHo6UlBTO/VJHApGaO3funABA2L9/vyAIgiCXy4VGjRoJ06ZNUyoXHx8vABBsbGyEnJwcxfbw8HABgNC6dWuhpKREsX3YsGGCrq6uUFhYqNjm6uoqABB+++03xbbc3FzBwcFB8PPzU2w7dOiQAEA4dOiQYtvo0aMFV1dXxeuMjAwBgDB37txK59StWzehW7dulbb/u44dO3YIAITFixcrtpWWlgpdunQRAAhr165VbO/Ro4fg4+OjdD5yuVzo2LGj4OnpWamtpzl79myluitUdd7dunUTAAgbN25UbIuNjRUACFKpVDh16pRi+19//VWp7vHjxwsODg5CZmamUltDhw4VzMzMhIKCgkptPU9FnN7e3kJxcbFi+7BhwwSJRCKEhoYqlQ8MDFR633NycgR9fX1h5syZSuXef/99wcjISHj06NEz238y5grvvvuuYGhoqHR9Xva9mzt3rgBA6Nu3r1Jb7733ngBAuHTpkmKbq6urMHr0aMXrqq5lr169lN6HCmvXrhUACPHx8Urb/11HcXGxYGtrK/j6+gpFRUWKcqtXrxYAKH3P//zzz4JUKhWOHTumVGdkZKQAQDhx4kSlOJ70tFgrfg88+T6NHj1aACAsXLhQsS07O1swMDAQJBKJsGnTJsX2ivf/yZ/bzz//XDAyMhJu3ryp1NasWbMELS0tISkp6ZmxUsPDYUFSexs2bICdnR1eeeUVAOVzKoYMGYJNmzZVGmYBgDfffBNmZmaK1xV3F7711lvQ1tZW2l5cXIwHDx4oHe/o6KjUy2NqaopRo0bhwoULSE1NFfXcnmfPnj3Q1tbGpEmTFNu0tLQwdepUpXJZWVk4ePAgBg8ejLy8PGRmZiIzMxMPHz5ESEgIbt26Vek8xWRsbIyhQ4cqXjdr1gzm5uZo0aKF0t2dFf+/e/cugPKeud9++w19+vSBIAiKuDMzMxESEoLc3FzExMQojj98+HC1e62A8p6uJ+fCBQQEQBAEjBs3TqlcQEAA7t27h9LSUgCAmZkZ+vXrh19++UXRXllZGTZv3oz+/fvDyMjome0aGBgo/l9xPbp06YKCgoJKw5ov+t49afLkyUqvK74/9uzZ88w4xXbu3Dmkp6dj4sSJShPlK4a1n7Rlyxa0aNECzZs3V7rur776KgDg0KFDosf39ttvK/5vbm6OZs2awcjICIMHD1Zsr3j/n3yft2zZgi5dusDCwkIp1qCgIJSVleHo0aOix0qqxWFBUmtlZWXYtGkTXnnlFcTHxyu2BwQEYOnSpYiKikJwcLDSMS4uLkqvK36pPzkk8eT27Oxspe1NmjSBRCJR2ta0aVMA5fM57O3tX+KMaiYxMREODg6VhjyaNWum9Pr27dsQBAGfffYZPvvssyrrSk9Ph5OTU63E2ahRo0rvmZmZ2XPf84yMDOTk5GD16tVYvXp1lXWnp6e/cFw1+V6Qy+XIzc2FlZUVgPLEbPPmzTh27Bi6du2KAwcOIC0tDSNHjnxuu9euXcPs2bNx8OBByGQypX25ublKr1/0vXuSp6en0msPDw9IpdJKc6RqW2JiYpXx6OjooHHjxkrbbt26hRs3bijmJv7by1z3qujr61dqy8zM7Knv/5Pv861bt3D58uU6i5VUj8kVqbWDBw8iJSUFmzZtwqZNmyrt37BhQ6XkSktLq8q6nra9Jj0hYpFIJFW2W1VPXHVUTKqdMWPGUydb1+bSFS/6nlfE/dZbb1U5jwyo+Tw1MeICgJCQENjZ2WH9+vXo2rUr1q9fD3t7e6VlK6qSk5ODbt26wdTUFAsWLICHhwf09fURExODmTNnVpoAXRvfr/9OFl7W0+p70e9XoPza+/j4YNmyZVXu/3dy+bJe5n2Wy+V47bXX8PHHH1dZtuKPL1IfTK5IrW3YsAG2traKu3qetG3bNmzfvh2RkZFKwzAvq6IX6MkPlJs3bwJAjVZgf9YHnIWFRZXDOxV/+VdwdXVFVFQUHj16pNR7FRcXp1SuoldAR0fnuR/+9YmNjQ1MTExQVlZW7+LW0tLC8OHDsW7dOvzf//0fduzYgQkTJjz1w7jC4cOH8fDhQ2zbtg1du3ZVbH+y51Vst27dUky+Bsq/h+VyeY2fGPC071kLCwsAqHSXalXfrxXxVAzvAUBJSQni4+OV7sD08PDApUuX0KNHjxdKBsVOIJ/Fw8MDjx49qnffo1R7OOeK1Nbjx4+xbds29O7dG4MGDar0NWXKFOTl5Ym61AAAJCcnY/v27YrXMpkMP/30E3x9fWs0JGhoaAig8gcSUP7LOjY2FhkZGYptly5dwokTJ5TKvf766ygtLcWqVasU28rKyvD1118rlbO1tUX37t3x7bffIiUlpVJ7T7ZTn2hpaWHgwIH47bffcPXq1Ur7/x13dZdiEMvIkSORnZ2Nd999F48ePapy3a9/q0i+nuz5KC4uxjfffFNrcf77j4+K74/Q0NAa1WNkZFRp2BIo/34FoDS3qKysrNJQbtu2bWFjY4PIyEgUFxcrtq9bt67Sz8HgwYPx4MEDfPfdd5Xae/z4MfLz818o1towePBgREdH46+//qq0LycnRzFXj9QHe65Ibe3cuRN5eXno27dvlfs7dOgAGxsbbNiwAUOGDBGt3aZNm2L8+PE4e/Ys7Ozs8MMPPyAtLQ1r166tUT0GBgbw8vLC5s2b0bRpU1haWsLb2xve3t4YN24cli1bhpCQEIwfPx7p6emIjIxEy5Ytlebo9OnTB506dcKsWbOQkJCgWHOrqg+VlStXonPnzvDx8cGECRPQuHFjpKWlITo6Gvfv38elS5de+r2pDYsWLcKhQ4cQEBCACRMmwMvLC1lZWYiJicGBAweU1iKr6VIML8vPzw/e3t6Kyddt2rR57jEdO3aEhYUFRo8ejffffx8SiQQ///xzrcYcHx+Pvn37omfPnoiOjsb69esxfPjwp67V9TT+/v7YvHkzwsLC0K5dOxgbG6NPnz5o2bIlOnTogPDwcGRlZcHS0hKbNm2qlFTo6OjgP//5D9599128+uqrGDJkCOLj47F27dpKc65GjhyJX3/9FRMnTsShQ4fQqVMnlJWVITY2Fr/++iv++usvtG3btsax1oaPPvoIO3fuRO/evTFmzBj4+/sjPz8fV65cwdatW5GQkABra+taaZtUgz1XpLY2bNgAfX39p64tJZVK0atXL+zduxcPHz4UrV1PT09s3rwZe/bswaxZs1BSUoLNmzc/d+HIqnz//fdwcnLC9OnTMWzYMMXikC1atMBPP/2E3NxchIWFYefOnfj5558rfXhLpVLs3LkTI0aMwPr16/Hpp5/CyckJP/74Y6W2vLy8cO7cOfTq1Qvr1q3D5MmTERkZCalUqli/qT6ys7PDmTNnMHbsWGzbtg1TpkzBl19+iaysLPzf//2fqsPDqFGjAKBaE9kBwMrKCrt27YKDgwNmz56N//73v3jttdewePHiWotx8+bN0NPTw6xZs7B7925MmTIFa9asqXE97733HoYPH461a9di+PDhSnelbtiwAR07dsSiRYuwcOFCvPLKK1i0aFGlOt555x188803SE5OxkcffYRjx45h586dleZQSaVS7NixA4sWLcKVK1cwY8YMzJ8/H2fPnsW0adOeO4/pWbGKzdDQEEeOHMFHH32Ew4cPY9q0aVi0aBFu3bqF+fPnV7oTkho+iaCK2bhEasrNzQ3e3t7YtWuXqkOheuLLL7/E9OnTkZCQUOnuQ1WbN28e5s+fj4yMDPacEImIPVdERLVEEASsWbMG3bp1q3eJFRHVHs65IiISWX5+Pnbu3IlDhw7hypUr+P3331UdEhHVISZXREQiy8jIwPDhw2Fubo5PPvnkqTdVEJF64pwrIiIiIhFxzhURERGRiJhcEREREYmIyRURERGRiJhcEREREYmIyRURERGRiJhcEREREYmIyRURERGRiJhcEZHGW7duHSQSCRISEhpEvfW1XSIqx+SKiKgBOnnyJObNm4ecnBxVh0JE/8LkioiolowcORKPHz+Gq6ur6HWfPHkS8+fPrzK5qs12iej5+GxBIiKR5efnw8jICFpaWtDS0qrz9lXVLhGVY88VEanM6tWr0aZNGxgaGkIikSh9NW7cuMpj5s2bB4lEgtjYWAwePBimpqawsrLCtGnTUFhYWKn8hQsXEBoaClNTUxgbG6NHjx44derUc2NLTEzEe++9h2bNmsHAwABWVlZ48803K81jqojn+vXrGD58OCwsLNC5c2cAynOfEhISKp3jk181aXfevHn46KOPAADu7u6KOirKVDXnqrrvQ8X53L59G2PGjIG5uTnMzMwwduxYFBQUPPd9IyL2XBGRikyfPh0rVqxAcHAwxo4di/v372P58uUoKSlB79694e/v/8zjBw8eDDc3N0RERODUqVP46quvkJ2djZ9++klR5tq1a+jSpQtMTU3x8ccfQ0dHB99++y26d++OI0eOICAg4Kn1nz17FidPnsTQoUPRqFEjJCQkYNWqVejevTuuX78OQ0NDpfJvvvkmPD09sXDhQgiCUKk+Gxsb/Pzzz0rbSkpKMH36dOjq6tao3QEDBuDmzZv45ZdfsHz5clhbWyvaqMqLvA+DBw+Gu7s7IiIiEBMTg++//x62trb4v//7v6e+Z0T0N4GIqI4dPXpUACBMmjRJafv8+fMFAMKZM2eeeuzcuXMFAELfvn2Vtr/33nsCAOHSpUuKbf379xd0dXWFO3fuKLYlJycLJiYmQteuXRXb1q5dKwAQ4uPjFdsKCgoqtR0dHS0AEH766adK8QwbNqxS+arq/XfMWlpawsGDB2vc7pIlS55a97/bre778OT5jBs3Tmn7G2+8IVhZWVV5HkSkjMOCRFTnli9fDktLSyxZskRpe7du3QAAN2/efG4dkydPVno9depUAMCePXsAAGVlZdi3bx/69++vNMTo4OCA4cOH4/jx45DJZE+t38DAQPH/kpISPHz4EE2aNIG5uTliYmIqlZ84ceJzY37STz/9hG+++QaLFy/GK6+88sLtPs+Lvg//Pp8uXbrg4cOHz3zPiKgckysiqlOlpaXYv38/QkNDYWRkpLSvuLgYAGBqavrcejw9PZVee3h4QCqVKuYZZWRkoKCgAM2aNat0bIsWLSCXy3Hv3r2n1v/48WPMmTMHzs7O0NPTg7W1NWxsbJCTk4Pc3NxK5d3d3Z8bc4WLFy9i4sSJGDZsGMLCwl6q3ed50ffBxcVF6bWFhQUAIDs7u8YxEGkazrkiojqVkJCAR48ewdvbu9K+8+fPAyj/0K+pJyeFi2Hq1KlYu3YtPvjgAwQGBsLMzAwSiQRDhw6FXC6vVP7JHqdnyc7OxsCBA9G0aVN8//33L91ubXna3YZCFfPJiEgZkysiqlN5eXkAoDSJGyj/0N6yZQtatmyJJk2aPLeeW7duKfUW3b59G3K5HG5ubgDKJ3cbGhoiLi6u0rGxsbGQSqVwdnZ+av1bt27F6NGjsXTpUsW2wsLCl1q0Uy6XY8SIEcjJycGBAwcqTYqvSbvVTSZf9n0goprjsCAR1amK4aYDBw4obV+xYgViYmIQHh5erXpWrlyp9Prrr78GAISGhgIo73kJDg7G77//rrQkQVpaGjZu3IjOnTs/c/hRS0urUi/N119/jbKysmrFV5X58+fjr7/+wi+//PLUYcTqtlsxpPq8ZO9l3wciqjn2XBFRnbKyskL//v2xY8cOjBgxAp06dcLx48fxyy+/4O2338aIESOqVU98fDz69u2Lnj17Ijo6GuvXr8fw4cPRunVrRZn//Oc/2L9/Pzp37oz33nsP2tra+Pbbb1FUVITFixc/s/7evXvj559/hpmZGby8vBAdHY0DBw7Aysrqhc77ypUr+Pzzz9G1a1ekp6dj/fr1SvvfeuutGrVbsVTFp59+iqFDh0JHRwd9+vSpNI/tZd8HIqo5JldEVOfWrl0Lc3Nz/P777/jtt9/g5eWFNWvWYNy4cdWuY/PmzZgzZw5mzZoFbW1tTJkypdLdhy1btsSxY8cQHh6OiIgIyOVyBAQEYP369c9c4woAvvzyS2hpaWHDhg0oLCxEp06dcODAAYSEhLzQOT98+BCCIODIkSM4cuRIpf0VyVV1223Xrh0+//xzREZGYu/evZDL5YiPj68yuXqZ94GIak4icHYiETUg8+bNw/z585GRkaFYPJOIqD7hnCsiIiIiETG5IiIiIhIRkysiIiIiEXHOFREREZGI2HNFREREJCIuxVDH5HI5kpOTYWJiIvrjOoiIiKh2CIKAvLw8ODo6Qip9dt8Uk6s6lpyczEdNEBERNVD37t1Do0aNnlmGyVUdMzExAVB+cfjICSIiooZBJpPB2dlZ8Tn+LEyu6ljFUKCpqSmTKyIiogamOlN6OKGdiIiISERMroiIiIhExOSKiIiISERMrtRIaZkcN1Jk4LqwREREqsPkSo18tPUyQr88hm+P3lV1KERERBqLyZUa2X7hAQBg5cHbKo6EiIhIczG5UkNlHBYkIiJSGSZXaqhMzuSKiIhIVeptcrVy5Uq4ublBX18fAQEBOHPmzDPLb9myBc2bN4e+vj58fHywZ88exb6SkhLMnDkTPj4+MDIygqOjI0aNGoXk5GSlOrKysjBixAiYmprC3Nwc48ePx6NHj5TKXL58GV26dIG+vj6cnZ2xePFi8U5aJOy4IiIiUp16mVxt3rwZYWFhmDt3LmJiYtC6dWuEhIQgPT29yvInT57EsGHDMH78eFy4cAH9+/dH//79cfXqVQBAQUEBYmJi8NlnnyEmJgbbtm1DXFwc+vbtq1TPiBEjcO3aNezfvx+7du3C0aNH8c477yj2y2QyBAcHw9XVFefPn8eSJUswb948rF69uvbejBcgZ3ZFRESkMhKhHt63HxAQgHbt2uF///sfAEAul8PZ2RlTp07FrFmzKpUfMmQI8vPzsWvXLsW2Dh06wNfXF5GRkVW2cfbsWbRv3x6JiYlwcXHBjRs34OXlhbNnz6Jt27YAgL179+L111/H/fv34ejoiFWrVuHTTz9FamoqdHV1AQCzZs3Cjh07EBsbW61zk8lkMDMzQ25uruiPv3GbtRsAIJEA8RG9RK2biIhIk9Xk87ve9VwVFxfj/PnzCAoKUmyTSqUICgpCdHR0lcdER0crlQeAkJCQp5YHgNzcXEgkEpibmyvqMDc3VyRWABAUFASpVIrTp08rynTt2lWRWFW0ExcXh+zs7CrbKSoqgkwmU/qqbfUvXSYiItIc9S65yszMRFlZGezs7JS229nZITU1tcpjUlNTa1S+sLAQM2fOxLBhwxTZZ2pqKmxtbZXKaWtrw9LSUlHP09qp2FeViIgImJmZKb6cnZ2rLEdERETqod4lV7WtpKQEgwcPhiAIWLVqVa23Fx4ejtzcXMXXvXv3ar1NIiIiUh1tVQfwb9bW1tDS0kJaWprS9rS0NNjb21d5jL29fbXKVyRWiYmJOHjwoNKYqb29faUJ86WlpcjKylLU87R2KvZVRU9PD3p6ek87XSIiIlIz9a7nSldXF/7+/oiKilJsk8vliIqKQmBgYJXHBAYGKpUHgP379yuVr0isbt26hQMHDsDKyqpSHTk5OTh//rxi28GDByGXyxEQEKAoc/ToUZSUlCi106xZM1hYWLz4SRMREZHaqHfJFQCEhYXhu+++w48//ogbN25g0qRJyM/Px9ixYwEAo0aNQnh4uKL8tGnTsHfvXixduhSxsbGYN28ezp07hylTpgAoT6wGDRqEc+fOYcOGDSgrK0NqaipSU1NRXFwMAGjRogV69uyJCRMm4MyZMzhx4gSmTJmCoUOHwtHREQAwfPhw6OrqYvz48bh27Ro2b96ML7/8EmFhYXX8DhEREVF9Ve+GBYHypRUyMjIwZ84cpKamwtfXF3v37lVMHk9KSoJU+k9e2LFjR2zcuBGzZ8/GJ598Ak9PT+zYsQPe3t4AgAcPHmDnzp0AAF9fX6W2Dh06hO7duwMANmzYgClTpqBHjx6QSqUYOHAgvvrqK0VZMzMz7Nu3D5MnT4a/vz+sra0xZ84cpbWwiIiISLPVy3Wu1FldrHMFAD+Oa49uTW1ErZ+IiEhTNeh1rkgcqw7fVnUIREREGonJlZris5uJiIhUg8mVmjoTn4WSMrmqwyAiItI4TK7U2KazXLCUiIiorjG5UmOf7biK4lL2XhEREdUlJldqLmTFUVWHQEREpFGYXKm5+Mx8DFx1kj1YREREdYTJlQY4n5iNprP/xLXkXFWHQkREpPaYXGmQXl8dx4BvTiCvsOT5hYmIiOiFMLlSI1LJ88vEJOXAZ94+zN5xBWVcDIuIiEh0TK7UiFRSjezqb+tPJcHjkz1Y8Md1FBSX1mJUREREmoXJlRqpSXJV4YcT8fCa8xcmb4xBuqywFqIiIiLSLNqqDoBE9JTcylBXCwXFZc88dPflFOy+nAJnSwP8MLodmtgaQ/ICyRoREZGmY3KlRgRBeQ7VgbCuaGJrongtlwu4n/0YW2Pu46uoW1XWcS/rMV5bXr421vIhrdHf14lJFhERUQ1IhH9/IlOtkslkMDMzQ25uLkxNTUWtu3H4bsUDm1s7m+P3yZ2eWT6/qBRrjsdj2f6bzyz3TtfGmPpqE5jo64gVKhERUYNSk89vJld1rDaTK7dZuxX/v/VFKHS0qj+l7uGjInwZdQs/RSc+tYyHjRG+HdkWTWyNXypOIiKihqYmn9+c0K6mapJYAYCVsR4W9PNGwqJe2De9K/xdLSqVuZORj6BlR+A2azfWnYjnqu9ERERVYHJFlTS1M8FvkzoiYVEvLBvcusoy8/64jqaz/8SoH87gfnZBHUdIRERUf3FCOz3TgDaNMKBNI+QWlOA/u69jy/n7SvuP3sxA5/87BAD4dqQ/gr3sOAGeiIg0Gnuu1IiVkW6t1W1mqIMlb7ZGfMTr+H1yJ9iZ6lUq8+7P5+EevgfTNl1Aeh7XzCIiIs3E5EqNeDmWT7AbEeBSa21IJBK0djbH6U+CcGNBT4SHNq9U5veLyWj/RRSaf/Ynjt/KrLVYiIiI6qN6m1ytXLkSbm5u0NfXR0BAAM6cOfPM8lu2bEHz5s2hr68PHx8f7NmzR2n/tm3bEBwcDCsrK0gkEly8eFFpf0JCAiQSSZVfW7ZsUZSrav+mTZtEO++XUXHfZzs3yzppz0BXC+9280DCol7YP70r2v+r3cISOd5acxpus3Zj9o4ryMovrpO4iIiIVKleJlebN29GWFgY5s6di5iYGLRu3RohISFIT0+vsvzJkycxbNgwjB8/HhcuXED//v3Rv39/XL16VVEmPz8fnTt3xv/93/9VWYezszNSUlKUvubPnw9jY2OEhoYqlV27dq1Suf79+4t27i9DQHl2pYopT552Jvh1YiBu/icUXw71rbR//akktPl8P9xm7cbxW5koLeOdhkREpJ7q5TpXAQEBaNeuHf73v/8BAORyOZydnTF16lTMmjWrUvkhQ4YgPz8fu3btUmzr0KEDfH19ERkZqVQ2ISEB7u7uuHDhAnx9fZ8Zh5+fH9q0aYM1a9YotkkkEmzfvv2FE6raXOdq+HencPLOQ3w51Bf9fJ1ErftFZOcXI/LoHXx75G6V+1/3sccX/X1gUYtzxYiIiMSgknWuSkpKcO/ePcTFxSErK+uF6ykuLsb58+cRFBSk2CaVShEUFITo6Ogqj4mOjlYqDwAhISFPLV8d58+fx8WLFzF+/PhK+yZPngxra2u0b98eP/zwQ6XHzjypqKgIMplM6UtTWBjpIjy0BRIW9cK29zrC3dpIaf+eK6nw+7s369dz9575PhIRETUUL5Vc5eXlYdWqVejWrRtMTU3h5uaGFi1awMbGBq6urpgwYQLOnj1bozozMzNRVlYGOzs7pe12dnZITU2t8pjU1NQala+ONWvWoEWLFujYsaPS9gULFuDXX3/F/v37MXDgQLz33nv4+uuvn1pPREQEzMzMFF/Ozs4vHNPzVOQm9XEphDYuFjg0oztufxGK2b1aVNr/8dbLcA/fg4GrTuJ2+iMVREhERCSOF17natmyZfjiiy/g4eGBPn364JNPPoGjoyMMDAyQlZWFq1ev4tixYwgODkZAQAC+/vpreHp6ihl7rXn8+DE2btyIzz77rNK+J7f5+fkhPz8fS5Yswfvvv19lXeHh4QgLC1O8lslktZZgKeZc1Urt4tDWkuLtLo3xdpfGSM0txEdbL+HYE3cUnk/MRtCyIwCA8NDmGBnoCkNdLsdGREQNR40+tbS0tFBWVgYAOHv2LI4ePYqWLVtWWbZ9+/YYN24cIiMjsXbtWhw7dqxayZW1tTW0tLSQlpamtD0tLQ329vZVHmNvb1+j8s+zdetWFBQUYNSoUc8tGxAQgM8//xxFRUXQ06u89pOenl6V22vDPz1XddLcS7M308fP4wMglws4cjMD4348iydHBiP+jEXEn7FwMNPHl0P90M7Nol72yhERET2pRsOCT86J+eWXX56aWD1JT08PEydOxLhx46rVhq6uLvz9/REVFaXYJpfLERUVhcDAwCqPCQwMVCoPAPv3739q+edZs2YN+vbtCxsbm+eWvXjxIiwsLOosgXqWiqsjqdd9V5VJpRK80twW8RG9cGNBT7zTtbHS/pTcQgz+Nhru4XswY8slZOQVqShSIiKi56tRz1Vd9RqEhYVh9OjRaNu2Ldq3b48VK1YgPz8fY8eOBQCMGjUKTk5OiIiIAABMmzYN3bp1w9KlS9GrVy9s2rQJ586dw+rVqxV1ZmVlISkpCcnJyQCAuLg4AOW9Xk/2cN2+fRtHjx6ttE4WAPzxxx9IS0tDhw4doK+vj/3792PhwoWYMWNGrb0XmsZAVwufvN4Cn7zeAncyHmHqxgu4nvLPTQBbz9/H1vP3IZEAK4b4ItTbAbra9XJFESIi0lD1cjLLkCFDkJGRgTlz5iA1NRW+vr7Yu3evYtJ6UlISpNJ/PlA7duyIjRs3Yvbs2fjkk0/g6emJHTt2wNvbW1Fm586diuQMAIYOHQoAmDt3LubNm6fY/sMPP6BRo0YIDg6uFJeOjg5WrlyJ6dOnQxAENGnSBMuWLcOECRPEfgteTAMbFnweDxtj7JnWBcWlcuy7noopGy8o9gkCMG3TRUzDRbR2Nsd/B7VCE1tjDhsSEZHK1WidqyfnXD3p8ePHEAQBhoaGAIDExERs374dXl5eVSYpmqw217l6M/IkziZkY9WINgj1cRC17voiK78Yy/ffxM+nEqvcP6ajGz4MbgoTfZ06joyIiNRZna9z1a9fP/z0008AgJycHAQEBGDp0qXo168fVq1aJUYTVA0NbUL7i7A00sXn/b2RsKgX/pjSGX4u5kr7151MgM+8fQhYeAC/X3yAMjnXziIiorolSnIVExODLl26ACi/087Ozg6JiYn46aef8NVXX4nRBFXDP2mEGmdXT/BpZIbt73VSPHLH8omV3tNkRZi26SI8PtmDYatP4XZ6HhcpJSKiOiHKnKuCggKYmJgAAPbt24cBAwZAKpWiQ4cOSEyseviGSCy62lL083VCP18nZOcXY9WRO1h99J9H7kTffYigZUehJZVgRIALPgxuBjMDDhsSEVHtEKXnqkmTJtixYwfu3buHv/76SzHPKj09XfR5RfR0FT0z6jws+DwWRrr45PXyR+7seb8LejS3Vewrkwv4KToRrefvQ5fFB7H+VCKHDYmISHSiJFdz5szBjBkz4ObmhoCAAMX6Uvv27YOfn58YTVA1/LPOFQGAl6Mp1oxphzsLX8e3I/3h5fBPon8v6zFm77iqGDY8l5DFYUMiIhKFKHcLAuXP90tJSUHr1q0VyyScOXMGpqamaN68uTjRqoHavFuw/8oTuHgvB9+NaovXvOyef4AGelRUio2nE/HdsfhKi5FKJcDANo0wLcgTjSwMVRQhERHVRzX5/BZtnat/L8YJlD8Ch+oOe66ez1hPG+909cA7XT1wP7sA3xy+gy3n7qGkTIBcALacv48t5+/DwUwfQ9u5YGxnN5hyWQciIqoB0ZKrqKgoREVFIT09HXK5XGnfDz/8IFYz9Cwc1qqRRhaGWPiGD77o740L93IQefgO9l0vf0ZlSm4hlh+4ieUHbqJVIzOM7eSG130coKetpeKoiYiovhMluZo/fz4WLFiAtm3bwsHBgatkqxjf/pqRSCRo42KB1aPaoqRMjoOx6fju6F2cS8wGAFy+n4vpmy9h+uZLeN3HHpO6NYFPIzMVR01ERPWVKMlVZGQk1q1bh5EjR4pRHb0gxbAgk6sXpqMlRUhLe4S0tMfj4jJsPJOEjacTcScjHwCw50oq9lxJhbWxHvr7OmJEB1e4WxupOGoiIqpPREmuiouL0bFjRzGqopegWKGds65EYaCrhfGd3TG+szsyHxXh+2Px2H0lGfeyHpe/Ph6P74/Ho7GNEd7wdcKoQDeYGXJ+FhGRphNlKYa3334bGzduFKMqegmC4snNqo1DHVkb62FWaHMc+/hVnP6kByZ0cYe1sR4A4G5GPpbuv4nWC/bhzciT+OVMEgpLqr6rloiI1J8oPVeFhYVYvXo1Dhw4gFatWkFHR/mv92XLlonRDD0H57PXDTtTfXzaywuf9vLC1Qe5WHsiAX9cTkZxqRxnE7JxNiEb4duuIKiFHd7wc0JISztoa4nydwwRETUAoiRXly9fhq+vLwDg6tWrSvs4ub3u8R2vO95OZlg6uDWWDGqFg7Hp2HbhPvZcSQUAHLiRhgM30qAtlaBXKwcMaeeMjh7WKo6YiIhqmyjJ1aFDh8Sohl6SYs4VE9o6J5VKEORlhyAvO5SUybH7cgq2XXiAozczUCoX8PvFZPx+MRkmetro7+eE4QEuaOHAR0MREamjl0qu5syZg379+sHf31+seOglcBHR+kFHS4r+fk7o7+eEx8Vl+PXcPey4+AAXknKQV1SKn08l4udTiXAyN0BfX0cMbOOEJrYmqg6biIhE8lLJ1f379xEaGgpdXV306dMHffv2RY8ePaCrqytWfFQDfHBz/WOgq4XRHd0wuqMbMvKKsPX8ffxxKRnXU2R4kPMYqw7fwarDd+BhY4Q+rR0xsE0jOFvy0TtERA3ZSz9bUC6X48SJE/jjjz/w+++/IyUlBa+99hr69euH3r17w9LSUvSgG7LafLZgzxVHEZuah5/Ht0cXTxtR6yZxJec8xq/n7mHnpWTc/XsNrQotHU3Rt7UjBvk3gtXfdyQSEZFq1eTzW7QHN1e4ceOGItE6f/482rdvj759+2LYsGFwcnKqblNqqy6Sq/XjA9DZkxOnG4qEzHxsPJOEPy4lIyW3UGlfOzcL9G7liDfaOPEZh0REKqTS5OpJGRkZ2LlzJ3bu3IkuXbpgxowZ1T5WXdVmchWy/Cji0vKw4e0AdGrC5KohunQvB9svPMCuyynIfFSktC+wsRV6t3ZAn9aOTLSIiOpYnSVX77zzDpYtWwZjY+MXj1bD1GZyFbz8CG6mPcLGtwPQkclVgyYIAs4mZGPnpQfYfTkF2QUlSvu7eFqjTytH9G7tAENd0Z6/TkRET1FnyZWWlhZSUlJga2v74tFqmNpMrl5bdgS30h9h44QArqekRgRBwKm7Wdh+4T52XU5BQbFy73FQC1u87uOAUG8HGOhqqShKIiL1VpPP75daNroGeVmNrVy5Em5ubtDX10dAQADOnDnzzPJbtmxB8+bNoa+vDx8fH+zZs0dp/7Zt2xAcHAwrKytIJBJcvHixUh3du3eHRCJR+po4caJSmaSkJPTq1QuGhoawtbXFRx99hNLS0pc+XzFwgXb1JJFIEOhhhcWDWuPa/BD8PL49BrRxgtHfidSBG+kI+/USWs7di7d/PIcdFx7w8TtERCpUL8cTNm/ejLCwMERGRiIgIAArVqxASEgI4uLiquwlO3nyJIYNG4aIiAj07t0bGzduRP/+/RETEwNvb28AQH5+Pjp37ozBgwdjwoQJT217woQJWLBggeK1oeE/t8WXlZWhV69esLe3x8mTJ5GSkoJRo0ZBR0cHCxcuFPEdeDGKpRi40pXakkgk6OJpgy6eNigdKMfRWxnYdTkFf11NRX5xmWJVeOmvQI8WdujdygEhLe2hr8MeLSKiuvLSw4KpqamwsRH3tv+AgAC0a9cO//vf/wCUL/fg7OyMqVOnYtasWZXKDxkyBPn5+di1a5diW4cOHeDr64vIyEilsgkJCXB3d8eFCxcUj+yp0L17d/j6+mLFihVVxvXnn3+id+/eSE5Ohp2dHQAgMjISM2fOREZGRpXrexUVFaGo6J+JyTKZDM7OzrUyLNhj6WHcycjHpnc6oENjK1HrpvqttEyOIzcz8MelZPx5NRVFpXLFPm2pBK82t0WvVg4I9rLn0CER0QuoybDgS/VcCYKAMWPGQE/v2WvxbNu2rdp1FhcX4/z58wgPD1dsk0qlCAoKQnR0dJXHREdHIywsTGlbSEgIduzYUe12K2zYsAHr16+Hvb09+vTpg88++0zRexUdHQ0fHx9FYlXRzqRJk3Dt2jX4+flVqi8iIgLz58+vcRwvgiu0ay5tLSl6tLBDjxZ2WCoXcDguHX9cSsa+62koKC7Dvutp2He9/DmH3ZvZolcre4S0tOdkeCKiWvBSv1lHjx4tVhwKmZmZKCsrU0pgAMDOzg6xsbFVHpOamlpl+dTU1Bq1PXz4cLi6usLR0RGXL1/GzJkzERcXp0gOn9ZOxb6qhIeHKyV+FT1XtYLPFiQAWlKJItGq6NHacyUV+66lIq+o9J+hQ8klvNLMFqE+DghpaQcTLu9ARCSKl0qu1q5dK1Yc9cI777yj+L+Pjw8cHBzQo0cP3LlzBx4eHi9Up56e3nN79sTCCe30b0/2aJWW+eDYrUzsvpKCvVdT8aioFFGx6YiKTcfHW4GuTW0Q6m2PUB8HrqNFRPQSapRc1ebdgRWsra2hpaWFtLQ0pe1paWmwt7ev8hh7e/sala+ugIAAAMDt27fh4eEBe3v7SnctVrT7sm2Jgc8WpGfR1pLilea2eKW5LRYPbIUTdzKx+3IK9l5LRU5BCQ7HZeBwXAZm/nYFnZpYoae3A3q2tIeNCR/BQ0RUEzVaiuHSpUv44osvsGzZMly6dKlWAtLV1YW/vz+ioqIU2+RyOaKiohAYGFjlMYGBgUrlAWD//v1PLV9dFcs1ODg4KNq5cuUK0tPTldoxNTWFl5fXS7UlJuZW9DxSafldh4sGtkLM7NewcUIARnZwVSRSJ24/xGc7rqLdFwcwaNVJrDkej9R/PZqHiIiqVqOeKx8fH/j4+KC4uBhRUVHYvXs3tLW10alTJ3Tq1Em0oMLCwjB69Gi0bdsW7du3x4oVK5Cfn4+xY8cCAEaNGgUnJydEREQAAKZNm4Zu3bph6dKl6NWrFzZt2oRz585h9erVijqzsrKQlJSE5ORkAEBcXByA8h4ne3t73LlzBxs3bsTrr78OKysrXL58GdOnT0fXrl3RqlUrAEBwcDC8vLwwcuRILF68GKmpqZg9ezYmT55cZ0N/z6KY0M7simpAKpWgo4c1OnpYY0G/lohJysbuy6nYezUFybmFOJeYjXOJ2fh813W0amSGnt72eN3bAW7WRqoOnYioXqrRUgxVEQQBx48fx4kTJyAIAtq0aYMePXpAW/vl7kL63//+hyVLliA1NRW+vr746quvFMN03bt3h5ubG9atW6cov2XLFsyePRsJCQnw9PTE4sWL8frrryv2r1u3TpGcPWnu3LmYN28e7t27h7feegtXr15Ffn4+nJ2d8cYbb2D27NlKt1wmJiZi0qRJOHz4MIyMjDB69GgsWrSo2udbmyu0d118CElZBfhtUkf4u1qIWjdppqsPcsvX0bqWivjMfKV9zexMEOJtj54t7eHlKO73MhFRfVNrj7+pjgsXLiAqKgrFxcVo3rw5QkNDYWBgIGYTDVptJlddFh/EvazHTK6oVtxMy8OfV1Kx91oqbqTIlPY5WxogxKt8MryfszmkUnafEpF6qbN1rqri5+enWO/pwoULCA0NhYWFBSIiItC8eXOxm6MnCIqlGFQbB6mnpnYmaGpngmlBnkjOeYxdl5Ox92oqYpJycC/rMb4/Ho/vj8fD2lgPwS3tENLSHh09rKCj9VJP2SIianBET65kMhl27tyJrVu34u7duwgNDcXAgQOZWNUh5lZU2xzNDfBOVw+809UD6XmF+OtaGvZdS8WJ25nIfFSEjaeTsPF0Eoz1tPFqc1sEt7TDK81sYaTHRUuJSP2J9puuoKAAgwYNQkpKCvr164cFCxYoJoJT3RC4iCipgK2JPkZ2cMXIDq54VFSKA9fT8Ne1VByOy8CjolLsvJSMnZeSoaMlQecm1ghuaY9gLztYGav+JhAiotogWnJlaGiIHj164MMPPxSrSnpBTK1IVYz1tNHfzwn9/ZxQWFKGY7cyse9aKvZdT0Pu4xIcisvAobgMhG+7An9XCwR72SG4pT3ceechEakR0Se007PV5oT2jhFRSM4txM4pndCqkbmodRO9DEEQcDYhG3uvpmL/jVTcy3qstL+JrTGCvcrnabVqZMbeVyKqd1Q6oZ1Uh1ky1VcSiQTt3S3R3t0Sn/VugdjUPOy/noZ911Nx9YEMt9Mf4Xb6I3xz+A6sjfXwmpcdgr3sEOhhBX0dLVWHT0RUIy+cXCUlJcHFxaXa5R88eAAnJ6cXbY6qQTHnigODVI9JJBK0cDBFCwdTvN/DE2myQuy9mooDN9IQfechMh8V4ZczSfjlTBIMdbXQ1dMGQV526NHcFhZGuqoOn4jouV74Hul27drh3XffxdmzZ59aJjc3F9999x28vb3x22+/vWhTVEMcUaGGxM5UH6M7uuHn8QGImfMavhrmhz6tHWGip42C4jLsvZaKGVsuoc1/9mNwZDRWH71TaUFTIqL65IV7rq5fv44vvvgCr732GvT19eHv7w9HR0fo6+sjOzsb169fx7Vr19CmTZtKq6VT7RA4MEgNnKm+Dvq2dkTf1o4oKZPj9N0sHLiRhv3X0/Ag5zHOJGThTEIWFu6JRWNrI0WPlr+rBbS5nhYR1RMvPaH98ePH2L17N44fP47ExEQ8fvwY1tbW8PPzQ0hICLy9vcWKVS3U5oT29l8cQHpeEXa/3xktHc1ErZtIlQRBwM20R9h3rXz48NL9XKX95oY6eKWZLV5tbovuzWxgoq+jokiJSF2p9PE39Gy1mVy1++IAMphckQZIlxUiKjYdB66n4fjtTBSVyhX7tKQSBLhbokcLO7zWwg4uVoYqjJSI1AXvFtRQnNBOmsLWVB/D2rtgWHsXFJaU4fitTETFpiPqRhrS84pw8s5DnLzzEJ/vug4PGyP0aGGHV5vboi2HD4moDjC5UkOc0E6aRF9HC0FedgjysoMgeONasgz7r6fhYGw6rjzIxZ2MfNzJuIvVR+/CVF8b3ZrZokdzW3RrasO7D4moVjC5UivlXVdMrkhTSSQSeDuZwdvJDNNfa4o0WSEOxabjwI10nLidCVlhKf64lIw/LiVDKgHauFjglea2CGphh6Z2xly8lIhEweRKjXBYkEiZnak+hrZ3wdC/hw9P3X2Ig7HpOBibjvvZj3EuMRvnErOx5K84OJrpo3tzW7zazBadPa25eCkRvTBRkqvHjx9DEAQYGpZPHE1MTMT27dvh5eWF4OBgMZqgauCdCURPp6+jhe7NbNG9mS0W9ANup+fhwI3yROtcQhaScwux8XQSNp5Ogo6WBIEe1ni1mQ1eaW4LVys++5CIqk+U5Kpfv34YMGAAJk6ciJycHAQEBEBHRweZmZlYtmwZJk2aJEYz9BwVN35yZIPo+ZrYmqCJrQkmdvNA7uMSHL2ZgUOx6ThyMwMP84tx9GYGjt7MwLw/rqOxtRFe+XuZh3ZuluzVIqJnEiW5iomJwfLlywEAW7duhZ2dHS5cuIDffvsNc+bMYXJVx5hbEdWMmYEO+rR2RJ/WjiiTC7h8PweH4jJwOC4dl+/n4m5mPu4ej8ea4/Ew1NVCRw+rv3vBbNDIgks9EJEyUZKrgoICmJiYAAD27duHAQMGQCqVokOHDkhMTBSjCaqGimFB9lwRvTgtqQR+Lhbwc7FA2GtNkZFXhMNx6Tj8d09WXmEpDtwonyQPAJ62xujezAbdm9mirZsF9LTZq0Wk6URJrpo0aYIdO3bgjTfewF9//YXp06cDANLT00VfKJOe7p/lYJldEYnFxkQPb7Z1xpttnVFaJkdMUg4Ox6XjUFwGbqTIcCv9EW6lP8J3x+JhpKuFQA9rdGtmg+5NbeBsyV4tIk0kSnI1Z84cDB8+HNOnT8err76KwMBAAOW9WH5+fmI0QdXAxfaJape2lhTt3S3R3t0SH/dsjnRZIQ7HZeDwzXQcvZmJR0WlOHAjDQdupAEAGtsYoXtTW3RrZoMAd87VItIUoixVPGjQICQlJeHcuXP466+/FNt79OihmItVUytXroSbmxv09fUREBCAM2fOPLP8li1b0Lx5c+jr68PHxwd79uxR2r9t2zYEBwfDysoKEokEFy9eVNqflZWFqVOnolmzZjAwMICLiwvef/995OYqP8NMIpFU+tq0adMLnaPYOCxIVLdsTfUxuJ0zvhnhj8tzg/Hru4F4r7sHWjqW99jfzcjHDyfiMfqHM2g9fx9G/3AGPxyPx52MR/xjiEiNibbOlb29Pezt7SEIAgRBgEQiQfv27V+ors2bNyMsLAyRkZEICAjAihUrEBISgri4ONja2lYqf/LkSQwbNgwRERHo3bs3Nm7ciP79+yMmJkbx4Oj8/Hx07twZgwcPxoQJEyrVkZycjOTkZPz3v/+Fl5cXEhMTMXHiRCQnJ2Pr1q1KZdeuXYuePXsqXpubm7/QeYpOsc4VEdU1qVSi1KuV+agIh+PK52kdvZWBnIISHLmZgSM3M4BdgJO5Abo2tUG3ptbo1MSaD5smUiOiPbh5zZo1WL58OW7dugUA8PT0xAcffIC33367xnUFBASgXbt2+N///gcAkMvlcHZ2xtSpUzFr1qxK5YcMGYL8/Hzs2rVLsa1Dhw7w9fVFZGSkUtmEhAS4u7vjwoUL8PX1fWYcW7ZswVtvvYX8/Hxoa5fnoRKJBNu3b0f//v1rfF5A7T642WfuX8grKsWhGd3hbs11eYjqC7lcwOUHuTh6s/wOxAv3cvDkb16pBPB3tUBXTxt0bWoDHyczSKX8M4moPqnzBzfPmTMHy5Ytw9SpUxXzraKjozF9+nQkJSVhwYIF1a6ruLgY58+fR3h4uGKbVCpFUFAQoqOjqzwmOjoaYWFhSttCQkKwY8eOmp/MEyrewIrEqsLkyZPx9ttvo3Hjxpg4cSLGjh371MdmFBUVoaioSPFaJpO9VEzPohgWrLUWiOhFSKUS+Dqbw9fZHO/38ETu4xIcv5Wp6NVKyS3E2YRsnE3IxtL9N2FuqIPOTazR1dMGXZpaw8HMQNWnQEQ1IEpytWrVKnz33XcYNmyYYlvfvn3RqlUrTJ06tUbJVWZmJsrKymBnZ6e03c7ODrGxsVUek5qaWmX51NTUGpxF5Tg+//xzvPPOO0rbFyxYgFdffRWGhobYt28f3nvvPTx69Ajvv/9+lfVERERg/vz5LxxHTXAOB1HDYGagg16tHNCrlQMEQcCdjEc4HJeBY7cyceruQ+QUlGDX5RTsupwCAGhia4wunuXJVofGVjDQ5cR4ovpMlOSqpKQEbdu2rbTd398fpaWlYjRRp2QyGXr16gUvLy/MmzdPad9nn32m+L+fnx/y8/OxZMmSpyZX4eHhSr1qMpkMzs7OtRI3J7QTNTwSiUSxWvzbXRqjqLQMZ+KzcPRmebIVm5qH2+mPcDv9EdaeSICOlgRtXS3Rpak1ujSxgZejKbQ4hEhUr4iSXI0cORKrVq3CsmXLlLavXr0aI0aMqFFd1tbW0NLSQlpamtL2tLQ02NvbV3mMvb19jco/S15eHnr27AkTExNs374dOjrPnmQaEBCAzz//HEVFRdDT06u0X09Pr8rttYEPbiZq+PS0tdDF0wZdPG0AAOl5hTh2MxPHbmXg+O1MZD4qRvTdh4i++xCLEQdLI1109LBCF09rdPa0gZM5hxCJVO2Fk6sne2MkEgm+//577Nu3Dx06dAAAnD59GklJSRg1alSN6tXV1YW/vz+ioqIUk8blcjmioqIwZcqUKo8JDAxEVFQUPvjgA8W2/fv3K+Z/VZdMJkNISAj09PSwc+dO6OvrP/eYixcvwsLCos4SqOpgzxWR+rA10cdA/0YY6N8IgiDgRkoejt8u79U6fTcLWfnFSkOIja2N0Nmz/A7EDo2tYGbAuxCJ6toLJ1cXLlxQeu3v7w8AuHPnDoDyHihra2tcu3atxnWHhYVh9OjRaNu2Ldq3b48VK1YgPz8fY8eOBQCMGjUKTk5OiIiIAABMmzYN3bp1w9KlS9GrVy9s2rQJ586dw+rVqxV1ZmVlISkpCcnJyQCAuLg4AP8sISGTyRAcHIyCggKsX78eMplMMfncxsYGWlpa+OOPP5CWloYOHTpAX18f+/fvx8KFCzFjxowan2NtEMA5V0TqTCKRwMvRFF6Opninq4diCPH4rUwcu5WJ6ymy8ucgZubjp+hEaEklaN3IDJ2alCdbbVwsoKstyvKGRPQML5xcHTp0SMw4lAwZMgQZGRmYM2cOUlNT4evri7179yomrSclJUEq/ecXRMeOHbFx40bMnj0bn3zyCTw9PbFjxw7FGlcAsHPnTkVyBgBDhw4FAMydOxfz5s1DTEwMTp8+DaD8cT5Pio+Ph5ubG3R0dLBy5UpMnz4dgiCgSZMmWLZsWZXrZqkC57MTaZYnhxDDATx8VIQTdx7i+K0MnLj9EA9yHiMmKQcxSTn4+uBtGOhoob27JTr/nWy1cDB56p3ORPTiRFvnCgCuX7+OpKQkFBcX/9OARII+ffqI1USDV5vrXDWd/SeKS+U4PvMVNLLgM82INN3djEc4fjsTx29lIvruQ+QVKt9gZGmkW96r5WGFTk2s0cjCgMkW0VPU+TpXd+/exRtvvIErV65AIpEolgSo+CEtKysToxl6nooJ7fzlSEQAGtsYo7GNMUYFuqG0TI4rD3Jx4nYmTtx+iPNJ2cjKL8Yfl5Lxx6Xy6RLOlgbo5GGNjk2s0dHDCtbG9WcuKVFDIkpyNW3aNLi7uyMqKgru7u44c+YMHj58iA8//BD//e9/xWiCaoCpFRH9m7aWFH4uFvBzscCUVz3xuLgM5xKzcPx2Jk7efogrD3JxL+sxNmXdw6az9wAAze1NEOhhhU4e1mjnbsnJ8UTVJEpyFR0djYMHD8La2hpSqRRSqRSdO3dGREQE3n///UqT36l2VExoZ8cVET2Pga7ykg+5BSWIvlveq3XiTibuZuQjNjUPsal5WHsiAVpSCbydzNDRwwodPazQzs0S+jpczJSoKqIkV2VlZTAxMQFQfpdgcnIymjVrBldXV8VdeVT7OKGdiF6UmaEOeno7oKe3AwAgTVaIk3fKk62TtzORnFuIS/dycOleDlYdvgNtqQRtXC3Q0cMKgY2t4OtiDj1tJltEgEjJlbe3Ny5dugR3d3cEBARg8eLF0NXVxerVq9G4cWMxmqBq+OfZguy6IqKXY2eqjzf8GuENv0YAgITMfJy88xAn75Q/oifzUTHOxGfhTHwWVuAW9HWkaOtqiUAPK3RobAUfJzMu+0AaS5Tkavbs2cjPzwdQ/uy93r17o0uXLrCyssLmzZvFaIKq4Z8bCVQcCBGpHTdrI7hZG2F4gAvkcgG30h8h+k75XYjRdx5CVlhafmfi7UwAgJGuFvzdLNHxiWSLj+khTSHqUgxPysrKgoWFBe9c+5faXIrBPXw3BAE480kP2Jo+f3V5IiIxCIKAa8kyRN8pfyzP2fgs5BUpL/tgpFu+xlZFz5aXgym0tdizRQ1HnS/FUBVLS8vaqpqeh/ksEdUhiaR8sru3kxkmdG2M0jI5rv6dbJ26+xDnErKQX1yGQ3EZOBSXAQAw0ddGOzdLBLhbokNjK3izZ4vUSK0lV1T3+OBmIqoPtLWk8HU2h6+zOSZ190DJ32tsnbr7EKfuZpX3bBWW4mBsOg7GpgP4p2croLEVAtwt4e1kBh32bFEDxeSKiIhqlY6WFG1cLNDGxQLvdQdKy+S4dL8i2XqImMTsSj1bFXO2AtzLv1o7mzPZogaDyZWaeHLqHKe5EVF9pq0lhb+rBfxdLTD5lSYoLpXjWnIuTt3Nwun48jlb+cVlOHozA0dvlidbetrlxwS4W6G9uyX8XMy5zhbVW0yu1MSTtyUwtyKihkRX+5/V4yd194BcLuDKg1ycjn+I03ezcDYhC7LC0r+XgnioOKZ1IzO0c7NEe3dLtHOzhJEeP9KofuB3ohriHZpE1JBJpRK0djZHa2dzvNPVA6VlcsSl5eFMfBZO383CmYQsZOUX42xCNs4mZOObw3cgkQDejmYIcC9Pttq6WcLSSFfVp0IaismVmnhyPQ2mVkSkTrS1pGjpaIaWjmYY28kdgiDgdvojnI7Pwun48gnyqbJCXHmQiysPcvH98XhIJICnrbFSz5aDmT7/+KQ6weRKTdTScmVERPWORCKBp50JPO1M8FYHVwiCgHtZj3EmoTzROpOQhfjMfNxMe4SbaY+w4XQSAMDJ3ECRaLV3t4CHjTGTLaoVTK7UhFLPFX9XEJEGkUgkcLEyhIuVIQb5lz+uJz2vEOcSssuHEuOzcCNFhgc5j7H9wgNsv/AAAGBhqIO2bpZo52aBtm6WaOloyucjkiiYXKkJ5QntzK6ISLPZmujjdR8HvO5T/iDqvMISnEvMxtn48gnyl+/nIrugBPuvp2H/9TQAgIGOFlr9PUm+rVv53Ywm+jqqPA1qoJhcqSPmVkRESkz0dfBKM1u80swWAFBYUoarD3L/nhSfhXN/35FYMY8LKB8FaG5virauFmj7d++WI+dtUTUwuVITArjOFRFRdenraKGtW/ldhZNQvvzDzfQ8nE3IxrmELJxLyMaDnMe4kSLDjRQZfj6VCABwMNOHv6vF3wmXJVo4mPKxPVQJkys1wfnsREQvTiqVoLm9KZrbm2JkB1cAQEruY5xLyMb5xPLerespMqTkFmLX5RTsupwCoHwosY2rOfxdLeHvagHfRuYwM+RQoqZjcqWG+DcUEdHLczAzQJ/WBujT2hEAkF9UipikbEXCdfFeDh4VleLE7Yc4cbt8cVOpBPC0NUGbv1eg93e1gJuVIYcSNUy9fVDTypUr4ebmBn19fQQEBODMmTPPLL9lyxY0b94c+vr68PHxwZ49e5T2b9u2DcHBwbCysoJEIsHFixcr1VFYWIjJkyfDysoKxsbGGDhwINLS0pTKJCUloVevXjA0NIStrS0++ugjlJaWvvT5viylCe38ISYiEp2Rnja6eNpg+mtNsf7tAFyc8xp2v98Zn/driX6+jnAyN4BcAOLS8vDLmSTM2HIJr/z3MNr+5wDe/vEcVh2+gzPxWSgoVv1nBtWuetlztXnzZoSFhSEyMhIBAQFYsWIFQkJCEBcXB1tb20rlT548iWHDhiEiIgK9e/fGxo0b0b9/f8TExMDb2xsAkJ+fj86dO2Pw4MGYMGFCle1Onz4du3fvxpYtW2BmZoYpU6ZgwIABOHHiBACgrKwMvXr1gr29PU6ePImUlBSMGjUKOjo6WLhwYe29IdWgNOdKhXEQEWmKJxc3HRnoBgBIkxUqerbOJ2XjRrIMD/OLceBGGg7cKP9jXUdLghYOpmjjYgEfJzN08LCCk7mBCs+ExCYR6uHqkwEBAWjXrh3+97//AQDkcjmcnZ0xdepUzJo1q1L5IUOGID8/H7t27VJs69ChA3x9fREZGalUNiEhAe7u7rhw4QJ8fX0V23Nzc2FjY4ONGzdi0KBBAIDY2Fi0aNEC0dHR6NChA/7880/07t0bycnJsLOzAwBERkZi5syZyMjIgK7u8x+1IJPJYGZmhtzcXJiamtb4vXmaguJSeM35CwBwfUEIDHXrZd5MRKRRHheX4WpyLs4nZiMmMRsxSdnIfFRcqZydqR7aulqiVSMztHe3RDN7E/4er2dq8vld765ccXExzp8/j/DwcMU2qVSKoKAgREdHV3lMdHQ0wsLClLaFhIRgx44d1W73/PnzKCkpQVBQkGJb8+bN4eLiokiuoqOj4ePjo0isKtqZNGkSrl27Bj8/v0r1FhUVoaioSPFaJpNVO6aaqH8pMhERGehqoZ1b+arwQPnTNJKyCnA+MRuX7+fi0v0cXL6fizRZEXZfScHuK+UT5bWlEvg0MkMrJzO0dbOEr7M5nC0NVXkqVAP1LrnKzMxEWVmZUgIDAHZ2doiNja3ymNTU1CrLp6amVrvd1NRU6Orqwtzc/Kn1PK2din1ViYiIwPz586sdx4tSfrYgBwaJiOojiUQCVysjuFoZYUCb8tXk8wpLcCMlD6fvPsTV5Fycic9CdkEJLiTl4EJSDn6MLl8GwtpYF21cLNDa2RxtXS3QwtEUplzktF6qd8mVugkPD1fqVZPJZHB2dha9nSdHdzmfnYio4TDR10F79/IHTFdIeliAmKTsv3u4cnAtWYbMR8XYdz0N+/5eUV4qAVo4mMLb0QwBjcvX3GpqZ8J1t+qBepdcWVtbQ0tLq9JdemlpabC3t6/yGHt7+xqVf1odxcXFyMnJUeq9erIee3v7SnctVrT7tLb09PSgp6dX7TheFEcFiYjUR8WzEvv7OQEoXwYiNjUPZ+KzcDU5F+cSspAmK8K1ZBmuJcuw+dw9AICJvjb8XCzQwsEEAe6WaNXIHNbGtf8ZRMrqXXKlq6sLf39/REVFoX///gDKJ7RHRUVhypQpVR4TGBiIqKgofPDBB4pt+/fvR2BgYLXb9ff3h46ODqKiojBw4EAAQFxcHJKSkhT1BAYG4osvvkB6errirsX9+/fD1NQUXl5eL3C2tYM9V0RE6sVIT1uxblaF1NxCXEjKxoV7ObiQlI0bKXnIKyzF0ZsZOHozA98euQsAsDfVR0BjS3g5mMLPxQKtnc34gOpaVu+SKwAICwvD6NGj0bZtW7Rv3x4rVqxAfn4+xo4dCwAYNWoUnJycEBERAQCYNm0aunXrhqVLl6JXr17YtGkTzp07h9WrVyvqzMrKQlJSEpKTkwGUJ05AeY+Tvb09zMzMMH78eISFhcHS0hKmpqaYOnUqAgMD0aFDBwBAcHAwvLy8MHLkSCxevBipqamYPXs2Jk+eXCe9U8/CCe1ERJrF3kwfoT4OCP374dRlcgHXknNx9YEMl+/n4HR8FhIe5iNVVojfLybj94vln3/aUgk8bIzR1q18/paXgykf4yOyeplcDRkyBBkZGZgzZw5SU1Ph6+uLvXv3KiaPJyUlQSr9Z/3Tjh07YuPGjZg9ezY++eQTeHp6YseOHYo1rgBg586diuQMAIYOHQoAmDt3LubNmwcAWL58OaRSKQYOHIiioiKEhITgm2++URyjpaWFXbt2YdKkSQgMDISRkRFGjx6NBQsW1ObbUT1PLiLKCe1ERBpHSypBq0bmaNXIHMMDXACUL9NzISkH5xOzcSNFhrMJWch8VIy4tDzEpeVhw+kkAICJnjZaOJiiVSOzv5+ZaAIXS64s/6Lq5TpX6qy21rnKKSiG74L9AIBbX4RCR6veLr5PREQqIggCUnILceVB+V2JN1JkuHw/F4+KKq8ab22sC19nczS3N0U7d0t4O5rCSoPnbzXoda7oxUgkEjSxNS7/v4pjISKi+kkikcDR3ACO5gYIaVl+I1Zpmby8Jys1D2cTyu9OvJX+CJmPinHgRjoO3EgHDpUf72Cmj7ZulmhubwI/Z3P4uphzsdMqsOeqjtVWzxUREZFYikrLcPl+Lm6kyHDxXvmwYlJWQaX5vRIJ4GJpCH/X8kf5eDmYorWzOfR11G/CfE0+v5lc1TEmV0RE1BDlF5Xi4r0cXLyXg9jUPJyNz0KqrLBSOW2pBO7WRmjVyBxtXM3R3N4ELR3NGnzCxeSqHmNyRURE6uLhoyJcTZbhXEIWYlPzcOleDtLziiqV09WSoqm9MZrZmaKdmwVaOJiimb1Jg0q4mFzVY0yuiIhIXQmCgOTcQsSlynA+MRtXHshwPVmGzEeVEy4dLQm8HM3Q1NYYvi7maN3IHB42xjDQrZ8JF5OreozJFRERaZKKh1XfSMnD1Qe5iEnKRlxqHh7mF1cqqyWVwMvB9O+1t0zg72oJTzvjetHDxeSqHmNyRUREmk4QBNzPfoxL93MQl5qnWIcru6CkUlmJBGhmZ4LWjczhaWcMPxdzeDmY1XkPF5OreozJFRERUdUe5DzGhaRs3EzNw6X7ubh0Pwc5T0m43KyM4OdsjmZ/T5j3dTGHsV7tLQvBda6IiIiowXEyN4CTuQHQqvy1IAjIyCtCTFI2rqfk4XqyDBeSsvEwvxjxmfmIz8yvdHyrRmbwd7XAWx1cVTacyJ6rOsaeKyIiopfz8FERrjzIxcV75QueXkzKwYOcx4r9JnrauDg3WNTnJbLnioiIiNSWlbEeujezRfdmtopt2fnlz0yMScpGYYlcpQ+iZnJFREREDZ6FkS46NLZCh8ZWqg4FfLovERERkYiYXBERERGJiMkVERERkYiYXBERERGJiBPa61jFyhcymUzFkRAREVF1VXxuV2cFKyZXdSwvLw8A4OzsrOJIiIiIqKby8vJgZmb2zDJcRLSOyeVyJCcnw8TEBBKJuGtwyGQyODs74969e2q5QKm6nx+g/ufI82v41P0c1f38APU/x9o6P0EQkJeXB0dHR0ilz55VxZ6rOiaVStGoUaNabcPU1FQtf2AqqPv5Aep/jjy/hk/dz1Hdzw9Q/3OsjfN7Xo9VBU5oJyIiIhIRkysiIiIiETG5UiN6enqYO3cu9PT0VB1KrVD38wPU/xx5fg2fup+jup8foP7nWB/OjxPaiYiIiETEnisiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhIRkysiIiIiETG5IiIiIhKRtqoD0DRyuRzJyckwMTGBRCJRdThERERUDYIgIC8vD46OjpBKn903xeSqjiUnJ8PZ2VnVYRAREdELuHfvHho1avTMMkyu6piJiQmA8otjamqq4miIiIioOmQyGZydnRWf48/C5KqOVQwFmpqaMrkiIiJqYKozpYcT2omIiIhExOSKiIiISERMroiIiIhExDlXRERieXgdOL0QKCl4YqNQzYOfNY+jvtRR3XpeNo6XOb622q7NY1/kuLps62nH1Ze2nihfMR9KxxgI/fE57dQeJldERGK58D/gxgZVR0FE+pYqbZ7JFRGRWEoLy/9t0h9w6/kSFQl4fg+SptTxMser6tgXPf5F26yPbT3R0yQI//QoKVSjV0vpuL//fbIeQVA+RnjiWC29asRYe5hcERGJzaED0PpdVUdBRCrCCe1ERKKp+MuZj7Yi0mRMroiIRFPdSeNEpM6YXBERiaVizgcfyk6k0ZhcERGJhsOCRMTkiohIfOy5ItJoTK6IiMQicM4VETG5IiISEYcFiYjJFRGReDihnYjA5IqIqBYwuSLSZEyuiIhEw54rImJyRUQkHk5oJyIwuaqxiIgItGvXDiYmJrC1tUX//v0RFxen6rCIqF7ghHYiYnJVY0eOHMHkyZNx6tQp7N+/HyUlJQgODkZ+fr6qQyOieoPJFZEm01Z1AA3N3r17lV6vW7cOtra2OH/+PLp27aqiqIioXuDdgkQEJlcvLTc3FwBgaWlZ5f6ioiIUFRUpXstksjqJi4hUgXOuiIjDgi9FLpfjgw8+QKdOneDt7V1lmYiICJiZmSm+nJ2d6zhKIqo7nHNFREyuXsrkyZNx9epVbNq06allwsPDkZubq/i6d+9eHUZIRCrBYUEijcZhwRc0ZcoU7Nq1C0ePHkWjRo2eWk5PTw96enp1GBkRqYzAnisiYnJVY4IgYOrUqdi+fTsOHz4Md3d3VYdERPUG51wREZOrGps8eTI2btyI33//HSYmJkhNTQUAmJmZwcDAQMXREZFK8W5BIgLnXNXYqlWrkJubi+7du8PBwUHxtXnzZlWHRkT1BpMrIk3GnqsaEvh4CyJ6KvZcERF7roiIxMM/vogITK6IiETEuwWJiMkVEVEtYHJFpMmYXBERiYV3CxIRmFwREYmIc66IiMkVEZGIOOeKiJhcERGJj8OCRBqNyRURkVj4bEEiApMrIiIRcc4VETG5IiISD+8WJCIwuSIiEhGHBYmIyRURkfjYc0Wk0ZhcERGJhc8WJCIwuSIiEhGHBYkI0FZ1AHWlpKQEqampKCgogI2NDSwtLVUdEhGpG05oJyKoec9VXl4eVq1ahW7dusHU1BRubm5o0aIFbGxs4OrqigkTJuDs2bOqDpOI1A6TKyJNprbJ1bJly+Dm5oa1a9ciKCgIO3bswMWLF3Hz5k1ER0dj7ty5KC0tRXBwMHr27Ilbt26pOmQiavA4LEhEajwsePbsWRw9ehQtW7ascn/79u0xbtw4REZGYu3atTh27Bg8PT3rOEoiUi+c0E5Eapxc/fLLL9Uqp6enh4kTJ9ZyNESkETjnioigxsOCRESqw+SKSJOpbc/VgAEDql1227ZttRgJEWkO9lwRkRr3XJmZmSm+TE1NERUVhXPnzin2nz9/HlFRUTAzM1NhlESkVriIKBFBjXuu1q5dq/j/zJkzMXjwYERGRkJLSwsAUFZWhvfeew+mpqaqCpGI1A7vFiQiNe65etIPP/yAGTNmKBIrANDS0kJYWBh++OEHFUZGRGqJw4JEGk0jkqvS0lLExsZW2h4bGwu5XK6CiIhILQnsuSIiNR4WfNLYsWMxfvx43LlzB+3btwcAnD59GosWLcLYsWNVHB0RqQ/OuSIiDUmu/vvf/8Le3h5Lly5FSkoKAMDBwQEfffQRPvzwQxVHR0Rqg+tcERE0JLmSSqX4+OOP8fHHH0MmkwEAJ7ITUS1ickWkyTQiuXoSkyoiqj2cc0VEGpRcbd26Fb/++iuSkpJQXFystC8mJkZFURGRWuE6V0QEDblb8KuvvsLYsWNhZ2eHCxcuoH379rCyssLdu3cRGhqq6vCISG1wzhURaUhy9c0332D16tX4+uuvoauri48//hj79+/H+++/j9zcXFWHR0Rqh8kVkSbTiOQqKSkJHTt2BAAYGBggLy8PADBy5Ej88ssvqgyNiNQKe66ISEOSK3t7e2RlZQEAXFxccOrUKQBAfHw8BM6RICKx8PcJEUFDkqtXX30VO3fuBFC+oOj06dPx2muvYciQIXjjjTdUHB0RqQ/eLUhEGnK34OrVqxWPuZk8eTKsrKxw8uRJ9O3bF++++66KoyMitcNhQSKNpvbJVWlpKRYuXIhx48ahUaNGAIChQ4di6NChKo6MiNQOny1IRNCAYUFtbW0sXrwYpaWlqg6FiNQe51wRkQYkVwDQo0cPHDlyRJS6jh49ij59+sDR0RESiQQ7duwQpV4iUgN8tiARQQOGBQEgNDQUs2bNwpUrV+Dv7w8jIyOl/X379q12Xfn5+WjdujXGjRuHAQMGiB0qETVoHBYkIg1Jrt577z0AwLJlyyrtk0gkKCsrq3ZdoaGhXNWdiJ6NPVdEGk0jkquKOwVVoaioCEVFRYrXMplMZbEQUS3jOldEBA2Zc6VKERERMDMzU3w5OzurOiQiqjUcFiQiNU6uNm3aVO2y9+7dw4kTJ2oljvDwcOTm5iq+7t27VyvtEFF9wOSKiNQ4uVq1ahVatGiBxYsX48aNG5X25+bmYs+ePRg+fDjatGmDhw8f1kocenp6MDU1VfoiIjXHOVdEGk1t51wdOXIEO3fuxNdff43w8HAYGRnBzs4O+vr6yM7ORmpqKqytrTFmzBhcvXoVdnZ2qg6ZiBo6LiJKRFDj5AooX2Khb9++yMzMxPHjx5GYmIjHjx/D2toafn5+8PPzg1Ras867R48e4fbt24rX8fHxuHjxIiwtLeHi4iL2KRBRg8IJ7USk5slVBWtra/Tv31+Uus6dO4dXXnlF8TosLAwAMHr0aKxbt06UNoiogeIiokQEDUmuxNS9e3cIvN2aiJ6JyRWRJlPbCe1ERHWPPVdExOSKiEg87NUmIjC5IiISEe8WJCINS66Ki4sRFxeH0tJSVYdCROqMw4JEGk0jkquCggKMHz8ehoaGaNmyJZKSkgAAU6dOxaJFi1QcHRGpDa5zRUTQkOQqPDwcly5dwuHDh6Gvr6/YHhQUhM2bN6swMiJSL5xzRUQashTDjh07sHnzZnTo0AGSJ7rrW7ZsiTt37qgwMiJSL+y5IiIN6bnKyMiAra1tpe35+flKyRYRkSj4e4VIo2lEctW2bVvs3r1b8boiofr+++8RGBioqrCISN1wzhURQUOGBRcuXIjQ0FBcv34dpaWl+PLLL3H9+nWcPHkSR44cUXV4RKQ2OOeKiDSk56pz5864ePEiSktL4ePjg3379sHW1hbR0dHw9/dXdXhEpC74bEEigob0XAGAh4cHvvvuO1WHQUQagckVkSbTiJ6rmJgYXLlyRfH6999/R//+/fHJJ5+guLhYhZERkXphzxURaUhy9e677+LmzZsAgLt372LIkCEwNDTEli1b8PHHH6s4OiJSG3y2IBFBQ5KrmzdvwtfXFwCwZcsWdOvWDRs3bsS6devw22+/qTY4IlIjvFuQiDQkuRIEAXK5HABw4MABvP766wAAZ2dnZGZmqjI0IlJHHBYk0mgakVy1bdsW//nPf/Dzzz/jyJEj6NWrFwAgPj4ednZ2Ko6OiNQG17kiImhIcrVixQrExMRgypQp+PTTT9GkSRMAwNatW9GxY0cVR0dE6oNzrohIQ5ZiaNWqldLdghWWLFkCLS0tFURERGqJ61wRETQkuXoafX19VYdARGqJyRWRJtOI5KqsrAzLly/Hr7/+iqSkpEprW2VlZakoMiJSL5xzRUQaMudq/vz5WLZsGYYMGYLc3FyEhYVhwIABkEqlmDdvnqrDIyK1wTlXRKQhydWGDRvw3Xff4cMPP4S2tjaGDRuG77//HnPmzMGpU6dUHR4RqQvOuSIiaEhylZqaCh8fHwCAsbExcnNzAQC9e/fG7t27VRkaEakVDgsSkYYkV40aNUJKSgqA8gc479u3DwBw9uxZ6OnpqTI0IlJH7Lki0mgakVy98cYbiIqKAgBMnToVn332GTw9PTFq1CiMGzdOxdERkdrgswWJCBpyt+CiRYsU/x8yZAhcXFwQHR0NT09P9OnTR4WREZF64bAgEWlIcvVvgYGBCAwMVHUYRKRuOKGdiKBBydWtW7dw6NAhpKenKx7iXGHOnDkqioqI1BOTKyJNphHJ1XfffYdJkybB2toa9vb2kDzxV6VEImFyRUQiYc8VEWlIcvWf//wHX3zxBWbOnKnqUIhInXFCOxFBQ+4WzM7OxptvvqnqMIhI7XFCOxFpSHL15ptvKta2IiKqfUyuiDSZ2g4LfvXVV4r/N2nSBJ999hlOnToFHx8f6OjoKJV9//336zo8IlJLnHNFRIBEENRzkoC7u3u1ykkkEty9e7eWo/mHTCaDmZkZcnNzYWpqWmftElEd+J8FUJQDjI0FLJupOhoiElFNPr/VtucqPj5e1SEQkcbhnCsi0pA5V08SBAFq2llHRPUFhwWJNJrGJFdr1qyBt7c39PX1oa+vD29vb3z//feqDouI1InAnisi0pDkas6cOZg2bRr69OmDLVu2YMuWLejTpw+mT5/+wguIrly5Em5ubtDX10dAQADOnDkjctRE1PCwV5yI1HjO1ZNWrVqF7777DsOGDVNs69u3L1q1aoWpU6diwYIFNapv8+bNCAsLQ2RkJAICArBixQqEhIQgLi4Otra2YodPRA0Fny1IRNCQnquSkhK0bdu20nZ/f3+UlpbWuL5ly5ZhwoQJGDt2LLy8vBAZGQlDQ0P88MMPlcoWFRVBJpMpfRGRumNyRaTJNCK5GjlyJFatWlVp++rVqzFixIga1VVcXIzz588jKChIsU0qlSIoKAjR0dGVykdERMDMzEzx5ezsXPMTIKIGgj1XRKQhw4JA+YT2ffv2oUOHDgCA06dPIykpCaNGjUJYWJii3LJly55ZT2ZmJsrKymBnZ6e03c7ODrGxsZXKh4eHK9Uvk8mYYBGpK96JTETQkOTq6tWraNOmDQDgzp07AABra2tYW1vj6tWrinKSWvhrU09PD3p6eqLXS0T1Ee8WJCINSa4OHTokWl3W1tbQ0tJCWlqa0va0tDTY29uL1g4RNWAcFiTSaBox50pMurq68Pf3R1RUlGKbXC5HVFQUAgMDVRgZEakee66ISEN6rsQWFhaG0aNHo23btmjfvj1WrFiB/Px8jB07VtWhEZEqcc4VEYHJ1QsZMmQIMjIyMGfOHKSmpsLX1xd79+6tNMmdiDQNe66IiMnVC5syZQqmTJmi6jCIqD7inCsijcY5V0REYuGzBYkITK6IiETEOVdExOSKiEg8fLYgEYHJFRGRiDgsSERMroiIxMeeKyKNxuSKiEgsXOeKiMDkiohIRBwWJCImV0RE4uOwIJFGY3JFRCQ6JldEmozJFRGRGJTmWzG5ItJkTK6IiETByexEVI7JFRGRGJ7sueKcKyKNxuSKiEh0TK6INBmTKyIiUbDniojKMbkiIhIDFxAlor8xuSIiEgXvFiSickyuiIjExmFBIo3G5IqISAxc54qI/sbkiohIFJxzRUTlmFwREYmB61wR0d+YXBERiY7JFZEmY3JFRCQK9lwRUTkmV0REouCcKyIqx+SKiEgMvFuQiP7G5IqISHRMrog0GZMrIiJRcM4VEZVjckVEJAY+W5CI/sbkiohIFJxzRUTlmFwREYmNw4JEGk1b1QGQSAqzgW+dVB0FkQZjzxURlWNypU5KH6s6AiKyaQ1I+auVSJPxN4C60DUF3o6vep8YQxSCUPtDHQ29jYZWd0OoU8z66qouYycOCxJpOCZX6kKqBZi5qToKIiIijccJ7UREREQiYnJFREREJCImV0REREQiYnJFREREJCJOaK9jwt+PyJDJZCqOhIiIiKqr4nNbqMajrphc1bG8vDwAgLOzs4ojISIioprKy8uDmZnZM8tIhOqkYCQauVyO5ORkmJiYQCLyWjgymQzOzs64d+8eTE1NRa27PlD38wPU/xx5fg2fup+jup8foP7nWFvnJwgC8vLy4OjoCKn02bOq2HNVx6RSKRo1alSrbZiamqrlD0wFdT8/QP3PkefX8Kn7Oar7+QHqf461cX7P67GqwAntRERERCJickVEREQkIiZXakRPTw9z586Fnp6eqkOpFep+foD6nyPPr+FT93NU9/MD1P8c68P5cUI7ERERkYjYc0VEREQkIiZXRERERCJickVEREQkIiZXRERERCJicqUmVq5cCTc3N+jr6yMgIABnzpxRdUgvLCIiAu3atYOJiQlsbW3Rv39/xMXFKZXp3r07JBKJ0tfEiRNVFHHNzJs3r1LszZs3V+wvLCzE5MmTYWVlBWNjYwwcOBBpaWkqjLhm3NzcKp2fRCLB5MmTATTMa3f06FH06dMHjo6OkEgk2LFjh9J+QRAwZ84cODg4wMDAAEFBQbh165ZSmaysLIwYMQKmpqYwNzfH+PHj8ejRozo8i6d71vmVlJRg5syZ8PHxgZGRERwdHTFq1CgkJycr1VHVdV+0aFEdn8nTPe8ajhkzplL8PXv2VCrTUK8hgCp/JiUSCZYsWaIoU5+vYXU+F6rzuzMpKQm9evWCoaEhbG1t8dFHH6G0tFT0eJlcqYHNmzcjLCwMc+fORUxMDFq3bo2QkBCkp6erOrQXcuTIEUyePBmnTp3C/v37UVJSguDgYOTn5yuVmzBhAlJSUhRfixcvVlHENdeyZUul2I8fP67YN336dPzxxx/YsmULjhw5guTkZAwYMECF0dbM2bNnlc5t//79AIA333xTUaahXbv8/Hy0bt0aK1eurHL/4sWL8dVXXyEyMhKnT5+GkZERQkJCUFhYqCgzYsQIXLt2Dfv378euXbtw9OhRvPPOO3V1Cs/0rPMrKChATEwMPvvsM8TExGDbtm2Ii4tD3759K5VdsGCB0nWdOnVqXYRfLc+7hgDQs2dPpfh/+eUXpf0N9RoCUDqvlJQU/PDDD5BIJBg4cKBSufp6DavzufC8351lZWXo1asXiouLcfLkSfz4449Yt24d5syZI37AAjV47du3FyZPnqx4XVZWJjg6OgoREREqjEo86enpAgDhyJEjim3dunUTpk2bprqgXsLcuXOF1q1bV7kvJydH0NHREbZs2aLYduPGDQGAEB0dXUcRimvatGmCh4eHIJfLBUFo2NdOEAQBgLB9+3bFa7lcLtjb2wtLlixRbMvJyRH09PSEX375RRAEQbh+/boAQDh79qyizJ9//ilIJBLhwYMHdRZ7dfz7/Kpy5swZAYCQmJio2Obq6iosX768doMTSVXnOHr0aKFfv35PPUbdrmG/fv2EV199VWlbQ7qG//5cqM7vzj179ghSqVRITU1VlFm1apVgamoqFBUViRofe64auOLiYpw/fx5BQUGKbVKpFEFBQYiOjlZhZOLJzc0FAFhaWipt37BhA6ytreHt7Y3w8HAUFBSoIrwXcuvWLTg6OqJx48YYMWIEkpKSAADnz59HSUmJ0vVs3rw5XFxcGuT1LC4uxvr16zFu3DilB5U35Gv3b/Hx8UhNTVW6ZmZmZggICFBcs+joaJibm6Nt27aKMkFBQZBKpTh9+nSdx/yycnNzIZFIYG5urrR90aJFsLKygp+fH5YsWVIrwy216fDhw7C1tUWzZs0wadIkPHz4ULFPna5hWloadu/ejfHjx1fa11Cu4b8/F6rzuzM6Oho+Pj6ws7NTlAkJCYFMJsO1a9dEjY8Pbm7gMjMzUVZWpvTNAgB2dnaIjY1VUVTikcvl+OCDD9CpUyd4e3srtg8fPhyurq5wdHTE5cuXMXPmTMTFxWHbtm0qjLZ6AgICsG7dOjRr1gwpKSmYP38+unTpgqtXryI1NRW6urqVPrTs7OyQmpqqmoBfwo4dO5CTk4MxY8YotjXka1eViutS1c9gxb7U1FTY2toq7dfW1oalpWWDu66FhYWYOXMmhg0bpvRQ3Pfffx9t2rSBpaUlTp48ifDwcKSkpGDZsmUqjLb6evbsiQEDBsDd3R137tzBJ598gtDQUERHR0NLS0utruGPP/4IExOTStMNGso1rOpzoTq/O1NTU6v8Oa3YJyYmV1SvTZ48GVevXlWakwRAaZ6Dj48PHBwc0KNHD9y5cwceHh51HWaNhIaGKv7fqlUrBAQEwNXVFb/++isMDAxUGJn41qxZg9DQUDg6Oiq2NeRrp+lKSkowePBgCIKAVatWKe0LCwtT/L9Vq1bQ1dXFu+++i4iIiAbxmJWhQ4cq/u/j44NWrVrBw8MDhw8fRo8ePVQYmfh++OEHjBgxAvr6+krbG8o1fNrnQn3CYcEGztraGlpaWpXuiEhLS4O9vb2KohLHlClTsGvXLhw6dAiNGjV6ZtmAgAAAwO3bt+siNFGZm5ujadOmuH37Nuzt7VFcXIycnBylMg3xeiYmJuLAgQN4++23n1muIV87AIrr8qyfQXt7+0o3mJSWliIrK6vBXNeKxCoxMRH79+9X6rWqSkBAAEpLS5GQkFA3AYqscePGsLa2VnxfqsM1BIBjx44hLi7uuT+XQP28hk/7XKjO7057e/sqf04r9omJyVUDp6urC39/f0RFRSm2yeVyREVFITAwUIWRvThBEDBlyhRs374dBw8ehLu7+3OPuXjxIgDAwcGhlqMT36NHj3Dnzh04ODjA398fOjo6StczLi4OSUlJDe56rl27Fra2tujVq9czyzXkawcA7u7usLe3V7pmMpkMp0+fVlyzwMBA5OTk4Pz584oyBw8ehFwuVySX9VlFYnXr1i0cOHAAVlZWzz3m4sWLkEqllYbSGor79+/j4cOHiu/Lhn4NK6xZswb+/v5o3br1c8vWp2v4vM+F6vzuDAwMxJUrV5SS5Io/FLy8vEQPmBq4TZs2CXp6esK6deuE69evC++8845gbm6udEdEQzJp0iTBzMxMOHz4sJCSkqL4KigoEARBEG7fvi0sWLBAOHfunBAfHy/8/vvvQuPGjYWuXbuqOPLq+fDDD4XDhw8L8fHxwokTJ4SgoCDB2tpaSE9PFwRBECZOnCi4uLgIBw8eFM6dOycEBgYKgYGBKo66ZsrKygQXFxdh5syZStsb6rXLy8sTLly4IFy4cEEAICxbtky4cOGC4m65RYsWCebm5sLvv/8uXL58WejXr5/g7u4uPH78WFFHz549BT8/P+H06dPC8ePHBU9PT2HYsGGqOiUlzzq/4uJioW/fvkKjRo2EixcvKv1MVtxhdfLkSWH58uXCxYsXhTt37gjr168XbGxshFGjRqn4zP7xrHPMy8sTZsyYIURHRwvx8fHCgQMHhDZt2gienp5CYWGhoo6Geg0r5ObmCoaGhsKqVasqHV/fr+HzPhcE4fm/O0tLSwVvb28hODhYuHjxorB3717BxsZGCA8PFz1eJldq4uuvvxZcXFwEXV1doX379sKpU6dUHdILA1Dl19q1awVBEISkpCSha9eugqWlpaCnpyc0adJE+Oijj4Tc3FzVBl5NQ4YMERwcHARdXV3ByclJGDJkiHD79m3F/sePHwvvvfeeYGFhIRgaGgpvvPGGkJKSosKIa+6vv/4SAAhxcXFK2xvqtTt06FCV35OjR48WBKF8OYbPPvtMsLOzE/T09IQePXpUOveHDx8Kw4YNE4yNjQVTU1Nh7NixQl5engrOprJnnV98fPxTfyYPHTokCIIgnD9/XggICBDMzMwEfX19oUWLFsLChQuVEhNVe9Y5FhQUCMHBwYKNjY2go6MjuLq6ChMmTKj0B2pDvYYVvv32W8HAwEDIycmpdHx9v4bP+1wQhOr97kxISBBCQ0MFAwMDwdraWvjwww+FkpIS0eOV/B00EREREYmAc66IiIiIRMTkioiIiEhETK6IiIiIRMTkioiIiEhETK6IiIiIRMTkioiIiEhETK6IiIiIRMTkioiIiEhETK6ISGMcPnwYEomk0sNd60pUVBRatGiBsrIyUeqbNWsWpk6dKkpdRCQertBORGqpe/fu8PX1xYoVKxTbiouLkZWVBTs7O0gkkjqPyd/fH2FhYRgxYoQo9WVmZqJx48a4ePEiGjduLEqdRPTy2HNFRBpDV1cX9vb2Kkmsjh8/jjt37mDgwIGi1WltbY2QkBCsWrVKtDqJ6OUxuSIitTNmzBgcOXIEX375JSQSCSQSCRISEioNC65btw7m5ubYtWsXmjVrBkNDQwwaNAgFBQX48ccf4ebmBgsLC7z//vtKQ3lFRUWYMWMGnJycYGRkhICAABw+fPiZMW3atAmvvfYa9PX1FdvmzZsHX19f/Pzzz3Bzc4OZmRmGDh2KvLw8RZmtW7fCx8cHBgYGsLKyQlBQEPLz8xX7+/Tpg02bNonzxhGRKJhcEZHa+fLLLxEYGIgJEyYgJSUFKSkpcHZ2rrJsQUEBvvrqK2zatAl79+7F4cOH8cYbb2DPnj3Ys2cPfv75Z3z77bfYunWr4pgpU6YgOjoamzZtwuXLl/Hmm2+iZ8+euHXr1lNjOnbsGNq2bVtp+507d7Bjxw7s2rULu3btwpEjR7Bo0SIAQEpKCoYNG4Zx48bhxo0bOHz4MAYMGIAnZ3O0b98e9+/fR0JCwgu+W0QkNm1VB0BEJDYzMzPo6urC0NAQ9vb2zyxbUlKCVatWwcPDAwAwaNAg/Pzzz0hLS4OxsTG8vLzwyiuv4NChQxgyZAiSkpKwdu1aJCUlwdHREQAwY8YM7N27F2vXrsXChQurbCcxMVFR/klyuRzr1q2DiYkJAGDkyJGIiorCF198gZSUFJSWlmLAgAFwdXUFAPj4+CgdX1FnYmIi3Nzcqv8mEVGtYXJFRBrN0NBQkVgBgJ2dHdzc3GBsbKy0LT09HQBw5coVlJWVoWnTpkr1FBUVwcrK6qntPH78WGlIsIKbm5sisQIABwcHRVutW7dGjx494OPjg5CQEAQHB2PQoEGwsLBQlDcwMABQ3gNHRPUDkysi0mg6OjpKryUSSZXb5HI5AODRo0fQ0tLC+fPnoaWlpVTuyYTs36ytrZGdnV2t9iva0tLSwv79+3Hy5Ens27cPX3/9NT799FOcPn0a7u7uAICsrCwAgI2NTXVOl4jqAOdcEZFa0tXVFW09qSf5+fmhrKwM6enpaNKkidLXs4Yg/fz8cP369Rq3J5FI0KlTJ8yfPx8XLlyArq4utm/frth/9epV6OjooGXLli90PkQkPvZcEZFacnNzw+nTp5GQkABjY2NYWlqKUm/Tpk0xYsQIjBo1CkuXLoWfnx8yMjIQFRWFVq1aoVevXlUeFxISgh9//LFGbZ0+fRpRUVEIDg6Gra0tTp8+jYyMDLRo0UJR5tixY+jSpYtieJCIVI89V0SklmbMmAEtLS14eXnBxsYGSUlJotW9du1ajBo1Ch9++CGaNWuG/v374+zZs3BxcXnqMSNGjMC1a9cQFxdX7XZMTU1x9OhRvP7662jatClmz56NpUuXIjQ0VFFm06ZNmDBhwkudDxGJiyu0ExHVkY8++ggymQzffvutKPX9+eef+PDDD3H58mVoa3Mggqi+YM8VEVEd+fTTT+Hq6qqYsP6y8vPzsXbtWiZWRPUMe66IiIiIRMSeKyIiIiIRMbkiIiIiEhGTKyIiIiIRMbkiIiIiEhGTKyIiIiIRMbkiIiIiEhGTKyIiIiIRMbkiIiIiEhGTKyIiIiIR/T+maCtWwpS53wAAAABJRU5ErkJggg==", "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.12.0" } }, "nbformat": 4, "nbformat_minor": 5 }