{ "cells": [ { "cell_type": "code", "execution_count": 12, "metadata": { "code_folding": [ 0 ] }, "outputs": [], "source": [ "# Configuration step (Please run it before the code!)\n", "import numpy as np\n", "import sys, obspy, os\n", "import matplotlib.pyplot as plt\n", "\n", "from obspy.core import UTCDateTime\n", "from obspy.clients.fdsn import Client\n", "from obspy.geodetics import gps2dist_azimuth as gps2DistAzimuth # depends on obspy version; this is for v1.1.0\n", "#from obspy.core.util import gps2DistAzimuth\n", "\n", "#from PIL import Image\n", "import requests\n", "from io import BytesIO\n", "\n", "%matplotlib inline\n", "import warnings\n", "warnings.filterwarnings('ignore')" ] }, { "cell_type": "code", "execution_count": 13, "metadata": {}, "outputs": [], "source": [ "# Functions\n", "# collection of functions used in noise correlation processing\n", "\n", "def normalize(tr, clip_factor=6, clip_weight=10, norm_win=None, norm_method=\"1bit\"): \n", " \n", " if norm_method == 'clipping':\n", " lim = clip_factor * np.std(tr.data)\n", " tr.data[tr.data > lim] = lim\n", " tr.data[tr.data < -lim] = -lim\n", "\n", " elif norm_method == \"clipping_iter\":\n", " lim = clip_factor * np.std(np.abs(tr.data))\n", " \n", " # as long as still values left above the waterlevel, clip_weight\n", " while tr.data[np.abs(tr.data) > lim] != []:\n", " tr.data[tr.data > lim] /= clip_weight\n", " tr.data[tr.data < -lim] /= clip_weight\n", "\n", " elif norm_method == 'ramn':\n", " lwin = tr.stats.sampling_rate * norm_win\n", " st = 0 # starting point\n", " N = lwin # ending point\n", "\n", " while N < tr.stats.npts:\n", " win = tr.data[st:N]\n", "\n", " w = np.mean(np.abs(win)) / (2. * lwin + 1)\n", " \n", " # weight center of window\n", " tr.data[st + lwin / 2] /= w\n", "\n", " # shift window\n", " st += 1\n", " N += 1\n", "\n", " # taper edges\n", " taper = get_window(tr.stats.npts)\n", " tr.data *= taper\n", "\n", " elif norm_method == \"1bit\":\n", " tr.data = np.sign(tr.data)\n", " tr.data = np.float32(tr.data)\n", "\n", " return tr\n", "\n", "\n", "def get_window(N, alpha=0.2):\n", "\n", " window = np.ones(N)\n", " x = np.linspace(-1., 1., N)\n", " ind1 = (abs(x) > 1 - alpha) * (x < 0)\n", " ind2 = (abs(x) > 1 - alpha) * (x > 0)\n", " window[ind1] = 0.5 * (1 - np.cos(np.pi * (x[ind1] + 1) / alpha))\n", " window[ind2] = 0.5 * (1 - np.cos(np.pi * (x[ind2] - 1) / alpha))\n", " return window\n", "\n", "\n", "def whiten(tr, freqmin, freqmax):\n", " \n", " nsamp = tr.stats.sampling_rate\n", " \n", " n = len(tr.data)\n", " if n == 1:\n", " return tr\n", " else: \n", " frange = float(freqmax) - float(freqmin)\n", " print(frange)\n", " nsmo = int(np.fix(min(0.01, 0.5 * (frange)) * float(n) / nsamp))\n", " print(nsmo)\n", " f = np.arange(n) * nsamp / (n - 1.)\n", " print(f)\n", " JJ = ((f > float(freqmin)) & (f" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# after detrend and cosine taper\n", "stp2[0].plot()" ] }, { "cell_type": "code", "execution_count": 24, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "YT.ST01..BHZ | 2010-09-13T07:26:34.000000Z - 2010-09-13T07:26:59.575000Z | 40.0 Hz, 1024 samples\n", "YT.STT1..BHZ | 2010-09-13T07:26:34.000000Z - 2010-09-13T07:26:59.575000Z | 40.0 Hz, 1024 samples\n", "done!\n" ] } ], "source": [ "# 1 bit normalzation \n", "# then whitening; devonovolution bythe mean spetrca \n", "\n", "stw = stp2.copy() # copy stream\n", "n = 1\n", "for tr in stw:\n", " tr = normalize(tr, norm_method=\"1bit\")\n", " #tr = whiten(tr, 0.1, 0.2)\n", " tr = whiten2(tr, 1, 5, 3)\n", "\n", " print(tr)\n", " #print(n)\n", " n += 1\n", " \n", "\n", "print ('done!')\n" ] }, { "cell_type": "code", "execution_count": 25, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAvQAAADtCAYAAADQmuFsAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAAPYQAAD2EBqD+naQAAIABJREFUeJzsvXl8VNX9//862ZeZCUlIAiTsmwYQtSJawEc/omjrRrW1Kh+11f5ara2t9eOKWv24a/l8/bTVVmurVcTla+tP1CIWtyKKKwgCspMQIAlkm0yWSTJzvn/ceV/O3Ll35s7GzCTv5+Mxj0zueu6dc899ndd5n3OElBIMwzAMwzAMw2QmWalOAMMwDMMwDMMwscOCnmEYhmEYhmEyGBb0DMMwDMMwDJPBsKBnGIZhGIZhmAyGBT3DMAzDMAzDZDAs6BmGYRiGYRgmg2FBzzAMwzAMwzAZDAt6hmEYhmEYhslgWNAzDMMwDMMwTAbDgp5hGIZhGIZhMhgW9AyT5ggh8oUQfxVC1Ash3EKItUKIkw3b3CyEOCiEaBVCPCSEEMq6PwkhdgghpBDiW4b9CoUQS4UQnYHjXxwhLbOEEBuEEN1CiPeFEGOVdaOFEK8LIdqFELuFEBdGOFa4dN0lhNgbSNd2IcQVgeW3CiE8gY9XCNGv/P+nwDbfDhy3SwjxqhCiVDmux/DxCyGut0jfD4UQXwTueZ0Q4mbDeqcQ4snAPW8XQjwT5lpNr8dkuxVCiIFw9y0ehBBLhBA7A+nYIIQ427D+h0KIhsA1PyWEyFPW3SmE2BS4Zz807JclhHgkcB+ahBDXRUjHRCHEmkA++kIIMVNZVyqEeE4IcUgIsV8IcW2EY4VL11VCiF2B66kXQtwSWL5IyQO9Qgif8v+KwDbh8vomQz7yCSF+b5G+s4QQHwohOgLX8z9CiFxlfY4Q4uHAfesQQrwb4VrVPO8xud7dgd/3ZSHEsDDHkoFnhI51q53rE0J8K3Cv1fXzrM7DMMwRQkrJH/7wJ40/AIoB3AFgDLRK+EUADgFwBNZ/B8BeABMBjACwEcCVyv5XAfgWgJ0AvmU49kMA3gTgAnASgDYAUy3SkR84z48BFAC4F8BqZf37AH4LIBfAcQAOAjgqzHWFS9dkAK7A9ykADgCYYdjmZgBPG5ZVAmgP3JMiAH8D8KzF+UcBGAAwPkz6vhm4nnEAvgZwibL+/wJYErh3uQCOC3Otdq5nIYAPAAwkMS/dGTh/FoD/CNyr8YF1MwK//ywAJQBWAbhb2fc/AZwB4N8Afmg47s8ArA/c/8kA9gGYHyYdnwC4K5CPrgawC0BOYN3fADwPoBDAhMC608IcK1y6xgMoD3wfCWATgLMM21wE4L1o8rrJtm3GPKysvxjAgsD1VABYDeBWZf3DgesdDiAbwDci/H5PWqz7DwCNAI4OpOlxAEvDHEsCqLGRZ4KuD9ozuyNZeZQ//OFPbJ+UJ4A//OFP9B8A++nFHxADtynrfgjgfZN9vjaKjoCwnKv8/zSAuyzOeYb6IocmmLsDoskBwI+AaA2s/zOAe2xcS0i6DOsnB9J5rmG5maD/KYBVyv/jAPQCKDQ57n9ZiTSLdNwH4PeB77UBoZkdw28Xcj0B0fgVgLlIoqA3ScuHAC4IfL9fFYsB4VZnss+bCBXOHwH4T+X/OwH8zeKcUwF4AOQry/YA+I/A90MAjlHWLQ4nTMOly7B+JLTK7rWG5WaC3jKvmxz3ewDqAAib9/wnAF4LfC+DVvF12dz3TlgL+t8C+B/l/2oAXgBFFtvbFfRB1wcW9PzhT1p+OOSGYTIMIcRkaEJgR2BRLYANyiYbAUyzcZxSaI6+3X2DziOl7Ibmrk8DIJSPfgo76QiTvpuFEF0AtkFzfFfZ2M2Yxj0A+qG1Xhi5FIAeJiOEmCuEaA9z7FOgObyA5mJvB/A3IUSL0MKgvmlI++tRXM/NAF4A0BDxChNE4PefDmBzYJFZPhojhHDYOFw0ebAWwDYppTfM9onMR5cIITqhVYKLoLWsRCJcXjdyKbQKhwycb0wg9GiMxbHVfDQDmqt+VyDEaIMQ4lxD2jcY9v9eIM+tE0Kcb7xcw/c8aBVIKz4RQuwTQjwthCi32Cbo+gKMFkI0Cy187A4hRHaYczAMcwRgQc8wGYQQohDAUgD3Syk7AosdANzKZu7AskjQNp029zWeR99eStkJYA2AO4UW838CNGev2EY6TJFSPhA450kA/g6gz8ZulmlUFwghjoEWeqKLOynlB1JK05hjIcSvoVWi/hZYVA0tjOJdaJWiRwAsp5hlKeUDUsqg+HSr6xFCjANwITSH9YgghMgC8BSAv0sptwQWm+UjWh6JaPJgpN/oTQCLhRCOQOX1SsSXj5ZJKZ3QxPPTCM7vVtjNR+UAvg3gWeV89VLKYVLKeuNBhRAXAJgP4H8Ci6qhVao6oIWA/QLAUiHERCXtxyiHeAnAUQCqANwE4CkhxImBdW8CuEQIMT1QTtwJzYW3unenABgL4Fho4UBPm6Q35PqgtajNhJbvz4OWd39lcQ6GYY4QLOgZJkMIdKT7v9Cc+f9WVnmgxXETrsCySNA2TrN9DR3jxpicx3iuRdDcwAYAj0JzvxssjmULqfExNLHzE5vXFC6NxKXQwh7COfIAtA6UAK6DFnvdE1jcA2CPlPIvUsp+KeUL0MJoZoc7lsX1/B8At0speyOlxSRtpp0jbfAYtDj5q5RlZvmIlkfCMg8KraOv2nky0m90LbRWlR0A/n8Az+FwPjIeyzZSyq+ghc3cEcP1GNNIXARgvZTy60gHFEL8B4A/AjhHStkcWNwD7VrvkVL2SSnfh9YX5TSLa9gspWyUUg5IKd8CsAyaqIaUchW0cuEVaCFMO6FVXkxbfaSUqwN59yC0e/4dIURBpOsLnP9rKaVfSrkZwN0AjC0FDMMcYVjQM0wGEHBUn4XmuF1uaP7eDM19JGbgcJO+JVLKNmjN/ab7SimnSSkdgU+98TxCiCJooSy0fZ2U8mwpZYWUcja0Tn6fWBwrWnIATLKxnTGN46B1WN2pLMsCcAmCXUdThBDnQev4eqaUcrey6itov4WK8f9wqNfzLQCPCiEaAXwKIFsI0SiEiBhmotxTOy46AEAI8RCAb0CL4VfDXszyUb2U0o6gt8yDUspvK+lcHdh2shAi32L7VinlIinlCCnlNGjvqU8sjhUtseajoLyucCns5aPZ0Nz1C6WUnymrvjLZPJp8BChhNlLKR6WUk6WUVdBagTyILoxLGP63dX0m+zEMc6RJdRA/f/jDn8gfaB1M3wdQYLLuLAD10EYEMRvlJg9ap8ut0MJECnC4g9vDAFZAc+lnI/IoNw0Argh8N45yUwstJCEfmhDYCaXjo8nxwqXr/wMwDIdHY3EDONuwf7hRbs6ExSg3AE4H0IzAqCph0jcfWgfNE03W5QLYDeByaCOTfA9a58ZhFseyvJ5AmkcEPrOgjbwzIlL6YsxHt0ETq+Um62YAaIUm9s1GuckN/EZvBa6nAEBWYN01ANZBG8XF7ig3vwnklasQPMrNJGjhTTnQRitqAFAV5ljh0nV54P4KAMdDG7nm54b9rUa5sczrgW0mQwubGh7hns8I5LdzLdavAXB74HrnQAu/mWSx7bmB3yYLwKmBbb8ZWFeIw/1ZxkIrL662OM40aGEz2QBKoTn9b9q5PmgV0NHKNhsA3JzovMof/vAnuk/KE8Af/vAn/CfwcpbQmuc9ymeess0t0MRnG7ShKIWy7r3A/upnXGBdIbSQBk9A7FwSIS2zAi/wHmjDBI5V1v0XgJbAsf4Fi4qBzXS9GjhWJzRX9Ccm+4cI+sDy70CrTHQDWA6g1LD+GQC/M9lvHgCP8v+70MS1es//pKw/Bpqj7gHwOYA5yrpbAaxQ/o94PYHtxiG5w1ZKaCOfqNe0SFn/Q2hivBNaTLU6Es3TJr/XtwLrsqD1I2iHJl5/HSEdk6AJ2R5oFYFjlXWXQGs56gKwFsDsCMcKl67HADQFrnMHtBFzhGH/EEEfKa8H1v83gOUm+1F42pjA/08B8Bnu+QrD9m8Hlm8BcJ6ybhGATcr/L0B7xjsDafuBsq4MmuPfBe1ZvsmQrhUIDJcJrTKwLbBtIzQXvtLm9V0fyCNd0Cq1/40kVD75wx/+RPchN4xhGIZhGIZhmAyEY+gZhmEYhmEYJoNhQc8wDMMwDMMwGQwLeoZhGIZhGIbJYFjQMwzDMAzDMEwGw4KeYRiGYRiGYTIYFvQMwzAMwzAMk8GwoGcYhmEYhmGYDIYFPcMwDMMwDMNkMCzoGYZhGIZhGCaDYUHPMAzDMAzDMBkMC3qGYRiGYRiGyWBY0DMMwzAMwzBMBsOCnmEYhmEYhmEyGBb0DMMwDMMwDJPBsKBnGIZhGIZhmAyGBT3DMAzDMAzDZDAs6BmGYRiGYRgmg2FBzzAMwzAMwzAZDAt6hmEYhmEYhslgWNAzDMMwDMMwTAbDgp5hGIZhGIZhMhgW9AzDMAzDMAyTwbCgZxiGYRiGYZgMhgU9wzAMwzAMw2QwLOgZhmEYhmEYJoNhQc8wDMMwDMMwGQwLeoZhGIZhGIbJYFjQMwzDMAzDMEwGw4KeYRiGYRiGYTIYFvQMwzAMwzAMk8GwoGcYhmEYhmGYDIYFPcMwDMMwDMNkMDmpTkCqEEIIAKMAdKY6LQzDMAzDMAwTwAlgv5RS2t1hyAp6aGK+IdWJYBiGYRiGYRgDNQD22d14KAv6TgDYu3cvXC5XqtPCMAzDMAzDDHHcbjdGjx4NRBlBMpQFPQDA5XKxoGcYhmEYhmEyliEv6O3S2dmJzz77DG63O9VJyWhcLhdOOOEEOJ3OVCeFYRiGYRhmUMCC3gZPPfUUHn/8cQwMDKQ6KYOCnJwc/OQnP8EVV1yR6qQwDMMwDMNkPCzoI/Dee+/h0UcfxaWXXoqLLroIlZWV0AbIYaJFSonm5ma8+OKLeOyxxzBhwgR861vfSnWyGIZhGIZhMhoehz4CK1asQG1tLX75y1+iqqqKxXwcCCFQVVWFa6+9FrW1tfjnP/+Z6iTFxcaNG+Hz+VKdjKTzxRdf4B//+Eeqk8EwDMMwjAVJE/RCiAohxBtCiC4hxFYhxHyL7QqFEEuFEJ1CiHohxMUW260QQgwYlu0RQnQLITyBz58SfR27d+/GjBkzEn3YIc8xxxyDPXv2pDoZMdPQ0IBjjjkGjzzySKqTknS+8Y1v4IILLkh1MhiGYRiGsSCZITePAmgEUAHgNAAvCSEmSylbDdvdBWA4gGoAtQBWCCG+kFJupQ2EEAuhDbJvxgIp5QcJT30An8+H3NzcZB1+yJKbm5vRfRJ27doFQBv2lGEYhmEYJpUkxaEXQjgALATwGyllt5RyOYCNAM4z2fxSAPdIKd1SyrUAXgVwiXKsAgD3ALg5zjTlCyFc9IF1BSEs3/3ud/Ff//VfQcseeeQRzJkzB9OmTYPD4YDD4UBWVhYKCwv1/1evXh1yrPfeew+zZ89GSUkJysvLMX/+fNTX1+PHP/6xvl9ubi7y8/P1/x966CEAwF/+8hdUV1fD5XLhxz/+Mfr7+/Xj3n777aitrUVWVhaWLl1q+9p27NgBIYR+rurqavzmN78JWl9QUBCy39y5c/HCCy8AAKZOnarvT5+cnBwsWLDAdjoygd27dwMAysvLU5wShmHC0dvbi+3bt6c6GQzDMEklWSE3kwF4pJTqTKwbAUxTNxJClAIYAWBDmO1uBvACrGd1fVkI0SSEeEUIMTZMmm4B0KF8Ypol9g9/+AOeeuoprF+/HoAWenH33XfjiSeewKZNm+DxeODxeDBlyhSsWLFC/3/evHlBx2lra8N3v/td3HrrrWhra0NdXR2uueYaZGVl4cknn9T3u/zyy3H77bfr/994441Yv349brjhBrz22muor6/Hzp07cd999+nHnjx5Mh555BGcfPLJUV9ffn6+fq6PPvoITz/9dFSx7lu3btX393g8WLduHRwOB2688cao05LOkEDo6+tLcUoYhgnHr371K0yZMiWjWwQZhmEikSxB7wBgHLDdHVhu3A4Ing1L304IMQ7AhQB+a3GeSwCMg1aBqAfwqhDC6pruB1CifGrCX4I51dXVuOeee/CTn/wEfr8fP//5z/Gzn/0M06ZNi7yzwtdff43i4mKcd955yMrKgsPhwPnnn4+amsjJeu6553DhhRfi+OOPx7Bhw7B48WI888wz+vrLLrsMCxYsQGFhYdTXpzJmzBjMmTMHW7ZsiWl/r9eLH/zgB/j5z3+O0047La60pBs7duwAAHR0dKQ4JclFbfkZCh2AmcEHmS91dXUpTgnDMEzySJag9wAwTr/qCiw3bgcEh7+o2/0fALdLKXvNTiKl/FBK2SuldAP4NYBJACZabOsNhPW4A9tHNaWuylVXXYW8vDycf/752LJlC2677baoj3HUUUfB4/Hgpz/9KVatWoXOTvvJ2bx5M4455hj9/xkzZmDXrl3o7TW9TTGza9curFmzBieeeGJM+19//fVwOBy46667EpqudKCxsRHA4Bf0TU1N+veurq4UpoRhYoPC4rZt25bilDAMw5iTCMMsWYJ+OwCHEKJaWTYDwCZ1IyllG7SOszMstvsWgEeFEI0APgWQLYRoFEKY2eEy8Dfp40oKIfDEE0/g1VdfxR/+8Afk5+dHfYzS0lKsXr0aXV1dWLRoEYYPH47LLrsMHo+xzhOKx+OBy3W4vkTf7ewbCa/Xi2HDhsHlcmHixIk46aSTMGfOnJD16mft2rUhx3nllVfwwgsvYNmyZcjOzo47XelGT08PgMEv6A8cOKB/T0T+YpgjDbVUchw9wzDpyJtvvomcnBzs378/ruMkRdBLKT3QOrfeFRiW8hxoQv1Vk82XArhNCOEUQsyG1nF2WWDdVADHBj7fAeALfN8qhBgjhDhZCJErhCgG8DCAOgA7k3FNRmprawEARx99dMzHmDFjBpYuXYqmpiZ88MEH+Pe//40HHngg4n4OhwNu9+GIJvrucBgjmqInPz8f7e3tcLvdaG5uRnNzM2655ZaQ9ernpJNOCjpGXV0drrzySjz99NO2Qogyke7ubgCDX9CrIp4deiYTaW3VBlZjh55hmHTk448/DvobK8mcWOpnAEYBaAGwBMAPpJStQohFQgjVqb8DQCuAAwBeBnANDVkppWyWUjZKKRsBHAwsa5RSDkAL03kCQBuAPdDi6M+VUmZkoO+sWbOwcOFCfPXVVxG3ra2txcaNG/X/N27ciAkTJpiOQBMPFRUVWLhwIVauXGl7n/7+flx00UX40Y9+hLPPPjuh6UknhopDr3b6ZYeeSSS///3vg1qAksWhQ4eC/jIMw6QTZMZu2LAhwpbhSZqgl1IelFJ+R0pZJKWcIqVcFVj+nJRymrJdj5RykZTSIaUcLaVcZnG8PVLKHOX/TVLKGYH9KqSU50opj4g7nwg2b96M//3f/9VfaFu3bsXrr79uK179kksuwUsvvYR169ahvb0d9957Ly677DJ9fX9/P3p7e+H3+/XvUsowRzSnra0Ny5cvj6oV4tZbb4Xf77fV0pDJDEVBzw49kyh6enpw7bXXYtGiRUk/Fwl5emYZhmHSiZaWFgDAZ599FtdxkunQMwbuvvtunHPOOQAAp9OJ1atX4/jjj0dxcTEWLFiA7373u7jhhhsiHue4447DQw89hLPOOgujR4/GhAkTgsJifvSjH6GwsBDvvvsurrjiChQWFmLNmjUhaQCABQsW6GPbA1qMPI0fP3nyZJSVldmeDdXn82HJkiVYv349SktLg8ainzlzpq1jZArd3d0oKytDe3t7qpOSVFjQM8mAKsLJni1aSsmCnmGYtIYEPQ22ESvJnCl20BPO9f76669Dlt1+++3699GjR+Pll1+OeI4nn3zSdPmVV16JK6+80nTd0qVLLSeUUtMAAG+99Zb+fdKkSWGvadKkSaYj6XzwweGJev1+v+X+g4menh5MmDBh0He045AbJhlQvx96kSWLzs5O9Pf3Iz8/nwU9wzBpCZkOat/IWGCHnmGixO/3w+v1oqysDAMDA4N6fHZ26JlkQA59vC+wSNDxR4wYoXdkZxiGSSdI0EczfLkZLOhtEEv8OROeTL6n1EoxbNgwAFqY0mDF6/VCCIHs7GwW9EzCUIV8MssCEvHl5eXs0DMMk5awoD9CFBcXo62tLdXJGHS0trYmZJjNVEAiobS0FMDgFvR9fX3Iy8tDcXExh9wwCUPtTJ7MfEUivqysjAU9wzBpSUtLC8aOHQuPxxNX2DIL+giceOKJ+OCDD5Ie6zmUaGlpwQcffBDzDLSphoQBOfR2Z+h1u9349NNPk5auZECCvqioiEMWmIShOvSJnuFahfLs8OHDWdAzDJOWtLW1YcyYMQDiMzi4U2wELrjgArzxxhtYtGgRFixYgIqKCmRlcT0oFvx+Pw4ePIh//etfyM/PxwUXXJDqJMWEUdDbdeivvvpqLFu2DF6vF3l5eUlLXyIhQZ+Xl4f+/v5UJ4cZJKgOfTIFPT2r5eXlXCFlGCbt8Pl86OvrQ1VVFQDN7HC5XDEdiwV9BEaOHIk///nPePbZZ/Hee+8N+nHHk43L5cIpp5yCSy+9FCNHjkx1cmIi1pCbgwcPAgC2bNmSMcN49vX1IT8/H/n5+UEdZBkmHlSHPpkha/SscsgNEwstLS2oq6vD8ccfn+qkMIMUKv8qKysBaHH0LOiTSE1NTdA478zQJlaHfvz48QCA9evXZ5SgJ4eeBT2TKI50yE15eTkGBgbQ39+P3NzcpJ2PGVwsXrwYjz/+ONatW4djjz021clJGfQcFRUVpTglgw/SE6qgjxWOHWGYKIlV0BNffvllwtOULFjQM8mgo6MDOTman3SkQm7U/xnGDkIIAMCzzz6b4pSklokTJ2LUqFFR7TMwMID9+/cnKUWDB6Ogj2coXxb0DBMl5FZEK+ip5t3c3JychCUBFvRMMnC73aioqACQfIc+JydHb8JmQc9EA/Ubinc4wUynsbEx6nDj++67D9XV1Tw6WgSo/GOHPkPw+XyDevKhoQaJgmhj6OlBzSRRQR14WdAziaSrqwvDhw8HkHyHvrCwEIWFhfr/mcytt96KDz/8MNXJGDLQ3BssSjVaW1ttb7t7924AwLvvvpus5AwKMsKhF0JUCCHeEEJ0CSG2CiHmW2xXKIRYKoToFELUCyEuVtadJITYIIToEEI0CyH+JoRwKOsnCiHWCCG6hRBfCCHSMjB53Lhx3KlmEBFryA0J+iM52sbJJ5+Mhx56KOb92aFnkoHX64162NdY6O7uRlFR0aAQ9B0dHbj//vtx+eWXpzopQwYS8jypnkY04aJUYX/rrbeSlZxBAZVJLpcLeXl5aevQPwqgEUAFgBsAvCSEKDPZ7i4AwwFUA7gQwGNCiKmBdTsBnA1gGIBxgfT+Rtn3eQCrAJQB+DOAV4QQadfRt6GhARs2bEh1MpgE0dPTAyEEnE4nAPuCnl4OR1LQr127FjfddFPM+7OgZ5JBb2/vEZlpmRx66syXyUNXfvzxxwC0eGbmyMAOPYImOtq7d6/t/drb26PeZyhChkZhYSGKi4vjKqOSIugDLvpCAL+RUnZLKZcD2AjgPJPNLwVwj5TSLaVcC+BVAJcAgJTyoJSyXmpzgwsAfgATA+eYCqAWwH1Syl4p5R8D1zPPIk35QggXfQA4E3nNzNDB6/UiPz8fBQUF+v92ONIhN9pjEx80bGVeXt6gnhGXObKogp4denusXbsWADBixIgUpyQx9Pf3Y+7cuWk9SAAL+uA5I6J5fkjQ81Df4aF7SsZD2gl6AJMBeKSUDcqyjQCmqRsJIUoBjACwwWo7IcQYIUQ7AA+ACwD8IbCqFsA2KaXXal8DtwDoUD4NFtslFHUynkQILCb1kKDPz8/X/7fDkQ65SUQzMTv0TDLo7e1FSUmJ/j1ZDCZBTyOGDJYOmg0NDVizZg3uvPPOVCfFEhLyQ1nQq3Hz0by7WNDbg8qkgoKCtBX0DgDGyH53YLlxOwDotNou4NAPA1AF4H4A9VGeg7gfQInyqYl4FQmgqalJ/z6UC4XBBAn6nJwcCCGiEvR5eXlHTFSoBXGslUkW9PHR0NCA5cuXpzoZaUdvby8KCwuRl5dnW9C73e6o8/Fg6hRLFfTBIuhpor10bvljh16bXIuI5vlpa2sDcFjYM+aoITfpKug9AIxTXbkCy43bAcHhL2bbQUrZDOBNAMuiPAft7w2E9billG4EVyKSxoEDB0y/M5kLiVwhBPLz8229kPx+PzweDyorKxPq0Pt8Pjz22GOmaVAFfaw95+laeabY2Jg3bx7OO88s0nBo09vbi4KCAhQUFNgS9D09PSgpKcHvfve7qM5DDn1eXh4AZHQepnIjnlEw0gkyu9L5N/F4PCgpKRnSgp7eI9nZ2ezQJwHVoU/LGHoA2wE4hBDVyrIZADapG0kp26B1nJ0RbjuFHACTAt83A5gshMi3uW9KYEE/+CCHHoBtQU8PaVVVVUIF/UcffYRrrrkGDz/8cMg6ckgAzSmOBXbo42PPnj0A0lu0pIJoBf2hQ4cAAE888URU5yGHPtMEfVVVFe67776gZVRuDBaHnubjSOdWk66uLlRVVQ3pUW5I0FdXV0ct6EePHo2Ojg4ONw5Db28vcnJykJOTg6KiorjyWlIEvZTSA61z612BYSnPgSa2XzXZfCmA24QQTiHEbGgdZ5cBgBDibCHEVKExCsC9AN4JnGMrgC0Abg50eL0KWqfZ1cm4plihZkUgsyYUYqyJRdDTS7iioiKhLzA695o1a0LWqQ59rC4Jj0OfGKIZv3ko0Nvbq3cst/P8UOV0x44dUZ2nu7sbhYWF+vOaKXm4ubkZixcvDlpGL/rB5tCnq9Hl8/nQ29uLqqoqeL3eoP5wQ4nOzk5kZWWhvLzc9rtLSon29naMGzcOPp8vrUeXqq+vT2n6yHQAkLYhNwDwMwCjALQAWALgB1LKViHEIiGE6qLfAaAVwAEALwO4JiDWAa3D7ApoYTSfAagD8FNl30sALADQHlh+vpRyIHmXFD09PT3IytJuczpnasY+sQh62qasrAz9/f0YGEhMNiWhuHHjRst1QOxSGJi/AAAgAElEQVR5jx362FEr82ocKhO9Q0+Cvq+vL8jt6+zsDCu0KFY/Ozs7qv4uqcTqfiTKoe/o6EgLc4kEfUNDg20Ht7+/P6hfWjKhChSNKjRUXXoKW4smHMTj8cDn82Hs2LEA0jfs5uuvv8bYsWNx3HHHpSwNGSHoA0NOfkdKWSSlnCKlXBVY/pyUcpqyXY+UcpGU0iGlHC2lXKase1JKOUFKWSylHCWlvFJK2aKs3yGlnCOlLJRSHielXJ+s64mV3t5eOJ1O5OTksKAfJMQi6EkM08geiXLpSbSbiQA15CbW82WaoF+/fn3apJPCRIzfhzo+nw/9/f0xCXrgcAfFlStXwuVy4Re/+IXlfl6vFwUFBRBCZEwetqr8dXV1weFwoLOzM64QhgsvvBBVVVUpD3UhYd7f32/7d7nzzjsxYsSIiGWulDJu04QEfFVVFYCh2zG2q6tLHynKbp6he1ddrUVdp2vH2Pp6bYyVbdu2pSwNZG4AaSzoGQ36seL9oZj0gUQuEL2gp7G3E5UXSOiYFbQdHR0oLy+3db433nhDj/dWUcehT4YY+uMf/4i6urqEHGv//v047rjj8MADDyTkePGi3nN26A9Dz0usgp7cvq+++grA4SnmzaDQHgAZ07HbqvLX3d2Nqqoq+P3+uMoPes6XLVsWfsMk09LSAiEEAPvlIU3QSJNsWXHDDTcgNzc3rvSRKK2srAQwdAV9d3c3iouLo9Iw9JzRvUtXhz4dyuWMcOgZDXKIWNAPHmJx6CksgAR9oh367u7uENfO6/WitLRUXx+Os88+G7NmzQpZnkyHvre3Fz/72c9w4YUXJuR4NMV4usxMqP7G6fDiSBdUQZ+fnx+zoCfXL1wohOp+ZYpDbyXou7q69PCPeOLoqQxK9XPS3d2N4cOHA7BfHo4ePRoA8O6774bdbsmSJQCCZzmNJX0AC3p1Lge7vxM945kk6FNVNrCgzyDYoR98JCLkJlF5QY2TNwojr9erD9kX7nw+nw+AuZDo6+tDbm5uUsQQjbxDE+bEy9tvvw0gWPylEiuH/tlnnzXtxDxUoHwar0NPy8LFlKvPaqbMdqw+h2plpbu7Wxf08cTRU9+OVFcye3p6dEFvtzyk337Lli2W26giPp5QD8oriTZhMg0KuYlGw9C9o983XUdmUp+BVHU255CbDIKafFnQDx7iEfTJCrkBQp1Kctcj5b1wI7D09/cnTdBT/GKiZgqlCsL27dsTcrx4oXteWloaJNIuu+wyzJ07N+bjtrS04OKLL06biku0xCroy8rKAIQK+nAv4kx36BsbGwFoIrWnp0e/B/E8M3T8VAv67u5uPSTQrlimsiqcUFdD+OIZXYryisulTXeTzBmN05l4Qm4ov6Zr60a6Cfri4uL0G7aSOQz9WIWFhSzoBwkUVw7Yj8tNZqdYqzh5VdCHO586Ggu59cTAwIAu6L1eb0LHEyZBn6h7oQ5tmA7jHtN1jRw50rTJOdby4K9//SteeOEFvPqq2SjA6Y8q6O0+P21tbRg3bhyA0JCbcC9iY+U70wQ9XRvlJRJIsbY0eL1e3S1NtaDv6emx3ceHIIEerjKrrovnGpNVZmca8YTc0H7p7NAXFRUBSJ2gV/vksUOf5gzFkJuOjg7cfPPNGfHyjAUamx2w7/oly6FvbW1FTU0NAHOHPj8/P2JlUh3CzjgBVX9/P3JycpCXlwcpZYjgjwcS9F1dXQkZBaGtrQ1Tp05Fd3d3QmM2V69ejfXrox9Ai+55RUWF/tuoI2/EckzgcMfAVI7MEA+qoM/NzbX1/LS3t6O6uhpZWVmmITdWFbhMdOjVZ4EEFOWleB16ErjV1dUpF/TxOPThBL367CfCoSdBP1Qd+nhCbvLz8+FwONLaoR8/fjyA1MX5q6ZDUVER+vv7Y57zgAV9klEF/VCp4S9evBgPPvggPv/881QnJSnEEpebrOZbj8ejx9UaC1uqeEQqiFWH3vgCHBgYQE5OTlIm5iFBD5gXpp2dnVEJ/ba2NkyePBlAYoeJPOWUU2Iap7inpwf5+flwuVz6C00VIrG6VhR//9lnn8W0f6pRBb1dkd3T04Pi4mK4XK4gQV9RUYH+/n7TZ3BgYAA+ny9I0GdCDH1vb6/emZ3eGVQhjFfQ07M+derUlE92FqtD73Q6wwp61WmNp9JCeWWoC3oKuYnGoadnOi8vL2MEfbo49EDsrUEs6JPMUHToaQSCdH2I48Uo6KNx6J1Op36MRNDV1WXZscxuDL3q0BsLEjXkho6ZKNra2vQWC2NeOXjwIFwuFxYuXGjrWAMDA+js7MSkSZMAJC6cQL1v0To4NEtpcXGxfn1qumKJlezr60NdXR0qKioGhUOfl5dny42icrSkpCRI0I8ZMwaA+ctYdQnpbyY49L29vSEdMY0OfazlB1V0p06dmlKHXkoZ5NDbeTf6/X60tbVh4sSJaG9vt2yVofwhhEiIQ19cXAwhxJAx5IxQyM1gdOhbW1vTzqEHWNCnLUNR0NO40Klu0k0WscTlJsuh7+7uRkVFBQDrkJtoHHq1IKEQGwq5Ua8jEbjdbowcORJAqKBftWoVANhu5SEnP9EOPYW3AMCHH35ous3DDz+Mf/zjHyHL6UXocDj030Z9JmIpD/bs2QO/349Zs2alfLIWt9sd07CAsTj0NPyvKujb29t1QW/W2qGeB8ickJtkOvR0n8aPHw+32x1z03680Iy/0YTcdHR0QEqJiRMnYmBgwLJC3NHRgfz8fIwaNSohMfT5+fm2O28PRtSQm97eXlvPvFHQp2sMvdvtRlVVFXJyctLKoY9VK7KgTzJDUdBTnPVgFfRqp9hoHXpyexLh0Esp0dXVpQv6WENu1I616ouV4r2T5dC73W6MGjUKQKigpyZ1qgBFgsQtOfTJEPT79u0LWd/T04Mbb7wRF1xwgek6EvSJcuh37twJAPjGN74Rs6BOBP39/SgpKcHdd98d9b6xCHoaLYwEfU9PD7xerz61vB2HPlNCbrxery7c6XmkvyT0Y70OOg71u0lV2A2VR8OGDYMQwta7kZ4hGoveKuymo6MDJSUlKCsrS4hDT/2QBqNDf8cdd2DRokVht6GQm2jCLtWQG6fTmbYOfWdnJ1wuF1wuV8oEPTv0GcRQG+XG6/XqD/NgFfTGTrF2Y+izsrL0ePREuD19fX3w+/2Wgt7uKDcej0efAERNF7l3qXDo6WV96NAhWyPW0PaVlZVwOp228159fT1Wrlxpub6pqQlVVVWoqKjQhxBUobHvySlWoZAbM0EfaW4AK3bu3Im8vDxMnz4dUsqUNRNv3LgRwOGWlGiIpVOsMeSGfu9wITeZ7NA7nU5kZWXpzy1dS7zx3HS8qqoqAKmLG6Z0kPNrR8BQOUt9hiIJerVlLBa8Xi+ysrKQnZ0d4tD39fWlvA9CIrj77ruxbNmysNdCLY3RvAMyIeSmr68PXq835YKeHfoMYqg59KrAGMyCPhaHnh7aRDXfqs3wQoiYQ246Ozv1SoGZQ58sQd/Z2akLemPayXHv6+uz1VxLL/fS0lKUl5fbcuh7enowduxYnHnmmZaVMmq9GDFiBJqamkLWk7ClCpHx+EVFRSEx9E6nE8OGDYtJbOzYsQPjx4/XW1RSFXazdu1aAObXbWTDhg2YN2+ennco7+fn58ccckPXbUfQRxse5/F4sGnTpojbJQvVBDIK+sLCQuTl5cUl6PPy8vSKQapCIag8iiY22yjorfI+Cfp4XXW1zDYea8GCBfozmMnQNaxYscJyGwq5ieYdoBpY6SroKe87nc6UhlSpeqK4uBhAGgp6IUSFEOINIUSXEGKrEGK+xXaFQoilQohOIUS9EOJiZV2VEOI1IUSzECLEphNC7BFCdAshPIHPn5J1PbEy1AS9+mJN5Egj6USsMfRUINqdjCoSlJ+sJv2gloRIrUMejwfl5eUhHb+ORMhNVVUVhBAhBb76slZj/K1QBf3w4cNt5T11iE61Y7BKS0sLysrKMGLECFOHnl4KZg6X6tCTeG9vb0dpaSmKi4tjKg82bdqE2tpavdNkqiaX+uKLLwDYq7Tfcccd+OCDD/QJf7xeL7Kzs/WKot1OsWrIDV03hY6Y3Ut6xqJ16H/xi19g+vTpKSuvwwl6Grs/1vKD8iR1zk+VoKfrKiwstC286R7YdejjFWnhTJj3338/6uM9+OCD+OSTT2JOTzI5cOCA6XKfzwev1xt1yI3aip2ugp60isvlSqmgz5RRbh4F0AigAsANAF4SQpSZbHcXgOEAqgFcCOAxIcTUwDo/gH8CuCzMeRZIKR2Bz1UJS32CGKqCfsyYMezQKyTToS8uLjYViHZHuens7DR1KZIZcuP3+9HZ2ak3jZuF3JB7b0ec0/4OhwPl5eW28p4q0K1eaK2trWEFPeV3K0FPMfRerxf9/f16zGZRUVFMDv2mTZswbdo0XdDbceillKZpjwf6Tcz6FRihPLV//379/2hFtlXIDYk7M4FrFnJjRwjTtdFoXUcaao1QK+KxzK5rRk9PT1oI+ngcegoXSragV8v5wsJC/VjqTLR2+7DQ3CxnnnlmzOlJNAMDA3o5aVVeqqFR0Ybc0L1L106xRoc+Vf1rzGLo08qhF0I4ACwE8BspZbeUcjmAjQDOM9n8UgD3SCndUsq1AF4FcAkASCkPSin/CCC2GVjSAKOgT4cZLJMJhdxMmDBhUAp6n88Hv98f0zj0ubm5ABLv0NNLMdaJpTweDxwOR4hTlkyHngS40+k0FfTt7e36iDV2HPru7m4UFBQgKyvLtkOvilwrwUshN1VVVaYhNyTo29vbQybdIvHkcDgAaBUwt9sNp9MZk0Pf1taGAwcOYNq0aXrnSDuC/osvvsDIkSPx4IMPRnW+cNBz3tDQELFMIyFELSKqoLcbQ28VcmPW90PdB4g+5IaGgY2lf0AioNYI9Xmka8nLy4tLfKSLoI/FoadrpucnkqBXRXgkPvnkE/zud78LWmY0YSiNe/fu1bexe/8++OADPW3pglpGWr2r6Z0SS8gNPXfp2ilWdegT1a8tFjLBoZ8MwCOlVKed3AhgmrqREKIUwAgAG8JtF4GXhRBNQohXhBBjrTYSQuQLIVz0AeCM4hwxQ7WvgoICSClTNkzYkYIektGjR6dlrTxe1Bcr/bVTwPX39yfVoTcT7WrITTgB4PF44HQ6Q16syXToKW+4XC5Lh55GrLEr6KkwjEbQ5+TkICsry9KhjxRyo4aYGcU1pYniIru6uuJy6GlCqWOPPVaPgbYTckMj4yxZsiSq84Wjvb0dFRUV6O3tjZgGum9mgt7O8zMwMICBgQE95MbtdqO1tVUXvVYCN9ZOsSRuzCpwRwKrkJuCggIIIeISH2q/DiFERjr0+fn5KC0tteXQ2xFGfX19mD17Nn75y1+GLDcrs8M981bQs+v3+9MmFJXy97BhwywFvRrWGa1Dr4bcpKMWSMcY+tzcXOTk5KSXQw/AAcDYS8kdWG7cDgA6I2xnxSUAxkGrQNQDeFUIYXVNtwDoUD4NFtslFLVwpv8HM+RAjBgxYlCGGMXq+iUz5KaoqAj5+fkhlUU6Z6QXW2dnJxwOR8h2aqdYut5ENUuq7oiVQ19ZWYnCwkJbow+ogj6akJtRo0ahsrIyrENfVlaGsrIydHZ26vdEvQ4aetN4TqND7/F4dIc+lhC8559/HtOmTcPUqVORnZ0Nl8tlS1DQtSXyeezo6MCECRMAhBc1Ukp9XopYBb0aC19SUgKfz4d9+/bprRRWLV6xDltJv2Oq3FQrQU/XkYiQGyGEqdA6cOAAbrjhhqSPBkR5kRx6O3lT7eQ8bNiwiJ1i7d6nLVu26N/V1iarTrFqeWS3D8uBAwf08mnXrl229kk21G/o6KOPjijo4wm5cTqd6O7uDmnBTDXpGEMPwPaoT2YkS9B7ABgHkHYFlhu3A4LdcrPtTJFSfiil7JVSugH8GsAkABMtNr8fQInyqbFzjnjw+Xzo7+/X4x6BwS/o3W637qAMBUEfSwx9MjrFmqWDmj0jFVaqQ69ul8yQG7UwLS4uDnGraRZZp9Npy90xOvQtLS0RQ0EaGxsxYsQIjBw50lTQ+/1+XdBTzLpR5LndbowbNw5AaBw9vdRUQU8Ovdk1h6O+vh4vv/wyLr/8cgghAGjOmh1BQddGLQSJoL29XR8PPFyFq7e3V88zVoKeJjCzwijoAW2CLRL0Vvnb6NDbrXzT7xjrCEJtbW348ssvY9oXCI6hNzr0QHzlBwl6AKbP1l/+8hf89re/xQsvvBBz+u2mA9DEi90QIjUf2HXo7bxv1WdaFVKqKLVy6O0K+ra2NtTW1gKwDu870tiZNdgs5MZuiCltT2VnqoaFtILyPplZqYihl1IGhScBiKu/ZbIE/XYADiFEtbJsBoCgscCklG3QOs7OCLedTejtLUxXSukNxOm7AxWApLcBqQXQUBH0HR0dekjBYBT06mQjgCZIfD5fRPch2Q69mftIzZ4FBQX6mPVGfD4furu7TWPoYwm5ef311229sNTmTqNDT+Orl5SUxCToy8vLMTAwEPEF0tjYqI8xbxbWQxM3RRL0NHSiUQDSb06CvrOzM2aH/k9/+hMcDgeuvvpqfZndSoGdzr/RQL9PuCEjCfpdXS6XLhqMgh4In69UZ5YEfV1dnf6bWAlc47CV0YbcxOrQz58/H8cee2zMLptZDL16z+IpP2iUG8Bc0FOla+nSpTEdP5p0CCH08snO9dgJuaEZZKOJoVefW/U3T6RD39raiilTpoQN7zvSUPkzevTouEJuXnjhhZAWVrUyRM9sqme2NuJ2u+FwOJCVlRV1GNubb74ZNOlgrNA7VnXoYx0BDUiSoJdSeqB1br0rMCzlOdCE+qsmmy8FcJsQwimEmA2t4+wyWimEKACQT9+FEPR9jBDiZCFErhCiGMDDAOoA7EzGNcWCcWQCIPbODpmC2+1GSUmJ3myUqpksk4WZQw9EFrrJcuiFEKYzblLNn2LoAfPKJAnCRHSK9fl8OOecc3DGGWdETLsx5EYVFl6vF36/Hw6HI2aHHog8Ok5bWxvKy8stz0FpLCkpsXwpud1uy46Z5LzQbLdutzvIoY+m0P7ss88wd+5cvXIA2Bf0Bw4c0N1BGmkmHjweD/x+vy2Hnl70o0aN0rczdooF7Al6K4fe6mXs9XohhNDPYUfQSynR2tqKysrKmARIU1MT1q1bBwB44403ot4fMA+5IdceQEI6xQLmgp5EijqSS7ysXr0aH374YUg6ioqK9PIrGoc+NzfXUtAbwyjsvG9VEW8l6I0OPY3AZTePUFlTVVWVNg49DSRQUVFh2aIZKeSmubkZF198Ma6//vqg/cwEfTp1CAYOj+4GWFeS29vbcdVVVwWVsz6fD9/+9rcxc+bMuMtTo54A0jPkBgB+BmAUgBYASwD8QErZKoRYJIRQHfg7ALQCOADgZQDXSCm3Kut7AOxRvtM6J4AnALQF1k8GcK6UMm0CtcwE/WB36GnEFBJXg+16zTrFAtEJerPC46STTsItt9wSVVpowg/qKKemwefzQUqpO2CA+W+hjjYT77CVNIShHeeCXrxmbrU6AkYyBb3b7TatUBBqp2OrYSI7Ozt1QW8shOk3pxea2+0OcujthtxIKbF+/Xocd9xxQcvtVgoaGxtx/PHHA0iMQ29nUicikqCPxqFXBX1TU5N+38N1iqWOpHSuSMKxo6MDPp8PEyZMiEmAfPrppwC0CjLNIhwNAwMD8Pl8EUNu4o2hB0IFvZQSGzduRFVVFfbt25ewEdlOOeUUzJkzJ2iZ2lJg93qo5UIIYSno6TeLNeRGzcvhHPqysjI4nc6oQm6oc326OPT0/igvL0dfX59peaSWgWbj0FNZQP1kCPXepaugp2gCwFrQL126FI8//jj+8Y9/6Mu2b9+uf6fKe6zQvTTG0KeVQw/oQ05+R0pZJKWcIqVcFVj+nJRymrJdj5RyUWAc+dFSymWG4wjDZ1xg+SYp5YzAfhVSynOllGnjzgPBL6Kh0im2q6tLn+gISGxHvHTArFMsEL1Dr+YDj8eDjz/+GA888EBUaenu7tZHUDG6j2o6wwl6NY7QyqHPyclBTk5OxOuk0VQobeFwu90oKChAbm5uyHlVVyjWkBsg8qRHkcJ61NYLM0E/MDCA7u5ufVxsM4eeKlR5eXlob2+PGEPv8XiCXhiAJsIPHjyIY489Nmi53UpBU1MTJk6ciOLi4oQ49PRiHjlyJLKzs20J+urqan2/aAW9MXaaGDtWG9QsXKdY1fmyE0NPeWbixInweDwhnaAj8dVXX8HlcmHBggXYsWNHVPtSmgHoQx0no1MsPSfGfN/V1QWPx4PZs2cnrL/Ftm3b9O/q8dR0RBNyQ/nGqlOsKujthtxYOfTGGHpV0LtcrrBx/EZaW1tRWlpqOVpWIrnvvvtw3XXXRcy79P6wCiekbQBrh56u31gOZYJDT5P8AdamAOVZtcVKNay2bt0ask80ZJJDP+RJtkP/4YcfBhWYR5pVq1bhuOOOC4o/plp/qgS93+/HggUL8MgjjyTl+IkIuTEWHqtXrw46pl2o8kTpUNOg1vzDhXupDn24kBuKdw13neroDZEmHCJhC4ROqx6roKd7QYLerkNvNU4yLSsuLtbTqooISldpaSlyc3MtHXpACwE4cOAA/H4/nE6npQg94YQTMGXKlKBl5H4Zl9sNuaEJvKqqqmwNARoJugfDhg2Dy+WKyqGXUpoK+nDD+aqx8A6HQ38JU2dkK4dXPQ+dy66gpxF8ou3I99VXX2H69OmYPHlyTIJevdYj3SmW8sbMmTMB2Js0LBL//ve/9e8fffSR/l2tgEcj6KmMtOvQ9/X1Rezf1NHRobfqWYXcqCKLyg27wsvr9aKnpwelpaWWHfABreP7ZZddFtc7c9OmTVi8eDEeeeSRiPMo0G8QbpK67u5u5OTkWIZdWgl6taNnOgt6tR+OWR78+uuvAQCff/65vmzDhg0YOXIkZs6cGbegzxiHnkmuoD948CDmzJmD008/PSHHi4XTTz8d69evDxr2i4RVqgT91q1b8a9//QvXXXddUpwQs06xQOSe/+FCbtav1+ZNoxetXajyROmwEvThWofCOfRqyI3ZOYyoDn2kwpteikCooI835KagoADFxcVhBb06U20kh764uBg5OTlwOp1BL71I8brGZmfqcOhyuUzvZXd3t/6CUMMdaHg5agkg7ITc+P1+PQzO7nCekVBFUzSCfmBgAL29vUFOa7Qx9AD0uFdy6MOF3KiV5Ly8PPj9/rACj0a4mThRGywt2jh6EvSTJk1CfX191MLb2Kqb6E6xdgQ9tQQlojWnrq4Oo0aNQlZWVlBYhpoOuxUU9fcsLS1Fb29vyH0wCnogctnc3t6OmpqaoP2BUEFPzxqVXXZbAEj0Rgq5ueiii/Dss8/qk1DFwvPPP49hw4Zh/PjxWLZsWdhtqcwMJ7iNphFgLujNOsWq77u8vLy0FvRWzxQJ+q+++kpftnv3bkyePBlTp05lh34w8Pvf/x633Xab7RkSkyHon3nmGQBarT4VHU/Va1eFs9GhP9KdgD/55BP9eyJ6oRtJRqdYKhTb29ujcgSNITfqMdWKh50Yehq6y2zYSitBf/DgQaxZs0bPf3v37tU7i0W6DoolBxLn0FOeAw4PXWmFx+OBlFJ36COF3AChzfyqoDd7uRsFfX19vb69WfiH2tqm3o+mpiZkZ2cHhZsA9kJu6F46nc6ECXq1VScaQQ9o9yyekBsAep6P5NCrFQfaTj2eGUaHPhoRMjAwgC1btuiC3u/3h8QWGzGW3VaCXr2WeBz6cKPcUAU4kQ59fX09xo0bh+HDhwe1DiXCoQdCR5kxE/SRjt3R0YHy8nI4HA5LQU+tYVJKXdDb7XRLlUTVoTdqB5/Ph88++wwAsHbt2ojHtGLFihU466yz8P3vfx9vv/12WI1i16Gn38ms8k3XZvwdjOFuNMNzomloaNANj2gxCnqzZ6q+vh4lJSVBfUoaGhpQU1OD8ePHx915nB36FCOlxLXXXot7770Xr732WthtzWLoEyVwN2/erH/ftCmWUT7tI6XEJ598ElQ4qDVyVdCn2qH/5JNPMGXKFBQVFWHjxo0JP34iYuiNLzC1ICXRZwfVoTcKRONU8UDkTrFWDr3VKCF33nkn5s6di0cffRSAJvAnT54MwJ6gT6RDr94LIPJsseqL3+l0oq+vL+Q3VB162tauQ+/z+eD3+4MEPYVgVFRU6PfS6plShXdzczMqKiqQlRVcXNsJuVFbYOzOoBsJSmdRUZEtQZ+Tk4OKigoA2n2PZ9hKALjgggsAQHdVo3HoI52rpaUF+fn5UY9iAgA7duxAX18fpk+frlcI9uzZY7n97bffjuzs7KA8oFZeUuXQ19TUwOl0JiQ8q76+HmPGjEFlZWWQ8FLTEUsMPQl6s1GncnNzo+q3Rn1paBZi9XyUf4qKivR5ZVRBH41DTzH0fX19IQK4vr4e/f39yMrKCgpNsgPln9bWVqxbtw6nnXYa5s2bh/3794etUFKZaVfQ04hRZg59a2trUMuXcWz1ZAh6n8+HWbNmoaqqKqbQYzsOfWdnJ44++mj09PTo6d+3bx9qamowatQoHDhwIK7O4+zQH2Ha2tqCfjBVJBqH4jKSTId++/bt+M53vgMgWNwngzvuuAOzZ8/GF198oS9TC2ejQx+LoJdSor29Pehef/bZZxEdLiO7du1CbW0tpk2bllRBH88oN0ZHsb29XRfC0bhi4TrF2o2hpxc6dXpSY5kjOfTUJE/Nks3NzXqoQqTCOxkx9Kqgj+RGq2KcWgqM5/F4PMjPz0d2djaA6Bx6Yz5RQ24qKyv1CZXUjmvhBD2N6KJiJ+RGrbAl0qEvKirSZ6uNJOidTmfQSD+q0I52lBsAuOuuu9DY2KhXNMN1ijXG0BvP5fP58NBDD+nX0NLSou/F9m8AACAASURBVLu1lH67ULP89OnTMWrUKGRnZ4d18Gjypvfeey/kWimGnn7fIyXoqX9HeXl5yERpsUCC3jjXg9GhtztsZSSHnjqf0nCYQGQTjQS9y+UKG3JD6Y5X0AOhk0uRIF24cGFQ6GIk+vr6MHnyZFx++eV48803IaXEqaeeim9+85sAEDZ8RzXfcnJyIobcAKEtwer9N5Z/qus8bNiwhAv6lStX6vcxUn8BM4wx9P39/UGVEr/fj66uLhx11FEAtPedlBINDQ2orq7GyJEj4fV64xpf38yhj3bSQRUW9GFYsmQJysrK8PTTT+vL1qxZg5ycHMyfPz9iSIdaOOfk5CA7Ozthgn7btm048cQT4XK5EjpmsBl/+MMfAAQP10SC3uFwhDj0sXSKvemmm1BaWooHH3wQgCYSZ82ahQkTJmDv3r2203rgwAGMHDkStbW1utBMJMnoFNve3h5TE7+dGHoSBoC1Q0+TaxjdF7VTLB1LTTeJQ2pVOHjwIGpqalBQUBC1Q68WpkZBT+Ex4TALuYnGoQdCBb3xZWZ8KYVz6I0FNQlaCp0xa9mJRdBH49DbEfSvvvpqxA7llGcA2BL0DocjaCz+aDvFGp+5rKysoP4E0XSKVY8HAB9//DFuuukmLFmyBECooI/mxbpu3Tp9orKcnBzU1NSELZtJKLz11lum11pYWKh36lQrQbGG3EgpQ0a56e7u1p+7gwcP6i0pZWVlcQt6v9+PvXv3RnTo6fezE8JK94CEmFHQHzp0SL8Guyaa2slV3dYYcgNo+SHaGHpjyA1gLujz8/Mxa9YsNDU1RTwm8dZbb2Hnzp145plnsGjRIsyZMwdjxoxBWVkZpk2bFlHQ07DHxtZH4zaE8T2j3n+rWXYBLT8lonVQ5cknn8TMmTNx3HHH6eFKdvH7/ejo6Aga5QYILhuoPJ46dSoAzWxraWmB1+tFTU2N/lvGMwwpO/RHkD//+c8AtI4mxObNmzFlyhSceOKJtgW9lbPi8/nwzjvvRD00mtvtRlNTEyZPnoyxY8eGbdaNF3XGTXUUE3JbZsyYYerQU2FtR9BLKfV7/OCDD6KjowOLFy/W71uk0CaVxsZGjBgxAhMmTIja3bdDrJ1i+/v7LUNu2traUFNTg6ysrKhq+3aGrbQTckPixejQR+oUSwV0XV0dpJR6aIjR6TLDGEMPHH4hGENu/H5/xHxk5tCHe4HYcei7urqCJnIyc+iFEHp+tyPoKXTGrCIYraC3E0OvOvTUryCccFq4cCGuu+66sMPxJVLQ2+kUq04oZIbdkBuzShSZFJS/WltbUV5eruclOw69z+dDW1sbVqxYgfnz5+vLx4wZE1bQU/40hnlQWtWKeCIcegrxUh164PA1qmK4rKws7tacpqYm9Pf364I+nEMvpQxbqQPsOfQHDx7UR6yxKvc+//zzoNBGakUyVspVl1nND9S6aDeGvq2tDYWFhSgoKNAdeqMI3LZtGyZNmqSPBmVX0C1fvhxTp07FZZddhtzcXNx44436urlz5+ojqJmh/gZWw4DaEfQUCmgs/9Rnr7q6OiF9Mvbt24e5c+fi9ttvx2uvvYYf//jHmDVrlj7/g11ocjw15AYIfo/T+4AE/f79+/VWVgq5oeV2oDBMFasYehb0Caa5uRlbt27F3Llz8c477+g/7ubNm1FbW4sZM2Zg3759YQWYsfZlLIh/+tOfYv78+XjllVeiShsJ+AkTJmDs2LEJdeiNvbYbGxv1TKgK5ObmZgghMG3atKDCiVxjK0F/4MCBkJj/DRs2oKGhAc888wx6enpw2WWX4ZVXXsHvfvc7nHrqqfjnP/9pK+0DAwNobm7GyJEjMX78eDQ3Nwe9kP1+P7Zs2RLX+Mperxc5OTl6IWYmzN55550QgWLWKZaEVXt7O8rKyixdEiuMIxCYdYq1E3JD4szKoQ8n6F0uF+rr69HZ2Ym+vj5UVlaGxKKaYXTo1fR1d3cjOzsbubm5lmJbpb+/H/39/VF1io3VoTcOW0mtG0a3zljxI0FPwtysIujxeJCbm4vs7Owgd5QqSkaKi4v1a7fC6NCrFXQjagucOpGKkXgEfSwx9CQOaIIoI3Y7xZqdi8LyqAwjh55+00gVpoGBAZx55pkoKyvDunXrcNZZZ+nrIpXNJHCNkxlRWtXnIhGdYtWKMoCQfJ9oh55EM4XchIuhByI76eo9KCwsRH5+vqmgp2uwapk84YQTdJEGHM6jxvezKkqpHCD3PFzIzfLly4PMr7a2Nr0CUlxcDKfTaerQT5kyRW95suvSb9++HTNnzsTTTz+Nnp4enHvuufq6efPm4euvv7bsNGoU9FYhN2q5auyrZTWxnjHkpqamRhfD8bB48WJ8+eWXuP/++zFt2jT88Ic/xLRp07B9+/aoBgdRh94FzPMgPRfDhw9HeXk59u3bp1dKYnHozzzzTH0yPsLKoedOsQmGeprfeuut8Pl8+jikJOgnTZoEAGHj3Xp7e4Mm5TG6eDQii3EimUioBWUiBf3SpUtx1FFH4Z133tGXUbjLUUcdFeLQ01TWVPBLKXXXODs7G/n5+SEZ86KLLsL06dODChkatvH888/H9ddfj+XLl+Poo4/GpZdeipNOOklfH4mDBw9CSqk79EBwp7RnnnkGtbW1OPvss6O4K8GYTVYDHH4RNzc3Y/78+bjiiiuC9jOG3Pj9fl0w0wQXVi6JFWqBbCxoVUFJabRy6OnFbhTs4TrF+v1+tLa24vjjj0dHR4eeh8mhN4q8zZs3Y9y4cVi8eDGA0Bh6IFjQU1OwHUFP+5k59FZuNLnrDocjbAx9OEGvVkoihdxQHwlKj5lb3NXVBafTGeKO0iyTRiht4Qp/o0MPWI/PT89ZTU0NPv7447DHjEbQFxcX6xXLWEa5UZ8dM6wErp1OsRT3TkZGS0uLfq8dDkdEQb9mzRqsWrUKM2fOxK9//Wu9wy6gjcITrpWQykAzQa869D09PREdep/PZ2tUEyC8oKc8kogYevU9VVlZiUOHDumiy1h2AZFbOY1lb2lpaUh5aRZyoz6XdH+ME/tZCXqjQ09CPJygP++88zB79mz9/9bW1qDn12zoShL05ODbFfR1dXUYO3YshBB6Xx/i1FNPBQD861//Mt1XFetWZpJaxgGh74je3l69smKMoVd/q5qaGjQ2NkYdjaAipcTrr7+OX//619i1axf+/e9/w+FwYMyYMfB6vVF14qaKIBktZu9Iei6cTidGjRqlO/TZ2dmoqqpCYWEhSkpKgn7L/v5+LF++PKRyIaXEqlWrQloprBz6WEexYkFvwfbt21FcXIwFCxaguLgYn376KVpbW9HY2Bgk6MNNHmKM4TQWAJQRVKFsh/r6euTm5mLEiBEYN26cHvIQL/fffz8ALY6WIEE/b968oIrDoUOHMHz48KCCwCiszGqa9BL7y1/+oi/bvn07qqurUVxcjHvuuQcrVqzAe++9h4KCAhx99NHYt2+freEcqbAlhx4IblV4/PHHAWiTnagTRQBawWhHTBudB6NIILH03HPPBe1ndOgB6HGj1DnHStC73W489thjIa0n4SaWUkNuqIOYVS9+q5CbgYEBCCGCWiPoHB0dHfD5fDj++OMBHJ54o7Ky0jTk5osvvkBdXR2eeuopSCnR0dERNuTGSnSYocbcE6NHj8bAwIBlcyidPysrKyT0gIjk0BtbGcwcevrNzzzzTACHHWGrkBuHwxHijlKnPSN0veFEp9rpmdw/q/kZ9u7di+LiYpx66qn48ssvLY8ZjaBX7yHli2g7xRqfOSNWITd2hq2k8ogqpOTQA1qFKVLIzdq1a+FwOPD5559jyZIlQSJm6tSp2L9/v2ne9Xq9+nJ1vfrcqs9FT09PkKAfGBgI6sC3dOlSnHbaaVixYoVlWo+0Q19XV6fPslxWVgafz6efKxaH3lhBo8mlrrzySnzve99DV1dXUKXE7Pc2Ovo0uhXNwxEphp7e2VaCno6vVppVhx5AyGyxvb29qKuri9qh9/l82Lt3rz58qxGa/GjlypWm69WQzXDvHrXsMRP05HJHCrnx+XxR9Q8wsnfvXrS0tOCEE07AmDFj9LKXXO9oRoij1oLq6moA4R16EvT79u1DQ0ODPkM2oN1jVdC/+OKLOO+88/R+h4Tan0/Nj1YOfaywoLdg9+7dGD9+PLKzs3H88cfj008/1UeTqa2tRWlpKcrKyiI69FaCvq+vT3/oYxH0FHM9duxYeDwe21NQW9HS0oLNmzfD6XQGvRQaGhrgcDgwadKkoEJKFaGdnZ1Bcc5USJgJekqnOszU9u3bdQdTCIEzzzxTf7EcffTRAOxNsUwP1ogRIzBixAgUFBTo9/bgwYNYu3Yt/va3v2HEiBFB/SJeeukljBs3Dt/73vcinsNYUBlDJ9RC0TiMFzndarxeV1cXBgYGwgr6X/ziF7jmmmtwzTXXBC23O7EUndNqplh6sVPIDVUO+/v7g+KW1XNQXiBBT0OtjRw50jTkhgry5uZmtLe3w+v1hjSNGx164PAY8NEK+traWgDWI0CpLyorYWyMoafrot81Goe+oqICJ554ol5ptgq5oRYDut6+vj709PSYCnq7Dj2NSBOpibi+vh6jR4/GzJkzsXHjRvh8Prz11lv4/ve/j3Xr1oWkE9CETVdXl+VkTepv6XK5dCGj5jkgfKdY4zNnJFyn2EgOfXd3NyorK/WKhiro7Tj0a9euxYknnhjijgKHO72aDalHbuLo0aNthdyo4stMqNIIZE899ZRlWiMJemMMfSIc+jFjxkAIETI0ojGGHjgspnp7ezF9+vQQ08Xo+lZUVGDt2rX461//ir///e948cUXgyolZvdJHWChtbU1aK4JOzH0qqA36xSrGkj0TNDIO4RxttidO3dCSokpU6Zg+PDhEELYEr4HDhzAwMCAPsGaGWeccQbeeust03AUOyE3HR0dYR16r9erX1ukkBsAcYXdUBlE7xwiFkFfV1eH3NxcvUyM5NBXV1frDj1dCwDduScouoHMQ0IV9GoetHLoY4UFvQW7du3Swzao08XmzZuRlZWlT8E+ceLEmB16emBnzpwZdedNevECh2dLjDfshsJ/rr76auzcuVPPaPv27UN1dTXKy8vR0dERFCZCY/fSZBvGcbuNgp468wLBoUqqoDdCsY7qbLRWUCFZWVmJrKysoCZvCiE45ZRTcPbZZ+sdbaWUuO+++wAAb7/9dsRCwfhSMYoEVZCrYsDKoVdj+cwE/aFDh/D8889jxowZePvtt4N+Z7NOsSTGjTHcVg69sVOslFJ/EQ0MDOjhYuo5KF2ANkRfbm4u3n//fT1cxMyhJyfU5/PpeY1iL8MJ+lgd+vHjx6OgoMBS0Ksvqry8PGRnZ4cIY7OQG+BwiES4sfTNCuqPP/4YN998c9Byo0NPMbbkDNN9pHOrqCNvWNHZ2anfQxpdx0rQ04gkM2fORE9PD7Zv3457770XL7/8Ml588cWgdKqCns5jhurElpSU6C8z9d4b74OReEJuIsXQd3V16eV8Q0MDOjs7gxz6SIL+yy+/DBEYBJVdZqNt0fMzYcKEoHtnFnLT3d0dVgAD0Ds/rly50jKswdiCqj5bfX196OjoCBL0tDxWSNADoZ1Ywzn0e/fuxaZNm3DXXXcFHc/Y4jJnzhx89NFHKCgowJQpU/Dqq6+ir6/PtqDfvn172In1zGLojSE3RpNENedIGxhD5owhN1ThmzJlCnJyclBeXh4S937TTTfprXwEvQsiCfqmpqaQFjcKj43UKdaOQ28l6NX35Lhx45CdnR31aDQqGzZsQHl5ud4ZlSgvL0dhYaEtQX/33Xfjueeew549ezB69Gi9Im6WV6wcelXQGx36t99+Gzk5Odi5c2dQJUoNB1Lf4V6vF9nZ2UGGAAv6JEAOPaAJ+j179uD999/HpEmT9B9/6tSputCk+G0V4wtFrdFTJjj55JNRX18fdjpyI2pBGY+gX7lypV7ofPTRR6ioqMDpp58Ov9+vH6+pqQlVVVX6S44K5I6ODl2E0v9GYWUU9HSuM844Q/8upQwr6B0OB0aNGmVrbN4DBw5g+PDh+ot7/PjxQYK+oqICY8eOxTnnnINt27Zh27ZtWL9+Pb788ku89NJLyM7ODttkDYQWVNRB1kzQq831xhh6OhYVGi6Xy7RQXb16Nfr7+7F06VIAh51wY0dQoxg3joNuFJyEMeSGjg2EF/QkasvKylBTU4Pdu3dj3Lhxely6UQg1NTXpzZs0f4OVQ28WchMu9MFM0GdnZ+Ooo46ynItADWOhkWrMHHozQU+/UTQhN0ashq10OBxwOBx6vlA77xqxE3Kjim8hRNip5/fu3as79ID2O1F+UytGZoLeKuzG6NDTS1cdxlMIEXfITTSdYtWXdnd3t17OU1luN+RmYGAA9fX1+twLRsjZU6eNJ+h3NTr0Xq9Xj4em+0ZlbriY861bt+LCCy9EZ2cnPv/8c2zfvj2kkhPOoacKBoWrkACNp+VXfYeqz45RTBqvh/qPGGOijS0up59+OgDtfbJw4UIsX74cAHSzK5KgP3jwYERBr5afQOSQG1XQk3NrDLkxOvTbtm1DSUlJUGXKeN8feughrFy5Muh8dHwqV82YM2cO8vPzQ4avpGGCw8XQU2ik0aFX76dZyA2NWGScWOrss8/Gz3/+c/zqV7+KOKLR448/jtdffz1o2a5duzBx4sSQzvFCCIwZMybi8NZdXV2444478J//+Z9Yvnx5UEXIStBnZ2ejoKAA1dXVaGxs1KMjCFXQt7e3o76+HhdddNH/Y+/Lo+MqzuxvtdRqtXbZsuRFXmXZ4A2DwWA2ExxsIGELOBAgDBAChCQMJAbMBEgY9jCTDGACyWQIEwhwAhOWSRz4DfsSdmyD8YYsWZstWXtray3d9fuj+3uqV/3WXiS1XfecPra6X7+ut1XdunW/78PAwIDOL9/a2qr1mSJXMxIrxDHHLVJG6Bljkxhjf2OM9TLGdjLGVpls52eMPckY62aM1THGviN9filjrIExFmCM/YExliV8VsEYe48x1scY+4wxdlgy2k4lu0VCDwBPPfWUTo2hpem6ujqUlpbG5G82Uujppqeb4IgjjkA4HLbMyBEKhXSKZ1VVlTaIlJaWIjs72zWh/5//+R+ceuqpWLlyJYaGhvDBBx/gmGOOiQn2JUJPHTy1U1To6W87hZ46sRNPPBH79u1DX18fmpub0dPTY0roAfvgMnH/FFQERNQv6mA//PBDHH300WCM4etf/zqys7Pxv//7v3j66adRUlKCc845B4cccohtKlKZ0AN6oitnQSEYEfpgMKhb8jUi9O+//z6mTZuGJUuWYOrUqRo5kO1NcodE7SFLgxmhly034netLDfU7pycHG1ySV5Oo99qbm7GsmXL4PP5NELvRKGn43Oi0Msd4erVq/HEE0/gqaeeivmOHOxlZA8zI/T0LLqx3Miw8tCLlhsrQu9EoRevLxCrKImglb+SkhJMmzYN//Ef/4GhoSGcfvrpuhWyRAg9DbriZEpW/WQ4sdwMDQ3F2AqcpK0UFXqatDi13NTX1yMUCmnfN8IJJ5ygyzNPoPNVXl4eY7mhuBd6LqjPtbKo9PX14dRTT0VBQQGOOeYYzJs3D7/+9a91vykTep/PB6/Xi0AgoBF6IpV0DuK13XDOTQk99VFmCj2p0zKhl/ve448/HuvXr8dvf/tbXXaX5cuXAxgRW0SS1tLSosthLxJ6cVIuV3r2eDzIzs7WWcaMCP2+ffs0gk2r0bLlZvLkyWhvb9faRQGxRFRlu5M4qaTVTSBCEjMyMgxX7wg+nw9z586NsX3JIoiR5WZgYABDQ0OuFXqzvu+WW27BqlWr8MADD+C5554zbXNLSwuuvvpqnHHGGTrrkXg/yZgxY4atQi8WAt25c6cu9sCM0Ofn54MxhqlTpyIUCmnpRQnUn3LOtT6Skm6I7o3W1laUl5djwoQJutUXIz4xXhX6hwE0AZgE4AYAf2aMxaZqAG4HUAJgGoBvA/gNY2w+ADDGFgP4NYBzAEyPvm4Vvvs0gFcBTADwnwCeZ4xlIkE0NTVhYGBA66jnzJmjdXAXX3yxtt3SpUvR29uL66+/HgDw0EMP6fYjE3qRNDQ3N8Pj8WDhwoXab5ph9erVKCoqwiuvvKLLQQ+MzE7d5qJ/7LHHkJeXh7179+LZZ5/Fhx9+iBUrVmD69Onwer0aoacc2HIHb6TQGxF6keSQ6kATpOrqai0YzY7QOzk+KipFIIU+HA7jo48+0jIP5OTkYPXq1bjzzjvxm9/8BhdddBEyMzOxePFiR4Re7qhE1cJIoaeMNkaWG9omNzfXkNDTdQGAhQsXaik/RUJNbQBGOlPy7FNAq5l9IBAIxCj0tA8rhV681tTJ0mTKLBh6ypQpqKio0DpWuqesCL3H49Ep1kYwUuiBSIaqQw89FBdddBFee+013WdyoKmZQi/noQfMFXp5yRmwJ/RmHnq6L+T0aiKceOjFFRhgZAAaHh7G22+/rZHJgYEBNDc3a5MzEisqKytx9tlno7q6WiMvbgi9WMiooKBAI6dWJEGGE8sNEOvDt0tbGQqFMDAwgJkzZ8Lj8WjPFokXdgo9iQVmJAMAzjjjDHz22WcxZIPO17Rp03TFncTJi1NCT5+XlZXh5z//ufYbb731lu435Sw3jDEtsJTIs6gSA/ET+ra2NnR3d8cQ+o6Ojpi+Sz4eaovsI5evp9frxT333IOysjKtjywvL9eJELIdq7e3V7NotLe3xyj09AwbPb+5ubloamrSsrhRMLa4Mt/b24spU6YgOzsb+/fvB+fcMCgWGBnzidATiouLdeddrNAu5lunAGCzdK6EefPmxRB6eawuKirS0qMSiODbZbmRFXrZ7kk48sgj8eqrr2LhwoUx/bEIqqAM6Iuu1dTUmE6enRD69957DxMnTsS1114LADjttNO0z6wIPTASywdAq8ILRDz0vb296O7uxvbt27UYQI/Ho8teSLEdcnVmo75t3BF6xlgegLMB/Jxz3sc5fwnAFwDOMtj8uwDu5JwHOOcfAHgRwIXRzy4E8D+c8485510A7gRwSfQ35gNYAOBuznmQc/5I9HhOMGmTjzFWQC8A+UbbASOBLdQZMcbw9ttv4yc/+QnWrFmjbbd06VIAkZzNkydPRk1NjY6QWRF6mrXTbN6M0G/evFkLtHjwwQe1m0TsACjTjRmeeeYZrF27VluiCwaDeOONN3Dbbbfh+OOPxyWXXIJAIIAVK1YgIyMDs2fP1maXsuXGSqEXPWfy8QKRzjwrKwuLFi0CEFkF+Oqrr8AYM122BiLXwQmhlxX6yspK9Pb24rXXXkNXV5culdhtt92Gzs5O9Pf3Y926dQCAJUuW4PPPP7fMGGSkFpop9DRYENEwstyIA4oRod+1a5c26Vu0aJGpQi+TFXniYZbblrIVifug9lop9H19ffB4PPD5fFi3bh1KSkqwevVq099qbm5GaWkp5s6di76+PkyYMEG3egAYW26A2BL1MswIfWFhITZt2oTly5fHqJVGCr1M6M089HSNxNSbRsv1dM6MYKYWp8JyIyr05eXl2LNnD9atW4eVK1dqQbq0PEx2hSuuuAIAsHbtWsydO1er+imfFycKveihJ4j/93q9tkGxVoTerLibXVAs3TcFBQWYNGlSjEJv56GvqakBY8zSw/zNb34TEydOxL/+67/q3g8EAlr6O2BkBUp8bum8kXputhpHffLEiRNx/fXXo7q6Grfffjs+/PBD3aqFrNDTd9ra2jTVUCb0TopLffnll3jooYd0Ci+NoUTAMjMzkZeXh87OTl0qVfF4ZIW+p6dHZ0U1UjMJHo8Hu3btiikwZEToc3NztYmMmeXGiJTm5ORg//79Mel2ZQtXbm4uysrK0NzcjL6+PgwNDek89HROqqqqNGVXHM9ly434bMo2DqP6FDKMCL3cZ9LzaFQJW3xW5fTIAwMD8Pv9yMrK0s6dUeYWESeffLIuNbaMjz/+GEcffTTmzJmjBcIODQ2hoaEhIYW+trYWlZWVuPfee1FfX4+1a9fqjktsO6AXLkTBcfHixdr/yc/f2NiI7du3Y9asWSgsLMSUKVN0AcA0zspZwdJFoa8E0MM5F0OavwCwUNyIMVYMYDKAz022W2Dw2YzohGEBgF2c8wGT78q4GUCX8DINtyblRVySWbBgAf793/9dp1iWlJTgggsuABDxuAGISUdlRug7OjpQVFRkmErujjvuwC9/+UtwzrWsLL/5zW/w8ssva+qmeIOJuei3bNmCVatWaVVu29ra8E//9E947rnntDZ+8MEH6O/vx+rVq7F+/XqEQiGUlZXhhBMic6GKigrs3r0boVAILS0tKC0t1XXw5K0rLCzUKfRi5ygfLzCSGWfy5MnIzc1FVVUVvvrqK0yfPl13nmTMmjULDQ0NtgFa+/bt0xF6IvD33HMPPB6PtjIAAMuWLcMnn3yC+vp6zRO3cOFCBAIBy0IRRg+g2MlR1VcAukwlACwVeiL0okoSDAbR1NSkEYa5c+eitrYWoVBIV+lU3Leo0MvqkpEC3d/frw0KsuVGVujF46QMO1RcrKWlRcsSZETo29raUFJSot23YuVTKhpkpNAD9oReVvxEMMZwyimn6FQuwFiht7PciJNXINZy44bQuw2KFScfBCIUdkGxokJ/9NFHY+fOnXjggQcAQIsZEXOGA8A555yDqqoq/OIXv9Dlxg6FQujp6dHa49ZyQxDbZKfQWxE5wDyPuVlQLG0n3jdlZWUxCr2d5aampgbl5eWWk42CggKsW7cOTz75pG5fNBmUg4rl4HnGmGOFfuLEiWCMYfbs2TjqqKPQ0dGhIzlmhL61tRXNzc3IycnRrgspynYKfV9fH9asWYNrr71Wl4XLaPWCBAt5nJCrdIqWBPGcyRM0GZWVlbr+H0ic0IvXlsY6mvAZ5bmnoSqOLgAAIABJREFUfrG0tBTNzc3a+RMV+jlz5sDr9WL79u2ora1FW1ubzsorW272798Pv9+PuXPn6jKqiBl9rDBv3jzU1tbq+icjyw1gTOjNFPpQKIShoSFkZ2frVijtViePPvpo1NTUmFYT37RpE5YuXYrDDz9cI/R1dXUIh8OWhL65udky9SkFtPr9fp0PHjAWWILBoO5Z2bhxI+6//37dmEiTs+rqamzfvl1T8uXUpETo01KhB5AHQO7hA9H35e0AoNtkO3k/AeF9p79BuAdAofAqN9kONTU1KC0t1Q06ZvjTn/6EhoYGbclPJINyByQT+uLiYmRnZ6OwsFC7+Js3b8Ztt92Gm266Cc8++yyee+45nHfeeRpZuummmzB9+nTdErxI6G+++Wa8/vrr+MlPfoL9+/fj8ccfBwBceeWVeOyxxzAwMID33nsPBQUFWLx4Mb7xjW/gpZdewltvvaVFWhOhb29vRzgcRllZGbKyspCXl4f29nYEg0EMDQ2hqKgI2dnZyMrKQldXF7q7u8EYMw2KpWNmjGHOnDnYvXs3tm3bpqV3M8OsWbPAObcMeuGco6mpSWe5KSsrw9y5c/HGG2/gpJNOirEtLFu2TBcxT2TTqtCXEw89qZw0WBilkKR90TY5OTkxnSodLxH6GTNmYGhoCM3NzTFWDCNCb1d9Tl5md2K5EcmQWfAOWa1opWN4eBjBYBD5+flahyx2qJQnn9rnltD39fUhMzNTt5ogYunSpdi3b59uCV9W6OUJDwXuiX2A1+tFTk6OFtgXCAR0KqNZpV4jGKnKREBlhT4nJ8fw2Dwej21VQVmhX7lypfb/n//859iyZQuam5u1e028LhUVFfB6vTrRgdpFkxs6P0aEXg5+pPOdn5+vy+qQqOXGaHLEOY95VimbhKzQ5+bmoqKiAsFgEBMmTNDueSeWGyv/POHcc8/FwMCAzmZA9588IRKfW3ouSKG3CyIlogmMKIdyddasrCzNhkffaWtr01Ziyb6RlZWF/Px8HbHknOOuu+7SZTx66qmnsG/fPqxfvx5/+tOftIw+mzdvxrRp03RE1ohEA+aWG0BP6O0mdkaQn0tS0MW2ULyC6KE3en6pX6d/jbINiQr9/v37NaVdPA+ZmZmorKzE9u3btRUFUWyithFIVKPUieL7tLpqhenTp4Nzruv/zAi9uEJstDooPqt0XmVCb2a5IRDpNcr+FAwGsX37dhx++OFYunQpNm/erMVjAOb2NhIiSBW/4YYb8Mgjj+i2kTPUiDASBeSV4tNOO01bzSdMnToV2dnZqKqq0hF6WqEhmFlu0kWh7wEgS0oF0ffl7QC9/UXcTt5PgfC+098AAHDOB6K2ngDnPAD9JEKH6upqS1+kCI/Hg2nTpmnKgEzorRR6esgnT56sXfy//OUvKCkpwSmnnILzzz8fDQ0NuOiiizBp0iScfPLJ6O/v11UjBCKEr62tDV999RVefvll3HvvveCc46GHHsIjjzyCtWvX4kc/+hG6u7vx1ltv4b333sOKFSu0jv2MM87QlcKeO3cudu/erR0LqanU+dNDTw+6qLzk5uZq+zUj9PQbVVVV+OKLL3RLWEaga2Flu+np6UFfX1+MQnPJJZcAAC688EKjr+kwZ86cGO+bDCeEnmxUThV6OmeyAkzHSytFYr5dmdDLHZJsuTFS6OVAOLeWGytCT8cH6P2aJ554Ig499FA8+OCDuu+IA0I8lhurTvDwww8HAF0udVmhl+9VmpDIx0j3el9fH8LhsE6hFwMz7QY1I0WIPMKk0ItFx8xgZwuRFfry8nJ85zvfwaOPPorvfve7ACLkq66uDiUlJYbnsbi4GF6vF83NzTGDPBXmMiL0lEZVttzI9qFkeejFwXhoaAic85iVP7PAbrIBko0SMD63bW1t2n1iFaQnorKyEnPmzMEbb7yhvUeTQdk2JT+3fr/fkUIv9h/ASJ8tE3r5+sqEXoSsFO/YsQO33HILLrjgAi3W6M9//jO+9rWv4Re/+AUmTZqk5d/+9NNPsWzZMt3+zBR62XJDVciBEVGEJmhWK7lGMFLoc3JytGMT+196him2AnBP6GWFnoi5XOn50EMPxfbt2/H+++9j+vTpunNvpNBPmjQpJue5U8uN0b1gZrmRC+cB5go9HbeZQm/W91mlot66dStCoRAOP/xwLFiwAJ2dndi/fz/27NkDj8ejjYEyxLGxvb0d//Zv/4ZrrrlGl2GqsbHRlNAbBVDL/M0IHo8HFRUV+OKLL1BTU2Or0MuWG6O+ze2kVdeeuL9pja8A5DHGxHxKiwF8KW7EOe9AJHB2scl22ww+q+Oc90Q/q2SM+Uy+Gzf27NnjmNATyPtqRejFAaKzs1NH6OlB3bRpE5YtW4YNGzZg8eLFuOKKK3DMMccAiASynnXWWfjRj36k+21ScB9++GF4PB5ceeWV+N73voc777wTu3fvxrXXXotFixZh1qxZ+MMf/oA333xTKwtthIqKCgwMDGi52+khoM5fHtQLCws1hV5UA80sN0DED/7OO++gpqbGltBPnz4djDHLTDd03kWFHgBuvfVWtLS04PLLL7f8DSDyIM2YMcOW0FsFxZIdQVT3zBR6IvQ0sMkqSW1tLRhj2vmngUQk9E4tN04UejvLjUyGzAi9mD+btgUiz8jixYuxbds2XZARtS9ey40doZ89ezaysrK0uJDBwUEEg0FLhV4MVhZBpEQe7MyyDInnT4R8roGRFb38/HxN3TarEksw8v6LkBV6IKKqXnXVVZg9ezZycnKwdetW1NbWaveXDMYYysrK0NTUZDjIm1WLlUkDfUeOUaGCZmaIx3IjqocizAK7idCLirtsuRkcHERJSYnmvXWq0AMRAidmvSCFns6NqG7KE3F6Tq0I/YQJE3TKOym3ototT5QBe0Iveuj/+te/IjMzE/n5+Xj88cfR1taG119/Heeddx58Ph8uv/xyPP744+jr63NF6D0eD7xer3bNurq6NFGEzj9N0BJV6GXLjTjhFc+r0YScSKNThd7McgNE0km+8847eOSRR7SsKIQJEyagv79fF1NQWlqqEXp6fpwq9EaE3igoFtBbbqjtZgq9HaE3m4Tn5uZi5syZhgr95s2b4fF4sGjRIm3lfufOnZq9zWwVlsbIuro6vPzyy9r7tAISCATQ3d1tSuiB2HvF6HkxQkVFBf72t7+Bc25I6Pv7+9Hb22touTHq2xhj+MEPfmD7u0ZICaGPEu4XAdweTUt5BiJk+0WDzZ8EcAtjLJ8xdjQigbOUZ+4pAOcyxpYxxgoB3ALgj9Hf2AlgO4D10YDXqwGEAbyTaPudKi8y5PyysqJg5KEHIh0ELXlv2rQJhx9+OObNm4fPP/9c88LTdi+88EJMACmlvHrggQdw/PHHo7i4GHfddRcuvvhiPPDAA1i+fDkYY7jsssvwzDPPoL+/XxcQIoP2T0USaPmWOn9SHaj9YkctqoF+v99UoT/ppJO0DsGO0GdlZWHatGmWCj2dd1mhB+AoEwChsrLSktAbqVwyScjJydH5oO2CYq0I/bRp07TvFRUVIS8vD/X19ejs7ER+fr5GGJ146J1abkSF3orQm5Foep9+Tx7AjWCl0JsVPSHYEXpSduj+MVpKlomxPODJbXFC6Cn9oFmbMjMzY0ioz+fTVce1I/RG11WErNDLbVi4cCG2bt2KXbt26VbpZNAqotG5c0roKUZHDOwDIvdVMoJi5ckREKt2iYO2aLmhYxftf7m5uQgGg1pgJllNNm7ciPb2drS0tDgm9GRjJJCHXg4Il61yRUVFmlhhFRQr2m2AyCSpqKgoRqGXCUpJSQna2trQ1NRkq9C/8847OPnkk3HFFVfgySefxHPPPYdwOIxzzjkHQMTW2dXVhdNOOw3t7e0xolFRUZGh5QbQx6D09PRobRFXLsRjdwo7y41Y9MmI0Iv3nRwwTOdSVuiNLDfyKtuFF16oWREvu+wy3WdyES6y3EyZMgX9/f3o6uoC59yxQk+k30qhLygoQGZmpm6b+vp6TJ482XSVVrwmbiw3QOQ5M1LoN23ahEMOOQQ5OTmoqKiAx+PRCL0VJ8vOzkZZWRnq6urwzjvv4NBDD0VFRYWm0JMVxw2hd6LQA5EVYOIftMJHEzq6TgBMs9wYnad7773X9neNkMq0ldcAmAqgDcC/Azifc97OGLuIMSaq6LcBaAewD8BzAH4YJevgnH8B4CcAXkIkiLUBkUw3hAsBrAbQCeAqAN/inBuXyXMIiqYWA2KdQs7x7NRyQxHaLS0taGxs1CwCTlFWVqYVg7nqqqsARDrLJ554QkvRRJ8tX74cN954o+XDMXv2bDDG8Nprr6G8vFwjddTBy57NeBR6ijnIysrS2m4FOXVlf3+/juTReTci9G5gR+iNyKMYLEoDhuiDlgcHWt6zU+j37Nmjy6BBKUpJoRcHCSceelnJbWlpQX5+vradkYc+EcuNTOitCmaIA4J8jktKSrRO0Qh2hB7Q3z9Gy+AyMRZXFUS4JfRWkDNGUH8hFtNK1HJjpNCLWLRoETZv3owdO3ZYxrKQQu+G0MtBmNOnT8eSJUu0gkCEVFhuiGTJBNbMcnPYYYdh48aN+Od//mdtW7r2tJ0YWE3WEjEntRUqKipQXV2tWbJIoZdXs+QVwKKiIq2/pbFEVobN7pHS0lKdQi9mHCJMnDgRnZ2daGxstCX0u3fvxvz583HZZZdpucJPOukk7Xtz5szBNddcg7fffhvLli3Dcccdp9tfcXGxJvxQfABBJPTd3d1aP55sQk+Em46tvb3dkKAbqcz029Qn2gXF9vf3o66uDnl5eTHKcllZGV566SW8/vrrOv88EJsyVLTcAJGCUlTF14lC7/V6Y/KfyylMPR4PysvLdXFqYgFLgrgaLSr0Rik/ra4VWY5kfPbZZxop9vl8mDNnDrZv347du3fbiqwzZ85EXV0dPvvsMxx55JG6rHDxEHqnCv33vvc9ABFOQ+PQ5MmTMTAwgK6uLl3RNqMsN3bjhBukjNBzzls456dzznM45/M4569G3/8T53yhsF0/5/wiznke53w65/wpaT+Pc86ncc7zOeeXilltOOdVnPPjOOd+zvnhnPPNiba7vr7eMpraCqR2EIwI/fDwMIaGhmIIfWNjo7Y85JbQA5Hg3F/+8pda1h0jlJWV4cMPP8R9991nua/s7GyUl5ejtbVVN7EhhV4m9ERyZDXQykPv9/vx2muvob6+XhcgZwYKogUineb8+fMxb948bf9NTU1agHEiqKysRFVVVUyRGoLRQy4TXcoWYWa5AaDlMBYJPRWxEBV6OSWeU0Jv5KGXldympiadwpMsy41M6M3IsQgrQk/3nRniJfTiMrhMjK0U+q6urqQQevF8Dg8PIxQKuVborSw3g4ODGBwctDzvxx57LDZv3ozm5mZbQt/c3JyQ5QaIKHCvvPKKbrtkFJYCYgPaAGvLjZz69bTTTtPd7/Q+Pcd1dXU45ZRTMGXKFNxxxx2OxQhgxMZIqxOyh97MckPPOGWUAkaeUzpeWUghTJo0yVahp+e/ublZlyAAiDx3RCo555rFaPHixfjhD38Iv9+PX/3qV7rv/PrXv8bf//53vPLKKzGrUzROUHpW8XOfz6dT6M0IfbI89KWlpWhtbUVra2uMQt/f32/YZ69duxbXX3+9Fo9lZ7kBInEHsn+ecMYZZ+BrX/tazPtGhJ4sN0BEvJLjn+xQWloaQ+h9Pp9u7BWTawDmhD5Ryw0QIfTV1dUx1+aTTz7RTQSXLl2K999/H59//rntszZjxgzs3r0bW7ZswbJly3R1WxoaGsAYi7HkiohXoZ8+fTpefvllXbEsuv7Nzc06Qu9UoY8XqVTo0xJ20dRWkAc2I0IPjChvRChmzpyJcDiMjRs3Ii8vzzInuxkWLFiAG264wbG1xA6kPImEUiT0ubm52o1ICr2sBubk5GBgYEBbspbJycknn6xLX2iFBQsW4Msvv9RSedbX16OlpQWPPfYYgJGUlYkef2VlJYLBYIwtgGBG6OXsL1YKPTAygImEngLbRIVeXikyI/RykKWRh14mfp988gmWLFmiOw7AWVCsE8sNdfBuFPpQKITBwUHdOaaJstkkywmhFwcrGijFgVaefJq1ma4PXVsxKBYYGeTcEnqRsIg5oROx3Mj5vo1ANhgAtpYbUugzMjJ058UNofd4PDHPqJO0lfFabuQBWRy0rdKdArGVeOvq6jBz5kysWrUK/f39WLJkiWOCSeMJ3YOk0Msqr5HlRm4jqdsiATYj9HYeejEIWPa8ix76pqYmBINBzWK0YcMGtLa26r4PRCYbp556aowFiI6FLDfyJJMEDkCv0NM9bGahsoOZ5aa0tBShUAi7d++2tNyIv0cTGLovrIJiidBt2bLFMekmiJabgYEBBAIBzXIDRBR6uq5OFHrAmNDL972cy90NoRfPsxPLzaGHHopQKKSLK3n//fcxPDysy8S1YsUKfPDBBwgGgzErGTJmzJiBN998EwMDAzjqqKNQUVGBxsZGBINBNDQ0aNn6zBCvQg8Aa9as0U2IxeJhorW1oKAA3d3dWhxE2ij06Yo9e/Zo1ga3sCP01BE0NTWBc6511vRbL774Ig477DBdcNNYgYicOFsmxUb2bIokRyb0wEjGEDlVoBssXrwY3d3dqK2txbPPPos1a9bghBNO0DJHyDno44Vd6kqjjtBIoTfyFMoKvUzogREVa2hoCI2Nja4VerFTFTvU3NxcDA4OYng44kgbGhrCRx99pLu+oxEUawY6X3TOZIU+HA6b5i12qtC3traip6cnYYVetNzIxXHcWm7ErET0HrWrs7NTq8hsBivLjZPYhXnz5mHx4sU4++yzdZM7GaTQd3V1oaCgQEfK7Sw3dtfGLig2HsuNUc51IPYellVKEbLlhgjOtddeixkzZuDb3/625XGJIDJG1kDqC+WUrfKx0r0gn0PZomJ0jSdNmqSzqhnF/4iB0LICKlpuKK+8GDPgNr1ecXExent70dHRYUjoiUgPDQ1h4sSJ8Hg8KbPckJC0c+dOxx56GTKhHx4exuDgoG7/O3bscC0OijUARDLo9/tRXFyMvXv3ulbo5dUao/5bJPThcNiQ0Is2QTNC70ShX7BgAYDIhIfw7LPPYsqUKdpnwIg1F7B3LojixNFHH62bRFulrBSPLR6F3ghi7Y7W1lZkZ2drsXXhcFj3vCuFPoWoqanRBSK6gVOFntRfenApC0dDQwOOP/74RJqfNNx444247rrrNH8YEOng+/r60NjYqCP0ZEOQyalI6Pv7+xEOhy3VQitQ4Oyrr76Kt956C2vXrsXKlSvx9ttva8vB8ayqyJg9ezays7O1tGwyjGbt1MkRYc7NzbUl9NR5GFUj7ejoQGNjI8LhcAyhnz59OlpbW9HY2OjKciPbYP7whz+gr69P1wnK+zAKig2HwwiFQinz0BupuqRCtbW1YWBgQFtGJVhNLgi00lFbW4v29nZkZWXpfoNWMEg5sfLQd3R0IBAIwOfzxaQiNUsbagQ7hb6zs1OryGwGK8sNrSJYEXrGGLZs2YLnn3/eUkggT2hdXV3MpNxOobdTuZwExTqx3DhR6J3GgQB6y01/fz/279+PGTNm4KijjkJtbS1uuOEGy+MSUVRUBJ/PpxF6cYVWzPBk5KGnbeRjtrPc0MopwajvYoxpKzPy/TphwgR0d3djaGjIsFCUW9CxNDY2xtyTtGIpVhsXJ6vJIPRiXQRxZdiph16GHNAs9l00IQHcnzOq+dLR0aERemovZbpxq9BT/ALBTKFvbGxES0sLrrvuOgSDQS3Dntg2o6BYcYXFybWaOHEi5s2bpxXK3Lp1K/74xz/iRz/6kU4sWL58Oe666y789re/tZ1AnnnmmZg6dSquu+46ZGRkaJPPmpqauAi9G4VeRkFBAXw+H5qamrSUlYwx7Tk1qgydDChCLyHeDDdA/ITe7/drpZ+tPPCjiVWrVuHXv/51TBESIKJeu1Ho+/r6dB11PCgvL8fkyZNx5ZVXwuPx4Oyzz8bKlSvR2tqK7du322bpcAqv14tly5bh/fffj/lMLpRDoE5O7NCdKvTyqgUtjco56AmkmGzatEmnMDvJcgNErkVvby9++tOf4tJLL8Xy5ctj9kHkSg6KFYmT2yw3WVlZpinHgFiFXq5mCURy+a5YsQKLFi3SqU3ypMgIdB737NmjK3BGyM3NRTgc1s4fTULkDr20tBRdXV2oq6uLqXYLxO+hFy0FXq8Xubm56OjocG256erq0kijE4UegCObGlkIdu3aFdOegoICw9UTo8mZEWTLzRtvvKELNHRquYlHobdqm2i5oT7bLLWnHci/u2/fPq2vIIIr9xVGlhv5d2WF3qhflccjM4Ly6aef6goZEei56+joQHV1NcrKymyfMyvQsdTX15sq9OI4IaYNTYaHXqwtYUTojTz0TiaSdB3EmAyPx6OdP6eZkEQUFxejvb1d6+eovVOmTNEIfUFBgWMyKGcKMxrHVq1aBc45Kisr8dBDD+G0007Tqq0TkmW5ASLxO++++y5qamqwevVqzJ8/PyYld0ZGBv7lX/4FV155pe0xejwe1NfXa3EdlNCjpqYG9fX1WipUMyRToWeMaTZFMb2ovOqnFPoUI54c9ITCwkIEAgFwzg0LYZgRegB45plncPPNNzsOshoL0E25Y8eOGIV+YGAAra2thgq9SOjjtdwwxnD55ZeDc46zzz4bEydOxIoVK5CZmYmXXnoJ+/fv1yZFiWLFihX4xz/+gVAopCMqVDjIzEMv2jTsqg6SoiETNkp9Sl5becmTYhs451rubHHfVllugMi1ePHFF9HT04Nbb71VR+acWG7oc7eWGztSaaXQ0722efNmrTiUWCTKiUI/ZcoUeL1e1NbW6oKzCbJfmsierFqTyvPJJ5/oAqySZbmh/qKoqAhNTU1aRWYziMHXTU1NKCoqwrHHHgtgJD2n1fedgpaQt27dGnPuzBR6pxMKkSR0dXXh5JNPxk9+8hPt80Sy3Fh56O0UenHwlYlVPCB1Vb4uYnpfs6BYWawQj8PIkw6MjEcEoyw3QOTeN7pHxODM3bt3x0VMRdB909DQYOqhF+M+RIU+GR56kXCL41C8lhuPx4OsrCzDAnrASHxCPPZdsjvRfUfWGrqHnKasJDgh9BUVFfjud7+LefPm4dNPP8VLL70Usx8nhJ7+tRJwAOD000/Hli1bMGfOHOTk5ODll1+Omx8QxBidjIwMzJw5E1999RWqqqpsM1IlU6EHRlL9NjY2apMJeZxRCn2KUVNTE1fKSiAysIXDYfT29hoqCvQAUREpsRNduHAh7r777qQFtaYCNNFpb2+PUeiBSBpEkeSIJFL2HMeDdevW4f7778d///d/A4g8HEceeSR+97vfAUDSCP03v/lNNDQ0IDMzE+Xl5VpgmJni6FShFzs48s3KhJ5m9bW1tSgrK4vpUMRlwyOOOEL7v1zW3ijLDRDpSN58800sWbIkZoCWCb1RUCx9bkWGMjIy4PP5dEGxdoTbCaF/4YUXtPfcEvqMjAxMnz4de/bs0aWqIxhl5jHaJ3XMH330kS4IKtGgWJmwFBcXa6s0Vgq9SKapENyePXu0TFpAcgg9KfRNTU0xQfsU6CUHLXd3d8Pv99tmsRI99CR2iBUe7VQsoyJdThV6p5YbuWZDPCCFXq7ybGW5oftKVhedKvT9/f3aiptbgkLPSFtbG6qrq+NK1iCCjrehoSHmnpYV+ry8PF1xvmRYbkTCLY6zNNFwS+jFdgOx48Mdd9wBADGBw05AefIptTC1jQj9/v37HdttgFj7lVnc0eOPP46PPvoIRxxxhGFRvKysLIRCIYRCIUuF3uv12sYCnnfeeTj33HPxta99Da+++mpSYuBkLFq0CH/9618RDAYts3gBsfasRBR6YCTVr2j3kccZpdCnED09Pdi3b1/cHRfNLgOBgKGiICv0yRhoRxOURxWAjsyIxyEuDSfTcgNEOrl169bpBuGzzjpLy0xk98A6xcqVK7F+/XqsXr0aPT09eOKJJwDYkwQxa4ad5YZsSoFAwJDQyznoRdA1EBV62r+ZQiz6snfs2IGFCxdCBmMMXq9XZ7kxUuhpNcLKriCqjm4UeqNz7PP5UFlZiVdeeQUVFRU48cQTXRN6YCR1pVOF3qjNYgVL8RkwUujtOmozDz0QeaacEnoaqMU6DTU1Nejs7ITH47E9905QWFioXRMjQs85j/HymynHMsTzQH2jXHDLTikVK40Czjz0dvdwVlYWMjMz0dvbqxF6o+wtTmFG6K0sN5QhjIL1CUQkh4eH0d/fb0roAWgTvngJfXt7u6uquGYQnzm5bxOzfgGp8dDLfcvFF1+MY445RksMIBL6gYEBeDwe28mouBIrK/RHH300OOcx6UCdQFToxVUhIvS1tbWulH9aRRfbatRn2gmKoqgTDAaRkZGBzMzMGIXeyXVijOG5557D66+/HreIaoelS5dqCS7k6uQyxGMwW413gylTpqCxsVFH6JVCP4qgQgdilLUbGBF6cUChz/fs2YOCggJH+dfHExhjGtkTs6OIFqVUEnojfP/73weAGKKfKO655x688sorOP3007VS0mZBfhQUKy7p2hH6oqIi1NfXg3MeQ+j7+vqwdetWU0L/9ttv4+GHH47pNEWyIhMDsfy3VREhUS01CooFIvd3OBy2PN9iGshEFXoAmtd/zZo1qKio0OVLdkvoGxsbY9Qg+j3qaM3aXFBQoJFUO0LvxnIjCwBFRUXYtWsXAGtVuLCwEAMDAxgYGEBNTY32fO3atUsLqE1G1izGmEYq5aVrurfkVK9uCD31K1QARvTQuz2X9P3MzMwYldGoCJwZGGOaStzS0oIJEyYYqpZOYUXozSw3F1xwAf7whz9o1Vjl46X71SwoFhixXhllubECEfra2tqEhC5CQUGBLhE5Bmu8AAAgAElEQVSECCOFXly5MJug2cGK0D/xxBN4//33tXsgKytS2ZkUeiekVCyo5DRmxAkoZWhTU1MMoR8cHMTmzZtdkWBR0KG2xtNOkdCLxF0m9MkkqYlAtDDbxb8Y3SuJKPTz5s3DF198gZaWlhiFXnnoRwGUPcNuJmcGO0Kfk5OD3Nxc7Ny5M0YhTBfQjS4Gy4gPiqgaJNNDb4aJEyeiq6sL999/f1L3S1i0aBF27twJwDwNnxPLDWNMN4ErLi7WSKm4wkGWpY8//tjUQnTYYYfhmmuuiXlfVnzlCURGRgZ27tyJlpYW0wBiudiRkeWGBgWnhD4eD708aaJJ9ve//31dRWZShp0Qx1mzZqG2thY1NTUxaqMYNExtNjs+0c9KSDQoVlYgxawUVtkZaKAOBALYs2cPjjvuOOTk5GDnzp22VWbdgiaYMpGga/P555/rgskTUejp31AohHA47Gq1AzAPaJNrRtiRGlKJW1paErLbAJFnu62tTfNFO7HcZGZm4tJLL42ZlBkRYBmJKvRZWVmYPHkyNm7cCCB+oUsEHbN8Dxl56I2KFSWi0NtNCiiFKAXFOiGlouXGLNVtPKCqrfX19boxlfqcjo4OV4SezjtN7pJB6MVnTCwMlmySmgiOP/54HHHEEfjZz35mK2wY3SuJKPSHHnqoZkOUFXqzytCJIn654QDEl19+idmzZ8f9QIodqFxwhlBWVobq6uq0s9sQXn31VXzwwQe6G11cphOPSyT0PT098Hg8SVEvZCR7kiBi/vz5qK2tRTAYtLTcGAXFytUf5fNEuYTloFiCXSENGbJCL3YUHo8HkyZNwjvvvAPA3J4kqqVmlhvyZltdS1mhd0rozYoh/fSnP8Vxxx2HpUuX4t1339VqOQwMDNiuFhBmzpypkSmZ0MvKidUk4dvf/jbuu+8+nW1Jrt4Zr4detNzQv1bnju79rq4u1NTU4LjjjkN9fT1qa2ttA2rd4u6770Z7e3uM1ausrAwTJkzA+eefDyAyGT3yyCNNCx7JEM8DKfQdHR267D12AXZGGSqMBmM5baWdD5kyrSSL0AORpAJiCju/36/zijsZ4Ol4rVY+EyX0QITE//3vfweQHEsjPaeyQi+mrfR6vcjKykJ2drY2qXWSRtIIVgq9EYigDw0NxU3okzHGzZw5E83NzVrWNYI4Ppit4BqB+oFEFXox05lM6N1abkYDkyZNwqeffupo22Qr9OLzQqKnUuhHEZs2bdLynccDO4UeGAkuS1eFftWqVfjZz35m+D6gJ/d07KTQy+W+0wHz5s0D5xxVVVUJBcXKg4NItERCL9oZ3BJ60U5g1FGUlpbi7bff1tknZMiWGyOFngh9si03wWAQXV1d8Hq9MW33+XxaBcEpU6ZgcHAQ7e3tjnLcE0RFSyYUbhT6e++9F8FgUJfyk9S9ZCn0tMxul2pNtFWQmjdt2jQ0NDQkXaFftGgR3n333Zj7nzGmI7sbNmwAYF7wSIZ4z9HKCxC5b5ymKxTvfSAyIJsp9E6DYgFoFZ+TSei3b9+us0KJlhuziYgMozSPMsR7g1LuuiX0NGmdOXNmUmIx1qxZAyA284tYaE+c6MgrF27tY24UerEd8Sj0ybTc0PlpaGjQnSuR0I+15UYk9HIe+vFiuXGDZCv0NOEqKCjQxluPx6MrJBcMBpWHPpkQ7QUffvihlvItHlBH1NXVZToQkXc3XQm9GTZu3BiTi5oxpg1WZpkYxjvI9rJr1y5ThYfIBM26/X4//H4/QqEQhoaGDBULM0Lv8XiwYcMGHHHEEa6DqawsN0BErejr68PMmTNNO3MnCr0Ty404GDu13AARj79ciVSGWALdzTK3E0Jv56EnGKkq8rJzImkrqcCcUX5wESQiNDc3o6OjA1OnTkV5eXlKCL0VbrjhBnz/+9/HpZdeqili8VhuAoGA1kd2d3c79k6LVhrAnBjLWU/sSE1JSQlaW1uxf//+pBH6L7/8MmYlU/SKO1EFRQIM2Ftu6JjdEhSyn8r5yOPFvffei6qqqphjJEIoTgJFshwMBuNSMt0q9CQsOL0OclBsdnZ2UmLjRPVdJPQ+nw/nnXcezj//fFdZ3WTLjdO4IxkyoRc99IODg+CcjyvLjRskW6H3eDz46KOPUFdXp5uIko0vHA5jYGAgoUlDzG8mbU9RMMaOYox9zhjrY4y9xRgzXRdijFUwxt6LbvsZY+ww4bMTo9/vYYy9KX1vFmOMRz+j10XxtJcCz7Zu3Yqenp6ECL3X60VOTo6lQk/KW7w+/fGKrKwsQ+sLKbXd3d0ptcakCpMmTYLP50N9fb2tQt/T06PlLheLlLgh9ADwwx/+0PEyodwOK4WY7j2rpXORXJkFxcZjuXGi0AORtIh29wkRo3379rki9OXl5bjmmmtw8cUXx5Azn88HxpjOcuN2wJNTt8VjuaHvUIVG0XZiBLp3KM5jypQpGqE3yuaTKnzve9/D7373OyxcuBC1tbXgnMcVFNvd3a1NZN0q9HJQrJ1CbxcUC4wUetu3b19c2UpElJSUIDMzEzt27NBZvmjy6yZVnhPLjd/vR2ZmJgKBgCMya4Tzzz8fd9xxB37/+9+7+p4ZMjMzDYNrnSj08RJ6yljiVKGnbFtOV0rEoNhkWUrFuBnZWvPss8/imWeesbWhicjLy4PH40mZQi9accaT5cYNkq3QA5FVdnl8p+reyfoNEUn10DPGfAD+AuB2AE8CuDX67wkmX3kawN8BrAJwGYDnGWPzOOfDAPoA/A7ANACnG3w3xDlPeA3wo48+wrHHHouNGzciJycHRx55ZEL7o7zQZp0H+RlPPPHEhH4nXUDErqurKy0JPWMM5eXlaGxs1Py2Rgo9LX+LgxHgjNAnaxCwynIDOJtMivYHs6BYp5YbGjycKPTU6dXX17si9ERYnZBvj8eDhx9+2PAzymjiJpBXhttMD0YrKqTkZGdn47e//a1toTk6Vzt27AAATaFvamqC3+8f9VidmTNnoru7Gx0dHXF56Lu7u7Xl6e7ubp0C6HQfgDMPvZNJW2lpKf7xj3+gubnZ1v5kB4/Hg8mTJ6OhoUGnrtIq5vDwMMLhsCuF3orQM8a0tKZmAf12mDBhAm655RZX34kHYh9qpNAnQuiBSH/oRHUV04G6sT4B8aveRiAlfmhoKCmxCx6PBwUFBUnPciMTesq4le6Wm2Qo9GagcSbeSbYVkq3QnwRggHP+e855EMBdAJYxxmJKrzLG5gNYAOBuznmQc/5ItD0nAADn/BPO+Z8A1CW5jTq8//774Jzj6aefxllnnZXwybUj9LR8tmLFioR+J11Ay8mjufyfbEybNg2NjY3o7e3VvNIi8vLyEA6HsX///rgIfbLiCqyCYgHg6quvxn333Ycbb7zRch9mlhs6BrdZbpwotUTM6+rqbAl9dnY2ioqKXCv0diDlBBgdhV5WhOR75Morr7S1Ovh8Pvh8Pi3lLin0nHNUV1e7Kj6TDJCtac+ePY4VenESKSr0biw3ThV6OW2lHamZNGkSqqur484nLoMEI5HQUx/pdDUCGDleSjZg9h0aj1JBHpIJJwp9PORKJJr0jFn1t9QOt9YnILkKPRBR4l944YWkqd1FRUXo6urC0NAQhoeHU6LQDwwMHBCWm1So5wSy3KQDoV8A4HP6g3PeB2A3gNgqNpFtd3HOB4T3vjDZ1ggZjLG9jLFaxth/MMYs707GmI8xVkAvAPkA8OKLL+Laa6/F1q1bccUVVzj8aXPYEfrbb78d27ZtS0u1Oh6ICn06E/qGhga0tbVh4sSJMQMCDUB79+6NIfRUpETu4Oj6k1c6GSCyEg6HDbM0HHLIIbjxxhstK/I5yUPvNijWCTkmQl9bW+vo2aDUlckk9KJC78QmJCPRoNh41aDCwkLs2LEDPp8PxcXFupSkckaaVIPsAbW1tY6DYmWFPhmWGyuFnjIj9ff3O1LoCYkq9ACwZMkSAPqUp0Rc3eRaFxX6/Px8U5JKFULThdAHAoGke+iByPPlxEZDnng3lptUKPSpABUzTCR4VywuaEboDwTLTSoVehob04HQ5wEISO8Fou8nsq2MVgDLAEwHcCKApQDus/nOzQC6hFcDEOmkN2zYgB/84Ac4+eSTHfy0NWRCb5St40Dzz1uBbt50VujJctPa2mqoeBKJb2xs1MionUI/Z84cXHTRRXjyySeT1k4iRqSwx7Ps6SQPfUdHBzwej+X+iaSEQiH09/c7Vuid2jRSRegTVejd5GKWPfTxDoJTp05FS0sLpkyZAsaYLuCXCORoYdKkSfD7/a4U+qysLITDYQwPD6OnpyfuoFg3WW6ckhq5qE+i+PGPf4xzzz0Xq1ev1t4jy00ihN4MNB6Z1XcYLygsLATnHPv27Uu6h56+70R1d+uhF9uYbIU+2SgqKkJnZ2dCfaZVUCwwotCnu+XmoFDoGWPvRoNRjV53AugBIMtrBdH3ZbjZVgfOeQ/n/DPOeYhzXgvgJgDfsvnaPQAKhVc5EPHQV1dX4ze/+Y3dzzoCeRZ7enrg9XpdBa4ciDgQCD1ZblpaWiwJvZFCb0bofT4fnnzySVe5hO1A6iMRm3gGQNFyY6XQ5+bmWi5d03V3OngUFhZq+3Oj0Hd2dupyeicCstyEQiEEg8GEPPTxZLmJVw2i1IKkxovXxU15+GSAMYZZs2bhq6++Qm9vb0xAmBHklZ+ioiItN7uZMCLDbR56saqzFURCnwz7UklJCZ577jldsHJOTg6Ghoa0Z8VtUKwTQj/eFXo6H/X19TEKPdWbGA2FXrTcHGgKPVlukqHQm1luzFak0wE+n0+LY6HnJRXHQePMmBN6zvnxnHNm8roFwDYAWiL3qA2mAsCXBrvbBqAyGkhLWGyyrRNYGpE55wOc8wC9AHQDEYuBnMIuEVAHarcUerBAJPROBvfxiJkzZyIYDGLbtm1JI/SpAJGVeIuwANZBsUTuidBbQSb0duTY4/FoEz43hL6trQ3FxcWu81MbgZSTeBWsRLLcJHKPEKEXA/pfe+01PPjgg2PS/8ycORP/+Mc/AFhXuSWIKUuByPOUn5+fMg+9XATOCnPmzEFRURGuuuqqlJ1LOn6a0LhNW2n1bMmWm/GqIBOhb29v1/WhlAbRqaddRjwKPVlu3HroxzuhT4blRs5mI+ahB9LfcgPo4y2SMa7IGDcKvQO8CcDPGLs8StR/BuBTznmNvCHnfCeA7QDWR/3tVwMIA3gHABhjHsZYNgAvAA9jLJsx5o1+tpwxVskimArgbgAvJflY4kJhYSECgQACgUBa5l1PNmg2ms4eeprwffrpp5aEPhQKjQtCTyQxGZYbUaFnjCErKwsdHR22gwERejeFn2hQd0voJ06caLu9E9Bk3Cq3txUS8dAnYrmhe1LMs3/yySfjxz/+cVz7SxSzZs3Cli1bADgj9HQvNTU1AYg8T3l5ea489G6y3AAj+bid5KHv6OjAo48+ansc8YLa2d7eDiA1Cj0d73iN3RJXLESFHkhM9R1ND326WG5S7aFPZ8sNMHKvpMI/D0Dr21Jh60kqoY8GuJ4D4DoAnYhkrLmYPmeMPcoYE3vGCwGsjm57FYBvRVNWAhFvfD+AP0b30w/gP6OfzQXwf4jYcz5CRO1fl8xjiReiQj9eO8/RRE5ODlpbWzE8PJz2hB4wXnYXB9SxJPQUFJsqyw197iRgVCb0TsixW0Lf19eHPXv2JI3Qk5pJ5Mft/RqP5SYZQbHf+c53sG7dOqxduzau7ycboo3MCaGne6m5uRmAsULv1nJj5aEHnAV2jxaIWLlV6MPhsE7RNoIoMDHGklLtNRUQn7Vk9qHxeugPRMsN9W+JEHqxvsqBGBQLjNwrqbKn5efno6enJyUKfVLz0AMA5/xjAIaRWJzzq6W/qwAcZ7LtmzCx0XDOnwLwVEINTRFky83Bjry8PNTVRTKPpiuhLygowMSJE9HW1mZI6MVO3GwwShbptEIyLTecc4RCoZgYENqnXcEiGiza2toAOCP0X34ZcdstXrzYZsuRXPRbt27F3Llzbbd3AhrwKC1nPISe6kw4VejFojfxDoJ5eXm4//774/puKiBOgJ30gfT8GCn0wWAQjDHbWCRxcgRYV4oFnBVHGy3Ea7kBgNbWVixYsMB0O4rp6urqQn5+fkosBMmA+KyZKfTxFElzS9KIoDudYB9sQbEZGRnIysrSCL1RUOyBQOhTrdB3d3enheXmoIci9HpMnTpVIznpSugBYPr06QCMFXqPx6N1jqQuZ2VlgTGmEfrRWIJMpuVmeDiyUCYr9BQHIVdalSH7op0MHpTd5KSTTrLdlgj9zp07k67QE6F3G/NBKnEoFEI4HHZE6IHISkgiCv14wymnnOJqeyJBZgp9dna2rX9drJIMwFSNHY8KfbyWGyBC6J1absZzDJM4YTMSRZLhoXdC0uKpFEuBu+NdoS8qKkJPT4+2AhlvW8U0qwey5SaVCr1I6JPZ7yddoT/YUVxcjFAohMbGxlHPMDEeIS65j+cBxQ7XX389/va3v5mSTQp0ocGIMaZ1fKPpoU9Wlhsi9LIyOnnyZNTU1NgSeiJpe/fuBeBsMvfee+9hYGDAkYpIhB7AuLPcOF0hEQPMElHoxxsmTJiAE044ARkZGY62N1PoA4GA44mOkeXGykM/ngh9vJYbwJ7QFxYWYnBwEPv370+b/pdW3OxqeTiBrNDbiWy0MuSmUiwALXPSeFfoAaCqqgperzdhQm8UFJvuWW4A5/aseJGfn68ljEh24K0i9EkGEZ3du3dr2ScOZoiEXgzaSzdccskluOSSS2y3EweMsSD0ybLcmCn0RKTtCD0NHnv27IHH43FkuXGT5zs/P1/z6SeT0Hd3d6O9vR0ZGRlxF5ZyukIiBpgNDAykDeFygrfeestxVhhx8uf1euH3+5Gfn4/GxkbHEx2joFgrhZ5WYcYDAUvEctPf32+r0AORCszpcn8ddthhAPTHmKygWLt+i2IOSJCxg0hkx7tCTylYt2/fjpKSkrizNhkp9PQcjeZ4l2yMpkIPAC0tLUn/DWW5STKow7BTTg4WiIR+PAyeqQJZS2gwAvR+zNEKinVDKI0gF6cyI/Rifm4jkN+1pqYGRUVFSU/5xxjT7E/JJPScczQ0NOjy4juF2/Mv5nRO10HQDG7OHZGgvXv3auedAsecKmWiQh8KhTA4OGjroc/Kyoq5v8cCsuXGyXMrPn9lZWWm29GzUV1dnTaEnq6JbLkZjaDYwsJChEIhDA8Pu5pYUeXP8Uzo6T7Ztm1bQn2mSOjp/Hq9XmRkZKCvr++AsNykWqEHIpxBEfpxDlEBUIQ+OeXS0wFELg855BDtvdFW6GkCkUrLDQ1YThX66urquILZnICCrU844YSk7I8IT21tbVzxHm4JvWy5OVA89G5BHvm9e/dq550Cx+Ih9FapLuXiaOMBouXGSbwAoBdKrDIJUVzKrl27xj2h37JlC3bs2KH9LSv08Twf4iqYE0uMeI6cVooFxleQtRloEmhWT8Up/H4/uru7EQqFdNckJydHKfQOIBL6ZPf5Yy9PHGAQZ74qbeXIQ/Ld7353jFuSWuzYsSNG7SNCPzg4OCodnDiBABKz3Jgp9OSLtiO8ROL37NmTMuvZHXfcgZqamqTtnwbzeO0JPp9PN6Fya7lJx0EwGWCMIScnB52dnZp/mgLHnJ4X0XJjlT1CJPTjhXyJpNDpAC+OLRSwbwQx1mS8j0dLluiT4yXDxuHxeOD1erViYk4qVhPcWG5aW1sBjI+YDDMUFBRoE99EFXpaTRKvCVkgRysJRLLhNt4iXpD9NBUKvSL0SYbX60VxcTE6OjqUQh9FR0fHuM1/nCwYdZBjpdAnQujtFHoi6k7IKnXwqcpudMsttyR1f6JCH08qTBosnaYjMyujfjCCCtDRNRDTVrpV6K2qy4qWm/FCvjIzM5GZmYn29va47gErhT43N1ebHI13hV4GXZ/e3t6E+lCxOrAbQu/GckMEd7xMEo3AGENpaSnq6+sTVuhpAiOO66Tch8PhtBQnZIXezlYaL0QPfbIn2cpyk0KkWweaKhQVFY0Lr+poY7QJvd/vRzgc1gqHxGu5sQqK/fGPf4wHHngAK1eutN0Xkf90SVdK7ayvr4+rzRQU65TQy2XU03EQTBaIaNF5z8/P1yqhOiX04XAYw8PDaafQAxEi6EahF2F3r5JKn27jUWZmJrKystDb25vQhNcNoRcJlhuFnuptjJdJohnIR5+oQk8xY6JoSatsQHxjz1hjtD30jY2NSRd9FaFPAchPt3r16jFuicJYYiwUemCkrH0qLDfZ2dm49tprHaXaIkKfKg99slFaWqqdM6tAQzMQcbBSiEXIZdTTcRBMFohcix56IEKUnAys8moHYK3Qt7W1jSvy5ff7XVenpGfQznPf3d0NAJg3b178DRwj0MpNogp9MBhEX1+fK4V+woQJtvuWA5rH0yTRCKQ6p5rQp7vlZjQ89H19fXGNM1Y4+GTTUcCGDRswNDQ0KtVBFcYv/H4/enp6Rm0Jkjqgrq4uMMYc5wEXYWe5cYN0U+gzMjIwe/Zs7Ny5E4sWLXL9fVKJiUA5qUoJjORuPpgtN0S0iFCJy9JO0t3KgzFgfP7pXhxvhD4/Px/Nzc2urBA7d+7UinFZgbb5xje+EXf7xgq5ubkIBAIJ9aE+nw9dXV3gnNsSbpGgOpkApZtCf9VVV6GjowPHH3983Pvw+/2aaCmnaW5paQHgrDL4eIMssKSqP87OzkZGRgZCoVDSbT2K0KcAP/zhD8e6CQrjAH6/Hw0NDQBGZwlSVOh9Pl9cqSLtFHo3kG0U6QA63ngJPTCS49yO0NPnfX19B71CL98rRAgaGhqwdOlS2+87Vej9fr+2cjae1NTJkyejqqrK1QA/d+5cR7Eeb775JhobG1OmOKYSubm5GllOhNCTgm5HuMWVRye/l26E/swzz8SZZ56Z0D7E+0gk7jk5Oaivr495P10gBlD39fWl7HlhjGHSpEloampSCr2CQrrA7/ePqqeQBpdAIBD3kmdWVpaWxxtIjNBTG8yq645H0HWKJ3NOIoR+cHBQKfQYUfzo3/b2dkfWB6cKPRCxGzQ0NIwr8kXpJe3SwcaDZKV1HQvk5uZqZDwRD32qCLdM6MfTJDFVEJ8r2XJjZMVJJ5BtsqenJ6XHMHv27JQQeuWhV1BIEczSe6Xy94CIQp8IoQdGUv8lYrnZsGEDNm3ahOOOOy7ufYw2/vjHP+Kmm26KKwsEDe5E6O0IiHi9gPQMJEsWrrzySgAjdRzEwdQtobeLYaD9OdnvaIEIfaoya6QrcnNzNZIYr2Kal5en2Y6cEPqFCxfi+uuvd7TvdMpykyyYKfRi9psDgdCncpWB7KjjntAzxo5ijH3OGOtjjL3FGJtpsW0FY+y96LafMcYOEz67mjFWzRgLMMbqGGM3S989jTFWxRjrZYy9yBhLj8g7hYMGfr9fU25Gw3YiWm7iJfRE4Ht6enR/x4Pp06c7skuMJyxcuBD33ntvXN8lUtnV1aUt31ohKysLHo/H8QTgQMaZZ56JYDCIU089FYCeKDgh3qLlxk6hp/3NmDEjoTYnEzSwK0KvR25uLvbt2wcg/iw9hYWF2Lt3r7Y/O2zduhW/+tWvHO07IyMDXq8XbW1t8Hq9CfWX6QJ6rnJzc3UWJXEyk46WGyDSh3d2diIcDqd0UkL3crKf96QSesaYD8BfADwIYAKAdwE8afGVpwG8Gt32PwE8zxijNf5XABzFOS8AcDSAixlj34j+Tmn0u9cCmASgM/qbCgrjBiKhGI1ML/EUqJFBxIgCOw9mkukWYo5zv99vG8PAGNMFmB3MCj2gP34xfWCyFXoKFh9PhJ7IULpkhBot5ObmamR8tAi9W2RnZ6O1tXVcWbhSCRpnZMIrEvp0VehF21AqJyUUo5XsVcJkK/QnARjgnP+ecx4EcBeAZYyx2fKGjLH5ABYAuJtzHuScPxJtzwkAwDmv4Zy3CV8JA6iI/v8cAJ9wzjdyzvsA/BzAWsaY6ZocY8zHGCugF4D0vOMU0gYioR9Nhb6lpSXuDlUm9Ac7yXQD0UPv1B4gEno1eRqBuMIUj4feaoWkt7cXwPgi9E7SwB6MyM3N1fqiRAi9uL9kIzs7G+3t7QeF3QYYIe5mhN7j8aRtX1ZQUIDGxkYAqSX0N954I/7v//4PlZWVSd1vsnuRBQA+pz+iZHs3AKMIswUAdnHOB4T3vhC3ZYxdyBjrBrAXQA6AZ01+Zw+AIYwQfiPcDKBLeDU4PSgFhXgwVgp9S0tL3J0RkSCl0LuHaLlxet6UQm8OembcWm7s8o2TnWz69OlJaGVyQIHjy5YtG9uGjDOI1zHeqpqjQegHBwfT1mbiFmQPCwQCuvdF5T6eDGvjAeJqTiqvZ2ZmJr7+9a8nfb/JJvR5AALSe4Ho+6635Zw/xTnPB7AYwOMAuuP4HcI9AAqFl3m9bAWFJED0x40GMabf6O7uVgr9GEAMinWj0CeaxeNAhRtCLyr0dgFt06ZNAzASiDoesHTpUnDOMX/+/LFuyrgCEXDGWNwES1wdTQWhp3ozyQ5wHK849thjASCmBoKZcp9OKCgo0Ah9Oh6HK0LPGHuXMcZNXncC6AEgT6MLou/LcLwt53wrgD4At7n9rrCPAc55gF4YmRwoKKQEogI4GoqFSAjjHfyI0JP6ko4V/8YKYtVIN4SeirHEq0AeqCAi72R1yw2hf/LJJ/HKK6/EVXhNYXRBBLygoCBuWxIp9D6fLyXXnCaIU6ZMSfq+xyPMJi5E6NN5paKgoEBbMU3H43D1hHDOj+ecM5PXLQC2IaKmAwAYYzmI2GC+NNjdNgCV0UBawmKTbYFIznyqoiH/ziwAXkTsPQoK4wLl5aO7COTxeDQCHq+6IGfKUd5e5yBC3tzc7IrQp3vu5lTh0Ucfxdlnn+0o/oQIfTAYtCX0JSUlWCy5LnsAABI/SURBVL16ddLaqZA6iIQ+XhChT1VfNnXqVAAHD6EHgKeffhovvPCC7j0a7yhrVzpCtGcd8ITeAd4E4GeMXR4l6j8D8CnnvEbekHO+E8B2AOujAatXIxL4+g4AMMb+iTFWyiI4AsCPALwe/frzAI5ijJ0anTTcDuBZznl/ko9HQSFuUEc/miBCHm9nRES0s7NT2W1cgryjgUDAsX0mJydHW7pWhF6PZcuW4fnnn3dExIj49fb2pjyHtMLogQhWIlU7aR+kpCcbB5tCDwAXXHABzjrrLN17Rx11FACgqalpLJqUFIgTx3TsQ5JK6KMBrucAuA6RVJInALiYPmeMPcoYe1T4yoUAVke3vQrAtzjnw9HPjkYkSLYbwJ8BPArg4ejv7I9+92EArQCKEUlhqaAwbjAWdhWzlGJuv9/R0aEIvUt4PB7tvLtR6EOhEID0HEDGC7KysuD1ehWhP8CweHFkIb6+vj7ufRChnz07JtleUkDWsFRU+U0nHAgTGiL0WVlZaWk3jb+uuwk45x8DWGLy2dXS31UADMtIcs6vAXCNxe9sBLAx/pYqKBx4IGU4UULf2dmpgjTjQGFhIQKBgCtCD0SUeuXpTgx5eXno6elBT0/PqKSJVUg9qDAdFQuLB5S5K1WEnuKjlD0RWLdu3bjKHuUWROjTtaZA0gm9goLCCD777LNRVQvpt5TlZmxQWFiI+vp6V2krAWW3SQby8vI0hX6041cUUoNk5HZfsmQJzj33XNx6661JaFEs1qxZAwBYuXJlSvafTrj//vvHugkJgVZz0lVcUYReQSGFOPzww0f192bMmIEvv/wyKZabVHlOD2S49fwqQp885ObmoqenB729vWmrsCnE4rPPPkvoevr9fjz33HNJbJEelZWV4JynbP8KowdS6NO1YrMi9AoKBxBImUyksFRGRgY6OzsxZ86cZDbtoABZPSg3tR0UoU8eRMuN8tAfOBhtUUTh4AUp86eddtoYtyQ+KNOXgsIBBMqsk5kZ/1ydSKay3LgHWW3mzZvnaHtF6JMH0XKjCL2CgoJbHHfccbj88stx++23j3VT4oJS6BUUDiAsWLAAABAOh+Peh9/vR09PjwqKjQNtbW0AnBN6ypCRDK/wwQ6y3ChCr6CgEA8KCgrwX//1X2PdjLihCL2CwgGEtWvXgjEWkyPYDZRCHz+oSFRlZaWj7ZcvXw4AqKurS1mbDhbk5eWhubkZfX19itArKCgcdFCWGwWFAwiMMaxduzahFGqK0MePI444AgAwefJkR9sfeeSRAIBt27alrE0HC/Ly8rQJlSL0CgoKBxsUoVdQUNCBCL2y3LjH7373O1RXV2u5qe2Qm5uLpUuX4t57701xyw585ObmoqGhAYC+hLuCgoLCwQBluVFQUNBBKfTxw+/3uy5gs2nTphS15uBCXl4eAoEAgAOjaqWCgoKCGyiFXkFBQQdF6BXSEaLNxqnlSUFBQeFAgSL0CgoKOijLjUI6gooPeTwelJSUjHFrFBQUFEYXitArKCjooBR6hXTEjBkzAERStqZr6XYFBQWFeJF0Qs8YO4ox9jljrI8x9hZjbKbFthWMsfei237GGDvMYJtCxlgTY+xV4b1ZjDHOGOsRXhcl+1gUFA5GKIVeIR1x1FFHjXUTFBQUFMYMSSX0jDEfgL8AeBDABADvAnjS4itPA3g1uu1/AnieMSYH6t4OoMrguyHOeZ7w+lPCB6CgoICqqsjjNn/+/DFuiYKCc0ybNm2sm6CgoKAwZki2Qn8SgAHO+e8550EAdwFYxhiLSfvAGJsPYAGAuznnQc75I9H2nCBsswjAsQAeS7RhjDEfY6yAXgBUrXUFBQMsWbIEAHD++eePcUsUFNxh/fr1ePjhh8e6GQoKCgqjjmQT+gUAPqc/OOd9AHYDWGiy7S7O+YDw3hfStg8C+CkAozr2GYyxvYyxWsbYfzDG7Gqn3wygS3g12B2MgsLBiA0bNqCnpwder3esm6Kg4Ar33HMPrrnmmrFuhoKCgsKoI9mEPg9AQHovEH3f1baMsQsA7Oecv2Pw3VYAywBMB3AigKUA7rNp2z0ACoVXuc32CgoHJTIzM7WMIQoKCgoKCgrjH64IPWPs3WgwqtHrTgA9AAqkrxVE35dhui1jLBcR7/wNRu3gnPdwzj/jnIc457UAbgLwLau2c84HOOcBegHotj1gBQUFBQUFBQUFhXEOV5ViOefHW33OGFsD4Brh7xwAFQC+NNh8G4BKxphPsN0sBvArAJUAZgP4OFpC3Q8gmzG2k3NuFqnnrNa6goKCgoKCgoKCwgGEZFtu3gTgZ4xdHs148zMAn3LOa+QNOec7AWwHsD4asHo1Il75dwBsBTADESvNUgC3AfgQ0YBZxthyxlgli2AqgLsBvJTkY1FQUFBQUFBQUFAY93Cl0NuBcz7AGDsHwH8BeBjAxwAups8ZY49Gt7s6+taFAP4bwHoAOwB8i3M+HP2sSfheF4BBzvn+6FtzESHxkwB0AHgekaBX1wgEZBu/goKCgoKCgoKCwugjXl7KOOdJbkp6gDE2DSrTjYKCgoKCgoKCwvhDOee80enGBzOhZwCmQgXHjnfkIzLxKoe6VgoRqHtCQYa6JxRkqHtCQUY63RP5APZyFyQ9qZabdEL0JDme+SiMDaJB0QDQHc1OpHCQQ90TCjLUPaEgQ90TCjLS7J5w3b5kB8UqKCgoKCgoKCgoKIwiFKFXUFBQUFBQUFBQSGMoQq8w3jGASJGxAbsNFQ4aqHtCQYa6JxRkqHtCQcYBfU8ctEGxCgoKCgoKCgoKCgcClEKvoKCgoKCgoKCgkMZQhF5BQUFBQUFBQUEhjaEIvYKCgoKCgoKCgkIaQxF6BQUFBQUFBQUFhTSGIvQKCgoKCgoKCgoKaQxF6NMIjDEfY+wxxlgdYyzAGPuAMbZC+Hw9Y6yFMdbOGPslE8qiMcYeZYxVMcY4Y+wkab9+xtiTjLHu6L6/Y9OOoxhjnzPG+hhjbzHGZgqfTWeM/ZUx1skYq2GMfdtmX1btup0xVh9t11eMscuj7/8LY6wn+hpgjA0Jfz8a3ea06H57GWMvMsaKhf32SK8wY+ynJu27lDH2WfR81zLG1kuf5zPGfh89552MsT9aHKvh8UQ/Y4yxuxhj+xhjHYyx/2WMTbU6d+MZjLFJjLG/Rc//TsbYquj7j0rnfpAx9oXJPuZHz0MLY6yVMfYX+ZxEr89X0X1tZ4xVOGjb3xljw9J7f4v+Thdj7EPxuVJIDtLwnljFGNscfV63McbOTOT4FWIx3u4JxthJ0fFA/O0ThM9VP5FCpOH9ML76CM65eqXJC0AugNsAzEBkMnYBgFYAeQBOB1APoALAZABfAPie8N2rAZwEYDeAk6T9/hLAywAKABwDoAPAfJM2+KK/cwWAbAB3AXhH+PwtAP8GwAvgcAAtAA6xOCardlUCKIj+fx6AfQAWS9usB/C49F4pgM7oOckB8N8AnjD5/akAhgHMtmjfsdHjmQVgB4ALhc+fBfDv0XPnBXC4xbGaHg+AcwHURa9tFoDHATw91vdcAvfqnwH8V/T8nwmgDcAEg+3+DuAXJvtYDuAyAMXR++5BAP9P+PwbALYAWACAAZgLoNimXWcDeBfAsPT+YgBZ0f+fCWA/oml91evguycAZCDSD/5TdD+nAui125d6pfc9gchYVGXRXtVPqPuBPht3fcSYX0D1SvACAnsBLAPwNIBbhPcvBfCWwfY7EEuc9wE4Xvj7cQC3m/zeGvEGjz54fQBmIzKxCCNKWqOf/yeAOx0cR0y7pM8ro+08U3rfiNBfBeBV4e9ZAIIA/Ab7XQdhQuKgnXcDeCj6/wUAqgFkxHHddMcD4KcA/iR8fjqALWN9f8V5T+YBGARQLrz3JoDLpO0mAxgCUOFwv/MAdAt/fwhglYt2/f92zjXEqiqK47+l46ihZg9LqNREo/IFQUplZVkiOiFFRKmDFdSHvkRRYFH2InpBQRBk2AcpwTKCzArM3kVU9KW0h5WaTmFF+Jqc6LX6sPYdj7dzz71TU3PP9P/Bhnv3Pvt51l13r73X3kOAjcBMqib0mWcsKXwHDu3rsewvoWwyARyRZKA1E/c9BQa7QvllgjoT+sxz0hP/c3loRh0hl5sSY2YTgcOBL4nJ5UeZ5I+BSQ2UcRjxA2k070H1uPt+YnV9EqHkKqG7ikbaUdC+pWb2E7AZ+AbY0EC26jZuIymAnGfbgW43GTObaWa7C8o+C9iUPp8KfAGsNLMfLVygTq9q+7oG+/MMcIKZHW9mQ4HLgPUN9LUZmQh0untHJi5Ppi4D3nf3ryoRFq5cC2uU2z32ZjYQOAWYbOHGtMXMbjELNzMzG2PhAjUmk38psBroqC445VlHGH7rgEfdfU+D/RX1KZVMuPuPxCJJu5kNNLM2oAv4tGfdFgU0q0wcZ2bfJ5eMZamMbqQn/jVKJQ/NqCNa+qpi8c9Ik74ngXvcfY+ZDQP2Zh7ZS1i89ag8s6/BvNX1dD/v7vvM7B3gdgtf8ynAxcD7DbQjF3e/18zuI7bRziUs+HoMI9yC/tLGbISZTSWs9zWZ+t4GRuYVambXEwbUyhR1DDCHcD+6gnCbWWtmE9x9t7vf24P+7CTGaQvwO2GQXNNAX5uRWjJyRFVcO/BYNsLdp+YVaGYTiN2RS1PU0YT+mkPI2UjCAPqacK/aTuY9mtk44BJCmY/Oq8Pd28ysFWgj3NtE71FGmXgaWAEsJ36nC9z959pdFD2k6WSC2CmeRiy4nEjIwE+EW2WlbumJf4cyykNT6Qit0JcQMxtETEK/BO5M0Z2EH3eFESmuHpVnhuflNbNNmcMgY3Lqqa5rEWFpdwCPEKvfHTXKaggP3iP83a9usE9FbazQDjzv7kUr8gCY2SLgOmC+u3el6C5gm7s/7u6/uvtqwo1mRlFZNfpzG7GzcBTxLt7kgOFQNuqOv5lNIvr7VL3C0oGm9cCt7v5qiq68g/uT8bSNUKrzahTzUMpfqGzd/Rd3fxa4MbVR9A6lkgkzOwl4AriQONMyD1hlZsfUa5tomKaTCXff6e6fufsf7v4JcBdwUc5z0hO9T6nkoRl1hCb0JcPMBhBC5MAST45bwCeExVlhCgdcQ2ri7ruI1eHcvO4+yd2HpbC9uh4zO4RwZak8/7W7t7n7KHefARxJWqHPKauntBAHWOpR3cZxxIHV7BbdAGAhMZaFmNkCwiKf6+5bM0kbifeQpfp7Edn+TANWu/sPyWBYAczuQVnNxBfAsCrFVi2P7cALSf5qYmZHEm5Jj7n78kp8yvctB4930djPAh4xs53AB8BAM9tZ8Gc8iHwXLfH3KJtMTAY2uvtb6c/8deLQ+vTibooe0IwykVt8QZr0RO9RNnloPh3RG474Cv9dIA6ZvgEMqYqfn4RpPPm33LQSB8A+J7abhpBO5wMPEKfGhxOry/VuuekArkyfq2+5OZnYOhtM/Pi+AgYX9KeoXVcR218DgHOI7be2qvxFt9zMpcYtN8D5xAGWljrjPZu4SWh6TtogYCtxyn0g4V70AzCyRlk1+wPcAbxCuPS0Ag8Cb/e1vP0DOV1DGCVDgQvI3FaQ+r+D2J4sKmME8CGxmpKXfjfhxzocOJbYHl1c49mj0u9iNHH24bf0uYU40D0/yd5g4FrCBW10X49jfwolk4mJwB7gtPTsTGKlMFcvKvQbmZgFHJc+TyRcH5em79ITkoesPDSdjujzF6jQg5cFYwlrsisJTiWcmdJvIiafu4irKC2T9/WUNxvGpbShwKpU1g4y1zLWaMepSbC7CNeQsZm0G9KPsBN4uZ5w12nXc6msfYSVfnVO/r9M6FP8PMKY2A+speoqKcIV6OGcfGcSB3Mq318j/uiz4/1oJn0qsbrXmZTIGZm0m4GXMt9r9ocwPB4Hvkvvb0NfKoZekNVRwItp/DcD52XSzk1yOign3yZgUfq8JMlDduyz76aVMHD3EEbmskxaxT1sTE4d4zj4RpPxwLuEgbULeAs4u6/HsL+FMslEiltMLDTsI1YPL+/rMexvodlkgrht7BvCT3or4dLaktKkJyQP3fKQ0ptKR1RWQoUQQgghhBAlRD70QgghhBBClBhN6IUQQgghhCgxmtALIYQQQghRYjShF0IIIYQQosRoQi+EEEIIIUSJ0YReCCGEEEKIEqMJvRBCCCGEECVGE3ohhBBCCCFKjCb0QgghhBBClBhN6IUQQgghhCgxmtALIYQQQghRYv4EOxYfLcia9wMAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# 1-bit normalized and whitened waveforms\n", "stw[0].plot()" ] }, { "cell_type": "code", "execution_count": 26, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "correlating stations ('ST01', 'STT1')\n", "2010-09-13T07:26:46.786000Z\n", "2010-09-13T07:26:34.003000Z\n", "# n = 1\n", "...done\n" ] } ], "source": [ "# Cross-correlate\n", "#xcorr_w, stack_w = correlateNoise(stw, ['MLAC','PHL'], 7200)\n", "#xcorr_w, stack_w = correlateNoise(stw, ['ST01','STT1'], 12.75)\n", "xcorr_w, stack_w = correlateNoise(stw, ['ST01','STT1'], 12.786)\n", "\n", "#xcorr_w, stack_w = correlateNoise(stw, ['ST01','STT1'], 100)" ] }, { "cell_type": "code", "execution_count": 27, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEWCAYAAABv+EDhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAIABJREFUeJzt3Xl8E2X+B/DPt3dLKRQoR2kF5C4ICkWOIiAKAiKguMsqyqEuoiiiLIgXivex/nAVV2RlvS9cQUHFA1ERBAUUUEBukKuc5Si0lDaf3x8zE9M0SdMmaZL2+369+moz88wzz0yT+c5zTYQklFJKVT0RwS6AUkqp4NAAoJRSVZQGAKWUqqI0ACilVBWlAUAppaooDQBKKVVFaQAIISLymog86qe8KCLN/JGXP/l6jCKSKyLn+rNMZr47ReRSf+erAsefn5eqSgOAF0Sku4j8ICLHReSoiCwTkU7mulEisjTYZayMRORbEbnJcRnJRJLbg1UmV8IheIjIYBFZIyInROSwiCwWkSYiMtMMqrkiUiAiZx1eLzS3PV9EVovIafP3+Q75thWRL8w8Q2ZSkYjEiMizIrLHPJadIvKcuS7X4ccmInkOr7d4cT5micgmc9tRQT1QH2kAKIWIJAH4BMALAGoBaAhgGoAzwSxXsIlIlDfLVPCZNcE3AEwEUANAEwAvAigiOdYMqokAHgfwvvWaZH8RiQHwMYC3ACQDeB3Ax+ZyADgLYA6AGyv2qEp1D4BMABcCqA6gF4CfAftNhHXMfwC4wmFZc0/nw8x7LYBbrfzCmQaA0rUAAJLvkiwimUfyS5LrRKQ1gJkAupp3CMcAQEQuF5FfzLut3SLykGOGDjWKY+b6Uc47FZHqIvKNiDwvhlgR+aeI/CEiB8w7t3iH9JNEZL+I7BORGzwdkIjUEpFXzbQ5IvKRw7q/i8hWs6YzX0RSHdZRRMaJyBYAWzwsayUiX5l5bBKRv7opR7KIfCIih8xyfCIiaea6xwBcBGCGeW5nOOyvmfl3DRF5w9x+l4jcLyIR5rpRIrLUPGc5IrJDRPq7KoeDTiKywUz/qojEOZR1oBh30MfM/107c/mbAM4BsMAs52QReV1EJprrG1rnyHzd1DwvEZ7yNdelisiH5vHtEJHxDuseEpE55vGfFJH1IpLp5rjOB7CD5Nc0nCT5Ick/SjkfgHHhjALwHMkzJJ8HIAB6AwDJTSRnA1jvRV4QkX+Z7/kTYtQmLvL2mETkAhH52Vz3PoA4lzsxdAIwj+Q+85h3knzDmzKWhuSLJL8GkO+P/IJJA0DpNgMoMj/U/UUk2VpBciOAsQCWm3cINc1VpwCMAFATwOUAbhGRIQAgIo0ALIRRo0iB8eFc47hDEakN4GsAy0iOp/G8jidhBKPzATSDUROZaqbvB+AfAPoAaA6gtOaINwEkAGgDoC6A6WY+vQE8AeCvABoA2AXgPadthwDoDCDD1TIRqQbgKwDvmHn/DcC/RSQDJUUAeBVAIxgX0TwAMwCA5H0Avgdwm3lub3Ox/Qsw7mjPBdATxjkf7bC+M4BNAOoAeBrAbBERD+dlOIDLADSFca7vN8/LBQD+C+BmALUBvAxgvojEkrwexe8inwbwHYwLJ8xybQfQw+H19yRtnvI1A8QCGHebDQFcAmCCiFzmUN5BMP4/NQHMt86dCz8DaCUi00XkYhFJ9HAOnLUBsI7Fnxmz1lxeHithvIdrwXiPfOAYaOHmmMSocXwE471bC8AHAIZ62M8KAHeJyK0icl4p//eqi6T+lPIDoDWA1wDsAVAI441Zz1w3CsDSUrZ/DsB08+97YNyZuEr3GowLwm8AJjksFxhBpanDsq4w7upgbvOkw7oWAAigmYt9NABgA5DsYt1sAE87vE6EUcVvbL4mgN5O2xRbBmAYjAucY5qXATzocIyPujn+8wHkOLz+FsBNLvbXDEAkgAIAGQ7rbgbwrcP/ZavDugRz2/pu9r0TwFiH1wMAbDP/fgnAI07pNwHo6bDtpQ7rmgLIgRHgZprl2mOuex3AXaXlCyN4/eG07h4Ar5p/PwRgkcO6DAB5Ht6DXWA01RyCcef6GoBEpzQPAXjLadkDAN5zWvY2gIecljUDwHJ8tnIAtC/tmGAE0H0AxGH9Dx7eS5EAxgFYBqO5dh+AkW7+75e6yaPE+XBavxTAqLIecyj9aA3ACyQ3khxFMg1AWwCpMC7qLolIZzGabw6JyHEYtYQ65up0ANs87O5yAPEwLhyWFBgXsNVmU8ExAJ+by2GWZ7dD+l0e8k8HcJRkjot1qY7bkswFcATGHahlt/NGTssaAehsldMs63AA9Z03EpEEEXnZbL45AWAJgJoiEumh/JY6AKJR/Fh3OZU12+FYTpt/err7dT6HVvNXIwATnY4p3WF9MSS3wQjY58NoxvoEwD4RaQnj4v6dF/k2ApDqtO5eAPVcHR+A0wDixE0/DMkVJP9KMsUsUw8A93k4F5ZcAElOy2oAOOnFtiWIyD9EZKMYAyqOmXnVcUji7phSAeyleeU1uX2f02iufZFkFozaxGMA/itGs60yaQAoI5K/w7h7amstcpHsHRi1hHSSNWBczK0q6G4Yd4ju/AfGxf0zszkFAA7DaB5pQ7Km+VODRkcVAOyHceGwnOMh/90AaolITRfr9sG48AAAzP3XBrDXIY2r43VcthvAdw7lrEmjaeQWF9tNBNASQGeSSfizmcQ6V55GlRyGUTtp5LDsHKeylpXzOdxn/r0bwGNOx5RA8l0P5fwOwNUAYkjuNV+PhNGRajX5ecp3N4wanuO66iQH+HB8RmHJlQDm4s/3sCfrAbRzakJpBy/b/B2Z7f2TYTQxJtNoMj2OP//fnuwH0NCpHJ7e53Y0+u1ehFHbcNUUWWVpACiFGB2aE+XPzsl0ANfAaGMEgAMA0uTPURGAMergKMl8EbkQwLUO694GcKmI/FVEokSktjgMqzPdBqMpYIGIxJO0wQgM00WkrlmOhg7twXMAjBKRDBFJAPCgu+MhuR9GH8S/xeiEjRYR68L7LoDRYgz7i4UxCuJHkju9PV8w7nZbiMj1Zt7RItLJzZ1XdRiB7ZiI1HJR7gMw2vddHUcRjON+TIwO80YA7oIxWqW8xolImlmW+wC8by7/D4CxZs1ORKSaGB391T2U8zsY/8cl5utvzddLzbKXlu9PAE6KyN0iEi8ikWIMuexU1oMSY9DB3x3eO61gtLWv8LylvdxFAMabfRPjYQS8xWZeYrbhx5iv48z3jivVYTShHgIQJSJTUbJ24c5yc9vx5nvqKhgjfFwSkQki0ss8d1EiMtLc/y9e7s8tMYaYxsEIXNHmMYfltTQsC13BTsJoj/1RRE7B+ND8BuPuFTA+COsBZIvIYXPZrQAeFpGTMDpq51iZ0Rh5McDc/iiMu8H2jjs0q7ljYPQ5fGy+2e4GsBXACrO5ZBGMu2eQXAijSWqxmWZxKcd0PYy7598BHAQwwcxnEYw23w9h3HE1hdGJ6zWSJwH0NbfbB6NK/xQAVxeF52A0dx2GcV4/d1r/LwBXizEq53kX298Oo6llO4z22Hdg9IeU1zsAvjTz2wbgUfOYVgH4O4wOyRwY53iUw3ZPALjfbKr5h7nsOxgXHCsALIXRjGe99pivGSQGwhzBA+McvQKjyaSsjsG44P8qIrkwzvM8GB3jHpEsgNHJP8LMZxSAIeZywKiB5eHPGkEejJsXV74w970ZRvNNPlw3Kborx1Xm/o/C6Gua62GT0wCehfH+OwyjP2Ao/TOH5EsYx9kNwCzz7x4etwhRUrxJTSmlVFWhNQCllKqiNAAopVQVpQFAKaWqKA0ASilVRYX0w7vq1KnDxo0bB7sYSikVNlavXn3YnPBXqpAOAI0bN8aqVauCXQyllAobIuLpSQDFaBOQUkpVURoAlFKqitIAoJRSVZQGAKWUqqI0ACilVBXllwAgIv8VkYMi8pub9SLGVxtuFZF1ItLBH/tVSilVfv6qAbwGoJ+H9f1hfFVhcxhPuXzJT/tVSilVTn4JACSXwHhEqzuDAbxBwwoY3/rUwB/7VkopVT4V1QfQEMWf+70Hxb+6z05ExojIKhFZdejQoQopnFJKVUUh1wlMchbJTJKZKSlezWZWSilVDhUVAPai+PetpsG3725VSinlo4oKAPMBjDBHA3UBcNz8blqllFJB4peHwYnIuwB6AagjIntgfLl3NACQnAngMxjfg7sVxnd1jvbHfpVSSpWfXwIAyWtKWU8YX8pcJoWFheUuk1JKKc9CrhPYkQYApZQKHA0ASilVRWkAUEqpKiqkA8DZs2eDXQSllKq0QjoAaA1AKaUCJ6QDgNYAlFIqcEI6AGgNQCmlAkcDgFJKVVEhHQC0CUgppQInpAOA1gCUUipwQj4A2Gy2YBdDKaUqpZAOAACQl5cX7CIopVSlFPIBID8/P9hFUEqpSkkDgFJKVVEaAJRSqooK+QCgfQBKKRUYIR8AtAaglFKBoQFAKaWqKA0ASilVRfklAIhIPxHZJCJbRWSKi/U1RGSBiKwVkfUi4vWXwmsAUEqpwPA5AIhIJIAXAfQHkAHgGhHJcEo2DsAGku0B9ALwrIjEeJO/BgCllAoMf9QALgSwleR2kgUA3gMw2CkNAVQXEQGQCOAoAK8e9KMBQCmlAsMfAaAhgN0Or/eYyxzNANAawD4AvwK4g6TLh/yIyBgRWSUiqwAdBqqUUoFSUZ3AlwFYAyAVwPkAZohIkquEJGeRzCSZCWgNQCmlAsUfAWAvgHSH12nmMkejAcylYSuAHQBaeZO5BgCllAoMfwSAlQCai0gTs2P3bwDmO6X5A8AlACAi9QC0BLDdm8w1ACilVGBE+ZoByUIRuQ3AFwAiAfyX5HoRGWuunwngEQCvicivAATA3SQPe5O/BgCllAoMnwMAAJD8DMBnTstmOvy9D0DfsuYbERGhAUAppQIkpGcCR0RE6CggpZQKkJAOACKiNQCllAqQkA4A2gSklFKBE9IBQGsASikVOCEdALQGoJRSgaMBQCmlqqiQDgDaBKSUUoET0gFAh4EqpVTghHwA0BqAUkoFRkgHAG0CUkqpwAnpAKA1AKWUChwNAEopVUWFdADQJiCllAqckA4AVg2AZLCLopRSlU7IBwAAOHPmTJBLopRSlU9IBwARAaBfCqOUUoEQ0gHAqgFoAFBKKf8L6QCgNQCllAqckA4AWgNQSqnA8UsAEJF+IrJJRLaKyBQ3aXqJyBoRWS8i33lVOA0ASikVMD5/KbyIRAJ4EUAfAHsArBSR+SQ3OKSpCeDfAPqR/ENE6nqTtwYApZQKHH/UAC4EsJXkdpIFAN4DMNgpzbUA5pL8AwBIHvQmY6sPQJ8IqpRS/uePANAQwG6H13vMZY5aAEgWkW9FZLWIjHCXmYiMEZFVIrLqxIkTALQGoJRSgVBRncBRADoCuBzAZQAeEJEWrhKSnEUyk2RmcnIyAA0ASikVCD73AQDYCyDd4XWauczRHgBHSJ4CcEpElgBoD2Czp4y1D0AppQLHHzWAlQCai0gTEYkB8DcA853SfAygu4hEiUgCgM4ANpZaOA0ASikVMD7XAEgWishtAL4AEAngvyTXi8hYc/1MkhtF5HMA6wDYALxC8rfS8tZOYKWUChx/NAGB5GcAPnNaNtPp9TMAnilLvlYNoKCgwMcSKqWUchbSM4GtGoA+DVQppfwvLAKA1gCUUsr/Qj4AREREaABQSqkACOkAAAAxMTHaBKSUUgEQ8gEgNjZWawBKKRUAIR8AtAaglFKBERYBQGsASinlfyEfAGJjY7UGoJRSARDyAUBrAEopFRghHwC0E1gppQIj5AOAdgIrpVRghHwA0BqAUkoFRsgHAK0BKKVUYIR8ANAagFJKBUbIBwCtASilVGCERQDQGoBSSvlfyAcAnQimlFKBEfIBQGsASikVGCEfALQTWCmlAsMvAUBE+onIJhHZKiJTPKTrJCKFInK1t3lrJ7BSSgWGzwFARCIBvAigP4AMANeISIabdE8B+LIs+WsNQCmlAsMfNYALAWwluZ1kAYD3AAx2ke52AB8COFiWzK0aAEnfS6qUUsrOHwGgIYDdDq/3mMvsRKQhgCsBvFRaZiIyRkRWiciqQ4cOISYmBiRRVFTkh6IqpZSyVFQn8HMA7iZpKy0hyVkkM0lmpqSkIDY2FgC0H0Appfwsyg957AWQ7vA6zVzmKBPAeyICAHUADBCRQpIflZZ5TEwMAKCgoADVqlXzQ3GVUkoB/gkAKwE0F5EmMC78fwNwrWMCkk2sv0XkNQCfeHPxB2CvAWhHsFJK+ZfPAYBkoYjcBuALAJEA/ktyvYiMNdfP9CV/qwagTUBKKeVf/qgBgORnAD5zWubywk9yVFny1hqAUkoFRsjPBNYagFJKBUbYBACtASillH+FfADQYaBKKRUYIR8AtAaglFKBEfIBQDuBlVIqMEI+AGgnsFJKBUbIBwCtASilVGCEfADQGoBSSgVG2AQArQEopZR/hXwA0GGgSikVGCEfALQGoJRSgRHyAUA7gZVSKjBCPgBoJ7BSSgVG2AQArQEopZR/hXwAiIyMRGRkpNYAlFLKz0I+AABGLUBrAEop5V9hEQBiY2O1BqCUUn4WFgFAawBKKeV/YREAYmNjNQAopZSf+SUAiEg/EdkkIltFZIqL9cNFZJ2I/CoiP4hI+7LkHxMTo01ASinlZz4HABGJBPAigP4AMgBcIyIZTsl2AOhJ8jwAjwCYVZZ9aA1AKaX8zx81gAsBbCW5nWQBgPcADHZMQPIHkjnmyxUA0sqyA60BKKWU//kjADQEsNvh9R5zmTs3AljobqWIjBGRVSKy6tChQwC0E1gppQKhQjuBReRiGAHgbndpSM4imUkyMyUlBYAOA1VKqUDwRwDYCyDd4XWauawYEWkH4BUAg0keKcsOtAaglFL+548AsBJAcxFpIiIxAP4GYL5jAhE5B8BcANeT3FzWHWgnsFJK+V+UrxmQLBSR2wB8ASASwH9JrheRseb6mQCmAqgN4N8iAgCFJDO93UdMTAzy8/N9LapSSikHfukDIPkZyRYkm5J8zFw207z4g+RNJJNJnm/+eH3xB4C4uDjtA1AB9c0336Bly5ZYsWJFsIuiVIUJm5nAGgBUoGzduhVDhw7F5s2bceONN+p7TVUZGgBUlXf33XeDJGbMmIENGzZg+vTpwS6SUhVCA4Cq0vLy8vD5559j+PDhGDduHC666CL873//C3axlKoQYREA4uLitBO4HJ5++mn069cPZ8+eDXZRQtaiRYtw+vRpDB5sTF7v0aMH1qxZg9zc3CCXTKnAC4sAYNUASAa7KGHj448/xt13340vvvgCs2fPDnZxQtbHH3+MpKQk9OzZEwCQlZWFoqIi/PTTT0EumVKBFzYBgCQKCwuDXZSwUFRUhJtuugkdO3ZEVlYWHnroIb2jdcFms2HBggUYMGCA/bunu3btChHBsmXLglw6pQIvLAJAXFwcAGgzkJfWrVuHw4cPY+LEiXjyySdx4MABbdd2YcOGDTh48CD69etnX1azZk20adNGA4CqEsIiAMTGxgKAdgR76fvvvwcAdO/eHVlZWahTpw6++eabIJcq9Fjn6aKLLiq2PCsrC8uXL0dRUVEwiqUqgM1mw7x583Ds2LFgFyWoNAB4YdGiRZg6dapXnalLly5Fv379MGHCBOzZs6cCSue6DOeccw7S09MhIujVqxe++eYb7UNxsnTpUqSmpqJJkybFlnfq1AknTpzArl27glQyFUg2mw233HILrrrqKowcObJKfy7CKgD4swlo2rRpGDp0KNauXesxnc1mw9ixY/HII49gyJAhHoPQ2bNncdNNN2H58uWYOXMmhgwZUuHPMCKJpUuXonv37vZlvXr1wu7du7Fjx44KLUuo+/7779G9e3eYjyexa926NQBg48aNwShWWCOJOXPm4MEHH8TLL78cks/weuGFFzBr1ix06dIF8+fPx7vvvhvsIpVLYWEhFi5ciG3btpU/iJEM2Z+OHTuSJN977z0C4Pr16+kvaWlpBMDY2Fju2bPHbbpPPvmEAHjVVVcRAF999VW3af/v//6PAPjJJ59w7ty5BMBJkyb5rcze2LZtGwHwpZdesi9bv349AfCVV16p0LKEsl27dhEAn3/++RLrjhw5QgB85plnglCy8FVUVMQJEyYQgP2nVatWXLduXbCLVsx5553HLl26sLCwkJ06dWKzZs0Cvs/ff/+dvXv35s8//+yX/Gw2G8eOHWs/z7fffrt9HYBV9PIaG/SLvKcfKwDMmzePAPx28o4ePUoAvPrqqwmA7733ntu0l112GVNTU3nmzBmmp6dzyJAhLtMVFhayXr167NOnD202G0ny5ptvJgB++eWX9nRnzpzhsmXLmJ+f75djcfbmm28SANeuXWtfZrPZWLduXQ4fPjwg+wxHb7/9NgHwl19+cbm+bt26vPHGGyu4VOFt5syZBMDx48ezsLCQn376KVNTU1mnTp1i78dgWrt2LQFwxowZJMnnn3+eALh9+/aA7fPYsWNs2bIlAbBr167264Mv/vWvf9kv/MOHD6eI2ANtpQsAn332GQFw+fLlPp84klyyZAkB8OOPP2ZcXBzvuusul+msQPHAAw+QJMeNG8eEhASePn26RNpvv/2WAPj+++/bl506dYqtW7dm/fr1OWfOHE6aNIk1atQgAF5zzTV+eSM4mzx5MmNiYlhQUFBs+dChQ9mkSRO/7y+YioqKePbs2XJtO3bsWCYlJbGwsNDl+p49e7Jbt26+FK9KOX36NFNTU5mVlVXsfb1lyxampaUxNTWV+/fvD2IJDZMnT2ZUVBQPHjxIkty4cSMB8OWXXw7YPu+66y5GRkbabwjfffddn/Lbu3cvq1WrxgEDBrCoqIhHjhxhzZo1OWDAAJKVMAAsXryYAPjNN9/4dOIsL774IgFw9+7d7NatG7Oyslym+/zzzwmAX3/9NUnyyy+/JAAuWLCgRNrbbruNcXFxPHnyZLHl69atszc3iQj/9re/ccyYMQFrkhkwYADbtWtXYvk///lPAgiJD6E/LF++nE2aNGHfvn3LFUjbtm3Lfv36uV1/8803Mzk5OSBBujKymj9dfUbXrl3L+Ph49ujRo9wB2x9sNhvPOeccXn755cWWNWzYkFdffXWp22/atImPPvooP/zwQ545c8br/bZv356XXnopCwsL2aZNG3bu3Llc5beMGDGCMTEx3LZtm33ZY489RgD8/fffK18AWLZsGQHw888/9+nEWcaOHcuaNWvSZrNxwoQJjI+PL3HHTJLTpk2jiPD48eMkjeabpKQk3nTTTcXSFRUVsUGDBrzyyitd7u/s2bNcvHgxN27cSNJoLurevTvT0tL8foE555xzeO2115ZYbp3DuXPn+nV/wbBt2zZGR0ezevXqBMD58+eXaXurjf/RRx91m+a5554jAB44cMDX4lZ6NpuNTZo0Yc+ePd2msZomg9mvsmbNGgLg7Nmziy0fNWoUk5OT3dYGLYMGDbK3uf/73//2ap85OTkUEU6bNo0k+cADDzAiIoJHjx4t1zH8+OOPBMC777672PK9e/dSRPjggw+WKQCE1Sggfw0D/fXXX3HeeedBRHDhhRciLy8P69evL5FuxYoVaNOmDZKSkgAYX0zTv39/LFiwADabzZ5u+fLl2L9/P/7yl7+43F9UVBQuvvhitGrVCgAQGRmJ6667Dnv27MGWLVv8ckwAcOLECfzxxx9o27ZtiXUdOnRAdHQ0li9f7rf9eUJ6Pyph8uTJ6NmzJ9577z2v0s+ZMwdnz57Fzz//jBYtWuDuu+8u05j9H374AQCKjZRyFqiRQEVFRXj33Xfx5ptvBmSewenTp/HII49gzJgxePbZZyvkOVDLly/Hjh07MHr0aLdphg8fjkGDBmHq1KnYvn17wMvkyieffAIAGDBgQLHlffr0QU5ODtasWeN22wMHDuDTTz/FHXfcgebNm+Pjjz/2ap8//PADSNrfa3369IHNZivXvBySmDBhAurVq4f77ruv2LrU1FT07Nmz7COavI0UwfixagC//fZbifb18rLZbKxRowZvvfVWkn+Ompk5c2aJdLVq1SrREWh1Hjr2R0yYMIGxsbH2moI3Nm/eXGK0jq9++OEHj3fEXbp0Yffu3f22P3f27t3LNm3acPDgwTxx4oTHtCtXriQAe9+I1dzmSWZmJi+88EKSf/4/lixZ4nX5Jk+ezOjoaJd9ORZrlJA//z/Hjh1jp06d7HeR7du392uT3NatW9m6dWuKCOvVq0cA7NixI/ft2+e3fbhy6623Mj4+vtT3/+7du5mYmOh2IEWgde7cmZ06dSqxfMeOHQTAF1980e22zz77LAFww4YNnDhxIqOjo736vE+ZMoVRUVE8deoUSbKgoICJiYkcO3Zsmcs/a9YslzUYy8svv2y9typXE9CWLVsIgG+88UaZT5qzP/74o1gVzmazsXbt2iUu9NYF+j//+U+x5Tk5OYyKiuKUKVNIGs0/aWlpHDRoUJnKYbPZmJ6e7lXbo7f+85//eBzRcOeddzIuLq7U9stTp07xscceY6NGjdi0aVNOmjTJ67bbo0ePslWrVkxISGBkZCQ7dOjgsnmNNM5Bjx49mJKSwuzsbNaoUYPXXXedx/x37txJAHzqqadIGhfVyMhI3nfffV6VjyS7du3KLl26eExjs9lYrVo13nHHHV7nW5pJkyYRAF9//XW+//77TEhIYOfOnZmXl+dz3idOnGBGRgZr1arFr776iiQ5d+5cJiQksEuXLgEbdVZQUMDatWtz2LBhXqV/+OGHCYArVqwISHncyc7OLtYU48gaJTdixAi327dr187edv/dd98RAD/44INS95uVlVWizf+KK65g06ZNy1T+xYsXMzo6mpdeeimLiopcpjl8+DCjo6MrXwCwLtrOF+PysNrCFy5caF/Wp08fdujQoVi6N954gwD466+/lsjjkksuYUZGBkmjM7K8wWnkyJGsXbu2239oWd1xxx2sVq2a2/w+/PDDUu+WHe9S+/bty4EDB9r/dnchd/TEE0/Y7+Stu3N3ox6sc/fCCy+QNPpmSruTnD59OgFwy5Yt9mXMzrRpAAAfKElEQVRZWVnMzMwstWykcaGMioriPffcU2rajh07sm/fvl7lW5qtW7cyJiaGo0aNsi+z/h/XXXedT31BNpuNw4YNY0RERIka1AcffEAA9hqvv1kDIz766COv0p88eZJ169blxRdfHJDyuGMN/HA37PeKK65gy5YtXa7bt29fsf6Ls2fPslatWrz++us97rOgoICxsbGcOHFiseXW0NMdO3Z4VfZ33nmH8fHxzMjIYE5Ojse05ue1YgMAgH4ANgHYCmCKi/UC4Hlz/ToAHbzJ1woABw4cKDZ21xfWh87xjWA1CTjeGY8bN46JiYkuO4ZmzJhhHw00evRoRkdHl/qPccUKMv4aI927d2+XVVzLsWPHPF78zpw5w27dujE6OrpYZ/FLL71EAHz66ac97r+oqIhNmjRhr1697K+bNWvGCy+80OUF7vbbb2dsbCyPHTtGklyxYkWpo6MGDhxY4oP68MMPU0R46NAhj+Uj/5zYt2jRolLTDh8+nOnp6aWm88bo0aMZHx/PvXv3Flv+yCOPEAAfeeSRYudo27ZtfOqpp/jkk0+W2iz21ltvEQAfe+wxl+vvuusur5vXyurWW291OzTaHauD/fvvv/dLGVauXMnevXvz3HPP5f33329vbnGUmZnJdu3auQ20jz76KAG4/By/++67BMBVq1bZlw0dOrTUu/gNGzYQAN98881iy9etW0cAfO211zxub7PZ7DWmiy66iNnZ2R7Tk/Ym0YoLAAAiAWwDcC6AGABrAWQ4pRkAYKEZCLoA+NGbvK0AcPz4cQLgP//5z1JPQGmsOwHHtldrprFjUOjYsSN79+7tMo/8/Hy2a9eOCQkJBOB2HkFprBm6b731Vrm2d2Q1ZTmPUHLWs2dPnn/++S7X3XPPPW4nxg0aNIgJCQncuXOn27ytYbOO27/wwgsEwGXLlhVLe/bsWdatW7dYE5g1msRdc1pRURGTk5NLHKM1MuKdd95xWzbLnXfeydjYWK8uWNZFwXlob1kdOnSIsbGxLtt9bTYbr7nmGgJgv379eM8997B79+72fgLrZ/DgwS5HJO3YsYM1atRgt27d3I5iOX36NJs3b84mTZr4fCzOZW/YsKHb0W/unDp1iikpKezfv7/PZVi6dCmjoqJYt25d9u3blyLC/v37F7uZ+/XXXwmA06dPd5vPV199VWLSpuXmm28uMWfECtye+rjmzJnjstZRVFTE2rVrF6sNumKO6OGIESO8qn2TRg2rogNAVwBfOLy+B8A9TmleBnCNw+tNABqUlrcVAPLz8z3e4ZTF1KlTKSLF2rSt9n6rc+X06dOlNhOsWbOGMTExHDJkSKnDx9w5c+YMo6KieO+995Zre0e7d+/2qpb05JNPEkCJjsFly5YxIiKCN9xwg8vtdu7cyfj4eI9v2quuuoopKSnFPnwnT55kUlJSibZ9K1g4D0v9+9//zqSkJJd9DtYH2fnOqbCwkMnJyW7L7qhdu3ZuA7szq7boeOdXHlaz2G+//eZyfWFhIadPn86kpCRGRkYyIyODjz/+OHft2sXjx4/z6aefZlxcHM8991z7UGLSOLft2rVjjRo1ijWJubJkyRKKSLFHBvjK6sB//fXXy7zt448/TsC32f379u1j/fr12axZMx45coTknx2lw4YNs38ub7nllmKTv1zJycmx18SctWzZstjcAZKcP3++yxsbR9aQT1d9PFdeeaXHiZlLly6liHDEiBFlbiKu6ABwNYBXHF5fD2CGU5pPAHR3eP01gEw3+Y0BsArAqnPOOYekcacBgFOnTi3TiXBlzJgxrFu3brFlRUVFrF69OseNG0fSOPkwZwp7kp2dXe6Lv6V169Z+GRVhNW2UVq22psI7NrOcPHmS5557Lhs3buyx/X38+PGMiorirl27Sqw7evQoY2JiOGHChBLrxo0bx9jYWB4+fNi+7Oqrr2ZycnKJzkmrNvbTTz+VyMdqitq6dWuJdUOGDOG5557rtuzkn02Jjz/+uMd0FquG5lyFL4uzZ88yPT3d66DjzooVK5iSksKEhAQ+88wznD9/Pjt06MCIiAiv58fcfvvthB8nVN5zzz2MjIws9n/11rFjx5iUlMS//OUv5dq3zWbjgAEDGB8fX6Kf7qmnniIAXn/99fYJat70gbRp04aXXXZZsWXO7f8Wa5SYp/kAV155pdt+BetRDq4+S/n5+WzatCkbN25c6ig6V8I6ADj+WDUAkoyNjeXkyZPLfDKcDRo0yOVM2Ysuuohdu3Yl+eesWW/a3Hw1dOhQtmjRwud8rOaK0oam2Ww2ZmRksHnz5jxz5gxtNhtHjx5NEeF3333ncds//viDUVFRLu8irTsvV3fLVpvns88+S9L48ERGRrr8f2ZnZxMAn3zyyRLrhg8fzgYNGrhsx7XalV19oCyvvvpqme7oz5w5U+YRRs6sWsS8efPKnYdl9+7d7Nu3r71ZqF69emWa2Jebm8sWLVowJSXF43nyhs1mY9OmTdmnT59y53HvvfdSRIrVarz1+uuvEwD/9a9/uVx/33332c/TRRdd5NXM3dtuu40JCQnF0lr9KytXriyW1mazsWbNmhwzZozb/Jo3b86hQ4e6XGdNSnNVe7KCwxdffFFqmV2pdE1AJJmUlOSXIXmdO3d2+aadNGkSo6OjmZ2dzauvvpqNGzf2eV/esKqJvg7Tu/rqq70eWmY9W+m+++6zt/t7e5G74YYbGBcXV6I9ukePHmzVqpXbTrZu3boxNTWVBw8e5JQpUxgREeG2P6Ft27Yl/kfWNH53d4zWB8rTaKzBgwczPT29TCNuWrZsyauuusrr9M569uzJRo0a+VxTtNhsNm7YsIFff/21vfO8LH7//XfWqFGDLVq04Jo1a0gage6jjz6y34xkZmZy1qxZHsu8atWqUjvsS3PgwAHGx8eXOprG2datW5mUlMTu3bt7bB7ZvXs3P/zwQ68HaFjBeunSpfZl11xzDevWretyPz179nT7WIfTp09TRPjQQw+5XF9UVMR69eqVCBC5ubmsW7cue/XqVe6RYRUdAKIAbAfQBH92ArdxSnM5incC/+RN3o4BICUlpVyTJ5w1atTI5RvOfIYGR40axZiYGN58880+78sb1giD0h6Ze+rUKY4aNYqPPvpoiZEkpOe7DVcuv/xy+x3SqFGjvH6zbdq0iSJSrH/Emqjn6dEKq1evZmxsLFu2bMmoqCiP8x/Gjx/P+Pj4Yh21Vj+Nu8k6Vgexu36A3NxcxsXFlbkNfPDgwWzVqlWZtrFYF0lrzkKoWLJkCRs0aMDIyEi2aNGCiYmJBMAGDRpw6NChbN++PQFw5MiRboPApEmTGBUVZW97L69JkyZRRLh69Wqv0p86dYoXXHABk5OTPQ5IKI/Dhw9TRPjwww+TNJrvkpOTOXLkSJfpx48fz4SEBJfnaPXq1QTA//3vf273d9ttt5WYQGr1F3nqWyhNhQYAY38YAGAzjNFA95nLxgIYa/4tAF401//qTfMPnQJAWlpaqb3mf/zxh8dRDjabjbGxsW6f0d+vXz8CYGJiYsBnT1qsNnlPj6Qm/2z/BsAWLVoU6yQ9efJksTeuN/Lz8/nll19y0aJFZX5A11//+lcmJSXZ236HDh3K6tWrl9oW/PLLL1NEOGrUKI8XDquD+LPPPrMvs0YTuWr/twwZMoSNGzd2Gcysu7vFixeXdnjFPPTQQxQRj++rdevWcfTo0bzzzju5YcMGksaFKiMjg/Xq1Sv3c18C6dChQ5w6dSqHDh3KW265hQsWLLC/D2w2G6dNm0YALgNmQUEB09PT/TKKJycnh3Xq1GHPnj1LvQk5cuQIu3XrxoiICJcPZPSH888/3z5H4fvvvycAzpkzx2Xa2bNnEwA3bdpUYp3VRPX777+73Zc1J8lqBsrJyWFycrL9qZ7lVeEBIFA/jgGgWbNmvOaaa9wetDXN/M4773Sb5tixYx6Hk37xxRcEwCeeeMJtHv6Wl5fHiIgIjx3cRUVFbNmyJTMzM+0Texw7JhcsWOBTm2FZrV27ltHR0ezRo4e90/bBBx/0altvps/n5eWxWrVqxTruvJk9aU2Fd1WbGjZsGGvXrl3mYGedW3ed6wsXLqSIMDEx0ZqFyYEDB7Jjx45uhxWGizvuuMPlBdC6uJX1IXzuWDc37kb5nTp1ijNnzmS9evUYExPDDz/80C/7dWXixImMiYnhrl277I+OdtfUZtXwXAWIf/zjH4yNjfX4frPZbGzUqBF79+5Nm83Ge++91+eRUWQlDQBt2rTx2BZrfbnLJZdc4jaN1czjadz9jz/+6LeZud5q0aKFx+abTz/9lAD49ttvs6ioiO3atWOLFi3sVc/Ro0ezRo0aZXpEra/eeecde42kSZMm5WqP9mTIkCH29vozZ84wMTGRt9xyi8dtrOn+zu2u2dnZjI6OdjlCqTTWKJDnnnuuxLq8vDw2bdqUrVq14pEjR3jgwAE+8MADrFu3Li+44IKAPmO+Ipw5c4adO3dm9erV7UG1sLCQLVu2ZPv27f32JFvHuRBPPPEE8/LyeOrUKc6bN48jRoywPyeqW7duPg/JLY013Llr166Mj4/nFVdc4TZtXl6e20EC/fr1czvfxtHTTz9NAMzKyiIAl0/yLatKGQA6dOhQYiyuxZpBmpCQ4HHmpvWlLdazUkLFkCFD2Lp1a7frR48ezVq1atkng1jNGbNnz2ZBQQFr1apV6jN0AmHevHn84IMPAhJ4XnnlFcKcRGPVzLx53ED37t1LjPKy2lU9Vcc9adCggct+I2ssezjf5Zdm9+7dTE1NZVpaGr/99lv7l5q4axYpr7y8PPvjlkWEIkIATE5O5ogRI/jdd99V2HczWP/Xxo0bl/o48IyMDA4cOLDE8rS0NK8+kzabzT4Q47rrrvPLM5sqZQDo1q2b27t762Jx0003EYDLqeAk+f777xNw/XyfYLr33nsZFRXl9kLaokWLYrNjbTYbu3TpwtTUVPtsw8rwnH9HBw4cYLVq1ditWze2atWKaWlpzM3NLXU7a9z35s2bSRp9HY0bN/bp2TMDBw4sEaDz8/OZkpLic3ttOFizZg2Tk5PtNb7JkycH7GK8ePFi3n///Zw6dSoXLVrk9QxYf8rPz+eDDz7o1Q3DtddeW+Km05pUVpbO/+3bt/vtnFbKAHDxxRe7fZSxNQbeapt092wda3ytN8+MqUjWWGNXX3pvTV5yfjNZk9UAsGbNmm6DXjiz+hcA757dQxp3rLGxsbzyyitps9l4//33E/Dty4QefPDBEh3B1hecVOa7f0cnT57k22+/XWH9TOHCmnTmOKjB+mx++umnQSlTWQJAFMJEbGwscnNzXa7Lzs5GcnIyzjvvPADA5s2b0a5dO5fpIiMjUatWrYCWtawyMjIAABs2bLD/bbG+vCQrK6vY8qysLLz00kvIy8vDoEGDkJCQUDGFrUDDhg1DdnY2AOCSSy7xapu0tDRMmzYNU6ZMwXXXXYf3338fI0eOxGWXXVbucmRmZoIkfvjhB/Tt2xcAMGPGDLRs2dLrcoW7xMREXHvttcEuRshp3749AGDdunXo1asXAOC3334DAJdfzBRqwioA5Ofnu1yXnZ2N+vXro3nz5gCMAODKwYMHUbduXUREhNYXobVs2RIigg0bNpRYt2zZMsTExKBjx44l1o0dO7YiihdUd9xxR5m3mThxIubPn4958+ahZ8+emD59uk9luOSSS1CzZk28+uqr6Nu3LxYvXowff/wRM2bMCLn3kqpYVgBYu3ZtsQBQvXp1pKenB7Fk3gmbd29cXJzbr4S0AkBiYiJSU1M9BoCUlJRAFrNcEhIS0KRJE7cBoFOnToiLiwtCycJTVFQUli5ditzcXHz99ddITk72Kb/4+HiMGDECc+fOxaFDhzBlyhSkp6fjxhtv9FOJVbiqX78+GjZsiGXLltmX/fLLL2jbti1EJIgl807YBIDY2FiPAaBevXoAgBYtWrgNAIcOHULdunUDVkZfZGRklAgAeXl5WLVqVYnmH1U6EfHr3fnf//53FBQUoEePHli5ciWmTZumQVkBAPr27YuvvvoKhYWFyMnJwYoVK9C7d+9gF8srlSYA1K9fH4DnAGA1AYWijIwMbNq0CQUFBfZlq1atwtmzZzUAhIC2bdvihhtuQEJCAm655RZcf/31wS6SChH9+/fHsWPH8NNPP+Hzzz9HUVERBg4cGOxieSVs+gDi4uJc9gHk5uYiNzfXHgAaNWqEI0eO4PTp0yU6Rg8dOhSSTUCA0dFYUFCAtWvXolOnTgBgr1Z269YtmEVTptmzZwe7CCoE9enTB5GRkVi4cCG2b9+OlJQU+2c41IV9DeDAgQMAYA8ADRo0AADs37+/WLr8/HycPHkyZGsAXbt2BfDnqB8AWLp0KVq1aoU6deoEq1hKqVLUrFkTXbp0wVtvvYVPP/0Ul19+OSIjI4NdLK+EXQAwhrn+yRomaAWA1NRUACUDwKFDhwAgZGsAaWlpSE9Px/LlywEANpsNP/zwgzb/KBUGbr75ZuTk5OD48eMYNmxYsIvjtbAJAHFxcbDZbCgsLCy23DkAuKsBHDx4EABCtgYAGE09Vg3g999/R05ODrp37x7kUimlSnP99dcjJycHp0+fRr9+/YJdHK+FTQCIjY0FgBLNQO5qAPv27SuWLtRrAIARAHbv3o29e/figw8+AAD72GKlVGgTEcTHxwe7GGVSKQJARESEvZ28du3aiI6ODssagNUPMHfuXLz00ksYMGAAGjduHNxCKaUqrbAZBWQFAOeRQNnZ2ahbt66900VEUL9+/RI1gHAIAB06dEDXrl0xfvx4ALD/VkqpQAibAGBNunFVA7AmgVlSU1NddgLHxMSgevXqgS2oDyIjIzF//nxkZWUhNjYWffr0CXaRlFKVWNgEAHdNQEePHkXt2rWLLWvQoAG2bNlSbJk1CSzUp2fXqVMHa9euRUFBgT5nRikVUGFzhfEUAJyf7pmamuqyEziUO4AdxcXFISkpKdjFUEpVcj4FABGpJSJficgW83eJp26JSLqIfCMiG0RkvYiU/fGOgL13/fTp08WW5+TklHjYV4MGDZCTk1OsvyCUHwOhlFLB4GsNYAqAr0k2B/C1+dpZIYCJJDMAdAEwTkQyXKTzyGq7d/xOAJJuAwDw5xBRILxqAEopVRF8DQCDAbxu/v06gCHOCUjuJ/mz+fdJABsBNCzrjhITEwEAJ0+etC/Ly8tDQUFBiQDgPBeAJA4cOKA1AKWUcuBrAKhH0hpukw2gnqfEItIYwAUAfvSQZoyIrBKRVdbkLcB1DSAnJwcA3NYArABw8OBB5OXl6Zh6pZRyUOooIBFZBKC+i1X3Ob4gSRGhi3RWPokAPgQwgeQJd+lIzgIwCwAyMzPt+bmqAbgLAI0aNQIA7Nq1CwCwc+dOANAAoJRSDkoNACQvdbdORA6ISAOS+0WkAYCDbtJFw7j4v01ybnkK6qoGcPToUQAoMQooOTkZtWrVsg8F3bFjBwCgSZMm5dm1UkpVSr42Ac0HMNL8eySAj50TiDHwfjaAjST/r7w7iouLQ0REhFc1AABo1qwZtm7dCkBrAEop5YqvAeBJAH1EZAuAS83XEJFUEfnMTJMF4HoAvUVkjfkzoKw7EhFUr17dqz4AAGjevLk9AOzYsQN16tSxNyMppZTycSYwySMALnGxfB+AAebfSwH4ZfptYmKi1wGgWbNmeOedd5Cfn48dO3Zo849SSjkJm5nAgNEP4NwEJCKoUaNGibTNmjUDSezYsQM7d+7U5h+llHISVgHAVQ2gRo0aLp+Z07x5cwDA5s2bsWvXLq0BKKWUk7AKAM41AFfPAbI0a9YMALBkyRIUFBRoAFBKKSdhFQBc1QBctf8DxhfDJCcnY+HChQB0BJBSSjkLuwDg3AfgLgAARjPQxo0bISJo3bp1RRRRKaXCRth8HwAAl8NA09PT3aZ/9tln8csvv+Cyyy6zzw5WSillCKsAUNYaQPfu3dG9e/eKKJpSSoWdsGoCql69Ok6fPo2ioiL7o6DddQIrpZTyLOxqAABw6tQpRERE4OzZsx5rAEoppdwLqwDg+EC4wsJCAK5nASullCpdWAUAx0dCHz9+HABQv76rJ1UrpZQqTdj1AQBGDWD/fuN7aKwvf1FKKVU2YRUAHGsAGgCUUso3YRUAHGsA2dnZEBH9nl+llCqnsAoAzjWAOnXqIDo6OsilUkqp8BRWAcC5D0Cbf5RSqvzCKgA41wA0ACilVPmFZQCw+gA0ACilVPmFVQCIiopCXFwcTpw4oQFAKaV85FMAEJFaIvKViGwxf7udlisikSLyi4h84ss+q1evji1btqCwsFADgFJK+cDXGsAUAF+TbA7ga/O1O3cA2Ojj/tC5c2f7l7zoLGCllCo/XwPAYACvm3+/DmCIq0QikgbgcgCv+Lg/DBs2DGfPngWgk8CUUsoXvgaAeiT3m39nA6jnJt1zACYDsJWWoYiMEZFVIrLq0KFDJdYPGjQIsbGxADQAKKWUL0oNACKySER+c/Ez2DEdSQKgi+0HAjhIcrU3BSI5i2QmycyUlJQS65OSktCvXz8AGgCUUsoXpT4NlOSl7taJyAERaUByv4g0AHDQRbIsAINEZACAOABJIvIWyevKW+iHHnoInTt3RkJCQnmzUEqpKk+MG/dybizyDIAjJJ8UkSkAapGc7CF9LwD/IDnQm/wzMzO5atWqcpdPKaWqGhFZTTLTm7S+9gE8CaCPiGwBcKn5GiKSKiKf+Zi3UkqpAPLpC2FIHgFwiYvl+wAMcLH8WwDf+rJPpZRS/hFWM4GVUkr5jwYApZSqojQAKKVUFaUBQCmlqigNAEopVUVpAFBKqSrKp4lggSYiJwFsCnY5AqQOgMPBLkQA6fGFNz2+8NWSZHVvEvo0D6ACbPJ2Rlu4EZFVlfXYAD2+cKfHF75ExOvHJ2gTkFJKVVEaAJRSqooK9QAwK9gFCKDKfGyAHl+40+MLX14fW0h3AiullAqcUK8BKKWUChANAEopVUWFZAAQkX4isklEtppfNFNpiMh/ReSgiPwW7LIEgoiki8g3IrJBRNaLyB3BLpM/iUiciPwkImvN45sW7DL5m4hEisgvIvJJsMvibyKyU0R+FZE1ZRkuGS5EpKaI/E9EfheRjSLS1WP6UOsDEJFIAJsB9AGwB8BKANeQ3BDUgvmJiPQAkAvgDZJtg10efzO/GrQByZ9FpDqA1QCGVKL/nwCoRjJXRKIBLAVwB8kVQS6a34jIXQAyASR5++194UJEdgLIJFkpJ4GJyOsAvif5iojEAEggecxd+lCsAVwIYCvJ7SQLALwHYHAp24QNkksAHA12OQKF5H6SP5t/nwSwEUDD4JbKf2jINV9Gmz+hdRflAxFJA3A5gFeCXRZVNiJSA0APALMBgGSBp4s/EJoBoCGA3Q6v96ASXUCqEhFpDOACAD8GtyT+ZTaRrAFwEMBXJCvT8T0HYDIAW7ALEiAEsEhEVovImGAXxs+aADgE4FWzCe8VEanmaYNQDACqEhCRRAAfAphA8kSwy+NPJItIng8gDcCFIlIpmvJEZCCAgyRXB7ssAdTd/N/1BzDObJKtLKIAdADwEskLAJwC4LEPNRQDwF4A6Q6v08xlKkyYbeMfAnib5NxglydQzOr1NwD6BbssfpIFYJDZTv4egN4i8lZwi+RfJPeavw8CmAejybmy2ANgj0ON9H8wAoJboRgAVgJoLiJNzE6MvwGYH+QyKS+ZnaSzAWwk+X/BLo+/iUiKiNQ0/46HMVjh9+CWyj9I3kMyjWRjGJ+7xSSvC3Kx/EZEqpkDE2A2jfQFUGlG45HMBrBbRFqaiy4B4HHwRcg9DZRkoYjcBuALAJEA/ktyfZCL5Tci8i6AXgDqiMgeAA+SnB3cUvlVFoDrAfxqtpMDwL0kPwtimfypAYDXzdFqEQDmkKx0wyUrqXoA5hn3KIgC8A7Jz4NbJL+7HcDb5s3zdgCjPSUOuWGgSimlKkYoNgEppZSqABoAlFKqitIAoJRSVZQGAKWUqqI0ACilVBWlAUAppaooDQCqyjAflXurw+tUEflfAPbzkIjsFZGHPaRpaj6SONddGqUCTecBqCrDfDjdJ4F+DLeIPAQgl+Q/vUibSzIxkOVRyh2tAaiq5EkA1p33MyLS2PpiHhEZJSIfichX5peG3CYid5lPVVwhIrXMdE1F5HPzaZLfi0ir0nYqIj3Nfa4x86se4ONUyish9ygIpQJoCoC25tMgrRqBo7YwHl8dB2ArgLtJXiAi0wGMgPGo5FkAxpLcIiKdAfwbQO9S9vsPAONILjOfkprvp+NRyicaAJT60zfml9icFJHjABaYy38F0M68eHcD8IH5PBkAiPUi32UA/k9E3gYwl+QeP5dbqXLRAKDUn844/G1zeG2D8VmJAHDMqkF4i+STIvIpgAEAlonIZSQrxRNEVXjTPgBVlZwEUO72d/OLbXaIyF8A49HXItK+tO1EpCnJX0k+BeNx56X2GyhVETQAqCqD5BEYd+C/icgz5cxmOIAbRWQtgPXw7vuqJ5j7XAfgLICF5dy3Un6lw0CV8jMdBqrChdYAlPK/XABjvJkIBuBAxRVLqeK0BqCUUlWU1gCUUqqK0gCglFJVlAYApZSqojQAKKVUFfX/5P8I4nI/qSEAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Plotting\n", "plotStack2(stw,stack_w,6)" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "anaconda-cloud": {}, "kernelspec": { "display_name": "Python 3", "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.6.4" } }, "nbformat": 4, "nbformat_minor": 1 }