{ "cells": [ { "cell_type": "markdown", "id": "aa903e1f", "metadata": {}, "source": [ "# Tutorial for Conductivity module" ] }, { "cell_type": "code", "execution_count": 1, "id": "55008c01", "metadata": {}, "outputs": [], "source": [ "# import the necessary packages along with HavNegpy\n", "import numpy as np\n", "import pandas as pd\n", "import matplotlib.pyplot as plt\n", "import os\n", "import HavNegpy as dd\n", "%matplotlib qt" ] }, { "cell_type": "code", "execution_count": 2, "id": "e27dd328", "metadata": {}, "outputs": [], "source": [ "# extract the data\n", "\n", "filename = 'cond_example_data.txt'\n", "col_names = ['Freq','Sig1','Sig2']\n", "df = pd.read_csv(filename, sep='\\t',index_col=False,usecols = [0,1,2],names=col_names,header=None,skiprows=10,encoding='unicode_escape',engine='python')\n", "\n" ] }, { "cell_type": "markdown", "id": "0018dcc4", "metadata": {}, "source": [ "## Fitting the real part of Conductivity" ] }, { "cell_type": "code", "execution_count": 3, "id": "ec2de83e", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'Example for Conductivity fitting')" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# plot the data \n", "\n", "x = np.log10(df['Freq'])\n", "y = np.log10(df['Sig1'])\n", "plt.scatter(x,y,label='example data')\n", "plt.xlabel('log f [Hz]')\n", "plt.ylabel('log ( $\\sigma´$)')\n", "plt.legend()\n", "plt.title('Example for Conductivity fitting')" ] }, { "attachments": { "data_example-2.png": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAABbp0lEQVR4nO3deVzU1f4/8NeAMAOKgzDCDIlCbohYgqaClWaimJLV/bpcwyTNEjP3sk2BumabersuZF5zSUtbrimZ5G6aCypaIl41wzCFULRBVBaH8/uD38xlnBkYYFbm9Xw85vG485nz+cx7gO68POdzzpEIIQSIiIiIyGW42bsAIiIiIrItBkAiIiIiF8MASERERORiGACJiIiIXAwDIBEREZGLYQAkIiIicjEMgEREREQuhgGQiIiIyMUwABIRERG5GAZAIiIiIhfDAEhERETkYhgAiYiIiFwMAyARERGRi2EAJCIiInIxDIBERERELoYBkIiIiMjFMAASERERuRgGQCIiIiIXwwBIRERE5GIYAImIiIhcDAMgERERkYthACQiIiJyMQyARERERC6GAZCIiIjIxTAAEhEREbkYBkAiIiIiF8MASERERORiGACJiIiIXAwDIBEREZGLYQAkIiIicjEMgEREREQuhgGQiIiIyMUwABIRERG5GAZAIiIiIhfDAEhERETkYhgAiYiIiFwMAyBRNatWrYJEIjH52LNnj71LtIiQkBAkJiZa9JpvvvkmWrdujSZNmsDX19ei1zZl3759GD58OO655x54enpCLpcjJiYGaWlpuHnzpk1qqE7793PhwgWrXP/y5ctISUnBiRMnDF5LSUmBRCKp13Xv/nuo6X0aasOGDejcuTO8vLwgkUhw4sQJo7UvXboUq1atMjjfWj8DIlfTxN4FEDmilStXIiwszOB4eHi4HapxfJs2bcLcuXPxxhtvYNCgQZBKpVZ/z+TkZLz11luIiYnB22+/jbZt2+LWrVs4cOAAUlJScPbsWSxcuNDqddjS5cuXkZqaipCQEHTt2lXvteeeew5xcXH1uu7GjRvRvHlzs96nIa5cuYLRo0cjLi4OS5cuhVQqRYcOHYzWvnTpUigUCoN/qFjrZ0DkahgAiYyIiIhA9+7d7V2G08jOzgYATJ48GQEBARa55q1bt+Dt7W30ta+++gpvvfUWxo0bh+XLl+v1+gwaNAivvPIKDh48aJE6nEWrVq3QqlWrep0bGRlp4WqMO3v2LCoqKpCQkIA+ffrojnt7e9e79uoa8jMgcjmCiHRWrlwpAIgjR47U2O6LL74QAMSiRYv0js+ZM0e4ubmJbdu26Y6lpKSIHj16iBYtWggfHx8RGRkp/v3vf4vKykq9c9u0aSMGDx4s0tPTRdeuXYVMJhNhYWEiPT1dV1tYWJjw9vYWDzzwgEGNY8aMEU2bNhXZ2dmiX79+wtvbWygUCvHiiy+KmzdvGrzXmDFj9I6p1WoxY8YMERISIjw8PERQUJCYMmWKKCkpqfFn0aZNGwFA75GcnCyEEEKj0Yj33ntPdOzYUXh6eoqWLVuK0aNHi4sXL+pdo0+fPqJz585i7969Ijo6Wnh5eYkRI0aYfM+IiAjRokULg89lyu3bt8Wrr76q99kmTpworl+/bvBZBg8eLLZu3SoiIyOFTCYTHTt2FCtWrDC45sGDB0VMTIyQSqVCpVKJV199VXzyyScCgMjNzdW1q/7zuPu97v4d/PHHH2L8+PGiVatWwsPDQ6hUKvG3v/1NFBQUiN27dxv8nKtfOzk5WVT/v/ShQ4eK1q1bC41GY/DePXr0EJGRkUZrqel91qxZIwCIAwcOGFwzNTVVNGnSRFy6dMngNSGq/j7vvmafPn2M1m7sb6pNmzZ1/hlor2Xu73Tfvn2iV69eQiqViqCgIPHmm2+K5cuXG/xOiRoDBkCiarQB8NChQ6KiokLvcefOHb22EyZMEJ6enrogtnPnTuHm5ibefPNNvXaJiYlixYoVYvv27WL79u3i7bffFl5eXiI1NVWvXZs2bUSrVq1ERESE+OKLL8T3338vevbsKTw8PMScOXNE7969xX/+8x+xceNG0aFDBxEYGChu3bqlO3/MmDHC09NTtG7dWsydO1ds27ZNpKSkiCZNmoghQ4YYvFf18HHz5k3RtWtXoVAoxIIFC8SOHTvERx99JORyuejXr59BWK0uKytLjBs3TgAQGRkZ4uDBg7qA9/zzzwsAYtKkSSIjI0N8/PHHomXLliI4OFhcuXJFd40+ffoIPz8/ERwcLBYtWiR2794t9u7da/T9Ll++LADUGBCrq6ysFAMHDhRNmjQRs2fPFtu2bRMffvihaNq0qYiMjBSlpaUGv4Pw8HCxZs0a8cMPP4hhw4YJAHr1nDp1Snh7e4vw8HDxxRdfiE2bNomBAweK1q1b1zsA/vHHH0KlUun9DjZs2CDGjh0rTp8+LdRqte7v88033xQHDx7U+1nfHX42bdokAIjt27frve/p06cFAPGvf/3LaC01vU9ZWZlQKpXi6aef1rtmRUWFCAoKEsOGDTP5e/j111/FkiVLBADxzjvviIMHD4pTp04ZrT0rK0vce++9IjIyUvf+WVlZdf4ZaD+bOb/Tn3/+WchkMnHfffeJ9evXi82bN4vHHntMhISEMABSo8QASFSN9svF2MPd3V2vbWlpqYiMjBShoaEiJydHBAYGij59+hgExeo0Go2oqKgQb731lvD399cLVm3atBFeXl7ijz/+0B07ceKEACBUKpVeb9e3334rAIjNmzfrjml7WD766CO995w7d64AIPbv36/3XtXDx7x584Sbm5tBr+LXX38tAIjvv/++xp+b9ou3eqjTBo2JEyfqtT18+LAAIF5//XXdsT59+ggAYufOnTW+jxBCHDp0SAAQr776aq1thRAiIyNDABDvv/++3vENGzYIAOKTTz7RHWvTpo2QyWTi999/1x27ffu28PPzEy+88ILu2IgRI4SXl5coKCjQHbtz544ICwurdwAcO3as8PDwEDk5OSY/y5EjRwQAsXLlSoPX7g4/FRUVIjAwUIwaNUqv3SuvvCI8PT3F1atXTdZS2/t4enqKP//8U3dM+7M0Fdq1tD14X331VY21CyFE586ddT2E1dXlZ6D9bOb8TocNGyaaNm2q9zes0WhEeHg4AyA1SpwFTGTEmjVrcOTIEb3H4cOH9dpIpVJ8+eWXKCoqQlRUFIQQ+OKLL+Du7q7XbteuXejfvz/kcjnc3d3h4eGBOXPmoKioCIWFhXptu3btinvuuUf3vFOnTgCAvn376t0Ppz3++++/G9T+9NNP6z0fNWoUAGD37t0mP+93332HiIgIdO3aFXfu3NE9Bg4cWO/Zz9r3u/sm/h49eqBTp07YuXOn3vEWLVqgX79+dX6f2uzatctoHcOGDUPTpk0N6ujatStat26tey6TydChQwe9n/Xu3bvx6KOPIjAwUHfM3d0dI0aMqHedW7duxSOPPKL73TZUkyZNkJCQgP/85z9Qq9UAAI1Gg88++wxDhw6Fv79/va6blJQEAFi+fLnu2OLFi9GlSxc8/PDDDS/cCsz5ne7duxf9+vWDQqHQHXNzc8Pw4cNtWiuRrTAAEhnRqVMndO/eXe/RrVs3g3bt2rXDQw89hNLSUjz99NNQqVR6r2dmZmLAgAEAqr4wf/rpJxw5cgRvvPEGAOD27dt67f38/PSee3p61ni8tLRU73iTJk0MvtiVSiUAoKioyOTn/fPPP/HLL7/Aw8ND7+Hj4wMhBK5evWryXFO073f3zwQAgoKCDOox1s4Y7Rd5bm6u2XU0adIELVu21DsukUigVCoN6jAWjKRSqd7vqqioSPdzrc7YMXNduXLF4hMYxo4di9LSUqxfvx4A8MMPPyA/Px/PPvtsva8ZGBiIESNGYNmyZdBoNPjll1+wb98+TJo0yVJlW5y5v9PqgV7L2DGixoCzgIka4N///je2bNmCHj16YPHixRgxYgR69uype339+vXw8PDAd999B5lMpjv+7bffWqWeO3fuoKioSO8Lr6CgAIDxL0EthUIBLy8vfPrppyZfryvt++Xn5xsEm8uXLxtc09z121QqFbp06YJt27bVOFO4eh137tzBlStX9EKgEAIFBQV44IEHzHrfu6+p/blWZ+yYVCpFWVmZwfG7g2fLli3xxx9/1LmWmoSHh6NHjx5YuXIlXnjhBaxcuRJBQUG6f5TU15QpU/DZZ59h06ZNyMjIgK+vr0HPs7Px9/fHn3/+aXDc2O+UqDFgDyBRPZ08eRKTJ0/GM888g3379uG+++7DiBEjcP36dV0biUSCJk2a6A0L3759G5999pnV6lq3bp3e888//xxA1TCyKUOGDMH58+fh7+9v0PPZvXt3hISE1LkO7XDu2rVr9Y4fOXIEp0+fxqOPPlrna2rNnj0b169fx+TJkyGEMHi9pKQE27ZtAwDd+9xdxzfffIObN2/Wq45HHnkEO3fu1AsMGo0GGzZsMGgbEhKCX375Re/Yrl27UFJSonds0KBB2L17N86cOWPyfbXrK97dc1yTZ599FocPH8b+/fuRnp6OMWPGGNymUNf36datG2JiYvDee+9h3bp1SExMRNOmTc2uyRx399CZW1t99enTB7t27dLr7a6srMRXX31l0fchchTsASQyIjs7G3fu3DE43rZtW7Rs2RI3b97E8OHDERoaiqVLl8LT0xNffvkloqKi8Oyzz+p6+AYPHowFCxZg1KhReP7551FUVIQPP/zQagsle3p6Yv78+SgpKcEDDzyAAwcO4B//+AcGDRqEBx980OR5U6dOxTfffIOHH34Y06ZNw3333YfKykrk5eVh27ZtmDFjhl7Ppjk6duyI559/HosWLYKbmxsGDRqECxcuYPbs2QgODsa0adPq/TmHDRuG2bNn4+2338Z///tfjBs3TrcQ9OHDh7Fs2TKMGDECAwYMQGxsLAYOHIhZs2ahuLgYvXv3xi+//ILk5GRERkZi9OjRdX7/N998E5s3b0a/fv0wZ84ceHt7Y8mSJUZ3Hxk9ejRmz56NOXPmoE+fPsjJycHixYshl8v12r311lvYunUrHn74Ybz++uvo0qUL/vrrL2RkZGD69OkICwtD27Zt4eXlhXXr1qFTp05o1qwZgoKCEBQUZLLWv//975g+fTr+/ve/o6yszKwdYMx5nylTpmDEiBGQSCSYOHGi+T88M3Xp0gXr16/Hhg0bcO+990Imk6FLly71+hmY44033kB6ejoeffRRvPHGG/Dy8sLHH3+s+526ubG/hBoZ+85BIXIsNc0CBiCWL18uhBAiISFBeHt765ax0Prqq68EALFw4ULdsU8//VR07NhRSKVSce+994p58+aJFStWGMws1K5XdjcA4sUXX9Q7lpubKwCIDz74QHdMuw7gL7/8Ivr27Su8vLyEn5+fSEpKMljLz9gadCUlJeLNN9/Urdknl8tFly5dxLRp0/RmuxpjbBawEP9bB7BDhw7Cw8NDKBQKkZCQYHIdwLrau3ev+L//+z+hUqmEh4eHaN68uYiOjhYffPCBKC4u1rW7ffu2mDVrlmjTpo1ufb2kpCST6wDerU+fPgYzUn/66SfdmnFKpVK8/PLLRtcBLCsrE6+88ooIDg4WXl5eok+fPuLEiRNGfwcXL14UY8eOFUqlUrde4fDhw/Vm3H7xxRciLCxMeHh41LoGntaoUaMEANG7d2+jrxurxdT7VP9cUqlUxMXFGb2mMXWZBXzhwgUxYMAA4ePjo1sHsLbaaloH8G7Gfqf79u0TPXv21PudvvfeewKA+Ouvv8z+nETOQCKEkfETInI6iYmJ+Prrrw2GFomsIT09HY8//ji2bNmCxx57zN7lWM2AAQNw4cIFnD171t6lEFkUh4CJiMhsOTk5+P333zFjxgx07doVgwYNsndJFjN9+nRERkYiODgY165dw7p167B9+3asWLHC3qURWRwDIBERmW3ixIn46aefEBUVhdWrV5s9e9sZaDQazJkzBwUFBZBIJAgPD8dnn32GhIQEe5dGZHEcAiYiIiJyMZzWRERERORiGACJiIiIXAwDIBEREZGLYQAkIiIicjGcBdwAlZWVuHz5Mnx8fBrVTDgiIqLGTAiBGzduICgoyGV3eWEAbIDLly8jODjY3mUQERFRPVy8eBGtWrWydxl2wQDYAD4+PgCq/oCaN29u52qIiIjIHMXFxQgODtZ9j7siBsAG0A77Nm/enAGQiIjIybjy7VuuOfBNRERE5MIYAImIiIhcDAMgERERkYvhPYBWJoTAnTt3oNFo7F0KORl3d3c0adLEpe9RISIi63D6ALhnzx488sgjRl/LzMzEAw88YPS1xMRErF69Wu9Yz549cejQIYvVVl5ejvz8fNy6dcti1yTX4u3tDZVKBU9PT3uXQkREjYjTB8CYmBjk5+frHZs9ezZ27NiB7t2713huXFwcVq5cqXtuyS/ZyspK5Obmwt3dHUFBQfD09GRPDplNCIHy8nJcuXIFubm5aN++vcsuVkpERJbn9AHQ09MTSqVS97yiogKbN2/GpEmTag1cUqlU71xLKi8vR2VlJYKDg+Ht7W2V96DGzcvLCx4eHvj9999RXl4OmUxm75KIiKiRaHRdCps3b8bVq1eRmJhYa9s9e/YgICAAHTp0wPjx41FYWFhj+7KyMhQXF+s9asNeG2oI/v0QEZE1NLpvlxUrVmDgwIG1btE2aNAgrFu3Drt27cL8+fNx5MgR9OvXD2VlZSbPmTdvHuRyue7BbeCIiIgsT1MpcPB8ETaduISD54ugqRT2LqnRcdgAmJKSAolEUuPj6NGjeuf88ccf+OGHHzBu3Lharz9ixAgMHjwYERERiI+Px9atW3H27Fls2bLF5DmvvfYa1Gq17nHx4sUGf06ynD179kAikeCvv/6y6HVTUlLQtWtXi16TiIiMy8jOx4Pv7cLflx/ClPUn8Pflh/Dge7uQkZ1f+8lkNoe9B3DSpEkYOXJkjW1CQkL0nq9cuRL+/v54/PHH6/x+KpUKbdq0wblz50y2kUqlkEqldb42uR6JRIKNGzfiiSeesHcpREROIyM7H0lrs3B3f1+BuhRJa7OQlhCFuAiVXWprbBw2ACoUCigUCrPbCyGwcuVKPPPMM/Dw8Kjz+xUVFeHixYtQqRzrD0tTKZCZew2FN0oR4CNDj1A/uLtxNjERETUumkqB1PQcg/AHAAKABEBqeg5iw5X8HrQAhx0Crqtdu3YhNzfX5PBvWFgYNm7cCAAoKSnBzJkzcfDgQVy4cAF79uxBfHw8FAoFnnzySVuWXSN7dIMLIfD+++/j3nvvhZeXF+6//358/fXXutf69++PuLg4CFH1n+hff/2F1q1b44033gAAaDQajBs3DqGhofDy8kLHjh3x0Ucf6b1HYmIinnjiCbzzzjsIDAyEr68vUlNTcefOHbz88svw8/NDq1at8Omnn+rOuXDhAiQSCdavX4+YmBjIZDJ07twZe/bsqfHzHDhwAA8//DC8vLwQHByMyZMn4+bNmzWe8+677yIwMBA+Pj4YN24cSktL9V4/cuQIYmNjoVAoIJfL0adPH2RlZele1/ZMP/nkk5BIJLrn58+fx9ChQxEYGIhmzZrhgQcewI4dO2qshYjIVWTmXkO+utTk6wJAvroUmbnXbFdUI9ZoAuCKFSsQExODTp06GX39zJkzUKvVAKp2WDh58iSGDh2KDh06YMyYMejQoQMOHjwIHx8fW5ZtkrYb/O7/GLTd4NYKgW+++SZWrlyJtLQ0nDp1CtOmTUNCQgL27t0LiUSC1atXIzMzE//6178AABMmTEBgYCBSUlIAVK1/2KpVK3z55ZfIycnBnDlz8Prrr+PLL7/Ue59du3bh8uXL+PHHH7FgwQKkpKRgyJAhaNGiBQ4fPowJEyZgwoQJBvdZvvzyy5gxYwaOHz+OmJgYPP744ygqKjL6WU6ePImBAwfiqaeewi+//IINGzZg//79mDRpksnP/+WXXyI5ORlz587F0aNHoVKpsHTpUr02N27cwJgxY7Bv3z4cOnQI7du3x2OPPYYbN24AqAqIQNUtCfn5+brnJSUleOyxx7Bjxw4cP34cAwcORHx8PPLy8sz87RARNV6FN0yHv/q0o5pJhLYrh+qsuLgYcrkcarUazZs313uttLQUubm5CA0NrfP6bZpKgQff22XyX0ISAEq5DPtn9bNoN/jNmzehUCiwa9cuREdH644/99xzuHXrFj7//HMAwFdffYXRo0dj+vTp+Oijj3D8+HF06NDB5HVffPFF/Pnnn7qexMTEROzZswe//fabbpmTsLAwBAQE4McffwRQ1ZMol8vx73//GyNHjsSFCxcQGhqKd999F7NmzQIA3LlzB6GhoXjppZfwyiuv6HaFuX79Onx9ffHMM8/Ay8sLy5Yt09Wyf/9+9OnTBzdv3jT6e4mJicH999+PtLQ03bFevXqhtLQUJ06cMPr5NBoNWrRogc8//xxDhgwBYP49gJ07d0ZSUpLJUNqQvyMiImdy8HwR/r689t24vhjfC9Ft/Rv0XjV9f7uKRtMD2JjYqxs8JycHpaWliI2NRbNmzXSPNWvW4Pz587p2w4YNw1NPPYV58+Zh/vz5BuHv448/Rvfu3dGyZUs0a9YMy5cvN+jl6ty5s94ad4GBgejSpYvuubu7O/z9/Q3WZqweTJs0aYLu3bvj9OnTRj/PsWPHsGrVKr3PMnDgQN0uLcacPn1a7z3ufk8AKCwsxIQJE9ChQwfdkkAlJSW19uTdvHkTr7zyCsLDw+Hr64tmzZrhv//9L3sAiYgA9Aj1g0oug6luDQkAlbzqXnhqOIedBOLK7NUNXllZCQDYsmUL7rnnHr3Xqs9+vnXrFo4dOwZ3d3eDWdNffvklpk2bhvnz5yM6Oho+Pj744IMPcPjwYb12d0/UkUgkRo9pa6qJqR1fKisr8cILL2Dy5MkGr7Vu3brW65qSmJiIK1eu4J///CfatGkDqVSK6OholJeX13jeyy+/jB9++AEffvgh2rVrBy8vL/zf//1frecREbkCdzcJkuPDkbQ2CxJAbzKI9v/lk+PDOQHEQhgAHVCAj3lDfea2M1d4eDikUiny8vLQp08fk+1mzJgBNzc3bN26FY899hgGDx6Mfv36AQD27duHmJgYTJw4Ude+eu9hQx06dAgPP/wwgKoh4GPHjpkcPo2KisKpU6fQrl07s6/fqVMnHDp0CM8884zee1a3b98+LF26FI899hgA4OLFi7h69apeGw8PD2g0GoPzEhMTdRONSkpKcOHCBbNrIyJq7OIiVEhLiEJqeo7eSJhSLkNyfDiXgLEgBkAHpO0GL1CXGp0Or70H0NLd4D4+Ppg5cyamTZuGyspKPPjggyguLsaBAwfQrFkzjBkzBlu2bMGnn36KgwcPIioqCq+++irGjBmDX375BS1atEC7du2wZs0a/PDDDwgNDcVnn32GI0eOIDQ01CI1LlmyBO3bt0enTp2wcOFCXL9+HWPHjjXadtasWejVqxdefPFFjB8/Hk2bNsXp06exfft2LFq0yOg5U6ZMwZgxY9C9e3c8+OCDWLduHU6dOoV7771X16Zdu3b47LPP0L17dxQXF+Pll1+Gl5eX3nVCQkKwc+dO9O7dG1KpVPez+c9//oP4+HhIJBLMnj3brB5OIiJXEhehQmy4kkugWRnvAXRA2m5wAAb3Qli7G/ztt9/GnDlzMG/ePHTq1AkDBw5Eeno6QkNDceXKFYwbNw4pKSmIioqqqiM5GUFBQZgwYQKAqlnBTz31FEaMGIGePXuiqKhIrzewod5991289957uP/++7Fv3z5s2rTJ5HqR9913H/bu3Ytz587hoYceQmRkJGbPnl3jWo8jRozAnDlzMGvWLHTr1g2///47kpKS9Np8+umnuH79OiIjIzF69GhMnjwZAQEBem3mz5+P7du3Izg4GJGRkQCAhQsXokWLFoiJiUF8fDwGDhyo+zkSEdH/uLtJEN3WH0O73oPotv4Mf1bAWcANYK1ZwFoZ2fkG3eAqF+0G184CPn78uEtty8ZZwERElsdZwBwCdmjsBiciIiJrYAB0cNpucCIiIiJLYQAkpxASEgLerUBERGQZDIBERERkM5pKwVubHAADIBEREdkEJzc6Di4DY2UctqSG4N8PETUWGdn5SFqbZbDVaYG6FElrs5CRnW+nylwTA6CVaLc1u3Xrlp0rIWem/fu5e5s8IiJnoqkUSE3PMbq5gfZYanoONJX8R6+tcAjYStzd3eHr64vCwkIAgLe3t8k9a4nuJoTArVu3UFhYCF9fX7i7u9u7JCKiesvMvWbQ81edAJCvLkVm7jWufGEjDIBWpFQqAUAXAonqytfXV/d3RETkrApvmA5/9WlHDccAaEUSiQQqlQoBAQGoqKiwdznkZDw8PNjzR0SNQoCPeTsZmduOGo4B0Abc3d35RU5ERC6rR6gfVHIZCtSlRu8DlABQyquWhCHb4CQQIiIisip3NwmS48MBVIW96rTPk+PDuR6gDTEAEhERkdXFRaiQlhAFpVx/mFcplyEtIYrrANoYh4CJiIjIJuIiVIgNV3InEAfAAEhEREQ24+4m4VIvDoBDwEREREQuhgGQiIiIyMUwABIRERG5GN4DSERERBahqRSc4OEkGACJiIiowTKy85GanqO3569KLkNyfDiXeHFAHAImIiKiBsnIzkfS2iy98AcABepSJK3NQkZ2vp0qI1MYAImIiKjeNJUCqek5Rrd40x5LTc+BptJYC7IXBkAiIiKqt8zcawY9f9UJAPnqUmTmXrNdUVQrBkAiIiKqt8IbpsNffdqRbTAAEhERUb0F+Mhqb1SHdmQbDIBERERUbz1C/aCSy2BqsRcJqmYD9wj1s2VZVAsGQCIiIqo3dzcJkuPDAcAgBGqfJ8eHcz1AB8MASERERA0SF6FCWkIUlHL9YV6lXIa0hCiuA+iAuBA0ERERNVhchAqx4UruBOIkGACJiIjIItzdJIhu62/vMsgMHAImIiIicjGNIgCePXsWQ4cOhUKhQPPmzdG7d2/s3r27xnOEEEhJSUFQUBC8vLzQt29fnDp1ykYVExEREdlPowiAgwcPxp07d7Br1y4cO3YMXbt2xZAhQ1BQUGDynPfffx8LFizA4sWLceTIESiVSsTGxuLGjRs2rJyIiIjI9iRCCKfenO/q1ato2bIlfvzxRzz00EMAgBs3bqB58+bYsWMHHn30UYNzhBAICgrC1KlTMWvWLABAWVkZAgMD8d577+GFF14w672Li4shl8uhVqvRvHlzy30oIiIiB6KpFI1qcge/vxvBJBB/f3906tQJa9asQVRUFKRSKZYtW4bAwEB069bN6Dm5ubkoKCjAgAEDdMekUin69OmDAwcOmAyAZWVlKCsr0z0vLi627IchIiJyMBnZ+UhNz9Hb71cllyE5PpzLuzgxpx8Clkgk2L59O44fPw4fHx/IZDIsXLgQGRkZ8PX1NXqOdmg4MDBQ73hgYGCNw8bz5s2DXC7XPYKDgy32OYiIiBxNRnY+ktZm6YU/AChQlyJpbRYysvPtVBk1lMMGwJSUFEgkkhofR48ehRACEydOREBAAPbt24fMzEwMHToUQ4YMQX5+zX+YEol+97UQwuBYda+99hrUarXucfHiRYt8ViIiIkejqRRITc+BsfvEtMdS03OgqXTqO8lclsMOAU+aNAkjR46ssU1ISAh27dqF7777DtevX9eN4y9duhTbt2/H6tWr8eqrrxqcp1QqAVT1BKpU/+u+LiwsNOgVrE4qlUIqldbn4xARETmVzNxrBj1/1QkA+epSZOZe49p/TshhA6BCoYBCoai13a1btwAAbm76nZlubm6orKw0ek5oaCiUSiW2b9+OyMhIAEB5eTn27t2L9957r4GVExEROb/CG6bDX33akWNx2CFgc0VHR6NFixYYM2YMfv75Z5w9exYvv/wycnNzMXjwYF27sLAwbNy4EUDV0O/UqVPxzjvvYOPGjcjOzkZiYiK8vb0xatQoe30UIiIihxHgI6u9UR3akWNx2B5AcykUCmRkZOCNN95Av379UFFRgc6dO2PTpk24//77de3OnDkDtVqte/7KK6/g9u3bmDhxIq5fv46ePXti27Zt8PHxscfHICIicig9Qv2gkstQoC41eh+gBIBSXrUkDDkfp18H0J64jhARETVm2lnAAPRCoHa6ZFpClFMuBcPv70YwBExERETWERehQlpCFJRy/WFepVzmtOGPqjj9EDARERFZT1yECrHhyka1EwgxABIREVEt3N0kXOqlkeEQMBEREZGLYQAkIiIicjEMgEREREQuhvcAEhERuTBNpeAEDxfEAEhEROSiMrLzkZqeo7fnr0ouQ3J8OJd4aeQ4BExEROSCtIs8Vw9/AFCgLkXS2ixkZOfbqTKyBQZAIiIiF6OpFEhNzzG6xZv2WGp6DjSV3CyssWIAJCIicjGZudcMev6qEwDy1aXIzL1mu6LIphgAiYiIXEzhDdPhrz7tyPkwABIREbmYAB9Z7Y3q0I6cDwMgERGRi+kR6geVXAZTi71IUDUbuEeony3LIhtiACQiInIx7m4SJMeHA4BBCNQ+T44P53qAjRgDIBERkQuKi1AhLSEKSrn+MK9SLkNaQhTXAWzkuBA0ERGRi4qLUCE2XMmdQFwQAyAREZELc3eTILqtv73LIBvjEDARERGRi2EAJCIiInIxHAImIiJqhDSVgvf2kUkMgERERI1MRnY+UtNz9LZ7U8llSI4P5+xeAsAhYCIiokYlIzsfSWuzDPb6LVCXImltFjKy8+1UGTkSBkAiIqJGQlMpkJqeA2HkNe2x1PQcaCqNtSBXwgBIRETUSGTmXjPo+atOAMhXlyIz95rtiiKHxABIRETUSBTeMB3+6tOOGi8GQCIiokYiwEdWe6M6tKPGiwGQiIiokegR6geVXAZTi71IUDUbuEeony3LIgfEAEhERNRIuLtJkBwfDgAGIVD7PDk+nOsBEgMgERFRYxIXoUJaQhSUcv1hXqVchrSEKK4DSAC4EDQREVGjExehQmy4kjuBkEkMgERERI2Qu5sE0W397V0GOSgOARMRERG5GAZAIiIiIhfDIWAiIiIno6kUvL+PGoQBkIiIyIlkZOcjNT1Hb8s3lVyG5PhwzvAlszWKIeCzZ89i6NChUCgUaN68OXr37o3du3fXeE5iYiIkEoneo1evXjaqmIiIqO4ysvORtDbLYL/fAnUpktZmISM7306VkbNpFAFw8ODBuHPnDnbt2oVjx46ha9euGDJkCAoKCmo8Ly4uDvn5+brH999/b6OKiYiI6kZTKZCangNh5DXtsdT0HGgqjbUg0uf0AfDq1av49ddf8eqrr+K+++5D+/bt8e677+LWrVs4depUjedKpVIolUrdw8+PW+MQEZFjysy9ZtDzV50AkK8uRWbuNdsVRU7L6QOgv78/OnXqhDVr1uDmzZu4c+cOli1bhsDAQHTr1q3Gc/fs2YOAgAB06NAB48ePR2FhoY2qJiIiqpvCG6bDX33akWtz+kkgEokE27dvx9ChQ+Hj4wM3NzcEBgYiIyMDvr6+Js8bNGgQhg0bhjZt2iA3NxezZ89Gv379cOzYMUilUqPnlJWVoaysTPe8uLjY0h+HiIjIqAAfWe2N6tCOXJvD9gCmpKQYTNK4+3H06FEIITBx4kQEBARg3759yMzMxNChQzFkyBDk55u+GXbEiBEYPHgwIiIiEB8fj61bt+Ls2bPYsmWLyXPmzZsHuVyuewQHB1vjoxMRERnoEeoHlVwGU4u9SFA1G7hHKG9notpJhBAOebfo1atXcfXq1RrbhISE4KeffsKAAQNw/fp1NG/eXPda+/btMW7cOLz66qtmv2f79u3x3HPPYdasWUZfN9YDGBwcDLVarffeRERE1qCdBQxAbzKINhSmJURxKRgzFBcXQy6Xu/T3t8MOASsUCigUilrb3bp1CwDg5qbfmenm5obKykqz36+oqAgXL16ESmX6PxypVGpyeJiIiMja4iJUSEuIMlgHUMl1AKmOHDYAmis6OhotWrTAmDFjMGfOHHh5eWH58uXIzc3F4MGDde3CwsIwb948PPnkkygpKUFKSgr+9re/QaVS4cKFC3j99dehUCjw5JNP2vHTEBERVTG120dchAqx4UruBEIN4vQBUKFQICMjA2+88Qb69euHiooKdO7cGZs2bcL999+va3fmzBmo1WoAgLu7O06ePIk1a9bgr7/+gkqlwiOPPIINGzbAx8fHXh+FiIgIQO27fbi7SRDd1t+OFZKzc9h7AJ0B7yEgIiJL097nd/eXM+/zsxx+fzvwLGAiIiJXw90+yFYYAImIiBzEofNF3O2DbIIBkIiIyAFkZOfjxc+zzGrL3T6ooZx+EggREZGzM3Xfnync7YMaigGQiIjIjmq67+9uElSt+cfdPqihOARMRERkR5m512q87+9uyfHhXPOPGowBkIiIyI7MvZ/P18uDS8CQxXAImIiIyI7MvZ9vydNR6N2u9i1SiczBHkAiIiI76hHqB5VcBlODuhJU7QLS617u/EGWwwBIRERkR+5uEiTHhwOAQQjUPud9f2RpDIBERER2FhehQlpCFJRy/eFgpVzG+/7IKngPIBERkQOIi1AhNlyJzNxrKLxRigCfquVe2PNH1sAASERE5CDc3SSIbst7/cj6OARMRERE5GLYA0hERGQDmkrB4V1yGAyAREREVpaRnY/U9By9HT9UchmS48M5wYPsgkPAREREVpSRnY+ktVkG270VqEuRtDYLGdn5dqqMXBkDIBERkZVoKgVS03MgjLymPZaangNNpbEWRNbDAEhERGQlmbnXDHr+qhMA8tWlyMy9ZruiiMAASEREZDWFN0yHv/q0I7IUBkAiIiIrCfCR1d6oDu2ILIUBkIiIyEp6hPpBJZcZ7PGrJUHVbOAeoX62LIuIAZCIiMha3N0kSI4PBwCDEKh9nhwfzvUAyeYYAImIiKwoLkKFtIQoKOX6w7xKuQxpCVFcB5DsggtBExERWVlchAqx4UruBEIOgwGQiIjIBtzdJIhu62/vMogAMAASERE1CPf4JWfEAEhERFRP3OOXnBUngRAREdUD9/glZ8YASEREVEfc45ecHQMgERFRHXGPX3J2DIBERER1xD1+ydkxABIREdUR9/glZ8cASEREVEfc45ecHQMgERFRHXGPX3J2DIBERET1wD1+yZlxIWgiIqJ64h6/5KwaRQ9gVlYWYmNj4evrC39/fzz//PMoKSmp8RwhBFJSUhAUFAQvLy/07dsXp06dslHFRETUWGj3+B3a9R5Et/Vn+COn4PQB8PLly+jfvz/atWuHw4cPIyMjA6dOnUJiYmKN573//vtYsGABFi9ejCNHjkCpVCI2NhY3btywTeFEROTwNJUCB88XYdOJSzh4vogLO1Oj4fRDwN999x08PDywZMkSuLlV5dklS5YgMjISv/76K9q1a2dwjhAC//znP/HGG2/gqaeeAgCsXr0agYGB+Pzzz/HCCy/Y9DMQEZHj4T6/1Jg5fQ9gWVkZPD09deEPALy8vAAA+/fvN3pObm4uCgoKMGDAAN0xqVSKPn364MCBA9YtmIiIHB73+aXGzukDYL9+/VBQUIAPPvgA5eXluH79Ol5//XUAQH6+8f9ACwoKAACBgYF6xwMDA3WvGVNWVobi4mK9BxERNS7c55dcgcMGwJSUFEgkkhofR48eRefOnbF69WrMnz8f3t7eUCqVuPfeexEYGAh3d/ca30Mi0b9RVwhhcKy6efPmQS6X6x7BwcEW+axEROQ4uM8vuQKHvQdw0qRJGDlyZI1tQkJCAACjRo3CqFGj8Oeff6Jp06aQSCRYsGABQkNDjZ6nVCoBVPUEqlT/u4+jsLDQoFewutdeew3Tp0/XPS8uLmYIJCJqZLjPL7kCswJgRUUFCgoKcOvWLbRs2RJ+ftbf2kahUEChUNTpHG14+/TTTyGTyRAbG2u0XWhoKJRKJbZv347IyEgAQHl5Ofbu3Yv33nvP5PWlUimkUmmdaiIiIufCfX7JFZgcAi4pKcGyZcvQt29fyOVyhISEIDw8HC1btkSbNm0wfvx4HDlyxJa1mrR48WJkZWXh7NmzWLJkCSZNmoR58+bB19dX1yYsLAwbN24EUDX0O3XqVLzzzjvYuHEjsrOzkZiYCG9vb4waNcpOn4KIiBwB9/klV2C0B3DhwoWYO3cuQkJC8Pjjj+PVV1/FPffcAy8vL1y7dg3Z2dnYt28fYmNj0atXLyxatAjt27e3de06mZmZSE5ORklJCcLCwrBs2TKMHj1ar82ZM2egVqt1z1955RXcvn0bEydOxPXr19GzZ09s27YNPj4+ti6fiIgciHaf36S1WZAAepNBuM8vNRYSIYTBNKZhw4Zhzpw56NKlS40nl5WVYcWKFfD09MRzzz1ntSIdVXFxMeRyOdRqNZo3b27vcoiIyIK4DmDjxe9vEwGQzMM/ICKixk1TKbjPbyPE728HngVMRERkTeaEO+0+v0SNjUEA3LlzJ3r27IlmzZrZox4iIiKr4/AuuTqDWcB5eXl46KGHatwRg4iIyFlxmzciIz2Azz77LLy9vREbG4uTJ0/aoyYiIiKrqG2bNwmqtnmLDVfyXj9q1IyuAzhixAh88803tq6FiIjIqrjNG1EVkwtBd+jQwZZ1EBERWR23eSOqYjIAahUXF9uiDiIiIqvjNm9EVWoNgC1atOBwMBERNQrc5o2oSq0BUAiBtLQ09OzZE7169cKkSZNw+PBhW9RGRERkUdpt3gAYhEBu80aupNYACAA///wzevTogb59++LMmTPo06cPpk2bZu3aiIiILC4uQoW0hCgo5frDvEq5DGkJUVwHkFyCWTuBfP7554iNjdU9P3nyJJ544gm0atUKM2bMsFpxRERE1hAXoUJsuJLbvJHLqrUH0N/fH8HBwXrHunTpgn/961/4+OOPrVYYERGRNWm3eRva9R5Et/Vn+COXUmsAvP/++7FixQqD4+3atcPFixetUhQREVF9aSoFDp4vwqYTl3DwfBE0lcaWfSZybbUOAf/jH//AI488gkuXLmHixIm47777cPv2bbzzzjsIDQ21RY1ERERm4R6/ROaptQewV69eOHToEC5fvoy+ffuiRYsWCAoKwtdff4358+fbokYiIqJacY9fIvNJhBBm940XFhbi2LFjqKysRM+ePaFQKKxZm8MrLi6GXC6HWq1G8+bN7V0OEZHL0lQKPPjeLpPbvElQNct3/6x+vNeP+P0NM2cBawUEBGDQoEHWqoWIiKhe6rLHb3Rbf9sVRuSgzFoHkIiIyJFxj1+iumEAJCIip8c9fonqhgGQiIicHvf4JaobBkAiInJ63OOXqG6MBsC8vLw6XeTSpUsWKYaIiKi+uMcvkfmMLgMTGBiIxx9/HOPHj0ePHj2MnqhWq/Hll1/io48+wgsvvICXXnrJ6sU6Gk4jJyJyPJpKwT1+qUb8/jaxDMzp06fxzjvvIC4uDh4eHujevTuCgoIgk8lw/fp15OTk4NSpU+jevTs++OADLg1DREQ2YU640+7xS0Sm1bgQdGlpKb7//nvs27cPFy5cwO3bt6FQKBAZGYmBAwciIiLClrU6HP4LgojIdrjNG1kKv7/ruBMI6eMfEBGRbWi3ebv7C0vb98d7/Kgu+P3NWcBEROTgNJUCqek5BuEPgO5YanoONJXszyAyFwMgERE5tLps80ZE5mEAJCIih8Zt3ogsjwGQiIgcGrd5I7I8BkAiInJImkqBg+eLUFBcCr+mHibbcZs3orozug5gddOnTzd6XCKRQCaToV27dhg6dCj8/PgfHhERWYaxJV+M4TZvRPVT6zIwjzzyCLKysqDRaNCxY0cIIXDu3Dm4u7sjLCwMZ86cgUQiwf79+xEeHm6ruh0Cp5ETEVmeqSVfjOE6gFQf/P42owdQ27u3cuVK3Q+puLgY48aNw4MPPojx48dj1KhRmDZtGn744QerF0xERI1XTUu+AFU9fn5NPfHm4E5Qyr24zRtRPdXaA3jPPfdg+/btBr17p06dwoABA3Dp0iVkZWVhwIABuHr1qlWLdTT8FwQRkWUdPF+Evy8/VGu7L8b34nZvVG/8/jZjEoharUZhYaHB8StXrqC4uBgA4Ovri/LycstXR0RELoVLvhDZRq0BcOjQoRg7diw2btyIP/74A5cuXcLGjRsxbtw4PPHEEwCAzMxMdOjQwdq1mpSVlYXY2Fj4+vrC398fzz//PEpKSmo8JzExERKJRO/Rq1cvG1VMRETVaWf8nvuz5v/v1uKSL0QNU+s9gMuWLcO0adMwcuRI3Llzp+qkJk0wZswYLFy4EAAQFhaGf//739at1ITLly+jf//+GDFiBBYvXozi4mJMnToViYmJ+Prrr2s8Ny4uDitXrtQ99/T0tHa5RER0F3Nn/AJV9wAqueQLUYPVGgCbNWuG5cuXY+HChfjtt98ghEDbtm3RrFkzXZuuXbtas8Yafffdd/Dw8MCSJUvg5lbVoblkyRJERkbi119/Rbt27UyeK5VKoVQqbVUqERHdpS4zfrnkC5HlmL0QdLNmzXDffffh/vvv1wt/9lZWVgZPT09d+AMALy8vAMD+/ftrPHfPnj0ICAhAhw4dMH78eKP3OhIRkXXUNuP3bkq5DGkJUVzyhcgCau0BBIC//voLK1aswOnTpyGRSNCpUyeMGzcOcrnc2vXVql+/fpg+fTo++OADTJkyBTdv3sTrr78OAMjPzzd53qBBgzBs2DC0adMGubm5mD17Nvr164djx45BKpUaPaesrAxlZWW659pJMEREVHeZudfMGvad9Ehb9G7Xkku+EFlQrT2AR48eRdu2bbFw4UJcu3YNV69excKFC9G2bVtkZWVZrbCUlBSDSRp3P44ePYrOnTtj9erVmD9/Pry9vaFUKnHvvfciMDAQ7u7uJq8/YsQIDB48GBEREYiPj8fWrVtx9uxZbNmyxeQ58+bNg1wu1z2Cg4Ot8dGJiFyCuTN52wf6ILqtP8MfkQXVug7gQw89hHbt2mH58uVo0qSqw/DOnTt47rnn8Ntvv+HHH3+0SmFXr16tdV3BkJAQyGT/mwn2559/omnTppBIJGjevDnWr1+PYcOGmf2e7du3x3PPPYdZs2YZfd1YD2BwcLBLryNERFRfXPOP7IXrAJoxBHz06FG98AdUzQJ+5ZVX0L17d6sVplAooFAo6nROYGAgAODTTz+FTCZDbGys2ecWFRXh4sWLUKlM31silUpNDg8TEVHd9Aj1g0ouQ4G61Oh9gJzxS2Q9tQ4BN2/eHHl5eQbHL168CB8fH6sUVVeLFy9GVlYWzp49iyVLlmDSpEmYN28efH19dW3CwsKwceNGAEBJSQlmzpyJgwcP4sKFC9izZw/i4+OhUCjw5JNP2ulTEBG5Fnc3CZLjq3aZuntwlzN+iayr1gA4YsQIjBs3Dhs2bMDFixfxxx9/YP369Xjuuefw97//3RY11iozMxOxsbHo0qULPvnkEyxbtgyTJ0/Wa3PmzBmo1WoAgLu7O06ePImhQ4eiQ4cOGDNmDDp06ICDBw86TKglInIFcREqpCVEQSnXX9iZM36JrKvWewDLy8vx8ssv4+OPP8adO3cghICnpyeSkpLw7rvvuvSQKO8hICKyDE2lQGbuNRTeKEWAj4wzfsmq+P1tRgDUunXrFs6fPw8hBNq1awdvb29r1+bw+AdERETkfPj9bWISyPTp082+wIIFCyxWDBERERFZn9EAePz4cbNOlkjYPU9ERDXj8C6R4zEaAHfv3m3rOoiIqBHKyM5HanqO3o4fKrkMyfHhnOBBZEdm7wVMRERUFxnZ+Uham2Ww3VuBuhRJa7OQkW16u04isi4GQCIisjhNpUBqeo7RBZ61x1LTc6CpNGseIhFZGAMgERFZXGbuNYOev+oEgHx1KTJzr9muKCLSYQAkIiKLK7xhOvzVpx0RWRYDIBERWVyAj6z2RnVoR0SWxQBIREQW1yPUDyq5zGCPXy0JqmYD9wj1s2VZRPT/MQASEZHFubtJkBwfDgAGIVD7PDk+nOsBEtkJAyAREVlFXIQKaQlRUMr1h3mVchnSEqK4DiCRHRldCJqIiMgS4iJUiA1XcicQIgfDAEhERFbl7iZBdFt/e5dBRNVwCJiIiIjIxTAAEhEREbkYBkAiIiIiF8N7AImIqF40lYKTO4icFAMgERHVWUZ2PlLTc/T2+1XJZUiOD+fyLkROgEPARERUJxnZ+Uham6UX/gCgQF2KpLVZyMjOt1NlRGQuBkAiIjKbplIgNT0Hwshr2mOp6TnQVBprQUSOggGQiIjMlpl7zaDnrzoBIF9diszca7YriojqjAGQiIjMVnjDdPirTzsisg8GQCIiMluAj6z2RnVoR0T2wQBIRERm6xHqB5VcBlOLvUhQNRu4R6ifLcsiojpiACQiIrO5u0mQHB8OAAYhUPs8OT6c6wESOTgGQCIiqpO4CBXSEqKglOsP8yrlMqQlRHEdQCInwIWgiYiozuIiVIgNV3InECInxQBIRET14u4mQXRbf3uXQUT1wCFgIiIiIhfDAEhERETkYhgAiYiIiFwM7wEkIiIDmkrBCR5EjRgDIBER6cnIzkdqeo7enr8quQzJ8eFc4oWokeAQMBER6WRk5yNpbZZe+AOAAnUpktZmISM7306VEZElMQASERGAqmHf1PQcCCOvaY+lpudAU2msBRE5EwZAIiICAGTmXjPo+atOAMhXlyIz95rtiiIiq3D4ADh37lzExMTA29sbvr6+Rtvk5eUhPj4eTZs2hUKhwOTJk1FeXl7jdcvKyvDSSy9BoVCgadOmePzxx/HHH39Y4RMQETmHwhumw1992hGR43L4AFheXo5hw4YhKSnJ6OsajQaDBw/GzZs3sX//fqxfvx7ffPMNZsyYUeN1p06dio0bN2L9+vXYv38/SkpKMGTIEGg0Gmt8DCIihxfgI6u9UR3aEZHjkgghnOJmjlWrVmHq1Kn466+/9I5v3boVQ4YMwcWLFxEUFAQAWL9+PRITE1FYWIjmzZsbXEutVqNly5b47LPPMGLECADA5cuXERwcjO+//x4DBw40q6bi4mLI5XKo1Wqj70NE5Ew0lQIPvrcLBepSo/cBSgAo5TLsn9WPS8KQU+P3txP0ANbm4MGDiIiI0IU/ABg4cCDKyspw7Ngxo+ccO3YMFRUVGDBggO5YUFAQIiIicODAAZPvVVZWhuLiYr0HEVFj4e4mQXJ8OICqsFed9nlyfDjDH1Ej4PQBsKCgAIGBgXrHWrRoAU9PTxQUFJg8x9PTEy1atNA7HhgYaPIcAJg3bx7kcrnuERwc3PAPQETkQOIiVEhLiIJSrj/Mq5TLkJYQxXUAiRoJuywEnZKSgtTU1BrbHDlyBN27dzfrehKJ4b9GhRBGj9ektnNee+01TJ8+Xfe8uLiYIZCIGp24CBViw5XcCYSoEbNLAJw0aRJGjhxZY5uQkBCzrqVUKnH48GG9Y9evX0dFRYVBz2D1c8rLy3H9+nW9XsDCwkLExMSYfC+pVAqpVGpWXUREzszdTYLotv72LoOIrMQuAVChUEChUFjkWtHR0Zg7dy7y8/OhUlUNTWzbtg1SqRTdunUzek63bt3g4eGB7du3Y/jw4QCA/Px8ZGdn4/3337dIXURERESOyuHvAczLy8OJEyeQl5cHjUaDEydO4MSJEygpKQEADBgwAOHh4Rg9ejSOHz+OnTt3YubMmRg/frxuZs+lS5cQFhaGzMxMAIBcLse4ceMwY8YM7Ny5E8ePH0dCQgK6dOmC/v372+2zEhEREdmCXXoA62LOnDlYvXq17nlkZCQAYPfu3ejbty/c3d2xZcsWTJw4Eb1794aXlxdGjRqFDz/8UHdORUUFzpw5g1u3bumOLVy4EE2aNMHw4cNx+/ZtPProo1i1ahXc3d1t9+GIiIiI7MBp1gF0RFxHiIicjaZScHIHuTx+fztBDyAREVlGRnY+UtNz9Pb7VcllSI4P5/IuRC7G4e8BJCKihsvIzkfS2iy98AcABepSJK3NQkZ2vp0qIyJ7YAAkImrkNJUCqek5Rrd30x5LTc+BppJ3BBG5CgZAIqJGLjP3mkHPX3UCQL66FJm512xXFBHZFQMgEVEjV3jDdPirTzsicn4MgEREjVyAj6z2RnVoR0TOjwGQiKiR6xHqB5VcBlOLvUhQNRu4R6ifLcsiIjtiACQiauTc3SRIjg8HAIMQqH2eHB/O9QCJXAgDIBGRC4iLUCEtIQpKuf4wr1IuQ1pCFNcBJHIxXAiaiMhFxEWoEBuu5E4gRMQASETkStzdJIhu62/vMojIzjgETERERORiGACJiIiIXAyHgImIGglNpeD9fURkFgZAIqJGICM7H6npOXpbvqnkMiTHh3OGLxEZ4BAwEZGTy8jOR9LaLIP9fgvUpUham4WM7Hw7VUZEjooBkIjIiWkqBVLTcyCMvKY9lpqeA02lsRZE5KoYAImInFhm7jWDnr/qBIB8dSkyc6/ZrigicngMgERETqzwhunwV592ROQaGACJiJxYgI+s9kZ1aEdEroEBkIjIifUI9YNKLoOpxV4kqJoN3CPUz5ZlEZGDYwAkInJi7m4SJMeHA4BBCNQ+T44P53qARKSHAZCIyMnFRaiQlhAFpVx/mFcplyEtIYrrABKRAS4ETUTUCMRFqBAbruROIERkFgZAIqJGwt1Ngui2/vYug4icAIeAiYiIiFwMewCJiJyAplJweJeILIYBkIjIwWVk5yM1PUdvxw+VXIbk+HBO8CCieuEQMBGRA8vIzkfS2iyD7d4K1KVIWpuFjOx8O1VGRM6MAZCIyEFpKgVS03MgjLymPZaangNNpbEWRESmMQASETmozNxrBj1/1QkA+epSZOZes11RRNQoMAASETkgTaXAT79eMatt4Q3TIZGIyBhOAiEicjDGJn3UJMBHVnsjIqJqGACJiByIdtKHOXf1SVC13VuPUD9rl0VEjQyHgImIHERNkz7upl0BMDk+nOsBElGdsQeQiMhB1Dbpozol1wEkogZgACQichDmTuaY9Eg7TIvtwJ4/Iqo3hx8Cnjt3LmJiYuDt7Q1fX1+jbfLy8hAfH4+mTZtCoVBg8uTJKC8vr/G6ffv2hUQi0XuMHDnSCp+AiMg85k7m6N1OwfBHRA3i8D2A5eXlGDZsGKKjo7FixQqD1zUaDQYPHoyWLVti//79KCoqwpgxYyCEwKJFi2q89vjx4/HWW2/pnnt5eVm8fiKimlTf41fRTAplcyn+LC4zeh8gJ30QkaU4fABMTU0FAKxatcro69u2bUNOTg4uXryIoKAgAMD8+fORmJiIuXPnonnz5iav7e3tDaVSafGaiYjMYWy5F19vDwhUhb3qIZCTPojIkhx+CLg2Bw8eREREhC78AcDAgQNRVlaGY8eO1XjuunXroFAo0LlzZ8ycORM3btyosX1ZWRmKi4v1HkRE9WFqj1/1rQoAgNzbQ++4Ui5DWkIUJ30QkUU4fA9gbQoKChAYGKh3rEWLFvD09ERBQYHJ855++mmEhoZCqVQiOzsbr732Gn7++Wds377d5Dnz5s3T9UgSEdVXbXv8SgB4ebhjybgoXL1ZhgCfqmFf9vwRkaXYpQcwJSXFYALG3Y+jR4+afT2JxPD/FIUQRo9rjR8/Hv3790dERARGjhyJr7/+Gjt27EBWVpbJc1577TWo1Wrd4+LFi2bXSESkZe4ev25uEgzteg+i2/oz/BGRRdmlB3DSpEm1zrgNCQkx61pKpRKHDx/WO3b9+nVUVFQY9AzWJCoqCh4eHjh37hyioqKMtpFKpZBKpWZfk4jIGHOXe+Eev0RkLXYJgAqFAgqFwiLXio6Oxty5c5Gfnw+VquremG3btkEqlaJbt25mX+fUqVOoqKjQXYOIyFrMXe6Fe/wSkbU4/CSQvLw8nDhxAnl5edBoNDhx4gROnDiBkpISAMCAAQMQHh6O0aNH4/jx49i5cydmzpyJ8ePH62YAX7p0CWFhYcjMzAQAnD9/Hm+99RaOHj2KCxcu4Pvvv8ewYcMQGRmJ3r172+2zEpFr6BHqB5VcBlODuhIAKi73QkRW5PABcM6cOYiMjERycjJKSkoQGRmJyMhI3T2C7u7u2LJlC2QyGXr37o3hw4fjiSeewIcffqi7RkVFBc6cOYNbt24BADw9PbFz504MHDgQHTt2xOTJkzFgwADs2LED7u7udvmcROQ63N0kSI4PBwCDEMjlXojIFiRCCHP2HScjiouLIZfLoVara1xvkIhcU/VFno3N5DW2DqCKe/wSWR2/vxvBMjBERI7InHAXF6FCbLiyxpBIRGQN7AFsAP4LgoiM0S7yfPf/uWpjHRd0JrIvfn87wT2ARETOpLZFngEgNT0Hmkr+25uI7IcBkIjIgsxd5Dkz95rtiiIiugsDIBGRBXGRZyJyBgyAREQWxEWeicgZMAASEVkQF3kmImfAAEhEVEeaSoGD54uw6cQlHDxfpDehg4s8E5Ez4DqARER1YO76fmkJUQbtlFzkmYgcBNcBbACuI0TkWuq6vl9tO4EQkX3w+5s9gEREZqltfT8Jqtb3iw1X6kKeu5sE0W39bVkmEZFZeA8gEZEZuL4fETUmDIBERGbg+n5E1JhwCJiIqBpT9+1xfT8iakwYAImI/r+aZvjGhiuhkstQoC41eh+gBFWzfLm+HxE5Aw4BExHhfzN8777Pr0BdiqS1WdieU8D1/Yio0WAAJCKXV9sMX+B/M3zTEqKglOsP8yrlMoMlYIiIHBmHgInI5dVlhm9chAqx4Uqu70dETo0BkIhcQk2LMtd1hi/X9yMiZ8cASESNXm3bt3GGLxG5Gt4DSESNWm2TOzKy89Ej1A8qucxgcoeWBFWBkTN8iaixYAAkokbL3MkdADjDl4hcCgMgETk9TaXAwfNF2HTiEg6eL4Kmsire1XVyB2f4EpGr4D2AROTUarq/r+xOpVnX0E7u4AxfInIVDIBE5LS09/fdPcSrvb9vav8OZl2n+uQOzvAlIlfAAEhEDqumpVtqu79PAmD9kTwom0vxZ3EZt28jIqqGAZCIHFJtS7eYe3/ftP7t8c8d5yAB9EIgJ3cQkSvjJBAisgtTEzcA85ZuMXfx5hBFU07uICK6C3sAicjmaurdiw1X1jq0m5qegw//736z3ivAR4botv6c3EFEVA0DIBFZXE337pkzccOcoV1IqkJjgbrUrPv7OLmDiOh/GACJyKIa2ru38kCuWe9ztaQMyfHhSFqbxfv7iIjqiPcAEpHZarpvD6j93r3Fu36ttXfvr1sVZtUS4CPj4s1ERPXEHkAiAlDzsC1Q+6xcc5ZlMbd3z9fLA+rbFWYN7XLxZiKiumMAJGrkagt2QO3hrrb79tISoiD38rRY796zvUPqtHQL7+8jIqobBkCiRqy2YKdtU1O4WzIqCm9vqX1W7isDO5pVkzm9e5P6tUdHpY9B7cq7aiciovphACRyYg2ZbZuWEGXWpIzZm7JRdLPcZA3aWbnXamhTnbm9exzaJSKyHgZAIgfU0GFbc9fS85F61DpsW1P4q86vmdSsZVnq0rvHoV0iIutw+AA4d+5cbNmyBSdOnICnpyf++usvgzZTpkzB/v37kZ2djU6dOuHEiRO1XresrAwzZ87EF198gdu3b+PRRx/F0qVL0apVK8t/CKJqGjrZQtvGEmvpHfztqsU+l7K5zOxlWdi7R0RkXw4fAMvLyzFs2DBER0djxYoVRtsIITB27FgcPnwYv/zyi1nXnTp1KtLT07F+/Xr4+/tjxowZGDJkCI4dOwZ3d3dLfgRyEbaabGHJtfT+F81q5tfUA9dv1j4r191NgrSEKPbuERE5OIcPgKmpqQCAVatWmWzzr3/9CwBw5coVswKgWq3GihUr8Nlnn6F///4AgLVr1yI4OBg7duzAwIEDG144NSq26LUzd7KFOcO25s62jW7rj2+y/qh12Hb24E548fPjZs3KZe8eEZHjc/gAaA3Hjh1DRUUFBgwYoDsWFBSEiIgIHDhwwGQALCsrQ1lZme55cXGx1Wul+jOnR86cdrbqtTN3soW5w7bmzLbtda+/WcO2cREqpLlJzJ6Vy949IiLH5pIBsKCgAJ6enmjRooXe8cDAQBQUFJg8b968eboeSao7cwKZrUKbue1s2Wtn7mQLc4dt6zLb1pxhW/bsERE1HnYJgCkpKbUGqSNHjqB79+42qqiKEAISiekvs9deew3Tp0/XPS8uLkZwcLDF66gp3FgqRFkqaJnbxtwhUluENu0WYZYId5butTOHucO2dZlta264Y88eEVHjYJcAOGnSJIwcObLGNiEhIVZ7f6VSifLycly/fl2vF7CwsBAxMTEmz5NKpZBKpVarC6g53ACwSIiyVNCqS5vaAhkAm4W21PQc9AsLtMiQrKV77cyZbGHusG1dZ9sy3BERuQ6JEMLYd43DWbVqFaZOnWp0GRitlJQUfPvtt7UuA6NWq9GyZUusXbsWw4cPBwDk5+ejVatW+P77782eBFJcXAy5XA61Wo3mzZub+1FMMhVu7v6Sv/s1ADWGqLq2qSlo1fVaseFKPPjeLpPDn9pQI4RAQXFZjW32vvwI+nywu8Zr+TX1NCuUzR7cCW9vOV1rO0tZ91xPzPzqZ7MnWwDGg5329wOYH+SJiEifpb+/nZHD3wOYl5eHa9euIS8vDxqNRhfu2rVrh2bNmgEAfv31V5SUlKCgoAC3b9/WtQkPD4enpycuXbqERx99FGvWrEGPHj0gl8sxbtw4zJgxA/7+/vDz88PMmTPRpUsX3axgW9NUihp7pEyp3qslhKi158ucNub0jpl7LXPufavp9eptPjt4wWL30f1+7ZZZ7cxhyV67uky24D15RERUXw4fAOfMmYPVq1frnkdGRgIAdu/ejb59+wIAnnvuOezdu9egTW5uLkJCQlBRUYEzZ87g1q3/fekvXLgQTZo0wfDhw3ULQa9atcpuawBm5l6rNQiZUpcQZU4bc4KWudey5L1vlgxtbfy8zWpnTrgzd4kUa0y24LAtERHVh8MHwFWrVtW4BiAA7Nmzp8bXQ0JCcPdIt0wmw6JFi7Bo0aIGVmgZhTfqF/6swZJBy9x738xhydA2OjoE/96fa5H176zRa8dgR0RE1uTwAdBVBPjI7F2CjrlByxzmzlgVQuDP4jKbhTbPJm4WHZJlrx0RETkTBkAH0SPUDyq5zGS4qUldQpQlg5Y51zL33jcAtbaxRmiz5JAsgx0RETkLp5kF7IisNQsYMAw3wsj/1j4H9GflGju/Lm2qzwK2xLW0n82Wy9MAll3rkIiIGg/OAmYAbBBr/AE1xnUAtWy9QDUREZExDIAMgA1irT+gxrgTCBERkaNgAGQAbBD+ARERETkffn8DbvYugIiIiIhsiwGQiIiIyMUwABIRERG5GAZAIiIiIhfDAEhERETkYhgAiYiIiFwMAyARERGRi2EAJCIiInIxDIBERERELqaJvQtwZtpNVIqLi+1cCREREZlL+73typuhMQA2wI0bNwAAwcHBdq6EiIiI6urGjRuQy+X2LsMuuBdwA1RWVuLy5cvw8fGBRCKx6LWLi4sRHByMixcvuuQ+ha7++QH+DFz98wP8GQD8Gbj65wes8zMQQuDGjRsICgqCm5tr3g3HHsAGcHNzQ6tWraz6Hs2bN3fZ/+gBfn6APwNX//wAfwYAfwau/vkBy/8MXLXnT8s1Yy8RERGRC2MAJCIiInIxDIAOSiqVIjk5GVKp1N6l2IWrf36APwNX//wAfwYAfwau/vkB/gyshZNAiIiIiFwMewCJiIiIXAwDIBEREZGLYQAkIiIicjEMgEREREQuhgHQCcydOxcxMTHw9vaGr6+vvcuxiaVLlyI0NBQymQzdunXDvn377F2Szfz444+Ij49HUFAQJBIJvv32W3uXZFPz5s3DAw88AB8fHwQEBOCJJ57AmTNn7F2WTaWlpeG+++7TLXwbHR2NrVu32rssu5k3bx4kEgmmTp1q71JsJiUlBRKJRO+hVCrtXZZNXbp0CQkJCfD394e3tze6du2KY8eO2busRoMB0AmUl5dj2LBhSEpKsncpNrFhwwZMnToVb7zxBo4fP46HHnoIgwYNQl5enr1Ls4mbN2/i/vvvx+LFi+1dil3s3bsXL774Ig4dOoTt27fjzp07GDBgAG7evGnv0mymVatWePfdd3H06FEcPXoU/fr1w9ChQ3Hq1Cl7l2ZzR44cwSeffIL77rvP3qXYXOfOnZGfn697nDx50t4l2cz169fRu3dveHh4YOvWrcjJycH8+fNdphPEFrgMjBNZtWoVpk6dir/++svepVhVz549ERUVhbS0NN2xTp064YknnsC8efPsWJntSSQSbNy4EU888YS9S7GbK1euICAgAHv37sXDDz9s73Lsxs/PDx988AHGjRtn71JspqSkBFFRUVi6dCn+8Y9/oGvXrvjnP/9p77JsIiUlBd9++y1OnDhh71Ls4tVXX8VPP/3kUqM/tsYeQHIo5eXlOHbsGAYMGKB3fMCAAThw4ICdqiJ7UqvVAKoCkCvSaDRYv349bt68iejoaHuXY1MvvvgiBg8ejP79+9u7FLs4d+4cgoKCEBoaipEjR+K3336zd0k2s3nzZnTv3h3Dhg1DQEAAIiMjsXz5cnuX1agwAJJDuXr1KjQaDQIDA/WOBwYGoqCgwE5Vkb0IITB9+nQ8+OCDiIiIsHc5NnXy5Ek0a9YMUqkUEyZMwMaNGxEeHm7vsmxm/fr1OHbsmMv1+mv17NkTa9aswQ8//IDly5ejoKAAMTExKCoqsndpNvHbb78hLS0N7du3xw8//IAJEyZg8uTJWLNmjb1LazQYAO3E2A2+dz+OHj1q7zLtRiKR6D0XQhgco8Zv0qRJ+OWXX/DFF1/YuxSb69ixI06cOIFDhw4hKSkJY8aMQU5Ojr3LsomLFy9iypQpWLduHWQymb3LsYtBgwbhb3/7G7p06YL+/ftjy5YtAIDVq1fbuTLbqKysRFRUFN555x1ERkbihRdewPjx4/VuDaKGaWLvAlzVpEmTMHLkyBrbhISE2KYYB6JQKODu7m7Q21dYWGjQK0iN20svvYTNmzfjxx9/RKtWrexdjs15enqiXbt2AIDu3bvjyJEj+Oijj7Bs2TI7V2Z9x44dQ2FhIbp166Y7ptFo8OOPP2Lx4sUoKyuDu7u7HSu0vaZNm6JLly44d+6cvUuxCZVKZdDj3alTJ3zzzTd2qqjxYQC0E4VCAYVCYe8yHI6npye6deuG7du348knn9Qd3759O4YOHWrHyshWhBB46aWXsHHjRuzZswehoaH2LskhCCFQVlZm7zJs4tFHHzWY8frss88iLCwMs2bNcrnwBwBlZWU4ffo0HnroIXuXYhO9e/c2WP7p7NmzaNOmjZ0qanwYAJ1AXl4erl27hry8PGg0Gt2ssHbt2qFZs2b2Lc4Kpk+fjtGjR6N79+6Ijo7GJ598gry8PEyYMMHepdlESUkJfv31V93z3NxcnDhxAn5+fmjdurUdK7ONF198EZ9//jk2bdoEHx8fXW+wXC6Hl5eXnauzjddffx2DBg1CcHAwbty4gfXr12PPnj3IyMiwd2k24ePjY3DPZ9OmTeHv7+8y94LOnDkT8fHxaN26NQoLC/GPf/wDxcXFGDNmjL1Ls4lp06YhJiYG77zzDoYPH47MzEx88skn+OSTT+xdWuMhyOGNGTNGADB47N69296lWc2SJUtEmzZthKenp4iKihJ79+61d0k2s3v3bqO/7zFjxti7NJsw9tkBiJUrV9q7NJsZO3as7u+/ZcuW4tFHHxXbtm2zd1l21adPHzFlyhR7l2EzI0aMECqVSnh4eIigoCDx1FNPiVOnTtm7LJtKT08XERERQiqVirCwMPHJJ5/Yu6RGhesAEhEREbkYzgImIiIicjEMgEREREQuhgGQiIiIyMUwABIRERG5GAZAIiIiIhfDAEhERETkYhgAiYiIiFwMAyARObS+ffti6tSpVn+f//73v+jVqxdkMhm6du1qshaJRAKJRKLbkcdSVq1apbu2LT4vEbk2BkAiIgDJyclo2rQpzpw5g507d5psN378eOTn5+u2JLtw4YLJQFiX8DpixAjk5+cjOjq6PuUTEdUJ9wImIgJw/vx5DB48uNbN5r29vaFUKi3+/l5eXvDy8oKnp6fFr01EdDf2ABKRU7l+/TqeeeYZtGjRAt7e3hg0aBDOnTun12b58uUIDg6Gt7c3nnzySSxYsAC+vr4mrymRSHDs2DG89dZbkEgkSElJsXjde/bs0Q3xVn8kJiZa/L2IiGrDAEhETiUxMRFHjx7F5s2bcfDgQQgh8Nhjj6GiogIA8NNPP2HChAmYMmUKTpw4gdjYWMydO7fGa+bn56Nz586YMWMG8vPzMXPmTIvXHRMTg/z8fN1j165dkMlkePjhhy3+XkREteEQMBE5jXPnzmHz5s346aefEBMTAwBYt24dgoOD8e2332LYsGFYtGgRBg0apAtxHTp0wIEDB/Ddd9+ZvK5SqUSTJk3QrFmzeg/vxsTEwM1N/9/Ut2/f1k0o8fT01F27qKgI48ePx9ixYzF27Nh6vR8RUUMwABKR0zh9+jSaNGmCnj176o75+/ujY8eOOH36NADgzJkzePLJJ/XO69GjR40B0BI2bNiATp066R17+umnDdpVVFTgb3/7G1q3bo2PPvrIqjUREZnCAEhETkMIYfK4RCIx+N+1nWdJwcHBaNeund4xLy8vg3ZJSUnIy8vDkSNH0KQJ/y+YiOyD9wASkdMIDw/HnTt3cPjwYd2xoqIinD17Vtf7FhYWhszMTL3zjh49atM6TVmwYAE2bNiAzZs3w9/f397lEJELYwAkIqfRvn17DB06FOPHj8f+/fvx888/IyEhAffccw+GDh0KAHjppZfw/fffY8GCBTh37hyWLVuGrVu3GvQK2tqOHTvwyiuv4MMPP4RCoUBBQQEKCgqgVqvtWhcRuSYGQCJyKitXrkS3bt0wZMgQREdHQwiB77//Hh4eHgCA3r174+OPP8aCBQtw//33IyMjA9OmTYNMJrNr3fv374dGo8GECROgUql0jylTpti1LiJyTRJhi5tjiIjsaPz48fjvf/+Lffv2Neg6ffv2RdeuXfHPf/7TMoXZ6T2IiNgDSESNzocffoiff/4Zv/76KxYtWoTVq1djzJgxFrn20qVL0axZM5w8edIi19Nat24dmjVr1uCQSkRkDvYAElGjM3z4cOzZswc3btzAvffei5deegkTJkxo8HUvXbqE27dvAwBat25t0W3bbty4gT///BMA4OvrC4VCYbFrExHdjQGQiIiIyMVwCJiIiIjIxTAAEhEREbkYBkAiIiIiF8MASERERORiGACJiIiIXAwDIBEREZGLYQAkIiIicjEMgEREREQuhgGQiIiIyMX8P1syGX1CZcpjAAAAAElFTkSuQmCC" } }, "cell_type": "markdown", "id": "f1d5e828", "metadata": {}, "source": [ "![data_example-2.png](attachment:data_example-2.png)" ] }, { "cell_type": "code", "execution_count": 4, "id": "39511da9", "metadata": {}, "outputs": [], "source": [ "# instantiate the HN module\n", "cond = dd.Conductivity()" ] }, { "cell_type": "markdown", "id": "648e6324-33ae-48b0-8b24-389481a63f58", "metadata": {}, "source": [ "**Select the region of interest (ROI) to fit the data using the select range method**.
\n", "The data in ROI is shown as image in the next cell" ] }, { "cell_type": "code", "execution_count": 5, "id": "75be026b", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "x_lower_limit -1.039386875850182 x_upper_limit 5.614181794555549\n" ] } ], "source": [ "#select range\n", "x1,y1 = cond.select_range(x,y)" ] }, { "cell_type": "markdown", "id": "1a5ad9db", "metadata": {}, "source": [ "**Plot of the ROI for fitting**" ] }, { "cell_type": "code", "execution_count": 6, "id": "6f00ed52", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'log ( $\\\\sigma´$)')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "plt.scatter(x1,y1)\n", "plt.xlabel('log f [Hz]')\n", "plt.ylabel('log ( $\\sigma´$)')" ] }, { "attachments": { "ROI.png": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAABDRUlEQVR4nO3de3RU1d3/8c8kQCYBDIQGJuAtBQKMiRLAxBStl8fUtJiiomLLrbZoRYpWfqJV0ZBaxRbpoz4qKt6WBmur1dYIBpcVtbWRAAFMDMilQRAmSU0gBJLhMjO/P9KkDDNJJmHu5/1ai9XOOTvnfGcb6qd7n72PyeVyuQQAAADDiAl1AQAAAAguAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAbTK9QFRLLjx4+rsbFRcXFxiokhSwMAEAmcTqeOHDmixMRE9eplzChkzG/tJ42Njdq1a1eoywAAAD1w9tlna9CgQaEuIyQIgKcgLi5OUusvUHx8vF+v7XA4tG3bNqWlpSk2Ntav145U9Ik7+sMTfeKJPvFEn3gyWp+0tLRo165d7f8eNyIC4Clom/aNj49XQkKCX6/tcDgkSQkJCYb4y+gL+sQd/eGJPvFEn3iiTzwZtU+M/PiWcb85AACAQREAAQAADIYACAAAYDAEQAAAAIMhAAIAABhM1AfAPXv26KabblJWVpZycnJ011136eDBgx2237x5s6677jplZmbqsssu0xtvvBHEagEAAAIv6gPg/PnzNWLECH366ad67733tG/fPj3yyCNe2zY2Nurmm2/WVVddpXXr1umhhx7S4sWL9fnnnwe5agAAgMCJ+gC4c+dOuVyu9j8mk6nDTZvff/99DRgwQNOmTVOvXr2Uk5Oj/Px8rVixIshVAwAABE7EbwRtt9tVW1vr9VxycrLmzZunpUuX6pVXXpHD4dDYsWN15513em2/fft2paWluR0bMWKE3nzzTb/XDQAAvHM4XSqrblBdk12D+5uVlZqk2BhTqMuKKhEfADdv3qyZM2d6PffUU0/JZDJpzpw5uvHGG7V//37Nnz9fDzzwgJYsWeLR/vDhwx6jg2azWc3NzZ3W4HA42ndR95e26/n7upGMPnFHf3iiTzzRJ57oE0/h1Cerv6jRr9/dopqDR9qPWU6L0wNXjtEV51j8co9w+J6hZnK5XK5QFxEolZWVmjFjhtatW6devVqz7oYNGzRt2jStX79e/fr1c2v/m9/8RnV1dXriiSfaj7366qv685//rL/85S8e129ubtaWLVsC+h0AADCKz762a0npgQ7PL8gZoAtON/vtfmPGjPH7q1wjRcSPAHbGZrPJ4XDI6XS2H+vdu7dMJpPXdx2mpaXp008/dTu2Y8cOjRw5stP7pKWlBeRdwBUVFcrIyDDUexk7Q5+4oz880See6BNP9ImncOgTh9Oluas/6vC8SVJRlV03/SD7lKeDm5ubtW3btlO6RqSL6gA4fvx4xcfH6+GHH9Y999yjQ4cOaenSpcrNzfW6ECQ3N1dLlizRyy+/rGnTpmnDhg0qLi7W008/3el9YmNjA/YXJpDXjlT0iTv6wxN94ok+8USfeApln5Ttqneb9j2ZS5Kt0a4NuxuVM3zQKd2Lf+5Rvgo4KSlJL7zwgnbt2qWLLrpIV111lc4++2w9/PDD7W0mTZqkZ555RpI0cOBAvfjiiyopKVF2drYWLlyohQsX6oILLgjVVwAAwBDqmux+bYfORfUIoCSlp6fr5Zdf7vD8ypUr3T5nZGTo9ddfD3BVAADgRIP7+/Zsn6/t0LmoHgEEAACRISs1SSmJZnX0dJ9JUkpi65YwOHUEQAAAEHKxMSYV5FslySMEtn0uyLeyH6CfEAABAEBYyEtP0bLp42RJdJ/mtSSatWz6OOWlp4SosugT9c8AAgCAyJGXnqJcq4U3gQQYARAAAISV2BjTKW/1gs4xBQwAAGAwBEAAAACDIQACAAAYDAEQAADAYAiAAAAABkMABAAAMBi2gQEAAEHhcLrY3y9MEAABAEDAlVTaVFhcJVujvf1YSqJZBflW3vARAkwBAwCAgCqptGlOUblb+JOkmka75hSVq6TSFqLKjIsACAAAAsbhdKmwuEouL+fajhUWV8nh9NYCgUIABAAAAVNW3eAx8ncilyRbo11l1Q3BKwoEQAAAEDh1TR2Hv560g38QAAEAQMAM7m/2azv4BwEQAAAETFZqklISzeposxeTWlcDZ6UmBbMswyMAAgCAgImNMakg3ypJHiGw7XNBvpX9AIOMAAgAAAIqLz1Fy6aPkyXRfZrXkmjWsunj2AcwBNgIGgAABFxeeopyrRbeBBImCIAAACAoYmNMyhk+KNRlQEwBAwAAGA4BEAAAwGAIgAAAAAZDAAQAADAYFoEAAIBT4nC6WN0bYQiAAACgx0oqbSosrpKt8b/v8k1JNKsg38r+fmGMKWAAANAjJZU2zSkqdwt/klTTaNeconKVVNpCVBm6QgAEAADd5nC6VFhcJZeXc23HCour5HB6a4FQIwACAIBuK6tu8Bj5O5FLkq3RrrLqhuAVBZ8RAAEAQLfVNXUc/nrSDsFFAAQAAN02uL/Zr+0QXARAAADQbVmpSUpJNKujzV5Mal0NnJWaFMyy4CMCIAAA6LbYGJMK8q2S5BEC2z4X5FvZDzBMEQABAECP5KWnaNn0cbIkuk/zWhLNWjZ9HPsAhjE2ggYAAD2Wl56iXKuFN4FEGAIgAAA4JbExJuUMHxTqMtANTAEDAAAYDAEQAADAYAiAAAAABkMABAAAMBgCIAAAgMGwChgAAHjlcLrY3iVKEQABAICHkkqbCourZGu0tx9LSTSrIN/KBs9RIOqngPfs2aObbrpJWVlZysnJ0V133aWDBw922L6goEDp6enKzMxs//PHP/4xiBUDABBaJZU2zSkqdwt/klTTaNeconKVVNpCVBn8JeoD4Pz58zVixAh9+umneu+997Rv3z498sgjHbavqKjQgw8+qI0bN7b/mTp1ahArBgAgdBxOlwqLq+Tycq7tWGFxlRxOby0QKaI+AO7cuVMul6v9j8lkUnx8vNe2R48e1bZt25Senh7kKgEACA9l1Q0eI38nckmyNdpVVt0QvKLgdxH/DKDdbldtba3Xc8nJyZo3b56WLl2qV155RQ6HQ2PHjtWdd97ptf3WrVt1/PhxPfHEE9qwYYP69++vKVOmaPbs2YqJ6TgrOxwOORwOv3yfE6954n+CPjkZ/eGJPvFEn3iiTzyd2Cc1jc0+/UxNY7McjgEBrCpw+GcvmVwuV0SP4a5du1YzZ870eu6pp57S119/rcOHD+vGG2/U/v37NX/+fJ155plasmSJR/tPP/1Uzz77rH7xi18oMzNTW7Zs0dy5czVr1izNnj3bo31zc7O2bNni9+8EAECoVNYdUcHH+7tsV3jxQKUPjgtCRYEzZswYJSQkhLqMkIj4ANiZyspKzZgxQ+vWrVOvXq2DnRs2bNC0adO0fv169evXr8trPP/881q1apXeeustj3NtATAtLc3vv0AOh0MVFRXKyMhQbGysX68dqegTd/SHJ/rEE33iiT7xdGKfyBSj7y75SLUHj3h9DtAkyZJo1sd3XhyxW8I0Nzdr27Zthg6AET8F3BmbzSaHwyGn09l+rHfv3jKZTF7/0n/wwQf65ptvdMMNN7QfO3r0qMxmc6f3iY2NDdj/iATy2pGKPnFHf3iiTzzRJ57oE09tfbLoh+doTlG5TJJbCGyLewX5VvXpHbkRgn/uUb4IZPz48YqPj9fDDz+sI0eOqL6+XkuXLlVubq7XhSAul0uLFy9WaWmpXC6XNm7cqFdeeYVVwAAAQ8lLT9Gy6eNkSXQfALEkmrVs+jj2AYwCkRvffZCUlKQXXnhBjz76qC666CLFxcXpsssu04IFC9rbTJo0Sfn5+brllluUm5ure+65R4sWLVJtba2+9a1vad68eZo8eXIIvwUAAIHjcLr02b/qtW53i+yn1euC4cmKjTEpLz1FuVYLbwKJUlEdACUpPT1dL7/8cofnV65c6fb5hhtucJsCBgAgWnm87WPtOre3fcTGmJQzfFBoi0RARPUUMAAA8I63fRgbARAAAIPhbR8gAAIAYDC87QMEQAAADKamscWndnVNHYdERDYCIAAABlJSadODK317i9Xg/p3vg4vIFfWrgAEAQKu2hR9dPdnX9raPrNSkYJSFEGAEEAAAA+hs4Yc3BflW9vyLYgRAAAAMoKuFH22S+vbmbR8GQAAEAMAAfF3Qcf+V5xD+DIAACACAAfi6oMNyGgs/jIAACACAAWSlJikl0azOnupLYeGHYRAAAQAwgNgYkwryrZLkEQLbPt8/aTQLPwyCAAgAgEHkpado2fRxsiS6T/NaEs1akDNAV5xjCVFlCDb2AQQAwEDy0lOUa7WorLpBdU12De5v1vgzE1Xx+eZQl4YgIgACAGAwsTEm5Qwf1P7Z4XCEsBqEAlPAAAAABkMABAAAMBgCIAAAgMHwDCAAAFHC4XS5Le7ISk1iWxd4RQAEACAKlFTaVFhc5fa+35REswryrbzaDR6YAgYAIMKVVNo0p6jcLfxJUk2jXXOKylVSaQtRZQhXBEAAACKYw+lSYXGVXF7OtR0rLK6Sw+mtBYyKAAgAQAQrq27wGPk7kUuSrdGusuqG4BWFsEcABAAggtU1dRz+etIOxkAABAAggg3ub+66UTfawRgIgAAARLCs1CSlJJrV0WYvJrWuBs5KTQpmWQhzBEAAACJYbIxJBflWSfIIgW2fC/Kt7AcINwRAAAAiXF56ipZNHydLovs0ryXRrGXTx7EPIDywETQAAFEgLz1FuVYLbwKBTwiAAABEidgYk3KGDwp1GYgATAEDAAAYDAEQAADAYAiAAAAABsMzgAAAhDmH08XiDvgVARAAgDBWUmlTYXGV2/t+UxLNKsi3sr0LeowpYAAAwlRJpU1zisrdwp8k1TTaNaeoXCWVthBVhkhHAAQAIAw5nC4VFlfJ5eVc27HC4io5nN5aAJ0jAAIAEIbKqhs8Rv5O5JJka7SrrLoheEUhahAAAQAIQ3VNHYe/nrQDTkQABAAgDA3ub+66UTfaASciAAIAEIayUpOUkmhWR5u9mNS6GjgrNSmYZSFKEAABAAhDsTEmFeRbJckjBLZ9Lsi3sh8geoQACABAmMpLT9Gy6eNkSXSf5rUkmrVs+jj2AUSPsRE0AABhLC89RblWC28CgV8RAAEACHOxMSblDB8U6jIQRaJ+Cri2tla33XabsrOzdeGFF2rx4sU6cuRIh+03b96s6667TpmZmbrsssv0xhtvBLFaAACAwIvqAOh0OnXrrbfq6NGjKikp0TvvvKOtW7dq0aJFXts3Njbq5ptv1lVXXaV169bpoYce0uLFi/X5558Ht3AAgCE4nC6V7qzXXzftVenOet7qgaCJ6ing6upqVVZWas2aNRo4cKAkaf78+Zo2bZruvfde9e/f3639+++/rwEDBmjatGmSpJycHOXn52vFihU699xzg14/ACB6lVTaVFhc5fa2j5REswryrSzuQMBFfAC02+2qra31es7pdEqS4uPj24+ZTCYdO3ZMe/bskdVqdWu/fft2paWluR0bMWKE3nzzTT9XDQAwspJKm+YUlXu857em0a45ReWs8EXARXwA3Lx5s2bOnOn13BNPPKGRI0dq8eLFWrhwoY4fP64nn3xSUmtwPNnhw4fdwqIkmc1mNTc3d1qDw+GQw+Ho4Tfo+Jon/ifok5PRH57oE0/0iadQ94nD6dKid77wCH9S6/t9TZIKi6t02ajkoK30DXWfBJtRvmdnIj4AZmdn68svv+zwvNVq1UMPPaTvfe97Sk5O1o033qiPP/5YiYmJHm3j4+PV1NTkdsxut6tv376d1rBt27aeFe+DioqKgF07UtEn7ugPT/SJJ/rEU6j6pLLuiGoOdrwY0SXJ1mjXHz4oU/rguOAVJn5PjCTiA2BnXC6XGhsb9dhjj8lsbt1E8+OPP1bfvn111llnebRPS0vTp59+6nZsx44dGjlyZKf3SUtLU0JCgv8KV+v/O6moqFBGRoZiY2P9eu1IRZ+4oz880See6BNPoe6T3Zv3SdrfZbvThpyhsecNDXxBCn2fBFtzc3NAB28iQVQHQJPJpLvuukt5eXn6xS9+od27d+vRRx/VjBkz1KuX51fPzc3VkiVL9PLLL2vatGnasGGDiouL9fTTT3d6n9jY2ID9hQnktSMVfeKO/vBEn3iiTzyFqk8sib4NGFgSE4Jen1F+T4zwHbsS1dvASNJjjz2mtWvXasKECZo5c6Zyc3N1++23t5+fNGmSnnnmGUnSwIED9eKLL6qkpETZ2dlauHChFi5cqAsuuCBU5QMAokxWapJSEs0e7/dtY1LrauCs1KRglgWDieoRQKl1enbFihUdnl+5cqXb54yMDL3++uuBLgsAYFCxMSYV5Fs1p6hcJsltMUhbKCzIt/KqNwRU1I8AAgAQbvLSU7Rs+jhZEs1uxy2JZraAQVBE/QggAADhKC89RblWi8qqG1TXZNfg/q3Tvoz8IRgIgAAAhEhsjEk5wweFugwYEFPAAAAABsMIIAAAfuZwupjaRVgjAAIA4EcllTYVFlfJ1vjfV46mJJpVkG9lcQfCBlPAAAD4SUmlTXOKyt3CnyTVNNo1p6hcJZW2EFUGuCMAAgDgBw6nS4XFVW77+rVpO1ZYXCWH01sLILgIgAAA+EFZdYPHyN+JXJJsjXaVVTcEryigAwRAAAD8oK6p4/DXk3ZAIBEAAQDwg8H9zV036kY7IJAIgAAA+EFWapJSEs3qaLMXk1pXA2elJgWzLMArAiAAAH4QG2NSQb5VkjxCYNvngnwr+wEiLBAAAQDwk7z0FC2bPk6WRPdpXkuiWcumj2MfQIQNNoIGAMCP8tJTlGu18CYQhDUCIAAAfhYbY1LO8EGhLgPoEFPAAAAABsMIIAAAPnI4XUztIioQAAEA8EFJpU2FxVVub/tISTSrIN/K4g5EHKaAAQDoQkmlTXOKyj1e9VbTaNeconKVVNpCVBnQMwRAAAA64XC6VFhcJZeXc23HCour5HB6awGEJwIgAACdKKtu8Bj5O5FLkq3RrrLqhuAVBZwiAiAAAJ2oa+o4/PWkHRAOCIAAAHRicH9z14260Q4IBwRAAAA6kZWapJREs8f7fduY1LoaOCs1KZhlAaeEAAgAQCdiY0wqyLdKkkcIbPtckG9lP0BEFAIgAABdyEtP0bLp42RJdJ/mtSSatWz6OPYBRMRhI2gAAHyQl56iXKuFN4EgKhAAAQDwUWyMSTnDB4W6DOCUMQUMAABgMARAAAAAgwn4FPCBAwcUGxur/v37B/pWAAD0iMPp4tk+GEpAAuAnn3yi1157TWvXrpXd3rozekJCgi666CJdd911mjhxYiBuCwBAt5VU2lRYXOX2ureURLMK8q2s7kXU8msA3Ldvn+666y7t379fP/jBDzRr1iwNGTJEDodDdXV12rBhgx566CENGjRIjzzyiIYNG+bP2wMA0C0llTbNKSqX66TjNY12zSkqZ4sXRC2/BsD58+fr9ttvV05Ojse5kSNHauLEibrtttv08ccfa/78+frjH//oz9sDAOAzh9OlwuIqj/AnSS61bvJcWFylXKuF6WBEHb8GwKKiIvXq1fUlL774YqaBAQAhtW5Xg9u078lckmyNdpVVN7D1C6KOX1cB+xL+etIWAAB/q2s64mO7jkMiEKlOOQCWlpbq8OHD/qgFAICgGdw/zsd25q4bARHmlAPgvn37NG3aNP373//2Rz0AAASMw+nSZ/+q1993t8jpkiynxamjp/tMal0NnJWaFMwSgaA45XnYKVOmKD4+Xj/96U9VXFzsj5oAAPA7j+1e1q7TgITe7Qs+TlwM0hYKC/KtLABBVPLLg3g/+MEPNGbMGH9cCgAAv+tou5fG5mOSpMSE3jrwn/8uSRb2AUSU89tKjNTUVH9dCgAAv/Flu5f43rF66mfj9M3hI7wJBIbg96W4hw4dUr9+/fx9WQAAeqSs2rftXmJiTJo8lhcUwBj8ug2MJJ1//vlavXq1vy8LAECP+LqNC9u9wEj8PgLocrn0hz/8Qc8//7xMJpPS09M1efJknXfeef6+FQAAXfJ1Gxe2e4GR+H0EUJK2bt2qc889V1lZWaqurtb06dP18MMPB+JWXaqtrdVtt92m7OxsXXjhhVq8eLGOHOl488+CggKlp6crMzOz/Q+vrAOAyORwuuR0uTQgvneHbdjuBUYUkNdxLF261O1Vb19++aXmzp0ri8Win/70p4G4pVdOp1O33nqrkpOTVVJSIpfLpTvuuEOLFi3S4sWLvf5MRUWFHnzwQV199dVBqxMA4H8e2754wXYvMCq/jwAOGDBAFovF7dioUaO0cOFCvf766/6+Xaeqq6tVWVmpBx54QAMHDlRSUpLmz5+v4uJiNTU1ebQ/evSotm3bpvT09KDWCQDwr7ZtXzoLf1Lrdi/Lpo9juxcYjt9HAEePHq0333xTd999t9vxM888Uzabzd+3k91uV21trddzTqdTkhQfH99+zGQy6dixY9qzZ4+sVqtb+61bt+r48eN64okntGHDBvXv319TpkzR7NmzFRMTkNlyAICfdbbtS5t+fUxaNm28vjNyMCN/MCS/B8Bf/vKXmjlzpurq6vSjH/1Io0ePlt1u17PPPqvTTz/d37fT5s2bNXPmTK/nnnjiCY0cOVKLFy/WwoULdfz4cT355JOSWoPjyZqampSVlaUZM2bo97//vbZs2aK5c+cqJiZGs2fP7rAGh8Mhh8Phny90wjVP/E/QJyejPzzRJ56M2Cef/au+y5G/Q0ddklySyykDdU2HjPZ7YpTv2RmTy+Xq7P8k9cjWrVv10EMPacOGDWq7fFxcnB5//HFdfPHF/r5dp/bs2aOHHnpImzZtUnJysm688Ubdc889WrVqlYYPH97lzz///PNatWqV3nrrLY9zzc3N2rJlSyDKBgD00N93t+ixtY1dtvtldqIuOjO+y3aIXmPGjFFCQkKoywiJgCwCGT16tF599VXV19ersrJSLpdL5557rpKSgrvCyuVyqbGxUY899pjM5tbl/R9//LH69u2rs846y6P9Bx98oG+++UY33HBD+7GjR4+2/2xH0tLS/P4L5HA4VFFRoYyMDMXGxvr12pGKPnFHf3iiTzwZsU/sp9VLa9d12W68dYTGjkgOQkXhz2i/J83Nzdq2bVuoywipgATANoMGDQr6iN+JTCaT7rrrLuXl5ekXv/iFdu/erUcffVQzZsxQr16eX93lcmnx4sU666yzdMEFF2jTpk165ZVXdM8993R6n9jY2ID9hQnktSMVfeKO/vBEn3gyUp9cMDxZKYlm1TTavT4HaJKUFB+j7G9/yzB94iuj/J4Y4Tt2JepXNjz22GNau3atJkyYoJkzZyo3N1e33357+/lJkybpmWeekSTl5ubqnnvu0aJFi5SZmakFCxZo3rx5mjx5cqjKBwB0U2yMSQX5rYv8Tl7e0fb5p2NPY/EHDC2gI4DhIC0tTStWrOjw/MqVK90+33DDDW5TwACAyJOXnqJl08d57ANoSTTr/kmjNeRYTQirA0LPrwHQ4XD4PKzanbYAAHRXXnqKcq0WlVU3qK7JrsH9//O2D5dTmzYRAGFsfp0CnjZtmkpLS7ts98knn2jatGn+vDUAAB5iY0zKGT5Ik8cOU87wQUz7Av/h1xHAJUuW6J577tFvfvMbXXnllcrMzNSQIUPkdDpVV1enDRs2qKSkRImJifrd737nz1sDAADAR34NgGeccYaKior00Ucf6Q9/+IOee+45tbS0SGp9G8eFF16oO++8U5dccok/bwsAAIBuCMgikEsuuUSXXHKJXC6X9u/fr5iYGA0YMCAQtwIAAEA3BXQVsMlkCvrmzwAAAOhc1G8DAwCIPg6ny2N1Lws8AN8RAAEAEaWk0uaxv19KolkF+VblpaeEsDIgckT9m0AAANGjpNKmOUXlbuFPkmoa7ZpTVK6SSluIKgMiCwEQABARHE6XCourvL7ft+1YYXGVHE5vLQCcKCBTwJdddplMJs9nMXr37q2kpCRdeuml+tnPfqaYGPInAMA3ZdUNHiN/J3JJsjXaVVbdoJzhg4JXGBCBAhIAr7/+ev3pT3/S7NmzdcYZZ2jv3r168cUX9Z3vfEff/va39dprr8lut2vevHmBuD0AIArVNXUc/nrSDjCygATA1atX69lnn9XIkSPbj2VlZenOO+/UAw88oEsuuUQzZswgAAIAfDa4v9mv7QAjC8gc7FdffaWzzz7b7dgZZ5yh6upqSdLpp5+ugwcPBuLWAIAolZWapJREszra7MWk1tXAWansPwt0JSABcPTo0Xr22Wfdjr344osaMWKEJOmTTz7RsGHDAnFrAECUio0xqSDfKkkeIbDtc0G+lf0AAR8EZAp44cKFuummm/TGG28oJSVF+/btk8vl0rJly1ReXq65c+fq8ccfD8StAQBRLC89Rcumj/PYB9DCPoBAtwQkAFqtVq1evVpr1qxRTU2NkpOTlZeXJ7PZrAMHDuijjz7SoEGs0AIAdF9eeopyrRbeBAKcgoAEwEOHDunuu+/W3/72N0mt7wR+55139Nhjj2nAgAGBuCUAwEBiY0xs9QKcgoA8A7h06VIdPnxYK1eu1ObNm/XXv/5VTqdTS5YsCcTtAAAA0A0BCYBr1qzR0qVLNXz4cMXFxSktLU1LlizRBx98EIjbAQAAoBsCEgBbWlrUv39/t2OnnXaanE5nIG4HAACAbghIADzvvPP0+OOPy+VqfR+jy+XS448/royMjEDcDgAAAN0QkEUg/+///T/NnDlT77zzjoYNG6a9e/fKZDLppZdeCsTtAAAA0A0BCYCjRo3S6tWr9cEHH6ihoUHDhg3TxRdfrH79+gXidgCAKOBwutjaBQgSvwbAffv2uX3+zne+0/7fDx48qIMHD2ro0KH+vCUAIAqUVNo8NndOYXNnIGD8GgAvu+wymUze/9+ay+WSyWTSli1b/HlLAECEK6m0aU5RuVwnHa9ptGtOUbmWTR9HCAT8zK8BsG3jZwAAfOFwulRYXOUR/iTJpdZ3/BYWVynXamE6GPAjvwbAYcOG+fNyAIAoV1bd4DbtezKXJFujXWXVDbz5A/CjgGwDAwCAL+qaOg5/PWkHwDcEQABAyAzub/ZrOwC+IQACAEImKzVJKYlmdfR0n0mtq4GzUpOCWRYQ9QiAAICQiY0xqSDfKkkeIbDtc0G+lQUggJ8RAAEAIZWXnqJl08fJkug+zWtJNLMFDBAgAXkTCAAA3ZGXnqJcq4U3gQBBQgAEAISF2BgTW70AQcIUMAAAgMEQAAEAAAyGAAgAAGAwBEAAAACDIQACAAAYDKuAAQAB5XC62N4FCDMEQABAwJRU2lRYXCVbo739WEqiWQX5VjZ4BkKIKWAAQECUVNo0p6jcLfxJUk2jXXOKylVSaQtRZQAIgAAAv3M4XSosrpLLy7m2Y4XFVXI4vbUAEGgEQACA35VVN3iM/J3IJcnWaFdZdUPwigLQjgAIAPC7uqaOw19P2gHwr6gKgC0tLZo6dareeustt+PV1dWaNWuWMjMzdeGFF+qZZ57p9Dpvv/22cnNzNXbsWF1zzTXauHFjIMsGgKgzuL/Zr+0A+FfUBMDt27dr2rRp2rRpk9vxY8eO6ZZbblFGRobWrl2r5557TitWrNB7773n9Tpr167Vgw8+qEceeUTr1q3TD3/4Q82ZM0ctLS1B+BYAEB2yUpOUkmhWR5u9mNS6GjgrNSmYZQH4j6gIgKWlpZo1a5auvvpqDR061O3cunXrVFdXp9tuu019+vSR1WrVjBkztGLFCq/XeuONNzRp0iSNHz9evXv31k9+8hMNHDhQq1atCsZXAYCoEBtjUkG+VZI8QmDb54J8K/sBAiESEfsA2u121dbWej2XnJys0aNHa82aNYqLi9NLL73kdn779u1KTU1Vnz592o+NGDFCzz33nNfr7dixQ1OmTHE7NmLECG3durXD+hwOhxwOh69fxydt1/P3dSMZfeKO/vBEn3gKZZ/kjhmsp348Vr9+d4tqDh5pP25JNOv+SaOVO2ZwSOri98ST0frEKN+zMxERADdv3qyZM2d6PffUU0/p8ssv7/BnDx8+rPj4eLdj8fHxam5u9rm92WzusL0kbdu2rcNzp6qioiJg145U9Ik7+sMTfeIpVH0yRNIT3xugLf8+qv12pwaaYzQmuY9ij9Vo06aakNTUht8TT/SJcUREAMzOztaXX37Zo59NSEjweH6vpaVFffv29do+Pj5edrv7qjS73a6BAwd2eI+0tDQlJCT0qL6OOBwOVVRUKCMjQ7GxsX69dqSiT9zRH57oE0/h0ifjQ3ZnT+HSJ+HEaH3S3Nwc0MGbSBARAfBUjBw5Urt27dLx48fVq1fr192xY4dGjhzZYfvt27e7HduxY4e++93vdniP2NjYgP2FCeS1IxV94o7+8ESfeKJPPNEnnozSJ0b4jl2JikUgncnOztbAgQO1dOlSHTlyRFu3btWrr76qa6+91mv7a6+9VsXFxfrss8907Ngxvfzyy6qvr1dubm6QKwcAAAiMqA+AvXr10osvvqht27Zp4sSJuvnmmzVjxgxdc801kqR9+/YpMzNT69evlyTl5OSooKBAixYtUlZWllauXKnly5drwIABIfwWAAAA/hN1U8Affvihx7GzzjpLL7zwgtf2Q4cO9djoefLkyZo8eXJA6gMAAAi1qB8BBAAAgLuoGwEEAASHw+lSWXWD6prsGty/9a0ebOwMRAYCIACg20oqbSosrpKt8b/bZqUkmlWQb1VeekoIKwPgC6aAAQDdUlJp05yicrfwJ0k1jXbNKSpXSaUtRJUB8BUBEADgM4fTpcLiKrm8nGs7VlhcJYfTWwsA4YIACADwWVl1g8fI34lckmyNdpVVNwSvKADdRgAEAPisrqnj8NeTdgBCgwAIAPDZ4P5mv7YDEBoEQACAz7JSk5SSaFZHm72Y1LoaOCs1KZhlAegmAiAAwGexMSYV5FslySMEtn0uyLeyHyAQ5giAAIBuyUtP0bLp42RJdJ/mtSSatWz6OPYBBCIAG0EDALotLz1FuVYLbwIBIhQBEADQI7ExJuUMHxTqMgD0AFPAAAAABkMABAAAMBgCIAAAgMEQAAEAAAyGAAgAAGAwrAIGALhxOF1s7wJEOQIgAKBdSaVNhcVVsjXa24+lJJpVkG9lg2cgijAFDACQ1Br+5hSVu4U/SapptGtOUblKKm0hqgyAvxEAAQByOF0qLK6Sy8u5tmOFxVVyOL21ABBpCIAAAJVVN3iM/J3IJcnWaFdZdUPwigIQMARAAIDqmjoOfz1pByC8EQABABrc3+zXdgDCGwEQAKCs1CSlJJrV0WYvJrWuBs5KTQpmWQAChAAIAFBsjEkF+VZJ8giBbZ8L8q3sBwhECQIgAECSlJeeomXTx8mS6D7Na0k0a9n0cewDCEQRNoIGALTLS09RrtXCm0CAKEcABAC4iY0xKWf4oFCXASCAmAIGAAAwGAIgAACAwRAAAQAADIZnAAHAIBxOF4s7AEgiAAKAIaz+okYPrtzq9r7flESzCvKtbO8CGBBTwAAQ5T772q65r21yC3+SVNNo15yicpVU2kJUGYBQIQACQBRzOF16cdNBubycaztWWFwlh9NbCwDRigAIAFFs3a4G1bc4OzzvkmRrtKusuiF4RQEIOQIgAEQph9Olf+6s96ltXZO960YAogaLQAAgCpVU2lRYXOXx3F9HBvc3d90IQNQgAAJAlCmptGlOUbnX5/5OZpJkSWzdEgaAcTAFDABRxOF0qbC4yufwJ0kF+Vb2AwQMhhFAAIgiZdUNPk/7WtgHEDAsAiAARBFfF3P84tIRuiM3jZE/wKCYAgaAKOLrYo6JI75F+AMMLKoCYEtLi6ZOnaq33nrL7Xh1dbVmzZqlzMxMXXjhhXrmmWc6vc73v/99nXfeecrMzGz/s3PnzkCWDgB+kZWapJREszqKdia1vgKORR+AsUVNANy+fbumTZumTZs2uR0/duyYbrnlFmVkZGjt2rV67rnntGLFCr333nter3Po0CFVV1dr1apV2rhxY/uf4cOHB+FbAMCpiY0xqSDfKkkeIZBFHwDaREUALC0t1axZs3T11Vdr6NChbufWrVunuro63XbbberTp4+sVqtmzJihFStWeL1WZWWlBgwYoGHDhgWjdAA4ZQ6nS6U76/XXTXtVurNeuVaLlk0fJ0ui+3SwJdGsZdPHsegDQGQsArHb7aqtrfV6Ljk5WaNHj9aaNWsUFxenl156ye389u3blZqaqj59+rQfGzFihJ577jmv16uoqFB8fLymT5+u7du3a9iwYZo3b54uvfRS/30hAPATbxs+p/xnde8/7r5Mn+38t9ZVbtP56Wm6YHgyI38AJEVIANy8ebNmzpzp9dxTTz2lyy+/vMOfPXz4sOLj492OxcfHq7m52Wt7k8mkjIwMzZ8/X0OHDlVJSYnmzZunoqIijR071uvPOBwOORwO376Mj9qu5+/rRjL6xB394clofbL6ixrNfW2Tx55/NY12zSkq11M/HqvLRyfLfDBeGWcNkFxOGaRrOmW03xNfGK1PjPI9OxMRATA7O1tffvllj342ISFBLS0tbsdaWlrUt29fr+1nz57t9vmHP/yh3n33Xa1evbrDALht27Ye1eaLioqKgF07UtEn7ugPT0boE4fLpftX/tvrhs9tx+5/+3N9a1KyYk0mQ/RJd9EnnugT44iIAHgqRo4cqV27dun48ePq1av16+7YsUMjR4702v6FF16Q1WpVTk5O+7GjR48qLi6uw3ukpaUpISHBr3U7HA5VVFQoIyNDsbGxfr12pKJP3NEfnozUJ5/9q171Ld4fjWlT3+LUkf6nK+HQXkP0ia+M9HviK6P1SXNzc0AHbyJB1AfA7OxsDRw4UEuXLtUvf/lLVVdX69VXX9Udd9zhtb3NZtMbb7yh5cuXKyUlRX/5y1+0ceNGFRYWdniP2NjYgP2FCeS1IxV94o7+8GSEPvnm8DGf250pY/RJd9EnnozSJ0b4jl2J+gDYq1cvvfjii/r1r3+tiRMnKiEhQTNmzNA111wjSdq3b58mTZqk5cuXa8KECbrrrrsUExOjH//4x2pqampfMHLWWWeF+JsAwH/5uuHz4P5x0sEAFwMg4kRdAPzwww89jp111ll64YUXvLYfOnSoNm7c2P65T58+uvfee3XvvfcGrEYAOFVtGz7XNNq9PgdoUuu2L+efnaSKz/cEuzwAYS4q9gEEAKNhw2cAp4IACAARKi89hQ2fAfRI1E0BA4CR5KWnKNdqUVl1g+qa7Brcv/U9v4z8AegMARAAwpTD6fIp2MXGmJQzfFAIKgQQqQiAABCGOnvFG1O7AE4VzwACQJgpqbRpTlG5W/iT/vuKt5JKW4gqAxAtCIAAEEYcTpcKi6s6fcVbYXGVHE5vLQDANwRAAAgjZdUNHiN/J3JJsjXaVVbdELyiAEQdAiAAhJG6po7DX0/aAYA3BEAACCO+v+LNt3YA4A0BEADCSNsr3jraxc+k1tXAWalJwSwLQJQhAAJAGOEVbwCCgQAIAEHmcLpUurNef920V6U76z1W9PKKNwCBxkbQABBEvm7wzCveAAQSARAAgqRtg+eTd/Br2+D55NE9XvEGIFCYAgaAIGCDZwDhhAAIAEHABs8AwgkBEACCgA2eAYQTAiAABAEbPAMIJywCAQA/cThdHa7abdvguabR7vU5QJNat3lhg2cAwUAABAA/6Gp7l7YNnucUlcskuYVANngGEGxMAQPAKWrb3uXkRR5t27uUVNokscEzgPDBCCAAnIKutncxqXV7l1yrRbExJjZ4BhAWCIAAcAq6s71L26bObPAMINQIgADQhc4Wd7C9C4BIRAAEgE50tbiD7V0ARCIWgQBAB3xZ3NG2vUtHT/CZ1BoY2d4FQDghAAKAF76+u1dq3b5FkkcIZHsXAOGKAAgAXnRncQfbuwCINDwDCMCw/Lm4g+1dAEQSAiAAQwrE4g62dwEQKZgCBmA4LO4AYHQEQABRx+F0qXRnvf66aa9Kd9bL4XS5nWNxBwCjYwoYQFT57Gu75q7+SDUHj7QfO3FqtyeLO06eKraccD0AiEQEQABRY/UXNVpSesDjeNvU7rLp43TkuNOna7G4A0A0IwACiBidrdp1OF369btbvP6cS63TtoXFVXr02vN8uheLOwBEMwIggIjQ1ardsuoGt2nfk7VN7crU+nM1jXavzwGa1DrFy+IOANGMRSAAQqqzBRttfFm16+u+fd8cOsLiDgCGxwgggJDpalRP6nrVbk+mdnOGD2JxBwBDIwACCIjOnteT/juqd3KwO3HBRndW7cokWU6L63Aa+OSpXRZ3ADAyAiAAv+tqZM/XUb1cq6VbU7sPXDlGt762SSbJ7dodTe2yuAOAUfEMIIBu6eqZPV+e1+vOXnzdeSXbFedYtCBngIacFud2zpJobh9RBAAwAgjgP7qaspX8N7J31xWjfKqprsmuK88d6vuqXZdTF5xu1k0/yNaG3Y1M7QJABwiAQJTzR7Bra9PVM3uJ8X18GtlrOHzUp9oH9zcrNsakgnyr5hSVdzm163C0fmZqFwA6xxQwEMVKKm268Lcf6kfLP9Ptr2/Sj5Z/pgt/+6FKKm1ubbqasvX1/bk1jS0+1ZXUL04piWaPbVjamNQaQE9csLFs+jhZEt2ng5naBYCeYQQQiFAOp0uVdUe0e/M+WRITerTKNtdq8WnKtn9cb7+O7FlOM/s8qteGVbsA4D9RNwLY0tKiqVOn6q233vI499VXXyk7O1tff/11p9d4++23lZubq7Fjx+qaa67Rxo0bA1Uu4MHXjZG/u+QjFXy8X3f86XOPkT1fR+w+21nvU7Ar/dc3PtXenZG9nozqtU3tTh47TDnDBxH+AKCHomoEcPv27br77rv1xRdfaOrUqW7n/va3v+n+++/XgQMHOr3G2rVr9eCDD2r58uU699xztWLFCs2ZM0dr1qxRfHx8AKtHtAvHZ/F8DXae78zwrrsje4zqAUBoRM0IYGlpqWbNmqWrr75aQ4cOdTv35JNP6ve//73uuOOOLq/zxhtvaNKkSRo/frx69+6tn/zkJxo4cKBWrVoVqNIR4XwdsQvHZ/F8DXY5wwcFbGSPUT0ACL6IGQG02+2qra31ei45OVmjR4/WmjVrFBcXp5deesnt/HXXXae5c+dq7969Xd5nx44dmjJlituxESNGaOvWrR3+jMPhkKNt+aGftF3P39eNZMHuE4fTpXW7GlTXdESD+8fp/LM9R6ZWf1GjX7+7xe3tE5bT4vTAlWN0xTmW9jZzX9vU4YjdUz8eq8vHDNGid77o8lm8vr1jfBrZ++aQb5snZ6cO1Jsb4lR78EinW6ycf9YA3T9ptOZ2ssny/ZNGSy6nHA4pd8xgXTYq2Wv/BfKfH39vPNEnnugTT0brE6N8z85ETADcvHmzZs6c6fXcU089pcsvv7zDnx0yZIjP9zl8+LDHVK/ZbFZzc3OHP7Nt2zafr99dFRUVAbt2pDrVPnG4XNry76Pab3dqoDlGY5L7KNbkHuw++9quFzcdVH2Ls/3YoPgY/XTsabrgdHN7myWlBzyuX3PwiG59bZMW5AzQ+cPidP/Kf3c6Ynf/25+r9uvTOnyFWVtbW6Ndf/lsi0/fsekbmwbFx7jVf7JB8TGKO7hbM86J15JS7/d2SZpuNavi880aIunOnAEe/ZL0n34ZcqxGmzbVuP28WdKZknRQqvh8j0+1+wN/bzzRJ57oE0/0iXFETADMzs7Wl19+GfD7xMfHy253Hz2x2+0aOHBghz+TlpamhIQEv9bhcDhUUVGhjIwMxcbG+vXa4cCX0bWT29Y2tqjp33t1/aXj1ae356+uP0fsHi31HLFraHHq0dID7SN2c1d/1OH3M0kqqrLr3NEjVN/ifeS6TX2LU9/EJEk60Gk7SbIMGSJt+VeX7bIzRmn0iGOa+9omSd5H7B68+lyNP8ei8ZlSaqpnv6QkmnX/pNHt/SJJY8dKN/3A9392wRbtf296gj7xRJ94MlqfNDc3B3TwJhJETAAMlpEjR2r79u1ux3bs2KHvfve7Hf5MbGysX//COJwurfuqXut2t8h+2gFdMDy5/V+wviwkCPd2vix06Kztso1/16IfnuPW1tfFE96mYmsPHtHc1za1b4vy4MqtnU7FPrhyqxLj43wasVu7a3+HbU5kMvn2OO53RiTrrY37unwrRtvvzLKYGI9+sXjp6x+cO0xXpA/16Z9xbKw0ceRgn+oNFX//nYwG9Ikn+sSTUfrECN+xKwTAk1x77bWaO3euvv/972v8+PFasWKF6uvrlZubG5T7ewSZtevag4wkn4KTrwErFO18WcHaVdvag0fc2oZivzt/r57NGT5Ify7/uutg9+1BAVtly9szAMA4omYVcE/t27dPmZmZWr9+vSQpJydHBQUFWrRokbKysrRy5UotX75cAwYMCHgtna0CvaWoXLd0sUK0q2uEup2vK1gdTpfPbY8ed4Zkvzt/r55tC3bernxysGtbZTvktDi3dqyyBQD4KipHAD/88EOvx08//XSP5wiHDh3qsdHz5MmTNXny5IDV540vgcebE0evLhs9xKdRrlC183V0ray6QfrPf++q7aulu6JmxK4t2PkyZZuXnqLLRiXrDx+U6bQhZ3h9EwgAAB2JygAYicqqGzoNMp3pbhgKVTtfQ1hdk+/98FVDx6uz3UVOsOvOlG364DiNPW8oz7MAALqFABgmuhN6OuJrGApVO19D2OD+5q4b/cdZSb6tvo60YMezeACAQCIAhonuhJ6O+BqGQtXO1xCWlZokqfXZuK7azsg5W8//o5pgBwBANxh+EUi4yEpN6nSxQGfaFhLMyDnbpwUHoWrXnYUOsTEmn9r26RXT7cUTvryiLC89Rf+4+zL94aYL9PgNY/WHmy7QP+6+zGOBhcQiCwBA5CEAhglfAk9n57oThkLVrichzJe2BDsAALrH5HK5Oltkik40Nzdry5YtGjNmjN/eBNLZXnpS5O8D2MbXjaVPbFvT2KyDtXv0o8uzOnwTiK/XjAYOh0ObNm3S2LFjWQTyH/SJJ/rEE33iyWh9Eoh/f0cangEMM23PlH22899aV7lN56enub0JxJfnzXx9Li1U7aTuPQ/X1tbhGKBNm+o6DHU8YwcAgG8IgGEoNsakC749SOaD8Rr7bfepR19DTri3AwAAocMzgAAAAAZDAAQAADAYAiAAAIDBEAABAAAMhgAIAABgMARAAAAAgyEAAgAAGAwBEAAAwGAIgAAAAAbDm0BOgdPplCS1tLT4/doOh0NS6/sKjfBeRl/QJ+7oD0/0iSf6xBN94slofdL27+22f48bkcnlcrlCXUSkqq+v165du0JdBgAA6IGzzz5bgwYZ8/WlBMBTcPz4cTU2NiouLk4xMcymAwAQCZxOp44cOaLExET16mXMyVACIAAAgMEwbAUAAGAwBEAAAACDIQBGgJaWFk2dOlVvvfVWqEsJifr6et16662aMGGCsrOz9dBDD+n48eOhLissNDQ0KDc3V2vXrg11KSG3detW3XjjjcrKytLEiRN11113qaGhIdRlhVRpaamuu+46jRs3ThMnTtSDDz4ou90e6rLCgsPh0IwZM/SrX/0q1KWE3KpVq2S1WpWZmdn+Z8GCBaEuCwFGAAxz27dv17Rp07Rp06ZQlxIyv/zlL5WQkKC///3vevPNN1VaWqqXX3451GWF3IYNGzR16lTt3r071KWEnN1u1+zZs5WZmal//OMfevfdd3XgwAHde++9oS4tZBoaGvTzn/9cP/rRj7R+/Xq9/fbbKisr03PPPRfq0sLCk08+qfXr14e6jLBQUVGhyZMna+PGje1/lixZEuqyEGAEwDBWWlqqWbNm6eqrr9bQoUNDXU5IfPXVVyorK9OCBQsUHx+vM844Q7feeqtWrFgR6tJC6u2339add96pO+64I9SlhIV9+/Zp9OjRmjt3rvr06aOBAwdq6tSpWrduXahLC5mkpCT985//1DXXXCOTyaQDBw7oyJEjSkpKCnVpIVdaWqr3339f3/ve90JdSlioqKhQenp6qMtAkBlz7XOYsNvtqq2t9XouOTlZo0eP1po1axQXF6eXXnopyNWFh+3bt2vAgAEaMmRI+7Hhw4dr3759OnjwoE477bQQVhc6F154ofLz89WrVy9CoKRvf/vbev75592OrV69Wuecc06IKgoP/fr1kyRdfPHFqq2t1YQJE3TNNdeEuKrQqq+v13333aenn36amQS1bofyxRdfKD4+Xs8//7wcDocuvvhi3XnnnUpMTAx1eQggRgBDaPPmzfre977n9c8///lPDRw4UHFxcaEuM6QOHz6s+Ph4t2Ntn5ubm0NRUlhITk427N5VXXG5XPrf//1frVmzRvfdd1+oywkL77//vj755BPFxMTotttuC3U5IeN0OrVgwQLdeOONGj16dKjLCQsNDQ2yWq264oortGrVKr3++uvatWsXzwAaAP8GCaHs7Gx9+eWXoS4jrCUkJHi8aq/tc9++fUNREsLYoUOHdM899+iLL75QUVGRRo0aFeqSwoLZbJbZbNaCBQt03XXXqbGx0ZCjO88++6z69OmjGTNmhLqUsPGtb33L7ZGa+Ph4LViwQNdff70OHTrUPoqM6MMIIMLayJEjdeDAAX3zzTftx3bu3CmLxaL+/fuHsDKEm927d2vKlCk6dOiQ3nzzTcOHv/LycuXl5eno0aPtx44eParevXt7jKobxV//+leVlZVpwoQJmjBhgt599129++67mjBhQqhLC5mtW7fq0Ucf1YnvhDh69KhiYmLUp0+fEFaGQCMAIqydffbZGj9+vB5++GEdOnRIe/bs0dNPP61rr7021KUhjDQ2NmrWrFkaN26cXnjhBRY6SBo1apTsdruWLl2qo0ePau/evfrtb3+ra6+91rD/Yi8pKVF5ebnWr1+v9evX68orr9SVV15p6NXAAwYM0IoVK/T888/r+PHj2rdvn5YsWaKrr77asL8nRkEARNh74okndPz4cf3P//yPrr/+el100UW69dZbQ10Wwshbb72lffv26b333tP48ePd9jMzqr59++r555/X9u3bNXHiRM2YMUPf+c53DL01DjxZLBY9++yz+tvf/qasrCxNmTJFGRkZeuCBB0JdGgKMdwEDAAAYDCOAAAAABkMABAAAMBgCIAAAgMEQAAEAAAyGAAgAAGAwBEAAAACDIQACAAAYDAEQAADAYAiAACLG119/rVGjRunrr7/2+7Wbm5v1s5/9TOedd56mTZvmcf5Xv/qVzjnnHGVmZqq+vl7/93//pxkzZni91qhRo7R27dpO7/fSSy8pMzPTp7YA4G8EQACQtGXLFv3jH//Qxx9/rBUrVnhtk5+fr40bN2rQoEGnfL8bb7xRGzduPOXrAEBP9Ap1AQDQU3v37tWSJUu0du1axcTE6IILLtDdd9+twYMHS5L++c9/6ne/+512796ttLQ0jR8/Xp9//rleffVVt+t88MEHmj9/viTp0ksv1b333qvrrrvulOurra1VXl6e2zG73a7MzEy99tprp3x9AOgpAiCAiHTs2DH99Kc/VXp6ut5//325XC4VFhbqlltu0Z/+9CfV1NTolltu0X333acpU6Zo06ZNuuWWWzRmzBiPa11++eVavny5Zs6c2a1RuQ0bNmjChAkdnh8yZIjb9T766CPNnz9fv/rVr7r3ZQHAzwiAACLS+vXrtWfPHv35z39Wv379JEmFhYXKyspSZWWlSktLNWbMGE2dOlWSNGHCBF1//fWqqKjwWw3jx4/3GE2UWp8BPNnnn3+u+fPn63e/+53OPfdcv9UAAD3BM4AAIlJ9fb0GDhzYHv4kqV+/fhowYID27t0rm82mYcOGuf3MGWecEewyJUlfffWVfv7zn+v222/X5ZdfHpIaAOBEBEAAEWnYsGHav3+/Dh061H6sqalJ+/fvV3JysoYNG6Z9+/a5/czJn4OhoaFBs2fP1qRJkzRr1qyg3x8AvCEAAohIGRkZGjFihAoKCtTU1KSmpiYtWrRIZ555psaNG6fJkydry5Yt+stf/iKHw6HNmzfrT3/6U1BrbG5u1s0336wxY8bo3nvvDeq9AaAzBEAAEalXr1569tlndfz4cV1xxRW69NJLdezYMb300kvq1auXLBaLnnjiCS1fvlwTJkzQb3/7W1144YXq3bt30Gp8//33VVFRobVr1+r8889XZmZm+x8ACCUWgQCIGKeffrq+/PLL9s8pKSl6/PHHvba12WxKTk7WypUr24898sgjHV47Ozvb7dpdmTdvXofnTrzOVVdd5fM1ASBYGAEEEJX279+vH//4x6qsrJQkbd26Ve+8844uvfTSEFcGAKFHAAQQlaxWq+677z7Nnz9fmZmZmjdvnm6++WZdeeWVPb5mcXFx+6vgTlXbq+AAIBRMLpfLFeoiAAAAEDyMAAIAABgMARAAAMBgCIAAAAAGQwAEAAAwGAIgAACAwRAAAQAADIYACAAAYDAEQAAAAIMhAAIAABgMARAAAMBgCIAAAAAGQwAEAAAwGAIgAACAwRAAAQAADIYACAAAYDAEQAAAAIMhAAIAABjM/weDQIG0yzRYzQAAAABJRU5ErkJggg==" } }, "cell_type": "markdown", "id": "8450c57f", "metadata": {}, "source": [ "![ROI.png](attachment:ROI.png)" ] }, { "cell_type": "code", "execution_count": 7, "id": "38443d0d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "enter the fc value:3\n", "enter the DC conductivity value:-10\n", "enter the exponent:0.8\n", "dumped_parameters {'fc': 1000.0, 'DC': -10.0, 'n': 0.8}\n" ] }, { "data": { "text/plain": [ "()" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#dump the initial guess parameters using dump parameters method (varies for each fn), which dumps the parameters in a json file'\n", "#this is required before performing the first fitting as it takes the initial guess from the json file created\n", "cond.dump_parameters()\n" ] }, { "cell_type": "code", "execution_count": 8, "id": "c1758bae", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "fit parameters:\n", " [274.40483973 -11.04255983 0.97224746]\n", " DC cond = -11.043\n", "fit parameters dumped for next iteration {'fc': 274.4048397299556, 'DC': -11.04255983252695, 'n': 0.9722474585471939}\n" ] }, { "data": { "text/plain": [ "()" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# perform least squares fitting\n", "# the plot is shown as a image in the next cell\n", "\n", "cond.fit(x1,y1)\n" ] }, { "attachments": { "final_fitting.png": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAHgCAYAAAA10dzkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAABfd0lEQVR4nO3deXhU5f3+8feZyWQlEAIBsiB7ApRoEAQR96/UVrS0gqJScKl1Q6wiUhcI2BStRaz7Am610mqDtm64/FRqa4sICDQgIUBZk5AQAkkg22Tm/P4YJslkkhBCJpPM3K/ryhXOOp95SJib55znOYZpmiYiIiIiEjQs/i5ARERERNqXAqCIiIhIkFEAFBEREQkyCoAiIiIiQUYBUERERCTIKACKiIiIBBkFQBEREZEgowAoIiIiEmQUAEVERESCjAKgiIiISJBRABQREREJMgqAIiIiIkFGAVBEREQkyCgAioiIiAQZBUARERGRIKMAKCIiIhJkFABFREREgowCoIiIiEiQUQAUERERCTIKgCIiIiJBRgFQREREJMgoAIqIiIgEGQVAERERkSCjACgiIiISZBQARURERIKMAqCIiIhIkFEAFBEREQkyCoAiIiIiQUYBUERERCTIKACKiIiIBBkFQBEREZEgowAoIiIiEmQUAEVERESCjAKgiIiISJBRABQREREJMgqAIiIiIkFGAVBEREQkyCgAioiIiAQZBUARERGRIKMAKCIiIhJkFABFREREgowCoIiIiEiQUQAUERERCTIKgCIiIiJBRgFQREREJMgoAIqIiIgEGQVAERERkSCjACgiIiISZBQARURERIKMAqCIiIhIkFEAFBEREQkyCoAiIiIiQUYBUERERCTIKACKiIiIBBkFQBEREZEgowAoIiIiEmRC/F1AZ1ZTU0NJSQlhYWFYLMrSIiIinYHT6aSqqopu3boREhKcUSg433UbKSkpYffu3f4uQ0RERFqhf//+9OjRw99l+IUC4CkICwsDXD9AERERbXpuh8NBTk4OycnJWK3WNj13Z6U28aT28KY28aY28aY28RZsbVJRUcHu3btrP8eDkQLgKXBf9o2IiCAyMrJNz+1wOACIjIwMil/GllCbeFJ7eFObeFObeFObeAvWNgnm27eC952LiIiIBCkFQBEREZEgowAoIiIiEmQUAEVERESCjAKgiIiISJBRABQREREJMgEfAPft28cvf/lLxowZw7hx45g7dy6lpaVN7r9p0yauuuoqRo4cycUXX0xmZmY7VisiIhJcMjNh2DBISjq5r2HDYMUKf1ffeQX8PICzZ89m9OjRPP/88xw7dow777yT3/3udzzyyCNe+5aUlHDLLbdw1113MXXqVNauXcvMmTNJSUnh9NNP92mdmZmQng5lZe41Fuz2VGy2E2f06GjIyIApU3xaooiISJtLT4fs7NYdO3++PvtaK+AD4M6dOxk1ahSmaWKaJoZhNPnUjs8++4yYmBimTZsGwLhx47jiiitYvny5zwOg9y+AAYS2+Hj9EoiISGfk7viwWCA+vmXH5OeD01m/00ROVqcPgJWVlRQUFDS6LS4ujlmzZrFkyRLeeOMNHA4HaWlpzJkzp9H9t2/fTnJysse6wYMHs+IEfcwOh6N2FvXWKiuzAAYWi1n7C2C327HZbM0e5/olMCgrM3E4nKdUQ1vYvXs3/fv398m53W18qm0dKNQe3tQm3tQm3tQm3vzbJq7Pv/h4kz17WvY51q+fhdxcA2jdZ5/+7gMgAG7atIkZM2Y0uu25557DMAxuv/12brzxRg4fPszs2bNJT09n8eLFXvsfO3bMq3cwPDyc8vLyZmvIyclp/Rs4zm5PBULp2dPOe+9ltfi4yy5LpbAwFLvdzsaNLT/uZKxYsYKtW7cyf/78Zvdbv349b7zxBk899ZRP6nDLyvLN++ys1B7e1Cbe1Cbe1Cbe/NEm7s+/k/kca80x4qnTB8CxY8eybdu2Rrdt3ryZ++67j7Vr1xISEkJkZCRz585l2rRpLFiwgC5dunjsHxERQVmD/uTKykqioqKarSE5OfmUnwXsvtfPZrORlpaGw+EgKyuL1NTUZp/L2PA4X/j666/Zt2/fCc+/a9cun9bR0jYJFmoPb2oTb2oTb2oTb/5sk9Z8jp3qZ195eXmbdN50Zp0+ADYnPz8fh8OB01nXPWyz2TAMo9Ef8OTkZP797397rNuxYwdDhgxp9nWsVmsb/sJ41tbyczf+nlrju+++Y9GiRezcuZOhQ4fSr18/DMPAYrGwbNkyPvjgA/Lz8zEMg/PPP59FixaxadMmHn74Yex2O6NHj+aTTz4B4NFHH+W///0vhw4domfPntx+++1MOcWbFdu2vTs/tYc3tYk3tYk3tYk3/7aJ5+dYxc4KIgY1fs9+U8e0lP7eA3wamFGjRhEREcEjjzxCVVUVhw4dYsmSJUyYMKHRgSATJkygqKiI119/HbvdzjfffMMHH3zA5MmT/VC9fxw+fJhbb72VSy+9lLVr13Lffffx+eefA/Dxxx/zxhtv8Mwzz7Bu3Treeustvv76az744APGjh3Lww8/TEJCAhs2bKB3797MmzcPm83GRx99xHfffcfPf/5zMjIyOHbsmJ/fpYiIdGTF/6+YNUPWsH3WdmrKavxdTkAK6AAYGxvLK6+8wu7duznvvPP46U9/Sv/+/T2mgJk4cSIvvvgiAN27d+fVV1/lk08+YezYscybN4958+Zx9tln++sttLt//OMfRERE8Mtf/hKbzcaoUaNqA/D555/PihUr6N+/P8XFxRw+fJiYmJgmB+H89re/ZcGCBdhsNvLy8oiKiqKyspKSkpL2fEsiItKJOI45yLk1B0zIfTaXgysO+rukgBTQl4ABRowYweuvv97k9o8++shjOTU1lbfeesvHVXVcBQUFxMfHYxhG7brTTjuNrVu3Ypomf/jDH1i1ahWxsbEMGzYMu92OaZqNnmvfvn38/ve/rx0Z3K9fPwCPS/IiIiL17Vqwi8pdlQB0O78bfa7v4+eKAlPAB0A5OX369CE3Nxen04nF4uogPnDgAACPP/44eXl5fPnll7UDaK644opGz2O327n11luZPXs21113HYZhsHnzZt5///32eSMiItLplK4rZf8f9gNghBmkLEvBsBgnOEpaI6AvAcvJu/jiizFNk2eeeYbq6mo2b95c+zi8o0ePEhYWhtVqpaqqildffZWcnBzsdjsAYWFhVFRUUFNTg91up7KykvDwcAzDIC8vr3bqHff+IiIiblbTybZfbIPjF4n6L+hPZPKpzbAhTVMAFA9du3bllVdeYfXq1YwZM4aHHnqISy+9FIC7776byspKzjnnHC6++GI2btzIpEmTaofSn3XWWfTo0YOzzjqLffv28cgjj/Dcc88xcuRIZsyYwfjx4+nZs2fQD70XERFvlx3dx7H/ugYJRp0RRd85ff1cUWDTJWDxMnTo0Cbvg3zzzTebPK5Xr158+OGHtcspKSlMmjTJY59bbrmlbYoUEZGAkUQ5Pyvd7VqwQMrLKVhs6qPyJQVAERER8RvDNJnDNkJxDSj8KDKJn/+0a7PH5Oe3R2WBTQGwg8nPh6QkAAt2e2rtbOfN7S8iItJZ9Y6qwXb8xr9cwnnm6ACqjrbs2OhoHxYW4BQAOwj3D7HTCbm5AAYQetLHi4iIdCb3P2Ijfd6ZnH0gl9yQKHqGt+wpHdHRkJHh4+ICmAJgB5GRAfPnQ92jiE3sdjs2mw1XGGyafglERKQjycyE9PT6n2knYvBZlySg7jPtFJ8aKiegANhBTJni+cPucDjZuDGLtLQ0PbNQREQ6lfR0yM5u/fHz5ysA+poCoIiIiLQpd8+fxQLx8d7bI5127ireworoAewI61a7Pj/fdStUy3sOpbUUAEVERMQn4uNh/37v9dk37+TAK0cYUb2BlGUpxP/ClRKTktz3wYuvaZIdERERaTeHvzzMgVdcjxi1drHS/dLufq4oOCkASos5HA727dvn7zJERKSTcpQ72HbLttrlgb8fSHhSuB8rCl4KgOJl5cqVjBs3jpSUFIYPH05eXh4A99xzD3//+9/9W5yIiHRauxfupnJnJQDdzu1Gwi0Jfq4oeCkAipfMzEwmTpzItm3b+P7770lIcP2CHj582M+ViYhIZ1X2XRn7lriuIhmhBsnLkjEszU9zJr6jQSDiYcqUKWzZsoW1a9eyfPlynE4nX3zxBS+88ALr1q1jw4YNbNmyhRdffNHfpYqISCfhtDvZ9ottHH/gB/3T+xM1NMq/RQU59QCKhxUrVjB69GhuvfVW/t//+3+16xctWlS7XuFPREROxv4n9nN0o+v5blGpUfS9r6+fKxL1APrBvif2se+JEw+mcAx0wCrPdVk/yaLsuxNPkNR3dl/6ztYvmIiI+FdNSQ17HtnjWrBAysspWELV/+RvCoB+UFNaQ3Vu9Qn3s8R6/4LYD9pbdGxNaU2rahMREWlLId1CGPn1SLbdvI1u47vRdUxXf5ckKAD6RUjXEEITQ0+4nyPG4bXOFmdr0bEhXfVXKyIiHUOX1C6c+Z8zcdqd/i5FjlNK8IOWXJ51OBxs3LjRa33q+6k+qurEQkNDKdPzeUREpBUMq6Fn23cgCoDSYj/96U9ZuHAhmzdv5s9//rO/yxERkQ7sIgr5d15PkpJafr9ffr4PCxIPCoDi5U9/+lPtn7dtq5ux/YorruCKK67wR0kiItKJnMdBbuV7dpmRLM4dylZO7r6/6GgfFSa1FABFRESkzdiP2Lm5fDsAAyhnWPdKSiNbHgCjoyEjw1fViZsCoIiIiLRIZiakp0Nzt4PfdPh/XFzumq3iu/AerImIA+qC3ZQp7VGpnIgCoIiIiLRIejpkZze9/QwOczGuG/mOYeXRyiEU5dU97m3+fAXAjkIBUERERFrE3fNnsUB8vOc2m+ng/oIcOD6DWWa3gYR1CScR1+AOp7P5nkNpXwqAIiIiclLi42H/fs91O+/fzb7HKgDoOr4rr/0zAeP4AOCkJMjNbecipVl6FouIiIickrINZex73PWIUyPUIGVZCobFOMFR4k8KgCIiItJqzhon227eVnvpt9+8fkQNi/JvUXJCCoAiIiLSek7oMbEHRohB1IgoTvv1af6uSFpA9wCKiIhIq1lCLQz4zQDirorDdJhYQtW31BkoAIqIiMgp65Laxd8lyElQTBcREZGT5jjm8HcJcgoUAEVEROSkxDiq+GbgN+x+eDfOaqe/y5FW0CVgEREROSkzjmzHXmln98LdOCudDHx0oL9LkpOkHkARERFpsXM5yJjKIgBscTb6zunr54qkNRQARUREpEUinXZ+xfba5cFPDcbWw+bHiqS1dAlYREREWuSakv/Rk2oANoTF8vM5veC+Ex+Xn+/jwuSkKQCKiIjICR3+x2EuLncluXKsPFqVzMG8k3vcW3S0LyqT1lAAFBERCXIrVsD99w/Hbm/8zjCb6eDRghz6HF9+2RjIISMcC2CxQNeuEBHR/GtER0NGRpuWLadAAVBERCTILVxoYffuphPczeyhDxUAbKYrfzcTME3XNqcTevWCrVvbo1JpKwqAIiIiQa6szPXdYjGJj/e8rGs1nZxfeAhqwI7BG71SSLC59snPdwVA9/HSeSgAioiICADx8bB/f8O1FhyVZ7L3kb1YIi18e39U7ZakJMjNbdcSpY0oAIqIiEizrOFWBvxmgL/LkDYU8PMA7tu3j1/+8peMGTOGcePGMXfuXEpLS5vcf8GCBYwYMYKRI0fWfr399tvtWLGIiIiIbwV8AJw9ezaDBw/m3//+Nx9//DF5eXn87ne/a3L/rKwsMjIy2LBhQ+3X1KlT27FiERER/zJNkx337uDY98f8XYr4SMAHwJ07d2KaZu2XYRhENDFWvbq6mpycHEaMGNHOVYqIiHQcB149wP4n9rNu5Dpyn9dNfoGo098DWFlZSUFBQaPb4uLimDVrFkuWLOGNN97A4XCQlpbGnDlzGt0/Ozubmpoann76adavX090dDSTJ0/m5ptvxmIJ+KwsIiJCVX4VO+7dAYBZbRIx6AQT/Emn1OkD4KZNm5gxY0aj25577jkMw+D222/nxhtv5PDhw8yePZv09HQWL17stX9ZWRljxoxh+vTpPPHEE2zdupWZM2disVi4+eabm6zB4XDgcDja7D25z1n/u6hNGlJ7eFObeFObeFObNKaukyPnzhwcJa626fXzXnS7pFszbWUBDMDE4XD6vMq2or97MEzTPZVj4Nm8eTPTp09n7dq1hIS4su769euZNm0a69ato0uXLic8x8svv8zKlSt59913vbaVl5ezVTNfiohIJ3fZZakUFoYysVs+c0q2AWB0N4jKjMKIafpxb+7jevWqZuXKrPYqt80MGzaMyMhIf5fhF52+B7A5+fn5OBwOnM66/5XYbDYMw8BqtXrt//nnn1NUVMQ111xTu666uprw8PBmXyc5ObnNf4AcDgdZWVmkpqY2WmswUpt4Unt4U5t4U5t4C9Y2WbHC9cSPxiZtLiqCKOzcULKrdt0zDOWb63sBrse4Pfywk8mTPY+z2SzHv9tIS0vzVeltrry8nJycHH+X4VcBHQBHjRpFREQEjzzyCA888ABHjx5lyZIlTJgwodGBIKZp8uijj9KvXz/OPvtsNm7cyBtvvMEDDzzQ7OtYrVaf/SPiy3N3VmoTT2oPb2oTb2oTb8HWJgsXQnZ209tv43/0pBqA1cTyzuFecLiu92/BAitXX93U0Y13rHRUnalWXwnoABgbG8srr7zC448/znnnnUdYWBgXX3wx9913X+0+EydO5IorruC2225jwoQJPPDAAyxcuJCCggJ69uzJrFmzmDRpkh/fhYiIyKmre9yb64kf9fXOK+ZyMx+ACsPKX3olkxiix70FsoAOgAAjRozg9ddfb3L7Rx995LF8zTXXeFwCFhERCSQNH/fmrHHydvh2OD4u4u2uAykO8b71KT/f9ei3huukcwr4ACgiIiJNs4RYWJ6YwtV7t1GCjTdKEjBLvPdzOpt+7m90tG9rlLanACgiIhKgMjMhPd11+dbdW9dYT15FRTe+sIwm2rRjGAYGrkvFXbvCkSOu8NfYpWNwhb+MDF+/E2lrCoAiIiIBKj3de+BH4z15BmClEiuYdfv16gUREa79G146ls5NAVBERCRA1R/4AZ49ef2ry9hj64JpuCZyttvt2Gw2wNDAjyCgACgiIhLg3Jdu3T15OavKWXf6Brqc2YWUl1MITw5n48Ys0tLSsFqtJCU1fb+fBAY94FZERCSIGKbJtl9uw1nppPQ/peS/rKG8wUgBUEREJIhcWJ5PyVeuYb7hA8IZ8JsBfq5I/EEBUEREJEj0pIprS3bWLicvTcYapadiBCMFQBERkWBgmvyKHCJN14zPfW7oQ+wlsX4uSvxFAVBERCQInFV5kHM5BICtl41BSwb5uSLxJ40CFhERCXBRTjvXH9leu/yEfQhrTrfV28OC3Z6KzebqF6o/abQEJgVAERGRAHddyU5inHYA/k0P3j0cB4fr72EAoV7HOZ11f9bj3gKLAqCIiEgAMzCpNlw9exUWK2/1SibRajTYq/GJoN2TRutxb4FHAVBERCSAmRj8MSaZ+z7pRXV+NVlXhXnt43A4G50IWo9/C1wKgCIiIkEg5twYf5cgHYgCoIiISAAyTRPXvX2uwRxJSc3t3fggEAlcCoAiIiKdTGYmpKdDWVnj262mkweLNjK2pjd/IwGn0zjBs30bHwSigR+BSwFQRESkk0lPh+zsprdPYx8plJJCKf05yhOkeGwPCYHeveuv8RwEAhr4EegUAEVERDoZd8+fe5RuffH2Y1xfuBsAJ7AmLoHE45177tG9vXt7Du5oOAhEAp8CoIiISCfVcJSu6TTZcP42SgtNAPrd15evf193Hdc9uldEj4ITEREJEHkv5FH671IAwgeF039hf/8WJB2WAqCIiEgAqNxbyf/u/1/tcsqyFKyRupwrjVMAFBER6eRM0yTn9hwcRx0AxN8cT/eLuvu5KunIFABFREQ6ucK/FFK8shiA0PhQBi4e6OeKpKNTABQREenETNNk3+J9tctDnh+CLcbmx4qkM9AoYBERkU7MMAzS/pHG/+7/H/ZiO3E/jfN3SdIJKACKiIh0ciHdQkh+IRnTYfq7FOkkdAlYREQkQBhWw98lSCehACgiItIJXUYesTWV/i5DOildAhYREelkhlSVcC85VBVayXt5EAk3J/i7JOlkFABFREQ6iMxMSE+ve9ZvY2ymg4VF27AAEaaDu29x8M6trucCd+0KERFNH5uf3+YlSyelACgiItJBpKdDdnbz+9zIXvpRDsBWonnHTMJpgtMJxcUte53o6BPvI4FNAVBERKSDcPf8WSwQH++9va/9KNcV7gWgBoMlRgoYBjibP66+6GjIyGjDoqVTUgAUERHpYOLjYf9+z3XOGicbxm2jrNA11cug+aex4zddAEhKgtzcxo8TaYxGAYuIiHQCuU/lUrbO1UUYOSySfg/183NF0pkpAIqIiHRwFf+rYNf8Xa4FA1JeTsESpo9waT399IiIiHRgpmmy7ZZtOCtcN/olzkyk2znd/FyVdHa6B1BERKQDMwyDpLuSKN9ajhFiMOCRAf4uSQKAAqCIiEgH1/MnPYm5IIbK3ZWEROujW06dLgGLiIh0AiHdQuhyRhd/lyEBQgFQRESkA6rKrcI0TX+XIQFKAVBERKSDiXLaWTdqHZt/spnK/ZX+LkcCkG4kEBER6WCuK9mJvdzOoQ8PYY22MvzPw/1dkgQY9QCKiIh0IKMo5oLyAwBYu1oZtHiQnyuSQKQeQBERER/LzIT09Lpn/TblcJ6DxeTULi8uG8TH/cPo2hUiIpo+Lj+/jQqVoBHwAbCgoIBFixaxZs0abDYbEydOZPbs2YSFhTW6/6ZNm/jtb3/Ljh076N69O7fffjtXXXVVO1ctIiKBJD0dsrNPvN8d7CIe1z1/G+nGh2Y8Zg0UF7fsdaKjT6FICSoBHQCdTid33HEHcXFxfPLJJ5imyT333MPChQt59NFHvfYvKSnhlltu4a677mLq1KmsXbuWmTNnkpKSwumnn+6HdyAiIoHA3fNnsUB8fOP7DKouYfLB/QBUYeEJIwXTNE54nFt0NGRktFXFEugCOgDu2rWLzZs3s2rVKrp37w7A7NmzmTZtGg8++CDRDf6r9NlnnxETE8O0adMAGDduHFdccQXLly9XABQRkVMWHw/793uvd1Y7WXfmNsoPupaHPdafvXMjSUqC3NymjxNprU4fACsrKykoKGh0m9Ppem5iRL0bJwzDwG63s2/fPoYP9xxVtX37dpKTkz3WDR48mBUrVrRx1SIiInX2PrqX8i3lAHQ5swtJs5P8XJEEuk4fADdt2sSMGTMa3fb0008zZMgQHn30UebNm0dNTQ3PPvss4AqODR07dswjLAKEh4dTXl7ebA0OhwOHw9HKd9D0Oet/F7VJQ2oPb2oTb2oTb/5pEwtgACYOh9Nji2malO84/jljhSFLh2Aa5vH6mj6uLQXbz0mwvM/mdPoAOHbsWLZt29bk9uHDh7No0SJ++MMfEhcXx4033shXX31Ft27dvPaNiIigrMEQrcrKSqKiopqtIScnp9ntpyIrK8tn5+6s1Cae1B7e1Cbe1Cbe2rNN7PZUIBS73c7GjY287t0QMTIC5z4nO9gBG1t4XBvTz0nw6PQBsDmmaVJSUsKTTz5JeHg4AF999RVRUVH069fPa//k5GT+/e9/e6zbsWMHQ4YMafZ1kpOTiYyMbLvCcf3vJCsri9TUVKxWa5ueu7NSm3hSe3hTm3hTm3jzR5vYbJbj322kpaU1vlMjq1t0XBsItp+T8vJyn3bedAYBHQANw2Du3Ln86Ec/4s4772Tv3r08/vjjTJ8+nZAQ77c+YcIEFi9ezOuvv860adNYv349H3zwAc8//3yzr2O1Wn32C+PLc3dWahNPag9vahNvahNv/mkTo5Wv2drjTk6w/JwEw3s8kYB/EsiTTz7JmjVrGD16NDNmzGDChAn86le/qt0+ceJEXnzxRQC6d+/Oq6++yieffMLYsWOZN28e8+bN4+yzz/ZX+SIiEoBM0+T7676n8O1CTNP0dzkShAK6BxBcl2eXL1/e5PaPPvrIYzk1NZW33nrL12WJiEgQO/DaAQr/UkjhXwpJ+CqB5OeTT3yQSBsK+B5AERGRjqQqr4ods3fULve4oocfq5FgpQAoIiLSXkyT7TO34yhxTUPSe3pvevxYAVDaX8BfAhYREfGVzEzXc34bzCDmJT/f9X1I/kGK/l4EQJnVxqGLBjPMxzWKNEYBUEREpJXS0yE7u2X7RmNnlrm9dvkJxxC+vsVG7/lNH+MOjiJtTQFQRESkldw9fxaL63m9TSkogDtqdhCLHYB/0ZN/EAc1rmf9nkiDR9eLnDIFQBERkVMUHw/79ze9/b15xXRb5Hpu/THDyjMMAdM4YXAEV/jLyGjDYkVQABQREfEp0zRJXPk/jh5fPnPZYEIXhEHuiYOjiK9oFLCIiIgPGYbB6R+fTtyUOGIujqHPTX38XZKIegBFRER8LbR3KD/I/AGOcgeGYfi7HBH1AIqIiLQXa6SeQSsdgwKgiIiIDxx89yD2Q3Z/lyHSKAVAERGRNnZ001G+n/o93w7/loPvHPR3OSJeFABFRETakLPGSfZN2Zg1JvZCO8c2H/N3SSJeFABFRETa0L7H93H0O9ekL1EjojjtgdP8XJGINwVAERGRNlK+rZzdC3e7FiyQ8koKllB91ErHo2lgREREGsjMdD3n1/2ot6a4n9Wbnw99E03mFWWTUm0CkHRPEl3HdPVxpSKtowAoIiLSQHo6ZGe3fH+nE87KyyWFUgByieCmPw+g+q3G93cHRxF/UQAUERFpwN3zd6Jn9RYUQE0N9KaCX/K/2vWLSWFX/onn/IuOPtVKRVpHAVBERKQJJ3pW74oVMH+eyR27c4iocgLwHglkWWJIbCY4giv8ZWS0YbEiJ0EBUEREpJWmTIHJk+HAH3ux854yDpRZWeoYeMLgKOJvCoAiIiKnwDAM4m+IJ/bSWCacUUX5wRC6+7sokRNQABQREWkDYfFh/C80zN9liLSIJicSERFpheqian+XINJqCoAiIiInqbqomrU/WMvWG7ZiL7b7uxyRk6ZLwCIiIidp5z07sRfaKfhjAYbVYOgrQ/1dkshJUQ+giIjISTi08hAFbxYAEBITwoDfDvBzRSInTwFQRESkhWpKa8i5Nad2edATgwiL18AP6XwUAEVERFpo59ydVO2vAqD7hO70uaGPnysSaR3dAygiIkEjMxPS0y0UF6diszXdB+J+Vm9+PiQluf48ynKYe/a5NliiLCQvTcYwDF+XLOITCoAiIhI00tMhO9sAQlu0v9MJubkQjoOr2Va7/rWQgVx3boTX/u7gKNLRKQCKiEjQKCtzfbdYTOLjARrvwSsogJqauuWb2EUilQD8l278sSQRs6Tp14mObpt6RXzF5wHwyJEjWK1WovXbICIiHUTPnnb27LFitVob3b5iBcyf7wqMhmnSL+8YAFVYeKNXCgm2pi/9RkdDRoZPyhZpMz4JgP/85z/585//zJo1a6isdP2PKTIykvPOO4+rrrqK8ePH++JlRURE2sSUKa4vF4OkxNM5I+8Acd2crCuI9GdpIm2iTQNgXl4ec+fO5fDhw1x22WVcf/319O7dG4fDQWFhIevXr2fRokX06NGD3/3udyQmJrbly4uIiPiGYbCSeBK7+LsQkbbRpgFw9uzZ/OpXv2LcuHFe24YMGcL48eO56667+Oqrr5g9ezZvv/12W768iIiIiLRAmwbAN998k5CQE5/yggsu0GVgERHpsJzVTrZctYXEOxOJnRDr73JE2lybBsCWhL/W7CsiInIyXPP91Y36dXNP01JUZKNfv8aPjY6Gx8/YS9T7hzj0/iEGLBoANLGzSCd1yils9erVnH766URFRbVFPSIiIqfMNd9f09udToPc3Ma39ecoYdl7XAtWiL00Fp5v+xpF/OmUHwWXl5fHtGnTOHjwYFvUIyIicsrq5vuDxMS6L4vFvd4kMdH02JaYCCGGk7lsIwQTgNPmnkb0KE1jJoHnlHsAJ0+eTEREBDfddBMffPBBW9QkIiLSJuLjYf/+uuWkJNeTPZqaB/CWbvsZVupKj5FDI+mXrku/Epja5Ea8yy67jGHDhrXFqURERPyiPKecKaW7AXACKa+mYA1vfKJokc7ulC8Buw0YMKCtTiUiItKuTKfJtl9sIxQnAJ9GJdFtXDc/VyXiO20+FPfo0aN06aKZMkVEpONqOAp4wtFcri9xPdw3l3CeKR/A8qS67e7RwyKBos16AN3OOussPv3007Y+rYiISJMyM2HYMNc9fklJdYEtP79uXVISFBS41rtGAbtGAh/MdXBFyZ7acy0mhQrTSm4utV9OV8cgeqy9BIo27wE0TZO//OUvvPzyyxiGwYgRI5g0aRJnnHFGW7+UiIgI0PS0L04nTU734laNlbutZ3K/dRs7qyPYRHcsFtcAkvqioyEjo+1qFvEnn8zGnJ2dzcSJE4mIiGDLli38/Oc/59prr+XBBx/0xcuJiEiQqz/tS3y8q+fP6aTRIAcmdrsdm81Gfr6B0wnOPhHcvvcMBvZ1Qp736GGRQOOTALhkyRKPR71t27aNmTNn0qdPH2666SZfvGSTCgoKWLRoEWvWrMFmszFx4kRmz55NWFhYo/svWLCAd955B5vNVrvu/vvvZ+rUqe1VsoiItJI7uLmne2ksyDkcTjZuzCItLY1+/ay1PYSGxcBuaNSvBIc2vwcwJiaGPn36eKxLSUlh3rx5vPXWW239cs1yOp3ccccdVFdX88knn/D++++TnZ3NwoULmzwmKyuLjIwMNmzYUPul8CciElhM0yTvxTzCnA5/lyLiF20eAIcOHcqKFSu81p922mnkt/Mwql27drF582bS09Pp3r07sbGxzJ49mw8++ICyhg+IBKqrq8nJyWHEiBHtWqeIiLSvmg9r2HnnTh4tXEsqR/xdjki7a/NLwHfffTczZsygsLCQa6+9lqFDh1JZWclLL71EUlLSiU9wkiorKylwD+tqwHl82FZERETtOsMwsNvt7Nu3j+HDh3vsn52dTU1NDU8//TTr168nOjqayZMnc/PNN2OxNJ2VHQ4HDkfb/i/Sfb62Pm9npjbxpPbwpjbxFjxtYgEMwMThcDayXKdiXwWVT1QC0MtRSRQ1FLfguEAWPD8nLsHyPpvT5gEwLS2Nv/71ryxatIgZM2Zgmq7nKYaFhfHUU0+19cuxadMmZsyY0ei2p59+miFDhvDoo48yb948ampqePbZZwFXcGyorKyMMWPGMH36dJ544gm2bt3KzJkzsVgs3HzzzU3WkJOT0zZvphFZWVk+O3dnpTbxpPbwpjbxFuhtYrenAqHY7XY2bsyqXc7Ph4SEeh/2psmckp2cVe1a/H/04ht6Ysk3SUhwUFRkOX4+13mCTaD/nEgdw3QnNB84dOgQmzdvxjRNTj/9dGJjY331Uk3at28fixYtYuPGjcTFxXHjjTfywAMPsHLlSgYNGnTC419++WVWrlzJu+++67WtvLycrVu3kpycTGRkZJvW7XA4yMrKIjU11etZlcFKbeJJ7eFNbeItWNqkXz8LubkGiYkme/Y4GTHCQna24bXfJRTwEFsBKMbGDYyhDJvXfkOHmmzeHFw9gMHwc+JWXl5OTk4Ow4YNa/PP787CJ6OA3Xr06MEFF1zgy5dolmmalJSU8OSTTxIeHg7AV199RVRUFP36eT/g+/PPP6eoqIhrrrmmdl11dXXtsU2xWr0fKN5WfHnuzkpt4knt4U1t4i142sTAarWSkQHz59dNDwPQzVHFXQXb4Xi3xxuxyRw7YoMG08W45vszgqS9PAXLz0kwvMcTafNBIB2JYRjMnTuXpUuX4nQ62b17N48//jjTp08nJMQ7+5qmyaOPPsrq1asxTZMNGzbwxhtvaBSwiEgnM2UKbN3qmgJm/37Yt8/k7bO3E23WABDywxDeKYytDX3u6WL273cdN2WKH4sXaQc+7QHsCJ588kkefvhhXn/9dbp06cKUKVO48847a7dPnDiRK664gttuu40JEybwwAMPsHDhQgoKCujZsyezZs1i0qRJfnwHIiJyqgrfLqTo70UA2OJshM4N9XNFIv7VpgHQ4XC0uFv1ZPY9FcnJySxfvrzJ7R999JHH8jXXXONxCVhERDq36oPVbL9ze+3yoGcHkRtzgufDiQS4Nr0EPG3aNFavXn3C/f75z38ybdq0tnxpERGRRtlibfR7sB+WcAtxU+KImxzn75JE/K5NewAXL17MAw88wG9/+1suv/xyRo4cSe/evXE6nRQWFrJ+/Xo++eQTunXrxu9///u2fGkRERHy812PgfNkAH3p07UH5V+HUNrPgt2eis1moZ2fTyDSYbRpAOzbty9vvvkm//jHP/jLX/7C0qVLqaioAFyTMZ977rnMmTOHCy+8sC1fVkREAlRmJqSne47mbYz7eQBOJ7XP9m0ol/rTfXjeAxgd3foaRTojnwwCufDCC7nwwgsxTZPDhw9jsViIiYnxxUuJiEgAS0+H7OzWHdvfWo69T2NzvJnY7XZsNhtgHJ/25VSqFOl8fDoK2DAMv0z+LCIigcHd81d/nr4Tyc+H85yFpDu+57RrkhiQMQBrRN2gQ4fDycaNWaSlpWk+OAlaAT8NjIiIdH7uefpaYnh8NXcf2I4F2L9kP93GdyPuZxr4IVJfQE8ELSIiwef6I9uJwQ5Azyt70vOnPf1ckUjHowAoIiIBo/CvhYytPAhAmSWE5OeTMQzvZwKLBDsFQBERCQjVhdVsn1k34fPr3ZIJ7a0nfog0xif3AF588cWN/o/LZrMRGxvLRRddxC9+8QssFuVPERE5daZpknNHDvYi16Xfr+jJmgjd9yfSFJ8EwKuvvpq//vWv3HzzzfTt25fc3FxeffVVzjnnHAYOHMif//xnKisrmTVrli9eXkREgszBvx6k6B3Xs35LLTaedCYTpUu/Ik3ySQD89NNPeemllxgyZEjtujFjxjBnzhzS09O58MILmT59ugKgiIicsuqCanJm5tQuvx4zhCPFoUT5sSaRjs4n12D37NlD//79Pdb17duXXbt2AZCUlERpaakvXlpERIKM46iDiAERAMRdFce3Eb38XJFIx+eTADh06FBeeuklj3WvvvoqgwcPBuCf//wniYmJvnhpEREJMhGDIhi5eiQDFw9kyHNDTnyAiPjmEvC8efP45S9/SWZmJvHx8eTl5WGaJi+88ALfffcdM2fO5KmnnvLFS4uISBCyhFg4bc5p/i5DpNPwSQAcPnw4n376KatWreLAgQPExcXxox/9iPDwcI4cOcI//vEPevTo4YuXFhGRIGCaJqbdxBKq2SREWsMnAfDo0aP8+te/5osvvgBczwR+//33efLJJ4mJifHFS4qISADLz4ekpLrlceWFTCrbw9LuQ/lfaFevfUWkeT75r9OSJUs4duwYH330EZs2beK9997D6XSyePFiX7yciIgEqOho13enE3JzXV8VuVVMP7ydpJpyFhz8jpDcY7XbcnNd+9Y/VkS8+SQArlq1iiVLljBo0CDCwsJITk5m8eLFfP755754ORERCVAZGTB0KCQmHv9KMJkXto2u1ADwbUQcNYlRdduPfw0d6jpWRBrnk0vAFRUVRDf4r1fXrl1xuv9bJiIiQS0zE9LToazs5I77ofMAI6uKAbD1tjFnSzIP6JZykZPmkwB4xhln8NRTTzFnzhwMw8A0TZ566ilSU1N98XIiItLJpKdDdvbJHdObSiazo3Y5ZVkKth62Nq5MJDj4JADee++9zJgxg/fff5/ExERyc3MxDIPXXnvNFy8nIiKdjLvnz2KB+PgT72+YJnfnZROFA4A+N/Sh5xU9fVihSGDzSQBMSUnh008/5fPPP6e4uJjExEQuuOACunTp4ouXExGRTio+HvbvP/F+uc/nsX3mEQCKrGGc++Rg3xYmEuDaNADm5eV5LJ9zzjm1fy4tLaW0tJSEhIS2fEkREQlw5TvK2XnfztrlZTEpTOnmk/4LkaDRpr9BF198MYZhNLrNNE0Mw2Dr1q1t+ZIiIhLgSleXYlabAPydBLaEx/q5IpHOr00DoHviZxERkbbSZ3ofokZEsfSc3bxUORAN+hU5dW0aABMTE9vydCIiIgBEj4zmDz1Sqcz1dyUigUEPURQREREJMgqAIiLS4ex5dA+7Fu7CadcDBER8QcOoRESkQynbUMbu9N2YNSbFK4sZ+e+RWGzqrxBpS/qNEhGRDsNZ5ST7+mzMGteo39hLYxX+RHxAv1UiItJh7H54N8eyjgEQdUYU/eb383NFIoFJAVBERDqEkv+UsPexvQAYNoNhbwzDEqqPKRFf0D2AIiLiN/n5kJQEYc4aHincSu/jYz7eiujPtMu6eO0rIm1DAVBERNpddLTru9MJublwLzvpTSUAm+nK0tLTcJY2f6yItJ4CoIiItLuMDJg/H8rKIK2iiMuLXd17lYaFV3sNIz6k8ceKRke7jhWRU6MAKCIibSIzE9LTXaHuZFx6rO7xHqe/NJjvfhnRxpWJSEMKgCIi0ibS0yE7++SPu5dUbmQXw7uUc8HN8W1fmIh4UQAUEZE24e75s1gg/iRyXH6+haXOQSR1NfmV0filXxFpWwqAIiLSpuLjYf/+lu+flOQaCGIq/Im0G02wJCIi7cp0mOTckcPRzUf9XYpI0FIAFBGRdrXvD/vIeyGP9aPWk/+qJvcT8QcFQBERaTdHs46y66FdAJh2k4hBGvEr4g8KgCIi0i6cVU62/nwrZrUJQN97+xJzQYx/ixIJUgqAIiLSLnYt2MWx/x4DIGpEFP0z+vu3IJEgFlABsKKigqlTp/Luu+96rN+1axfXX389I0eO5Nxzz+XFF19s9jx/+9vfmDBhAmlpaVx55ZVs2LDBl2WLiAS8I18fYd/v9wFg2AyGvTkMa7jVz1WJBK+ACYDbt29n2rRpbNy40WO93W7ntttuIzU1lTVr1rB06VKWL1/Oxx9/3Oh51qxZQ0ZGBr/73e9Yu3YtP/nJT7j99tupqKhoh3chIhJ4akpryJ6RDa4rvwzIGECXM7r4tyiRIBcQAXD16tVcf/31/OxnPyMhIcFj29q1ayksLOSuu+4iNDSU4cOHM336dJYvX97ouTIzM5k4cSKjRo3CZrNxww030L17d1auXNkeb0VEJOBsv3M7lbsqAeh2bjf6zunr54pEpFNMBF1ZWUlBQUGj2+Li4hg6dCirVq0iLCyM1157zWP79u3bGTBgAKGhobXrBg8ezNKlSxs9344dO5g8ebLHusGDB5PdzPONHA4HDoejpW+nRdzna+vzdmZqE09qD29qE2/t2yYWwABMHA4nAMe+P0bBn13/flu7Wkl+PRknTnA0f5wv6efEW7C1SbC8z+Z0igC4adMmZsyY0ei25557jksuuaTJY48dO0ZEhOc0AxEREZSXl7d4//Dw8Cb3B8jJyWly26nKysry2bk7K7WJJ7WHN7WJt/ZoE7s9FQjFbrezcWPd60Uui6RifgW2221kH8mGjS07ztf0c+JNbRI8OkUAHDt2LNu2bWvVsZGRkV7371VUVBAVFdXo/hEREVRWVnqsq6yspHv37k2+RnJyMpGRka2qrykOh4OsrCxSU1OxWnWjNKhNGlJ7eFObeGvPNrHZXHcVFRXZmDTpTM9tpgP7C1Z4wfu4oiL38TbS0tJ8WiPo56QxwdYm5eXlPu286Qw6RQA8FUOGDGH37t3U1NQQEuJ6uzt27GDIkCFN7r99+3aPdTt27OD8889v8jWsVqvPfmF8ee7OSm3iSe3hTW3irT3aJDra9d3pNMjNbbj1xB830dFGu/696efEW7C0STC8xxMJiEEgzRk7dizdu3dnyZIlVFVVkZ2dzZ/+9CemTJnS6P5Tpkzhgw8+4JtvvsFut/P6669z6NAhJkyY0M6Vi4h0LhkZMHQonB9Xwu1d95CUYJKYSIu+hg51HS8i7SPgewBDQkJ49dVX+c1vfsP48eOJjIxk+vTpXHnllQDk5eUxceJEli1bxujRoxk3bhwLFixg4cKFFBQUMHjwYJYtW0ZMTIx/34iISAc3ZQr89NIa1qVtpfJgJb9MK2bEuyOw9bD5uzQRaSDgAuCXX37pta5fv3688sorje6fkJDgNdHzpEmTmDRpkk/qExHpbDIzIT0dyspOvO8th7dzfrnrPup1q01+Pt7Kw791hUMR6TgCLgCKiEjbSk+HZmbCqnURhZyPa8qXY1hZYB9G/jYL8+crAIp0NAqAIiLSLHfPn8UC8fGN79OjppJ7C7fVPu3jje5DKCiJAGfLeg5FpH0pAIqISIvEx8P+/d7rTYfJxou3UlLgmly317W9eGt5b/r2pZHRwCLSEQT8KGAREfGtvY/tpeSfJQCE9QtjyPNDMAzDz1WJSHMUAEVEpNVKvy1l94LdrgULDPvTMGwxGvUr0tEpAIqISKvt/d1ezBrXjX/9HuxHzHkx/i1IRFpEAVBERFpt2PJhJNyWQNezu9IvvZ+/yxGRFtIgEBERaTVrhJXkF5JxlDuw2NSnINJZ6LdVREROmTVSz1YV6UwUAEVEpMUc5Q6+v/Z7yreV+7sUETkFCoAiItJiO361g8K3Clk3ah1FHxT5uxwRaSUFQBERaZFx5QXkv5zvWjAhYkiEfwsSkVbTIBARETmhRMq58XBO7fLzYcl8fUlUs8fk5/u6KhFpLQVAERFpVvcoJ7P4nkhcj3r7lN68fbgPHG7Z8dHRPixORFpFAVBERJr1+0E7icg5CkBeSAR/jRtCYgtvIIqOhowMHxYnIq2iACgiEqQyMyE9HcrKmt5nVMVB7inOBaAag2dif0CVJaQ22E2Z0k7FikibUgAUEQlS6emQnd309t5UcjPbapefZTDfFHapXZ4/XwFQpLNSABQRCVLunj+LBeLjvbeff+wwUUdqAFgTHsd3sQkkGq7BHU5n8z2HItKxKQCKiAS5+HjYv7/RLRz+Ipxd6buY/VEyv44xAEhKgtzcdi1RRNqYAqCIiDSp+/91J+biGAzD8HcpItKGNBG0iIjUctqdXusU/kQCjwKgiIgAruf8fnf2d+xdvBfTNP1djoj4kC4Bi4gIpmmSc3sOR787ytHvjlJdUM3gxwf7uywR8RH1AIqICHkv5VHwRgEAligL8Tc1MixYRAKGAqCISJAbVF3Kjrt21C4PfWUoUcObf86viHRuugQsIhLEulHNXcVbMB2ue/6S7k6i19Refq5KRHxNPYAiIkHKME3m8z09HFUAdDu3GwN/P9DPVYlIe1AAFBEJUleV7mIURwAI7RPK8L8Ox2LTx4JIMNBvuohIECp6r4ifHN0LgAMY/tfhhMWH+bcoEWk3ugdQRCQIRQyOID8kgviaCl4yBvGfa2NafGx+vu/qEpH2oQAoIhJgMjMhPR3KygAs2O2p2Lwu7UZRZo5iAvm8YyZBK57tGx3dBsWKiF8oAIqIBJj0dMjOdi8ZQGgTe4bwDn2914ZA797Nv0Z0NGRknEKRIuJXCoAiIgHG1fMHFgvEx5vY7XZsNhsjKg+zLawbdsPa6HH5+eB0usLf/v3tWLCItDsFQBGRABUfD3v2ONm4MYt+Rf3YPHEz0aOjGfHuCMISvAd8JCVBbisuBYtI56NRwCIiAc6530n2tdnggLI1ZRz44wF/lyQifqYeQBGRAFZTVkPF7Aqch50A9LiiB6f9+jQ/VyUi/qYeQBGRTi4zE4YNc13CTUqqm6blQJ7JM/E5OP/nCn+5IZFcuW4YfU8zavcdNgxWrPBj8SLiF+oBFBHp5DxH/da5yfwfoysPAVBGCL+uGUFuvvc/+/Pnw5Qpvq5SRDoSBUARkU7Oc9Svqwfwh858rmMf4HrSx3M9hkF4JIn1jnOP+nUfLyLBQwFQRCRAxMe7pm+5JO4Is4tyatdH/jqMjxbFYG0w+4tG/YoEL90DKCISYC46locNE4CEOxMIvaqpiaBFJFgpAIqIBJiXug/lHRLZFBbLwMcH+rscEemAdAlYRCTAOA0LzzKE03o4uTPE9Hc5ItIBqQdQRKSzM026YPda7TD0T7yINE49gCIinURmpmvKl4ajdi/I38tPyWV+XipJSdG18wDm5kJoqAU4E4sFunaFiIi649z75ed7zh8oIoEvoP57WFFRwdSpU3n33Xc91u/atYvrr7+ekSNHcu655/Liiy82e54f//jHnHHGGYwcObL2a+fOnb4sXUTkhNzz/eXm1n2l5ubzS3MXcVTzB3MDVblVOJ11xzidBk6nQU2NQXGx57Hu/ZxOz+Xo6PZ/byLSvgKmB3D79u38+te/ZsuWLUydOrV2vd1u57bbbmPChAksW7aMHTt2cOutt9KvXz9+/OMfe53n6NGj7Nq1iy+++ILExESv7SIi/tJwvr+RFUXcXbytdvubRn/CEsKIrYDSUnegM49/N2qPc3PPA1h/fXQ0ZGS00xsSEb8JiAC4evVq7r33Xm6//XYOHz7ssW3t2rUUFhZy1113ERoayvDhw5k+fTrLly9vNABu3ryZmJgYhT8R6bDi42HLX0vY9H/f4+7syySJr+L7sn+/574Oh5OEBAeFhaG18wS6uecBbLheRAJfpwiAlZWVFBQUNLotLi6OoUOHsmrVKsLCwnjttdc8tm/fvp0BAwYQGlo3D9bgwYNZunRpo+fLysoiIiKCn//852zfvp3ExERmzZrFRRdd1HZvSETkFCTaj5F1eRbOSlf8+09EL16oGESCYfi5MhHpLDpFANy0aRMzZsxodNtzzz3HJZdc0uSxx44dI6L+Xc9AREQE5eXlje5vGAapqanMnj2bhIQEPvnkE2bNmsWbb75JWlpao8c4HA4cDkfL3kwLuc/X1uftzNQmntQe3gK/TSzEUcXcQ/+lxlEDQMwlMbz0fQpmhUF+vklSkvcxRUWu270bbncN+jAAE4fD2fDAgBX4PycnL9jaJFjeZ3M6RQAcO3Ys27ZtO/GOjYiMjKSiosJjXUVFBVFRUY3uf/PNN3ss/+QnP+HDDz/k008/bTIA5uTkNLq+LWRlZfns3J2V2sST2sNboLZJWNVQfssWejiqALAMt1AzvwbL9dVABE6n0eyj3ZrabrNVsnHj974pugML1J+TU6E2CR6dIgCeiiFDhrB7925qamoICXG93R07djBkyJBG93/llVcYPnw448aNq11XXV1NWFhYk6+RnJxMZGRkm9btcDjIysoiNTUVa8MHeAYptYkntYe3QG+TMTWF9Md19SJiSASnf346ob1CeewxWLDA9Joexs012MPAYjE9BoGAa9DHww+HNvkf3EAU6D8nrRFsbVJeXu7TzpvOIOAD4NixY+nevTtLlizh7rvvZteuXfzpT3/innvuaXT//Px8MjMzWbZsGfHx8fz9739nw4YNPPzww02+htVq9dkvjC/P3VmpTTypPbwFapv8Kyqeo0dMbrTs5sefnk5EvOv2lquvdn01xuFw1BsEYjQx2CPw2qolAvXn5FQES5sEw3s8kYAPgCEhIbz66qv85je/Yfz48URGRjJ9+nSuvPJKAPLy8pg4cSLLli1j9OjRzJ07F4vFwnXXXUdZWVntgJF+/fr5+Z2IiMCHJPB9795cOUAfYCLSegEXAL/88kuvdf369eOVV15pdP+EhAQ2bNhQuxwaGsqDDz7Igw8+6LMaRURaqmJnBRGDPAeyVVkU/kTk1ATUk0BERDqjzEwYNsw1L1/9rzu67eHfg79lcmyhx6Pa3I9uS0pyHbdihX/rF5HOJ+B6AEVEOhv3I97qm8FurmY3ADMPf893hyNx0gWoe3Sb2/z5MGVKOxUrIgFBPYAiIn5W/xFviYkwK3oXNx4PfwB/6ToIe2IXLBbP/dzLTY3+FRFpinoARUQ6iPg+Jv/6xW72/GZP7bpBSwaxbHZfoO7RbQ25Lwk3rm4iaBERNwVAEZEOwWRKWYPw94dB9L27b+1ydLTre8NLwA2XPdU9Hs59vIiIAqCIiL+ZJr9gFz8t21u7avBTg0m6y7NbLyPDdb+f+5Kva4Jn16XghhM81zs5drud2FgbGRl6VrCIuCgAioj42ZVlu7mSeuHvmcEk3el9TXfKFM/BHu5LwvHxNDHBMzgcTjZuzCItLU2T34pILd0YIiLiZ9lhMVQe/+d4yLNDGg1/IiJtST2AIiJ+9n1Ydx4ilR/EVPDGzAR/lyMiQUABUETERzIzXXP8eU3TYppg1N2Pl58PuXRnY2l3vjze+Rcd7brnT/P7iYgvKACKiPhIYxM8W3EymxwKCOcN+nts0wTPItJeFABFRHyk/gTP8fEQ4azhruItpFYdBly9fH/v2t9rNK97WRM8i4ivaBCIiIiPxcfDjtWV/LnnhtrwZ4QZpC+LZP/+uilc3KN5m57SRUSkbagHUETEx06rLuO7s7OozqsGICQ2hBHvjSDm3Bj/FiYiQUsBUETEh87iEPOLvqfadAAQPiic01eeTmRypNe+7ke65ed7LjfFvZ+IyMlSABQR8ZELj+VxIzlYTddy17O7MuL9EYTGhXrs17pHvHkfLyLSUgqAIiInqcnpXer5v6O53FyyvXb52/CevLB3GOHnW72md2ndI95c3NPFiIicDAVAEZGT1Nj0Lg19RE8uYy+9qeJtknipchBmngF53tO7tOYRbyIip0IBUETkJDWc3qVxYTyQdzojzCN8ZEkkQdO7iEgHogAoItJK9XvojvzrCBGDIghLCKvdnpQUxQe5USQe38/dsyci4m+aB1BE5BSYDpPdv93Nxgs3svXnWzEdpr9LEhE5IfUAiogc15LBHVA3/UplXhXPRW1lRNURAI6sOsKBNw4Qf6NmchaRjk0BUETkuJYM7nA7k8M8ZH5PbJUdAAewsmd/LpjRx3cFioi0EQVAEZHjWjK4w2I6mZC3h2nsqb2HpohQfsswisK6s8TqfYwmeBaRjkYBUESkgaamX6nKreL7676nJK+kdl3sj2K5Y+NQth4IJbHB/prgWUQ6KgVAEZEWqC6sZl3aOuxFdZd8M7sO5IWP+lJ6mtHoMZrgWUQ6KgVAEQl4Jzu4o/4lWncQmzIllJ6Te5L/Uj5F1jAWOoZzJLobLzYzl4ImeBaRjkoBUEQC3skM7gD3JVoTC5CLUfvkjsF/GIw10sptf+nHtgM2r0u+IiKdhQKgiAS8lj25o+4S7WlGOXNDc/hndQ/+avatPd4aYWXwE4M5+lff1ywi4ksKgCLSaXlf2rVgt6dis3lel21sVG3dpd26dQMSHVyYt5frzL3YqkwGGaWsIg4I99VbEBHxCwVAEem0vC/tGkBok/s3HH3rvrQLUPx5MY8WbKcPFbXbyyw2ejqqKW0iAGp6FxHprBQARaTT8r60a2K327HZbLjCoEvD0bfu5bIyqDpQxc57d1L450LcUzjXYDDw/r7c9Md+7Mq3anoXEQk4CoAi0uGc7Khdt+houOGGfcyZ0x+rtW5G5oajb5OS4FBuDZeX5fHt0D04Shy1+/6XbizvlcyaR6Oo/lPjr6vpXUSks1MAFJEOp3WjdgEMXnwxgTlzTnzMZRzgmtL/4Y5+IbEhPO8cxFtH+pBga3xePzdN7yIinZ0CoIi0m9b07Fksrq+uXSEiwns/d88buP5cXt7Is9ga8RHxTLfsIQY7fW7ow8DHBnJNWijmkZa/HxGRzkoBUETazcn27IEr1Dmd0KsXbN3qua1+zxs0fv+d45iDiWW5OKjhEwbWrq/CytLuQ1nx7wgiUyI9jtHgDhEJdAqAInJKWtqrB6fWs9eS89fnOOYg9/lc9i3ex7WlduwYrK1JoP6ULv8N70FkSt0xGtwhIsFCAVBETklrevXg1Hv2mtLTUcnuBQc5sPQA9oP22vVWTEZUHQaaHqWhwR0iEiwUAEWkUa29X8/93Vc9e42xmk5GVh7iTvIZc6iYfYvqbTRgdXgvXqzohyMqqtnzaHCHiAQLBUCRINPSYFdQADU1J3dup7Pue1v37DXnruItjKo85LnSCnFT4uif3p+f/zCK3Fz07F4RkeMUAEUChC+DXX0hIdC7d91yc5Ms+4LT7sQIMTCMuqla1oX3rA2AhZYwzlrQh4RfJBCWGOZxrAZ3iIi4KACKBIjW3ovn1tJg17u35+XQxiZZbquevdracDKi8gg5M4s4uOIgP8j8ATHnx9RuXxPRixFVh3mvog97e0SR+1CIx0TQGtwhIuJJAVCkg2u6Z8+C3Z6Kzea6oa6l9+J1pGDXnCjsjCsv5szKIlIppsshB3nPu7blLc3zCIDVFisvxA4nNxd6GdVe59LgDhERTwqAIn5y6pdsDSC00WOauxfP38GuOZX7Kyn6exG/LjrEMI4Qctj02scIdV3+NU0T9/N+64ffoiIb/fq17PU0uENEgpUCoEgba7978UxiY+3YbDbAaPd78Xwh99lc9j22j9QG60sJ4b8RPbjljR7EXhpLSLTrn676l3bdnE6jxYFWl3ZFJFgpAIq0UEcbZNG7N7z3XhZpaWlYrdYO3bMHYJomPWqqGEYJZx8+QnJ1Cb9lGMeoS2E9f9qTfY/tAyCfcDZF9WB9RE9WFXWjT6yFeVM8z1n/0m5dO5nHL+s2/zxfXdoVkWAWcAGwoqKCG264galTp3LllVd6bNuzZw9XX30177zzDknNDP3729/+xvPPP8/BgwcZOHAg8+fPZ+TIkb4uXfykIwa7jnQvXmuFOR2UrjlG2foySr4uoeRfJTxVUOXaWO76lsoRvqkXALuO6crAxwYydUks3xRGYakwoAKcnHjErlvPnnb27LF6DAIRERFPARUAt2/fzq9//Wu2bNnC1KlTPbZ98cUXzJ8/nyNHjjR7jjVr1pCRkcGyZcs4/fTTWb58ObfffjurVq0iouGsttKhKdi1k3r34rktIouz8w/x3dlNH+YAeuI5YMOwGJw29zSOvAYUNry027L2iYx0AAp/IiLNCZgAuHr1au69915uv/12Dh8+7LHt2Wef5eOPP+aee+5h3rx5zZ4nMzOTiRMnMmrUKABuuOEG3n77bVauXMnkyZN9Vn8ga2kQq6iA0tK6D33vEayuUa81NZYT7OeiYNe2Ip12kqnkBxUV7H2skvJt5RzLOsZvD5rcyGiPfY9hxdLgeEukhf86urKuqhubjRi2Ek25GYKlkZ69igpX+7sfFwctHbFrcsMNeUD/U3mrIiIBr9MEwMrKSgoKChrdFhcXx9ChQ1m1ahVhYWG89tprHtuvuuoqZs6cSW4LPoV37NjhFfQGDx5MdjMTrDkcDhwORwvehacVK2DhQkujwaiiwkJJyZmuzhU8R0LWDzwNQ1PD7XXn889+4A5izd+P1RinE4qL669pfNSr937NCwkxGwl2Ru29Y+7l3r1N9uype4P9+lnIzTWIj3etdy+DicNRryGwHK/Vvb7hclvt51L3s3cq5wPDNHE46n7WTqsu52b20q+wgq+7V7L0yPE0XQz/u7/ubKcBYTgAS+35sunKQFsF50yPJCo1iq5ndyUqLYpbR4aQnW14/Ei3tGcvOdlk82Zns/s4HA6yso606vcxULnbQm1SR23iLdjaJFjeZ3M6TQDctGkTM2bMaHTbc889xyWXXNLksb3rf9qfwLFjx7wu9YaHh1NeXt7kMTk5OS0+f3333z+c3btP/rLyiQJPSwORv/YDsFpNevSwe6wrKrLhdLqCiYsrlPXseWr7ube5l2Nj7bz3XlbtfpddlkphYSg9e7rWu5ftdjsbN9btZ7enAnXrGy77Yz+ArKysZvfrWnmU9X/ZillsYhaZXH6kF2HUkFBUyb/PPMIT+Ta6Yef5o4PZuLFuPhWjpj8XchDsUHOkib9IAw5YwunhqCY/P5yEBAdFRRbeJYm/OxLp+aEdPnTvbFJZ6cBqtWKa7v88NP5311BkpIMbbshj48amCvHkbhOpozbxpjbxpjYJHp0mAI4dO5Zt27b5/HUiIiKorKz0WFdZWUn37t2bPCY5OZnIyMiTfi273dXzUjdqsY67FwrM2kl9wfNDs35vlXu/htubOl977Vd/X4vFPP4eDPr0gT17PO/T6tfP1ROUePyBre7LqO79XL07WUyadGaz+zU8n3ube9lms5GWlla7n3siZff6hssdYb+eVDLcWclp+adxrv0wlTiIr3JQtayKHmE9cBxx8OvSbDYTzWe2/h7n+82RzZRPr6o93zXsdf2hCpz/hV64tvUwHB6vWxzm6vFzAFEDwlmXF87uqgiOdQ1nwYthhA8MJ3JYJGeMtZGXbYATCgvremidTsNjuSnJybB584nu2bPSksu67p+R1NRUDQI5Tm3iTW3iLdjapLy8vNWdN4Gi0wTA9jJkyBC2b9/usW7Hjh2cf/75TR5jtZ7aiMP4eMNrMtqkJLM25Ozfb9Rb7w41hsf9ZYmJnssNz+mv/RrWDO7LfUYzbWZ4/Lkt9vPcVrdsOk2spkkYTsKdJs5Sky4OB90xiXWYOEus2GJtHmcbVF1G2ddOflDpJBGT3hUmxX83cVY5cVY5ufiYk2OY7KrpgdVa18Pbx17Bjjv24axw4ih3MKfIiYmT6EIH36U5WXLAgRUnEXkODHM8lpC65H8V+7m6aD9broBZ7pUlUP0S5OPqtTsTACefNXi/Ry02ejrqAmBDZZYQDjrDqDI8H59WZrEwjbEcNMLoVW0h3+4ajWs5Ch/cV3d8/fv14GTv2YOMjOb+7lrnVH8nA5HaxJvaxFuwtEkwvMcTUQBsYMqUKcycOZMf//jHjBo1iuXLl3Po0CEmTJjQrnWcXnmIn3OQyMOw9QZqr3beWuyaQSOyGLZeD33tfcmlS+1xSfajzGCfa/t0jj8tAW4vhorjx31/HdxRDPMYRv0QdeBPB7iz+BCVQEQxbLkKZh06ftwhk81TXHV0OaMLDXtjZhzZThcqiSgy+e/lx+s9/jW3yKQKCC+CL6PiyaVX7XHVB6vZMmULmDD/INgxCT0Ihgk1mIQWwvqzXEFt6F+GerzmWRWFXFW6CwcQesDkm4FgOkxwur4/V2BiAmUFocBZHsduuWYLB98+CMAf3Svz4d+x8GJtg8D/fh1PyrIUj2MfKNrIxgudPOBeUQxb6s1Pd9Px70/Zw4C6ABjttJO/tO4Sa5r7D3Y4thlqb1QwwVnhxBJdFwDLWziqtSvel1LXhvfknJu7YuthIzQ+lJnzQ8kpDuWwEUZYfCj7Dlhcwa4Uvqo3GKOg0KCGCDBb9/zc5GTPJ5CIiEjHEfQBMC8vj4kTJ7Js2TJGjx7NuHHjWLBgAQsXLqSgoIDBgwezbNkyYmJi2rWuvvZyfkQBlEPBH+vWn+f+QwUUvAHde/SCegEwxlHNDylwbX+z7rjx9Y4r/AucA8Awj9c8uukoZ1ccrN3v4AoY695YCUXvuP7oKPW+eXZY1WH6Ug5VUPyR57bT3X+ogqwwz0vpZrVJyT9LAKiNWfVnBrFD2TrXH50VngMAIpwOEmoqjhcFlbs8X7eb+zWcjQxA8X7CWKPMGu8dHVhw9YU1L8T03KfGaHogjCXCQkmVhXKnFWeIBbPa83U3042VXZK4/V4rj/whhPxSK+Hdrdy3YDdD0oYQ1iuMsy62kXMgBPP4qFr349FePNaf996vO1dBKdSAqw3y6ta3NNg1HB3dGE2yLCLSsQVkAPzyyy8bXZ+UlOR1H2FCQgIbNmzwWDdp0iQmTZrks/oCU8tG+XrtdYLDnIDFCoZheGUuu2HhmBFCjemaP65HT9d3rK7v+/IM7A44FuJ9H1rE4Aiix0RjhBh8u86gvNrAFmZwwUUGn//D4FilQVg4XDeqi9ex/69LIrNuc/LMiwaHyyxEdTW4f56BJcyCJczCvfcbFB6xcCi0q8dxuSFRjN44GkuEBUuEhbQxVnYdsBCXYGF/rlF3Wb03XNvD8zXXE8uBbrH8fiF8/DLklkJipEnIufvpltbNNbAiBhwHXI3Wmh67kwl2U6Y0v5+IiHRsARkAA8E/ovrwcWlPevcy+eqfdSnpgvOhoBB694Kv/gXbLgzzOC4ntBvTGEOf3vCvr13rDMNg/Hg4UOD6gP/PajhnHJgNZtXp91A/Jr+ZVLvft2thzBjIPwDxfWDtOlcdRphR7/qly2/j0ijIN+nTB7I2H6/XcH394AeQl28QHw/VhgVK644LjQ/lvIrzMCwGAwbA/jxISDAwgdw84/g9kK59HQ4HbKw7dnVkb1ZH9naFpni87j28yh2o4uCWBu07cNFAWOT680/c+/WEWR/DZe7lHjD3Tu+/mxVdB/Dk72HFnyG3DBKj4fl698T9KwNyj0Big9+uaov1+OVzl6PW452dLZwhJ79Bz15+vmsEs3uwx6nfi6dgJyISLBQAO6hjFhu5hOK0mUSm1CWEfBvkAk4bRCZDVYPZdqstVvKIxAiByMF16wtDIB+whEDEADjYyN+8rbuNQyE2CgFbCIT3hWIrFAFhVghLbL7eEqCLFWwNeq/KLXAMqGg4MzCucGoNPz7K93gnn3nyUwZ2ao0Fu/oTI7unv2zYk3cyo2x1L56IiNSnACjiI60Pdi29ZGsSG2vHZrPRVDei7sUTEZHGKACKnCTfB7uW3Yv38MNOBg3KIi0tTVMaiIjISVEAFDmuowS7lt6L53DAxo0n3k9ERKQhBUAJeJ0t2ImIiPiaAqCckhOFK/c68AxTrdvPgt2eSlFRy86nYCciItI4BcAO4ERhqP7tXe6pPRoGL/d+Dbc3db5T3e9Uw1Xr9jOAxke9tvR8oGAnIiKiAOhH0dGu742HF6P2u9PZcFvjx9Tfr7lA1Nb71ddYuKqogNJSz3npunaFiIiT3c/EbrdTU2OjtNQ44fkao2AnIiKiAOhXGRkwfz6UlXlvq6gwKSkB16N8Paf4qB94Goamhtvrzte2+zXG1+HK4XCycaNGvYqIiJwqBUA/mjKl6bDkCjsbFXZERESkzTXybAYRERERCWQKgCIiIiJBRgFQREREJMgoAIqIiIgEGQVAERERkSCjACgiIiISZDQNzClwHp8sr6Kios3P7XA4ACgvL9c0MMepTTypPbypTbypTbypTbwFW5u4P7edjT1pIUgYpumaalhO3qFDh9i9e7e/yxAREZFW6N+/Pz169PB3GX6hAHgKampqKCkpISwsDItFV9NFREQ6A6fTSVVVFd26dSMkJDgvhioAioiIiAQZdVuJiIiIBBkFQBEREZEgowDYCVRUVDB16lTeffddf5fiF4cOHeKOO+5g9OjRjB07lkWLFlFTU+PvsjqE4uJiJkyYwJo1a/xdit9lZ2dz4403MmbMGMaPH8/cuXMpLi72d1l+tXr1aq666irOPPNMxo8fT0ZGBpWVlf4uq0NwOBxMnz6d+++/39+l+N3KlSsZPnw4I0eOrP267777/F2W+JgCYAe3fft2pk2bxsaNG/1dit/cfffdREZG8q9//YsVK1awevVqXn/9dX+X5Xfr169n6tSp7N2719+l+F1lZSU333wzI0eO5Ouvv+bDDz/kyJEjPPjgg/4uzW+Ki4u59dZbufbaa1m3bh1/+9vf+Pbbb1m6dKm/S+sQnn32WdatW+fvMjqErKwsJk2axIYNG2q/Fi9e7O+yxMcUADuw1atXc/311/Ozn/2MhIQEf5fjF3v27OHbb7/lvvvuIyIigr59+3LHHXewfPlyf5fmV3/729+YM2cO99xzj79L6RDy8vIYOnQoM2fOJDQ0lO7duzN16lTWrl3r79L8JjY2lv/85z9ceeWVGIbBkSNHqKqqIjY21t+l+d3q1av57LPP+OEPf+jvUjqErKwsRowY4e8ypJ0F59jnDqKyspKCgoJGt8XFxTF06FBWrVpFWFgYr732WjtX1zFs376dmJgYevfuXbtu0KBB5OXlUVpaSteuXf1Ynf+ce+65XHHFFYSEhCgEAgMHDuTll1/2WPfpp5/ygx/8wE8VdQxdunQB4IILLqCgoIDRo0dz5ZVX+rkq/zp06BAPPfQQzz//vK4k4JoOZcuWLURERPDyyy/jcDi44IILmDNnDt26dfN3eeJD6gH0o02bNvHDH/6w0a///Oc/dO/enbCwMH+X6VfHjh0jIiLCY517uby83B8ldQhxcXFBO3fViZimyR/+8AdWrVrFQw895O9yOoTPPvuMf/7zn1gsFu666y5/l+M3TqeT++67jxtvvJGhQ4f6u5wOobi4mOHDh3PppZeycuVK3nrrLXbv3q17AIOAPkH8aOzYsWzbts3fZXRokZGRXo/acy9HRUX5oyTpwI4ePcoDDzzAli1bePPNN0lJSfF3SR1CeHg44eHh3HfffVx11VWUlJQEZe/OSy+9RGhoKNOnT/d3KR1Gz549PW6piYiI4L777uPqq6/m6NGjtb3IEnjUAygd2pAhQzhy5AhFRUW163bu3EmfPn2Ijo72Y2XS0ezdu5fJkydz9OhRVqxYEfTh77vvvuNHP/oR1dXVteuqq6ux2WxeverB4r333uPbb79l9OjRjB49mg8//JAPP/yQ0aNH+7s0v8nOzubxxx+n/jMhqqursVgshIaG+rEy8TUFQOnQ+vfvz6hRo3jkkUc4evQo+/bt4/nnn2fKlCn+Lk06kJKSEq6//nrOPPNMXnnlFQ10AFJSUqisrGTJkiVUV1eTm5vLY489xpQpU4L2g/2TTz7hu+++Y926daxbt47LL7+cyy+/PKhHA8fExLB8+XJefvllampqyMvLY/HixfzsZz8L2p+TYKEAKB3e008/TU1NDf/3f//H1VdfzXnnnccdd9zh77KkA3n33XfJy8vj448/ZtSoUR7zmQWrqKgoXn75ZbZv38748eOZPn0655xzTlBPjSPe+vTpw0svvcQXX3zBmDFjmDx5MqmpqaSnp/u7NPExPQtYREREJMioB1BEREQkyCgAioiIiAQZBUARERGRIKMAKCIiIhJkFABFREREgowCoIiIiEiQUQAUERERCTIKgCIiIiJBRgFQRDqc/fv3k5KSwv79+332GhkZGXz44YeNblu5ciXjxo1j1KhRrFq1qtHa0tLSWLRoEeB67NqaNWu8zvPMM88wffr0E9YyduxYRowYUbvv+vXr9bQbEfEpBUARCTqrV6/m+++/5/LLL290e2ZmJhMnTmT9+vVcdNFFje7z4Ycf8tBDD7VJPWvWrOHWW2+tXR41ahSRkZGsWLGiTc4vItKQAqCIdHi5ubncfffdjBs3jvHjx3PvvfdSWFhYu/0///kPP/3pTznzzDO55pprWLx4cbM9b0uWLGly+5QpU/jmm2946623uOSSS9rsPfzmN7/xeEbx6aefTkpKCqtXr250/+nTp/PMM89QXV3dZjWIiLgpAIpIh2a327npppuwWq189tlnfPzxxwDcdttt1NTUsH//fm677TauvfZavv32W+bMmcPbb7/d5Pn++9//snPnTi6++OJGt69YsYLRo0dz66238vnnn7e4zttuu43Ro0d7fC1durR2e3p6Ohs2bGDDhg2sWbOGM844gx//+MecffbZjZ7vjDPOwGaz8eWXX7a4BhGRlgrxdwEiIs1Zt24d+/bt45133qFLly4APPzww4wZM4bNmzezevVqhg0bxtSpUwEYPXo0V199NVlZWY2e75tvvmHYsGGEh4e3aZ0vvvgiY8eO9Vj3zDPP8O2333qsM02TuXPnYrfbeeyxxzAMo8lzpqWlsXr1an70ox+1aa0iIuoBFJEO7dChQ3Tv3r02/AF06dKFmJgYcnNzyc/PJzEx0eOYvn37Nnm+/Px8evfu7bN6T+TRRx/l+++/5/nnnycsLKzZffv06cOBAwfaqTIRCSYKgCLSoSUmJnL48GGOHj1au66srIzDhw8TFxdHYmIieXl5Hsc0XK7PYrHgdDp9Vm9zXnvtNd577z2WLl1KbGzsCfd3OBxYLPpnWkTanv5lEZEOLTU1lcGDB7NgwQLKysooKytj4cKFnHbaaZx55plMmjSJrVu38ve//x2Hw8GmTZv461//2uT5EhISKCgoaMd34LJy5UqefvppXnjhBfr379+iYwoLC0lISPBtYSISlBQARaRDCwkJ4aWXXqKmpoZLL72Uiy66CLvdzmuvvUZISAh9+vTh6aefZtmyZYwePZrHHnuMc889F5vN1uj5xo8fz5YtW6iqqmrX97FkyRIcDge33367x2jgF198scljvvvuO84777x2rFJEgoUGgYhIh5OUlMS2bdtql+Pj43nqqaca3Tc/P5+4uDg++uij2nW/+93vmjz30KFDGTJkCF988QWXXXZZo/v86U9/Oql669da36xZs2r//MUXX5zUOTds2IBpmpxzzjkndZyISEuoB1BEOrXDhw9z3XXXsXnzZgCys7N5//33m5zAGeDee+/lj3/8Y3uV2Cqvv/46s2bNIjQ01N+liEgAUgAUkU5t+PDhPPTQQ8yePZuRI0cya9Ysbrnlliaf8gGuy8DDhg3j/fffb/XrXn755bWPgjtVY8eO5aWXXqpdXrduHVVVVUyePLlNzi8i0pBhmqbp7yJEREREpP2oB1BEREQkyCgAioiIiAQZBUARERGRIKMAKCIiIhJkFABFREREgowCoIiIiEiQUQAUERERCTIKgCIiIiJBRgFQREREJMgoAIqIiIgEGQVAERERkSCjACgiIiISZBQARURERIKMAqCIiIhIkFEAFBEREQkyCoAiIiIiQUYBUERERCTI/H/VZAQToqd/PwAAAABJRU5ErkJggg==" } }, "cell_type": "markdown", "id": "e0f12880", "metadata": {}, "source": [ "**Plot of the final fit**\n", "![final_fitting.png](attachment:final_fitting.png)\n" ] }, { "cell_type": "code", "execution_count": 9, "id": "a01483e5", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Do you want to use an existing file to save fit results? \n", " eg: existing file to save HN parameters, y or n:y\n", "Enter the analysis_file_name:cond_fit.TXT\n", "file exists\n" ] }, { "data": { "text/plain": [ "()" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# before saving fit results an analysis file has to be created using create_analysis file method\n", "\n", "cond.create_analysis_file()\n" ] }, { "cell_type": "code", "execution_count": 10, "id": "3acb6038", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "()" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#save the fit results using save_fit method of the corresponding fit function\n", "#takes one argument, read more on the documentation\n", "\n", "cond.save_fit(1)" ] }, { "cell_type": "code", "execution_count": null, "id": "f7cd2fd4", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.9.7" } }, "nbformat": 4, "nbformat_minor": 5 }