{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEWCAYAAABxMXBSAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/av/WaAAAACXBIWXMAAAsTAAALEwEAmpwYAAB9BUlEQVR4nO2dd1xV9fvA3x+GIogLxa2AosgQEHGvMlfDVWZmw8r0m+1s2K8y03ZZWllmORpmOdJMLcuVmhMciLgVBQcIKjJk3uf3x73c2Bsu4Of9ep3Xveecz3jO4XKe8/k8n+d5lIig0Wg0Gk1RsbK0ABqNRqOpWmjFodFoNJpioRWHRqPRaIqFVhwajUajKRZacWg0Go2mWGjFodFoNJpioRWHRqPRaIqFVhwaTT4opfYopdoppdyUUvuyHP9RKXVRKXVdKXVcKTU+R73xSqmTSqkEpdSfSqlmWc7dopTarJSKU0qF59Gnn1Jqm+l8pFLqjSznxprazNySlFKilAoop1ug0eSJVhwaTR4opWyB1sAJIADYl+X0e4CLiNQBhgJvZz68lVL9gHeBYUAD4AywJEvdRGAB8FI+Xf8EbDXV7QtMUkoNBRCRxSJSO3MDJgGnc8im0ZQ7WnFoNHnjDYSJMbRCZ7I8nEXksIikZO6atjam/TuBZaYyqcAMoI9Sqo2p7h4R+QHjAz8vXIDFIpIhIqeA7YBXPmUfBr4XHf5BU8FoxaHRZEEp9YhS6hrwL9Dd9H0y8IFS6ppSytVU7kulVBJwFLgIrMvaTB7fvYsowizgIaWUrVKqPdAd2JCHnK2BPsD3RWxXoykztOLQaLIgIgtFpB4QDHQDOgKhQB0RqSciZ0zlJgGOQG/gVyBzBPIncK9SqqNSqhYwFeOIxL6IIqwB7gFuYFRK80Vkbx7lHgK2Zcqj0VQkWnFoNCaUUg1Mo4o4oAewBTgGtAeuKqWey1reNJ20HWgBPGE6tgF4E1gBhJu2eCCyKP1jVDzTATugJTBIKTUpj+IPAd8V9xo1mrJAKw6NxoSIXDGNNiYC35q+/wncZRptzMqnqg3/2TgQkTki4i4ijTEqEBuMo5bCcAMyROR7EUkXkUjgZ+D2rIWUUj2BZsDy4lyfRlNWaMWh0eQm6yoqf4zTVgAopZyVUvcppWorpayVUoOAMcBG03k7pZS3MtIKmAfMFpGrpvNWSik7wNa4q+yUUjVMzR83HbvfVK4JMBoIySHfw8AKEYkvl6vXaApBKw6NJjcBwD6llBPGEcDVLOcE47RUJHAV+Bh4TkRWm87bYVxSmwDsAXYCb2Sp3wej/WId0Mr0/S8AEbkOjASeN7V9AONI5e3Myialcy96mkpjQZReyafRaDSa4qBHHBqNRqMpFhZVHEqplqbwC2FKqcNKqWdNx6cppc4rpQ6YttsLa0uj0Wg0FYNFp6qUUk2BpiKyTynliNEIORzjHG6CiHxsMeE0Go1Gkyc2luxcRC5i9LpFROKVUkeA5paUSaPRaDQFU2mM40opF4zB3byBF4BxwHUgCJicY2VLZp0JwAQABweHAA8Pj4oSV6PRaKoFwcHBMSLSqDh1KoXiUErVBv4B3hGRX5VSjYEYjEsfZ2Ccznq0oDY6d+4sQUFB5S+sRqPRVCOUUsEi0rk4dSy+qsoUvnoFxoigvwKISJQpnIMB+AboYkkZNRqNRvMfll5VpYD5wBER+STL8aZZio2gaOEaNBqNRlMBWNQ4DvQEHgQOKaUOmI79HzBGKeWHcaoqHGPsII1Go9FUAiy9qmo72XMXZLIuj2OaakhaWhqRkZEkJydbWhRNKbCzs6NFixbY2tpaWhRNBWDpEQdKqZYYk9E0xjjCmCcis00hpn/BmBEtHLg3r5VVmqpNZGQkjo6OuLi4YJy51FQ1RITY2FgiIyNxdXW1tDiaCsDixnEgHeNyW0+MiXOeVEp5AlOAjSLijjHy6BQLyqgpJ5KTk3FyctJKowqjlMLJyUmPGm8iLK44ROSiiOwzfY8HMp0Ah/FfBNDvMHqUa6ohWmlUffTf8ObC4oojKyYnQH9gN9DY5FkOcAnjVFbO8hOUUkFKqaDLly9XnKAajUZTBJRSKKW4du0alcFnrqwoseIwpdksbKtXjPZqY/TneM6Ul8CMGO94rrsuIvNEpLOIdG7UqFiOjxqNRlNh1K9fHysrK+bPn29pUcqE0hjHL5i2gsao1hiT1RRIXk6AQJRSqqmIXDT5dUSXQlaNRqOxOBMnTqRDhw706NHD0qKUitJMVR0RETcRcc1vA2ILayQ/J0BgNcYUmZg+fyuFrBpNgdSuXdvSIuRLpmzXrl3jyy+/LPf+Xn75ZT7++L/A1CJCvXr1SExMLPe+qxMiglKKJ598kkuXLnHu3Dnc3Ny46667OHbsmKXFKxWlURzdy6hMphPgrTnyb7wPDFBKnQBuM+1rNDctFaU4Dh06RMeOHc37Z86coVGjRjg4OJR739WJ5ORkRITmzZvTuHFjWrZsyZ9//omNjQ3Dhg0jLS3N0iKWmNIojplKqV4FFRCRQtfnich2EVEi0lFE/EzbOhGJFZH+IuIuIreJyJVSyKrRFJtPPvkEb29vvL29mTVrVoHHw8PD8fDwYOzYsXTo0IF77rmHpKSkXG1OmTKFOXPmmPenTZtmfrvPr7+sdU+dOoWfnx8vvfQSAMOHDycgIAAvLy/mzZtnLjtjxgzat29Pr169GDNmjLmPH3/8kS5duuDn58fEiRPJyMjI1U9OxRESEpJtX1M0Mv/+WRWum5sb8+fP59ixY7z99tv5Va30lMbGcRz4yGR/WAosEZH9xW1EKbUAuBOIFhFv07FpwONA5lKp/xMR7U1ezXnuuec4cOBAmbbp5+eX50O4MIKDg1m4cCG7d+9GROjatSt9+/bFYDDkebx+/focO3aM+fPn07NnTx599FG+/PJLXnzxxWztjh49mueee44nn3wSgKVLl7J+/fp8+/P39zfXff/99wkNDc12jxYsWECDBg24ceMGgYGB3H333Zw+fZoVK1Zw8OBB0tLS6NSpEwEBARw5coRffvmFf//9F1tbWyZNmsTixYt56KGHzO1dvXqVtLQ0mjRpYj6mFUfJyFQc9vb22Y7feeed3HfffUyfPh1PT09Gjx5tCfFKRYlHHCIyW0S6A30x2jIWKKWOKqXeVEq1K0ZTi4DBeRz/NOsIpKRyajQlYfv27YwYMQIHBwdq167NyJEj2bZtW77HAVq2bEnPnj0BeOCBB9i+fXuudv39/YmOjubChQscPHiQ+vXr07JlywLbLYjPPvsMX19funXrRkREBCdOnODff/9l2LBh2NnZ4ejoyF133QXAxo0bCQ4OJjAwED8/PzZu3Mjp06eztXfo0CF8fHyyHcupOMaNG1ese3mzkmkTymuKb9GiRXTr1o2nn36auLi4ihat1JQ65IiInAU+AD5QSvkDC4CpGFdUFaX+VpP/huYmpyQjg8pETie4/JziRo0axfLly7l06VKp3ja3bNnChg0b2LlzJ/b29vTr169A720R4eGHH+a9997Lt8yxY8do06aNed9gMPDvv//y9ddfA8a3aG3rKBr5jTgAatasyRdffEFgYCDt2rVj+fLl9O7du6JFLDGldgBUStkope5SSi0G/gCOASNLLRk8pZQKUUotUErVz6dv7QCoKRd69+7NqlWrSEpKIjExkZUrV9K7d+98jwOcO3eOnTt3AvDTTz/Rq1feJsDRo0fz888/s3z5ckaNGlVgf1lxdHQkPj7evB8XF0f9+vWxt7fn6NGj7Nq1C4CePXvy+++/k5ycTEJCAmvWrAGgf//+LF++nOho48r2K1eucPbs2Wx9tG7dmr1795oV0LvvvkufPn1wcnICYN++fXTq1KnkN/YmoiDFARAQEMBjjz1GdHQ0ffr0ISEhoSLFKxWlcQAcYLJPRGK0R6wF2ojIfSJS2qWzXwFtAD+MOcln5lVIOwBqyoqkpCRatGhh3rZs2cK4cePo0qULXbt2Zfz48fj7+9OpU6c8jwO0b9+eOXPm0KFDB65evcoTTzyRZ19eXl7Ex8fTvHlzmjY1pp4pqN1MnJyc6NmzJ97e3rz00ksMHjyY9PR0OnTowJQpU+jWrRsAgYGBDB06lI4dOzJkyBB8fHyoW7cunp6evP322wwcOJCOHTsyYMAALl68mK2PgQMHcsstt+Dh4UH79u05ceIEc+fONZ/fu3cvgYGBZXbfqzOZU1X5KQ4wTjVmjjqnT59eIXKVCSJSog3YBIwH6pe0jSxtuQChxT2XdQsICBBN1SMsLMzSIpQJZ86cES8vL0uLYSY+Pl5ERBITEyUgIECCg4PLpN3HHntM0tPT8zxXXf6WZcXPP/8sgISGhhZa9tFHHxUbG5silS1rgCAp5jO7NMbxW0XkWymHUOc6A6BGUzomTJiAn58fnTp14u677y6z6aVvv/0Wa+simS9vejKnAVu0aFFo2Q8++ABHR0eeeuqpKhHTqsTGcaXUPhEp8NdYxDJLgH5AQ6VUJPAm0K+qZwC8cOECu3fv5q+//uLMmTM0a9aMF154AW9vb0uLpikHXFxcCA2tPO83P/30k6VFuOmJiIigbt261K1bt9CyDRs25L333uN///sf33//PQ8//HChdSxJaVZVdVBKhRRwXgGF3zG4gXEF1jH5z49jJf8lcaoFVPpA/xcvXuTzzz/PtWKldu3auLq6sn37dlauXElYWJh5Xluj0VRfzp49W6TRRiaPP/44ixcv5umnn6Zbt260b9++HKUrHaVRHB5FKJPbLTU3i4AvMGYBzCQzidP7Sqkppv1Xii1hOXPkyBFWrlzJypUrCQoKynbu008/pWvXrvj7+2NnZ8fx48fp2LEjTz75JCtWrND5CzSaakxCQgK///57sepYWVmxePFiOnXqxMiRI9m9e3fljaFWXKNIeWzkMIBjXNLb1PS9KcbRSKUwjp86dUrefvtt8fb2zgz1Ll26dJF3331XQkJCJCUlJd+677//vgCydOnSCpG1KhAWFiYGg8HSYmhKicFg0MbxLPz0008CSGBgYLHrbty4UaysrGTUqFHlIFluKIFx3OJKQ/JWHNeyfFdZ93PUmwAEAUGtWrUqy3uZjeTkZOnXr5/cdtttZmXRs2dP+eyzz+T8+fNFbictLU0CAgLEyclJDh06VG7yViVOnz4tly9f1sqjCmMwGOTy5cty+vRpS4tSach8TiQkJJSo/vTp0wWQtWvXlrFkuSmJ4lBSCSz4Js/xNfKfjeOaiNTLcv6qiOTpBJhJ586dJed0UWlJT0/nhx9+4NFHHzUfe/PNN3nkkUdo3bp1ido8ceIEffv25erVqwwbNgwXFxc6d+5M+/btad++PTVq1Cgr8asEaWlpREZGVli+6syIpLa2tsWum5GRQVpaGnZ2dkWuYzAYzFngiktqaipWVlbmrTJjZ2dHixYtCryvGRkZWFlZVeg07S233MKWLVsA4/0syd+9uJw7d878fCjp8zU1NRV/f3/CwsJISUkp1+eCUipYRDoXq1JxNU1+G7AQ+Bx4BPAGo1IqYl0XKtFUlcFgkBUrVoiHh4cA0qpVK3F3dy9wGqo4REZGysSJE6VJkybmNxNAatSoIV5eXnL//ffLkiVLJCkpqUz6u5k4e/aszJs3T55++mmJjo6Wfv36SYsWLeSbb74Rkf/eBPP6W964cUMyMjJyHd+1a5d899135rp//vmnXLp0yTwdAcjw4cPl4YcflkuXLsm3334rAQEB8sEHH2T7+zZu3FiUUvLBBx/Ihx9+KNOnT8/2Rpqeni4ff/yxvPvuu3LPPfdkq1vUbciQITJjxgyZO3eu3HvvvTJ9+nQ5f/68rFu3TtatWyenTp2S8PBwSUtLkytXrsjMmTPNdVu2bJnte3BwsGRkZMitt94q06ZNk2XLlklUVJQkJCTIrbfeKps3bzbLvnr1amnUqJH0799fRo8eLU888YTUr19fdu/eLS+88IL4+fmZ2545c6bEx8fLyZMn5e2335Y6deqYZe/QoYOsX79ekpOT5fDhw/Lnn3/KbbfdJkOHDhVAmjZtKr6+vgJI586dZezYsTJ79mzZsmWL+Pv7S6NGjeS5556Tvn375nl/AgIC5LHHHpOLFy9KRERE2f8ARaRRo0YCyKpVq0rVzqJFi8yzG/n5zpQFVORUFVArr2MY82s8B8wrRls5FcdHwBTT9ynAh4W1URaKw2AwyOrVq6Vr164CiIeHh/z666/lOo0SFxcn+/fvlx9++EGef/55GThwoPkfqUGDBnL//ffLZ599JteuXSs3GSqCkJAQ8fb2lubNm8tzzz0nERER8sMPP8jUqVPlzTfflDlz5sjPP/8sS5Yskc8//1wGDRokgYGB8sgjj8jhw4dly5YtUrNmzWwPAX9/f5k1a1aRHqi2trYSHx9v3m/Tpo2sW7dORo4cWaIHdFluHh4eMmzYMIvLUV23sLAwqVWrVr7ne/bsKffff7+MHTtWOnbsKJs2bZJTp05JSEiIbNu2TaKiosy/41OnTsn48ePlm2++Mb98xMbGytdffy3//vuv3HfffeZ283oJKQ4Gg0FGjBghgDzwwAPlpjyoYMVxFePIYBnwBjCMEniRA0swhhVJwxi+5DHACdgInAA2AA0Ka6e0iiMqKkruuOMOAePb1rfffitpaWmlarOkpKamyvr16+X++++Xpk2bmn+Ibdu2lZ49e8odd9whzzzzjGzatEl2794tW7dulVOnTuXZ1pUrVyQlJcWs/GJiYiQ6OloMBoMkJydLcnKybNmyReLi4iQlJSXbD3/58uUiYvzHSE5OlpSUFJk8ebL5YZf1n2/atGnStm1bue+++6Rdu3YCyKhRo2TixIkV+pBo2LCh3HnnnTJmzBjx8fGRTz75RObNm1cmbQ8fPlzWrFmT7djcuXPlf//7n9jY2MgXX3wh9evXlyZNmkifPn0EkLfffluioqIkNDRUnn32Wdm6dav0799fJk+eLLfccos88MAD0q5dO2nQoIE0adJEevbsKS1btpRatWrJ8ePHzX/HhIQEOXr0qERGRsq+ffskJiZGrly5IlFRUbJmzRo5e/asbNiwQR5//HFxdXWV++67Tzp16pTt72RnZyc2Njbm/cwR76ZNmyQ9PV22bt0q33//vVy+fFn27dsnDg4O5rLjx4/Pdt0BAQHSvHnzbMcGDBggTzzxhLzyyisCyIIFC+S5556TLl26yGeffSbHjx+Xf//9V0aNGmWu8+2330pkZKScP39eli5dKvfee6/4+/ubz0+dOlWOHDki7777rhw/flwmTZokf/31l6SkpMj58+clISFBfv/9d5k7d65MnTpVrl27JqtWrZJt27Zle+EzGAzy2muvSePGjSvkd7h48eIyex68/fbbAsi4ceNKrYzyokIVh7E/XE0K4w2TArlgUgR1S9NuSbbSKI5Zs2aJnZ2d1KxZU2bPni2pqaklbqs8+Oeff2TGjBly1113iZeXV6E/WltbW2nSpEk2pdO8eXNp1qxZhT7Ec25vvfWWLFiwQIYOHSpt27YVPz8/qV27tsyYMUP8/Pxk9uzZsnTpUpk8ebJ8++23smDBAgGkX79+8vrrr8vKlStl3759kpqaKgaDQU6cOCHvvvuu/PLLL5KSkpKvoj969Gg2OTZu3CiAdOjQQTp06CC7d++WdevWycaNG8VgMIjBYJCffvpJFixYIOfOnRMRkX379pmnDoODg/MdAWZ9WF28eLFcpxiKisFgKLEcKSkpEh0dXWCZtLS0YreflpZWKf7PYmJiZP78+XLp0iU5cOCAzJo1SyZOnCj9+/cXZ2dnqVevnvl/KvP3s379ennyySez/aZuueUW81Tf66+/XuZyvvnmmwLI/fffX2ZT5pmURHGUqXHclD/8NYzBDh8pg/bCgXiM/iDpUoABpyTG8ePHj/Pyyy/z22+/MXjwYGbOnImnp2epZC5vRITg4GAcHBzYv38/x44d4+jRo9ja2rJz505Onz5N48aNSUxMzObN3LBhQ9q3b0+DBg3o1KkTM2fORERo2bIlx44do1GjRnTs2BFnZ2c++ugjdu3axQ8//MCqVasA6NChA126dKF58+Y0adKEwMBA7Ozs+PHHH+nbty9ffvklgYGBXLlyhXHjxtG5c2fS0tIIDQ2lVatW5uiqluLYsWMsX76cp59+mjp16lhUFk3V5MCBA9SsWZMOHTpUeN8iwvvvv8///d//0bdvX5YsWVJmjsQlMY6Xy6oqpdQRESn13TUpjs4iElNY2eIqjuPHj9OjRw/S09N58cUXefXVV3UMHo1GU6lZsmQJjz76KA4ODnz77bcMHz681G2WRHGUJqz6C0qp25RSzjmO1wSKvl7RAnz11Vd07NgRpRR79uzh9ddf10pDo9FUesaMGcP+/ftp3bo1I0aMYMyYMbmyOFYEpVkc3hh4AQhWSl1USv2llPoK+BdYUSbSGecP/1JKBSulJuQ8WZJEThs2bGDSpEn06dOHXbt20a5dcbLcajQajWXx8PBg586dvPnmm/z22294eHjw9NNPExUVVWEylMlUlVKqHuADtAfOisjfpW7U2G5zETlvGtX8DTwtIlvzKlvYVJWImJ2o2rdvz759+wpMsKLRaDSVnQsXLjB9+nS+/fZbMjIyeOihh3j99ddxd3cvchsVOlWVFRG5JiLbxJifo0yUhqnd86bPaGAl0KWkbd19993m71u3btVKQ6PRVHmaNWvG3LlzOXLkCPfddx/ff/897dq1w8fHh2nTphEcHIzBYCjzfitFyJG8UEo5AFYiEm/6/jcwXUT+zKt8QSOOf/75h+HDh3Pt2jUuXrxIkyZNyk9wjUajsRAXL17kxx9/ZPXq1fz777+ICHXq1KFr165069aNbt260bVr12yrHCvNqqqyQCnlhnGUAcbw7z+JyDv5lc9PcSQlJdGyZUvq1KnDDz/8QK9evcpHYI1Go6lEREVFsW7dOvbs2cOuXbsICQkxjz5cXFzw9PSkfv36LF68uPoojuKSn+KYN28eEydOZOvWrfTu3dsCkmk0Go3lSUhIIDg4mB07dnDw4EEOHTrE+fPniYuLq16KQyk1GJiNMUPgtyLyfn5lGzVqJDt37qRt27bmYykpKeZIpplRSjUajUbzHyWZqipNBsByRSllDcwBBmCMYbVXKbVaRMLyKh8TE4O7uzuPPvooHTp0ICUlheDgYAAeeughrTQ0Go2mjKi0Iw6lVHdgmogMMu2/CiAi7+VTPt8LOXbsmPbX0Gg0mjyoViMOoDkQkWU/EuiatYDJKXACgJOTEy4uLnk2dP/995ePhBqNptphMBhIS0sjLS2N9PR082d+W0HLXa2trbGxsSl0s7a2Npe1QNKuTsWtUJkVR6GIyDxgHpRPBkCNRlP1MRgMXL16laioKKKiooiOji7w+40bN/Jsp3bt2jRs2DDX5uTklO/xisg4WFqUUnlfcAFUZsVxHmiZZb+F6ZhGo9EQHx/PhQsXuHjxovnz4sWLuZTB5cuXSU9Pz1Xf2tqaRo0a0bhxY5ydnWnXrh3Ozs40btzYfCzzs2HDhsVKGVzdqcyKYy/grpRyxagw7gP0nJNGU40REeLi4nIpg7y+JyYm5qpvZ2dHkyZNcHZ2plWrVnTu3DlPRdC4cWMaNGhQ6XO5V1YqreIQkXSl1FPAeozLcReIyGELi6XRaEpIcnIy58+fJzIy0rxduHAhl1JITk7OVdfBwYFmzZrRtGlTAgICzN8zPzO/16lTR6+gLD5FixCbhUqrOABEZB2wztJyaDSagklMTDQrhYiIiGzKIXOLicmdVqdOnTrmB3+3bt3yVAZNmzbF0dHRAld101BovqOcVGrFodFoLE9iYiLnzp3LVyFERkZy9erVXPWcnJxo0aIFLVq0oGvXrrRo0YKWLVuajzVv3pzatWtb4Io0pUUrDo3mJkZEuHr1KmfPns13y2uk4OzsTIsWLXB1daV3797ZFEKmUqhVq5YFrkhTEVQLxaGUGhwQEGBpMTSaSofBYCAqKorw8PB8FUNCQkK2OrVq1aJ169a0bt2agIAA8/dWrVrRokULmjVrRs2aNS10RZryQCnVAPgFcAHCgXtFJPcwMrN8ZfUcLyqm0CTHAwIC3LQfh+Zmw2AwcOHCBU6fPs2ZM2dyKYVz586RmpqarU79+vXNyiCvrWHDhtrAfBOhlDqM0ZZ8RUTeV0pNAeqLyCv51SnWiMOklQrDICLXitNuKekCnATcKrBPjabCuH79OmfOnOH06dPmLXP/zJkzuRRDkyZNaN26NZ06dWLEiBFmheDi4kKrVq2oU6eOha5EU0mpAQwD+pn2vwO2AGWjOIALpq2g1xFroFVBjRQ16q1S6m5gORAoIvkNJ3KGJtFoqhTp6elERETkUg6ZiiGnjaFu3bq0adMGb29vhg4dipubG25ubri6utKyZUvtqKYpLgmAm4hcNO1fAhoXVKG4iuOIiPgXVEAptb+Q80WKequUcgSeBXYXU0aNptIRHx/PyZMnOXnyZC7lcPbsWTIyMsxlbWxsaN26NW5ubtx99924urqalYObmxv169e34JVoqiHZgm2JiBQUNBaKrzi6l0GZLsBJETkNoJT6GeMwKWe49BnAB8BLhbSXMzSJRmMRMpXDiRMnOHHihPn7yZMniYqKyla2YcOGuLm50aVLF+677z7ziMHNzY0WLVpgY1Mt1q1oqg5RSqmmInJRKdUUiC6ocHF/nbkDvuRARHK7fWanKFFvOwEtRWStUipfxZElOq7P5cvFdn7UaIpNVuWQU0nkVA5NmzbF3d2dO+64A3d3d9q2bUvbtm1p06aNdmjTVDZWAw8D75s+fyuocHEVxx5KEIK3OCilrIBPgHGFlc2MjquUur1Ro0Zry1Muzc1DQkJCnqOGEydOFEk5uLu706ZNG+3cpqkq1MWoMJYqpR4DzgL3FlShuIqjLNboFRb11hHwBraYlgQ2AVYrpYbmZyAXkXWdOxcrD4nmJsdgMBAREcGxY8c4duwYR48eNX+eP589CHNlUQ5paWlERkbmGctJoykMOzs7WrRokVeo9zgRiQX6F7Wt4iqORkqpF/I7KSKfFKGNAqPeikgc0DBzXym1BXixgFVVGk2+JCQkcPz4cbNiyFQOx48fz5Z3oW7durRv355bb70VDw8P2rVrV+lGDpGRkTg6OuLi4qL9LDTFQkSIjY0lMjISV1fXUrdXXMVhDdSmFCOP/KLeKqWmA0EisrqkbWtuToo6erCyssLFxQUPDw9uvfVW2rdvj4eHB+3bt6dx48aV/mGcnJyslYamRCilcHJyoqxswcVVHJdEZHppO80r6q2ITM2nbL/S9qepHqSlpXHy5EnCwsLM25EjRwodPWQqiLZt21b5UBlaaWhKSln+doqrOKp2fBJNlSAlJYUTJ06YlcPhw4cJCwvj+PHj2TK5ubq6mkcPmQqiqoweNJqqTEmmqgpEKbVPRMp15ZWmepCcnMyxY8fMiiFzO3nypNkhzsrKijZt2uDp6cmwYcPw9PTE09OT9u3b4+DgYOEr0GhuToqrONyVUiEFnFcYl3ZpNGYSExM5evRoNuUQFhbG6dOnMRiMTqvW1ta4u7vj5eXFqFGjsikIHUJDo6lcFFdxdChCmYzCi2iqIykpKRw7doxDhw4RGhrKoUOHOHz4MOHh4eYytra2tG/fnk6dOvHAAw+YFYS7uzs1atSwnPCaQomNjaV/f+OKzUuXLmFtbU2jRo0A2LNnT6n/fj169GDHjh25jk+bNo3atWvz4osv5lu3du3aucLDW4rPPvuMr776ik6dOnHmzJk8rykr+cle0HV//fXXTJs2jcaNG5OQkIC3tzdLly6tsP+hYikOETlbXoJoqg4Gg4EzZ86YlUOmoshqg7C1tcXDw4Nu3brx2GOPmRVEmzZt8lpHrqkCODk5ceDAAaBoD/OciAgigpWVVZ7nC3vAVhW+/PJLNmzYQIsWLcqtj0OHDvHuu+/yyCOPYDAYaNeuHSEhIVSUP1vef0GNBuM/elRUFBs2bODTTz/lscceo0uXLjg6OtK2bVuGDx/OG2+8QVBQEG3btuWVV15hyZIlhIaGkpiYSEhICEuWLOH1119n5MiReHh4aKVRjfnkk0/w9vbG29ubWbNmARAeHk779u156KGH8Pb2JiIiIs9yQDZ/mXfeeYd27drRq1cvjh07Zj4+fPhwAgIC8PLyYt68eQXKEx4ejoeHB+PGjaNdu3aMHTuWDRs20LNnT9zd3dmzZ0+h7c6YMYP27dvTq1cvxowZw8cff1xgn//73/84ffo0Q4YM4dNPP812TT/++CNdunTBz8+PiRMnZgtsWdh15yQkJAR/f2O82ZMnTyIitGvXrkDZyhKLRFIrLKy6yclwPMbYWJeBR/Vop3yJj48nNDTUPIrI/Mwa0tvZ2RkfHx8ef/xxfHx88Pb2xsvLq9I4yN1MPPfcc+a3/7LCz88v24O8OAQHB7Nw4UJ2796NiNC1a1f69u1L/fr1OXHiBN999x3dunXLt1zmQzCzrZ9//pkDBw6Qnp5Op06dyMzwuWDBAho0aMCNGzcIDAzk7rvvxsnJKV+5Tp48ybJly1iwYAGBgYH89NNPbN++ndWrV/Puu++yatWqfNs9ffo0K1as4ODBg6SlpWWTIz/mzp3Ln3/+yebNm2nYsCFvvPEGAEeOHOGXX37h33//xdbWlkmTJrF48WIeeuihIl13Tg4fPsxDDz1EWloa58+fZ82aNdSpU4erV69WSPTkClccRQyrvh/oLCJJSqkngA+B0RUta3UkIyODU6dOcfDgQfN26NAhzp79Ty87ODjg7e3N8OHD8fb2NisJZ2dnC0quqcxs376dESNGmFe6jRw5km3btjF06FBat25Nt27dCiyXVXFs27aNESNGYG9vD8DQoUPN5z777DNWrlwJQEREBCdOnChQcbi6uuLj4wOAl5cX/fv3RymFj49PNttbXu3u2rWLYcOGYWdnh52dHXfddVeJ78/GjRsJDg4mMDAQgBs3buT6fyrourMSERGBs7MzISHGdUrff/89M2bM4O+//+b5559n0aJFJZazqJRacSilJovITNP39iKS//jKSKFh1UVkc5byu4AHSivnzcj169cJCQnJpiRCQ0NJSkoCjCuZ2rdvT/fu3ZkwYYJZSbRu3TrfeWhN5aCkIwNLUFbLprds2cKGDRvYuXMn9vb29OvXr9C4XVkdPq2srMz7VlZWZntcSdotLiLCww8/zHvvvVfqtg4dOoSnp6d539fXl5kzZ/Lnn39y9OhRPvroI156qbBsFKWjxE8HpVQ9pdRC4B6l1CSlVC9gShGq5hVWvXkB5R8D/shHhglKqSClVNDNHFbdYDBw6tQpfv31V6ZNm8aIESNwc3Ojbt269O7dm6eeeorly5djb2/P448/zoIFCwgODiYhIYHDhw+zZMkS/u///o+hQ4fi6uqqlYam2PTu3ZtVq1aRlJREYmIiK1eupHfv3iUq16dPH1atWsWNGzeIj4/n999/ByAuLo769etjb2/P0aNH2bVrV5nInl+7PXv25Pfffyc5OZmEhATWrFlT4j769+/P8uXLiY42prm4cuVKtlE+5H/dOQkJCaFDB+MCVxHhu+++47bbbqNhw4Y88MAD5a40oPg5x22AjsBxU17xR5RSg4AY0/Ffy1I4pdQDQGegb17nM8OqA3Tu3Pmm8GpPTEzk0KFDuaaa4uPjAeOblLu7O4GBgYwfPx5fX198fX1p3ry59qbWlBudOnVi3LhxdOnSBYDx48fj7++fbTqooHI5y4wePRpfX1+cnZ3N0zuDBw9m7ty5dOjQgfbt25unv0pLfu0GBgYydOhQOnbsSOPGjfHx8aFuXaOb2u233863335Ls2bNitSHp6cnb7/9NgMHDsRgMGBra8ucOXNo3bp1odedk0OHDvHPP/+wdu1arKys6Nq1Kx9//DE///wzvr6+pbwbRUOJFP15q5T6FfAx7T4CTAUaARuAN4qQxAmlVHdgmogMMu2/CiAi7+UodxvwOdBXRArMRgVGxREUVH0C6IoI586dy6YgQkJCzCsoAOrUqWNWDL6+vnTs2BFvb2/zHKmmenHkyBHzm6amYkhISKB27dokJSXRp08f5s2bR6dOlTMwxurVq1mxYgVTpkzJ93eS129IKRUsIsVax1tcG4cP0A5wx5jU6SkgBGPGqI9N+4VRYFh1AKWUP/A1MLgoSqOqk5GRwbFjx9i/fz/79u1j//79HDhwgKtXr5rLtG3blo4dO/LAAw+YFUXr1q31KEKjKUcmTJhAWFgYycnJPPzww5VWaYDRmJ6fQb2sKa7iiBfj6+5xpdQFEfkRQCn1IkZFUihFDKv+Ecbw7ctMD8ZzIlIxd6ScSU5OJjQ01Kwg9u/fT0hIiDm6q52dHT4+PowaNQp/f398fX3x8fHRS141Ggvw008/WVqESklxFUcTpdRDwAEgNfOgiIgp5WuRKCysuojcVky5KiVxcXEcOHDArCD2799PWFiY2fGnbt26ZmegTp064e/vj4eHBzY2FnGv0Wg0miJR3CfUNCAQ40qnFkqpwxiX0R7BaOu4abl06VI2BbF//35OnTplPt+kSRP8/f0ZOnQo/v7++Pv74+rqqqeaNBpNlaO4saqy+fgrpVpgtHt0BLaWoVyVlkyjdXBwsHm6ad++fVy6dMlcxs3NjU6dOvHoo4+alUSTJk0sKLWmuiAi+mVDUyKKsxCqMEo1JyIikRj9MPL0s6jqiAgREREEBwcTFBRk/oyNjQWMDnQdOnRg4MCBZgXh5+dnXrKn0ZQldnZ2xMbG4uTkpJWHplhk5hwvqxQFejLdhIgQGRmZS0lkxmqytrbG29ubYcOG0blzZwICAvDx8aFWrVoWllxzs9CiRQsiIyPLLG+05ubCzs6uzCL23pSKQ0Q4f/58LiWR+Q9pbW2Nl5cXd911l1lJdOzYUSsJjUWxtbXF1dXV0mJoNJU2Om5N4HsgAIgFRotIeEn7u3btGkFBQezZs4fdu3ezZ88es03C2toaT09P7rjjDrOS8PX11UpCo9Fo8qGyRsd9DLgqIm2VUvcBH1DE6LhpaWkcPHjQrCD27NnD0aNHzefbt2/PgAEDCAwMpHPnzvj6+mpPa41GoykGlhhxFBod17Q/zfR9OfCFUkpJAcsCzp8/z8CBA9mxYweJiYmAMX9E165deeCBB+jSpQuBgYHUq1ev7K9Io9FobiIsoTjyio7bNb8yJk/zOMAJYzBFM0qpCcAEMKa1vHLlCh4eHtkaunDhAitXrjTH2tdoNBpNNoodR6VKG8dzRsetTkEONRqNpiJQSt0obh1LJF44D7TMst/CdCzPMkqpGkBrYFFFCKfRaDSagrGE4jBHxzUphfuA1TnKrMYYcRfgW3IrFo1Go9FYiApXHCKSjjH8+nqMMa6WZkbHVUplRsCdDzgppc5gNJT/X0XLqdFoNDcJxfYotYiNowjRcZOBUUqp5cB7gGNe7WQ1jrdq1arc5NVoKgsiQmpqKra2tgCkpqYSExND/fr1OX/+PDt27CAiIgJ/f38uX75MTEwMP//8MwaDgR49elC3bl1SU1MxGAw0bNiQQ4cOcc899zB8+HBOnz5NYmIi7dq1w8rKin///Zc9e4zZEp5//nmSkpK4fPkyDg4O5oySIkJaWho1atSw5G3RlI6Ywotkp1gZACsSpdSdwO0iMkkp1Q94UUTuzK+8No5rKhvx8fE4OmZ/5xERjh8/ztWrVwkLC+PChQvs2LGDuLg4kpOTcXFx4fTp0xw4cABXV1cmTZrE+vXrSU1NZevW/+KI2tnZoZTCwcHBHBanImnZsiXp6elcvHgRMGaj9PPzo0uXLnh7e3P9+nUefPBBLl++TPPmzbWvVCWmJBkAK7PieA94EEgH7IA6wK8i8kBe5bXi0FQUFy9epGHDhua3/gsXLrBp0ybc3NxYvnw5MTEx/PHHH8TExNC1a1d2796Np6cnp06dIj093ZyPJZMOHTrQtGlTQkNDAYiJiaFhw4ZER/+X/NLd3Z0GDRoQHh5O7969sbGx4cKFCzg7O3P58mXi4+Px8PCgZcuWtGnThoYNGxIfH0+jRo1wc3MjNTUVKysrXFxcUEpRq1Ytrly5wsGDB3n00Uc5e/Ysn376KfPmzSM+Pp777ruPHj160KtXL9577z0SExO5cuUKIsLp06cJDw+nf//+2NjYkJaWxqZNm7JlrMyJjY0NHTp04M477+TMmTOICI0bN+bpp5+mTZs2OmijBbG44lBKNShCMYOIXCtmu/0oYMShlBocEBDwh1YcmrLAYDDwzjvv8M8//9C2bVsyMjIIDg6mc+fODBw4kFGjRpWo3f79++Po6Ej79u1JT0/Hzc2NUaNG0ahR7lQ2GRkZbNu2jRs3btC2bVvc3d1Le1nlisFg4MKFC1y8eJHr168D8Pnnn9OrVy+ioqJYtmwZZ8+eLbCNW265BXd3d/M0mre3N3fddRd16tSpiEu4aVFKBQMDgV8AFyAcuFdE8n0TKGvFkQxcAAp6fbAWkWIZJApSHKYQJscDAgLctOLQ5Ed8fDzbtm3D09OTZs2a8c033xATE0NYWBj169dnyZIl5gdeUfD19cXJyYnIyEjS09NxcHBg6NChPPPMM9SvX59Dhw7RsWNHtm/fTqdOnfTDD1i5ciVNmzalffv2ODo6sm/fPjZv3sx3333HkSNH8q03aNAgAIYNG4aLiwtdunTBycmposSu9pgS8q0DrojI+0qpKUB9EXkl3zplrDj2i4h/acsUs8/uwLSAgICBWnFo8uLq1at4eXmZ5+MLw8/Pj169etGvXz/s7e2JiorCw8ODV155hTZt2vDpp5/qnCvlxI0bN7CysuLGjRv8/fff/PXXX2zdupUTJ06YExFZW1vTpUsX2rdvT7du3XB1dcXR0ZFu3brpKa8SoJQ6AQjQT0QuKqWaAltEpH2+dcpYcdiZVkSVqkwx+7wHGBwQEPCYVhw3N9HR0SQkJJCamkpQUBB//vknN27c4MqVK2zZsoX//e9/zJ0711x+/fr1uLu7U6tWLRo3bkxGRobO915JuXz5MlOnTqVbt24cPnyYP//8k8OHD2MwGLKV69SpE507d+aOO+7grrvu0oqkCCil9gNuIlLPtK8wBpmtl28lESn3DegFzClmHTtgD3AQOAy8lU+5e4BvAwICRHPzYDAY5Mcff5Rbb71VGjZsKBjfmLJtderUMX+fPHmyiIhMnTpVBg4cKAaDwcJXoCktGRkZcujQIdm4caPMmjVLOnXqlOs38NZbb0lSUpKIiMTGxkpGRoaFpa58AEHANcn+XL0qBTyfy21VlVLKH7gfGAWcwbgi6vNi1FeAg4gkKKVsge3AsyKyK0c5PVVVjYmJieHIkSMEBwdz5swZPvvss2zn3d3dqVGjBjdu3OChhx6iTp06ZGRk4OnpaV71IyJ6JHGTkJiYSFJSEt999x0vvfRSvuUGDRrEkCFDsLW1xdvbm65du1KzZs0KlLTyYDKOO2LBqap2wBjTFoPRSv+iiLQuZbv2GBXHEyKyO8c5G4zGcVetOKo+6enpLFu2jLCwMI4cOcLatWtJTs49szlkyBBGjRrFQw89hLW1tQUk1VQFVq9ezbBhw4pU9q677iI1NZXbbruN7t27ExAQUGY5uiszJsWxGYiV/4zjDUTk5XzrlLHiMADbgMdE5KTp2GkRcSthe9ZAMNAW41TXKznOZ3qO123VqlXbwpb7aSofIsLBgwd57bXXuHbtGocPHyYuLg4ApRT9+vXjgQceYMiQITRp0sR8XKMpKSJCXFwcW7duNb+ULFq0iPDw8Fyru5o0acLYsWNp27Yt9evXp1u3brRuXar34EqHUuok0A1YCrQCzmJcjnsl3zplrDiGYwxa2BP4E/gZY2rYUiVKVkrVA1YCT4tIaF5ltANg1SAxMZEjR45w48YNUlNTSUtLM5+ztramVq1a1KhRA1tbW2rWrKmVRBXGzs6OFi1amB0lqwIXLlwgJSWFoKAg1qxZw4oVK0hPTyclJQUAKysrXn/9df7v//6PuLg4nJ2dLSxx6bG4A2AWQRwwBiccA9yKMX/4ShH5qxRtTgWSROTjvM5rxVH5SE1NZeXKlRw4cID4+Hh27NjBuHHj6NKlC7a2tlhZWWFlZUWzZs2oU6fOTTvHXB0REWJjY4mPj8fVtVTvjRYnIyODY8eO8ffff7N8+XK2b99uPtejRw/s7e356aef8nTkrApUGsWRrQOl6mM0kI8Wkf7FqNcISBORa0qpWsBfwAcisiav8lpxWB4RYevWrfzwww/Ex8ezbt06EhISAHBwcCAxMZE///yTdu3a0apVK22bqOaICEePHqVDhw6WFqXMSE1NZebMmVy7do2DBw9y6NAhLly4gKOjIz179qRPnz7cf//9VWo6qySKo7yW3wbkcezOYrbREdgPhAChwNSCyuvluJbBYDDI8uXLpXfv3tK4ceNsSyHHjBkja9eulbS0NHP5sLAwC0qrqWhuhr93UFCQjBs3Tpo0aWL+7e/atcvSYhUZIEiK+YwvrzWK3yilHhKTPUIpNQZ4DshztJAXIhIClJmHuabsiIyMZOvWrYwdO9Z8zMHBAT8/PyZPnsywYcNo1qwZtWvXtqCUGk3FEBAQwMKFCwGYOnUqM2bM4JFHHuHgwYNVyr5THMpLcdwDLFdK3Q/0Bh7CGESryCilWmK0jTTGqMXnicjsshZUUzRCQkL45JNPOHDgAAcPHsx27u233+bZZ5/VikJz0zN9+nSuX7/O7NmzWbBgARMnTrS0SOVCeToAtgNWAeeAESJSrIToJieUpiKyTynliHFZ7nARCcurvLZxlD3//PMP8+fPZ+PGjVy4cAE7OzuaNm1K3bp1mTJlCr1796ZZs2bFavPIkSPVas5bUzA34987PT0dW1tb3N3dOXr0KFZWlsjQXXRKYuMo0ytSSh1SSoUopUKA5UADwBXYbTpWZETkoojsM32Px5hmtnlZyqvJjohw4sQJNmzYwB133EG/fv349ddf6dq1K7Nnz+bcuXOcPn2a/fv3M3r06GIrjcqEtbU1fn5+eHt7M2rUKJKSkgBu+lHTtWvX+PLLL/M9//XXX1fbt+iywsbGhjlz5nDixAk++ugjS4tTLpS1KrwTuCvL1hXjFFXmfolQSrlgtHfk9BqfoJQKUkoFXb5c7LS5GhOhoaE8/fTTNGzYkHbt2jFgwAC2bt3K66+/TnR0NL/++ivPPPNMlV1umBe1atXiwIEDhIaGUqNGjWzBD6sqIpIt6F/O/aJQmOLIDBevKZhM+99XX31lYUnKhzJVHCJytqCtJG0qpWoDK4DnRCRbwgQRmScinUWkc3V6qJU3IsKuXbt48803ad26NT4+PnzxxRcEBgbywQcf8NFHH3H27FlmzJhxU6T87N27NydPnsx2bPjw4QQEBODl5cW8efMAo/PiHXfcga+vL97e3vzyyy8AhIeH4+Hhwbhx42jXrh1jx45lw4YN9OzZE3d3d3Pe7rzazIvvv/+ejh074uvry4MPPkh4eDje3t7m8x9//DHTpk0z992+fXseeughvL292bZtW7b9iIgIfvzxR7p06YKfnx8TJ04kIyOD8PBwOnTowOOPP46XlxcDBw7kxo0bTJkyhVOnTuHn55dnrKeQkBB8fHxKdb9vBurWrcvYsWNzZXusLpSpcVwptU9EOpW2TJaythiVxmIR+bUsZLyZOXLkCH/++SeLFy8mODgYgM6dO3PvvfcyadKkCnfUeu655zhw4ECZtunn58esWbOKXD49PZ0//viDwYMHZzu+YMECGjRowI0bNwgMDOTuu+9my5YtNGvWjLVr1wKYQ6MAnDx5kmXLlrFgwQICAwP56aef2L59O6tXr+bdd99l1apVebaZMyHR4cOHefvtt9mxYwcNGzbkypUrhSaYOnHiBN999x3dunUjPDw82/6RI0f45Zdf+Pfff7G1tWXSpEksXryYPn36cOLECZYsWcI333zDvffey4oVK3j//fcJDQ3N9+8SGhqqRxxFxMvLi8WLF5OQkFDtpkDLelVVh0JsGQooUgYcU3Tc+cAREfmkLIS7GYmJieGjjz7ir7/+Mj8M2rZty5w5c7j33ntp2LChZQW0EDdu3MDPzw8wjjgee+yxbOc/++wzVq5cCUBERAQnTpzAx8eHyZMn88orr3DnnXfSu3dvc3lXV1fzm7iXlxf9+/dHKYWPjw/h4eH5tplTcWzatIlRo0aZ/y4NGjQoVHG0bt2abt265bm/ceNGgoODCQwMNF+3s7Mzffr0wdXV1XwPAgICCA8Pp1evXvn2ExERgYODA/Xq1StQHo2RFi1aAMYc9ZU99W9xKWvF4VGEMkUdu/UEHgQOKaUOmI79n4isK4lgNxOpqaksXbqUJUuWsG6d8Xa5u7vz8ccfM3jwYDw8PCqF13ZxRgZlTaaNIy+2bNnChg0b2LlzJ/b29vTr14/k5GTatWvHvn37WLduHa+//jr9+/dn6tSpANnCpVhZWZn3raysSE9Pz7fNomBjY5PNVpGznoODQ777IsLDDz/Me++9l61MeHh4Npmtra25caPghY9Z7RvR0dEMHjyYQYMGcezYMZYvX17pVw9VNJmLR86fP1/tFEeF2jhMW2QR29ouIkpEOoqIn2nTSqMATp48yWOPPYazszMPPvggBw4c4Pbbb2fNmjUcOnSIyZMn4+XlVSmURmUmLi6O+vXrY29vz9GjR9m1y5gC5sKFC9jb2/PAAw/w0ksvsW/fvlK3mZNbb72VZcuWERsbC8CVK1do3Lgx0dHRxMbGkpKSwpo1RfajpX///ixfvpzo6GhzewVFkXZ0dCQ+Pj7Pc1ntG3v37mXMmDG89957ODs7m+XV/EdmNOeoqCgLS1L26Ow21YBt27Yxd+5cfv75Z6ysrOjbty+TJ09m8ODBOrpsCRg8eDBz586lQ4cO5rzWYHzjfumll7CyssLW1rZYK2byazMnXl5evPbaa/Tt2xdra2v8/f1ZtGgRU6dOpUuXLjRv3hwPj6IM7I14enry9ttvM3DgQAwGA7a2tsyZM8f8UMuJk5MTPXv2xNvbmyFDhmRbTnro0CHuuOMOwKg4evToARiVol6ckpvM6cbqqFTzdQBUSq0uQv0rIjKuTCX6r/8FGJf3RouId2HlbzYHwLi4OObNm8fHH39MdHQ0tWvXZsKECbz44os0bdrU0uLly83oEFYdGTNmDM2aNSMpKYl77rmH/v3zjl96M/+909LSqFGjBm+99ZZ5SrMyUhIHwIJGHB2A8QX1B8wpomBOIlJctbsI+AJj2BGNiStXrvDWW2/x66+/EhkZSa9evXj++eeZOHEi9evXt7R4mpsEGxsbZs6caWkxKjW2trbUrVuXmJgYS4tS5hSkOF4TkX8KqqyUequwDpRSK4BopVQdjEmdNhdFMBHZanL802Bc0fLRRx/x448/Eh8fT69evczLKjWaiuaHH36wtAhVgri4OD7//HNGjx5Nz549LS1OmZGvcVxElhZWuShlgKMi8oSIjMUY/LDMuBk8x69evcr48eNxc3Nj7ty5DBo0iN27d7N582atNDSaKsKQIUPyXXRQFSnRqipTru+iMlgpNVkpdRuQVJL+8qM6e45fvnyZyZMn07x5cxYuXMj//vc/s8NWp05F8p/UaDQWRkRYunQp8fHxuZxMqzIlXY6b51IdpdRXSqknlFK9TFNTAIMxBijsCTRXSn1Xwj5vCs6ePcsrr7yCq6srs2bNYtSoUezfv5/PP/+8SmUV02g0RkaNGsWoUaPYsWMHP//8s6XFKRNKtBxXRL7O59R+wAe4D/BWSiVgzOB3CAgRkUJtIjcjBoOBv/76iy+//JI1a9aglOLee+9l6tSp1XJFiojoZcI3AeWVsqEqsmDBAkJCQhgzZgxBQUG888472RwwqxqFKg6lVJ7ryERkeh7HskVuU0q1wKhIOgJ3AEVWt0qpJUA/oKFSKhJ4U0TmF7V+VSA2NpaFCxcyd+5cTp06RePGjXnttdeYMGECLVu2tLR45YKdnR2xsbE4OTlp5VGNERFiY2Oxs7OztCiVgtq1axMSEsJzzz3HzJkz+eWXX/jggw8YM2ZMlfw/KDSRk1JqcpZdO4y+FUdE5NHyFKy4VCU/jri4OGbMmMGcOXNITk6md+/eTJo0iZEjR1KjRg1Li1eupKWlERkZWeRwG5qqi52dHS1atKi26VNLyvr165kyZQoHDhzA29ub+++/n5deegkbG8v4Y5fEj6PYGQCVUjWB9SLSr1gVy5mqojjCw8O55ZZbOHv2LA899BAvvPCCjjaq0dxkZGRkMH/+fBYvXszWrVvp1q0bgwcPpnHjxri6utKiRQtatmzJX3/9ha+vL87Ozjg6OpZLPLCKUhz1gb0i0rZYFUuAUmowMBuwxugD8n5+ZauC4jh79ix9+/bl+vXrrFmzxhyyQaPR3Ly88cYbLFy4kAsXLhRoF1JKYWdnR926dalXrx729vbUqlULa2trbGxsqFWrFvXq1cPKyooaNWpga2tLjRo1qFGjBlZWVhgMBuzs7LCxsTHv29jY8Prrr5e94lBKHQIyC1kDjYDpIvJFcToqLkopa+A4MACIBPYCY/LLOa6UkjfeeIPp03OZXioFMTEx9OvXj/Pnz7Nx40a9pFaj0WQjKSmJq1evcvLkScLDwzl37hy1a9emfv36xMbGEhcXR1JSEnFxccTFxZGYmEhycjIZGRlkZGSQmJhIXFwcBoOB1NRU0tLSSEtLIyUlBYPBgJWVFcnJyXkpp3JRHFnXgKYDUSKSXpxOSoJSqjswTUQGmfZfBRCR9/Ipb74QFxcXlFKICJGRkYSGhtK+ffs8+8nIyEBESExMJDExkdTUVGJjYzEYDCQkJJCQkMD58+f56quvePrppwFISUkhJSWF1NTUbJ8XLlxgyZIlDBw4kObNm5uPx8bGsmvXLnPSoPzi+mg0Gk15k/nMU0qRkZFBzZo1y3+qqqJQSt0DDBaR8ab9B4GuIvJUljITgAkATk5OAS4uLpYQVaPRaKoswcHBIiLFMp6UyIyvlFojIneWpG5ZYlr+Ow+qho2jJFy4cIFr167RokUL6tSpU3gFjUajKQZKqYIzeOVBSU30j5ewXnE4D2R1ZmhhOnZTcPz4cYYOHUqLFi3w8vKibt26jBkzhitXrlhaNI1Gc5NTIsUhIheLUk4pNVgpdUwpdVIpNSWP832UUvuUUummqanM437AJ8AtSqkjSqn7MXqjFyVHSJVn8+bNdOrUiW3btvH666/z448/8tJLL7F8+XJ8fX3Zv3+/pUXUaDQ3MUUxjrsD7wGeGB0AARARt0LqFboqyhQ2vQ7wIrBaRJabjrfDuJLLHfgcaA28JyJv5NdfdZmqOnz4MN27d6dVq1asX7+e5s2bm88FBwczfPhwrl+/zsyZM3nssceqpNepRqOpPCilzoqIS3HqFGXEsRD4CuOKqlswJlb6sQj1ugAnReS0iKRiDDcyLGsBEQkXkRDAkOP4cRE5ISLrRKQNcJibIKFTRkYGDz74ILVq1WLdunXZlAZAQEAA27Ztw8vLi8cff5y77rqLsLA8VydrNBpNUSl2pqmiKI5aIrIR4+jkrIhMwxh3qjCaAxFZ9iNNx4qFUqoLUAM4lce5apWPY/HixeZIuK1atcqzjIuLC9u3b2fmzJls3boVLy8vxo0bx4YNG8jIyKhgiTUazc1IURRHilLKCjihlHpKKTUCqF3OcgGglGoK/AA8IiKGnOerUz4OEeH999+nY8eOjBo1qsCyVlZWvPDCC5w+fZphw4axYsUKBgwYQKNGjRgwYADz588nISGhgiTXaDQ3G0VRHM8C9sAzQADwAPBwEeqValWUKZ/HWowpbHcVUrbKZ0jZtWsXR44c4fnnny+y3aJhw4asWrWK6Oholi1bxtChQ9m7dy/jx4/H0dERFxcXOnfuzNy5czl9+rQOc63RaPJEKdVAKfW3UuqE6bN+geXL62GilLLBaBzvj1Fh7AXuF5HDeZRdBKzJYhyvAfwB/C4iswrpxxo4HhAQ4FaVjePPPvssX3/9NVFRUdStW7fE7aSkpPDHH3+wfPlyjh49SnBwsPmco6MjHTp0wNfXl8aNG1OnTh2GDh1KmzZtLBaZU6PRWBal1GFgHXBFRN43rYCtLyKv5FsnP8WhlJpmsmcU1GGBZZRStwOzMMa4WiAi7yilpgNBIrJaKRUIrATqA8nAJRHxUko9gNEon1XJjBORA3n00R2YFhAQMLCqKo6MjAxatGhB9+7d+fXXX8u07fT0dA4cOEBQUBBhYWGEhIQQEhLCtWvXzCOQmjVr0rBhQ+6//358fX1xcXHB3d2dOnXq6HwKGk01Ryl1AuMq1n4ictFkItgiInnHaaJgxRGJ0Zci37rA4yLiUQqZS01maJKAgIDHqqri+Oeff+jXrx9Lly4t1L5RVhgMBrZv386xY8c4evQo69at49ixY9mms5RSNGjQADs7Ozw8PKhfvz5NmzbF3d0dZ2dnEhMTOXLkCH379qVXr17Uq1evQmTXaDRlh1JqP+AmIvVM+wq4mrmfFwXNT3wDOBbS5zeFCFRgWHRTbo/vMdpOYoHRIhJumqr6GuiMcanusyKypRBZqiwbN27EysqKQYMGVVifVlZW9OnThz59+gAwc+ZMbty4wfbt24mNjSUqKorLly9z5coVzp07R2RkJCdPniQqKipXEqaPP/4YADc3N5o0aYJSCl9fX1JTU2nTpg22trbY2dnRpEkTYmNjSU5OxtHRkbp165Keno6Dg4M5JLTBYKBhw4Y4ODhgZ2dHcnIy9erV08mANJryI6c7hGQNGpsX+SqO0uYHN9ke5pDFAVAptTpHWPTHMGq2tkqp+4APgNGYQpqIiI9Syhn4QykVmNfKKnIb4ascW7dupVOnThaPRVWrVi0GDBhQYBkR4dKlS1y9epXo6GjatWtHaGgoQUFBHDp0iKioKCIiIvjmm29IS0srM9nq1KmDn58ftWrV4uTJkzRu3JgOHTowduxYevToUaz8zZmjqop0nkxLS8Pa2hqgXJLxaDSlJEop1TTLVFV0QYXL0zheaFh0pdR6U5mdJmP6JYz5Pr4AdonID6ZyG4FXRWRPHv3YYDSOu1bFqaqUlBTq1q3Lk08+ycyZMy0tTpkhIqSkpJCcnMz58+e5cuUKdevWpX79+tSoUYOYmBiSk5MxGAzmEPbXr19HRDh8+DBJSUnUq1ePWrVqcerUKT777DMAPD09cXR0JCUlhSNHjpCSkgKAra0tderUoV27dly5coWoqCg8PDxo2LChud3mzZuzb98+Ll26REJCAs7OzqSlpdGkSROaNm2Kk5MTwcHBXLt2jcTERPz8/Khfvz4pKSlERERQp04dkpKSCAgIYPDgwXh5eZGUlMSlS5c4deoUJ06coGHDhtjb29OoUSOioqLYu3cvR44c4eTJk4DRnlS3bl1SU1O55ZZbePHFFwtM6LVr1y6WLl2KUorhw4dz9epVcz7vo0eP0qhRI5KSkrhy5QqHDx/Gw8ODBg0aMGDAAHr27KmVlKZQlFLBwGYgNotxvIGIvJxvnXJUHEUJix5qKhNp2j8FdAVGYhypjME4mtgPPCYiK3L0kRlWvW6rVq3anj17tlyupTzZvn07vXv3ZtWqVQwbNqzwChozN27c4LfffuPUqVPExcURERHBpUuXaNSoEYmJicTExBAREUHLli1JTU3lwIEDtGzZkjvvvJN69eoRHR1NWloa33//PdbW1rRp04akpCSaN29OgwYNCA8Px87Ojpo1a7J//358fHxwdnZm165d+QabtLa2zuaI2apVK9zd3Wnbti316tUjLS2NS5cuUatWLX7//Xeio6O55557mDhxIh06dOCPP/7A2tqakJAQ/vrrL8LCwqhRowZKKbOSzA8rKyvs7e3NPjwODg44OzvTqFEjWrdujbOzM4MHD8bd3Z1mzZrh6GiciY6KimLPnj2kpaVhb2/Pvn37sLW15fbbb8fT01OHtanmKKVOAt2ApUAr4Cxwr4jkG1G1KLGqnEQktgTClEZxXAM+whji5CxgC8wTkVX59VdVY1W9++67vPbaa8TExODk5GRpcao1mclrSktGRgarV68mOjqaZs2a4ezsjK2tLX5+fgDExcVx9epVDAYDbdq0ybfPhIQE3n33XebOncvVq1eznatRowZ9+/bl9ttv57HHHsNgMLB69Wpz6tDY2Fi6d++Ok5MT9vb25mkwgNjYWNasWcOBAweIjo5mx44d2NracuLEiWJfa4MGDWjZsiW1atWiadOmNG7cmGbNmmFra8v169dp2rQpLi4udO3aFWdn52K3r7E85ZJz3LRU6wDG5bF/SBGHKKWZqsrZh1JqBzA+v7SxUHUVx5AhQ4iIiCA0NNTSomgsxI0bN1i/fj2XLl3C1taWrl270qZNG2rVqlWm/SQnJ/PHH39w4cIFEhMTCQ4OJjw8nNtvv51+/fpha2tLcnIy/v7+xMTEsHr1ao4fP87Zs2c5d+4c169f5/z5gn14HRwcSE9Pp1GjRuYFDx06dKBNmzacOXOGEydOEB8fT+/evRkyZAh9+/aldu0KCUShyYfyUhwKuA14FAjEOJxZJCLHC6lXqAOgUupJwEdE/mcyjo8UkXuVUvYm2RKVUgOAN0SkT0H9VUXFISI0atSI4cOH8+2331paHI2mSCQkJJCcnEzdunWJiYlh7969bN68GSsrKy5cuEDdunWJj48nJSWFkydPcvDgQQDs7e3p0qUL6enp7N271zz1Vrt2berWrYuNjQ3Ozs44OjrSqFEj82iuSZMmNG/e3LzYItPOZWVlRe3atfWKu1JSEsVRqLuw6e3/b+BvpdQtGCPjTlJKHQSmiMjOfOqlK6WeAtbznwPg4awOgMB84AfTHNsVjDk3AJyB9UopA0al82BxLqqqcPbsWWJjY+ncuVh/M43GotSuXds8SmjatClDhw5l6NChBdZJTEykZs2a5ggFSUlJ7Ny5k507d3L+/HkSExPJyMjgypUrxMfHExQUxMWLF0lKSipUniZNmuDh4UGHDh3w9vbGwcGBAQMGULNmTc6dO8eOHTvYvn07Dg4OODk54ebmRseOHenYsSMODg6lvyE3IYUqDqWUE8b4VA8CUcDTGBMq+QHLANf86orIOoyu7FmPTc3yPRnI5fEmIuFAvl6L1YXMEVJlUhxpaWlERkbm8tXQaMqaZs2acffddxdYJiMjA4PBQHp6Okop0tPTzfvJycls3bqVEydOcPbsWebPn09qamqe7dSsWZNatWpx7do18zEbGxuGDBlC//79ueWWW/D29tar0IpIUQIU7cQYoXZ4phHbRJBSam75iHVzEBQURI0aNfDx8bG0KGYiIyPNARL1ahpNZSVzSbKbmxuursZ319TUVKKjozl27BghISHs3bsXe3t7HnvsMfz9/alZsybp6elERkayb98+li5dSlBQEL///jtgXAjQunVrGjVqRN26dc3TZnZ2dkRHRxMYGIiVlRVOTk507tyZxo0bW/IWWJSi2DjuFZGlOY6NEpFlhTZeQs9x07mOGL3H62D0bAw0jVDypCraOG677TauXbtGZZL7yJEjeHh4aKWhqfSICEePHqVDhw6laufcuXNs2rSJzZs3ExISQs2aNYmLiyMqKoqEhIR8HVnbtGlDjx49GDBgAE2bNsXDw4NmzZqZR+2pqak0adKEevXqVer/p/Iyju8TkU6FHcujXlFSx04COmYxjo8QkdEmw/o+4EEROWiaLrsmIvlmKqpqikNEaNCgAaNHj2bu3MozcDty5Eip/xE1moqiIn6vmYrg5MmTKKVIS0sz2022bNmSq3xOXx6ADh06MHPmTIYMGVKuspaEMjWOK6WGALcDzZVSn2U5VQdjGtnCMKeONbWXmTo265LaYcA00/flwBemVVwDgRAROQhQEj+Sys7p06e5du1apbJvaDSa3Nja2uLq6mqeEgPMCiA5OZljx45x9uxZQkNDiYuLIyYmho4dO5KSksK5c+c4ePAgkZGR3H777cyfP59HH33UUpdSdohInhvgizFh01nTZ+Y2EmOs9nzrmurfg3F6KnP/QeCLHGVCgRZZ9k8BDYHnMNpV1mMcebycTx8TgCAgqFWrVlKV+PnnnwWQffv2WVqUbISFhVlaBBEROXPmjCxcuNC8f/XqVWnQoIEYDAYREdmxY4cAEhERISIi165dk/r160tGRob07dtX9u7dW+S+Nm/eLHfccUeZyl9Srl69KnPmzLG0GAVy5swZ8fLyEhGRvXv3ytNPP11g+f3798vatWvN+7/99pu89957ZSJLZfm9FsaFCxcEY+hyefPNNy0tTjYwrnIt8Hmec8t3CYGIHBSR74A2IvJdlu1XEbmaX70ywgboBYw1fY5QSvXPQ8Yqmzo2ODiYmjVr4uXlZWlRKh1fffUVQ4YM4Y033qBfv35cunSJevXq0bRpU44cOQLAjh078Pf3Z8eOHYAxplOXLl2q/KqYa9eu8eWXX+Z5Lj29KAP9klOSnPWdO3c2xxHLjwMHDrBu3X+LK4cOHcqUKVOK3VdVpmnTpqSlpTFu3DjeeustXnkl3xxJVYJ8/8uUUpkG8f1KqZAs2yGlVEgR2i5K6lhzGZNdoy5GI3kksFVEYkQkCeOS3gJtKlWNoKAgfH19qVGjhqVFqVTEx8fz5ptvsnjxYmbMmMGiRYvMa+179OhhVhQ7duzg+eefz7bfs2dPczvLli2jS5cutGvXjm3btgHGaYVHHnkEHx8f/P392bx5c67+ExMTefTRR+nSpQv+/v789ttvucpcvHiRPn364Ofnh7e3t7n92rVr8/zzz+Pl5UX//v25fPkyAKdOnWLw4MEEBATQu3dvjh49ChhjRI0YMQJfX198fX3ZsWMHU6ZM4dSpU/j5+fHSSy+xZcsWevfuzdChQ/H09CQ8PBxvb2+zLB9//DHTpk0DoF+/fjz//PN07tyZDh06sHfvXkaOHIm7uzuvv/56nve7du3aTJ48GV9fX3bu3Mn06dMJDAzE29ubCRMmmCMJBwcHm+WcM2eOuf6WLVu48847AdizZw/du3fH39+fHj16cOzYMVJTU5k6dSq//PILfn5+/PLLLyxatIinnjJGHgoPD+fWW2+lY8eO9O/fn3PnzgEwbtw4nnnmGXr06IGbmxvLly/P5xdTdbCxseHrr7/G0dGRDz/8kF9++cXSIpWYgpbjPmv6vLOEbe8F3JVSrhgVxH3A/TnKrMY4/bUT49TWJhERUyiSl00e5KlAX+DTEspR6TAYDAQHBzN27FhLi1Igzz33HAcOHCjTNv38/Jg1a1a+562srFBKmYMIuri4mM/17NmTf/75h/Hjx3P69GlGjRrF119/DWB+6GaSnp7Onj17WLduHW+99RYbNmxgzpw5KKU4dOgQR48eZeDAgRw/nj0AwjvvvMOtt97KggULuHbtGl26dOG2227L5ij2008/MWjQIF577TUyMjLMTmqJiYl07tyZTz/9lOnTp/PWW2/xxRdfMGHCBObOnYu7uzu7d+9m0qRJbNq0iWeeeYa+ffuycuVKMjIySEhI4P333yc0NNR837ds2cK+ffsIDQ3F1dWV8PDwAu9vjRo1CAoKYvbs2QwbNozg4GAaNGhAmzZteP7553PFQ0tMTKRr167myMyenp5MnWp0tXrwwQdZs2YNd911F4888ghffPEFffr04aWXXsqzbw8PD7Zt24aNjQ0bNmzg//7v/1ixYgXTp08nKCiIL774AoBFixaZ6zz99NM8/PDDPPzwwyxYsIBnnnmGVatWAUYFvX37do4ePcrQoUO55557Crz2qkCNGjW4dOkSTk5OPProo/j4+ODp6WlpsYpNQfk4Lpq+xgA3RMSglGoHeGDMB14gUgrPcRG5qpT6BKPyEWCdiKwt8VVWMk6dOsX169e1YTwPHBwc+Oabb3j11Ve5dOkSoaGhTJ8+HXt7e3r06MF7773HmTNncHFxwc7ODhEhISGB4OBgunbtam5n5MiRAAQEBJgfttu3b+fpp58GjA+51q1b51Icf/31F6tXrzYnp0pOTubcuXPZVu4EBgby6KOPkpaWxvDhw83BDa2srBg9ejQADzzwACNHjiQhIYEdO3Zky+yYGWpj06ZNfP/994BxJU7dunVzBTsE6NKlSzbDbEFkenD7+Pjg5eVF06ZNAWOSrYiIiFyKw9raOpsT3ubNm/nwww/Nodq9vLzo3bs3165dMyf9evDBB/njj9yPgLi4OB5++GFOnDhhXn1UGDt37jSnS37wwQd5+eX/InkPHz4cKysrPD09iYqKKtL1VwXs7e0JCwvD39+fYcOGsXPnTho2bGhpsYpFURwAtwK9lVL1gb8wPsxHY7Q/FIiU0HPcdO5HjOFNqh2V0WM8LwoaGZQnQ4cOpWPHjvz+++8EBQUxc+ZM3njjDdzd3bl27Rq///473bt3B4yKYeHChbi4uGQLlpeZ2Mna2rpYtgERYcWKFbRvn3/ggj59+rB161bWrl3LuHHjeOGFF3jooYdylVNKYTAYqFevXqlGbllHOzY2NhgM/+Uzy+nhn3ndVlZW2ZJbWVlZ5Xkf7OzszJF1k5OTmTRpEkFBQbRs2ZJp06YVK4LAG2+8wS233MLKlSsJDw+nX79+Ra6bF1nlz5wyqy64urry9ddfc99993HnnXeycePGKhX+pCiWRGWyM4wEvhSRUUCRLLpKqcFKqWNKqZOm5CA5z9dUSv1iOr9bKeViOt5FKXXAtB1USo0oxjVVeoKCgrCzs6uSQ9TyJiEhgcy8Ko6OjnTo0IH4+Hjz+W7dujF79myz4ujevTuzZs3KZt/Ij969e7N48WIAjh8/zrlz53IpiEGDBvH555+bH1T79+/P1c7Zs2dp3Lgxjz/+OOPHj2ffvn2AcQoycy7+p59+olevXtSpUwdXV1eWLTP6y4qIOehf//79+eqrrwCjYTouLg5HR8ds15uTxo0bEx0dTWxsLCkpKaxZs6bQ6y4qmUqiYcOGJCQkmK+lXr161KtXj+3btwOY72FO4uLiaN68OZB9Oqqga+rRowc///yzud3evXuXybVUBUaPHs3KlSvZu3cvo0aNKtOMmeVNkRSHKUT6WCBzusi6gPKZlTJTxw4BPIExSqmcT0pz6liMNowPTMdDgc4i4gcMBr42Gc+rBcHBwfj5+ZkDvmn+Iy0tjYkTJzJq1Cjeeust/vjjD5599lnz+Z49exIREWEerXXv3p3Tp08XmEUvk0mTJmEwGPDx8WH06NEsWrQoV8rZN954g7S0NDp27IiXlxdvvPFGrna2bNmCr68v/v7+/PLLL2b5HBwc2LNnD97e3mzatMlsK1i8eDHz58/H19cXLy8vs8F99uzZbN68GR8fHwICAggLC8PJyYmePXvi7e2dpy3B1taWqVOn0qVLFwYMGICHh0cR72zh1KtXj8cffxxvb28GDRpEYGCg+dzChQt58skn8fPzy/ft/+WXX+bVV1/F398/2+jmlltuISwszGwcz8rnn3/OwoUL6dixIz/88AOzZ88us+upCgwfPpy5c+fyxx9/mPOuVAkKW68L9MFoxH7FtO8GfFaEet2B9Vn2X8WY/jVrmfVAd9N3G4z2FJWjjCvG4Io2BfUXEBBQ3OXLFiEjI0Nq164tTz31lKVFyZPKsi4+px9HVcDBwcHSItx0VJbfa2mZMWOGADJ48GCJi4ur0L4pgR9HUcKqb8Vo58jcPw08UwSd1ByIyLIfiTG7X55lxGhMjwOcgBilVFdgAdAaY+iRXBO0WVLH0qpVqyKIZHmOHz9OQkICAQEBlhalUlOvXj2z0Vmjqe689tprXLt2jZkzZ1K3bl0uXrxIkyZNLC1WvhQ6VaWUaqeUmqeU+ksptSlzK2/BRGS3iHhhTB71qlLKLo8yVc4BMDg4GKj8hnFLUxUVR2aub42muCil+Pjjj/nmm28AaN26Na+88kqJnDIrgqLYOJYB+4HXgZeybIVRGgdAMyJyBEgAvKkGbNy4kXr16pXp3HRZI9VsBYumelIdf6fjx48nNDSUdu3a8eGHH+Lj48Pq1asr3bUWRXGki8hXIrJHRIIztyLUMzsAKqVqYPTRWJ2jTKYDIGR3AHTNNIYrpVpj9B0JL8oFVWYMBgNr165lyJAhldYwbmdnR2xsbKX7oWo0WRFTPg47u1wTEVUeLy8vQkJCePXVVzl69CjDhg3DysqKefPmVZoRSFHCqk8DooGVQErmcRG5UmjjSt0OzOI/B8B3sjoAmqaffgD8MTkAishppdSDwBQgDWMujukisqqgvqpCWPVdu3bRvXt3fvrpJ8aMGWNpcfJEZwDUVBXs7Oxo0aJFtc45npqayrRp03jvvfcAY7Kp9u3bM3z4cJ5++mlq1apV6j7KKx/HmTwOi4i4Faej8qYqKI5XX32Vjz76iMuXL1O/fn1Li6PRaKoIIsKvv/7KN998w/r167Odu//++83xyQICAiiuvbdcFEdpKGUGwFcx+nlkAM+ISPa7lQNvb28JDQ0t82soLpn302AwYGVlhcFgIC0tja1btzJo0CBuueUWNm0q97UFGo2mmpKQkMDvv//Ovn372Lx5s3nBTVY8PDxwd3c35wOZNm0aTk5OxMfHExgYSK1atWjRogXOzs7Y29uXy4jDHngBaCUiE5RS7kB7ESnQZbWUGQA9gSUYk0E1AzYA7aSADIBKKRkyZAjR0dHUr1+f8PBwEhIScHJyYuTIkRw8eJCjR48SGRlJt27diImJ4cqVK0RGRtK9e3dSUlLYt28fvr6+XLx4EScnJ44dO2Z2yLG3t8fZ2Zn09HTS0tKyxc7JK+NXfhw+fFh7jGs0mjIjPT3dHBjzgw8+4OjRo1hbW+Pp6cmhQ4eK0kS5KI5fgGDgIRHxNimSHWL06i6oXndgmogMMu2/CiAi72Ups95UZqfJGH4JaITRvmEum7VcAf3leSFKKUSEWrVqcePGDcA4N5p1Dr9+/frm4HIBAQFcvXqV06dPm8/Xrl2bOnXq0L9/f2xsbLCxsWHHjh1cvnwZNzc3c4gGpRSbN2/G0dGRxo0bk5KSQtOmTYmIiGDKlCmMHj36pk5wr9FoKp7r168TFxfHmjVrcHBw4NSpU4SEhBAZGcmtt97Khx9+WC6KI0hEOiul9ouIv+nYQRHxLaTePcBgERlv2n8Q6CoiT2UpE2oqE2naP4XRSXAasEuMgQ5RSs0H/hCR5Tn6MDsAOjk5BWQNwa3RaDSawgkODhYRKVYGtKKsCU1VStXCGN4cpVQbsqyusiQiMg+YB1XDOK7RaDSVDaXUjeLWKYrimAb8CbRUSi0GegKPFKFecRwAI3M4ABalrkaj0WgsQKHDExH5C2NI9XEYDdadRSR3zs3clNgB0HT8PlPYdVfAHdhThD41Go1GU84UOuJQSm0Ukf78F1I967F8kdJlADxsynkeBqQDTxa0okqj0Wg0JeZycSvkaxw3eXXbA5uBfoAynaoD/CkilSrYUl42Du0FrcnKzeBprNEUl5I4ABY04pgIPIfRjyKY/xTHdeCLkghY0URGRuLo6IiLiwtKqcIraKotmbGNIiMji5y/W6PR5E2+ikNEZgOzlVJPi8jnFShTmZGcnKyVhgYw+vM4OTlx+XKxR+UajSYHRTGOf66U6qGUul8p9VDmVppOC8tFbipzr1IqTCl1WCn1U2HtFXCuNKJqqhH6t6DR5I1SqoFS6m+l1AnTZ4HB9IqSyOkH4GOgF8akSoFAibMQFSUXuSmsyatAT1Myp+eK0J5Go9Foio8dxmgdG0XEHdho2s+XovhxdAY8peyiIXYBTppS0KKU+hkYhnEFVSaPA3NE5CqAiEQX1h7GXOgaTZVl165dNGrUiMaNGxMXF0ezZs30KElTEdTA+AzuZ9r/DtgCvJJfhaK4mYcCZZn8Nq9c5M1zlGkHtFNK/auU2pXfVJQp5MiPQCc9d62pysyYMYPu3bvTtm1bHB0dadGihTm6skZTziQAjUXkomn/ElBgUL2iKI6GQJhSar1SanXmVkpBC8MGo9NfP2AM8I1Sql7OQqaQI68Av1WVnONZqV27drHrXLt2jS+//LIcpMmbadOm8fHHH5dpmz169AAq/loqmri4uCKVO3fuHFOnTs3z3KpVqwBITExk9erVWpFoyoNsPyrT7FKBM0xFURzTgOHAu8DMLFtJKUo4kUhgtYikicgZjOHZ3YvYXrWmoIetiFSJB8uOHTuA6qE4goOD8wypv2zZMurVq2e+1oJ4//338z139913k5yczIwZMxg2bBjW1tZs2bKl1MrcYDCwdOlS3nnnHU6dOlWqtjTVgiilVFMA02dB5gHjw6YiN4yjidOAK8a5tYOAV44yg4HvTN8bYpzaciqovYCAAMlJWFhYrmOW4IcffpDAwEDx9fWVCRMmSHp6uoiIODg4FHheROS7774THx8f6dixozzwwAMyevRosbOzE19fX3nxxRflzJkz0q5dO3nwwQfF09NTwsPDZebMmeLl5SVeXl7y6aefiojImTNnxMPDQ8aPHy+enp4yYMAASUpKylPet99+W9zd3aVnz55y3333yUcffZSvnAW1m5CQILfffrt07NhRvLy85Oeff8523TmvJVNWEZH/+7//k1mzZuUp34gRI+S1116T3r17S8uWLeXvv/8u8t+ipL+JK1euSHJysnk/LS1NgoODBZCnnnpK0tLSRERk79690qdPn8w3Nvnmm28KbNdgMEjfvn0FkPDwcFm1apW0atXKXB+Q3r17S6NGjbIdAyQjI6NE1yIiMm7cuGxtpaamlrgtTdUGCAI+AqYYd5kCfCgFPcfzPQHxGJ39cm7xwPWCGi1sA27HOIo4BbxmOjYdGGr6roBPMBrMD2HMRV5ge4UpjmeffVb69u1bptuzzz5b6B8lLCxM7rzzTvM/5hNPPCHfffediBgfoAWdDw0NFXd3d7l8+bKIiMTGxsqZM2fEy8vL3P6ZM2dEKSU7d+4UEZGgoCDx9vaWhIQEiY+PF09PT9m3b5+cOXNGrK2tZf/+/SIiMmrUKPnhhx9yyZtZPzExUeLi4qRNmzby0Ucf5StnQe0uX75cxo8fb2772rVr5uvOlD3zWs6cOSP+/v4iIpKRkSFubm4SExOT5z1t27atWZn9+uuvMm7cuEL/DpkUR3GkpKTIu+++KwkJCQJIv379RMSoRDw8PHI9yA0GQ65jK1asKLCP9evXC5BLSeZsJ6/t5MmTYjAYiqxA4uLiJCwsTH7//fdcbX3wwQd51vn777+zKXRN9QM4AThhXE11AmPivAZSwDO3IAdAx/zOlRYRWQesy3FsapbvgjHr4AtFba9z5xKvEC5XNm7cSHBwMIGBgQDcuHEDZ2fnIp3ftGkTo0aNomHDhoAxUf3169dz9dG6dWu6desGwPbt2xkxYgQODg4AjBw5km3btjF06FBcXV3x8/MDjAmrwsPDc7W1bds2RowYgb29PQBDhw4tUM4+ffrk266Pjw+TJ0/mlVde4c4776R379753icXFxecnJzYv38/UVFR+Pv74+TklKtcUlIScXFxPP/884AxrEy9evXybDM2NpYrV67Qtm1bAM6ePZtv+Jn09HRsbP77d/j888955plnAJg/fz4AW7ZsQUS45557OHr0aK427r777lzHTp48me81A/z555/UqFGD//3vf9mOv/DCC3zyyScF1s28Lj8/P+666y5ee+01bGxsSE5OJioqigsXLgDG1VorVqxg165d+bb1yiuvYGNjwyeffML+/ftp1KgRP/74Iw8++CAAzz//PMOGDTPbXAAyMjK4du0aH374Ibt27aJHjx689957udqOi4ujVq1a1KhRI9vxmJgYIiIi8Pf3L/A6NeVOnIjEAgXGH8xKUZbjVgtmzZplkX5FhIcffjjPf6iinC8KmUqiMGrWrGn+bm1tbc6IWBTykzM8PDzfdtu1a8e+fftYt24dr7/+Ov3798/XCAwwfvx4Fi1axKVLl3j00UfzLBMWFkZAQADW1tYAhISE4O3tTVpaGm+++SZJSUkYDAY+++wzzpw5A8CFCxewt7cnJiaGmJgYLl68yNKlS3F2diYqKoqvv/46myLw8vLi8OHD5v2sNgArq/zNgitXrsx17JVXXmH27Nm88cYbGAwGJk2alK3dTz/9FE9Pz2z3EGDmzJlMmzaNVatW8dBDD9G9e3dOnjyZp+f7gQMHOHDgADNmzMhXtrwYMWJENpknT54MGG0ueSmt3377DaUULVu2JCIiItf5rVu3mu01bm5uLFy4kAULFvDdd9+ZywwbNoyAgIA8fwetWrXi559/pk2bNpw/f55///2XO+64gxYtWpCWloaNjQ22trakpKRgZ2dXrGvVlDEFDUeq0lZZbRyHDx+Wtm3bSlRUlIgYp5vCw8NFxDhlU9D5zKmqzCmb2NhYiYmJkVatWpnbzzl1FRwcLD4+PpKYmCgJCQni5eVlnqrKWu6jjz6SN998M5e8mfWTkpLk+vXr5mmh/OQsqN3z58/LjRs3RETk999/l2HDhpmvW0RyXUtKSoq0a9dOXF1ds9l5srJgwQKZMmWKef/OO++UzZs3y+effy4TJ06UiRMnytChQyUsLEz27t2ba/vjjz+KNA1U2NaqVStp3ry5PPHEE/Lqq6+aj7/wwgsyefJkiYiIyLeum5ub/P777+Lg4CCA/Prrr3leaybx8fFiMBjEYDDIli1b5MqVK/Lnn3+WSO6ffvpJ7r//fjl16pSIiKSmppbJ/aiozdPTUwC544475LnnnhODwSAxMTH5/l40hYMxWnmxnrelHnEopWYDdUTkEaXUQDHm7yiszmBgNsZw69+KSJ7LSpRSdwPLgUARqZLp/Tw9PXn77bcZOHAgBoMBW1tb5syZQ+vWrQs97+XlxWuvvUbfvn2xtrbG39+fRYsW0bNnT7y9vRkyZAhPPvlktv46derEuHHj6NKlC2B8i/f3989zWiovOnXqxOjRo/H19cXZ2dk8NZWfnE2a5O/ic+jQIV566SWsrKywtbXlq6++ynbeyckp27V89NFH3HLLLdSrV888osirza5du5r3Q0JCsLGxYcOGDbz88svm6ZDExMQiXW9+zJkzh19++YWtW7cCcPnyZTKXfK9evZq77rorW/m3334712ikdevWnD17NtcI5vTp0+b6W7ZsoW/fvgXKknXZdmbZQYMGkZ6ejlIKKysrvvzyS4KDgzlx4gQJCQm0bduWHj16YGNjQ8+ePbNNB40ZM8b83dbWFoPBwNy5c7njjjuYNGkSa9eu5bPPPqNNmzbcdtttiAhnzpzBxcWFY8eOYW1tzeTJk+natSt9+/Y1jwBDQkI4cuQIZ8+exc3NjeDgYF599VUaNGjAuXPn8PD4L6D2sGHDmDZtGq1bt2by5MksXLiw8D8KxhEnwNq1xiwPO3bsYM8eY6qeWrVq8euvv3L9+nW6deuGs7OzHpmUF8XVNDk34FNgqun7+0Uob43RKO7Gf6uqPPMo5whsBXZhTB5VJUccmqKTkZEhvr6+cvz48SKVj4+PN48kZs6cKQMHDpQHH3xQZs+eLXv37pXz589LWlqaREdHS3R0tFy9etU84li1apWMHDlSpk6dKk888YSkp6eLwWDI1cfWrVtl7969IiKyfft2uXjxYrGuJ/NN+PTp07J8+fJsb88rV64sclsVRUZGhqSkpJRL2xs2bMi2Mi1nvxkZGRISEiJpaWkSHx8vIiJRUVFy7tw52bFjR65RR2Hb7bffXi7XUd2gBCOOfPNxFBWl1Dv8F+vkOxG5v5Dy3YFpIjLItP+qSYG9l6PcLOBv4CXgRSlkxJFXPo4jR47QoUOHYl2PxjKEhYVx5513MmLECGbOzO0mJCIcO3aM5s2b4+hoXLcRERFBVFRUtnLNmjUjIyODRo0a5XrbFBEOHjyIj49PviMaTeXl9OnTuLi4YGVlRUpKCqtXr+bee+/FxcUl3xF1aZ9vNwMlycdRFAfAwngT4wjiS6DAKLYmCg05opTqBLQUkbUUgFJqglIqSCkVpEOOVD0yMjLMznOenp6cPn2amTNncuXKlVyrn2JjY0lISODYsWPmYwkJCVhbW9OpUydcXV1p3bo1zZo1o2XLlnlOUSilqFmzplYaVRQ3NzfzdGDNmjUZNWqUeRptzZo1jB07FoCXXnrJXKdx48bs37/fIvJWZ0pk41BK3QqMBa5hjGW1G5gvIimlFUgpZYXRh2NcYWXFGHJkHhhHHKXtW1M+xMbGcvHiRRwcHMz5UY4dO0Z8fDwA9vb2ODo60qxZM65evWp+e/Tw8CAtLQ1ra+tsb5RBQUHY2dmZlYuVlVWeS3c1Nw933HEHgwcP5r333qNly5aEhobyxx9/EB0dTdeuXUlNTbW0iNWKkhrHF2AMdW4LdMQYksQLaFuEuoWFHHEEvIEtpsigTYDVSqmhhU1XaSoP8fHxZGRkoJQyL4tNTk6mVq1axMbGZlsKnJSURFJSUq5pp7x8JTLR6YA1ObG2tqZlS+OjZc6cObi5GQNmp6WlAcYQN3Xq1Mm1iCE+Pt48/akpGiVVHGdFZJXp+7Ji1t0LuCulXDEqjPsAs11EROIwhhkBQCm1hSLYOPJDRHRo6jwQETIyMsxOb5n76enp1KxZM997ZjAY8vRlEJFcU0l5ERkZWXrhs9CpU6cil9Xz3TcPrq6uLFu2jFGjRgHZk3g5OjqaR7ujR4/ml19+AeCRRx7BwcGBQYMGERwczNixY2nbti2pqam5nBdvdkpkHFdKzQCuALOkBA0opW4HZmFcYbVARN5RSk3HaN1fnaPsFkpoHD9z5gyOjo44OTlp5YHxwXnu3Dnq1q1boEdzgwYNcHZ25tKlS4Dx7d7R0dHsfJY13HfWKaPi4OjoSJs2bbhy5Qrnzp3D3d2dGzducP36dRo3bkytWrU4d+4ciYmJ1K5dm9TUVNzc3IiIiKBmzZpERUXRqVOnAh3ycl57bGws8fHxOuf4TYKIsGzZMkaPHl1mbTo4OLB//37c3fOLuVr1KIlxvKSKYwXgA9QBgoEDwAERKe7oo8zIS3GkpaURGRl5009rpKamEh8fT0JCQoX1WaNGDWrWrEm9evWK/HAvb+zs7GjRogW2traWFkVTgVy7do309HSuX7/O/PnzcXd3Z+3atTRt2pQjR46glOLvv/8uVpubNm3i2LFjLFmyhIkTJ3L//QUuJq3UlERxlNaHoxYQgNGQ/XFp2irtlpcfx83K5s2bzZFot2/fXuA691GjRomtra189dVXEh8fL19++aUopbL5Gixfvlzuvfde6d27d57tbN68WXvuaqodkZGRMm7cONmyZYuMHDlSbr311nz/lzp37myOkFzVwBJ+HJWFvEYc1ZGDBw9ia2tLnz59+Oeff8xeyU8++ST//PMPdevWNScQysvo1759e86cOUNKSsEL4E6fPo2joyM5E2RFR0cTFxeHu7s7169fJz4+nubNcyZw1GiqJ5cuXaJp06b5ns8rqkBlp8JHHCXdMObbOIYxV/iUPM6/gDGkegjGUL+tC2vzZhhxbNiwIdebTs68CvltmfGQ0tLSJCEhwcJXotFUbcLDw8VgMMhzzz2X63+tqkFVGHEopawx5uIYgNH5by8wRkTCspS5BdgtIklKqSeAfiJSoIWruo841q5dy5133lmiuikpKXpViEZTTly6dIkePXqYl50vW7aMe+65x8JSFR1LeY4Xly7ASRE5LSKpwM/AsKwFRGSziCSZdndh9PW4afnkk09KrDTefvttrTQ0mnKkSZMmZH1pfeSRRywoTcVgCcVRaMiRHDwG/JHXieoeciQmJgYXFxdznoRMatSoYV57DjB27Fieeuoppk2blq3c999/z2uvvVYRomo0NzUNGjRgw4YNgDEUTln7K1U2KnUiJ6XUA0BnIM+401INQ45cvXqVl19+mbCwMHbs2GE+HhISQrt27bCysjL7VbRs2RIrK6tsYcbffPNNS4it0dz09O/fn9mzZ/Pss89yzz33sH379mxZJasTlriqwkKOAKCUug14DegrZRADq7KTkZHB7t276dmzp/mYu7s7X375Jf369cv2A8z0Q+jevXuFy6nRaPLnmWeeISMjgxdeeIGXX36ZBx98sFqmxrWEcdwGo3G8P0aFsRe4X0QOZynjjzGB02AROVGUdquycXzTpk0MHz7cHAZh3LhxTJkyhfbt21tYMo1GU1zS0tKwt7cnPT0dMKYuLmgJr6WpEsZxEUkHngLWA0eApSJyWCk1XSk11FTsI6A2sEwpdUAptTqf5qo0IsKkSZPo378/8fHx2NnZcfLkSRYuXKiVhkZTRbG1tWXx4sXm/RYtWpCUlFRAjaqHdgC0ECLC2LFjWbJkCQBz585lwoQJOqaWRlNNMBgMtG3bljNnzrBnzx5zGubKRpUYcWiM/PPPPyxZsoTRo0eTnJzMxIkTtdLQaKoRVlZW5tzoBaUIqIpoxWEhpk2bRtOmTVm4cCE1a9a0tDgajaYcaNu2LTY2NoSGhlpalDLFIopDKTVYKXVMKXVSKTUlj/M1lVK/mM7vVkq5WEDMcmPevHn8888/TJkyhVq1allaHI1GU07Y2toSGBjIhx9+yB9/5OmOli8vvvgiU6dOLSfJSkdlDTkyCegoIv9TSt0HjCgs5EhAQIAEBweXo+TFR0RIT083b2fPnmXmzJksWrSIAQMGsHbtWh3iW6Op5oSFheHn50daWhqtW7fGw8ODevXq0bBhQ5ydnXF0dMTBwQF7e3vzp729PX369AHKPwFZheXjKA1Kqe7ANBEZZNp/FUBE3stSZr2pzE7T8t1LQCMpQFhra2tp1aoVycnJZGRkZLaDlZVVts+8jmWS2XzWbgo7ZjAYzIohLS0tm6LITHaUFTs7OyZMmMCHH36op6g0mpuEqKgovv32W8LCwjh69CgJCQlcvnyZq1evFlpXKw5AKXUPRv+M8ab9B4GuIvJUljKhpjKRpv1TpjIxOdqaAEww7XoD1Wsi0bI0BGIKLaUpKvp+li36fpYd7UWkWEnXq7Q/fNaQI0qpoOJqTU3+6PtZtuj7Wbbo+1l2KKWK7cdgCeN4UUKOmMuYpqrqArEVIp1Go9FoCsQSimMv4K6UclVK1QDuA3J6hq8GHjZ9vwfYVJB9Q6PRaDQVR4VPVYlIulIqM+SINbAgM+QIxkxUq4H5wA9KqZPAFYzKpTDmlZvQNyf6fpYt+n6WLfp+lh3FvpfVJuSIRqPRaCoG7Tmu0Wg0mmKhFYdGo9FoikW1UByFhTDRFA+lVLhS6pAppH3VCTlcSVBKLVBKRZv8kTKPNVBK/a2UOmH6rG9JGasK+dzLaUqp86bf5wGl1O2WlLEqoZRqqZTarJQKU0odVko9azperN9nlVccphAmc4AhgCcwRinlaVmpqgW3iIifXitfIhYBg3McmwJsFBF3YKNpX1M4i8h9LwE+Nf0+/URkXQXLVJVJByaLiCfQDXjS9Lws1u+zyisOoAtwUkROi0gq8DMwzMIyaW5iRGQrxtWAWRkGfGf6/h0wvCJlqqrkcy81JURELorIPtP3eIzJ9JpTzN9ndVAczYGILPuRpmOakiPAX0qpYFNYF03paSwiF03fLwGNLSlMNeAppVSIaSpLT/uVAFPUcX9gN8X8fVYHxaEpe3qJSCeM039PKqX6WFqg6oTJmVWvgy85XwFtAD/gIjDTotJUQZRStYEVwHMicj3ruaL8PquD4ihKCBNNMRCR86bPaGAlxulATemIUko1BTB9RltYniqLiESJSIaIGIBv0L/PYqGUssWoNBaLyK+mw8X6fVYHxVGUECaaIqKUclBKOWZ+Bwaiow6XBVnD6DwM/GZBWao0mQ84EyPQv88io4x5JOYDR0TkkyynivX7rBae46bleLP4L4TJO5aVqOqilHLDOMoAY0ian/T9LB5KqSVAP4yhv6OAN4FVwFKgFXAWuFdEtNG3EPK5l/0wTlMJEA5MzDI/rykApVQvYBtwCMhMGPR/GO0cRf59VgvFodFoNJqKozpMVWk0Go2mAtGKQ6PRaDTFQisOjUaj0RQLrTg0Go1GUyy04tBoNBpNsdCKQ6PRaDTFQisOzU2PUsopS4juS1lCdicopb4sh/4WKaXOKKX+VwZtfWSS+cWykE2jKQoVnnNco6lsiEgsRocylFLTgAQR+bicu31JRJaXthEReUkplVgWAmk0RUWPODSafFBK9VNKrTF9n6aU+k4ptU0pdVYpNVIp9aEp4dWfpvg/KKUClFL/mCILr88RHiO/fhYppb5SSu1SSp029btAKXVEKbXIVMbaVC7U1Ofz5XrxGk0BaMWh0RSdNsCtwFDgR2CziPgAN4A7TMrjc+AeEQkAFgBFDddSH+gOPI8xbtCngBfgo5Tywzgiai4i3qY+F5bVRWk0xUVPVWk0RecPEUlTSh3CGBftT9PxQ4AL0B7wBv42xpLDGmPY76Lwu4iIqe0oETkEoJQ6bGr7H8BNKfU5sBb4q0yuSKMpAVpxaDRFJwVARAxKqTT5L9CbAeP/kgIOi0j3krZtaisly3EDYCMiV5VSvsAg4H/AvcCjJehHoyk1eqpKoyk7jgGNlFLdwZj3QCnlVRYNK6UaAlYisgJ4HehUFu1qNCVBjzg0mjJCRFKVUvcAnyml6mL8/5oFHC6D5psDC5VSmS97r5ZBmxpNidBh1TWaCsa0UmpNWSzHNbU3jYpZQqzRAHqqSqOxBHHAjLJyAAQeALQvh6bC0CMOjUaj0RQLPeLQaDQaTbHQikOj0Wg0xUIrDo1Go9EUC604NBqNRlMs/h9Zp5gUsoaffQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 5 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import numpy as np\n",
    "from matplotlib.pyplot import *\n",
    "from pylab import *\n",
    "import matplotlib.image as mpimg\n",
    "#from urllib  import urlopen #python 2.7 \n",
    "from urllib.request  import urlopen #python 3.0 \n",
    "\n",
    "ShotNo = 39187\n",
    "\n",
    "\n",
    "#Create a path to data\n",
    "baseURL = \"http://golem.fjfi.cvut.cz/shots/\" #remote access\n",
    "#baseURL = \"file:///golem/database/operation/shots/\" #local access\n",
    "#dataURL = urlopen(baseURL+ str(ShotNo) + '/' + diagnSPEC + '.npz')\n",
    "\n",
    "\n",
    "fig=figure(1);subplots_adjust(hspace=0.001)\n",
    "\n",
    "sbp1=subplot(511)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/Diagnostics/BasicDiagnostics/U_Loop.csv')\n",
    "obj1=np.loadtxt(dataURL, delimiter=',')\n",
    "ylim(0,20);yticks(arange(0, 30, 5))\n",
    "title('#'+str(ShotNo));ylabel('$U_l$ [V]')\n",
    "setp(xticklabels, visible=False)\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1], 'k-',label='Loop voltage $U_l$' );legend(loc=0)\n",
    "\n",
    "sbp1=subplot(512, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/Diagnostics/BasicDiagnostics/U_IntBtCoil.csv')\n",
    "obj1=np.loadtxt(dataURL, delimiter=',')\n",
    "yticks(arange(0, 0.6 , 0.2));ylim(0,0.6)\n",
    "ylabel('$B_t$ [T]')\n",
    "setp(xticklabels, visible=False)\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1], 'k-',label='Toroidal mag. field $B_t$');legend(loc=0)\n",
    "\n",
    "sbp1=subplot(513, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/Diagnostics/BasicDiagnostics/U_IntRogCoil.csv')\n",
    "obj1=np.loadtxt(dataURL, delimiter=',')\n",
    "yticks(arange(0, 4.5, 1));ylim(0,4.5)\n",
    "ylabel('$I_p$ [kA]')\n",
    "setp(xticklabels, visible=False)\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1]/1000, 'k-',label='Plasma current $I_p$');legend(loc=0)\n",
    "\n",
    "sbp1=subplot(514, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/Diagnostics/BasicDiagnostics/U_LeybPhot.csv')\n",
    "obj1=np.loadtxt(dataURL, delimiter=',')\n",
    "yticks(arange(0, 0.14 , 0.03));ylim(0,0.14)\n",
    "ylabel('Intensity [a.u.]')\n",
    "setp(xticklabels, visible=False)\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1], 'k-',label='Whole spectrum radiation');legend(loc=0)\n",
    "\n",
    "sbp1=subplot(515, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/Diagnostics/Interferometry/ne_lav.csv')\n",
    "obj2=np.loadtxt(dataURL, delimiter=',')\n",
    "yticks(arange(0, 0.8 , 0.2));ylim(0,0.8)\n",
    "ylabel('$n_e$')\n",
    "xticks(arange(0, 30, 5));xlim(0,20)\n",
    "xlabel('Time [ms]')\n",
    "plt.plot(obj2[:,0], obj2[:,1]*1e-19, 'k-',label='electron density $n_e$');legend(loc=0)\n",
    "\n",
    "xticklabels = sbp1.get_xticklabels()\n",
    "setp(xticklabels, visible=True)\n",
    "\n",
    "savefig('basicgraph.pdf')\n",
    "savefig('basicgraph.jpg')\n",
    "\n",
    "show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "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"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}
