Tujuan: Memperkenalkan cara penentuan episenter yang sederhana yaitu menggunakan metode lingkaran
\n",
"
Keluaran: Peserta dapat memahami dasar secara sekilas kemudian mengimplementasikan dalam bentuk kode Python
\n",
"
Sesi Materi
\n",
"
Waktu/Tempat: Sabtu, 25 September 2021/ Zoom
\n",
"
\n",
"
"
]
},
{
"cell_type": "markdown",
"id": "4479c0c8",
"metadata": {},
"source": [
"## Sekilas Dasar\n",
"\n",
"Penentuan episenter metode lingkaran memanfaatkan minimal 3 stasiun seismik. Data yang dibutuhkan dari setiap stasiun adalah waktu tiba gelombang P dan gelombang S yang kemudian diselisihkan untuk mendapatkan waktu S-P. Selisih waktu S-P ini dipakai untuk memperkirakan jarak yang sudah ditempuh yang pada contoh ini akan mengaplikasikan metode Tau P. Metode ini sering juga disebut *circle and chord* karena memanfaatkan perpotongan lingkaran dan garis perpotongannya.\n",
"\n",
"![circle_chord](figures/circle_and_chord.png)\n",
"\n",
"Dalam gambar (Havskov et al, 2002) di atas S merupakan stasiun seismik dan d merupakan jarak yang didapatkan dari pencocokan nilai S-P dengan jarak."
]
},
{
"cell_type": "markdown",
"id": "9a5691e7",
"metadata": {},
"source": [
"### Data\n",
"\n",
"Data yang akan kita gunakan adalah data waktu tiba gelombang P dan gelombang S dari gempa Mamuju M6,2 yang direkam di stasiun TATO: Taiwan, PMG: Port Moresby, dan COCO: Cocos Island. Data waktu tiba ini dikonversi kedalam format `UTCDateTime` `obspy`:"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "ffcdaf13",
"metadata": {},
"outputs": [],
"source": [
"from obspy import UTCDateTime\n",
"\n",
"TATO_P = UTCDateTime(\"2021-01-14T18:34:08\")\n",
"TATO_S = UTCDateTime(\"2021-01-14T18:38:53\")\n",
"\n",
"PMG_P = UTCDateTime(\"2021-01-14T18:34:15\")\n",
"PMG_S = UTCDateTime(\"2021-01-14T18:39:04\")\n",
"\n",
"COCO_P = UTCDateTime(\"2021-01-14T18:33:28\")\n",
"COCO_S = UTCDateTime(\"2021-01-14T18:37:43\")"
]
},
{
"cell_type": "markdown",
"id": "6be7700a",
"metadata": {},
"source": [
"### Selisih waktu tiba S dan P\n",
"Selisih S-P kemudian dapat kita cari dengan cara yang sederhana yaitu:"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "fa048e6e",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"255.0 sekon\n"
]
}
],
"source": [
"TATO_SP = TATO_S-TATO_P\n",
"PMG_SP = PMG_S-PMG_P\n",
"COCO_SP = COCO_S-COCO_P\n",
"\n",
"print(COCO_SP,\"sekon\")"
]
},
{
"cell_type": "markdown",
"id": "2ef8f2ea",
"metadata": {},
"source": [
"### Mencari jarak dengan *rule of thumb*\n",
"\n",
"Dari nilai SP sebenarnya kita langsung dapat mencari radius lingkaran/jarak dengan rumus Havskov et al (2002): \n",
"\n",
"$$jarak=(SP-2)\\times10,$$\n",
"\n",
"dengan satuan kilometer, apabila akan dijadikan derajat maka kita bagi dengan angka 111,26 km:"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "feba269b",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[25.43591587273054, 25.795434118281502, 22.73952903109833]"
]
},
"execution_count": 7,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"jarak=[(SP-2)*10/111.26 for SP in [TATO_SP,PMG_SP,COCO_SP]]\n",
"jarak"
]
},
{
"cell_type": "markdown",
"id": "05c962fa",
"metadata": {},
"source": [
"Kemudian dapat kita plot:"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "e3e5a60c",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAATsAAAEICAYAAAAgMlPEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABOJElEQVR4nO2deVhV1f6H38WkaOaAoCAOaQ45IChOqWkO6XVA05yzrExTy6zMTLMsLavr7ddw85bmzSm10kpwHsrKBq/gjPMYCuKAEwIKnO/vj30gQEDQc9jnHNb7PPuBs/bea332Pvt89pqXEhE0Go3G1XEzW4BGo9EUBdrsNBpNsUCbnUajKRZos9NoNMUCbXYajaZYoM1Oo9EUC+xqdkqpNUqpx+2ZhquilKqhlBKllIeN441WSrW3QTztlVKn7lyRxtFRSk1SSn1hto475ZZmp5Q6oZRKVkolZtn+XZDIReQfIjL/zmXmq0+UUvfaIc6zWY1GKeVpDStQx0Sl1DCl1BZb6iooSqmpSqlFue0TkQYisrmIJRUpyuBZpdRupVSSUuqMUmqzUmqg2dpuhVJqnvX565Uj/P+s4cOKWpOIvCMiw2/n3MJcj1LKVym1WCl1WSl1USn1VR5x+imlliilYq3H/qaUanErLQXN2fUUkbuybM8W8DyH5ha5povAP7J8/oc1TGNDbJ1ztfIxMA54CfABqgCvAV3tkJY9OAQ8lvHBeo/6A0dNU3RnFPR6vgPOANUAP2BmHvHdBWwDmgIVgPnAKqXUXfmqEJF8N+AE0CmPfcOA34B/A5eBA0DHLPs3A8Ot/98L/Gw97jzwdZbj6gEbgATgINA/y755wKfAKuAqsBWoZd33CyDANSARGGAN7wHsBC4BvwNBOa7nFWA3cB3wyOW6BOPH8W2WsGXAZOOWZYaVBeYCccBpYDrgDtwHpADpVl2Xshy/ADgHnLSm4Wbd547x5Z4HjgFjrDo88ksrj+9lKrDoVt+n9bhvrJquAtFAaJZjmwA7rPu+Bb4Gplv3tQdOZTl2LLAPCASaA39Y73+c9fnwynF/xwCHgePWsI+AGOAKEAW0zXE9eerMcX11rPc91/23+u5yPNf/Z72GY8D91vAY4CzweI5n9DOMZ/gqxnNePcv+PK8tF13zrM9BPFA+y/O8BtgCDLOG1QJ+BC5Yn5mvgHI5vueXMZ7za9ZrrWSN5yqwMUv82b7LPJ6TRQU59g6u5yFrPLk+0wXwqStA0/yOsUWdXQsMh64IvAF8p5SqkMtx04D1QHmMH8QnAEqp0hgPyWIMNx8IzFJK1c9y7kDgTeu5R4C3AUTkAev+xtYc59dKqRDgv8BIjLf650C4UqpElvgGAd0xHo60PK7rB+ABpVQ5pVR5oC2wIscx84A0DCMPwfjChovIfuAZ4A+rrnLW4z/B+JHVBNphvO2esO57GuMhCAFCgUcKklYe2gtDGLAUKAeEYxgTSikv4HtruhWAJcDDuUWglHodwwjaicgpDLN5AeOZaAV0BEbnOK03xrOT8T1vA4KtaS0GvlVKlbyVzlzoAMSISGR+F82t72cLDKPwsepZCjSzHv8o8O8cOYkhGM94RYwXbdYi2K2uLScpGM9aRrH7MQyjz4oCZgABGC/XqhimlJW+QGeMF0BPDIOZBPhilOrG5qPBlhTkelpiZHTmK6UuKKW2KaXaFSRypVQw4IXhDXlTAMc8gTV3kmV7OssbMBZQWY7/HzBUbs7ZLQBmA4E54h8A/Joj7HPgjSxvhi+y7OsGHMiRS7g3y+f/ANNyxHcQ44eYcT1P3uKaBeOh/gLDNJ8B5ljDxHpMJYycoXeW8wYBP2W5N1uy7HMHbgD1s4SNBDZb//8ReCbLvoesOjxulVYu+qdS8Jzdxiz76gPJ1v8fwMjxZP1ut5A9Z3ca+MAaXjaf+zkO+D7H/e1wi+/gIsZLLF+duZz3GvBnjrBTGM9tClC9gN/d4Sz7Glk1V8oSdgEIzvKMLs2y7y4Mw696q2vLZd88jFxmG4zccTmMXJE3WXJCuZzXG9iR43sekuXzcuA/WT4/B/yQ5bu0Z87ulteD4Q0CPAV4YhjjJaDiLZ6Tu4E9wKv5HSciFLS+pLeIbMxj3+mMX7+Vkxhvm5xMwHjz/U8pdRH4l4j8F+Pha6GUupTlWA9gYZbPZ7L8n4TxMOVFdeBxpdRzWcK8cmiKyef8rCzAeHsqjKJvznQ8gTilVEaYWz5xV7QefzJL2EmM+iSs+mJy7LvdtApDzntb0lqnEsDN323O9MoBIzCqDy5nBCql6mCYYChQCuP7jMpxbra4lFLjMR70AIyH/m6Me5avTrk5Z34B8M8aICKB1mtKxfguC3I/47P8n2yNJ2dY1ucw81wRSVRKJVivJaYA13YTIrJFKeWLUXWyUkSSs2hFKVUJo3jcFihj1X8xRzQ59ean367c6nqsek6IyFzr56VKqclAa24uUQGglPIGIjBebjNupcEWxdgqKrvqahi5vWyIyBkReVpEAjByNLOsragxwM8iUi7LdpeIjLpNPTHA2zniKyUiS7LKKWBcv2L8cCphvIVypnMd482Tkc7dItIgjzTOY/zYqmcJq4aROwKj7qhqjn0FTcsexHHzd1s1xzEXMYreXyqlWmcJ/w9G/W1tEbkbo+ikcpybeX+UUm0xXob9Mep1ymHU7eY8pyD8CAQqpULzOcYe9zPz3liLtxWA2Du8tkUYjSw5i3wA72Dcw0bWe/xoAePMjWsYL6UM/e4YRd07PTYn+V3Pbm7+zeT5O7VWS/2AkWsfWZDEbWF2fsBYa9eMfhj1B6tzEddPKRVo/XgR40IswEqgjlJqqDUOT6VUM6XUfQVMPx6jDiyDOcAzSqkW1i4IpZVS3ZVSZQp7YdZcTU8gLEcOBxGJw6iD/JdS6m6llJtSqlaWeoZ4jB+dl/X4dIxK9reVUmWUUtWBFzEeAKz7xiqlAq11hBMLkVZuuCmlSmbZSuRzbG78gVEUe1Yp5WHtOtA8l3u0GaO+6julVMb+MhgVxolKqXrArV5cZTDqz84BHtY6wLsLqTdDz0GMapClSqnOSilv6w/y/izH3M79vBXdlFJtrN/3NIzcRswdXtvHGHVuv+SyrwxG9dJlpVQVjMaI2+UQRk65u1LKE6MqIK/npTDH5iS/6/keKK+Uelwp5a6UegSjbv83yOxOtdn6vydGg2EyRkORpSCJF9TsIlT2fnbfZ9m3FaiNkXN5G3hERC7kEkczYKtSKhGjgvl5ETkmIlcx6qcGYuQIzwDvUfAbOBWjUvOSUqq/GBXTT2NUYF/EqLQcVsC4bkJEokUkOo/dj2EUkfdZ01rG30WoHzFaDc8opc5bw57DeDMew8gpLsZoTAHDpNcBu4DtGM3wBU0rNwZhPAwZW6G6LYjIDaAPRvHrEkbOYSVGjijnsRuAJzGekybAeGAwRqvfHIxW3PxYB6zF+CGdxKhbu5Mi+hiMH9YHGC38pzAMaADwl/WYwt7PW7EYo4EuAaNLxKPW8Nu+NhFJEJFNOV+0Vt7EaC2/jNFTIefzUmCsVRCjMeqoT2M8o7l2GC/Msbmcm+f1iEgCRiPUeIxrmgj0EpGM305VrMaH8eLqgeEbl7L4Utv80le538eCoYwOgcNFpM1tR6JxGpRSW4HPRORLs7U4EkqpeRiV9q+ZrcUeKKXewmhYfNJEDTsxurXllpEqEHpsrCZPlFLtlFKVrcXYx4EgjFyKpphgrbOtDxw3U4eIBN+J0YEdzE4Z/dKWKaUOKKX2K6VaKaUqKKU2KKUOW/+Wt3W6GrtQF6NYfQmjYvkRa32XpviwHaPubI7ZQu6UOyrG5hqhUvMx+s19Ya2sLYXRGpcgIu8qpSZitErl7Mqh0Wg0dsOmZqeUKovRe7xm1kpIpdRBoL2IxCml/DE60ta1WcIajUZzC2w9CPsejCb2L5VSjTE6kj6P0fM8o/hzBqPf2k0opUZgdFKldOnSTevVq2djeRqNpiiIioo6LyIF7X9XJNg6ZxcK/Am0FpGtSqmPMPpbPSd/jw9FKXVRRPKttwsNDZXIyFsNb9RoNI6IUipKRPLr2F3k2LqB4hRGE/xW6+dlGH2B4q3FV6x/z9o4XY1Go8kXm5qdiJzBGAuYUR/XEaPTZjjwuDXscfIY66bRaDT2wh4TJz4HfGVtiT2GMYWRG/CNUuopjF7k/e2Qrkaj0eSJzc1ORHZizHaRk462TsuVSE1N5dSpU6SkpJgtRaMpMCVLliQwMBBPT0+zpdwSe+TsNLfBqVOnKFOmDDVq1CD7RCMajWMiIly4cIFTp05xzz33mC3nlujhYg5CSkoKPj4+2ug0ToNSCh8fH6cpjWizcyC00WmcDWd6ZrXZaTSaYoE2O00m7u7uBAcHZ27vvvtuvsd/9tlnLFiQ26Szt88777xj0/huxbx583j2WZdYGRQRoUOHDly5cgWAu+4yZl0/ceIE3t7eBAcHU79+fR577DFSU1MB2Lx5Mz169AAgPj6eHj160LhxY+rXr0+3bt0AOHfuHF27OssqlHmjGyg0mXh7e7Nz584CH//MM8/YXMM777zDpEmTCnVOeno67u7uNtdiD9LS0vDwsM/PbvXq1TRu3Ji77755IuRatWqxc+dO0tPT6dy5M9988w1DhgzJdszrr79O586def755wHYvXs3AL6+vvj7+/Pbb7/RunXrm+J2FnTOzln56iuoUQPc3Iy/X+W6eLpNqFGjBhMmTKBRo0Y0b96cI0eMFeumTp3KzJnGOsYff/wx9evXJygoiIEDjRXzrl27xpNPPknz5s0JCQlhxQqjL/m8efPo06cPXbt2pXbt2kyYMAGAiRMnkpycTHBwcOYPcdGiRTRv3pzg4GBGjhxJeno6YORaXnrpJRo3bswff/yRTW9htADExMTQvn17ateuzZtvvpkZnl/akydPpnHjxrRs2ZL4eGMdm4iICFq0aEFISAidOnXKDJ86dSpDhw6ldevWDB06lKlTp/Lkk0/Svn17atasyccffwwYZvPhhx9mpj958mQ++ugjEhMT6dixI02aNKFRo0bZtGflq6++olevXvl+l+7u7jRv3pzTp0/ftC8uLo7AwMDMz0FBQZn/9+7dm6/s+IwVCbdafsysrWnTplKc2LdvX8EPXrRIpFQpEfh7K1XKCL8D3NzcpHHjxpnb0qVLRUSkevXqMn36dBERmT9/vnTv3l1ERN544w355z//KSIi/v7+kpKSIiIiFy9eFBGRV199VRYuXJgZVrt2bUlMTJQvv/xS7rnnHrl06ZIkJydLtWrV5K+//hIRkdKlS2fq2bdvn/To0UNu3LghIiKjRo2S+fPni4gIIF9//XWu11FYLZUrV5bz589LUlKSNGjQQLZt23bLtMPDw0VE5OWXX5Zp06aJiEhCQoJYLBYREZkzZ468+OKLmfepSZMmkpSUlPm5VatWkpKSIufOnZMKFSrIjRs35Pjx4xISEiIiIunp6VKzZk05f/68pKamyuXLl0VE5Ny5c1KrVq3MdLJSrVo1uXLlSubnjHt5/PhxadCggYiIJCcnS/v27WXXrl0iIvLTTz9lfp9r166VsmXLSvv27WX69Oly+vTpzLhOnTolDRs2zPV+5/bsApHiAD6SddPFWGdk8mRISsoelpRkhOcomhSG/IqxgwYNyvz7wgsv3LQ/KCiIIUOG0Lt3b3r37g3A+vXrCQ8Pz8z9paSk8NdfxhIQHTt2pGzZsgDUr1+fkydPUrVq9sXLNm3aRFRUFM2aNQMgOTkZPz8/wMih9O3bN1ethdXSuXNnfHx8AOjTpw9btmzBw8Mjz7S9vLwy67maNm3Khg0bAKOv5IABA4iLi+PGjRvZ+p6FhYXh7e2d+bl79+6UKFGCEiVK4OfnR3x8PDVq1MDHx4cdO3YQHx9PSEgIPj4+pKamMmnSJH755Rfc3Nw4ffo08fHxVK5cOdt1JyQkUKZM7utKHT16lODgYI4fP0737t2z5doy6NKlC8eOHWPt2rWsWbOGkJAQ9u7di6+vL35+fsTG3rRooFOhzc4Z+euvwoXbAJV9zdKb9q9atYpffvmFiIgI3n77bfbs2YOIsHz5curWzT514datWylR4u/1lNzd3UlLy7n8q1HqePzxx5kx4+YlQUuWLJlnPV1hteS8HqVUvml7enpmnpNV+3PPPceLL75IWFgYmzdvZurUqZnnlC5dOlsceV3/8OHDmTdvHmfOnOHJJ40lH7766ivOnTtHVFQUnp6e1KhRI9e+bR4eHlgsFtzcbq6dyqizO3/+PK1btyY8PJywsLCbjqtQoQKDBw9m8ODB9OjRg19++YW+ffuSkpKSzaydEV1n54xUq1a4cBvw9ddfZ/5t1apVtn0Wi4WYmBgefPBB3nvvPS5fvkxiYiJdunThk08+QazTiO3YseOW6Xh6ema2FHbs2JFly5Zx9qwxSU5CQgInT57M7/Tb0rJhwwYSEhJITk7mhx9+oHXr1reV9uXLl6lSxVjzfP78+be81tx4+OGHWbt2Ldu2baNLly6Z8fr5+eHp6clPP/2Up466dety7NixfOOvWLEi7777bq4m/uOPP5JkLTFcvXqVo0ePUs36TB06dIiGDRve1jU5Cjpn54y8/TaMGJG9KFuqlBF+B2Q0DmTQtWvXzO4nFy9eJCgoiBIlSrBkyZJs56Wnp/Poo49y+fJlRISxY8dSrlw5pkyZwrhx4wgKCsJisXDPPfewcuXKfDWMGDGCoKAgmjRpwldffcX06dN56KGHsFgseHp68umnn1K9evU8z78dLc2bN6dv376cOnWKRx99lNBQY2h3YdOeOnUq/fr1o3z58nTo0IHjxwu/Ro2XlxcPPvgg5cqVy8y5DhkyhJ49e9KoUSNCQ0PJa1Lb7t27s3nzZu6999580+jduzdTp07l119/zRYeFRXFs88+m5lDHD58eGYx/qeffqJ79+6Fvh6HwuxKw7w23UBxCxYtEqleXUQp4+8dNk7kR/Xq1eXcuXN2i1/zN+np6dK4cWM5dOhQoc+NjY2VTp062UGVSNu2bSUhISHXfc7SQKGLsc7KkCFw4gRYLMbfO2iY0DgG+/bt495776Vjx47Url270Of7+/vz9NNPZ3YqthXnzp3jxRdfpHx5514UUBdjNbfkxIkT+e63WCykpqaSmpqa7U2qlEIphZubG56entkq9jU3U79+/VvWud2K/v1tP1Wkr69vZqu2M6PNzoHIMAhHRERISUnh2rVrJCUlkZKSwo0bN0hNTc2s0/L09MxsCcxo0RQR0tPTSU1NJT09HQ8PDzw9PfHy8sLb25vSpUtTqlQpbYROiohtl2K1J9rsHISSJUty4cIFh5nmSURITEzk0qVLJCYmkpycjJeXF6VKlaJUqVKULVs207Tc3d0LpNlisZCWlkZqaio3btwgKSmJs2fPZrYAli5dmjJlylCuXDlKlixp70vU3CEixnx2zvJdabNzEAIDAzl16hTnzp0zTYPFYiE5OTlz8/DwwNvbm5IlS1KiRAnc3Ny4fv06169f5+LFizZLt0SJEqSlpXHlyhXi4+NJTk7Gzc0Nb29vvL29KVGihEO8ADQ3kzFTsTOgzc5B8PT0NG2216ioKGbNmsXy5cu5//77CQsLo0ePHqY9xBaLhe3btxMeHk54eDgJCQmMHDmS4cOHU6lSrksOazS3RLfGFlOSk5OZP38+LVq0oG/fvtSuXZtDhw6xevVqnnnmGVPf1m5uboSGhvLWW2+xc+dOVqxYwcmTJ6lXrx6DBg3i119/daq6Io1joM2umHH9+nU++eQTatasyddff82UKVM4evQoEydOzBz76WiEhIQwe/Zsjh8/TqtWrXj66adp2bIlP/30k9nSNE6ENrtiQnp6OosWLaJevXqsWbOGNWvWsHr1anr06OE0c8GVK1eOsWPHsm/fPsaNG8fw4cPp2rVrgYahaTTa7IoBmzZtIiQkhE8//ZR58+axevXqbMPCnA03NzcGDRrE/v37CQsLo1u3bgwaNIiYmBizpWkcGG12LsyVK1cYOXIkTz75JG+++Sa///477dq1M1uWzfDy8mL06NEcPnyYunXr0qRJE+bOnavr8zS5os3ORdmwYQNBQUGICHv27OHhhx922e4bd911F1OnTmXTpk18+umndOvWTefyNDehzc7FSEpKyuymMXv2bGbPnp3rmgSuSFBQEFu3buX++++nSZMmtz3NksY10WbnQsTExNC2bVuuXbvGnj17eOihh8yWVOR4enoyZcoUNm3axIwZM3juuecy58fTFG+02bkIv//+Oy1atGDQoEEsXLiw2OTm8iIoKIg///yTI0eO0LVrVy5cuGC2JI3JaLNzAf773//Su3dvvvjiC8aPH++ydXOFpVy5cqxcuZImTZrQvHlzoqOjzZakMRE9XMyJERHeeOMNlixZwi+//JLnDLbFGXd3d/75z38SFBTEgw8+yPLly2nbtq3ZsjQmoM3OSRERXnnlFdatW8dvv/3msKMfHIWhQ4fi7+9P3759WbJkCR07djRbkqaI0cVYJ0REePnll9m0aRM//vijNroC0qlTJ5YtW8agQYP48ccfzZajKWK02Tkhr7/+Ohs3bmTDhg2Z651qCsYDDzzAt99+y4ABA9iyZYvZcjRFiDY7J2PWrFksW7aMDRs2UKFCBbPlOCXt2rVj8eLF9O3bl4MHD5otR1NEaLNzIn766SfefPNNVq5cia+vr9lynJrOnTvzzjvvEBYWxqVLl8yWoykC7GJ2Sil3pdQOpdRK6+d7lFJblVJHlFJfK6W87JGuK3Ps2DEGDRrE4sWLqVWrltlyXIKnnnqKrl27MmjQINLT082Wo7Ez9srZPQ/sz/L5PeD/RORe4CLwlJ3SdUmuXr1Kr169eO2113Qroo3517/+RVpaGhMnTjRbisbO2NzslFKBQHfgC+tnBXQAllkPmQ/0tnW6rkzGZJVjxowxW4rL4eHhwddff83333/PN998Y7YcjR2xRz+7D4EJQBnrZx/gkoikWT+fAqrkdqJSagQwAqBatWp2kOZ8fPvtt+zcuZMdO3bokRF2okKFCixevJiwsDDatWun17lwUWyas1NK9QDOikjU7ZwvIrNFJFREQnUFPJw9e5bnnnuOefPm4e3tbbYcl6Z58+YMGzaM0aNH6/nwXBRbF2NbA2FKqRPAUozi60dAOaVURi4yEDht43RdkmeffZahQ4fSsmVLs6UUC6ZOncr+/ft1cdZFsanZicirIhIoIjWAgcCPIjIE+Al4xHrY48AKW6brinz77bfs3r2bt956y2wpxYaSJUsyb948nn/+eeLj482Wo7ExRdXP7hXgRaXUEYw6vLlFlK5TkpKSwksvvcScOXN08bWIad68OYMGDWLq1KlmS9HYGOWo9ROhoaESGRlptgxT+OCDD/j5559ZsUJngM3gwoUL1K1blz/++IPatWubLccpUUpFiUio2TqyokdQOBiXL1/m3Xff5Z133jFbSrHFx8eHF154gSlTppgtRWNDtNk5GDNnzqRbt240aNDAbCnFmnHjxvHLL7+wfft2s6VobISez86BOHv2LLNmzdI/MAegdOnSvPbaa0yaNIm1a9eaLUdjA3TOzoGYPXs2ffr0oXr16mZL0QDDhw9n165d7N+//9YHaxwebXYOQlpaGp9//rkeEuZAeHl5MXz4cP7zn/+YLUVjA7TZOQirVq2iatWqBAcHmy1Fk4URI0awaNEiEhMTzZaiuUO02TkIs2bNYvTo0WbL0OSgatWqmZN9apwbbXYOwJEjR9ixYwePPPLIrQ/WFDmjR49m1qxZZsvQ3CHa7ByA5cuX88gjj1CyZEmzpWhyoWPHjsTHx3P48GGzpWjuAG12DkB4eDi9evUyW4YmD9zc3OjZsycRERFmS9HcAdrsTObs2bNER0fTvn17s6Vo8iEsLIzw8HCzZWjuAG12JrNq1So6d+5MiRIlzJaiyYeOHTuyfft2EhISzJaiuU202ZlMREQEYWFhZsvQ3AJvb286dOjA6tWrzZaiuU202ZmIiPDzzz/rRXSchE6dOvHLL7+YLUNzm2izM5GTJ09SsmRJAgICzJaiKQBNmzaluE475gposzORyMhImjZtarYMTQFp3LgxBw4c4Pr162ZL0dwG2uxMJCoqSpudE1GqVClq1arFnj17zJaiuQ202ZmINjvnIzQ0lKio21o8T2My2uxMZO/evTRu3NhsGZpC0LhxY52zc1K02ZlEWloa586dw9/f32wpduHChQsEBwcTHBxM5cqVqVKlSubns2fP4unpyWeffQbAmDFjCA4Opn79+nh7e2cet2zZMkSE6dOnU7t2berUqcODDz5IdHS0addVpUoVYmNjTUtfc/vomYpNIj4+Hh8fHzw8XPMr8PHxYefOnYCxHutdd93F+PHjAfjPf/5Dy5YtWbJkCc888wyffvopACdOnKBHjx6Z5wH8+9//5vfff2fXrl2UKlWK9evXExYWRnR0tCljiQMCAoiLiyvydDV3js7ZmURcXFyx7XKyZMkS/vWvf3H69GlOnTqV77Hvvfce//73vylVqhQADz30EPfffz9fffVVUUi9CX9/f212Too2O5OIi4tz2SJsfsTExBAXF0fz5s3p378/X3/9dZ7HXrlyhWvXrlGzZs1s4aGhoaYVZTPMzlGXINXkjTY7k4iLi6Ny5cpmyyhyvv76a/r37w/AwIEDWbJkicmKCoe3tzelSpXSY2SdEG12JpGUlMRdd91ltgybEhgISt28ffDBi5nHLFmyhHnz5lGjRg3CwsLYvXt3nvPE3X333ZQuXZpjx45lC4+KijJ1qcnSpUuTlJRkWvqa20ObnUmkpaXh6elptgybEhYGXl7Zw7y8oG7dAwAcOnSIxMRETp8+zYkTJzhx4gSvvvpqvrm7l19+mbFjx5KcnAzAxo0b2bJlC4MHD7bbddwKT09P0tLSTEtfc3toszOJtLQ03N3dzZZhU6ZMAbccT5S7O7RrZwyeX7JkCQ8//HC2/X379s3X7J577jmaNWtGo0aNqFu3LtOmTWPFihV4e3vbXH9BcXd312bnhLhmvwcnwMPDg/T0dLNl2BR/f3jiCZg7F27cMHJ1TzwBM2eOz/OcoKCgzHVZa9Sowd69e7PtV0rxxhtv8MYbb9hVe2FIT0932S5DrozO2ZmEh4eHS+YOsubu3N2Nz65GWlqaNjsnRJudSXh7e3Pt2jWzZdicjNydm5vx1xUbnK9du2ZqMVpze2izM4nKlSu7bOfUKVOgRg3XzNWlpKRw7do1fHx8zJaiKSQ6L24SrjzsyN8fjh41W4V9OHPmDJUqVUIpZbYUTSHROTuT0MOOnJPY2NhiOfLFFdBmZxKVKlXi3LlzLtci6+oU5zHNzo5NzU4pVVUp9ZNSap9SKlop9bw1vIJSaoNS6rD1b3lbpuuMeHp6UqFCBc6cOWO2FE0hOH36tM7ZOSm2ztmlAS+JSH2gJTBGKVUfmAhsEpHawCbr52JPw4YN2b17t9kyNIVg9+7dNGzY0GwZmtvApmYnInEist36/1VgP1AF6AXMtx42H+hty3SdlaZNm+opvp0MvUiS82K3OjulVA0gBNgKVBKRjNr4M0ClPM4ZoZSKVEpFnjt3zl7SHAZtds5FSkoKBw8eJCgoyGwpmtvALmanlLoLWA6ME5ErWfeJMRFYrpOBichsEQkVkVBfX197SHMo9OItzsXu3bupU6eO7lDspNjc7JRSnhhG95WIfGcNjldK+Vv3+wNnbZ2uM1KzZk2uXr1KfHy82VI0BSAqKorQ0FCzZWhuE1u3xipgLrBfRD7IsisceNz6/+PAClum66wopWjTpg0//vij2VI0BWDTpk20bdvWbBma28TWObvWwFCgg1Jqp3XrBrwLdFZKHQY6WT9rgLCwMMLDw82WobkFKSkpbNiwge7du5stRXOb2HS4mIhsAfIaR9PRlmm5Cj169GDChAmkpqa63GSersTmzZtp1KgRxaEu2VXRIyhMxt/fnzp16vDrr7+aLUWTDxEREYSFhZktQ3MH6IkAHICMomyHDh3MlqLJQVpaGvHx8Sxbtozp06cTERFBWloaIoKHhwceHh6UK1cOf39//P39TVnLVlMwlKMuCRcaGiqRkZFmyygSDhw4wIMPPsjJkyfxyrmIg6ZIOH/+PFFRUURFRbF9+3aOHTtGXFwc58+fp0yZMiQnJ/PAAw/g5eWVOXFneno6qampXLp0idjYWM6cOUPp0qXx9/enatWqBAcHExoaStOmTalRo0axmilFKRUlIg7VdK1zdg5AvXr1qFevHj/88EPmMoMa+3LmzBlWrVrF2rVr2bZtG5cuXSIkJISmTZvSt29fateujb+/P5UqVWLIkCG0bduWZ599Nt84RYQLFy4QFxfHiRMn2LFjBwsXLmTcuHEkJyfTpEkTOnbsSFhYGPfdd1+xMj9HQOfsHIRvvvmGWbNmsXnzZrOluCQiQnR0NOHh4YSHh3Pw4EG6dOlC9+7dadmyJbVq1cIt52pBGLOc1K9fnxMnTlC2bNnbTj8+Pp5t27axbt06wsPD8fT0pGfPnoSFhdGmTRuXa5xyxJydNjsH4caNG1SvXp2NGzeauiaqq3Ht2jUWL17MrFmzSEhIoHfv3oSFhdG2bdsCVRlMmzaN06dP89lnn9lMk4iwe/fuTOM9ffo0Tz/9NE8//TSBgYE2S8dMHNHsEBGH3Jo2bSrFjSlTpsjo0aPNluES7N+/X8aOHSsVKlSQXr16ybp16yQ9Pb1QcaSmpkpgYKDs3LnTTioN9u7dK2PGjJHy5ctLnz59ZOPGjWKxWOyapr0BIsUBfCTrZrqAvLbiaHanTp2SChUqyOnTp82W4rRER0dLr169pFKlSjJ58mQ5efLkbcc1d+5cad++vQ3V5c+VK1fkP//5jzRq1EgaNWokERERTmt62uy02d2S8ePHy8iRI82W4XScPHlShg0bJr6+vjJz5kxJTk6+o/iSkpKkatWq8vvvv9tIYcGxWCzyww8/SP369aVNmzby22+/FbmGO8URzU53KnYwJk6cyLJlyzh8+LDZUpyCy5cv89JLLxESEkKVKlU4fPgwL7300h33d5s1axZNmjShVatWNlJacJRS9OrVi927dzN8+HAGDRpEr169OHjwYJFrcSW02TkYPj4+vPjii7z22mtmS3F41q1bR6NGjbhy5QrR0dFMnz79jlpMM7h8+TLvvfceb7/9tg1U3j7u7u48/vjjHDx4kHbt2tGmTRtmzpyp1y25XczOWua1FddirIhIYmKi+Pv7y//+9z+zpTgkly5dkqeeekqqV68uGzZssHn8r776qgwbNszm8d4pR48elfbt20urVq3kwIEDZsvJF3QxVlMQSpcuzbvvvsvw4cO5ceOG2XIcig0bNhAUFISHhwd79uyhU6dONo1/165dzJkzh+nTp9s0XltQs2ZNNm3axJAhQ2jTpg3/+te/sFgsZstyHsx227y24pyzEzEqqbt37y6vv/662VIcAovFIu+++64EBATI+vXr7ZLGjRs3JDg4WP773//aJX5bcvToUWnVqpX06dNHrl69aracm8ABc3a6U7GdERHi4uI4ffo0sbGxxMXFERcXx4ULF0hNTSUtLQ13d3c8PDwyx1UGBATg7++Pm5sbffr0Yf369YSEhJh9KaaRnJzM8OHDOXjwID/88IPdOt6+9dZb/Pnnn6xatcophnJdv36d0aNHExkZyYoVK6hRo4bZkjJxxE7FemysjYmJiWHr1q2Zg8qjoqLw8PAgMDAw08QCAgKoW7cunp6eeHh4YLFYSEtL4+rVq/z1119s3bo10xivXLnC/fffz+DBg2nWrBlNmzYlJCQkczC6q3Pq1Ckefvhhateuza+//mq39R927drFJ598wo4dO5zC6ABKlCjBF198wSeffEKrVq1YunQp7dq1M1uW42J21jKvzVmKsenp6bJ161aZPHmyNGrUSCpWrCg9e/aUqVOnSkREhMTGxt5R/NeuXZM2bdrIgw8+KE888YQ0bNhQfHx85NFHH5VvvvlGLl++bKMrcTz27NkjVapUkRkzZti1c+3ly5elfv368uWXX9otDXuzfv168fPzk4ULF5otRUQcsxhruoC8Nkc3u4zhSJUrV5b69evLxIkT5bfffpO0tDSbp5WQkCC1a9fOrEv666+/ZNasWdK1a1cpU6aMdO3aVb777jtJTU21edpmsX37dqlUqZIsXrzYrumkp6dLz549ZeTIkU47WiGD6OhoCQgIkC+++MJsKdrsCrM5otnduHFDli1bJh06dMgcjnTo0KEiSXvfvn3i6+t7U4/+K1euyMKFC+X++++XwMBAeeutt+44N2k2UVFR4ufnJ8uWLbN7WpMmTZIHHnhArl+/bve0ioJDhw5J1apV5fPPPzdVhzY7JzW71NRUmT17tgQGBkqbNm1k8eLFkpKSUuQ6Vq1aJQEBARITE5Pr/h07dsiIESOkXLlyMnz4cDl16lQRK7xz9uzZI5UqVZLly5fbPa2lS5dKjRo15OzZs3ZPqyg5fPiwVKlSRebPn2+aBm12TmZ2FotFvv32W6lTp460b99e/vzzT7MlyXvvvSfBwcFy8eLFPI+5cOGCvPLKK1KhQgWZMGGCXLhwoegE3gEnT56UgIAAuxddRUS2bNkiFStWtPuMJmaxf/9+8ff3l/DwcFPS12bnRGb3v//9T5o1ayYhISGydu1ah6nPsVgsMnbsWGnZsqVcuXIl32NPnTolI0aMkIoVK8oHH3xgl/pEW5GYmCjBwcEyc+ZMu6e1detW8fX1lXXr1tk9LTP5888/pWLFirJ3794iT1ubnROYXUpKirz66qvi5+cnCxYsKPQcaEWBxWKRESNGSJs2bQrUGnvgwAFp166dww4zslgs0q9fP3nsscfs/lLZtm2b+Pn5mZbjKWrmz58vNWvWlPPnzxdputrsHNzstm3bJg0aNJDevXtLXFxckadfGNLT02XkyJHSokWLfIu0WY//5JNPxMfHR2bOnOlQubxp06ZJixYt7nhapluxZcsW8fX1LTZGl8H48eOlQ4cOcuPGjSJLU5udg5qdxWKRmTNnip+fnyxevNhhiqy3wmKxyPPPPy+NGjWS48ePF+ico0ePSrt27aR9+/Zy7tw5+wosAOHh4RIYGGj3FuTvv/9eKlasKGvXrrVrOo5IWlqadO3aVcaNG1dkaWqzc0CzS05Olscee0yCg4PvaFZbs7BYLPJ///d/UrlyZfn5558LdE5aWpq8+uqrcs8998iuXbvsrDBvzp07J/7+/vLrr7/aLQ2LxSLTpk2TwMDAYj2LTEJCggQGBsqmTZuKJD1tdg5mdrGxsdKiRQvp16+fJCYm2j09e7Ju3Trx8/MrVP+qxYsXS8WKFeW7776zo7K8GThwoLzwwgt2iz8xMVH69esnLVq0cPq+h7Zg1apVUqNGjSKZOECbnQOZ3b59+6Rq1aoybdo0pym23opDhw5JvXr1ZNSoUZKUlFSgc7Zt2yaBgYFF0gqaleXLl0vt2rXl2rVrdon/8OHDEhwcLEOHDrV7XaAzMWzYMBk1apTd09Fm5yBmt2vXLqlcubIsWLDAbmmYxaVLl2TAgAFSp06dAq9dEBMTI7Vr15Zp06bZWZ3BuXPnpHLlyrJlyxabx52eni4fffSR+Pj4yEcffeQyLzJbcfHixSIpzmqzcwCzy+ihv3TpUrvE7ygsW7ZMKleuLC+99FKBcnmxsbFy3333FYnhDR8+XJ5//nmbx3v48GFp27attG7dusiG8TkjK1eulHvvvdeurbPa7Ew2uwMHDoi/v3+R9NB3BM6dO5eZyyvImzw2NlZq165t1yJtxhjfgnSXKSgpKSnyr3/9S3x8fBy+87Sj0LlzZ5k1a5bd4tdmZ6LZXbhwQe69916ZM2eOTeN1Br777jupWbOmPPTQQxIVFZXvsTExMVK1alW7jU19+OGH5f3337dJXGlpabJw4UKpUaOG/OMf/5CDBw/aJN7iQFRUlPj7+9utYU6bnUlml5qaKp06dZIXX3zRZnE6G9evX5dPP/1UKleuLAMGDJDDhw/neWxkZKRdxo3+8ccfEhgYWODGk7ywWCyycuVKadSokbRq1arAXW402RkwYIBMnz7dLnFrszPJ7J5//nnp0qWLS833drtcvXpVpk+fLj4+PjJgwAD5+eefc63EX7JkiU1nBLFYLNKuXbs7ylknJSXJl19+KaGhoVK/fn354YcfdAPEHXDo0CHx8fGxy1AybXYmmN28efOkTp06Nq0jcgUuXbokH3/8sdSrV08aNGggn3766U3jbCdNmiTt2rWzSR3Yb7/9JrVq1bqtF87hw4flpZdekooVK0q3bt1k5cqVul7ORgwfPlymTp1q83iLvdkBXYGDwBFgYn7H2sLsTpw4IRUrVpQ9e/bccVyuisVikR9//FEeeeQRKVeunAwYMEAWLVokFy5ckPT0dHnwwQdtUsc2ZMgQ+eCDDwp8/MGDB+Wf//yntG3bVnx9fWXChAly9OjRO9ahyc6uXbukSpUqNm+ZLdZmB7gDR4GagBewC6if1/F3anYWi0U6deok77zzzh3FU5w4c+aMfPHFFxIWFiZlypSR9u3by+TJk6Vs2bJ3NKwsPj5eypUrJwkJCXkec/HiRdm4caOMHz9e6tatKwEBATJy5EhZtWqV7hRsZ9q0aWPzBilHNLuiXKKqOXBERI4BKKWWAr2AffZIbM6cOVy6dImXX37ZHtG7JJUqVeKpp57iqaeeIikpiU2bNhEREUGpUqUICQkhJCSE0NBQmjZtSsOGDQkICKBy5cqUKFEi33j/+9//0qdPH8qXL8/Vq1eJjY0lJiaGHTt2EBkZSVRUFPHx8TRu3JiOHTuyaNEimjRpgpubXsO9KBg9ejSzZs2iT58+ZkuxK0W2bqxS6hGgq4gMt34eCrQQkWezHDMCGAFQrVq1pidPnryttGJiYmjSpAmbN2+mQYMGdy6+mGOxWGjfvj1BQUHUqVOHqKgo9u/fT1xcHPHx8dx9992ZxleyZEk8PDxwc3MjPT2d1NRUNm7ciK+vLxcvXkRECAgIoEqVKgQHB9O0aVOaNm1K3bp1cXd3N/tSiyXXr1+nevXqbN68mXr16tkkTkdcN9ahzC4rd7JI9hNPPEFAQABvv/32bevVZOfw4cO0atWKgwcP4uPjkxlusVg4f/48cXFxnDlzhuvXr5OWlkZ6ejoeHh5ER0ezYMECIiIi8Pf3p0yZMk6zLmtxYuLEiSilmDFjhk3ic0SzK8pi7GmgapbPgdYwmxIdHc2qVas4fPiwraMu1tSuXZt+/foxY8YMZs6cmRnu5uaGn58ffn5+NG7c+KbzNm/ezLBhw6hbt25RytUUkj59+vDEE0/YzOwckaKsFNkG1FZK3aOU8gIGAuG2TmTy5Mm88sorlC1b1tZRF3umTJnCl19+SUxMjBHw1VdQowa4uRl/v/oq2/EiQnh4OD179ixyrZrCERoaSkJCAkePHjVbit0oMrMTkTTgWWAdsB/4RkSibZnGH3/8wfbt2xkzZowto9VYCQgIYOTIkUydOtUwthEj4ORJEDH+jhiRzfCio42vt2HDhiYp1hQUNzc3evToQUREhNlS7EaRNneJyGoRqSMitUTE5hVqM2bMYPLkyZQsWdLWUWusTJgwge+//57Tr7wCSUnZdyYlweTJmR/Dw8MJCwvTdXROQlhYGOHhNi9sOQwu07Z/4sQJfvvtNx599FGzpbg05cqVY+DAgXxxOo/q1r/+yvx31apV9OjRo4iUae6Ujh07sm3bNq5cuWK2FLvgMmY3e/Zshg4dSunSpc2W4vKMGjWK2e7upOa2s1o1ANLS0ti5cyctW7YsUm2a26dUqVI0atSIHTt2mC3FLriE2V2/fp25c+cyatQos6UUCxo1akSte+9lhZdX9h2lSoG1u8++ffuoVq0aZcqUMUGh5nZp2rQpUVFRZsuwCy5hditWrKBhw4a6e0MRMnrqVD6vUweqVweljL+zZ8OQIQBERUXRtGlTk1VqCos2Owfn+++/Z+DAgWbLKFb07NmTrSdPcnHHDrBY4MSJTKMDbXbOijY7B+bGjRusXbtWV4QXMaVLl6Zdu3asXbs21/3a7JyT++67j5iYGJdspHB6s/v111+pU6cO/v7+ZkspduTXVeHIkSO6WsEJ8fDwoFatWhw/ftxsKTbH6c0uoy+Xpujp0aMHa9euJTU1e7vsjRs3uHz5Mr6+viYpKxrc3d0JDg6mYcOG9OvXjyRrv0OlVLYuUGlpafj6+mYrfaxdu5bmzZtTr149goODGTBgAH9l6bZjJv7+/sTGxpotw+Y4vdlt3LiRrl27mi2jWOLv788999zDtm3bsoWfOXOGSpUqufwUTd7e3uzcuZO9e/fi5eXFZ599BhhF/L1795KcnAzAhg0bqFKlSuZ5e/fu5bnnnmP+/PkcOHCAnTt3MmTIEE6cOGHGZdxEQEAAcXFxZsuwOU79NCYmJnLixAmCgoLMllJsadas2U0V2rGxscWuWqFt27YcOXIk83O3bt1YtWoVAEuWLGHQoEGZ+9577z0mTZrEfffdlxkWFhbGAw88UHSC80Hn7ByQnTt30qBBAzw9Pc2WUmzJrfUuLi6OgIAAkxQVPWlpaaxZs4ZGjRplhg0cOJClS5eSkpLC7t27adGiRea+6OhomjRpYobUAqFzdg6IbvEzn9zMLj4+Hj8/P5MU2Y/AQKNLYcZ27VoiSkGpUheoVq0aTz31VOaxQUFBnDhxgiVLltCtW7c847xw4QLBwcHUqVMn29RZZuLn50d8fLzZMmyONjvNHdGwYUOOHj2aWTkPRgPFraZqd0bCwiDnoBEvL3j66Up88skneOXYGRYWxvjx47MVYQEaNGjA9u3bAfDx8WHnzp2MGDGCxMREu+ovKCVKlODGjRtmy7A5Tm12hw8fzlbv4cycOXOGgQMHUqtWLZo2bUq3bt04dOgQ0dHRdOjQgbp161K7dm2mTZtG1tml16xZQ2hoKPXr1yckJISXXnopc9/s2bOpV68e9erVo3nz5mzZssXmukuUKEH16tWzdVVIS0vDw6Mo54UtGqZMMabuy4q7uxGeG08++SRvvPFGtuItGDPHvP322+zfvz8zLCnnDDIm4uHhQVpamtkybI5Tm52r1A2JCA8//DDt27fn6NGjREVFMWPGDOLj4wkLC2PixIkcPHiQXbt28fvvvzNr1izAaNV79tlnWbRoEfv27SMyMpJ7770XgJUrV/L555+zZcsWDhw4wGeffcbgwYM5c+aMzfXnrOOxWCwu2RLr7w9PPJE1d5fCE09A5cq5Hx8YGMjYsWNvCm/UqBEfffQRjz32GHXr1qV169bs37+fwYMH2017YXBzc8NisZgtw/aYvbxZXtutllK0WCzi5eXlEsvsbdq0Sdq2bXtT+BdffCFDhw7NFnbkyBEJDAwUEZGhQ4fK3Llzc42zTZs2smnTpmxhr732mrz22ms2Uv03Q4YMkfnz52d+/r//+z95/vnnbZ6OIxAbK1KypLEIqbe3SFyc2YpsT0REhHTv3v2O4sABl1J02tdvQkICpUuXdomJOvfu3Ztr3WN0dPRN4bVq1SIxMZErV67keV5e54aGhmbOHmxLcubsXLUYBH/n7tzcyDdX58y4ajWE015RXFycU/flCgyEv+e/NIo6H34IVarAqVNmqbo9/P39s3WILVOmDJcvXzZPkJ2ZMgXWrcu7rs7ZuXz5sktOzeW0ObvExESn/kLyatnr1evvz/Xr17+pW8exY8e46667uPvuu2nQoEGeM1Tkdm5UVJRd1tEtU6YM165dy/zsqv20MvD3h6NHXTNXB65TF54TpzU7Z89q59ayJ5KWmVvYvXs3devWZcuWLWzcuBGA5ORkxo4dy4QJEwB4+eWXeeeddzh06BBgNAxkDFmaMGECr7zyChcuXACMDtjz5s1j9OjRNr+WnMVWV+2BX1xw1REwTusW6enpTt3il1H3M3cu3LgBXl5C1aobad16DCVLlqRGjRp8+OGHrFixgueee44xY8aQnp7O0KFDefZZY13xoKAgPvzwQwYNGkRSUhJKqczB5mFhYZw+fZr7778fpRRlypRh0aJFdnmI3dzcSE9Pz/zs6jk7VycuLo42bdqYLcPmOK3ZeXh4ZPuBOSNTpsCXXxr/u7srtmzpSuXKN6/buXnz5jzj6NGjR55z+Y0aNapIpqpPT0/PlssuX748ycnJJCcn4+3tbff0NbbFVXN2Tps18vDwcPpe3q7SspeamprN7JRSVKlS5e/FtDVORUxMTLZZWlwFpzU7Hx+fzPooZ2bKFKhRw7lb9s6fP4+Pj0+2sKCgIHbu3GmOIM1tc+HCBS5fvkyNGjXMlmJznNbs/P39iYuLyzZ0yhlxhZa93LoBufJaBq7M9u3bCQkJcer68Lxw2isqXbo0Xl5eLt2fy1nIrauCNjvnxJUn13BaswPdxcFRyK1CO8PsnD3nXdyIjIzUZueIVKlShVPONtzABTl16tRNOTs/Pz/KlCnD0aM3ty5rHBeds3NQGjZsyK5du8yWUaxJSEjg0qVLuVZot2vXLrNDtMbxOXLkCNevX6d27dpmS7ELTm12ul7IfLZv305wcHCuFdo9e/YkIiLCBFWa2yEiIoIePXq4ZOMEOLnZhYaGarMzmaioKEJDQ3Pd16VLF3799VeHmYFXkz+uviypU5td3bp1OXPmDJcuXTJbSrElvwrtsmXL0rJlSzZs2FDEqjSFJSEhge3bt9OxY0ezpdgNpzY7d3d3mjRpwu+//262lGKJiPDHH3/QvHnzPI8JCwsjPDy8CFVpboc1a9bw4IMPuvTwPqc2O4Du3buzcuVKs2UUS3bs2EGpUqUyp4LPjd69exMeHp5tCiiN47Fo0SL69u1rtgy7YjOzU0r9Uyl1QCm1Wyn1vVKqXJZ9ryqljiilDiqlutgqTTByDhEREbo/lwlEREQQFhaGUirPYwIDA2nbti2LFy8uQmWawnD06FEiIyPp16+f2VLsii1zdhuAhiISBBwCXgVQStUHBgINgK7ALKWUu60SrVu3Lt7e3nocpgkUtEJ79OjRzJo1S7+QHJTPPvuMJ554wiWWOMgPm5mdiKwXkYwZHP8EAq3/9wKWish1ETkOHAHyruQpJEopwsLCWLFiha2i1BSAmJgYTp48yf3333/LYzt16kRiYiJ//vlnESjTFIbk5GTmzZvHyJEjzZZid+xVZ/cksMb6fxUg61w/p6xhN6GUGqGUilRKRZ47d67AiQ0YMIAFCxY4/fx2zsT8+fPp06dPgWaLdnNzY9SoUZlLQGoch2+++YZmzZpRq1Yts6XYnUKZnVJqo1Jqby5bryzHTAbSgK8KK0ZEZotIqIiE+vr6Fvi80NBQKlSowLp16wqbpOY2SEtL4/PPPy/UFO/Dhg1jzZo1HDt2zI7KNIUhPT2d999/P9e1bV2RQpmdiHQSkYa5bCsAlFLDgB7AEPm7guY0UDVLNIHWMJuhlMqsF9LYn5UrV1KtWjWCg4MLfE6FChUYO3YsU5x54j4XY8GCBfj4+NCli03bDB0XWy1Ai9H4sA/wzRHeANgFlADuAY4B7reK71aLZOfk2rVr4uPjI8eOHSvUeZrC07lzZ1m0aFGhz7ty5YpUqlRJduzYYXtRmkKRnJwsVatWld9++80u8ePii2T/GygDbFBK7VRKfWY102jgG6sRrgXGiIjNK9dKlSrFsGHD+PDDD20dtSYLu3fvZvfu3TzyyCOFPrdMmTJMnjyZSZMm2UGZpjDMmjWLkJCQAjUwuQxmu21eW2FzdiIiZ86ckQoVKsjx48cLfa6mYHTv3l0++uij2z7/+vXrcs8998hPP/1kO1GaQnHp0iXx9fWVvXv32i0NXDxnZzqVKlVizJgxvPHGG2ZLcUl+/fVXoqOj76ibgpeXF++//z6jRo0iJSXFhuo0BeXll1+mT58+dlkw3ZFxKbMDGD9+PGvXrmXPnj1mS3EpRISJEyfy1ltvUaJEiTuK65FHHqFhw4b6pWQC69evZ/369bz//vtmSylyXM7s7r77biZOnMiECRMyGkg0NuC7777jypUrDB482Cbxffrpp8yfP193NC5CLl++zNNPP82cOXO4++67zZZT5Lic2QGMGTOG2NhYFi5caLYUl+DChQs899xzzJo1C3d324z08/Pz4+OPP+aJJ57Qxdki4uWXX6ZLly507tzZbCmm4JJm5+Xlxbx58xg/frxekMcGPP/88/Tv35+2bdvaNN7+/fvTsGFDXnnlFZvGq7mZ8PBw1q1bx8yZM82WYhouaXYAISEhjBo1ipEjR+ri7B2wYsUKtm7dyjvvvGOX+GfPns2aNWv48ssv7RK/BqKjoxk+fDjffPNNsSy+ZmJ2c3Be2+10PcnJ9evXJSgoSObMmXPHcRVH4uLiJCAgQH755Re7prN//37x9fW1WwfX4sz58+elVq1asmDBgiJNF931pGjx8vJi6dKlTJo0iT/++MNsOU7F9evXefjhhxkxYoTNi685qVevHvPnz6dfv37ExMTc+gRNgUhNTaV///706dOHoUOHmi3HdFza7ADuu+8+vvzySx555BG9xmwBERGeeeYZAgMDi2ws6z/+8Q/GjRtH7969uXr1apGk6cqICOPGjaNEiRLMmDHDbDkOgcubHRhTt48dO5bevXuTnJxsthyH58MPP2THjh3MmzevSJfVGz9+PM2aNaNHjx4kJSUVWbquhohklmaWLFlisxZ0Z6dYmB3AhAkTqFevHoMHDyY1NdVsOQ7L8uXLef/991mxYgWlS5cu0rSVUsyaNYsaNWrQq1cvbXi3gYjw5ptvsnLlStavX0/ZsmXNluQwFBuzU0oxd+5cbty4wdChQ0lLS7v1ScWMlStXMnr0aFavXk316tVN0eDm5sbcuXOpVKkSPXr00GvOFgIR4bXXXmPZsmVs2LCBihUrmi3JoSg2ZgdQokQJli9fTkJCAkOHDtU5vCysWLGCJ598koiICEJCQkzV4uHhwfz586lZsyZdunTh/PnzpupxBtLT0xk3bhyrV69m8+bNVK5c2WxJDkexMjuAkiVLEh4ezpUrV+jfv7+uwwOWLl3KiBEjWL16db5rwBYl7u7uzJ49mwceeIDmzZvrsc75cOnSJXr06EF0dDQ//vijztHlQbEzOzAM77vvvqNUqVK0bdu22LbSWiwWXn/9dV555RU2bNhAaGio2ZKy4ebmxowZM5g+fTodOnTghx9+MFuSw3Hw4EFatGhBnTp1WLt2LeXLlzdbkuNidke/vDZbdCq+FRaLRd59910JCAiQP/74w+7pORJXr16V3r17S+vWreXMmTNmy7kl//vf/yQwMFDeeustsVgsZstxCNasWSO+vr7yxRdfmC3lJnDATsWmC8hrKwqzyyAiIkJ8fX1l7ty5xeKHdPjwYWnUqJE89dRTkpKSYracAhMbGyutWrWSjh07yokTJ8yWYxrXrl2TcePGSUBAgGzZssVsObniiGZXLIuxOenRowebN2/mww8/pHfv3sTFxZktyS5YLBY++ugjWrZsyTPPPMOcOXPueG66osTf359ffvmFzp07Exoayueff268sYsRW7ZsoXHjxpw9e5bdu3fTunVrsyU5D2a7bV5bUebsMrh+/bq89tpr4uvrKwsXLnSpXN7hw4elbdu20rp1azl06JDZcu6Y6OhoadasmXTs2LFYTMOfkZvz9/eX77//3mw5twQHzNmZLiCvzQyzyyAyMlIaNmwo3bp1k+joaNN02ILExESZPn26+Pj4yAcffCBpaWlmS7IZqampMmPGDKlQoYJMnDhREhISzJZkc1JTU+WLL76QwMBAGTx4sJw/f95sSQVCm52TmJ2ISEpKisycOVN8fX1l2LBhcvLkSVP1FJYbN27IrFmzxN/fX/r37y+HDx82W5LdiImJkaeeekoqVqwo7733niQlJZkt6Y6xWCyyfPlyqVevnrRr187pGtC02TmR2WVw8eJFmTRpkpQvX17GjRvn8KaXkpIi8+fPl1q1aknnzp0lMjLSbElFxr59+6RPnz5SpUoV+fe//y1XrlwxW1KhSUtLk4iICGnRooU0btxYVq9e7ZTVKdrsnNDsMoiNjZVx48ZJhQoVJCwsTNauXSvp6elmy8rk+PHjMnHiRPHz85POnTvLxo0bzZZkGn/++af06dNHypcvL2PGjLHrkoG24uzZszJjxgypXr26NG/eXBYvXuxQz1dh0WbnxGaXQWJiosyePVuCg4Pl3nvvlWnTpsnu3btNefueP39eFixYIN27dxcfHx954YUX5ODBg0Wuw1GJiYmR119/Xfz9/aVdu3aycOFCh6rXu379uqxfv16GDBki5cqVkyeffFK2bdtmtiyb4IhmpwxdjkdoaKhERkaaLSNPRIQ///yTr7/+mhUrViAihIWF0bNnT1q1asVdd91l8zTT09M5cOAAa9euJTw8nJ07d9KhQwd69+5Nv379KFWqlM3TdAVSU1P54YcfWLBgAT///DOhoaGZ31WtWrWKVEtCQgJr1qwhPDyc9evXU69ePfr378/jjz9OhQoVilSLPVFKRYmIQw3J0WZnA0SE6OhowsPDWblyJbt27aJatWqEhobStGlTGjduTGBgIAEBAQWaNiktLY34+Hji4uI4cOAAUVFRREZGsnPnTvz8/OjUqRNhYWF06NABb2/vIrhC1yEpKYmNGzdmflfly5enRYsW2b4rW700Ml5OkZGRmd/h3r176dChA2FhYXTv3p1KlSrZJC1HQ5tdIXAms8tJamoq+/btIyoqiqioKHbv3k1sbCyxsbF4enoSEBBAxYoV8fT0xN3dHYvFQlpaGlevXiUuLo6EhAR8fHzw9/fn3nvvzfwhNmnSRI99tCEWi4WdO3cSGRmZaUj79++nVq1a1K5dG39/fwICAjL/+vr64uXlhYeHB0op0tLSSEtL4+LFi8TGxhIXF5f59+TJk+zevRt/f//M769p06a0bNmyWLygtNkVAmc2u7wQES5fvkxcXBznz5/P/LG4u7vj4eFB6dKl8ff3x8/PDw8PD7PlFkuuX7/O3r17OX78eDbziouL4+zZs5nfmcViyXxZlS9fPpsx+vv7U7VqVYKCgihXrpzZl2QK2uwKgSuanUZTXHBEs9NjYzUaTbFAm51GoykWaLPTaDTFAm12Go2mWKDNTqPRFAu02Wk0mmKBXcxOKfWSUkqUUhWtn5VS6mOl1BGl1G6lVBN7pKvRaDR5YXOzU0pVBR4C/soS/A+gtnUbAfzH1ulqNBpNftgjZ/d/wAQga2/lXsAC64QIfwLllFL+dkhbo9FocsWmY5KUUr2A0yKySymVdVcVICbL51PWsGwr2yilRmDk/AASlVIHbanPSkXAmZaYdya9zqQVnEuvM2kFqGu2gJwU2uyUUhuByrnsmgxMwijC3hYiMhuYfbvnFwSlVKSjDWPJD2fS60xawbn0OpNWMPSarSEnhTY7EemUW7hSqhFwD5CRqwsEtiulmgOngapZDg+0hmk0Gk2RYLM6OxHZIyJ+IlJDRGpgFFWbiMgZIBx4zNoq2xK4LCKuuTirRqNxSIpqHqHVQDfgCJAEPFFE6eaGXYvJdsCZ9DqTVnAuvc6kFRxQr8NO8aTRaDS2RI+g0Gg0xQJtdhqNpljg0manlHpBKRWtlNqrlFqilCqplLpHKbXVOnTta6WUl4n6/quUOquU2pslrIJSaoNS6rD1b3lruOlD7vLQ+0+l1AGrpu+VUuWy7HvVqvegUqqL2Vqz7HO44Yx56VVKPWe9v9FKqfezhDvUvVVKBSul/lRK7VRKRVp7YTjEvc3E7LUc7bVhdFo+DnhbP38DDLP+HWgN+wwYZaLGB4AmwN4sYe8DE63/TwTes/7fDVgDKKAlsNVB9D4EeFj/fy+L3vrALqAERpeko4C7mVqt4VWBdcBJoKKD39sHgY1ACetnP0e9t8B64B9Z7udmR7m3GZtL5+wwWpu9lVIeQCmMERsdgGXW/fOB3uZIAxH5BUjIEdwLQxdk12f6kLvc9IrIehFJs378E6MPJRh6l4rIdRE5jtES39xMrVYccjhjHnpHAe+KyHXrMWet4Y54bwW42/p/WSDW+r/p9zYDlzU7ETkNzMSYkCAOuAxEAZey/Dgzhq05EpXk7z6IZ4CMhUXzGnLnSDyJ8RYHB9SbdThjjl0Op9VKHaCttdrlZ6VUM2u4I+odB/xTKRWD8bt71RruMFpd1uysdV29MLL5AUBpoKupogqJGOUAp+gbpJSaDKQBX5mtJTeUUqUwhjO+braWQuABVMAo/r0MfKNyDDp3IEYBL4hIVeAFYK7Jem7CZc0O6AQcF5FzIpIKfAe0xshGZ3SmdsRha/EZ2Xzr34yii8MOuVNKDQN6AEOsBg2Op7cWfw9nPMHfwxkr43haMzgFfGctAv4PsGBMCOCIeh/H+I0BfMvfxWqH0erKZvcX0FIpVcr6NuwI7AN+Ah6xHvM4sMIkfXkRjqELsutzyCF3SqmuGHVgYSKSlGVXODBQKVVCKXUPxlyG/zNDIzjtcMYfMBopUErVAbwwZj5xqHtrJRZoZ/2/A3DY+r/j3FuzWkaKYgPeBA4Ae4GFGK1XNTEejCMYb6ASJupbglGfmIrx43sK8AE2YTwsG4EK1mMV8ClGy9seINRB9B7BqJPZad0+y3L8ZKveg1hb6szUmmP/Cf5ujXXUe+sFLLI+v9uBDo56b4E2GHXiu4CtQFNHubcZmx4uptFoigWuXIzVaDSaTLTZaTSaYoE2O41GUyzQZqfRaIoF2uw0Gk2xQJudRqMpFmiz02g0xYL/B53a3dv68+vvAAAAAElFTkSuQmCC\n",
"text/plain": [
"