{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Example of Signal Recovery 7270 driver usage"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-13T11:47:53.439509Z",
     "start_time": "2022-01-13T11:47:53.429509Z"
    }
   },
   "outputs": [],
   "source": [
    "\"\"\"\n",
    "General driver notes:\n",
    "\n",
    "With the SR_7270 read commands end with the echo \"\\n\\x00\",\n",
    "and write commands return only \"\\x00\".\n",
    "\n",
    "ask_raw function has been rewritten to read echo to remove from buffer.\n",
    "write_raw function has been rewritten to also read after writing using ask_raw.\n",
    "\n",
    "Get commands ending with . as in 'x.' are known as floating point in manual.\n",
    "\n",
    "Never change noise mode as TC values will not be correct.\n",
    "(If you want to change this then the driver will need to updated.)\n",
    "\n",
    "Change I_mode for choosing to measure current vs voltage using:\n",
    "'CURRENT_MODE_OFF', 'CURRENT_MODE_ON_HIGH_BW', 'CURRENT_MODE_ON_LOW_BW\n",
    "\n",
    "Change V_mode for measuring A, -B or A - B voltages using\n",
    "'INPUTS_GNDED', 'A_INPUT_ONLY', '-B_INPUT_ONLY', 'A_B_DIFFERENTIAL'\n",
    "\n",
    "Reference frequency is set by INT, EXT_REAR or EXT_FRONT.\n",
    "If in mode INT then oscillator frequency changes reference frequency.\n",
    "Otherwise, reference frequency controlled by input.\n",
    "\n",
    "\"\"\""
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:20.347509Z",
     "start_time": "2022-01-14T09:36:18.707509Z"
    }
   },
   "outputs": [],
   "source": [
    "import qcodes as qc\n",
    "from qcodes import Station, load_or_create_experiment, \\\n",
    "    initialise_database, Measurement\n",
    "from qcodes.dataset.plotting import plot_by_id\n",
    "from qcodes.utils.dataset.doNd import do1d"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Connect to device"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:25.398509Z",
     "start_time": "2022-01-14T09:36:20.893509Z"
    }
   },
   "outputs": [],
   "source": [
    "import qcodes_contrib_drivers.drivers.Ametek.SR_7270 as SR_7270\n",
    "SR = SR_7270.Signalrecovery7270('SR', address='USB0::0x0A2D::0x001B::13086379::RAW')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:25.458509Z",
     "start_time": "2022-01-14T09:36:25.445509Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SR:\n",
      "\tparameter    value\n",
      "--------------------------------------------------------------------------------\n",
      "IDN           :\t{'vendor': 'Ametek', 'model': '7270', 'serial': None, 'firmwa...\n",
      "I_mode        :\tCURRENT_MODE_OFF \n",
      "V_mode        :\tA_INPUT_ONLY \n",
      "frequency     :\t77 (Hz)\n",
      "noise_mode    :\tOFF \n",
      "osc_amplitude :\t0 (V)\n",
      "osc_frequency :\t77 (Hz)\n",
      "phase         :\t0 (Degrees)\n",
      "r             :\t0 (V)\n",
      "reference     :\tINT \n",
      "sensitivity   :\t0.1 (V)\n",
      "timeconstant  :\t0.1 (s)\n",
      "timeout       :\t5 (s)\n",
      "x             :\t0 (V)\n",
      "xy            :\t0j (V)\n",
      "y             :\t0 (V)\n"
     ]
    }
   ],
   "source": [
    "SR.print_readable_snapshot(update=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Setup db and station for test measurerment"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:36.275509Z",
     "start_time": "2022-01-14T09:36:34.595509Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Upgrading database; v0 -> v1: : 0it [00:00, ?it/s]\n",
      "Upgrading database; v1 -> v2: 100%|█████████████| 1/1 [00:00<00:00, 125.00it/s]\n",
      "Upgrading database; v2 -> v3: : 0it [00:00, ?it/s]\n",
      "Upgrading database; v3 -> v4: : 0it [00:00, ?it/s]\n",
      "Upgrading database; v4 -> v5: 100%|█████████████| 1/1 [00:00<00:00, 142.85it/s]\n",
      "Upgrading database; v5 -> v6: : 0it [00:00, ?it/s]\n",
      "Upgrading database; v6 -> v7: 100%|██████████████| 1/1 [00:00<00:00,  9.52it/s]\n",
      "Upgrading database; v7 -> v8: 100%|█████████████| 1/1 [00:00<00:00, 111.11it/s]\n",
      "Upgrading database; v8 -> v9: 100%|█████████████| 1/1 [00:00<00:00, 333.30it/s]\n"
     ]
    }
   ],
   "source": [
    "# create an empty database based on the config file\n",
    "qc.initialise_or_create_database_at('./test_signal_recovery_7270.db')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:36.469509Z",
     "start_time": "2022-01-14T09:36:36.297509Z"
    }
   },
   "outputs": [],
   "source": [
    "exp = load_or_create_experiment(experiment_name='test_signal_recovery_7270',\n",
    "                          sample_name=\"output2input\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:36.509509Z",
     "start_time": "2022-01-14T09:36:36.497509Z"
    }
   },
   "outputs": [],
   "source": [
    "station = qc.Station(SR)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:37.829509Z",
     "start_time": "2022-01-14T09:36:37.803509Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Lockin x value is 0.0 V\n",
      "Lockin y value is 0.0 V\n",
      "Lockin xy value is 0j V\n",
      "Lockin R value is 0.0 V\n",
      "Lockin Phase value is -180.0 deg\n",
      "Lockin Frequency value is 77.0 Hz\n"
     ]
    }
   ],
   "source": [
    "print('Lockin x value is ' + str(SR.x()) + ' V' )\n",
    "print('Lockin y value is ' + str(SR.y()) + ' V' )\n",
    "print('Lockin xy value is ' + str(SR.xy()) + ' V' )\n",
    "print('Lockin R value is ' + str(SR.r()) + ' V' )\n",
    "print('Lockin Phase value is ' + str(SR.phase()) + ' deg' )\n",
    "print('Lockin Frequency value is ' + str(SR.frequency()) + ' Hz' )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:39.313509Z",
     "start_time": "2022-01-14T09:36:39.124509Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Reference frequency is 77.0 Hz\n",
      "Oscillator frequency is 77.0 Hz\n"
     ]
    }
   ],
   "source": [
    "# Output is directly connect to A input\n",
    "SR.reference('INT')\n",
    "SR.I_mode('CURRENT_MODE_OFF')\n",
    "SR.V_mode('A_INPUT_ONLY')\n",
    "SR.osc_amplitude(0)\n",
    "SR.osc_frequency(77)\n",
    "print('Reference frequency is ' + str(SR.frequency()) + ' Hz' )\n",
    "print('Oscillator frequency is ' + str(SR.osc_frequency()) + ' Hz')\n",
    "SR.sensitivity(.1)\n",
    "SR.timeconstant(.1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:40.426509Z",
     "start_time": "2022-01-14T09:36:40.414509Z"
    }
   },
   "outputs": [],
   "source": [
    "#Set up lockin for polar measurement of output\n",
    "In_R = qc.ScaledParameter(SR.r, gain=1, name='R', unit='V')\n",
    "In_phase = qc.ScaledParameter(SR.phase, gain=1, name='phase', unit='degrees')\n",
    "\n",
    "output = qc.ScaledParameter(SR.osc_amplitude, gain=1, name='output', unit='V')  "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:54.188509Z",
     "start_time": "2022-01-14T09:36:41.692509Z"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Starting experimental run with id: 1. Using 'qcodes.utils.dataset.doNd.do1d'\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "(output_varied #1@C:\\Users\\g-gre-gre050342\\Qcodes_contrib_drivers\\docs\\examples\\test_signal_recovery_7270.db\n",
       " -----------------------------------------------------------------------------------------------------------\n",
       " output - numeric\n",
       " R - numeric\n",
       " phase - numeric,\n",
       " [<AxesSubplot:title={'center':'Run #1, Experiment test_signal_recovery_7270 (output2input)'}, xlabel='output (mV)', ylabel='R (mV)'>,\n",
       "  <AxesSubplot:title={'center':'Run #1, Experiment test_signal_recovery_7270 (output2input)'}, xlabel='output (mV)', ylabel='phase (degrees)'>],\n",
       " [None, None])"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZoAAAEWCAYAAABfdFHAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAAzMElEQVR4nO3dd5wV1fnH8c9D772XpSPSLCDYu0axgEqMHSsxTX+JERALKjZMTNREY4jdWKJUFTuiYFQULLtLFZDee2fb8/tjZuN13WULe3d27/2+Xy9e7L3Tnrkzd545Z849x9wdERGReKkUdQAiIpLYlGhERCSulGhERCSulGhERCSulGhERCSulGhERCSulGgSgJmlmNlOM6scdSzlnZm9bWZDymA7z5rZPfHejhSfmXU3s1lmZlHHcqDM7Akzuz2ibT9kZr8qyrxlmmjMbKmZ7QkvimvDL2OdMtjuKjOraWYnm9mEPNNGm1mamWWZ2Z3FXO9HZrY33J/cf2+UavBF4O7L3b2Ou2eX9bbNzM2s836mX2lmn5TStpaa2akHsg53P9PdnyuNeKTozGxOnu9JVu53xcy6mtlkM9tgZpvN7F0zOyhm2SfyLLvPzHbETG9kZhPNbJeZLTOzSwoJZzTwZy+FHxEWdv4f6PrM7Egzez/8XDaY2Wtm1jJ3urtf7+6jS2v7+4krv+/xn4GRZlatsOWjKNGc4+51gEOBw4Bb4rkxM2sLbHL3PUAf4Ks8sywChgFTSriJ34YX+dx/5xxAuMVmZlXKcnsSP+X9WB5IfO7eI/c7AtQFVgCvhZMbAK8DBwHNgS+AyTHLXh/7HQNejlkW4DEgI1z2UuAfZtajgH1oCZwETCrpvpSxhsBYoD3QDtgBPBNlQLncfQ0wHzi3KDOX2T9gKXBqzOsHgSnh3ycCKwuaH7gTeBV4nuDDngP0LcI2zwOeCf/+DzCggPn+DdxZzP35CLi2gGnDgZlAlfD1r8KYaxCcNA4MBVYDa4A/xixbCRgBLAY2hfvdKJyWu+w1wHJgesx7VWLiugf4FNgJvAE0Bl4EtgNfAu1jttcNeB/YDCwALoyZ9izBF3lK+LnPBDqF06aH290VbucXeT6Dg4G9QHY4fWv4fnWCu6HlwDrgCaBmOK0J8CawNYxnRvh5vADkAHvCdQ3bz3GpER7PTeF6vgSa5z1mQGXgIWAj8D3w23w+x9HAf8N9fw9oErOd14C1wLbws+iR53O7p5Dz50RgJcG5sjbcxwKPfbjMseFx3Upwsb4yfL8+wXdjA7AMuC1cV/Vw3p4x62gafo7NwtdnA9+E830K9M7zHRwOpAL7gJuB8Xn241HgkWJ8b04IP8/aBUxvFB6HxvlMqx0ue0LM6wyga8w8LwAPFLDuK4AP8jlPPwr3fw5wbkHfceBK4JOCzv+YYzoyPK+WApeWdH35xH84sCO/8yxm2zcB6wmuK1flmfcJgu/6DuBjoF2e60qVvLFSwPc4nOdWwuvr/v5F9ozGzNoAZxKUKIrqXOAVfrgD+vt+1j/KzLYS3P38Ivx7MPCymW0tg+cZfyL4Yt5mZl2A+4DL3H1vzDwnAV2A04HhMdVCvwMGEXwhWwFbCC72sU4gOAF+VsD2LwIuB1oDnYDPCO6EGgHzgFEAZlab4MR7CWgWLve4mXXPs667CO6uFgH3Arj78eH0Qzy42/xPbADuPg+4HvgsnN4gnPQA0JWgVNs5jPGOcNpNBF+WpgR3qCODVfnlBInpnHBdDxaw3wBDCC68bQkS7PUEF9a8riM4Bw8l+AIPymeeS4CrCD6basAfY6a9TXD8mhGUlF/cT0wFaUFwTNoR3HgUeOzNrF24zb8RfD6HEiQIwvfqAx3DZa8guMjsAyYAF8ds80LgY3dfb2aHAU8DvyT4rP4JvG5m1WPmvxg4i+B792/gDDNrEMZUheD8eL4Y+zyEIFntKmD68cBad9+Uz7QLCJLp9PB1VyDL3RfGzPMtkG+JBuhFcDMFgJlVJbgRe4/gOP4OeDG26q4g+zn/WxDcMLUm2NexB7i+WMcTJMOCtCA4D1oT3Iw+ZmYNY6ZfSnDz1ITg3Cn0nN3P9xiCa8khha0jikQzKaxfXUGQdUcVY9lP3P0tD55FvMB+dtDd7yL4ML8HOgADgHfcvb67N/DSe57xaJi4cv+NDrefQ/Blv4EgKT7o7l/nWfYud9/l7mkESSD3YnA9cKu7rwwvFHcCg/NUXdwZLpvfBRSCu4zF7r6N4OK02N0/cPcsgjvxw8L5zgaWuvsz7p4Vxjge+HnMuia6+xfhsi8SXOBKJHwAOxT4vbtvdvcdBEn4onCWTKAlwZ1WprvP8PDWqRgyCS6and09291nu/v2fOa7kOBOfKW7byFIgHk94+4Lw8/5VWL23d2fdvcdMcfoEDOrX8xYc4BR7r4v3Mb+jv0lBHfjL4efzSZ3/ya8aboIuCWMZylBSe3ycBsv8cPnS7iel8K/hwL/dPeZ4Wf1HMEN0pEx8z/q7ivcfY8H1SXT+eH8OAPY6O6zi7KzZlaL4Ibv2QKmtyFIrH8oYBVDgOdjzok6BKX0WNsIqufy04Dgbj7XkeE6HnD3DHf/kKBEfXE+yxbH7eEx/ZigNuDCA1wfZtab4Ibs5v3MlgncHZ4fbxGUQGKT3BR3nx6eW7cCR4WPF0pqB8Fnul9RJJpB7l6XoJjXjSAZFNXamL93AzXyqzc2s0PDEswWgjvmBcA04MQwGZxfwtjzc0OYuHL//a8FSPiFn0ZQLM1bIoEg2eZaRnAHC8Hd7cTc5EVw15BNcIef37L5WRfz9558Xuc2wmgH9I9NlgR3PS1i5s/7uR9IA46mQC1gdsz23gnfh6AkuAh4z8yWmNmIEmzjBeBd4BUzW21mD4Z3rnm14sefY36fab77bmaVzewBM1tsZtsJqkigeOczwIY8pdz9Hfu2BFVqeTUBqhKcQ7mWEdzVQnAO1jKz/mbWniBZTozZ3k15jn9bfjgX4aefy3PAZeHflxF83kV1PkGV6Md5J5hZU4KSxePu/nI+01MIrhuxpaedQL08s9bjx8kk1hZ+nIRaASvCG8NcsZ9dSWzJU1qL/W6XSNhA4G3gRnefsZ9ZN4U3hLnyfl//dyzdfSfBsTiQ2OoSVDnuV2RVZ2Gmf5agrh6CeslaudPDu7SmP12ySOv+Jize3QvcEf49l6BI2sDdJ+xv+dJiZmcBRwFTCS6gecXeSaQQPK+B4GQ4M08Cq+Huq2LmL61ut1cQVKPEbquOuxep2WIR5I1zI0Gi6xGzvfoePOQlvCO/yd07ElSV/sHMTilgXflvMLibu8vduwNHE5Tarshn1jVAm5jXxbmzuwQYCJxKUFXRPny/uE1m8+7T/o79CoJq0Lw2EtzJtot5LwVYBRCW3l8luEu/GHgzLEnmbu/ePNurledCnzfGSUBvM+tJ8NkWp8owb4kEgLB65z3gdXe/t4BlLwf+6+5LYt5bCFQJq6dzHULB1UupBNVtuVYDbc0s9lr4v8+OPNclfnwDVpCGYZV07Ppyv9vFXl9YZfoBMNrdi5PU8/O/c9yCFr+NwthyE2NBsRX03TuYoKpyv6L+Hc3DwGlmdgjBCVPDzM4K7z5vI3iQeSD6AF+Fze9auftPngeZWVUzq0HwWVQxsxq5z2/MrH3Y3LB9cTdsZk2AJwkepg0BzjGzAXlmu93MaoUtZK4iaKwAwQO7e8MTDDNramYDixtDEb0JdDWzy8PPoqqZHWFmBxdx+XUEzwX2N71NbhPI8M7xX8BfzawZgJm1NrOfhX+fbWadwyq2bQR38zkx69rftgjXcZKZ9QqP43aCi3BOPrO+CtwYbr8BwUPvoqpLUMW0ieDLeV8xlt2f/R37F4FTzexCM6tiZo3N7NCYRHKvmdUNl/0DwfOUXC8RPKy+lB+qzSA4FteHpR0zs9rhd7CgqifCEti4cD1fuPvyouxYWC12EkGJKPb9egQl0P+6+/5KsFeQp8otLDlMAO4OYz+G4AagoAvy+8Dh4XcegsYtu4Fh4bl/InAOwbNgCJ5jnB9+TzsTPPeIVdA5eZeZVTOz4wiScW4ruWKtz8xaAx8Cf3f3JwrYp+IYYGbHht/H0cDnYbXoBoLkellYWr+aH9/U/Oh7HOMEgpLWfkWaaMKde56g1LEN+DXBxXkVQYZdeYCbyG3O3AtIL2CefxHcYV9MUGe5hx/qttsSFHtX5b8oAH+3H7fxz62rHgtM9uCZ0iaCE+pJM2scs+zHBNVEUwna9b8Xvv8IwXOd9yx4nvU50L+oO10c4Z3t6QR1+KsJqorGUPQkfyfwXFjtkl899IcEd5drzWxj+N5wgv3+PKx2+oAf6pG7hK93EjRgeNzdp4XT7idoXLHVzGIfyufVguBCuJ2g6ulj8r/w/IvgLjoV+Bp4C8giSG6FeZ4fzo25BMeoNBR47MML+gCCBhObCS5auc8pf0fwnVkCfEKQBJ7OXam7zwyntyLmwuDuswgaRfydoFppEUFLqMI8R/C9Ks4d9uUED5TzVv+dBxwBXJXnu5SSO4OZHUVQ+nyNn/o1UJPgme/LwK/cPd8SjbuvIzgnB4avMwgSy5kEJcPHgSvcfX64yF8JWrWtC/c5b+ntTn56/q8l+CxXh/NffwDru5Yg8dwZ+9nkt29F9BLBc/HNBNfHy2KmXUfw/GcTQWOKT2Om/eR7bEFT8e4Uoam4Ff85a/Iws9sI6tD/WcrrbU/QSKFqnvpUiZCZnQk84e7tCp05yYVJYD7QooCGFuWWBS0qnwP65a3CK4V1nwj8293bFDJrmTOzZwl+QnJbKa3vIYJGRo8XNm+5/oFY1NxdXYgkMDOrSVCV8x7Bw/ZR/PCQXAoQPs/4A/BKRUsyAO4+l6AEJQfA3W8q6rxRP6MRKREzuzRPNUvuv/39xuAnqyH4fdAWgqqzefzwe57SinNkAXEWWq9dHlnwkHs7cBp5fppQwH7uDJ9TSBJT1ZmIiMSVSjQiIhJXFfoZTZMmTbx9+/ZRhyEiUqHMnj17o7uX6HeKJVGhE0379u2ZNWtW1GGIiFQoZras8LlKj6rOREQkrpRoREQkrpRoREQkrpRoREQkrpRoREQkruKWaMzsaTNbb2bpMe81MrP3zey78P+G4ftmZo+a2SIzSzWzw+MVl4iIlK14lmieJRh9L9YIYKq7dyHosTi3S/AzCXrt7UIw4t8/4hiXiIiUobglGnefTtAVdayB/DAWxXP8MEb7QMLBkNz9c6BB2AW1iIjE2J2Rxf1vz2Pllt1Rh1JkZf2Mpnk45jgEYzbkDk3cmh8PF7uSAoZSNbOhZjbLzGZt2LAhfpGKiJQzHy1Yz2l/mc4/P17CtAUV5/oXWc8A7u5mVuwePd19LMGgYvTt21c9gopIwtu0cx93vzmXyd+splPT2rz6y6Po16FR1GEVWVknmnVm1tLd14RVY+vD91fx4/Ha27D/US1FRBKeuzP+q1XcM2Uuu/ZlccMpXfjNSZ2oXqVy1KEVS1knmteBIcAD4f+TY97/rZm9QjBs7baYKjYRkaSzbNMubp2YzieLNtKnXUPuP78XXZvXjTqsEolbojGzl4ETgSZmtpJgkKQHgFfN7BqC8dZzx9h+i2As9EXAbuCqeMUlIlKeZWXn8OQn3/PwBwupUqkSowf15NJ+KVSqZFGHVmJxSzTufnEBk07JZ14HfhOvWEREKoK0ldsYPj6VuWu2c3r35tw9sCct6teIOqwDVqGHCRARSQS7M7L4y3sLefq/39OkTnWeuOxwzuiZOL/wUKIREYnQRwvWc+vEdFZt3cMl/VMYfkY36tesGnVYpUqJRkQkAhW9yXJxKNGIiJShRGmyXBxKNCIiZSSRmiwXhxKNiEic/aTJ8sAeXNq/XYVuslwcSjQiInGUunIrI8anMXfNdk7r3py7B/agZf2aUYdVppRoRETiYHdGFg+9t5BnYpos/6xHC8ySoxQTS4lGRKSUJUOT5eJQohERKSXJ1GS5OJRoREQOUDI2WS4OJRoRkQOQrE2Wi0OJRkSkBJK9yXJxKNGIiBTTnNXbGDYulTmrk7fJcnEo0YiIFNG+rGwe+3ARj3+0mAa1qvL4pYdzZs/kbLJcHEo0IiJF8PXyLQwbl8p363dy/mGtuf3s7jSsXS3qsCoEJRoRkf3Yk5HNXz9YyJMzltC8Xg2eufIITurWLOqwKhQlGhGRAsxcsonh41NZumk3F/dL4ZYB3ahXI3l/eFlSSjQiInns3JfFmLfn88Lny2jbqCYvXdufozs3iTqsCkuJRkQkxvSFG7hlQhqrt+3h6mM68MefdaVWNV0qD4Q+PRERYNvuTO6ZMpfXZq+kU9PajLv+KPq0U/cxpUGJRkSS3ntz1nLrpHQ278rg1yd24oZTulCjqrqPKS1KNCKStDbt3Meo1+fwZuoaurWoyzNXHkHP1vWjDivhKNGISNJxd95IXcOdr89hx95M/nBaV64/oRPVqlSKOrSEpEQjIkll3fa93DoxnQ/mreOQNvV5cPCRHNRCnWDGkxKNiCQFd+e1WSsZPWUuGVk5jBzQjauP6UCVyirFxJsSjYgkvJVbdnPLhDRmfLeRfh0aMeaC3nRoUjvqsJKGEo2IJKycHOffM5cx5u35AOrKPyJKNCKSkL7fuIvh41L5YulmjuvShPvP70WbhrWiDispKdGISELJznGe+mQJD723kOpVKvHg4N78vE8bdeUfoUgSjZn9HrgWcCANuApoCbwCNAZmA5e7e0YU8YlIxbRw3Q5uHpfKtyu2clr35twzqCfN69WIOqykV+bNLcysNXAD0NfdewKVgYuAMcBf3b0zsAW4pqxjE5GKKTM7h0enfsdZj85gxebd/O3iwxh7eR8lmXIiqqqzKkBNM8sEagFrgJOBS8LpzwF3Av+IJDoRqTDSV23j5nGpzFuznXMPacWoc7rTuE71qMOSGGWeaNx9lZn9GVgO7AHeI6gq2+ruWeFsK4HW+S1vZkOBoQApKSnxD1hEyqW9mdk8OvU7/jl9CY1rV+NfV/TltO7Now5L8lHmicbMGgIDgQ7AVuA14IyiLu/uY4GxAH379vU4hCgi5dzsZZsZNi6VxRt2cWHfNtw6oDv1a2lAsvIqiqqzU4Hv3X0DgJlNAI4BGphZlbBU0wZYFUFsIlKO7c7I4k/vLuDZT5fSqn5Nnr+6H8d3bRp1WFKIKBLNcuBIM6tFUHV2CjALmAYMJmh5NgSYHEFsIlJOfbpoI8MnpLJi8x6uOKodw87oRp3q+oVGRRDFM5qZZjYO+ArIAr4mqAqbArxiZveE7z1V1rGJSPmzfW8m9781n5e/WE77xrX4z9Aj6d+xcdRhSTFEcjvg7qOAUXneXgL0iyAcESmnPpy/jpET0lm/Yy+/PL4jvz+tqwYkq4BU7hSRcmfLrgxGvzmXCV+vomvzOjxx+TEc2rZB1GFJCSnRiEi58nbaGm6fnM7W3ZnccEoXfnNSJ6pXUSmmIlOiEZFyYcOOfdwxOZ2309fSs3U9nr+6P91b1Ys6LCkFSjQiEil3Z9I3q7jrjbnszshm2BkHMfS4jhqQLIEo0YhIZFZv3cOtE9OYtmADfdo1ZMwFvencrE7UYUkpU6IRkTLn7rz8xQrue2se2TnOqHO6c8VR7amsAckSkhKNiJSp5Zt2M3x8Kp8t2cTRnRrzwPm9SWmsAckSmRKNiJSJ7BznuU+X8qd3F1ClknH/+b246Ii2GpAsCSjRiEjcLVq/k2HjvuWr5Vs5uVsz7j2vJy3r14w6LCkjSjQiEjdZ2Tn8c/oSHpn6HbWqVebhXxzKwENbqRSTZJRoRCQu5q7ezrDx35K+ajsDerXgrnN70rSuBiRLRko0IlKq9mVl89iHi3j8o8U0qFWNf1x6OGf2ahl1WBIhJRoRKTXfrNjKsHHfsnDdTs4/vDV3nN2dBrWqRR2WREyJRkQO2N7MbP7y/kKenLGE5vVq8MyVR3BSt2ZRhyXlhBKNiByQmUs2MXx8Kks37eaS/inccmY36tbQsMryAyUaESmRXfuyGPPOfJ7/bBkpjWrx0nX9ObpTk6jDknJIiUZEim3GdxsYMT6N1dv2cPUxHfjjz7pSq5ouJ5I/nRkiUmTb9mRy35R5/GfWCjo2rc2464+iT7tGUYcl5ZwSjYgUydR56xg5MY2NOzP41YmduPGULhpWWYpEiUZE9mvzrgzufmMOk75ZTbcWdfnXFX3p3aZB1GFJBaJEIyIFmpK6hjsmp7NtTyb/d2oXfn1iZ6pV0YBkUjxKNCLyE+t37GXU5Dm8nb6WXq3r8+9r+3NwSw2rLCWjRCMi/+PuTPw6GFZ5T2Y2I87sxrXHdtCwynJAlGhEBIA12/YwcsIPwyo/OLg3nZpqWGU5cEo0Ikkud1jl+9+aR5aGVZY4UKIRSWLLN+1mxIRUPl28iaM6NmbMBRpWWUqfEo1IEsrJcZ7/bClj3llA5UrGfef14uJ+GlZZ4kOJRiTJLN6wk+HjUpm1bAsnHtSU+87rRasGGlZZ4keJRiRJZGXn8OQn3/OX9xdSs2pl/nLhIZx3WGuVYiTulGhEksD8tdsZNi6V1JXb+FmP5owe1JNmdWtEHZYkCSUakQSWkZXDPz5azN+nfUe9GlX5+yWHcVavlirFSJmKJNGYWQPgSaAn4MDVwALgP0B7YClwobtviSI+kUSQtnIbN4/7lvlrdzDw0FaMOqcHjWprWGUpe1H93PcR4B137wYcAswDRgBT3b0LMDV8LSLFtDczmwffmc+gx//L5l0Z/OuKvjxy0WFKMhKZMi/RmFl94HjgSgB3zwAyzGwgcGI423PAR8Dwso5PpCKbvWwLw8Z9y+INu7iwbxtuPas79WtqWGWJVhRVZx2ADcAzZnYIMBu4EWju7mvCedYCzfNb2MyGAkMBUlJS4h+tSAWwOyOLP7+7kGc+/Z5W9Wvy/NX9OL5r06jDEgGiSTRVgMOB37n7TDN7hDzVZO7uZub5LezuY4GxAH379s13HpFk8unijYwYn8byzbu5/Mh2DD+zG3Wqq52PlB9RnI0rgZXuPjN8PY4g0awzs5buvsbMWgLrI4hNpMLYsTeT+9+ez0szl9O+cS1eGXokR3ZsHHVYIj9R5onG3dea2QozO8jdFwCnAHPDf0OAB8L/J5d1bCIVxbQF6xk5IY112/dy3XEd+MNpB1GzmoZVlvIpqvL174AXzawasAS4iqAF3Ktmdg2wDLgwothEyq2tuzO4+825TPhqFV2a1eHxXx3NYSkNow5LZL+KlGjMrBlwDNAK2AOkA7PcPackG3X3b4C++Uw6pSTrE0kG76Sv5bZJ6WzZncHvTu7Mb0/uTPUqKsVI+bffRGNmJxE8P2kEfE3w3KQGMAjoZGbjgIfcfXuc4xRJWht37mPU5DlMSVtD95b1ePaqI+jZun7UYYkUWWElmgHAde6+PO8EM6sCnA2cBoyPQ2wiSc3def3b1dz5+hx27cvm5p8dxNDjO1JVwypLBVNYonnI3dfmN8Hds4BJpR6RiLB2215um5TGB/PWc2jbBvxpcG+6NK8bdVgiJVJYovnGzNKBl4Hx7r41/iGJJC9359VZK7hnyjwys3O47ayDueqYDhpWWSq0whJNa+BU4CLgPjP7nCDpTHb3PfEOTiSZrNi8m1smpPHJoo3079CIMRf0pn2T2lGHJXLA9pto3D0beBd4N2yKfCZB0nnYzKa6+6VlEKNIQsvJcV74fBlj3pmPAfcM6skl/VKopFKMJIgi/47G3TPMbC5BT8t9gIPjFpVIkliyYSfDx6fy5dItHN+1Kfef34vWGlZZEkyhicbM2hKUYi4GahNUnZ3r7vPjHJtIwsrKzuGpcFjl6lUq8afBvRncp40GJJOEVNjvaD4leE7zKkEz59llEpVIAluwdgfDxn3Ltyu3cVr35tw7qCfN6mlYZUlchZVoRgAz3F29JIscoMzsYFjlv334HXVrVOVvFx/G2b01rLIkvsIaA0wHMLMOBP2TtY9dxt3PjWdwIokifdU2bh6Xyrw12znnkFbceU53GtepHnVYImWiqI0BJgFPAW8AJerfTCQZ7c3M5tGp3/HP6UtoXLsaYy/vw+k9WkQdlkiZKmqi2evuj8Y1EpEEM3vZZoaNS/1hWOUB3alfS8MqS/IpaqJ5xMxGAe8B+3LfdPev4hKVSAWmYZVFfqyoiaYXcDlwMj9UnXn4WkRC0xas57aJ6azaukfDKouEivoN+DnQ0d0z4hmMSEW1Ycc+7n5zLm98u5rOzerw2vVHcUT7RlGHJVIuFDXRpAMNCMajEZFQbieY906Zx97MHH5/aleuP7GjBiQTiVHURNMAmG9mX/LjZzRq3ixJa/GGnYyckMbM7zfTr0Mj7juvF52b1Yk6LJFyp6iJZlRcoxCpQPZlZfPER0t4bNoialStxAPn9+LCvm3VCaZIAQrrgsY88HFh85R+aCLlz5dLN3PLhDQWrd/J2b1bcsc53WlWV93HiOxPYSWaaWY2nmD8mf8N5xwOGXAsMASYBjwbtwhFyoFtezIZ8858Xpq5nNYNavLMlUdwUrdmUYclUiEUlmjOAK4GXg67odkK1AQqEfym5mF3/zquEYpEyN15O30to16fw6ad+7j22A78/rSu1FaTZZEiK6yvs73A48DjZlYVaALs0ZDOkgxWb93DHZPT+WDeenq0qsfTQ46gV5v6UYclUuEUZ+CzTGBNHGMRKReyc5znPl3KQ+8tIMfh1gEHc9Ux7alSuVLUoYlUSCr/i8SYs3obIyek8e3KbZzQtSn3DOpJ20a1og5LpEJTohEB9mRk8/DUhTw543sa1KzKIxcdyrmHtNJYMSKloESJxswqARe7+4ulHI9ImZu+cAO3TkpjxeY9/KJvW24Z0I0GtapFHZZIwijsdzT1gN8QDOf8OvA+8FvgJuBbQIlGKqxNO/dxz5R5TPx6FR2b1Obl647kqE6Now5LJOEUVqJ5AdgCfAZcC4wEDBjk7t/ENzSR+HB3xn+1inumzGXXvixuOLkzvz6pMzWqqn8ykXgoLNF0dPdeAGb2JEGrs5Sw2bNIhfP9xl3cOjGNTxdvok+7htx/fi+6Nq8bdVgiCa2wRJOZ+4e7Z5vZytJKMmZWGZgFrHL3s8MfhL4CNAZmA5drWAIpLRlZOfxrxhIemfod1StX4p5BPbmkX4r6JxMpA4UlmkPMbHv4twE1w9cGuLvXO4Bt3wjMA3LXMQb4q7u/YmZPANcA/ziA9YsAMHvZFkZOSGPBuh2c2bMFd57bg+b11D+ZSFkprGeAuFRam1kb4CzgXuAPFrQhPRm4JJzlOeBOlGjkAOzYm8mf3l3AC58vo3ndGoy9vA+n92gRdVgiSSeq39E8DAwDcivHGwNb3T0rfL2SoKWbSIl8OH8dIyeks27HXoYc1Z6bTu9K3RpVow5LJCmVeaIxs7OB9e4+28xOLMHyQ4GhACkpKaUbnFR4m3dlcPcbc5j0zWq6NKvD45cdzeEpDaMOSySpRVGiOQY418wGADUIntE8AjQwsyphqaYNsCq/hd19LDAWoG/fvhoHR4CgyfKUtDWMmjyHbXsyueGULvzmpE4aUlmkHCjzROPutwC3AIQlmj+6+6Vm9howmKDl2RBgclnHJhXTuu17uX1SOu/NXUev1vX597X9ObjlgbRTEZHSVJ76OhsOvGJm9wBfA09FHI+Uc+7Oa7NWMnrKXDKychhxZjeuPbaDelkWKWciTTTu/hHwUfj3EqBflPFIxbFi825GTkxjxncbOaJ9Q8Zc0JuOTetEHZaI5KM8lWhECpWT4zz/2VIefHcBBtw9sAeX9W+nH16KlGNKNFJhLN6wk+HjUpm1bAvHd23Kfef1pE1DjRUjUt4p0Ui5l5Wdw9gZS3j4g++oUaUSf/75IVxweGuNFSNSQSjRSLk2Z/U2ho9PJX3Vds7o0YK7B/WgWV11HyNSkSjRSLm0Lyubv01dxBMfL6ZBrao8funhDOjVMuqwRKQElGik3Plq+RaGjUtl0fqdnH94a24/qzsNa2vES5GKSolGyo3dGVn8+d2FPPPp97SoV4NnrjqCkw5qFnVYInKAlGikXPh00UZGTEhj+ebdXHZkCsPP6KZOMEUShBKNRGr73kzuf2seL3+xgvaNa/HK0CM5smPjqMMSkVKkRCORmTpvHbdOTGf9jr0MPb4jvz+1KzWrqRNMkUSjRCNlbvOuDO56Yw6Tv1nNQc3r8sTlfTi0bYOowxKROFGikTLj7ryRuoY7X5/Djr2Z3HhKF35zUmeqVVEnmCKJTIlGysS67Xu5dWI6H8xbR+829XlwcH+6tVBX/iLJQIlG4ipvV/4jB3Tj6mPUlb9IMlGikbhZsXk3t0xI45NFG+nXoRFjLuhNhya1ow5LRMqYEo2UuuzcrvzfWUAlg9GDenJpvxR15S+SpJRopFQtWr+D4ePTmL1sCyd0bcp95/eidYOaUYclIhFSopFSkZmdw9jpS3jkg++oVb0yf7nwEM47TF35i4gSjZSC9FXbGDYulblrtjOgVwvuOrcnTetWjzosESknlGikxPZmZvPo1O/45/QlNKpdjScuO5wzeqorfxH5MSUaKZFZSzczbHwqSzbs4ud92nDbWd2pX0udYIrITynRSLHs2pfFn95dwHOfLaVV/Zo8f3U/ju/aNOqwRKQcU6KRIpvx3QZGjE9j9bY9DDmqPTf/7CBqV9cpJCL7p6uEFGrb7kzumTKX12avpGPT2rz6y6M4on2jqMMSkQpCiUb26905a7ltUjqbd2XwqxM7ceMpXahRVV35i0jRKdFIvjbs2Medr89hStoaDm5Zj2euPIKeretHHZaIVEBKNPIj7s6kb1Zx1xtz2b0vmz+e3pVfntCJquoEU0RKSIlG/mf11j3cOjGNaQs2cHhKAx4c3JvOzepGHZaIVHBKNEJOjvPSF8t54O35ZOc4d5zdnSFHt6eyOsEUkVKgRJPklm7cxfDxqcz8fjPHdG7MA+f3pm2jWlGHJSIJRIkmSWVl5/D0f7/nofcWUq1KJcZc0IsL+7ZVJ5giUurKPNGYWVvgeaA54MBYd3/EzBoB/wHaA0uBC919S1nHlwzmr93O8HGpfLtyG6ce3Jx7z+tJ83o1og5LRBJUFCWaLOAmd//KzOoCs83sfeBKYKq7P2BmI4ARwPAI4ktYGVk5PDZtEY9/tIh6Naryt4sP4+zeLVWKEZG4KvNE4+5rgDXh3zvMbB7QGhgInBjO9hzwEUo0pebbFVsZNi6VBet2MPDQVow6pweNaleLOiwRSQKRPqMxs/bAYcBMoHmYhADWElSt5bfMUGAoQEpKShlEWbHtycjmrx8s5MkZS2hWtwZPDenLKQfn+9GKiMRFZInGzOoA44H/c/ftsdU37u5m5vkt5+5jgbEAffv2zXceCXy+ZBMjxqeydNNuLu6Xwi0DulGvhrryF5GyFUmiMbOqBEnmRXefEL69zsxauvsaM2sJrI8itkSwY28mD7w9nxdnLielUS1euq4/R3dqEnVYIpKkomh1ZsBTwDx3/0vMpNeBIcAD4f+Tyzq2RDBtwXpunZDGmu17uebYDtx0eldqVVMrdhGJThRXoGOAy4E0M/smfG8kQYJ51cyuAZYBF0YQW4W1ZVcGo9+cy4SvV9GlWR3G/+poDk9pGHVYIiKRtDr7BCioPe0pZRlLongrbQ13TE5n6+5MfndyZ357cmeqV1FX/iJSPqhOpQJbv2Mvd0yawztz1tKzdT2ev7o/3VvVizosEZEfUaKpgNyd8V+tYvSbc9mTmc3wM7px3XEdqKKu/EWkHFKiqWBWbtnNyInpTF+4gb7tGjJmcG86Na0TdVgiIgVSoqkgcnKcf89cxpi35+PAXef24PIj21FJXfmLSDmnRFMBLNmwkxHj0/hi6WaO69KE+87rpa78RaTCUKIpx7Kyc3jyk+/56/sLqV6lEg8O7s3P+7RRJ5giUqEo0ZRT89ZsZ9i4VNJWbeP07s25Z1BPmqkrfxGpgJRoypl9Wdk8Nm0xj09bRP2aVXnsksMZ0KuFSjEiUmEp0ZQjXy/fwvDxqSxct5PzDmvN7Wd3V1f+IlLhKdGUA3sysnnovQU8/d/vaV6vBk9f2ZeTu6krfxFJDEo0Efts8SZGTEhl2abdXNI/hVvO7EZddeUvIglEiSYiO/Zmcv/b83lp5nLaNVZX/iKSuJRoIjBt/npGTkxj3fa9XHtsB246/SBqVlMnmCKSmJRoytCWXRnc/eZcJoZd+T/+q6M5TF35i0iCU6IpA+7OlLQ1jJo8h217Mrnh5M78Rl35i0iSUKKJs3Xb93L7pHTem7uOXq3r8+9r+3NwS3XlLyLJQ4kmTtyd12atZPSUuWRk5XDLmd245lh15S8iyUeJJg5WbN7NLRPS+GTRRvq1b8QDF/Sio7ryF5EkpURTirJznOc/W8qD7yygksHoQT25tF+KuvIXkaSmRFNKFq3fwbBxqXy1fCsndG3Kfef3onWDmlGHJSISOSWaA5SZncM/P17Mo1MXUat6Zf5y4SGcd1hrdYIpIhJSojkA6au2cfO4VOat2c5ZvVpy57k9aFq3etRhiYiUK0o0JbA3M5tHpn7H2OlLaFS7Gk9c1oczeraIOiwRkXJJiaaYvly6meHjUlmycRc/79OG287qTv1a6gRTRKQgSjRFtHNfFg++M5/nP1tGm4Y1eeGafhzXpWnUYYmIlHtKNEXw8cINjJyQxupte7jy6Pbc/LODqF1dH52ISFHoarkfW3dnMPrNeYz/aiWdmtbmtV8eRd/2jaIOS0SkQlGiKcDbaWu4ffIctuzO4DcndeJ3J3ehRlV1gikiUlxKNHms37GXUZPn8Hb6Wnq0qsdzVx9Bj1b1ow5LRKTCUqIJuTvjv1rF6Dfnsiczm2FnHMR1x3WkqjrBFBE5IEo0wMotuxk5MZ3pCzfQt11DHrigN52bqRNMEZHSUK4SjZmdATwCVAaedPcH4rm9nBznhc+XMead+QDcdW4PLj+ynTrBFBEpReUm0ZhZZeAx4DRgJfClmb3u7nPjsb3FG3YyYnwqXy7dwnFdmnDfeb1o26hWPDYlIpLUyk2iAfoBi9x9CYCZvQIMBEo90bz65Qpum5xOjSqV+NPg3gzu00adYIqIxEl5SjStgRUxr1cC/fPOZGZDgaEAKSkpJdpQh6a1OaVbM+4a2INmdWuUaB0iIlI05SnRFIm7jwXGAvTt29dLso4j2jfiCP3wUkSkTJSntrurgLYxr9uE74mISAVWnhLNl0AXM+tgZtWAi4DXI45JREQOULmpOnP3LDP7LfAuQfPmp919TsRhiYjIASo3iQbA3d8C3oo6DhERKT3lqepMREQSkBKNiIjElRKNiIjElRKNiIjElbmX6DeP5YKZbQCWlXDxJsDGUgynItA+Jwftc3I4kH1u5+5NSzOY/anQieZAmNksd+8bdRxlSfucHLTPyaEi7bOqzkREJK6UaEREJK6SOdGMjTqACGifk4P2OTlUmH1O2mc0IiJSNpK5RCMiImVAiUZEROIqKRONmZ1hZgvMbJGZjYg6nngws7ZmNs3M5prZHDO7MXy/kZm9b2bfhf83jDrW0mRmlc3sazN7M3zdwcxmhsf6P+EQFAnDzBqY2Tgzm29m88zsqCQ4xr8Pz+l0M3vZzGok2nE2s6fNbL2Zpce8l+9xtcCj4b6nmtnh0UWev6RLNGZWGXgMOBPoDlxsZt2jjSousoCb3L07cCTwm3A/RwBT3b0LMDV8nUhuBObFvB4D/NXdOwNbgGsiiSp+HgHecfduwCEE+56wx9jMWgM3AH3dvSfBkCIXkXjH+VngjDzvFXRczwS6hP+GAv8ooxiLLOkSDdAPWOTuS9w9A3gFGBhxTKXO3de4+1fh3zsILkCtCfb1uXC254BBkQQYB2bWBjgLeDJ8bcDJwLhwlkTb3/rA8cBTAO6e4e5bSeBjHKoC1DSzKkAtYA0JdpzdfTqwOc/bBR3XgcDzHvgcaGBmLcsk0CJKxkTTGlgR83pl+F7CMrP2wGHATKC5u68JJ60FmkcVVxw8DAwDcsLXjYGt7p4Vvk60Y90B2AA8E1YXPmlmtUngY+zuq4A/A8sJEsw2YDaJfZxzFXRcy/01LRkTTVIxszrAeOD/3H177DQP2rYnRPt2MzsbWO/us6OOpQxVAQ4H/uHuhwG7yFNNlkjHGCB8LjGQIMm2Amrz0yqmhFfRjmsyJppVQNuY123C9xKOmVUlSDIvuvuE8O11ucXq8P/1UcVXyo4BzjWzpQTVoScTPL9oEFaxQOId65XASnefGb4eR5B4EvUYA5wKfO/uG9w9E5hAcOwT+TjnKui4lvtrWjImmi+BLmErlWoEDxJfjzimUhc+n3gKmOfuf4mZ9DowJPx7CDC5rGOLB3e/xd3buHt7gmP6obtfCkwDBoezJcz+Arj7WmCFmR0UvnUKMJcEPcah5cCRZlYrPMdz9zlhj3OMgo7r68AVYeuzI4FtMVVs5UJS9gxgZgMI6vMrA0+7+73RRlT6zOxYYAaQxg/PLEYSPKd5FUghGGLhQnfP+9CxQjOzE4E/uvvZZtaRoITTCPgauMzd90UYXqkys0MJGj9UA5YAVxHcQCbsMTazu4BfELSs/Bq4luCZRMIcZzN7GTiRYCiAdcAoYBL5HNcw4f6doApxN3CVu8+KIOwCJWWiERGRspOMVWciIlKGlGhERCSulGhERCSulGhERCSulGhERCSulGhE8jCzK82s1QEs397MLtnP9Ja5vUsXcX1Dwuause81MbMNZlbdzF4xsy4ljVck3pRoRH7qSoLuTUqqPVBgogH+APyrGOubCJxmZrVi3hsMvBH+VuQfBH28iZRLSjSS8MzsD+HYJelm9n/he+3zjPXxRzO708wGA32BF83sGzOraWZLzexBM0szsy/MrHO4zLPh/Lnr2Bn++QBwXLj87/MJ6QLgnXCZK81sUji+yFIz+20Y79dm9rmZNQr7qPsYOCdmHRcBuaWcGcCpMV2wiJQrSjSS0MysD8Gv5fsTjMtznZkdVtD87j4OmAVc6u6HuvuecNI2d+9F8AvshwvZ7AhgRrj8X/PE0wHYkudX6z2B84EjgHuB3WEnmZ8BV4TzvEyQXAir9boCH4Yx5wCLCMajESl3lGgk0R0LTHT3Xe6+k6ATxuNKsJ6XY/4/6gDiaUnQtX+sae6+w903EHR7/0b4fhpBNRzAFOAYM6sHXAiMd/fsmHWs58Cq+0TiRolGklUWPz7/axQyv+fz9//WYWaVCPobK8yefLYVW7rJiXmdQzAUAGHJ6h3gPH5cbZarRrhukXJHiUYS3QxgUNjbb22CC/UMgo4Km5lZYzOrDpwds8wOoG6e9fwi5v/Pwr+XAn3Cv88Fqu5n+VwL+aGUUlwvEzQkaB4TQ66uQPpPlhApB5RoJKGFw1k/C3xB0HP1k+7+dTiWyd3h++8D82MWexZ4IrcxQPheQzNLBW4Ech/w/ws4wcy+JahO2xW+nwpkm9m3eRsDuPsuYHFug4Jiep+geuw/HtMbrpk1B/aEwwaIlDvqvVmkEOFgan3dfWMpre88oI+731ZK6/s9sN3dnyqN9YmUNjWHFClj7j7RzBqX4iq3Ai+U4vpESpVKNCIiEld6RiMiInGlRCMiInGlRCMiInGlRCMiInGlRCMiInH1/63nBqePoh9vAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZcAAAEWCAYAAACqitpwAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAAApBklEQVR4nO3deZhcVZ3/8feHhBD2AEmALJCERY0LEaLouOG4ASKg46iAbCqMOzNuoOMIyuC4/nQcNgERUAQ3UFQUcENcUIIgO5ouIAGEVCdB0p1ASPL9/XFOhZuil+ruW1Wdrs/refrpuvu5deve7z3LPVcRgZmZWZk2aXcCzMxs7HFwMTOz0jm4mJlZ6RxczMysdA4uZmZWOgcXMzMrnYPLRkrSLpJ6JI1rd1pGO0k/lXR0C7ZzgaT/bvZ2bOgkzZW0QJLanZaRknS2pP9q07a/KOldjczb9OAi6V5Jq/KF8KF8Am7Vgu0+IGlzSf8s6bK6aadKulXSGkmnDHG9v5b0WN6f2t+PSk18AyJiUURsFRFrW71tSSFp9wGmHyPptyVt615JrxzJOiLigIi4sIz0WOMk3V53nqypnSuS9pT0Q0lVScskXSXpaYVlz65b9nFJKwrTt5d0uaReSfdJOnyQ5JwKfCFKeLBvsN//SNcn6QWSrsnfS1XSdyXtXJseEe+MiFPL2v4A6errPP4C8DFJEwZbvlU5l9dFxFbAPOC5wEebuTFJM4GlEbEK2Af4c90sC4GPAD8Z5ibemy/stb/XjSC5QyZpfCu3Z80z2o/lSNIXEc+snSPA1sBi4Lt58iTgCuBpwI7An4AfFpZ9Z/EcAy4pLAtwBrA6L3sEcJakZ/azDzsDLwd+MNx9abHtgHOAWcCuwArg6+1MUE1E/B24Czi4kZmb+gfcC7yyMPw54Cf5837A/f3ND5wCfAe4iPQF3w7Mb2Cbrwe+nj9/Gziwn/m+CZwyxP35NfCOfqadCPwRGJ+H35XTPJH0QwngeOBB4O/AhwrLbgKcBHQBS/N+b5+n1ZZ9O7AI+E1h3PhCuv4b+D3QA/wI2AG4GHgUuAGYVdje04FrgGXA3cCbCtMuIJ28P8nf+x+B3fK03+Tt9ubtvLnuO3gG8BiwNk9/JI/fjHTXswh4GDgb2DxPmwz8GHgkp+e6/H18A1gHrMrr+sgAx2ViPp5L83puAHasP2bAOOCLQDdwD/DePr7HU4Hf5X2/Gphc2M53gYeAf+Tv4pl139t/D/L72Q+4n/RbeSjvY7/HPi/z4nxcHyFdoI/J47clnRtV4D7g43ldm+V5n1VYx5T8PU7NwwcBN+f5fg88p+4cPBG4BXgc+DDw/br9+Arwv0M4b16Wv88t+5m+fT4OO/Qxbcu87MsKw6uBPQvzfAP4TD/rPgr4eR+/01/n/b8dOLi/cxw4Bvhtf7//wjH9WP5d3QscMdz19ZH+vYEVff3OCtv+ILCEdF05tm7es0nn+grgWmDXuuvK+Pq00s95nOf5T/L1daC/lta5SJoBHEDKOTTqYOBSnrzTOX2A9Z8s6RHSXc6b8+c3ApdIeqQF9ROfJ52MH5e0B/Bp4K0R8VhhnpcDewCvBk4sFPm8DziUdBJOA5aTLvBFLyMd9Nf0s/23AEcC04HdgD+Q7ni2B+4ETgaQtCXpx/YtYGpe7kxJc+vW9UnSXdRC4DSAiHhpnr5XpLvKbxcTEBF3Au8E/pCnT8qTPgPsScq97p7T+Ik87YOkE2QK6U70Y2lVcSQpGL0ur+tz/ew3wNGki+1MUlB9J+liWu840m9wHumkPbSPeQ4HjiV9NxOADxWm/ZR0/KaScsQXD5Cm/uxEOia7km42+j32knbN2/w/0vczjxQUyOO2BebkZY8iXVgeBy4DDits803AtRGxRNJzgfOBfyN9V18FrpC0WWH+w4DXks67bwL7S5qU0zSe9Pu4aAj7fDQpQPX2M/2lwEMRsbSPaf9CCqC/ycN7Amsi4q+Fef4C9JlzAZ5NuoECQNKmpJuvq0nH8X3AxcViuf4M8PvfiXSTNJ20r+eMcH1FLyUFwP7sRPodTCfdgJ4habvC9CNIN0yTSb+dQX+zA5zHkK4lew22jlblXHpIUTOAXwCTilG3j/mLOZefF6bNBVYNsr3xeed3BP6JnEvqZ97h5lxWku54an+nFqbPIt193wl8tG58AE8vjPsc8LX8+U7gFYVpOwNP5P2pLTunj/UV77j/szD9i8BPC8OvA27On98MXFe3X18FTi7c7ZxXmHYgcFdhOIDdB/iOjiHfmeVhke7MdiuMeyFwT/78KVKRyFPWSV3Od4Btvo26O/C6Y1bLufwS+LfCtFf28T1+vDD93cDP+tnmpLzstoXvrZGcy2pgYmHcQMf+o8DlfaxnXF7P3MK4fwN+XdivrsK03wFH5c9nUfjN5nF382TO4F7gbXXTfwoclz8fBNwxhHNmC1Lueb9+ps8AHgAO62f6Lyicp8BLSIGoOM9xtX3vY/lzKeRqassDmxTGXVLbBgPkNPr6/edjuoZCroyU+/yv4ayvLu3PIV1PXlIYt/53lre9ig1zH0uAFxTmvbQwbStSbmQmA+Rc+kpnYZ5XAZXBjnurci6HRsTWpC/i6aQI2qiHCp9XAhP7KgeWNC/nVJaT7ozvBn4F7JdzLW8YZtr78v6ImFT4W99yIyLuzdudxVNzHpCKNWruI92pQrqLvTyn9RHSBWctKUj2tWxfHi58XtXHcK0hxa7AvrVt5e0dQboDqqn/3kfSCGMK6QJzY2F7P8vjIeX4FgJXS6pIOmkY2/gGcBVwqaQHJX0u36HWm8aG32Nf32mf+y5pnKTPSOqS9CjpIgxD+z0DVGPD3OxAx34mqbis3mRgU9JvqOY+0t0rpN/gFpL2lTSLlOO5vLC9D9Yd/5k8+VuEp34vFwJvzZ/fSvq+G/UG0gXy2voJkqaQchBnRsQlfUzfhXTdKOaSeoBt6mbdhnQD25flpDqfmmnA4ohYVxhX/O6GY3lsmCsrntvDkiv5fwqcEBHXDTDr0ohYUxiuP1/XH8uI6CEdi5GkbWvSTfWAWlosFhHXkiLpF/KoXtJFB0gnL09ecIa67psjZd1OAz6RP99Bym5OiojLBlq+LJJeS7or/wXpollvZuHzLqT6F0g/gAPqgtbEiHigMH+UlMzFpCKS4ra2ioiGmhg2oD6d3aTg9szC9raNVFFLRKyIiA9GxBxSMegHJL2in3X1vcGIJyLikxExl5RjPYhUTFTv76Q75ZqZfczTn8OBQ0i5gm1JNxCQcmZDUb9PAx37xaQiznrdpNzNroVxu5ByAERqRfgdUvHWYcCPI6J28V0MnFa3vS3qLu71afwB8BxJzyJ9t0MpDjwauCjybW9NLrq5GrgiIk7rZ9kjgd9FRKUw7q/A+Fz0XLMX/Rcd3UIqSqt5EJgpqXj9W//dUXddYsObrv5sl4ubi+urndtDXl8uDv05KYc5lEDel/W/8dxSd/uctlow7C9t/Z17zyAVQw6oHc+5fBl4laS9SD+SiZJem+8yP06qjByJfYA/56Zy0yLiKfU7kjaVNJG0/+MlTazVx0ialZsGzhrqhiVNBs4jVYgdDbxO0oF1s/2XpC1yy5ZjSQ0OIFW6nZZ/VEiaIumQoaahQT8G9pR0ZP4uNpX0PEnPaHD5h0nl/ANNn1FrrpjvEM8FviRpKoCk6ZJekz8fJGl3SSJVlK8lVeQ3si3yOl4u6dn5OD5KuvCu62PW7wAn5O1PIlVcN2prUp3aUtIJ+ekhLDuQgY79xcArJb1J0nhJO0iaVwgep0naOi/7AVJRb823SEWgR+TPNecC78y5GknaMp+Dxbv7DeSc1vfyev4UEYsa2bFcz/pyUs6nOH4bUk7zdxExUE71KNINaTEtvaQ6pU/ltL+IFPT7uwhfA+ydz3lIDVRWAh/Jv/39SMXGl+bpNwNvyOfp7qR6jKL+fpOflDRB0ktIAbjWum1I65M0nVR8e3pEnN3PPg3FgZJenM/HU4HrI2JxRFRJAfWtOVf+Nja8kdngPC54GSlHNaCWB5e8QxeRchf/IJVpn0fayV5Sxe5I1JoePxu4rZ95ziXdSR9GavmwinSHBCnK38eTdzF9OV0btsG/MY8/B/hhRFwZqWLy7cB5knYoLHstqQjoF6R291fn8f9LarBwtVJ7/uuBfRvd6aHId7CvJlXKPkgqBvosjQf2U4ALc5HKm/qY/kvSXeRDkrrzuBNJ+319LlL6OakZKqQK8p+Tijv+QCoi+VWe9j+kBhKPSCpWrNfbiXTxe5RUrHQtfV9sziXdLd8C3ARcSSovb+R5oYt48rdxB+kYlaHfY58v4geSGj0sI12oapWp7yOdMxXgt6QL//m1lUbEH/P0aRQuBhGxgFRHcTqpyGghqXx9MBeSzquh3EkfSaoUri/aez3wPODYunNpl9oMkl5IymV+l6d6N7A5qX7hEuBdEdFnziUiHib9Jg/Jw6tJweQAUg7wTFJ91F15kS+R6rMezvtcn0s7haf+/h8ifZcP5vnfOYL1vYMUbE4pfjd97VuDvkVqzLOMdH18a2HacaTWgEtJDSJ+X5j2lPNYqVn3XBpo1q26nGrHk/RxUpn4V0te7yxS09dN68pHrY0kHQCcHRG7Djpzh8sX/ruAnSLi0XanZyiUWkJeCDy/vniuhHXvB3wzImYMMmvLSbqA1Gjq4yWt74ukhiJnDjbvqH6Aqx0iwt13jGGSNicV01xNqjA/mScruq0fuX7iA6SWRxtVYAGIiDtIOSUbgYj4YKPzum8x22hIOqKuCKX2N9AzAE9ZDen5neWkYrE7efJ5m7LS+bF+0jloOfVopFRR/SipCerJddP62s+eXO9gHczFYmZmVjrnXMzMrHQbXZ3L5MmTY9asWe1OhpnZRuXGG2/sjohhPUc4HBtdcJk1axYLFixodzLMzDYqku4bfK7yuFjMzMxK5+BiZmalc3AxM7PSObiYmVnpHFzMzKx0Di5mZlY6BxczMyvdRvecy8aoUu3hl3ctAaDY204U3sWz4fi+xvU9b1GtK5/R0KNPf0noK23Rx9x9z9fgCs3GqJ223ZzD991l8BlHAQeXFvjsz+7iqtsfHnxGGxYN9T2QZhupeTMnObjYkxYu6eGVz9iRL715r/XjVLgiav04CuMK0/u4eA4272i43qqfq35fY/vex9GwF2Y2HKMiuEjan/Q2vnHAeRHxmTYnqTRPrF3HfUtX8ppn7sTWEzdtd3LMzFqi7RX6+Z3nZ5BeOToXOCy/NW5MWLxsJWvWBXOmbNXupJiZtUzbgwvwfGBhRFTyu60vJb/reizoqvYCsNuULducEjOz1hkNwWU6sLgwfH8et56k4yUtkLSgWq22NHEjVan2ADjnYmYdZTQEl0FFxDkRMT8i5k+Z0rLXEZSiUu1l8labse3mrm8xs84xGoLLA8DMwvCMPG5M6Kr2MMdFYmbWYUZDcLkB2EPSbEkTgLcAV7Q5TaWpdPeym4vEzKzDtL0pckSskfRe4CpSU+TzI+L2NierFMt7V7Osd7Ur882s47Q9uABExJXAle1OR9kq3bXKfAcXM+sso6FYbMzqWlJrhuxiMTPrLA4uTdTV3cOEcZswY7st2p0UM7OWcnBpokq1l1mTt2DcJu4jy8w6i4NLE3VVe5gz2UViZtZ5HFya5Im161i0dCW7TXVlvpl1HgeXJllU67DSORcz60AOLk1SyR1WuhmymXUiB5cm6XKHlWbWwRxcmqRS7XGHlWbWsRxcmqRS7XW3L2bWsRxcmiT1huwiMTPrTA4uTbCsdzXLVz7hnIuZdSwHlyaovX3SfYqZWadycGkCN0M2s07n4NIEXVV3WGlmnc3BpQm63GGlmXU4B5cmqHT3uL7FzDqag0vJah1Wur7FzDqZg0vJah1WOudiZp3MwaVkXUvcp5iZmYNLySrdboZsZubgUrKuJT1M2XoztpnoDivNrHM5uJSs0t3LnMnOtZhZZ3NwKVml2sNuU13fYmadzcGlRLUOK51zMbNO5+BSovUdVjrnYmYdzsGlRLVXG+822cHFzDqbg0uJKtVeJozfhOnbbd7upJiZtZWDS4m6qj3M3mFLd1hpZh3PwaVElWqvH540M8PBpTRPrF3HomUr3aeYmRkOLqW5b2nqsNI5FzMzB5fS1Johu8NKMzMHl9J0Vd1hpZlZTUuCi6QjJN0i6VZJv5e0V2Ha/pLulrRQ0kmtSE8zVKrusNLMrKZVOZd7gJdFxLOBU4FzACSNA84ADgDmAodJmtuiNJWqq9rDbs61mJkBLQouEfH7iFieB68HZuTPzwcWRkQlIlYDlwKHtCJNZat097q+xcwsa0edy9uBn+bP04HFhWn353EbkHS8pAWSFlSr1RYkcWiW9a7mkZVPuBmymVnW0uAi6eWk4HLiUJaLiHMiYn5EzJ8yZUpzEjcCXetbirlYzMwMmhhcJL1H0s35b5qk5wDnAYdExNI82wPAzMJiM/K4jUrFHVaamW2gacElIs6IiHkRMQ8YD1wGHBkRfy3MdgOwh6TZkiYAbwGuaFaamqXLHVaamW1gfIu28wlgB+BMSQBrcjHXGknvBa4CxgHnR8TtLUpTaSrusNLMbAMtCS4R8Q7gHf1MuxK4shXpaJauai/P2HnrdifDzGzU8BP6I7R6Teqwco7rW8zM1nNwGaFFy1aydl2w21S3FDMzq3FwGaH1zZCdczEzW8/BZYQq7rDSzOwpHFxGqKvaw9StN2Nrd1hpZraeg8sIVao9zrWYmdVxcBmBiKCr2us+xczM6gz6nIukqcCLgGnAKuA2YEFErGty2ka9Zb2r+ceqJ9wbsplZnX6DS+5k8iRge+AmYAkwETgU2E3S94AvRsSjLUjnqFR7+6Tf42JmtqGBci4HAsdFxKL6CZLGAwcBrwK+36S0jXrrO6x0zsXMbAP9BpeI+PAA09YAP2hGgjYmle7UYeW0Se6w0sysaNAKfUknSNpGydck/VnSq1uRuNGua0kPcya7w0ozs3qNtBZ7W65XeTWwHXAk8JmmpmojkV5t7PoWM7N6jQSX2m35gcA3cpf4HX+rXuuw0vUtZmZP1UhwuVHS1aTgcpWkrYGOb4a8aFkva9eFcy5mZn1o5H0ubwfmAZWIWClpB+DYpqZqI7BwSa0ZsnMuZmb1Gsm5BDAXeH8e3pL0vEtHq3SnZsizJzvnYmZWr5HgcibwQuCwPLwCOKNpKdpIVKq97rDSzKwfjRSL7RsRe0u6CSAilkua0OR0jXpd1R4XiZmZ9aORnMsTksaRiseQNIUOr9CPCCpVN0M2M+tPI8HlK8DlwI6STgN+C3y6qaka5ZbmDiudczEz69ugxWIRcbGkG4FXkJ5vOTQi7mx6ykYxv33SzGxgjb7PZTKwMiJOB7olzW5imka9LndYaWY2oEb6FjsZOBH4aB61KfDNZiZqtKtUe9hs/CZMd4eVZmZ9aiTn8nrgYKAXICIeBLZuZqJGu0q1l9mTt2QTd1hpZtanRoLL6ogInmwt1vEVDW6GbGY2sEaCy3ckfRWYJOk44OfAuc1N1uj1+Jq1LF6+ypX5ZmYDGLC1mCQB3waeDjwKPA34RERc04K0jUqLlq5k7bpwzsXMbAADBpeICElXRsSzgY4NKEVdboZsZjaoRorF/izpeU1PyUai1gx5jnMuZmb9aqhvMeAISfeRWoyJlKl5TlNTNkpVqr3suM1mbLVZI1+dmVlnauQK+Zqmp2IjUunuYc5k51rMzAbSSLHYij7+HmxmokariKBrSQ+7TXV9i5nZQBqqcwGqwF+Bv+XP90r6s6R9mpm40WZp72oefWyNcy5mZoNoJLhcAxwYEZMjYgfgAODHwLtJLxJrmKTnSVoj6Y2FcUdL+lv+O3oo62u1riW5T7GpDi5mZgNpJLi8ICKuqg1ExNXACyPiemCzRjeU3wnzWeDqwrjtgZNJjQaeD5wsabtG19lqle7cDNmvNjYzG1AjweXvkk6UtGv++wjwcA4WQ3lp2PuA7wNLCuNeA1wTEcsiYjkpl7T/ENbZUl1L3GGlmVkjGgkuhwMzgB+QXho2M48bB7ypkY1Imk7qAPOsuknTgcWF4fvzuPrlj5e0QNKCarXayCabotLtDivNzBrRyMvCuoH3SdoyInrrJi9scDtfBk6MiHWpR5mhiYhzgHMA5s+fH0NeQUkq1R6eOW3bdm3ezGyj0cj7XP5J0h3AnXl4L0mDVuRLeo+kmyXdDMwHLpV0L/BG4ExJhwIPkHJCNTPyuFHn8TVrWbRsJbu52xczs0E18hDll0h1I1cARMRfJL10sIUi4gzgjPrxki4AfhwRP8gV+p8uVOK/midfSjaqLFq6knXhbl/MzBrRUB8mEbG4rjhrbRkbj4hlkk4FbsijPhURy8pYd9n8amMzs8Y1ElwWS/onICRtCpxALiIbjog4pm74fOD84a6vVWq9Ic92sZiZ2aAaaS32TuA9pFZcDwDz8nBH6ar2uMNKM7MGNdpa7IgWpGVUq1R7XSRmZtagfoOLpP8D+m32GxHvb0qKRqGIoFLt4eB509qdFDOzjcJAxWILgBuBicDepE4r/0YqFpvQ9JSNIt09qcNK51zMzBrTb84lIi4EkPQu4MURsSYPnw1c15rkjQ4Vv33SzGxIGqnQ3w7YpjC8VR7XMWotxfwApZlZYxpp+vQZ4CZJvyK94vilwCnNTNRoU6n2MHHTTZi2rTusNDNrRCOtxb4u6aekbvEh9RH2UHOTNbp0VXuYtYM7rDQza9RArcVmRcS9ADmY/LBuuoDpEXF/U1M4ClS6e3nWdHdYaWbWqIFyLp+XtAkpqNxIer3xRGB34OXAK0gv+hrTweXxNWtZvGwlh+zlZshmZo0aqLXYv0qaS3qA8m3AzsBKUtcvVwKnRcRjLUllG92XO6z0q43NzBo3YJ1LRNwB/GeL0jIqrW+GPNnBxcysUY00Re5otWbIc9wM2cysYQ4ug+iq9rDTNhPZ0h1Wmpk1zMFlEF3VXudazMyGqJHXHEvSWyV9Ig/vIun5zU9a+9U6rHSfYmZmQ9NIzuVM4IXAYXl4BX28vngs6u5ZzYrH1jjnYmY2RI1UJOwbEXtLugkgIpZL6ohekf1qYzOz4Wkk5/KEpHHkd7tImgKsa2qqRomKW4qZmQ1LI8HlK8DlwFRJpwG/BT7d1FSNEl3usNLMbFga6bjyYkk3krp7EXBoRNzZ9JSNApVqD7Mnb+UOK83MhqiR1mK7AfdExBnAbcCrJE1qdsJGAzdDNjMbnkaKxb4PrJW0O/BVYCbwraamahR4fM1a7l++0pX5ZmbD0EhwWZdfcfwG4PSI+DCpE8sxbX2Hlc65mJkNWaOtxQ4DjgJ+nMdt2rwkjQ5dS9wM2cxsuBoJLseSHqI8LSLukTQb+EZzk9V+le7UDHn2ZOdczMyGqpHWYncA7y8M3wN8tpmJGg26lvSw87busNLMbDgGvXJK2gP4H2Au6U2UAETEnCamq+26ut1SzMxsuBopFvs6cBawhvR644uAbzYzUe0WEVSW9PgFYWZmw9RIcNk8In4BKCLui4hTgNc2N1ntVe15nBWPr3FLMTOzYWqkQuFxSZsAf5P0XuABYEzf0nctqfUpNqZ308ysaRrJuZwAbEGq1N8HOBI4upmJardKd26GPNXBxcxsOBppLXZD/thDapY85lWqvUzcdBN23mbi4DObmdlTNNK32J6SzpV0taRf1v6GuiFJ+0m6WdLtkq4tjN9f0t2SFko6aajrbYauaqrMd4eVZmbD00idy3eBs4FzgbXD2Uju6PJMYP+IWCRpah4/jvRWy1cB9wM3SLoiP1vTNpVqL8+ZsW07k2BmtlFrJLisiYizRridw4HLImIRQEQsyeOfDyyMiAqApEuBQ4C2BZfHnljL4uUrOfS509uVBDOzjV6/xWKStpe0PfAjSe+WtHNtXB4/FHsC20n6taQbJR2Vx08HFhfmuz+Pa5v7lq4k3GGlmdmIDJRzuZH0auNaxcOHC9MCGMoT+uNJLc1eAWwO/EHS9Y0uLOl44HiAXXbZZQibHbquqjusNDMbqX6DS0TMHsmKJb0HOC4Pfge4KiJ6gV5JvwH2IuVUZhYWm0F6jqY+LecA5wDMnz8/RpKuwVRycHGHlWZmw9dIa7GJkj4g6TJJ35f075IGbaMbEWdExLyImAdcDrxY0nhJWwD7AncCNwB7SJotaQLwFuCKEe3RCFWqve6w0sxshBq5gl4ErAD+Lw8fTupy/18b3UhE3CnpZ8AtwDrgvIi4DSA/9X8VMA44PyJubzz55euq9rjDSjOzEWokuDwrIuYWhn8lacituSLi88Dn+xh/JXDlUNfXDBFBpdrL6/d2SzEzs5FopPuXP0t6QW1A0r7AguYlqX2qK1KHlXNc32JmNiKN5Fz2AX4vaVEe3gW4W9KtQETEc5qWuhbrqqYOK92nmJnZyDQSXPZveipGiVozZPeGbGY2Mo10XHlfKxIyGlSqvWy+6Th3WGlmNkKN1Ll0jEp3D7Mnb+kOK83MRsjBpcDNkM3MyuHgkj32xFruX77K3b6YmZXAwSW7d2kvETjnYmZWAgeXrFJrhuyci5nZiDm4ZF1Las2QnXMxMxspB5es0t3LtG0nssUEd1hpZjZSDi5ZpdrjhyfNzEri4ELqsLKr2usiMTOzkji4kDqs7Hl8jSvzzcxK4uACLKy6Mt/MrEwOLrgZsplZ2RxcSN2+bL7pOHZyh5VmZqVwcCHlXOZMcYeVZmZlcXAh9YbsZshmZuXp+OBS67DSrzY2MytPxweXWoeVfrWxmVl5Oj64dC1JLcWcczEzK0/HB5eKn3ExMytdxweXrmqPO6w0MytZxweXSnev61vMzErW0cElItIzLq5vMTMrVUcHlyW5w0o/42JmVq6ODi5duTLffYqZmZWrw4NLbobslmJmZqXq6OBSqfawxQR3WGlmVraODi5d1V5mT3aHlWZmZevo4FKp9ri+xcysCTo2uDz2xFoeeGSV61vMzJqgY4PLPd2pw0o3QzYzK19LgoukbSX9SNJfJN0u6djCtKMl/S3/Hd2K9EDx1cbOuZiZla1VHWq9B7gjIl4naQpwt6SLga2Ak4H5QAA3SroiIpY3O0G1Z1xm++l8M7PStapYLICtJYkUUJYBa4DXANdExLIcUK4B9m9FgirVHqZP2twdVpqZNUGrgsvpwDOAB4FbgRMiYh0wHVhcmO/+PG4Dko6XtEDSgmq1WkqCuqq9rsw3M2uSVgWX1wA3A9OAecDpkrZpdOGIOCci5kfE/ClTpow4ManDyh53WGlm1iRNCy6S3iPpZkk3k+pcLotkIXAP8HTgAWBmYbEZeVxTLVnxOL2r17qrfTOzJmlacImIMyJiXkTMA+4CXgEgaUfgaUAFuAp4taTtJG0HvDqPa6quJfntk5MdXMzMmqFVtdmnAhdIuhUQcGJEdANIOhW4Ic/3qYhY1uzEdHXnZshTXSxmZtYMLQkuEfEgKVfS17TzgfNbkY6ariXusNLMrJk68gn9SndqKZZaRpuZWdk6Mrh0LelxfYuZWRN1XHB57Im1PPgPd1hpZtZMHRdcah1Wuqt9M7Pm6bjgUutTzDkXM7Pm6bjgUusN2XUuZmbN03HBpSt3WLn5hHHtToqZ2ZjVccGl4g4rzcyarqOCS63DSlfmm5k1V0cFl4cfTR1WOudiZtZcHRVcai3FnHMxM2uujgouFTdDNjNriY4KLl3VXndYaWbWAh0WXHrcYaWZWQt0VHCpVHtd32Jm1gIdE1xWrV7LA4+s8pP5ZmYt0DHBpXf1Gg7eaxp77zqp3UkxMxvzWvWa47abvNVmfOWw57Y7GWZmHaFjci5mZtY6Di5mZlY6BxczMyudg4uZmZXOwcXMzErn4GJmZqVzcDEzs9I5uJiZWekUEe1Ow5BIqgL3jWAVk4HukpKzsei0fe60/QXvc6cYyT7vGhFTykzMQDa64DJSkhZExPx2p6OVOm2fO21/wfvcKTamfXaxmJmZlc7BxczMSteJweWcdiegDTptnzttf8H73Ck2mn3uuDoXMzNrvk7MuZiZWZM5uJiZWek6JrhI2l/S3ZIWSjqp3elpBkkzJf1K0h2Sbpd0Qh6/vaRrJP0t/9+u3Wktm6Rxkm6S9OM8PFvSH/Px/rakCe1OY5kkTZL0PUl3SbpT0gvH+nGW9B/5d32bpEskTRxrx1nS+ZKWSLqtMK7P46rkK3nfb5G0d/tS/lQdEVwkjQPOAA4A5gKHSZrb3lQ1xRrggxExF3gB8J68nycBv4iIPYBf5OGx5gTgzsLwZ4EvRcTuwHLg7W1JVfP8L/CziHg6sBdp38fscZY0HXg/MD8ingWMA97C2DvOFwD7143r77geAOyR/44HzmpRGhvSEcEFeD6wMCIqEbEauBQ4pM1pKl1E/D0i/pw/ryBdcKaT9vXCPNuFwKFtSWCTSJoBvBY4Lw8L+Gfge3mWMbXPkrYFXgp8DSAiVkfEI4zx40x6LfvmksYDWwB/Z4wd54j4DbCsbnR/x/UQ4KJIrgcmSdq5JQltQKcEl+nA4sLw/XncmCVpFvBc4I/AjhHx9zzpIWDHdqWrSb4MfARYl4d3AB6JiDV5eKwd79lAFfh6Lgo8T9KWjOHjHBEPAF8AFpGCyj+AGxnbx7mmv+M6qq9rnRJcOoqkrYDvA/8eEY8Wp0Vqez5m2p9LOghYEhE3tjstLTQe2Bs4KyKeC/RSVwQ2Bo/zdqQ79dnANGBLnlp8NOZtTMe1U4LLA8DMwvCMPG7MkbQpKbBcHBGX5dEP17LL+f+SdqWvCV4EHCzpXlJx5z+T6iMm5eITGHvH+37g/oj4Yx7+HinYjOXj/ErgnoioRsQTwGWkYz+Wj3NNf8d1VF/XOiW43ADskVuWTCBVBF7R5jSVLtc1fA24MyL+X2HSFcDR+fPRwA9bnbZmiYiPRsSMiJhFOq6/jIgjgF8Bb8yzjbV9fghYLOlpedQrgDsYw8eZVBz2Aklb5N95bZ/H7HEu6O+4XgEclVuNvQD4R6H4rO065gl9SQeSyubHAedHxGntTVH5JL0YuA64lSfrHz5Gqnf5DrAL6XUFb4qI+krDjZ6k/YAPRcRBkuaQcjLbAzcBb42Ix9uYvFJJmkdqwDABqADHkm4Wx+xxlvRJ4M2kVpE3Ae8g1TGMmeMs6RJgP1LX+g8DJwM/oI/jmoPs6aTiwZXAsRGxoA3J7lPHBBczM2udTikWMzOzFnJwMTOz0jm4mJlZ6RxczMysdA4uZmZWOgcXs0zSMZKmjWD5WZIOH2D6zrVemxtc39G5aWpx3GRJVUmbSbpU0h7DTa9ZMzm4mD3pGFLXIsM1C+g3uAAfAM4dwvouB14laYvCuDcCP8rPcpxF6lPNbNRxcLExS9IH8rs/bpP073ncrLp3ZXxI0imS3gjMBy6WdLOkzSXdK+lzkm6V9CdJu+dlLsjz19bRkz9+BnhJXv4/+kjSvwA/y8scI+kH+f0c90p6b07vTZKul7R97hfuWuB1hXW8BajlZq4DXlno/sRs1HBwsTFJ0j6kp9b3Jb3b5jhJz+1v/oj4HrAAOCIi5kXEqjzpHxHxbNKT0F8eZLMnAdfl5b9Ul57ZwPK6p8efBbwBeB5wGrAyd0T5B+CoPM8lpIBCLrLbE/hlTvM6YCHpfS5mo4qDi41VLwYuj4jeiOghdXT4kmGs55LC/xeOID07k7rJL/pVRKyIiCqpC/kf5fG3korYAH4CvEjSNsCbgO9HxNrCOpYwsqI8s6ZwcLFOs4YNf/cTB5k/+vi8fh2SNiH17zWYVX1sq5iLWVcYXkfqVp+cg/oZ8Ho2LBKrmZjXbTaqOLjYWHUdcGjuRXdL0sX5OlJngFMl7SBpM+CgwjIrgK3r1vPmwv8/5M/3AvvkzwcDmw6wfM1feTI3MlSXkBoD7FhIQ82ewG1PWcKszRxcbEzKr3u+APgTqVfo8yLipvwukE/l8dcAdxUWuwA4u1ahn8dtJ+kW4ASgVkl/LvAySX8hFZX15vG3AGsl/aW+Qj8ieoGuWqOAIbqGVPT17Sj0NCtpR2BV7oLfbFRxr8hm/cgvIJsfEd0lre/1wD4R8fGS1vcfwKMR8bUy1mdWJjdhNGuRiLhc0g4lrvIR4Bslrs+sNM65mJlZ6VznYmZmpXNwMTOz0jm4mJlZ6RxczMysdA4uZmZWuv8PxggsNhIjwX0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "do1d(output, 0, .1, 11, 1, In_R, In_phase, do_plot = True,measurement_name='output_varied')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2022-01-14T09:36:54.222509Z",
     "start_time": "2022-01-14T09:36:54.218509Z"
    }
   },
   "outputs": [],
   "source": [
    "SR.osc_amplitude(0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "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.8.10"
  },
  "nbsphinx": {
    "execute": "never"
   },
  "toc": {
   "base_numbering": 1,
   "nav_menu": {},
   "number_sections": true,
   "sideBar": true,
   "skip_h1_title": false,
   "title_cell": "Table of Contents",
   "title_sidebar": "Contents",
   "toc_cell": false,
   "toc_position": {},
   "toc_section_display": true,
   "toc_window_display": false
  },
  "varInspector": {
   "cols": {
    "lenName": 16,
    "lenType": 16,
    "lenVar": 40
   },
   "kernels_config": {
    "python": {
     "delete_cmd_postfix": "",
     "delete_cmd_prefix": "del ",
     "library": "var_list.py",
     "varRefreshCmd": "print(var_dic_list())"
    },
    "r": {
     "delete_cmd_postfix": ") ",
     "delete_cmd_prefix": "rm(",
     "library": "var_list.r",
     "varRefreshCmd": "cat(var_dic_list()) "
    }
   },
   "types_to_exclude": [
    "module",
    "function",
    "builtin_function_or_method",
    "instance",
    "_Feature"
   ],
   "window_display": true
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}