{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# 6 - Sísmica de reflexão: Correção de Normal Moveout (NMO) e empilhamento\n",
"\n",
"Na aula passada, vimos que uma seção CMP é crucial para conseguirmos uma estimativa da velocidade. Nessa prática, veremos como utilizar a correção de NMO para fazer a **análise de velocidades**. Também veremos a técnica do empilhamento e como ela melhora a razão sinal-ruído do dado.\n",
"\n",
"\n",
"Utilizaremos as simulações de ondas da biblioteca [Fatiando a Terra](http://www.fatiando.org). Essas simulações utilizam o [método de diferenças finitas](http://en.wikipedia.org/wiki/Finite_difference_method) para calcular soluções da equação da onda."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Objetivos\n",
"\n",
"* Utilizar a correção de NMO para estimar velocidades.\n",
"* Ver como o empilhamento melhora a qualidade dos dados."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Questão para entregar\n",
"\n",
"
\n",
"\n",
"Explique como é feita a análise de velocidades. As velocidades determinadas são as velocidades verdadeiras das camadas? Por que? Qual é vantagem de realizar o empilhamento? \n",
"\n",
"
\n",
"\n",
"### Regras para a resposta\n",
"\n",
"* Coloque **nome, data e o número da prática** em sua resposta. \n",
"* A resposta pode ter no **máximo 1 página** (não uma folha).\n",
"* **Execute o notebook** antes de responder. As simulações abaixo foram feitas para te ajudar.\n",
"* **Pense e organize** sua resposta andtes de começar a escrever."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Instruções\n",
"\n",
"Esse documento é um [Jupyter notebook](http://jupyter.org/), um documento interativo que mistura texto (como esse), código (como abaixo), e o resultado de executar o código (números, texto, figuras, videos, etc).\n",
"\n",
"O notebook te fornecerá exemplos interativos que trabalham os temas abordados no questionário. Utilize esses exemplos para responder as perguntas.\n",
"\n",
"As células com números ao lado, como `In [1]:`, são código [Python](http://python.org/). Algumas dessas células não produzem resultado e servem de preparação para os exemplos interativos. Outras, produzem gráficos interativos. **Você deve executar todas as células, uma de cada vez**, mesmo as que não produzem gráficos.\n",
"\n",
"**Para executar uma célula**, clique em cima dela e aperte `Shift + Enter`. O foco (contorno verde ou cinza em torno da célula) deverá passar para a célula abaixo. Para rodá-la, aperte `Shift + Enter` novamente e assim por diante. Você pode executar células de texto que não acontecerá nada."
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Setup\n",
"\n",
"Rode as células abaixo para carregar os módulos necessários para essa prática."
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/home/leo/miniconda3/envs/geofisica2/lib/python2.7/site-packages/h5py/__init__.py:36: FutureWarning: Conversion of the second argument of issubdtype from `float` to `np.floating` is deprecated. In future, it will be treated as `np.float64 == np.dtype(float).type`.\n",
" from ._conv import register_converters as _register_converters\n"
]
}
],
"source": [
"%matplotlib inline\n",
"from __future__ import division, print_function\n",
"import math\n",
"from multiprocessing import Pool, cpu_count\n",
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import ipywidgets as ipw\n",
"from fatiando.seismic import RickerWavelet, FDAcoustic2D\n",
"from fatiando import utils\n",
"from fatiando.vis import anim_to_html\n",
"from fatiando.vis.mpl import seismic_image, seismic_wiggle"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Simulação de um CMP para um modelo de duas camadas\n",
"\n",
"Rode as células abaixo para simular uma seção CMP para um modelo de duas camadas. A célula abaixo cria nosso modelo."
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"shape = (150, 200)\n",
"spacing = 10\n",
"extent = [0, shape[1]*spacing, shape[0]*spacing, 0]\n",
"\n",
"# Velocidades\n",
"v1, v2 = 4000, 5000\n",
"\n",
"densidade = np.ones(shape)*1600\n",
"velocidade = np.ones(shape)*v1\n",
"l1 = 40\n",
"densidade[l1:,:] = 1800\n",
"velocidade[l1:,:] = v2\n",
"l2 = 100\n",
"densidade[l2:,:] = 2000\n",
"velocidade[l2:,:] = 8000"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Em seguida, precisamos definir onde serão localizadas as fontes e os receptores em nossa simulação. Vamos aproveitar e calcular também os espassamentos (offsets) dos receptores. Lembre-se: offset é a distância da conte ao receptor."
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Utilizando 14 fontes e 14 receptores.\n",
"Fontes: [940 910 880 850 820 790 760 730 700 670 640 610 580 550]\n",
"Receptores: [1060 1090 1120 1150 1180 1210 1240 1270 1300 1330 1360 1390 1420 1450]\n",
"Offsets: [120 180 240 300 360 420 480 540 600 660 720 780 840 900]\n"
]
}
],
"source": [
"fontes = np.array(list(reversed(range(55, shape[1]//2 - 3, 3))))\n",
"recep = np.array([shape[1] - s for s in fontes])\n",
"offsets = (recep - fontes)*spacing\n",
"print(\"Utilizando {} fontes e {} receptores.\".format(len(fontes), len(recep)))\n",
"print('Fontes: {}'.format(fontes*spacing))\n",
"print('Receptores: {}'.format(recep*spacing))\n",
"print('Offsets: {}'.format(offsets))"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Vamos rodar as simulações que precisamos (uma por fonte). **A barra de progresso não irá aparecer** pois vamos rodar as simulações em paralelo para agilizar o processo."
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [],
"source": [
"def simul_shot(fonte, its=1200, verbose=False):\n",
" shot = FDAcoustic2D(velocidade, densidade, spacing=spacing, taper=0.005, padding=50, verbose=verbose)\n",
" shot.add_point_source((0, fonte), RickerWavelet(1, 120))\n",
" shot.run(its)\n",
" return shot"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Simulando... Pronto.\n",
"CPU times: user 32.2 ms, sys: 34.3 ms, total: 66.5 ms\n",
"Wall time: 29.9 s\n"
]
}
],
"source": [
"%%time\n",
"print('Simulando...', end=' ')\n",
"pool = Pool(processes=cpu_count())\n",
"shots = pool.map(simul_shot, fontes)\n",
"pool.close()\n",
"print('Pronto.')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Vamos dar uma olhada na simulação de 1 tiro para ter uuma ideia do que acontece."
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
" "
],
"text/plain": [
""
]
},
"execution_count": 6,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"shots[0].animate(embed=True, dpi=60, fps=8, every=20, cutoff=0.5, cmap='Greys')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Vamos extrair os dados CMP da nossa simulação."
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAsgAAAIHCAYAAACCFaoXAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMi4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvhp/UCwAAIABJREFUeJzt3XmcZWV97/vvr6onuht6ZmpmwQCiojIZCaJBr0YP6FUjaqLJTUJUfKl4PA4Z8BxyY+ToK1FPOB4RiYmJGgVJOoIiiBiHqzI0YXJuFJpGgW7oqXru5/5Ru822abqqnl17/fa31uf9evGSqt5V/XWx9lrf9axnPTtKKQIAAAAwaig7AAAAADBIKMgAAABAFwoyAAAA0IWCDAAAAHShIAMAAABdKMgAAABAFwoyAAAA0IWCDAADKiJeHRE3R8SGiHggIr4YEadHxH+PiBIRb97t9W/tfP+/d74+MyJ2dn5+fUT8ICJ+P+X/DAAYoSADwACKiLdJ+qCk90o6QNJhkv63pHM6L/mhpNft9mOv7Xy/26pSylxJ+0l6p6SPRcTx/coNAFMBBRkABkxEzJN0kaTzSymfL6VsLKVsK6X8Wynlv3VedpOk2RHxpM7PPEnSPp3vP0YZ9S+SHpFEQQaAvaAgA8DgeaakWZKuGuN1n9ToqLE0Opr8D4/3wogYioiXSpov6Y7JCAkAUxUFGQAGzyJJD5dSto/xun+U9KqImC7p3M7Xuzs4Ih6V9LCk90j63VLKDyY1LQBMMdOyAwAAHmO1pMURMW1vJbmUcm9E/Fij85R/VEq5LyJ2f9mqUsohfcwKAFMOI8gAMHj+P0mbJb1kHK/9B0n/VXuZXgEAmBhGkAFgwJRS1kbEhZIuiYjtkr4saZuksyQ9R9JI18v/WdJKSd9sPCgATFGMIAPAACql/LWkt0n6M0kPSbpP0psk/ctur9tUSrm+lLKp+ZQAMDVFKSU7AwAAADAwGEEGAAAAuvS1IEfECzofbfrjiHjXHv58ZkT8c+fPvxMRR/QzDwAAADCWvhXkiBiWdImkF2r0U5tetYePN/0DSY+UUo6W9DeSLu5XHgAAAGA8+jmCfIqkH5dSVpRStkr6jKRzdnvNOZL+vvPvV0j6zdjDIp4AAABAU/pZkJdq9KnrXVZ2vrfH13QWw1+r0U+QAgAAAFL0cx3kPY0E775kxnheo4g4T9J5kjRjxoxnHHjggb2nAwAAwJS0evVqbdiwoXpWQj8L8kpJh3Z9fYikVY/zmpURMU3SPElrdv9FpZRLJV0qSYcffnh55zvf2ZfAAAAA8Hfxxb091tbPKRY3STomIo6MiBmSzpW0bLfXLJP0us6/v1zSDYWFmQEAAJCobyPIpZTtEfEmSddKGpZ0eSnlroi4SNLNpZRlkj4u6ZMR8WONjhyf2688AAAAwHj0c4qFSinXSLpmt+9d2PXvmyW9op8ZAAAAgIngk/QAAACALn0dQe4X16WSmV4NAAAw+BhBBgAAALpQkAEAAIAuFGQAAACgi+UcZObyAgAAoF8sC7IrHi4EAAAYfEyxAAAAALowgtwgRmIBAAAGn2VBZqoCAAAA+oUpFgAAAEAXCjIAAADQxXKKhSumhgAAAAw+RpABAACALpYjyK4jmq4jyAAAAG3CCDIAAADQxXIE2RUj381z3eYAACAPI8gAAABAFwoyAAAA0IUpFhgT0xQAAECbUJAb5DqXl4IMAADahCkWAAAAQBfLEWTXkVgAAAAMPsuCjGY5X5AwPQQAAEyUZUF2LWyUNQAAgMFnV5AjQtOnT8+OUWXr1q3ZEQAAADAGy4I8Y8aM7BhVKMgAAACDz64gT5s2TQsWLMiOUWX9+vXZEaq4TmmRfLMzHQcAgDx2BXn69Ok6+OCDs2NUuffee7MjVHEua64FGQAA5LEsyAcccEB2DAAAAExRlgV56dKl2TFgwnn0GwAA5OCT9AAAAIAuFGQAAACgi90UCzTP+UE3plgAAICJYgQZAAAA6MIIcoOcR2JduW5zRr4BAMjDCDIAAADQhRHkBg0PD2dHqLJjx47sCAAAAI2hIDdo2jTPzU1Bbp7r1BCJ6SEAAH+ejc3UzJkzsyNU2bJlS3YEAACAxlCQG7Rw4cLsCFXWrl2bHaF1nEeQAQBwZ1eQf/SjH+mFL3xhdowqF198cXaEKvfcc092BAAAgMbYFeR169bpS1/6UnaMKn/3d3+XHQEmmMcLAEAeu4Ls7MADD8yOAAAAgDGwDjIAAADQhRFkTGmuD7sxxQIAgDyMIAMAAABdGEEGBpDryLfE6DcAwB8FGWNyLmsAAAATRUHGmJwLMqOZzWN/AQC4oyBjTENDvlPVd+zYkR0BAACYoSBjTDNmzMiOUG3Tpk3ZEQAAgBm7gnzsscfq8ssvz47RKrNmzcqOUM21IDvf6neeYgEAgGRYkOfMmaNnPvOZ2TFaZcmSJdkRqq1ZsyY7AgAAMGNXkJ25jqx973vfy45Q7Qc/+EF2hNZxHv0GAECiIGMcDjjggOwIAAAAjaEgY0wLFizIjgAAANAY3/W7AAAAgD5gBBkYQK7z1SXmIAMA/FGQAUwqyj0AwB1TLAAAAIAujCBjSnMezUTzXPcXRr4BYHLZFeSdO3dq48aN2TFggsIDAAAmyrIgb968OTsGTAwPD2dHqLJ9+/bsCAAAtJZlQXYdQf6zP/uz7AitM3PmzOwIVZwLsvPot+sdBwDA5LIryKUUbdmyJTtGlXe/+93ZEVpn9uzZ2RGquF4EAgAwFdgVZOcpFq5l7Tvf+U52hGr7779/doQqDz74YHaEVnIe/QYATB67guw8guzqtNNOy45Q7Z577smOUOXOO+/MjgAAQGtZFmTn+Zlo1pIlS7IjAAAAM5YFedu2bdkxYGLOnDnZEQAAgBnLgswIMgAAAPrFriBLow/qARhMLJXWPB4uBIDJZVeQmWIBAACAfrIsyK4jyEwNAdAPzqP2jH4DGER2BVmSduzYkR2hyoYNG7IjwIRz4QEAwF1fC3JEvEDShyQNS7qslPK+3f78DEkflPQUSeeWUq4Y63c6jyCvX78+OwJMDA8PZ0eo5noBCwDALn0ryBExLOkSSc+TtFLSTRGxrJRyd9fL7pX0e5LePpHf7XoC5gNOMF5DQ0PZEaq5vj8BANilnyPIp0j6cSllhSRFxGcknSPplwW5lPLTzp+Ne0i4lGI7Z21kZCQ7QpXly5dnR2gd148ll6S1a9dmR6jmemxhSg4ATK5+FuSlku7r+nqlpFNrflFEnCfpPElavHix7UnM9SG9pz/96dkRWmfu3LnZEao5F2RXrsdEABhU/SzIexrSqDqKl1IulXSpJD3hCU8orrdwmWLRvOc85znZEaq8733vG/tFA2rlypXZEQAA6Ek/C/JKSYd2fX2IpFW9/lLnKRas39y8G2+8MTtClcMPPzw7QrVbbrklOwIAAD3pZ0G+SdIxEXGkpPslnSvp1ZPxi10LsuvIN5q3ePHi7AgAALRW3wpyKWV7RLxJ0rUaXebt8lLKXRFxkaSbSynLIuJkSVdJWiDpv0TE/yilPGkcv7tfsfvKdQ4ymjdtmuUS5QAATAl9PQuXUq6RdM1u37uw699v0ujUi4n+3t7DJaAgA+gH51UsXI/nAKY2hqkaREEGBptz0QQATB67guz8SXquuTdt2pQdAcAU5XxRwug3MHXZFWTJ96Dk+pAey9MBAIA2sSzIrlxHkNetW5cdAUacRwQBAJBMC7LrCLJr7o0bN2ZHgJHh4eHsCNVc7/IAACaXZUF25TqC7DzF4txzz82O0DrOS9RRkAEAkmlBdh2JpSA377LLLsuO0Dr77rtvdoRqrvu66zERAAaVZUFGs5yXp5szZ052hCrO01r222+/7AjVHn744ewIAIABYFeQnZd5cx3l2bZtW3aE1pk7d252hGq33XZbdoRqP/nJT7IjAAAGgF1BdkZBRhscccQR2RGqffOb38yOAAAYABTkBrmOfDtPsUDz5s2blx0BAICeWBZk15FY19wUZAAA0CaWBRnNch35BjD4nD9YxnXQA8DYKMgNcj2YOq8Nu2nTpuwIQN85l0xnrtvd9VwENImCjDE5jyA7L5cGAAByWBZk16tfcjePEWQAADBRlgUZzXKeYsEIMibC9ZY5AGByUZAb5DpVwXkE2bUg77///tkRWmn69OnZEaqwVjkATC4KMsbkXJBdi8OHP/zh7AitNGPGjOwIVVz3c8n7+AJg6qIgY0zOJ7CtW7dmR6jyyle+MjtCK7l+yInrnRIAGFSWBdm1sLnmdp0aInnPn3Z16aWXZkeo9sQnPjE7QpX7778/OwIATCmWBRkYLz4FsHl//Md/nB2h2g9/+MPsCFXuvvvu7AgAMKVQkDEm15FvyXtuJpp36KGHZkcAAAwACnKDXIuma26JEWRMzKxZs7IjAAAGAAUZY6IgN895mwNt4bpuNscXYGwUZExprieCzZs3Z0cAMEW5FnvJ95gOPxRkTGmuK3CMjIxkRwAa4VzWAExddgW5lGJ7Bel6InAtmZJvdta1BQAgj11BduZa7J25bvMtW7ZkRwAAoLUoyBiTa8mUGEFGOwwNDWVHqOZ8fAEwdVGQMaW5nnyd5yC//e1vz47QOjNmzMiOUI27Jc1zPS4CTaIgY0pzHUF2XZ5Okt7whjdkR2gd5/WbnQsyRROYuijImNJcT2Bbt27NjlDtqKOOyo5QbcWKFdkRqixcuDA7QrW1a9dmRwCAx6AgY0yuJdPZjh07siO00jnnnJMdocrHPvax7AjVOL4AGEQU5Aa5nghcl6eTfLc5BTnHnXfemR2hytKlS7MjVLv11luzIwDAY1CQMaW5FmTnOcho3qGHHpodAQCmFAoyMIBcHy4EAGAqoCBjTK6jsJJv0dy2bVt2hGobNmzIjgBgL5ynzTlzPpe2EQUZGECuxV7yXrYLmAiKJjB1WRZk1/Lg/GlXaJbzSIPzh5wAACCZFmRgvFyLpmtuSdq8eXN2BAAAekJBBgaQc0FmigUmgmkKAAYRBRkYQM4FeePGjdkRYGR4eDg7QjXWKwemLgoyxuRc1lyzO594nadYXH/99dkRWmfWrFnZEaq5Xgy6HheBJlGQAUwq5yXqTj/99OwIrbPPPvtkR6jmWpABjI2C3CBW38B4OY/wOBdk19HMiy66KDtCtTPOOCM7QrWHHnooOwKAPqEgA5hUztNDXL3nPe/JjlBt+fLl2REA4DEoyJjSXEdiXXNLvndKkOOwww7LjlDtW9/6VnYEAH1CQQYwqSjImIiFCxdmRwCAx6AgAwPIuWRu3749O0K1TZs2ZUcAsBesm9085zuavaAgAwPI+YDkXJDXrVuXHQEAMAAoyAAmlXO551MA0QaMwgJjoyADmFQUZACAOwoyxuQ8H5Y1nJtHQQYAuKMgAwPIuWQ627p1a3YEGGGqAjB1UZABTCrnOw4jIyPZEaqcffbZ2RFaafr06dkRqjh/2qUzBj68UJABTCrngrx58+bsCFUuueSS7AitNHv27OwIVdauXZsdoRolE02hIAMDiJNADtcpFoccckh2hGrOa0/PmTMnO0IV54IMNIWCDGBSOZd75zWcXR1//PHZEap9+tOfzo5Q5f7778+OAAw8CjKASeU8xcK53Lv66U9/mh2h2oIFC7IjAOgTCjIAdFCQMRHOU1sA7B0FGQA6duzYkR2hivM8Xmeuc5CBiWjrcoYUZExprrf7h4eHsyO0kuv+snHjxuwIAMbQ1qLpyrIgu94G5c2B8XLdxyXv7K4P6a1fvz47AtAIzqNoimVBBjC4nAuy6wgyH/wAAJOLggwAHRRkAIBEQQYAe64fcAIAg4qCDADmKMiYiKGhoewI1ZyncMELBblBrm9sHopAW7i+R52Xefvwhz+cHaF1ZsyYkR2h2pYtW7IjtI7rcbFXFGQAk8r5YOqa3bk0vPzlL8+O0DqzZs3KjlDNeV93Pb60FQUZADpcT2CbN2/OjlDtoIMOyo5Q7Stf+Up2hCoLFy7MjlBt7dq12RHQEhRkAOhwXcXC9RMA3Z111lnZEaosX748O0K1FStWZEdAS1gWZNc5sa6jUwAGG8cWTITzCDLQFMuC7PoErusoj/PJ1/ViynUkU+JjsjExzvu6q6VLl2ZHAAZeXwtyRLxA0ockDUu6rJTyvt3+/G2S/lDSdkkPSfp/Sik/G+N3UpAB9IXrxaBzyVyzZk12hNbhIhYT4TrQ1Ku+FeSIGJZ0iaTnSVop6aaIWFZKubvrZcslnVRKGYmIN0j6n5JeOdbvdi3IrjuZa2lADvaX5jkX5I0bN2ZHAIDH6OcI8imSflxKWSFJEfEZSedI+mVBLqV8tev135b0O2P90qGhIdslarZv354doQoj38Bgcz22SNLIyEh2BKARroNkbdXPgrxU0n1dX6+UdOpeXv8Hkr441i+NCM2cObPHaDlcTwQUZLSF6+i38wiyc7lH81zfoxIF2U0/C/Ke9oQ97tkR8TuSTpL07Mf58/MknSdJ+++/v+bNmzdZGRvleiuRExgw2JxLAxfgmAhKJprSz4K8UtKhXV8fImnV7i+KiLMk/amkZ5dS9vgROaWUSyVdKknHHXdc2XfffSc/bQP22Wef7AhVnD+EwPXky0kgh3PRdLVt27bsCDDCsRFN6WdBvknSMRFxpKT7JZ0r6dXdL4iIp0n6qKQXlFIeHM8vHRoa0pw5cyY7ayNc505zQALQL84F+dBDDx37RZhUzitwuA7YtFXfCnIpZXtEvEnStRpd5u3yUspdEXGRpJtLKcskvV/SXEmf65Swe0spZ+/t9w4PD9tOsXDNvXr16uwI1ZgegjZwHvl2fTZDkq688srsCK3jOtAk+U6zdD6+9KKv6yCXUq6RdM1u37uw698n/DmdQ0NDmj179iSka57ryPf06dOzI1Rznh4CtIHze/SMM87IjtA6rud/ybcgt5XdJ+kNDQ3ZzuWdP39+doQqrsVektavX58doXWcV1RwvX3rfOt2y5Y9PnpiYcaMGdkRqrzxjW/MjlDtla8c86MSBtaDD45rJikGhF1Bjgjbk5jrFAvX3JL00EMPZUeo4lx4nLneSnTNLXlfULn6yEc+kh2h2lve8pbsCNXuuuuu7AiYALuCLPmeDFzXb3a+peX6qYsUZAB4rEWLFmVHQEvYFeSdO3fazlljakjzXB/ocH6yH5gI51VymFPavMWLF2dHQEvYFeRt27Zp1arHLKds4eijj86OUMW5ILtelDB3Gm3hfLfk4Ycfzo4AoE/sCvLWrVt17733ZseocuSRR2ZHqOJaMiXv6SFAGzjPQXZeog4YL+e7PL2wK8jr1q3T9ddfnx2jylOe8pTsCFWc53y5jn67PogqeY8IunJ9LkPyXqvcdbofcri+TynIJkZGRrR8+fLsGFVcl3jZf//9syNUc12izvXhQomCjIlxHkF2LTwAxmZXkNevX68bbrghO0aVb37zm9kRqhx//PHZEaotWLAgO0IV1/VVJR4wzOBc1MiOtmjrSKwru4Ls7Mtf/nJ2hCqvec1rsiNUc13D2XnuNE/2YyKcSwMXg5gI1zuDbb0QpCA36Fvf+lZ2hCpr1qzJjlBt6dKl2RGqOBdkoC22bt2aHaHK+eefnx2hlaZN86xcbb0Q9PyvZcp1SaBbbrklO0I115VDXB8ulKT77rsvOwLQCNe7Je95z3uyI7SS6zMxjz76aHaEFBRkjOnWW2/NjlDt5S9/eXaEKq4HUsn7lrnrA4bOq544377dsmVLdoQqS5YsyY5QbcWKFdkRqrke1ynIwOP4/Oc/nx2hmutIycKFC7MjVHMua64F2ZnzKhauBdnZySefnB2h2tVXX50docrKlSuzI6SwK8gzZszQIYcckh2jiuuV7w9/+MPsCNVcl9ZzLsiuH+8t+c4pdR6FdS7IztvdlfMzMQcffHB2hCrOd5F7YVeQ58yZY3sF6VqQnble+R544IHZEarNnDkzOwLQCNfpRBT7HM4DH21kWZBPPfXU7BhVPvvZz2ZHqOJ8MP3iF7+YHaHKr//6r2dHqOa6tJ4kPfTQQ9kRYMR12S7XB8bdzZ07NzsCJsCuIM+ePVtPfepTs2NU+a3f+q3sCFVc501Jvhclf/EXf5EdoZrzChyuhceZ6yis5Ds9hAtBYGx2BXloaEj77rtvdowqxx13XHaEKs4F2XX+tPO6k67vT8m3rDnf5XHmWpA3bdqUHQEYeJYF2XWO4ymnnJIdASZWr16dHaGa8zw71xFk59U3nMu964Xs9u3bsyPAiOt7tNfcdgW5lGJ71X7AAQdkR6ji/MDYz3/+8+wIVVwfLpR875RIvp90RUHO4Xoucr1TAjTJ7mxQSrG9PbTffvtlR6jyohe9KDtCtY9//OPZEapcc8012RGqPeMZz8iOUM11iTrXY6LkvW62a0EGJsL1gqrX3HYFeefOnRoZGcmOUWWfffbJjlDltNNOy45QzbUgf/rTn86OUO3CCy/MjlDN9SL2kUceyY7QSq5TclynhkhMVcwwffr07AhVWlmQXT+9yHXe1+GHH54doXW+//3vZ0doJdePggUmwvUDcSTpggsuyI7QOhRkE6UUbd68OTtGFdeR70WLFmVHqOb6wJjzp0U5j0657i+uI5nuXG89b9iwITtCtZe97GXZEaq5dhfX9Zt7PS5aFmTXeV+uV+2uDy5Jsv1QGdcPOJGk9evXZ0eo5rpEnWtRc+f6gKHrXVjJdzRTkj7xiU9kR6hyxBFHZEeo0uvzDZbNx/Wg5PogjfMT8meccUZ2hCrOBfn222/PjlDN9QFD5xFk1wEPyffY6DpY4+68887LjlBl+fLl2RGq9Dq4Z1mQXQ+orreeXYu9JJ144onZEVrnc5/7XHaEas9+9rOzI1RxHlVzve0s+Q7WcMcBE7F06dLsCFVmzJjR08/bFeRSiu1ByTW382jDggULsiO0zg033JAdoZprcaAgYyJc93PJe8DG1fz587MjVGnlFAvX21qu6306H5BcH7pyXRJQYgWODK5zpyXvOeuunAvyAw88kB2hdVwvwFu5ioXrcmmuBdl1aojkOx3npS99aXaEap/61KeyI7SO64WgJK1atSo7Quu43s2UpNWrV2dHQEtYFmTXEWTXg5JryZR8y/2xxx6bHaGVXG/3uy7DhBwc04Gx2RVkybdour6xe53onsl1eojragrufvKTn2RHqLJ48eLsCNWcb/e7cj0XAU2yLMiuXIs9ow3Nc75l7uzWW2/NjlDlJS95SXaEVnI9prvehZW8lzSEFwpyg1wPSq5zviXfBfGXLFmSHaGVPvnJT2ZHqPK7v/u72RGquT6bIfkeG12LPdAkCnKDXA9KricByTe7674iSQcddFB2hGrXXXdddoTWmTlzZnaEaq7HF+dpLc7Ljv7Jn/xJdgRMwIQKckTMkbS5lOI5FIoqzmXNdXqI64lXkp71rGdlR6h2xRVXZEdoHeeCvHHjxuwIVZynKTgX5Ne85jXZETABey3IETEk6VxJr5F0sqQtkmZGxEOSrpF0aSnlR31PuRvXwuZ61e5aMiXfoum6moIk/cZv/EZ2hGoU5OY5f5jPmjVrsiO0zrp167IjVDvrrLOyI1S5+eabsyNUGRkZ6ennxxpB/qqk6yW9W9KdpZSdkhQRCyU9R9L7IuKqUso/9pRigijIzXIuyMz7bt5RRx2VHQFG5syZkx2hmusx3ZnzCLKriy++ODtClZ///Oc9/fxYBfmsUspjlgEopayRdKWkKyPC8yNWElDsm+dakJ1HkA8++ODsCK3zs5/9LDtCNecl6ly5Hhcl35WJnLX1ztpeC/KuchwRT5C0spSyJSLOlPQUSf9QSnl0TwW6n0optkUTzXMdiXU+CTjPKXW1bNmy7AjVXvva12ZHaB3nu4LO86edj+ttNN6H9K6UdFJEHC3p45KWSfqUpN/qVzAMDucLEtfsziM8rtvc2de+9rXsCNVe//rXZ0eAEeeC7Dx/uo3GW5B3llK2R8RLJX2wlPK/ImJ5P4MBk8G1rLmOfEvecwRPOumk7AhVrrzyyuwI1aZP952l57qGs/MFuHNBXrVqVXYETMB4C/K2iHiVpNdJ+i+d7/ke1TAhriVT8r2V6Hwrznn+9Ite9KLsCFVcnzJ3N2vWrOwIVXp9uj+T6zFdkh599NHsCJiA8Rbk35f0ekl/WUq5JyKOlNToyhXI41yQXTmfBJxHv5/4xCdmR4AR1/n2FOQcnEu9jKsgl1LulvTmrq/vkfS+foUCJovrAcn5FqjzFIujjz46OwKMzJ8/PztClUceeSQ7QjXnu2vO00PaaKwPCvk3SZdK+tLuq1VExFGSfk/ST0spl/ct4R64XkG6vjmcl3lz5VrsJe8T2OzZs7MjtI7zHYe5c+dmR6jifHxxHjxwnbPeVmONIP+RpLdJ+mBErJH0kKRZko6Q9BNJf1tK+de+JgR64HoicM0teRfkadPGO+sMk+WOO+7IjlBtyZIl2RGqOB9fyN68ffbZJztClV6fhxlrHeSfS3qHpHdExBGSDpK0SdIPSym+k5jQGq4HJNfckndBds7u6vLLG70BOane+973ZkdoHdc7sZLv9LPPfvaz2RGqXHDBBT39/LiHS0opP5X0057+tpZjakjzmB7SPOdb5q7ZDz/88OwI1a666qrsCNU++MEPZkeo4nxcdD4fua7wc/rpp2dHqNLrFCjuJwIDyPViSvKeI+h6AnvFK16RHaHaBz7wgewI1VznlDqXTGfr16/PjlDF9WHUXt+fFGRgADlPsXDO7jrF4rjjjsuOACPOc+2dL8Bdp1i01YTfJRGxQNKhpZTb+5AHgLxLJiew5j3pSU/KjgAjrh9wIkkbNmzIjlDNdQpXW42rIEfEjZLO7rz+NkkPRcTXSilv62M2DAjn2/2utxKd5wg6l3vXgux6CxQ59ttvv+wI1datW5cdoZrr+ch1akiv3WW8I8jzSinrIuIPJf1dKeU9EcEIMoDHcC7IriM8ride5JgzZ052hGrOxxfXqS1r1qzJjlCl1+P5eP9rTYuIgyT9tqQ/7elvBDAm55OA8x0H1xFk523u7O67786OUMV1/WZJ+v73v58doZrrncEHHnggO0KVXp/HvdfGAAAgAElEQVQpGW9BvkjStZK+WUq5qfMpej/q6W8GMCU5l3vXEWTX1TfcffSjH82OUOU973lPdoRqriVT8j02bty4MTtClUamWJRSPifpc11fr5D0sp7+ZqABriNrrstHuXPdX1xX35CkM888MztCNdcPUHj/+9+fHaGac0F2vQB33ua9GO9Dek+U9BFJB5RSToiIp0g6u5Ty//Y13eNwvQpr606GiXPdxyXv7K4rcDiPIL/4xS/OjlDtxhtvzI5QZcaMGdkRqjmfR10vZF3nTve6r4z3//XHJP03SR+VpFLK7RHxKUkpBdmVa3FwPiABE+H6HnUdmZKk448/PjsCjDifj1zvULX1juZ4C/LsUsp3d9sxfY/ImBDX0iB5H0xdOe8vrtmdR5APO+yw7Agw4jqaKUlbtmzJjlDF9bjYq/HuaQ9HxBMkFUmKiJdL8nysETDgOtIgeY82uE6xcF19Q5Jmz56dHQFG9tlnn+wI1Vzfp64X4E2tg3y+pEslHRsR90u6R9Lv9PQ3Aw1wvfJ1Hvl23ebOnKdYABPhvIaz64ecjIyMZEeo0tQqFisknRURcyQNlVI8P1YFAPbCtdy7PvwjUe4zrF27NjtCtUWLFmVHqPbzn/88O0IV14/3bqQgR8R8Sa+VdIRGPzREklRKeXNPf3sl15OYK+fRTGAiXI8tzlNynMu9q9tv9/0gXOcPOXE9vrhODel1e493isU1kr4t6Q5JvkdiVHF9U0u+5d55mzuXNdf9xXkU1vX2rSSddNJJ2RGq/NVf/VV2hGqf+tSnsiNUcz02uh4XmyrIs0opb+vpbwIA9IXzBZXrA0CS9Nu//dvZEaq84x3vyI5Qbf78+dkRqrk+BDxz5szsCFWaWgf5kxHxR5K+IOmX65SUUtb09LcDfeZaHFyv2N257i+uuSXvKRZPfvKTsyMAfee6MlFTBXmrpPdL+lN1lnrr/O9RPf3tAKYc57LmynVkSvIeQT7qKE6BGD/XNZyHhoayI6QY73+tt0k6upTycD/DABjlXDKdR79d5wg67y/OI8jOa/Kiea4F2fW42Kvx/te6S9JAPElRSrE+GThyLjxAGzgfE50LcltH1lDH9YJq06ZN2RGqNPWQ3g5Jt0XEV/Wrc5BTlnlDs5xPvpT75rG/NM95mzsXZOfVQ9C8efPmZUeo4voBJ00V5H/p/AMAU5Zz0XTlfPvWudy7cp6zPnfu3OwIVVyfcWikIJdS/r6nvwVI4lp4XEcykcN1P5d8T76StHHjxuwIrfO5z30uO0K1U089NTsCJmCvBTkiPltK+e2IuEP/uXrFL5VSntK3ZACAKc95BNl1NPP888/PjlDtggsuyI5QbcWKFdkRqrR1rv1YI8hv6fzvi2t+eUS8QNKHJA1LuqyU8r7d/vz1ks7X6BznDZLOK6XcXfN3AVOJ84ggo9/Ncy6Zzvu660fwnnLKKdkRql1yySXZEaq5TrFwLch9XQe5lPJA51/fWEp5525/8cWS3vnYn/rlnw9LukTS8yStlHRTRCzbrQB/qpTyfzqvP1vSX0t6wYT/XwDAJHAta84XJc5TLLZs2TL2iwaQ60dku2tr0XQ13of0nqfHluEX7uF73U6R9ONSygpJiojPSDpH0i8Lciml+9HIOdrDNA4AAAaRa0Heb7/9siMAA2+sOchvkPRGSUdFxO1df7SvpG+O8buXSrqv6+uVkh4zQz0iztfoB5HMkPTcx8lxnqTzJGnBggVj/LUAMrmOwkq+IyXO29w5u+sUi5kzZ2ZHgBHX42KvxhpB/pSkL0r6K0nv6vr++lLKmjF+dk9bdE8P+l0i6ZKIeLWkP5P0uj285lJJl0rSYYcd5ns0BQD8CueC7LrMm/O0FjTP9YKq1yktY81BXitpraRXVfzulZIO7fr6EEmr9vL6z0j6SMXfAwAw5VyQXR+OdC32yOFakPv6kF6PbpJ0TEQcKel+SedKenX3CyLimFLKjzpfvkjSjzQOrgfUtt6mAICpyPWT9Fw/Otid64XJ4sWLsyNUmTatt4rbt4JcStkeEW+SdK1Gl3m7vJRyV0RcJOnmUsoySW+KiLMkbZP0iPYwvWIqodgDg831PYocrvuL6/rN7u64447sCFVcPyJ7eHi4p5/v5wiySinXSLpmt+9d2PXvb3nMDwEAWsN1moLkm915BPm8887LjlDtne/c28Jfg+uTn/xkdoQq06dP7+nn+1qQMTW4jpJIjH4D6B8KcvPOPffc7AjVnvvcPS7UNfD233//7AhVKMgAAFvOF+Cu2V3nwkrSUUcdlR2hdVw/4KRXFGRMaa4nMEa+gcHnenxxHkFetGhRdgS0BAUZAIAKrgXZ9QNOpPaOZqJ5FGQAMOda1CTvuyWu2935g0Jc533DDwUZGECuJ17Ju/Cgec77uivngjwyMpIdAS1BQQYAoEWcC7Lz9BB4oSADAFDB9Xa/a25J2rhxY3aE1vnFL36RHaFKr6u1UJABAGgR5xHkDRs2ZEeo9oY3vCE7QpULL7xw7BcNoFWrVvX08xRkAAAquM6fds0teU+xOOuss7IjVHnZy16WHSEFBRkAgBZxLsibN2/OjlDttNNOy46ACaAgAwAAC1u2bMmOUG3evHnZETABFGQAAFrEeQTZeYrF9OnTsyNgAijIAIA0zmXNdc1v523e68oEmbZv354dARNAQQYAoIJz0XTlvESdc7lvIwoyAACw4DwKyxrOXijIAAC0iPMorHN25zWc24iCDAAALDhPa9m0aVN2hCoHHHBAdoQqq1ev7unnKcjAAHJ9+AcA+sl5BHlkZCQ7QpX3vve92RGqXHTRRT39PAUZGEDOoySUe2CwOR9fnLO7LlF3+umnZ0eoMnfu3J5+noIMAAAsOBdk1ykWRx55ZHaEKjNnzuzp5ynIAADAgnNBdv0UQNcPOOn1biYFGQAAoM9YB9kLBRnApHIe4WH+NCbCeV935bzNXQuy68h3r/sKBRkAOlxPvhR7tIXzKhau2V0fLqQgA0DLuRZ7iXKP9tixY0d2hCrr16/PjlCl1+1NQQYApKHcYyKc9xfXEeR169ZlR6hCQQYAAOPmXDKdL0pcC/LmzZuzI1TpdXtTkAEAqOBaNJ1LpjPXgrxx48bsCFUoyMBecAJrnus2l7y3O9AGzscX1+yuH3DSyoLsupMB48U+noPt3jwuSprnvJ877y+u2911igWrWAAAbLmWBsm7rLly3l9cp1i4FuRWjiADAJDNtaxR7DER27dvz45QhRFkAACAAed6QeW6fnOvKMgAALSIa1Fz57rdGUEGAAAAurgW+15RkAEAALBHrlMsGEEGAAAYcK4jsa65e0VBBgAAwB65Lk/HCDIAAGgF59FM1+X1XAtyr4ayAwAAAACDhIIMAAAAdGGKBQAAQJ+5Tg9xzd0rRpABAACALowgAwPI9WEOYKLaOjqFOhwbm9fW96hlQeYNAgBTA8dzAIOIKRYAAABAF8sRZAAAACfcLfHCCDIAAADQhYIMAAAAdGGKBQB0cAsUACCZFmROYgAAAOgXplgAAAAAXSxHkAEAAJxw99sLBRkAkIbSAGAQUZABoIOyBgC/qq3HRQoyxtTWNwcAAGgnu4IcERQ2AAAqcQ4FxmZXkAEAGAQUTWDqoiADgDmKGoB+aevxxbIgu/7HKqVkR4AJ130cAPqJYyOaYlmQXfHGBgAAGHyWBdm1aDKCDABTh+u5CJiIoaF2fuiyZUF2/Y+1c+fO7AhVKPYA+oWSibZw7S5tZVeQnZd5c31zuBZ7iXKPiXE9tgBt4XoehR+7giz5vkFci6br9pZ8tzkA9IvzMR3Ncx046DW3ZUF2/Y/lelDasWNHdgQAGDiux3QAY7MryBGhadPsYkvyHc1kmkLzXPcVYKIomWgL18G9trJrmhGh4eHh7BhVXIum85vaOTswXpRMtAXH9Oa1dZvbFWTJ92TgOvLtWuwl3+yu+ziAwdfWwgNMhGVjc31zu5Ye1xF7yXeqgmuxl7yzAwAgGRZk52XeXDkXZNcHDF2LvURBRntwLsJEsL94sSvIpRTbE7Drm8N15FvyLfeuxV7y3c8l33LvvM2BiWBfb15bt7ldQZZ8R9dcy5prbol9JYNryQQAYBe7glxK0fbt27NjVHG9CnMeQXbN7pobwOBzPRcBTeprQY6IF0j6kKRhSZeVUt73OK97uaTPSTq5lHLz3n5nKUVbtmyZ9KxNcD0oOZc115FY19yS9/7iPLUFzXM9piMH+4uXvhXkiBiWdImk50laKemmiFhWSrl7t9ftK+nNkr4znt9bStG2bdsmO24jXEvP9OnTsyNUcz0gOZdM5+yuU3Kcp7W4vkcBTG39HEE+RdKPSykrJCkiPiPpHEl37/a6v5D0PyW9fTy/dOfOndq8efNk5myM64nAuSC7XpS4rpktyfYCVvJ9jwJtwXu0eW3d5v08Cy+VdF/X1yslndr9goh4mqRDSylfiIjHLcgRcZ6k8yRp/vz52rp1ax/i9p/rTuZckF2zOxdk5+yuzzc4jyADwCDq55lsT23wl0fxiBiS9DeSfm+sX1RKuVTSpZK0dOnS4jpC5Xrr2fWCRPK9Ze56MSX57ueSb3bX/RwABlU/C/JKSYd2fX2IpFVdX+8r6QRJN3bKwIGSlkXE2Xt7UK+UYvsgjWvRdJ3SIvlmnzFjRnaEaq7TWiTfgux8QcXod/Oc9xegKf0syDdJOiYijpR0v6RzJb161x+WUtZKWrzr64i4UdLbx1rFwplrsXcdsZd8b5m77iuS9xQL13Lv/B51RtHERLC/eOnbmayUsj0i3iTpWo0u83Z5KeWuiLhI0s2llGU1vzcibEd5XG+Duo58S9LIyEh2hCpz5szJjlDNdd63xAgyAGBUX4d6SinXSLpmt+9d+DivPXM8vzMibG8/uxZN59FM1zWzXaeGSNKsWbOyI1RzHf12HfmWfO/yAJja7M4GEWE7QuVaNF1HviXfW8+uF1OS95xSRpCBwca+jqbYFeShoSHts88+2TGquBYH12Iv+Y5Obdq0KTtCNedyP3v27OwIVVxHviXf96jke0wHJqKtAwd2R9WIsL2F63owdS7IrkXTdWqI5D09xLUgO0+xcD35St7HRgB7Z1mQmWLRLNftLfmWNdepIZL3CLLrdCLngswtcwCDyLIgu95OnDlzZnaEKs5lzXUk1rXYS76j9pLvRaxzQXYeQXbFRQkwNsum6frmdi32rquGSL6j384F2fmCyjW7634ueRdk13MRgLFZNjbXubyuJwLXkW/Jd8kx1/WbJe9y7zr67bxutvPoN4Cpy7IgM0+wWc6jU64rnjhflDiXe9eC7DpoIPneWZN8Bz1cz6FAk+yOTKUU23mCrrfjXE8Ckm+5d80teZc11znrzkulOR9fXI/pAMZmV5Al36tf12LvfAJzLZquI9+StH79+uwI1VynhzivHOI8PcT12Oh6LgKaZFeQSym2Bdl1tMF1e0u+JzDnByOdb5m7jsQ6F+S5c+dmR6jmOm3O9WFUoEmWZzLXW7iuV+2uxV7yLZquDxdK3tnXrl2bHaGK87zvhQsXZkeo5lqQAYzNsiA7FzZHziPIrtmdT7zO00PWrVuXHaGKc0F2HbWXfO+WuN5Zk3yP6fBj9+6OCNs3t+sb2zW35Dtq71yQnVfgcD22uD5cKPm+RyXfZxwYZALGZleQnbmefJ0Lsmt212lEkvcUC9fsziPIrkvrSb4Xg67nIsn7ggpeLAuy69Wvc+lx5Xowdc0teY9+u04PcV45xLncL1iwIDtCFef3KA8YNq+tA02WBdm1aLoWe+fRBtc3tnNBdp2XKfk+1Ol6bJF8l9Zz5nxMB5rieyYz5Frs0TzXYi957+euI8jOFyXOUyxcL2Sd9xfncu98XG8ju3eJ8zrIzm9sV67b3HUfl7yzuz505Vx4nNdwdr3d7zzFAmiK5VHVdYTKtTi4lkzkcB1Vk3wL8uzZs7MjVFu9enV2hGquo9/z58/PjlDN+Xzk2gHayrIgu6LYN8/5YOrKuSC7zkF2/rjmhx9+ODtCtQ0bNmRHqMKHs+RwXvO7jSjIGJNrsXfmXOydL6hcs7sWe8l7X3def9qV8wOp8GJXkCPC9grS9eTrXJBdt7lzaXDmur84F2Tn7Bs3bsyOUMX5Lo/r+R9+7AqyRHlomvPBFJgI133duWQ6Z3ddw9n5Vr/zA6munAfJemG3p0WE7QHV9Yln11E1yfeN7bzNnbHdm+f8gOGjjz6aHaGK89rTc+fOzY5QzXVwz3XgoFd2BXloaMh2rVLmTjXP9aLEmetJQPItyK4XgpJ3QXbd7s6fvLjffvtlR6hGB/BiWZBnzpyZHaOK68HUtTRIvle+rrmRw/k96rq0nuQ7H9Z5BNm5ZDoPHrSRZUF2vcXC3KnmuRYH19ySd3ZXznNKXafMSbIdrHF9uNCd6wVVW9k1tuHhYdtbLBTk5rkWB9fcEqPfGZy3ufOIoOv0ENe505L3sdH1bonrB+L0yq6xOc9Bdp1i4XxAch3hcV5f1bmsoXmux0XJ9/ji/GyGc3amWHixLMj77rtvdowqjCA3z7WsORdk5xOY6/QQ55Lpus0l3xFk5/3FeXqI6ycYOt/l6YVdYxsaGtKsWbOyY7SK8wiy68Mozvu469qwkndZc+V6ESv53jJ3LjzOx5dFixZlR6jiOvLd635uV5Ajwvag5LqTOXMt984jyM4jPK77i3Phcb4ocZ1i4ZpbkjZs2JAdoZrrQ3rOx5de2BVkZ65TLFxvI0q+hcd15FvyPoG5Xpi47ueS9wiya3FwPqY7H19cLwZdu0vrRpBLKbYHVNedzPWqV/KdquD6IKrknd31PepckJ25zuV1HkF2XoHDtbu09e633dlgx44dWrt2bXaMKq4PFzqvU+pakF2XMpS8R3geeeSR7AhVtm7dmh2hmmtpkHwvTObMmZMdoZrzQ8CudwZdp7W2bgR5+/btWr16dXaMKq634+bNm5cdoZrriKBrsZe8b9+6jn47X5Q4F2TXW+bOgx7OXAuyawdoXUEeGRnRrbfemh2jyrHHHpsdocrhhx+eHaHa/PnzsyNUcb4F6nqnRPIduXe+7ew6Civ5lnvnC3DXQQ9JWr9+fXaEKm1dns5uT9uwYYO+8Y1vZMeo4noScx35lnxHSlw/Tl3yLsiuIyXOo/auD0ZKvgXZ+ZjuepdH8i3Izs8h9cKuIG/atEm33357dowqrjvZkiVLsiNU23///bMjVHEuPK7z1STf7e6aW5LtMyWSb0F25nx3ra0f2ezKriCPjIzolltuyY5R5eGHH86OUMV5tMF1YXbXW/2S9wiya3bn/cX1mRLJ94Ex52ktznfXXPd11/2l11Vm7Aqys5/97GfZEapce+212RGqLViwIDtClYMPPjg7QrVjjjkmO0I119u3rlNDJO/Rb9eHrlwfLpS89xfXouk6DYqCjL679957syNUu/vuu7MjVHG9mJK8H+p0LcjOpcF1m0u+68O6FjXJewqXq7ZeCFKQMaW5jn4fffTR2RGqPfGJT8yOUO3II4/MjlDFeYqFc3bXaXOuU0Mk71UsXMv9unXrsiNU6fUZAd89DRgH1xPBnXfemR2h2v33358dodqhhx6aHaGK64lX8h79dt7urpyfiXFdVWlkZCQ7QpVWjiC73tZynveFZl1//fXZEaqdeuqp2RGqua5V7rzSjOta5ZLvA2POHyzjfB513V9cR5BbV5D33XdfnXTSSdkxqixfvjw7QhXX9ZuR47rrrsuOUO35z39+doQqixcvzo5QzXmKhWv2hx56KDtCNeel9VwLsuv+0ropFrNmzdJxxx2XHaOK64MRX//617MjwMhtt92WHaHa97///ewIVZ72tKdlR6jmPE3BdXqI8zZ3fWBM8l3D2bW79MquIM+ePVtPf/rTs2NUcV1j1fmK/Vvf+lZ2hNbZunVrdoRqX/3qV7MjVDnrrLOyI1Q75JBDsiNUc50e4lrsJd/5sJLvR3y7Tmtp3TJvM2bM0NKlS7NjVHEtmq6Lm0vSXXfdlR2hivOnizn77ne/mx2hygMPPJAdoZpzQXYtmq65Jd9nkCTfT9N13ea9PtBpV5CHh4dtPx3Ndb1P16tHyXf+9BVXXJEdoZVWrFiRHaHKv//7v2dHqPaUpzwlO0K1OXPmZEeo4no3U/ItmZLvChyuI9+tK8gRYfsGcT0ouT7ZL/k+feu8zJvrPF5nX/7yl7MjVHvlK1+ZHaGa67KArlNDJN+yJvlOP3N9uLDXrmhXkEsptmvbuj4Y4Xw7zvVT3U488cTsCNUoyM27+eabsyNUc54ecthhh2VHqOJ8THd90E2SNm3alB2hiuvd79aNIO/cuVMbN27MjlHFdSdzPiC5jvCcfPLJ2RGqOd/uX7VqVXaEKs7r2ro+GCnJ9oHxefPmZUeo5rq0nuQ75c91KlHrRpCdC7Lr/CPXKS2Sb7l3HfmWvJcccy3IzpynE7mOCDqPILsONDlzPY+2bgR5x44dWr9+fXaMKq5TQ5y5zvt2XalFkk477bTsCNWuvvrq7Ait86//+q/ZEaq9613vyo5Q5clPfnJ2hGrO86ddP2radXCvdQV5586dtguFuy7z5sx19Nv1oQjJe1WCM888MztClRtvvDE7QjXX47kkfeMb38iOUOWEE07IjlDNuSC7jty7TuFqXUGWel/8OYvrCLLzCcx1UXnn24gLFy7MjlDt+OOPz45QxbkgO7v99tuzI1RxPRdJ3sdGCrIXu4LsvMyb6wiyc0F2/cAN55OA8+i360NXrideyfciVpKuuuqq7AhV/vzP/zw7QrUDDjggO0I11wcMH3zwwewIKewKsuQ7H8a12Dt/UIhruXct9pL3OqXHHHNMdoQqL3zhC7MjVLvyyiuzI1R7+OGHsyNUcX4w0vmTFxcsWJAdoYrrJ+n1yrIgu3It9s5cR+2dR9W2bNmSHaGaa7k/+uijsyPAyG233ZYdodrzn//87AjVXJdLc/0Mh9bNQS6l2JYe1xFk19ySb0F2HfmWfD+9UPKdquC8tJ7zdCLXZd6WLVuWHaHam9/85uwI1VynWLhOm2vdOsiS9y1/R84F2fXWkPM+7loaJGn79u3ZEao4z8t86lOfmh2h2re//e3sCFW++93vZkeo5jz9bPHixdkRqrgOHPR6/rcryM4fNe1cNF2xzZvnPD3E9Wlt56WvXJfWk3wLsvMF+F133ZUdodrznve87AhVXI8vrRtBdp5i4bo8HZrnXOxd35+S7+iU661byXflEOS45pprsiNUO+uss7IjVHFdunPatN4qrl1BlnyLputVu3PhcZ1i4VyQnbk+YLh169bsCNX233//7AjVjjrqqOwIVVasWJEdoZrrqL3k+z51fbiwdVMsJN/VIFyLvXNBdt1XXIu95F3uXU9g69evz45QbcmSJdkRqrkur3fJJZdkR6jmXJBXr16dHaGK69zpVhZkV65lzbXYS77l3vVug+RdkF05Pxjp7AlPeEJ2BBi5//77syNUOeigg7IjVGndMm+Sb9FE81zLveuDqJL36LdrdueC7Lyk4QknnJAdAUauvfba7AhVTj755OwIKewKckRQkDFuriOxriPfknd219Fv1+XpJN+VQyRp0aJF2RGqvOhFL8qOUO3qq6/OjlDtiiuuyI5Q5d3vfnd2hBR2BbmUYjsqSLFvHtu8ec5lzfUTo5w5LwvoOn/6uOOOy45Qzbkg/8d//Ed2hCqud6h6HSCzK8iS721zV5TM5jnv466j9pJvuXcu9q4nX8n3oc5nPOMZ2RFg5Gc/+1l2hCq9rkpEQcaURrlvnvP7k4LcPNdtLkkbN27MjlDlsMMOy45Qbd68edkRqrmus37ddddlR6iybt26nn7esiC7ci4OQBu4jn4z7zuH6/QQ109Gk7wfGLv++uuzI1S5+eabsyNU6fX92deCHBEvkPQhScOSLiulvG+3P/89Se+XtGvtk78tpVzWz0yYOIp981xXU3BHQW4eBbl5zsf0M844IztCNdeC/E//9E/ZEVL0rSBHxLCkSyQ9T9JKSTdFxLJSyt27vfSfSylv6lcOtJvricA1N3I4F2Rnrssxus6dlnw/vRB++jmCfIqkH5dSVkhSRHxG0jmSdi/IE+Y6r5TSg/Fy3ldc35/OnAuyc3bXfd15ab3DDz88OwJaop8Feamk+7q+Xinp1D287mURcYakH0q6oJRy3x5eMyW4Hkydy5pzdjTP9T3qOjVE8n5Iz3V6iPNHk7uuPS1Jz33uc7MjVLnhhhuyI6ToZ0He05lm97byb5I+XUrZEhGvl/T3kh6zB0XEeZLOk6QFCxZMdk5MYa6Fh2Kfw3V/ceY8guz6PnWeYuF8MfjkJz85O0IVCvLkWynp0K6vD5G0qvsFpZTVXV9+TNLFe/pFpZRLJV0qSYcddlhxPSgBbeD8/nTN7lzsnQuya3bXudOS97rZzitwtFE/C/JNko6JiCM1ukrFuZJe3f2CiDiolPJA58uzJX1vPL/Y9STmyvnki+Y5j/C47uuuuSXv/cW1ILtODZG8C/KBBx6YHQET0LeCXErZHhFvknStRpd5u7yUcldEXCTp5lLKMklvjoizJW2XtEbS743j91ofUB2x5BgmwvkClsLTPOfjuetIrPP+4vphG5J0xBFHZEeo4rpudq9z7fu6DnIp5RpJ1+z2vQu7/v3dkt7dzwzonWtpkLxH1tA857KG5rkW5BkzZmRHqLZ58+bsCNVcjy9nn312doQqV199dU8/zyfpNch5ZM0V2xxt4HrilbzvULlud+eVQ5z3ly1btmRHqHLMMcdkR6gyc+bMnn7esiBTeprF9kZbuO7rrrnduW5315FviekhGU444YTsCFX22Wefnn7esiC7cj2YOmObow2c93PXUVhnzgV5+vTp2RGqrVu3LjtCFdcPZ5k2rbeKa1mQnU8GwHiwj6MtnPd112ccnLS4hD0AAAx5SURBVC9KnLO7rj/tPK2lFxRkAJOK92fznLe5c3bX4uC8zZ1Hv133F9di3+t+blmQAQDI5lo0XXNL3qsquW73kZGR7AhVer3bQEHGmFzf1MBEsa9jIthfmudckF2nh7jOne51X6EgAwBQwbUgu86dlnxLpuSbnRFk4HE4H0xdT2DO2F8wEewvaAvX0W/XYt/KOcjOB1RHzicB9pXmsb80z3mbO3PdX5w5b3PXguyau1cU5Aa5nsRct7fku82dOe8vQBs4Hxedjy/O272NLAuyK+c3tivXbc6BFG3hvK+7Hl+csc2b5/we7YVlQeYNAqAfXI8tzicw5+yuXOeUSr7vUWeu+0vr5iBHhO1i20AbuB5MgbZwPoc6X1C5lnvnbd4Lu4Is+e5kQBsMDw9nR6jW1hNBJucLKtf9xfkc6rrNJd/t7rzNe0FBBgYQ+3iOtp4IMjmPZjqXe1fO29x1X2/r+YiCDAwg9vEcrtudYp/DtfAgh+vxxTV3r+wKMnOQgcHW1oMp2sf1wsQ1t8TxJYPrNu81t11Blnz/YwFtwAUsMNicP/jB+fzvemHS1mO6ZUFu63+sLM4HJDSP/aV5ztvctTQ4c95fnLnu666dq3UjyEyxwERwIsBEuO4vHBMxEc4PujlzLciux8Ve2RVkqb3/sTBx7CvNo6w1j/0cE8FSjGiD1o0gS5wMmkbhaZ7zPu6c3ZXzNuf4golg9Lt5bb0osSvITLFonvPJ15XzNnfO7npscd7maJ7z/uKc3ZVrQW7lCLLrSQzNcz2YOu/jrttc8s7uynmbO2dH81yLZltZFmQOSs1yLmuunPdxsjfPNbfkfXxx3u7AVNe6EWSmWDSPk0DznLc52ZvnfEx03eaSb3bX3O5cR5Bd95fWFWTJ9z+WK+ft7ZzdlfM2d32633mbO2d3vTBx3uZAUyjIGJPz9nbN7nrilXy3ueSb3TW3xL6ewTW35J0dzWIEGX3H9sZEOO8vrtkpmTlcs7vmlryzO79P28iuIEeE9RvEEW/q5jnv487ZXfd1521O9ua55pZ836OS93ZvI7uC7Mz5je3K9YDkmhs5nPcXsjfP+Vzkus0l7+yOWjnFwvnN7cj5Te2a3TW35P3+dN7uaJ7rvu68n5O9ea77ea8sC7LrTubKeXs7Z0fzXPcX19ySd3ZgIlyLZlvfo5YFGZjqXA+kUnsPpgD6j+MLmhJuC1dHxEOSfpadY0AslvRwdogphO05ediWk4vtOXnYlpOL7Tm52J6T59dKKfvW/rDdCHIpZUl2hkERETeXUk7KzjFVsD0nD9tycrE9Jw/bcnKxPScX23PyRMTNvfy8731cAAAAoA8oyAAAAEAXCrK3S7MDTDFsz8nDtpxcbM/Jw7acXGzPycX2nDw9bUu7h/QAAACAfmIEGQAAAOhCQR5QEXFoRHw1Ir4XEXdFxFs6318YEddFxI86/7ug8/2IiA9HxI8j4vaIeHru/4PBEhGzIuK7EfEfne35PzrfPzIivtPZnv8cETM635/Z+frHnT8/IjP/IIqI4YhYHhFf6HzNtqwUET+NiDsi4rZdT17zXq8XEfMj4oqI+H7nGPpMtufERcSvdfbJXf+si4i3si3rRcQFnXPQnRHx6c65iWNnhYh4S2c73hURb+18b9L2TQry4Nou6b+WUo6TdJqk8yPieEnvkvSVUsoxkr7S+VqSXijpmM4/50n6SPORB9oWSc8tpTxV0omSXhARp0m6WNLfdLbnI5L+oPP6P5D0SCnlaEl/03kdftVbJH2v62u2ZW+eU0o5sWuJJ97r9T4k6UullGMlPVWj+ynbc4JKKT/o7JMnSnqGpBFJV4ltWSUilkp6s6STSiknSBqWdK44dk5YRJwg6Y8knaLR9/iLI+IYTea+WUrhH4N/JP2rpOdJ+oGkgzrfO0jSDzr//lFJr+p6/S9fxz+P2ZazJd0q6VSNLsg+rfP9Z0q6tvPv10p6Zuffp3VeF9nZB+UfSYd0Dj7PlfQFScG27Gl7/lTS4t2+x3u9blvuJ+me3fcxtmfP2/X5kr7JtuxpGy6VdJ+khZ1j4Rck/V8cO6u25SskXdb19Z9Lesdk7puMIBvo3FZ5mqTvSDqglPKAJHX+d//Oy3a98XZZ2fkeOjpTAm6T9KCk6yT9RNKjpZTtnZd0b7Nfbs/On6+VtKjZxAPtgxo9GO3sfL1IbMteFElfjohbIuK8zvd4r9c5StJDkv6uMwXosoiYI7Znr86V9OnOv7MtK5RS7pf0AUn3SnpAo8fCW8Sxs8adks6IiEURMVvSb0k6VJO4b1KQB1xEzJV0paS3llLW7e2le/geS5R0KaXsKKO3Cg/R6G2Z4/b0ss7/sj0fR0S8WNKDpZRbur+9h5eyLcfvWaWUp2v0NuD5EXHGXl7L9ty7aZKeLukjpZSnSdqo/7zNuidszzF05sSeLelzY710D99jW3Z05sOeI+lISQdLmqPR9/zuOHaOoZTyPY1OOblO0pck/YdGp6Y+nglvSwryAIuI6Rotx/9USvl859u/iIiDOn9+kEZHQ6XRq6FDu378EEmrmsrqpJTyqKQbNTq3e35E7PrI9e5t9svt2fnzeZLWNJt0YD1L0tkR8VNJn9HoNIsPim1ZrZSyqvO/D2p0jucp4r1ea6WklaWU73S+vkKjhZntWe+Fkm4tpfyi8zXbss5Zku4ppTxUStkm6fOSfl0cO6uUUj5eSnl6KeUMjW6XH2kS900K8oCKiJD0cUnfK6X8ddcfLZP0us6/v06jc5N3ff+1nSc1T5O0dtdtBkgRsSQi5nf+fR+NHqi+J+mrkl7eednu23PXdn65pBtKZ+JS25VS3l1KOaSUcoRGb7veUEp5jdiWVSJiTkTsu+vfNTrX807xXq9SSvm5pPsi4tc63/pNSXeL7dmLV+k/p1dIbMta90o6LSJmd87xu/ZNjp0VImL/zv8eJun/1ug+Omn7Jh8UMqAi4nRJX5d0h/5znuefaHQe8mclHabRN9srSilrOm+2v5X0Ao0+afz7pZSbGw8+oCLiKZL+XqNPDQ9J+mwp5aKIOEqjo6ALJS2X9DullC0RMUvSJzU693uNpHNLKSty0g+uiDhT0ttLKS9mW9bpbLerOl9Ok/SpUspfRsQi8V6vEhEnSrpM0gxJKyT9vjrve7E9J6Qzv/M+SUeVUtZ2vse+WSlGlxh9pUanAyyX9IcanQvLsXOCIuLrGp2TvU3S20opX5nMfZOCDAAAAHRhigUAAADQhYIMAAAAdKEgAwAAAF0oyAAAAEAXCjIAAADQhYIMAAAAdKEgA0CiiHhzRHwvIv4pImZGxPURcVtEvHKCv+fMiPj1vfz5SyLiwgn+zus7H48LAK0ybeyXAAD66I2SXlhKuafzCU/TSyknVvyeMyVtkPStx/nzd0g6e4K/85OdfH9ZkQcAbDGCDAANiIi3RcSdnX/e2vne/5F0lKRlEfFOSf8o6cTOCPITIuJ9EXF3RNweER/o/MySiLgyIm7q/POsiDhC0uslXdD52d/Y7e9+oqQtpZSHO19/IiI+EhFfjYgVEfHsiLi8M5L9ia4fXabRjxkGgFZhBBkA+iwinqHRjzs+VVJI+k5EfK2U8vqIeIGk55RSHo6I7+g/P7p7oaSXSjq2lFIiYn7n131I0t+UUr4REYdJuraUclynbG8opXxgDxGeJenW3b63QNJzNTqq/G+d1/yhpJsi4sRSym2llEc60z4WlVJWT+pGAYABRkEGgP47XdJVpZSNkhQRn5f0G5KW7+Vn1knaLOmyiLha0hc63z9L0vERset1+0XEvmP8/QdJemi37/1bp3jfIekXpZQ7OtnuknSEpNs6r3tQ0sGSKMgAWoOCDAD9F2O/5FeVUrZHxCmSflPSuZLepNER3yFJzyylbPqVvyD2+ldskjRvt+9t6fzvzq5/3/V197lhVufnAaA1mIMMAP3375JeEhGzI2KORqdOfH1vPxARcyXNK6VcI+mtknY9uPdljZblXa/b9f31kh5vJPl7ko6eaOgYbd0HSvrpRH8WAJxRkAGgz0opt0r6hKTvSvqOpMtKKXubXiGNlt0vRMTtkr4m6YLO998s6aTOg3t3a/ThPGl0HvFL9/SQnkYL+tNijGHmPXiGpG+XUrZP8OcAwFqUUrIzAAD6LCI+pNF5x9dP8GeWlVK+0r9kADB4GEEGgHZ4r6TZE/yZOynHANqIEWQAAACgCyPIAAAAQBcKMgAAANCFggwAAAB0oSADAAAAXSjIAAAAQJf/HxyhT/1VfqHxAAAAAElFTkSuQmCC\n",
"text/plain": [
"