Skip to content

Instantly share code, notes, and snippets.

@christianpanton
Last active September 29, 2022 09:26
Show Gist options
  • Save christianpanton/cd98a4c2e26687d60acbe7694bc7c693 to your computer and use it in GitHub Desktop.
Save christianpanton/cd98a4c2e26687d60acbe7694bc7c693 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "be671641",
"metadata": {},
"source": [
"Closest station is Bornholm with name \"BSD\"\n",
"\n",
"Download event 1\n",
"\n",
"```\n",
"curl -o HHE-1.mseed http://geofon.gfz-potsdam.de/fdsnws/dataselect/1/query?net=DK&sta=BSD&cha=HHE&starttime=2022-09-26T00:02:00Z&endtime=2022-09-26T00:05:00Z\n",
"curl -o HHN-1.mseed http://geofon.gfz-potsdam.de/fdsnws/dataselect/1/query?net=DK&sta=BSD&cha=HHN&starttime=2022-09-26T00:02:00Z&endtime=2022-09-26T00:05:00Z\n",
"curl -o HHZ-1.mseed http://geofon.gfz-potsdam.de/fdsnws/dataselect/1/query?net=DK&sta=BSD&cha=HHZ&starttime=2022-09-26T00:02:00Z&endtime=2022-09-26T00:05:00Z\n",
"```\n",
"\n",
"Download event 2\n",
"\n",
"```\n",
"curl -o HHE-2.mseed http://geofon.gfz-potsdam.de/fdsnws/dataselect/1/query?net=DK&sta=BSD&cha=HHE&starttime=2022-09-26T17:03:00Z&endtime=2022-09-26T17:06:00Z\n",
"curl -o HHN-2.mseed http://geofon.gfz-potsdam.de/fdsnws/dataselect/1/query?net=DK&sta=BSD&cha=HHN&starttime=2022-09-26T17:03:00Z&endtime=2022-09-26T17:06:00Z\n",
"curl -o HHZ-2.mseed http://geofon.gfz-potsdam.de/fdsnws/dataselect/1/query?net=DK&sta=BSD&cha=HHZ&starttime=2022-09-26T17:03:00Z&endtime=2022-09-26T17:06:00Z\n",
"```\n",
"\n",
"Download and build this tool to convert the mSEED files: https://github.com/iris-edu/mseed2ascii\n",
"\n",
"Convert from mSEED to ASCII\n",
"`./mseed2ascii *.mseed`"
]
},
{
"cell_type": "code",
"execution_count": 77,
"id": "3726fba3",
"metadata": {},
"outputs": [],
"source": [
"from matplotlib import pyplot as plt\n",
"import pandas as pd\n",
"import numpy as np\n",
"from scipy.signal import butter, filtfilt\n",
"import dateutil.parser \n",
"import datetime\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 78,
"id": "ccd3702a",
"metadata": {},
"outputs": [],
"source": [
"def sample_start_time(filename):\n",
" header = open(filename).readline().strip().split(\", \")\n",
" return dateutil.parser.parse(header[3])\n",
"\n",
"def samples_aligned(filename, ts):\n",
" samples_per_second = 100\n",
" sts = sample_start_time(filename)\n",
" offset = ts - sts\n",
" skip = int(offset.total_seconds() * samples_per_second)\n",
" # add one to skip, to skip header as well\n",
" return np.loadtxt(filename, skiprows=skip+1)\n",
"\n",
"def make_dataset(files, ts):\n",
" ts = dateutil.parser.parse(ts)\n",
" data_n = samples_aligned(files[0], ts)\n",
" data_e = samples_aligned(files[1], ts)\n",
" data_z = samples_aligned(files[2], ts)\n",
" rows = np.min([np.size(data_n), np.size(data_e), np.size(data_z)])\n",
" data = np.stack((data_n[:rows], data_e[:rows], data_z[:rows])).transpose()\n",
" \n",
" # also generate a nice x axis for the dataset\n",
" # converting to local danish time at the time of incident (UTC + 2h)\n",
" xaxis = np.array([ts + datetime.timedelta(hours=2) + datetime.timedelta(milliseconds=i*10) for i in range(rows)])\n",
" return xaxis, data"
]
},
{
"cell_type": "code",
"execution_count": 79,
"id": "3aa74dac",
"metadata": {},
"outputs": [],
"source": [
"event1_ts = \"2022-09-26T000230\"\n",
"\n",
"event1_files = [\n",
" 'DK.BSD.01.HHN.D.2022-09-26T000159.730000.txt',\n",
" 'DK.BSD.01.HHE.D.2022-09-26T000157.850000.txt',\n",
" 'DK.BSD.01.HHZ.D.2022-09-26T000157.730000.txt']\n",
"\n",
"event2_ts = \"2022-09-26T170300\"\n",
"\n",
"event2_files = [\n",
" 'DK.BSD.01.HHN.D.2022-09-26T170256.870000.txt',\n",
" 'DK.BSD.01.HHE.D.2022-09-26T170258.890000.txt',\n",
" 'DK.BSD.01.HHZ.D.2022-09-26T170257.250000.txt']\n",
" \n",
"x1, event1 = make_dataset(event1_files, event1_ts)\n",
"x2, event2 = make_dataset(event2_files, event2_ts)"
]
},
{
"cell_type": "code",
"execution_count": 81,
"id": "65613414",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"[<matplotlib.lines.Line2D at 0x291aede50>]"
]
},
"execution_count": 81,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAjkAAAGdCAYAAADwjmIIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAACa1klEQVR4nO3dd3wT9f8H8Fd294Iu6ICyyihbStlopWAdKKKioijOL7hQFH6ioKggKg5ExAUuBFFEBQTKXmUVCpRRZlmlLasLOpP7/ZHm8rnLZbVJk6bv5+PRB01yuVyPy9373p/35/ORcRzHgRBCCCHEw8hdvQGEEEIIIc5AQQ4hhBBCPBIFOYQQQgjxSBTkEEIIIcQjUZBDCCGEEI9EQQ4hhBBCPBIFOYQQQgjxSBTkEEIIIcQjKV29AfVBp9MhNzcX/v7+kMlkrt4cQgghhNiA4ziUlJSgWbNmkMvtz8s0iiAnNzcX0dHRrt4MQgghhNTC+fPnERUVZff7GkWQ4+/vD0C/kwICAly8NYQQQgixRXFxMaKjo/nruL0aRZBjaKIKCAigIIcQQghpYGpbakKFx4QQQgjxSBTkEEIIIcQjUZBDCCGEEI9EQQ4hxGMdOF+I15YeQH5xuas3hRDiAo2i8JgQ0jjdM3c7ACCvqBy/PJXo4q0hhNQ3yuQQQjzeiYISV28CIcQFKMghhHiMaq0OryzJxM87zwqe5zgXbRAhxKUoyCGEeIzVh/Pw1/6LeGt5luB5inEIaZwoyCGEeIyL18v436+WVvC/UyaHkMaJghxCiMeo1hmjmUe+28W8QlEOIY0RBTmEEI8hZ4Z+P5ZnLDbWUYxDSKNEQQ4hxGOYm97m2o3K+t0QQohboCCHEEIIIR6JghxCiMfQUYUxIYRBQQ4hxGN0bBbo6k0ghLgRCnIIIR5DJTcW5bQO83PhlhBC3AEFOYQQj8H2opKbKUImhDQeFOQQQjwGx4yHQ93GCSEU5BBCPAYb2OgoyiGk0aMghxDiMdjeVdTTihBCQQ4hxGNwTGBDIQ4hhIIcQojHYJM3lMghhFCQQwjxGIKaHIpyCGn0KMghhHgMCmwIISwKcgghHkNQk0PxDiGNnlODnGnTpkEmkwl+4uPj+dfLy8sxbtw4NGnSBH5+fhgxYgTy8/MF6zh37hxSU1Ph4+ODsLAwTJw4EdXV1c7cbEJIA1WlZYMcinIIaeyUzv6Ajh07Yt26dcYPVBo/8pVXXsHKlSuxdOlSBAYGYvz48bjvvvuwfft2AIBWq0VqaioiIiKwY8cOXLp0CY899hhUKhU++OADZ286IaQBKKvUQi4HNEoFqnU6/vmOzQORW1Tuwi0jhLia04McpVKJiIgIk+eLiorw/fffY9GiRbj11lsBAAsWLED79u2xc+dO9O7dG2vXrsWRI0ewbt06hIeHo2vXrpg+fTreeOMNTJs2DWq12tmbTwhxY+VVWrR/ezUCvVU4MHWIIJMTHqBx4ZYRQtyB02tyTpw4gWbNmiEuLg6PPPIIzp07BwDIyMhAVVUVkpOT+WXj4+MRExOD9PR0AEB6ejoSEhIQHh7OL5OSkoLi4mIcPnzY7GdWVFSguLhY8EMI8Txnr94EABSVVYHjOMEoxzTgMSHEqUFOYmIiFi5ciNWrV2PevHk4c+YM+vfvj5KSEuTl5UGtViMoKEjwnvDwcOTl5QEA8vLyBAGO4XXDa+bMmDEDgYGB/E90dLRj/zBCiFuQMZNwanWcILChmhxCiFObq4YNG8b/3rlzZyQmJiI2Nha///47vL29nfa5kydPxoQJE/jHxcXFFOgQ4uG0HCec1kFnYWFCSKNQr13Ig4KC0LZtW5w8eRIRERGorKxEYWGhYJn8/Hy+hiciIsKkt5XhsVSdj4FGo0FAQIDghxDieX5Kz+F/5zhh9obGzCGE1GuQU1pailOnTiEyMhI9evSASqXC+vXr+dezs7Nx7tw5JCUlAQCSkpJw6NAhFBQU8MukpaUhICAAHTp0qM9NJ4S4oV92nuN/13EctBZqcqj5ipDGx6nNVa+99hruuusuxMbGIjc3F1OnToVCocCoUaMQGBiIsWPHYsKECQgJCUFAQABeeOEFJCUloXfv3gCAIUOGoEOHDhg9ejRmzZqFvLw8TJkyBePGjYNGQz0nCCFG1mpydBygkIEQ0og4Nci5cOECRo0ahatXryI0NBT9+vXDzp07ERoaCgD49NNPIZfLMWLECFRUVCAlJQVfffUV/36FQoEVK1bg+eefR1JSEnx9ffH444/j3XffdeZmE0IaIB0nbKISN1fpOA4KUJRDSGPi1CBn8eLFFl/38vLC3LlzMXfuXLPLxMbGYtWqVY7eNEKIh9HpOMFUDuLmKqrRIaTxobmrCCEeQcdx0FrI5FCMQ0jjQ0EOIcQjiLuQa3WmzVWEkMaFghxCiEfQdyE3PmaneABoBGRCGiMKcgghHmHf2euC7I1WNBogZXIIaXwoyCGEeITnf90nCGSqRakbjkZAJqTRoSCHEOIxKqqNkUyVVhjVaCmTQ0ijQ0EOIcRjXC6p4H+vNqnJoSCHkMaGghxCSIMV5KMSPGYDmSrqXUVIo0dBDiHE7ZRWVOO1pQew4Vi+yWvL919EnxnrkXWxCEM6hAteY+OYalFzVW1jnCO5xfhs3XGUVlTXbgWEEJehIIcQ4nZ+Tj+LPzIu4MmFe01ee3lJJnKLyjF+0T6T1yyNk3PqcikKisvt3pb3Vx3BZ+tO4Iv1J+x+LyHEtSjIIYS4nbNXb/C/v7fiCP7OvAhAmJ0prdBClKwRjZMjfPHhb3eh1wfr7d6W7SevAtB3USeENCxOnbuKEEJqg83IfLftDADgnq7N8dRPxsyOXCY107j5LuR1lVtY5tD1EUKcjzI5hBC3Y65+ZlP2Zf53mcy0W7iwJsdykHOxsMwkSLIkt8j+pi5CiGtRJocQ4nZsScIo5XKLM41X66RH/9PqOLz9dxZ+3XUOAJAzM7XW20kIcW+UySGEuBzHcYJCYQ62ZVh0oiiHXYe5TE6VVscHOIQQz0ZBDiHE5cYt2ocBszairFILACiv0lp9z8XCMhSXVwmeE4yTI65KtvI8IcTzUJBDCKl3FwvL8N6KI7hw/SYAYNWhPFwsLMOzv2Twj61p4qvG1hNXBM+xxcbmCo+t1eoQQjwHBTmEkHr35II9+G7bGTz2/W7B81uOXzbzDlMJUYEmz9nUXGWmVudySQUWbD+DoptVkq8TQhoeCnIIIfUuO78EAHD6yg0rSxpdLa0wea5v6yaCx8JpHaSDGangJ+PsNQyfux3v/HsEE37PtHmbCCHujYIcQohbmb7iiOTzgz/eJHgsg+moxoLiZTOtUheuC8e72Xn6KkbMS8fFmnFw1h8rkHxfgJd7dEYtLq/Csz/vxX+HLvHPbT1xGb0/WI9N2dLbTkhjRUEOIcSltp4QNlF9XzP4n1hxuXDuKJlMBnGyxkzyRuDBb9IFjzfb2EQmk8lsWs7ZvtxwEmsO5+P5X43TWoz+fjfyissxZsEeF24ZIe6HghxCSJ1xHIezV2+YDK7HcRwuXL9pcdC90aK6HHuIZxYXDw4oRbxIjoUmM0Nvr7ooKCnH0r3nbeoxZov8Wsy/RUhjRUEOIaTOPlydjYEfbcKcDScF80vNTjuOfh9uxNKMC075XHFQI26+ssW5azfNvpbHBBRFZVWYvOwQThaUSC57/UYlPlh1FMfyigXP93p/PSb+cRC3iprbasuWQZrLKrU4kS+9nYQ0JhTkEELq7OvNpwDog5r2b6/G7jPXAABzNpwEALz7r3SdTV3IYDoycm2CHCkv/rZfcq6q33afw71f7ZB8T78PN+CbLacx9LOtkq87aloIW/7C2z7ZhNs/3YLVWZesL0yIB6MghxDiUFVazqSHkqOCD5ZMZnnEY1tJveWfA7l4eXGm5PIlotogAMgrKscNG5u2LpdU4JUlmdiTc82ezeTZMt+WIaB67pd9VpYkxLNRkEMIcTpxs5Ja4YhTj8wkqBHX6NjCR62QfP54QQm0FiqZLxWVITtP3yT00Zpsmz+vz8z1+Gv/RYz8Ot36whJoKENCbOcefSIJIR5NnH1QKWSoa02vTGYa1Jgb5dgSP430abDwZhUszQCRNGMDAGDbG4NRKVqwvEoLL5V08FRV1xGXRW+vpmkqCDGLghxCiNMZMi4bjxVgeeZFm5t2LJFLBDknC0rtXo+l7I+55q8bFcYmq4yz101etzXYKq2oNhtkmSOevFQcYEkpq9Si/durAQCnPrgDCrl7dIcnxNmouYoQ4nSGa/4TC/fg78xch63TEaU+loIcc691nLqG//2n9LMmr9tab3PcQg+oLccv46kf96KguBw5V25gw7F8AKa9qyqqrAc5Tyw0dtM/nFtk07bVl6KbVfhl51kU3qy0aXmO4/DlhhNYnWV9fjNCKMghhDhcLUpj7KbTcSaFx7VhqVjZlkJmqUzOwfO2BRI3K8xntB77YTfWHc3HW39nYdDHm/Dkwr3YflI4Ien1G5WoqBYGOVKFyYHeKuZ189tTXqXFHZ9vxZt/HbK43SXlVUievRnv/HvY4nJiHMfh3wO5gm74czedxJTlWXj7b9vWtfvMNXy89jieq5nM1RFqM4bRpaIyFFgZs+hmZTVun70ZU//Oqu2mkTqiIIcQ4nAXJbpfl1aY9kqqCx3H1arQWGznafNZF1sGFwSAvCLh37tkzzkM/GgjzliZm6tMdHFNP3UVaw4LMxTsOnaevor/mAzGrjNXUVwunFBUqvmK/TMMr58sKDG5SK89ko8jl4rx665zqKw2nyFacfASThaUYsH2HLPLSPnnQC5e+G0/kmdv4Z/7Zstp/jVbXL1hW8bHVnM3nkSnqWuQcdb23m4V1VokzdiAXh+sR5WF5sL/DuXhREEpfkw/22C685dXaSXniWuoKMghhNRJ2pF8yedbTFopeNzt3bUO/VwdZ3sQUuvPsDFTVFwmDOByi8px9upNjF9kuQt3RbUxyNHpOIz6diee/TkD55kBCo/nG+uMxEHd5dJKfLbuuOA5Q48vdr3RIT784xsV1SgoLkfy7C3o9cF6wbIbmXm7xv5ofoqIa7UMNL7bKj1lhxR2tOn84nI+Q8VWE4mnBLHV99vO4P55O1BcXoWP1mSjWsfhreW2Z6UKmZnq2RGor5RW4K3lWXxtGHt8vrLkgGAdhy4U4akf9+LUZfvryJyF4zjEv7UaPd5b5zGBDgU5Hqhaq0OLSSuRNGO99YUZZ67cEJx0CbGkSqtDcXkVnv5pr43LOzYg0XGcTXNV1UWJjdknc8HW4Vzh6MfiKRlm/ncM9361HTcqqhH3f6v453/YLh0MrDokzPK8tTwLW08Im7DEmYVDF4sEhdJllVqcv24Moq4zActf+y/yv289cQU/peegxaSVfJajpLwKOh0n6DJvz9hEveNCbFru6KVidHl3LaavOIKf0nOQ+MF6fmDJn3caa6A+Xnvc5L1/Z17ELztN66QMOI7D9BVHsPfsdfy26xz/vD3F2MVlxiCHDcYm/XkIP+88ywe3Cma+M6Vo/f/31yGsO5qP0d/t4rdr+f6LFuu0nG0tc8Oy49RVwWtllVqHNA/XNwpyPERBcTl/cjMMAHapqNzmJoIPVx/D4I83YfKfltvi3VVltQ4bjuULTvDlVVoU3ayy8C5TOh2H/rM2mMx47Uzrj+bjf79mOOQEIq7H+HnnWUz68yD//B8ZF7DyoPm0ufh44TjpuhetjkObN/9D52mOzc7Yg+NqNy6OPb63MfNga68ucRPQhetl2H+u0OSGJM/M6MhSzV/VouBRnFXKzisR9Ga7UanFjpPGCxg7uOGwThGC9xrqZEbMS8f5azeRMG0t/verMDuVddG0/qi8SouXFu/ni6UN1ErhJUccIBmO0zkbTqCyWofvt53ht2F2mj6gYS++UUHeuFRUxr/vcG4RXlqciSnLswTZMNYlZt+yH2/uWMq5cgNzN57EjlNXsPO0/rMnLTOeJ28y+3bdUf3fe6wmm8Y2R7aPDMCKg7kYt2gfblZW41DNfjMM3PjZuhN4eUkmhnxqbMrT6ji0mLQS8W/9B61OHwRZ+v5aotNxSD9l2rzJevZnY50TW+t17UYlen2wzmJ2z11RkONm1hzOs3sk1L0519Drg/Vo8+Z/OHihkP+iAcB/h2z7QszbpB+WfxlzJ+euRn+/Cy0mrRScID9YdRRPLtyLL9afAKC/44x/azW6vLvW7AVDyi+7zuL8tTKcuXIDf+13znxLLI7jMPbHvVh1KA/TV9o39cHV0grBXeSC7WfQcvIqtJi0EtdvVCKvqBxvLc/C4j3nsfn4ZeQWluG1pQcwbtE+FJWZnui+WH8CnaauETQz3TlnG+L+b5WgGcPw2a6m4zinjKTMOiRxAa8Lc9srnmH9Pws9h2KYpicAeDgxRvBY3PzR1F8tyOTcrKzGJ2nGDMiBC4X873GhvmY/1zCK9WpRzdBvu8+ZLLtwRw7+zszFkwv1Wb6Kai2Ky6uglBsvOdtOXDGpSSqv6SkmzliZs/LQJSTN2ICWk1ehWqtD6hfb+NcMTUoXrt/EA1+n8+dCtqnJz8vYfd/c/03qF1vx0ZpsPPztLjz0zU7kFpYJis0/W3ccs9OO45stpwTvu1FRLQgg1Uo5xi/aj5UHL/HnWwOO4/B5zbmL9fC3OwHo98vCHTl4eUkmxi3aZ9IkaUnG2Wu4WFiGpRnnMerbnXj8h91YvPuc2WZmg8slxu/4+qP5KCmvxsbsy9ibc63WgZYrUJBTRzodZ1KZf+1Gpc0HoU7H8dmH/eeu49mfMzDy63TodBw4jsMvO89KXpBY9zMjp9795XbBa7amPr1UxkPBluwPx+nvMFpMWsn3lLhYWIYWk1YihbkTqY2dp69i6t9Z/In5taUH0GLSSly7UYmisio+Pf+/X413HQt35AAwzpXU78ON/GuLdplPXZdVavHh6mP8/yF797cnx7TXDGtTdgFaTFppV1CacfYav98A4DRzZy41V5JBXlG5oGdNzpUb6PHeOrR/ezW/7e8w80P9ezAXRy4ZL9AHzhcJjkmpjMBs5sJXrdVBp+P45pYpy4W9Q2Qy14+zoi88du5nWPvu2cuWMW2sEQ/+F+KrFjx+b+VR4WdW6wTH1g1Rj66Xl2Tit93nkHH2usXu6Oa+D9U6Dsv2XYBOx+FSURlm/ncMe0XL3vPldnSethY3K43nlke/3yV4DOj3t6VpK15besAk22Tw6lJhzcuF6zfxf38dwkuLM7E75xqe/3Uftp64LPjMUokpOgDw39GKaq3JmE7XRV3dN2ZfxhfrT+CDVccEz8/ffEqQNWGbPs9eFWaZxIX6hn1w5JKxuTOdyWCNWbAbm7IL+HNvtVYnCGQrqrXgOP3/y4h56eg7cwPeqMnS7z9XiEnLDuHpn/bi9z3n+b9VfB0zFNNXVusw8Y+D/PP3f52OcYv2CTJ41VqdTdONuEKDCXLmzp2LFi1awMvLC4mJidi9e7f1NzmZVsch7v9WIf6t1Xy2gOM4dJ+ehpTPtmDfOeMXnU35V1br0P6t1Zj050HE/d8qtHnzP9yoqBZM/Dd52SG0nLwKU5Znocs7wiaBymodTuSX2HRQfbv1DK6WVqDFpJVIqBnbQ6fjMG/TKbSYtJL/4rCjs65isj/ZeSV4+NudJk0WYxYY05aGnhJDZm/Wvye/BFodx18opZzIL0GLSSvR8711Jq899M1O/Jh+Fi/8th9llVr8UTODdffpaYKT9ZrD+jsRcQ1ClVZncnH6OT0HT/24F0U3qwRNMEM+24x5m06h8zumzS4ZNSfqBdvPYPvJK+A4Dmeu3MDve86D4zh+H4z8Ot3kZA3oCwvZGamLy6swYp4xIL1cUiHItLHp4YlLD6DztDX89vaesR6PfLcLX23SB3HvrzJeyHaduWayD0rKq5FzxXgivVxazqfQAWNtyNDPtuD+eTtMjqWCkgpcYbI11UzxS9bFIuRctdxrqD7odM5vrnI0R4ztIp7o81KR+eAY0B9X7P89W48D6M9jk5cdwoh5O2o1n9YfGRcw4fcD+HD1MSTN2ICvN58SZJNPFpTwn/+tqPlvjWh/FJdXCeqCpD7LXJZLPP7S87/uw6Jd5wRZl9Hf7xZkV9hsxbG8EhTerMT8zcYsS7spq00+x9bu5meu3sQvzBhKpUzAI775FN9YGjJc93Zrzj/H7lO5TIYxC/YgefYW6HQc3vjzEHq+tw45V27gmy2n0G7KarScvAoTfhcGfmKv/2kMXv7IuAC2bCi/uAILtp9B2yn/Sb7335recJdLKtDvw42Cpi530iBGPF6yZAkmTJiAr7/+GomJifjss8+QkpKC7OxshIWFuWy73mLGPnj2lww8mhiDc0wm4L6vdiBnZqog/f/3uL74Yv0JlFVpsXjPef753/cafweAJaLHqw5dwh0JkQDAH3RP9m2Jt+/qgH6tm2KbaPwMVo+aQKKkohrnr91E/1nGLMeob3diz5vJghTumSs3sPvMNTww33hBjvu/VciZmQpAn+3YfNy0VwN7x3PwQiH+9+s+XCoqh0ymr5/o27oJFj7RCyqFHLfXZHuulFbgi/UnMH5wa1wurUAT5q50w7ECfpRWA6lsxzlR2/venOtQyI3zGh3PL8UXNRmeXh+s44OJ6cM74fw1490Kx3GCO6rs/BIs2XNOkCExEN89d3hbH0CemXEHsvNLBDNRp0++FYt2neOzTAaHLhYKCicNf8e7/x7B0prArsu7a5ExJZlf5uf0s+gd10SQaj6cW4QQH+H25BeXQ85kW85fKxOc3C9eL8PGYwX8xWeLqHj12o1KXGAuhpdLKnC1tAIfr83Gb7uFx6arOKoLeX1is2WOcl1Ud6YfCdr4uFgU8C/aZdq8ZHDggm3Nc52aByDrorCoen5NV3Axtru4mPi7O6SOWWBbfMDcIJy7JgzWu76bZvX9f+67CLVCbjUrd6OiGr4aJZ/BYWuBrpQKs0HimeuLyqrgo1aazSTeYG6q4t9azW/LrDXHbG7qE8stLIO/l4r/zKs3KiTPfQbzt5zGn/su8H9L3hHbywLqk4xz1xwTIzExEbfccgu+/PJLAIBOp0N0dDReeOEFTJo0yer7i4uLERgYiKKiIgQEBDhsu8RdZG0V7KMyOTHZ4q07O2D6CuFBlzMzFf0+3IAL1y3fzRlEBnoJvmxSHkmMwa8WToRi4QEa/Dw20eYT1OjesYIeEqxQfw1/dyU+WQPA1Ls6mHzxvn2sp6CHj0IuQ9twfxytSfWyAY8l/ds0RUW1DrvPWL+bjW3iY5JyBoBZIzoL7o4AoHNUIA5KXDyeGRDHjxECABqlHMemD0XLyasEyz3aOwa/7DT//5HSMRx3dWmG8Yv2888N6xSB4vIqbK8pMI0L9YVaIeeDmocTY5B+6irfbHVrfBg2MHU3C564BTtOXjG583YHXaODkHm+EN1jgnCioFRyVnB31TsuxOK4PLUxuF0oNmYbbzo6RAYImjneGBqPD1cfk3qrS/lrlHh9WDzeWu6ZA+Wldo7E9RuVJr2UbNWiiQ9yJM4xgP5cIR4E0hHYc27Lpr5Wx3kSO/7eMJPi8rqq6/Xb7ZurKisrkZGRgeRk492sXC5HcnIy0tOlZ/GtqKhAcXGx4Med1CbAAWAS4AD6QMvWAAeA1QAHgF0BDqBPa9pzB2YuwAGE6WOpuGS5xJQAn6wVzgCt1XGC9dhanLr1xBWbi5SlAhwAJgEOAJMAR1NzEjAEOG3D/SCX6ZsV2O02sBTgAPp2e0OGS6XQZ2+ulFYIanAuXC8TbPOiXecEJzC2aRUAPl6TbbHZwJX+N6gVgJppHRpYl1ZHBzgATL7/4qbTGw4ehNFRSiqqUejggf3cycn80jrVdJkLcAA4JcABhOdcewMcQN/hw924fZBz5coVaLVahIeHC54PDw9HXp50Wm7GjBkIDAzkf6Kjo52ybak1zUfEuZoFevG/HzhfaPK6ITvRiukZcqWWvX/E6fO6+L874k2eW/JMb7SL8Bc8F+yjRpi//m8UD85mSXJ7fVPt+Wtl/HZ3iQoCoC8SZVPildU6k54srEJR4H04t9gkpe4OQv01UNYEclw9FB47WmJL28aJsccJUfd1caHsrjO1yyTUh/0S32dPkZ1fYjJOkqdzx++j2wc5tTF58mQUFRXxP+fPO6eG4MuHu2HNywNsWtZQz2INW39hzt1dmtm0rvrQ1E9YC9KxmXQ6cVSvaJx8f5jJ8wFe1svCyqq08GYKowHgswe74rZ4YT3WY0ktrK6rrj4e2cXmZZ8Z0MrkuV4tQzBucGvBc2EBXmjqrzZZ9renewsem/SiGZ7A/27I9iREBZqsh+05Zw0bUBpEBJg+5yo+agXfs4sd8Xhg21BXbhbP3PFvsMuGplBbBPuo4G9m9vKbosyNtV6CrrRBNDRBQ+LoZhl30TzI22TgQlvMH90Dof4aJ2xR3bj9/1LTpk2hUCiQny/s05+fn4+ICOmuhBqNBgEBAYIfZ5DJZGgb7mfyfJfoIJvef/qDO/DTk70EzzXxMz1Izsy4Q3AH+MWobpj7cHfBMsemD7XpMw2kgjO2kt8gZ2Yqtr0xGEM6hGP2A6YX+MhAb/73359Nwhejukl+XoCXCkqF8HAL9ddg7SsD+cdfP9odPz7ZC2/d2UGwXN/WTZHQXHjxbtHUF4/3aSF4rlNz0wu8lPAAjSBQbB7kjaS4JoJlhnQwZg7ZYmhbAsx24f7Y/eZtAPT1NCyZTIZwUdAwsG2oSRFnZKAXbmkRLHhu5+TbTP4Osc6iIOf+HlH82COAaVB5QhR4Pi8KwADgnm7Gv/mJvi3w61OJJsvUF7lMxhdU6ziO7xVWm5OyI4gvdIktm5hZ0vGfKxXQAqaZHGLdI6Kxhh5PirX6nnhRRtYSqe+qFHuPY0dlBvu1bsr/nldcjmqJlMxzA4U3bYufEd6ENQ/yhjty+yBHrVajR48eWL/emMbX6XRYv349kpKSXLhleuLxQnJmpuLvcX2x7Y3BJstmvzcUvz3dGwFeSnz5cDfI5TLJC/OU1PYmn7Hk2STsmHQrnw1J7RyJDa8OxP8GtcLpD+4QdAGXkpoQiXbh+i/lS7e1QbsIf5x8f5igSSXEVy2ZcYoK9sE3j/XEfd2j+L8P0J8I2GkgYpv4IJYZqOzTB41B0b3dTQOoPW8mIyLQC1tfH4wlz/TG0E6RGNg2FGP7tcQ/4/vyy308sgs+EQVYMSE+6N+mqeC57jFBGNrRGPhqRBcgwwXpoVti8NHIzvzza14ZgGcGxAmWZdu8f306ETkzU5EzM1VwUevTqgmW1+yLP5/vg+n3dMRDt0RjzSsD+OanOzsbA4S7agKkMNHdTnSw8OQw75Hu2P7GrVAq5HyNzaO9Y0wuqFJj1QzpECGYcXpsv5ZoH2kM8kN81fzJedWL/aFSyBHko19eIZehs+h47Cbap7fFh6Nva+F+r08ygO/mqtUZBwO0Z0h+RxJflLzV9XNKLSqrcsjfbC4bMbidMTM2tl9Lu9YpPn85gziza0mwjwqnPrgDfmYyX4D+HGcQ28QHbwwzbWoWk+pMAMDkXAIA0cE+EksCUcHeiA4xfv+fGRCHqXd1kFzWwHDz89qQtljyrOk18PA7KSZBmyVtwvyQzXRpF9cwju4di0PThphMydE7rglyZqbiq0e6Y0pqe5tvMutbg+hCPmHCBDz++OPo2bMnevXqhc8++ww3btzAE0884epNA6APbG5WVsNHbdyd7JfGEPFqlAoktWqCg9NS+NfETRAAcHfXZvyAXoaLHAA0E0XKcaF+eH2o+S9jUz8NX5tyV5dmGCoaREupkKNnC+OBa9iWb0b3wDM/Z+C/l/pLrrdLdBAfDLE9zAwZip+e7IXsvBIM79oc3aKDkVdcjvgI/YX2wNtD8NbfWXg5uQ3/vugQH8EEggDQOSpIEHCJXxfvt/gIf8hkMrxye1t+RNaKah0WPZ2Ih7/dhWcHxGF0UiyOXirBbfFhkMtlgvUPjg/DswPiMH/LaaR2jhQ0BbQOFWbrdky6FSsPXsKYvi2gUsj59fSIFWZeAAgCjNCaLF1TUbYuPMALpz+4A+uPFaBXyxBBkHLi/TsEy94WH4b1xwr4E+l/L/XHsM+N3U99NcJup23C/BAZ6MX3NOsd1wQzR3QWrHPn5Nuw5nAe+rcJRRDz2f+O74e4UF/4apToERsMrY5Dv5rA8sGe0SbDHNQHmcw4HxA7GKBS4Zog56Yoa6JR2n7xlTJpWDxm/me9J1R5lc5kzqouUYE2dwE3+PSBrhi3aB96xgZjb814MiqFDH5exuPA10JwICXIx/Sc1jrMz+LUF4bj2pzuMUF4684O/FhiIb5qyZnuDVITIrGyZgyqdRMGQiGX4d8X+uHtv7Ow/eQVjBvcWjCkQyTTTPvirW3go1bi5PvD8L9f96Fni2Dc0iIEWh0nGHhV7Lene6NHbDAqqrWCXpOAPpgx7N+PR3bBR2uO4XJJBf4Z3w9fbz7FL3+zUosn+upvTB76Zqfk5wyOD8PS5/rwjz8ckYAVBy/hnbs7wkulgK9GiffvTUCYvxfmbjqJMH8NX5z+/eM9MfZH4Vxz4o4afVs34XtlAjVNo14qwfWtU3Pjee0ON69NbRBBzoMPPojLly/j7bffRl5eHrp27YrVq1ebFCO7EnsAGOTMTIVOx0Fu4x3XE31bAADC/L0Q6K0fr+D4e6Z1LLaIa+qL11La8fPM9DKT1mRrCAwDVA3pGGFzDdFvT/fGqG93IorJRgxoG4oBNTUSLZr6okVTY0FwoI/KbJOWNYPahWJTtnB8HvF2skW96yYMROswP2S/N5S/+ESZuaMCgMl3tMfkO/R3oVodh7/2X8T/3dHepJmtWZA3npa4W5MS6K3i/y9fuk0f2Invnpv6ayCXy3B7B+vH83eP98SxvBK0CdMHXu0jA5AzMxUcx0lmdpQKOd65uyNf+yBuCgT0A0He09WYaZP6v//z+T6CxxOGtMWBC4WCQebqg1wm4/9O9o5TIXePpLQ4e2iv5wa2MhvkDGwbKjk+lYG1YOThxBiTMXJSO0eif9shUCvkiH9LPyZVlZaDr9oYrJn7m9iA7J27OyK3qAytmvoJlmePpfdWHMF3286gqZ/apKg9tXOkSZAT4KXkp7v436DWgnoPSwEOAMx9pDumlVTgZmU1XwLQsqkvfh6rb2qtrNYJgpxAH2NQZxj4UqmQ45vHegrW2zzIm//sEF81hnWK4HujqpVyqJVyySancCaICvZRYcUL/VGp1SHEVy3IwhuO6d5xTfDv+H74bc85jOwRhfGL9vOfe1M0YvWDt8TgwVtMMzcvJbfBc4PioNVx+HLDSdzZuZmg80FiyxDsOnMNj/aOxcnLpfyx8dXDPdDlXePgqP41AS97rWjRxPz0H+6mQQQ5ADB+/HiMHz/e1ZthN1sCHKmLyoGpQ+z+LH8vJT9myIu3tUFy+3ColXLEhvhIZowA4Z3n9OGd7P7MpFZNbA6I6mrhE72sLwTT/Vmbu2uFXIbsWgaYYtb+Ly2l0cVkMpkgO8Q+b2DI7nxbc4KODvHBgalDIJPZf1duTniAF1a/PMDusaJm3d8ZrzNDxNtLX5Oj/52tG3BVTY6YOCB2pLu7NKtTkPPG0HhBkPNGTRY4gMnaAPoatGDmfHGzshqThsVjzvoTmHpXR36YhMhAL6x+uT92nb6G0b1j+XNdRbUWvVqGoKcoszlhSFtEBnljSIdwzPzvGJ9p0b/HtEs0O59XebXWJANqjmH+LX1QJP0etVIuCFg4jkOvFiE4nFuE5Pbmbza2T7oVD3+7EztOXcXy//VFiJ+aD3KaBekDGblchogAL+TVjCr+9aPdBbO2q5VyQcC2gJlxfhDTTJgQFYiEKH3ngl4tQ/ghHXq0MM0Ym2M49xky/myW99vHe+LA+UL0bdUUf2Rc4I8Nfy8lEpoH8vO2Gd7jq1HCW6VAWZXWphsyd9FgghxiXY/YYD7TkdIxAmqlHNnTh1qdZ+jE+8NQXqXlI3bSsBmyOyy2Cay+zLwvQTBbMwD0rmNhrkxmvHFgu/vbWp9yd5dm+OeA6VhLjuLoYKtZoBc/jYO3Whisi0c6V4i+58v+1wf3MVPFBHqr8OfzSfzUIn4WejaGMgHFvwcuYcvrg/F0/zjB3G43K7WIjwjgm6INNEoFfpeoFfFRK/n6nsl3xGNV1iVwHPDira0FhbkpHcOx5nA+Fj2diM/STuBYXjFuaREiyHj8b1ArfMVMcrn4md74evMptAnzw9h+tmVZt0+6lQ/SO0QG4penElFWpbX6XVkk6vX4z/i+KC6rFnTC2DRxEP7LuoSrpZVI6RiB33afx6nLN/j9w2rRxJcPKKRaBAD9/jpz5QZ6xgZjUB16EgZ6q7DptUHwUikQ4KVC/zb6dbHHllwuw/djeqLX+/o6WDZrvf7VgTh4oRApHaU7/bgjCnI8yAf3JqDPzA1oHebHH7S2TKSoUsihcuIdKDG1/63b8e6KI/zdtCd6qFeMSZCjUtYtCJAxmRx2rHaVjTU5zi5QZpvQ2DoUQ9OAwZ/P98GIeTtM3i82qlcMP2M4OxRAXKgvPn2wK2553zj3m3jwtlZNTXt+do02ZgF2nLyC0b1NexFdvVEpCCANgZtCLhM0Pddl4LeoYB9kTx+GQxeL0CUqEAq5DK8kt0XnqEAMaheK6zerEOKrRrfoYJRXafnM0tLnkpB5rhBP9W8pCHJ6xzVB7zj7A+iMKckoLq9GRE1zUm26hXeuGZuK5aVS4N5uUfzjeY9256d+EWe4fn82iZ++xtxxHObvxXdyqCv2/9BgSMdw9I4LQfeYYP7zpDL0zYK8TWpD3R0FOR6kWZA3Dr+TYlfPA+Iawb5qfPpgV1dvRr3zqmNhrlwmHbgrbazJcfTk6eIpQxbtNjYHdY4K5IOcd+7pKJifiG2iZKcyAYy1Z2wPJwBQK4z7Lrl9uMmYJNmiSR8VEhdMNsizNB1G77gQLNyRAwC4vaN004S8jjtTrZQLivVfYjojGJrXvdUKQZbhlhYhuKWFsL6wiZmmeFs08dNIDtvhaD5qpdlmfW+1AiO6R+H8tZuSAVN90CgVWPyM63srOwPdvnsYX43S5kJnQupqwu1tAZjvNvydqHAzyKduzWbsODksWwt+63phFntcNADl8K7GIQN2nb6Gf8b3xfeP9zRp0mEL9R+6RTgi++cPdsOHIxLw2UPdENPEWCjPZhmkMlJ9WgkzGWzT2aB2pk0cE1PaCR4bejwueioRYcxYTq8ktxUsN/2ejhjYNhQP9bK9m7IzbJ90K6YP74Rd/3eb9YXd3CcPdMHvzyV57ACDrkSZHEJIrb14Wxs8dEs0wgK8+GEPWMlMgeIryW0hk8ngp1GitJbzKcllprUnTf00Ngf2tY3/m/iqcVVinqW7ukTiB6ZwtHWYsYnoYmGZ4M6cnVzWV6PE7jdvg1ohN5nAMdBHxfeWSekYgeT2YRjQNlRwAVRJ/CELnrgF7aas5h+zgRD7GX8+3wcXC8tMBi19ObktXq4JaCqqtWgT5odWoX4mY3CNTmqB0fUwurg1zYO8JZvbCGFR2EgIqRPDXb94FG6xiEB9s0Bd7rxlMplJk1N4gMbmZii5xPttep+Z6Ejc081SV3ZxD8cwfy8E+ei7IU+/p6NgAEwDL5UC3z1+Cx5LaiGs16j5I14fqs/GTB4WD7Work4hk/E9hZYxQwD0iA22OnK3RqnA2lcG4OvRPSwuR4i7o0wOIcQhUjtHYtwi86/LoL8w16Ubu0xm2uSkVMj5dVt/f+1SOeYyQOJtYR/5iHpDJdU0J0WK5geTyWQ2ZUbY6TnWZOVhwu1t8b9BrfFo71iTbuCAPjD77vGeJs/bqrb7ihB3QkEOIcSpOjUPQNbFYgwWTahaG3KZDOJkiUous7kZqrbNVeZqecRPZ5wzToa59hXh/HBh/l7Y99btJsGPrbowc1XdrDI290kFOLUZ84oQT0RBDiHEqZb/ry/KHDQOU15RuUnAoZDb3gRVmy7k3WKCkF8zVg1Lo5SbbIuO6WklNd6KuUE5bcEONBgn0T2cpXbRNBeEuBuqySGEOIx48lFAf3F21ECTpRXVJg1TKoVpsGGOeDnxTO8A8NH9nbHptUH84/aRAfyAfADw89heeDgxBjsn32ayLc2ZXlP2jGRtL9M5ooXaRZiOik1IY0RBDiHEYRw9zod4njOpcXKUCtOKHHMjD8tkwrqZxyRqYUb2jBYMmCYO3LrFBOODexMQ7Ks2ySDd1j4cSrkM/Vo3dWpNC8dJhzmrX+6Pb0b3QFdRzylCGisKcgghDtOrpW3z6nxwbwKS21uv0ekdF4IP7k3gH0s1TSnlMpOAwlyzlFwmE2RBLGVEPh7ZBcntw/gZ3/l1M58l/tzmQd7YMfnWOhX82sJMjIP4iAAMaUBD7hPibFSTQwhxmCf6toSvRok+rZpaXO7hxBg8nBhjdYJPpVyOXi2NI9yWV+kksjZyk8DHfJBjPkAQu79HFO7vEWXyvLXBlcP8vSwvUAcjukfhz30XMP7W1k77DEI8CQU5hBCHUSnkeCTRcQO0KRUysPmW0opqiS7kMkEXcpnEgIEG4vcO6RCOtuF+OJ5favM2mZvVftaIzjavo7Y+HtkZU1LbC2YJJ4SYR81VhBC3pZTLTDIv4vhFX3hsfKzvZm6uJkf4vJdKgbWvDERnpnt2bUlNfOhoMpmMAhxC7EBBDiHEZeY90h3jBrcy+7pUs5O4zFghlwl6P8ll5sfDUZg547GTbNYWO6UDIcQ9UJBDCHGZYQmRmJgSL3iOnRBSPFUBIJXJkeE3ZvZvc5N4Gl6TUpcgZ++UZGx9fXCdxsAhhDgHBTmEELfCjjUjE/WG0j8nfKyUyzGmTwvB6+a6b5t7XmdrNbKEpn4aRIf4WF+QEFLvKMghhLiV8iqt4LFKlM2R6i7eNtzf+Fgmw5XSCsl1m2vGckRzFSHE/VCQQwhxK7d30M+c3bGZftTeFk2EWRJxoLLz9FUcvFDIvG5+ED5zr9UhkUMIcWPUhZwQ4lbCA7xwYOoQ+NZMZCnO3IgLj4/llQjqYcRxjFwGGBI15jI5wb5q4MoNq9vmrard5JqEENegTA4hxO0EeqsEE1KypJIx7HPi7uNsDy1zXcvfv7cT2ob74atHuku+/kBP/aCAS59LsrTZhBA3Q5kcQkiDIhXksM1NUrOUV2k5ydcM4iMCsPaVgWY/88MRnTH1ro7wdeKkm4QQx6NMDiGkQRE3V0UFe1sMcpTMPAxyGRBaM+Fmh0jbZ+qWyWQU4BDSAFGQQwhpUMTJmBdvayN4LG6REo+G/MvYRIzsEYX5o3s4aQsJIe6CghxCiMu9f28nAMBH91uf/0k8pk2Al0rwuKBE2H2cre2RyWRoF+GPj0Z2obFtCGkEKP9KCHG5RxJjcW+35vBRWz8lhfl7YUDbUGw5frnmGU6yTsdAUHhsYTlCiOehTA4hxC3YEuAYfPlwN/53a2PcsDOSWxpDhxDieSjIIYQ0OOJgpV2Ev5klbetCTgjxTBTkEEIaHDY7wwF4ZkCc+WWpuYqQRouCHEJIg8P0CgfHAV5K40jEXaODcH+PKP6xMMihKIeQxoSCHEJIgyMXZHI4k+ClCTPNA2VyCGm8KMghhLi9RxJjAABP9G0BQNRcxQEy0ZmMne+KXVY8DxYhxLNRF3JCiNt7M7U9EuOaYGDbUADCAmIOUlM5sL9TcxUhjRUFOYQQt+ejVuLuLs0kXyu6WSkxyrHxCaWCmqsIaayouYoQ0qDtybkuqtERNkuxrykoyiGkUXFakNOiRQvIZDLBz8yZMwXLHDx4EP3794eXlxeio6Mxa9Ysk/UsXboU8fHx8PLyQkJCAlatWuWsTSaENECdowJNRjxWmAlsqCaHkMbFqZmcd999F5cuXeJ/XnjhBf614uJiDBkyBLGxscjIyMBHH32EadOm4ZtvvuGX2bFjB0aNGoWxY8di//79GD58OIYPH46srCxnbjYhpAF4flArtA7zw4O3RJvU2rAJG+pdRUjj5dSaHH9/f0REREi+9uuvv6KyshI//PAD1Go1OnbsiMzMTMyePRvPPPMMAODzzz/H0KFDMXHiRADA9OnTkZaWhi+//BJff/21MzedEOLm3hgajzeGxgMAtDrh3A5sYTJN60BI4+XUTM7MmTPRpEkTdOvWDR999BGqq6v519LT0zFgwACo1cbxLFJSUpCdnY3r16/zyyQnJwvWmZKSgvT0dIufW1FRgeLiYsEPIcRziTM0bCxDhceENF5Oy+S8+OKL6N69O0JCQrBjxw5MnjwZly5dwuzZswEAeXl5aNmypeA94eHh/GvBwcHIy8vjn2OXycvLs/jZM2bMwDvvvOPAv4YQ4s7EtTZs9ubi9TKzyxFCPJtdmZxJkyaZFBOLf44dOwYAmDBhAgYNGoTOnTvjueeewyeffII5c+agoqLCKX8Ia/LkySgqKuJ/zp8/7/TPJIS4D7ZZ6vSVG5LPE0I8n12ZnFdffRVjxoyxuExcnPREeYmJiaiurkZOTg7atWuHiIgI5OfnC5YxPDbU8Zhbxlydj4FGo4FGo7G4DCHEc5mLZRQ0aAYhjYpdQU5oaChCQ0Nr9UGZmZmQy+UICwsDACQlJeHNN99EVVUVVCoVACAtLQ3t2rVDcHAwv8z69evx8ssv8+tJS0tDUlJSrbaBENI4sBmbUb2i8dtufTaXmqsIaVyccl+Tnp6Ozz77DAcOHMDp06fx66+/4pVXXsGjjz7KBzAPP/ww1Go1xo4di8OHD2PJkiX4/PPPMWHCBH49L730ElavXo1PPvkEx44dw7Rp07B3716MHz/eGZtNCPEQbLdxJTNleUWVzhWbQwhxEacEORqNBosXL8bAgQPRsWNHvP/++3jllVcEY+AEBgZi7dq1OHPmDHr06IFXX30Vb7/9Nt99HAD69OmDRYsW4ZtvvkGXLl3wxx9/YPny5ejUqZMzNpsQ4iHMjZNzLI96WhLSmDild1X37t2xc+dOq8t17twZW7dutbjMyJEjMXLkSEdtGiGkEWCbpVQK6TFzCCGej8rwCCEeR25mAEA5DZRDSKNCQQ4hxOOY60WlpCCHkEaFghxCiGfhOGEvKuZXGieHkMaFghxCiMdhg5nLJcYBSBWUySGkUaEghxDicdhYZsvxK/zvdyREumBrCCGuQkEOIcTjsBkbtndVRKCXKzaHEOIiFOQQQjwOW5PzVH/pqWYIIZ7PabOQE0KIK2g5TtBc1TsuBI8lxSIy0Nt1G0UIcQkKcgghHiXrYrGg8Fghl+Hde2iUdEIaI2quIoR4HHODARJCGhcKcgghHkcuk/6dENK4UJBDCPE4wuwNRTmENFYU5BBCPI6cObNRJoeQxouCHEKIx6GaHEIIQEEOIcQDUZBDCAEoyCGEeCA2sKEYh5DGi4IcQojHEfSuoqIcQhotCnIIIR4hsmZeqm4xQYLAhmIcQhovGvGYEOIRljyThF92ncXYfi1x9upN/nkZdSEnpNGiTA4hxCPENPHB/93RHuEBXjQYICEEAAU5hBAPJBMUHlOUQ0hjRUEOIcTjUCaHEAJQkEMI8XA0Tg4hjRcFOYQQj8Mxv1OQQ0jjRUEOIcTjaHXGMEetpNMcIY0VffsJIR6HDXJUCsrkENJYUZBDCPE4OibIUSroNEdIY0WDARJCPE732GC0aOKD8AAvV28KIcSFKMghhHgcL5UC618dRGMdE9LIUZBDCPFIChogh5BGjxqrCSGEEOKRKMghhBBCiEeiIIcQQgghHomCHEIIIYR4pEZReMxx+jEziouLXbwlhBBCCLGV4bptuI7bq1EEOSUlJQCA6OhoF28JIYQQQuxVUlKCwMBAu98n42obHjUgOp0Oubm58Pf3h8yBk/UVFxcjOjoa58+fR0BAgMPWS4RoP9cP2s/OR/u4ftB+rh/1sZ85jkNJSQmaNWsGudz+CptGkcmRy+WIiopy2voDAgLoi1QPaD/XD9rPzkf7uH7Qfq4fzt7PtcngGFDhMSGEEEI8EgU5hBBCCPFIFOTUgUajwdSpU6HRaFy9KR6N9nP9oP3sfLSP6wft5/rREPZzoyg8JoQQQkjjQ5kcQgghhHgkCnIIIYQQ4pEoyCGEEEKIR6IghxBCCCEeyeOCnLlz56JFixbw8vJCYmIidu/eDQC4du0aXnjhBbRr1w7e3t6IiYnBiy++iKKiIqvrXLp0KeLj4+Hl5YWEhASsWrWKf62qqgpvvPEGEhIS4Ovri2bNmuGxxx5Dbm6uxXVu27YNffv2RZMmTeDt7Y34+Hh8+umnNv89rlbf+xkApk2bhvj4ePj6+iI4OBjJycnYtWuX1fVu2rQJ3bt3h0ajQevWrbFw4UKb/x5Xayj7uSEfz67Yx6znnnsOMpkMn332mdX10rEs5Iz93JCPZcA1+3nMmDGQyWSCn6FDh1pdb70cz5wHWbx4MadWq7kffviBO3z4MPf0009zQUFBXH5+Pnfo0CHuvvvu4/755x/u5MmT3Pr167k2bdpwI0aMsLjO7du3cwqFgps1axZ35MgRbsqUKZxKpeIOHTrEcRzHFRYWcsnJydySJUu4Y8eOcenp6VyvXr24Hj16WFzvvn37uEWLFnFZWVncmTNnuJ9//pnz8fHh5s+fb9Pf40qu2M8cx3G//vorl5aWxp06dYrLysrixo4dywUEBHAFBQVm13v69GnOx8eHmzBhAnfkyBFuzpw5nEKh4FavXm3T3+NKDWk/N9Tj2VX72GDZsmVcly5duGbNmnGffvqpxfXSsSzkrP3cUI9la9vlzP38+OOPc0OHDuUuXbrE/1y7ds3ieuvrePaoIKdXr17cuHHj+MdarZZr1qwZN2PGDMnlf//9d06tVnNVVVVm1/nAAw9wqampgucSExO5Z5991ux7du/ezQHgzp49a9f233vvvdyjjz7KP7b376kv7rKfi4qKOADcunXrzC7z+uuvcx07dhQ89+CDD3IpKSm1/nvqS0Paz1IawvHsyn184cIFrnnz5lxWVhYXGxtr9eJLx7KQs/azlIZwLHOc6/bz448/zt1zzz12bWt9Hc8e01xVWVmJjIwMJCcn88/J5XIkJycjPT1d8j1FRUUICAiAUmmcwqtFixaYNm0a/zg9PV2wTgBISUkxu07DemUyGYKCgvjnBg0ahDFjxph9z/79+7Fjxw4MHDiw1n9PfXCX/VxZWYlvvvkGgYGB6NKlC/+8eD9bWy/tZ8fsZ7GGcDy7ch/rdDqMHj0aEydORMeOHSU/i45l1+xnsYZwLNd2uxx5zti0aRPCwsLQrl07PP/887h69argdVcdzx4T5Fy5cgVarRbh4eGC58PDw5GXlye5/PTp0/HMM88Inm/VqhWaNm3KP87Ly7N5nQBQXl6ON954A6NGjRJMWBYTE4PIyEiT5aOioqDRaNCzZ0+MGzcOTz31VK3+nvri6v28YsUK+Pn5wcvLC59++inS0tIE6xHvZ3PrLS4uRllZGe1nM+u0dz8bNKTj2ZX7+MMPP4RSqcSLL75odvvoWHbNfjZoSMcy4Nr9PHToUPz0009Yv349PvzwQ2zevBnDhg2DVqvll3HV8dwoZiEXKy4uRmpqKjp06CCIWAFg/fr1tV5vVVUVHnjgAXAch3nz5gle++mnnyTfs3XrVpSWlmLnzp2YNGkSWrdujVGjRtV6G9yJM/bz4MGDkZmZiStXruDbb7/FAw88gF27diEsLAyA+f3sydxpP3vq8ezIfZyRkYHPP/8c+/btg0wmM7scHcvTBK/V93721GMZcPw546GHHuJ/T0hIQOfOndGqVSts2rQJt912GwDXHc8ek8lp2rQpFAoF8vPzBc/n5+cjIiKCf1xSUoKhQ4fC398ff/31F1QqlcX1RkREWF0nYAxwzp49i7S0NJunnW/ZsiUSEhLw9NNP45VXXuEPOFv/nvrm6v3s6+uL1q1bo3fv3vj++++hVCrx/fff273egIAAeHt7036WWCdg/342aEjHs6v28datW1FQUICYmBgolUoolUqcPXsWr776Klq0aGH3eulYll5nbfezQUM6lgHXnzNYcXFxaNq0KU6ePGn3eh19PHtMkKNWq9GjRw9BFKrT6bB+/XokJSUB0EevQ4YMgVqtxj///AMvLy+r601KSjKJbNPS0vh1AsYA58SJE1i3bh2aNGlSq79Bp9OhoqLC5r/HFVy5n6Ww+6w266X97Jj9bO097rifXbWPR48ejYMHDyIzM5P/adasGSZOnIg1a9bUer3uuI9t3S532s9S3P1YtnW76uucceHCBVy9elWyGdDW9TpsP9tVDu3mFi9ezGk0Gm7hwoXckSNHuGeeeYYLCgri8vLyuKKiIi4xMZFLSEjgTp48KejqVl1dza/j1ltv5ebMmcM/3r59O6dUKrmPP/6YO3r0KDd16lRB97nKykru7rvv5qKiorjMzEzBeisqKvj1jB49mps0aRL/+Msvv+T++ecf7vjx49zx48e57777jvP39+fefPNNm/4eV3LFfi4tLeUmT57Mpaenczk5OdzevXu5J554gtNoNFxWVha/HvF+NnRTnDhxInf06FFu7ty5kt0UaT/XbT831OPZFftYilSvHzqWXbOfG+qxbG27nLWfS0pKuNdee41LT0/nzpw5w61bt47r3r0716ZNG668vJxfj6uOZ48KcjiO4+bMmcPFxMRwarWa69WrF7dz506O4zhu48aNHADJnzNnzvDvj42N5aZOnSpY5++//861bduWU6vVXMeOHbmVK1fyr505c8bsejdu3MgvN3DgQO7xxx/nH3/xxRdcx44dOR8fHy4gIIDr1q0b99VXX3Fardamv8fV6ns/l5WVcffeey/XrFkzTq1Wc5GRkdzdd9/N7d69W7AO8X42bFPXrl05tVrNxcXFcQsWLLD573G1hrKfG/LxXN/7WIrUxZeOZdfs54Z8LHNc/e/nmzdvckOGDOFCQ0M5lUrFxcbGck8//bRJIOKq41nGcRxne96HEEIIIaRh8JiaHEIIIYQQFgU5hBBCCPFIFOQQQgghxCNRkEMIIYQQj0RBDiGEEEI8EgU5hBBCCPFIFOQQQgghxCNRkEMIIYQQj0RBDiGEEEI8EgU5hBBCCPFISldvQH3Q6XTIzc2Fv78/ZDKZqzeHEEIIITbgOA4lJSVo1qwZ5HL78zKNIsjJzc1FdHS0qzeDEEIIIbVw/vx5REVF2f2+RhHk+Pv7A9DvpICAABdvDSGEEEJsUVxcjOjoaP46bq9GEeQYmqgCAgIoyCGEEEIamNqWmlDhMSGEEEI8EgU5hBBCCPFIFOQQQgghxCNRkENIA6TVcSiv0rp6MwghxK1RkENIAzTs8y3oNHUNblZWu3pTCCHEbVGQQ0gDdDy/FNU6DpnnCl29KYQQ4rYoyCGkAdNxrt4CQghxXxTkENKAcaAohxBCzKEghxBCCCEeiYIcQgghhHgkCnIIacBkqN1Q54QQ0hhQkEMIIYQQj0RBDiENGBUeE0KIeRTkEEIIIcQjUZBDSANGNTmEEGIeBTmEEEII8UgU5BBCCCHEI9VbkDNz5kzIZDK8/PLL/HPl5eUYN24cmjRpAj8/P4wYMQL5+fmC9507dw6pqanw8fFBWFgYJk6ciOpqmpSQEGfjOCpqJoQ0bPUS5OzZswfz589H586dBc+/8sor+Pfff7F06VJs3rwZubm5uO+++/jXtVotUlNTUVlZiR07duDHH3/EwoUL8fbbb9fHZhPSaP2deRE931uHjLPXXL0phBBSa04PckpLS/HII4/g22+/RXBwMP98UVERvv/+e8yePRu33norevTogQULFmDHjh3YuXMnAGDt2rU4cuQIfvnlF3Tt2hXDhg3D9OnTMXfuXFRWVjp70wlxS2yG5fU/DjjlM15anImrNyrx9E8ZTlk/IYTUB6cHOePGjUNqaiqSk5MFz2dkZKCqqkrwfHx8PGJiYpCeng4ASE9PR0JCAsLDw/llUlJSUFxcjMOHD5v9zIqKChQXFwt+CPEUbCtSblG5Uz9LV8smq4pqrYO3hBBC7OfUIGfx4sXYt28fZsyYYfJaXl4e1Go1goKCBM+Hh4cjLy+PX4YNcAyvG14zZ8aMGQgMDOR/oqOj6/iXEOI+6rNSpjYxzrYTV9BuymrM3XjS8RtECCF2cFqQc/78ebz00kv49ddf4eXl5ayPkTR58mQUFRXxP+fPn6/XzyfEmdy9IHjSsoMAgI/WZLt4SwghjZ3TgpyMjAwUFBSge/fuUCqVUCqV2Lx5M7744gsolUqEh4ejsrIShYWFgvfl5+cjIiICABAREWHS28rw2LCMFI1Gg4CAAMEPIZ6ifjM57h1QEUKIJU4Lcm677TYcOnQImZmZ/E/Pnj3xyCOP8L+rVCqsX7+ef092djbOnTuHpKQkAEBSUhIOHTqEgoICfpm0tDQEBASgQ4cOztp0QtwaxR2EEGIbpbNW7O/vj06dOgme8/X1RZMmTfjnx44diwkTJiAkJAQBAQF44YUXkJSUhN69ewMAhgwZgg4dOmD06NGYNWsW8vLyMGXKFIwbNw4ajcZZm04IqVHXeOrU5VK0CvVzyLZYotNxkMtpigtCiJBLRzz+9NNPceedd2LEiBEYMGAAIiIisGzZMv51hUKBFStWQKFQICkpCY8++igee+wxvPvuuy7cakJcqz5nHi8pr8by/Rdr/f5jl0ocuDXS3v47C31mbkDhTRpWghAi5LRMjpRNmzYJHnt5eWHu3LmYO3eu2ffExsZi1apVTt4yQhqO+m6uenlJJoZ3a27z8jKZ9O/O8lP6WQDAot3n8L9BrZ3/gYSQBoPmriKEeASakZ0QIkZBDiENTEMqPKawgxDiShTkENLAOKsm53BuEf7IuEDdxgkhHqNea3IIIe4r9YttAIAgb5WLt4QQQhyDMjmENDDOTrQcuWQ619vWE5fx8uL9dvdgqo/CYzGtjkN+sXPn9CKENAwU5BDSwDi7MUlqUs7R3+/G8sxcfLhaOFWDOzZtPfPTXiR+sB7bTlxx9aYQQlyMghxCGhhXBhYXC8v431//4wCSZ2/G+Ws38c6/h3HqcqnEO+o/lbP+mH6E9B+2n6n3zyaEuBeqySGkgXGX3Mnvey8AAPrP2ggAWLLnPI68O1Ry2ayLRQjz1yAswHmT9YoLst0xy0QIqV8U5BDSwIiv3RzHQVZPxS+WAoeblVrJ508WlODOOfqi5pyZqU7ZLsB0nBwKcQgh1FxFSENjEuTU30dfKbV/6oT95wodvyESxHEeJXIIIRTkENLA1ee1/KhEzysxYXDBQaVwzWnGFT27CCHuhZqrCGlgxLUnOo6Doh4LfHMLy/DdVvNFvReuG4uTdRygVNTPthXerMKlojLrCxJCGg0KcghpYExrcuq2vsKbldh37rrNyz+5cA+O5UnPLl50s0rwWMdxUNRTSuXrzafw9eZT9fJZhJCGgYIcQhoYcUxT12keRszbgVOXb9i8vLkABwBuVFYLHus4wJmJHEuF0NRaRQihmhxCGphqnU7wuK6ZHHsCHHvpe345bfXYlH3ZeSt3kp93nsWc9SdcvRmENAoU5BDSgOQVlaPX++sFz7lTLyJxQDNleRYOXChy2uedKDCfVXJXby3Pwidpx5FzxXnBJSFEj4IcQhqQBTtMC34dPSu5eLwZe4gDrpLyaszbVLc6mfIqLTLOXoNWZ/p3/peVV6d1u1JxeZX1hQghdUJBDiENSEWVzuS5XaevuWBL6s8zP2dgxLx0fLPltMlr9TUGj6OwNUTVEkEbIcSxKMghpAFZdeiSyXNPLNyDaq1p8OMpthzX1938uCPHtRviAGymi6adIMT5nBrkzJgxA7fccgv8/f0RFhaG4cOHIztbOItxeXk5xo0bhyZNmsDPzw8jRoxAfn6+YJlz584hNTUVPj4+CAsLw8SJE1FdLezFQUhjMLRThOTzbFbAlRdPW4uMl++/iNVZpgGbs7y34gjmbjxZb59nDvs/QzEOIc7n1CBn8+bNGDduHHbu3Im0tDRUVVVhyJAhuHHDWHD3yiuv4N9//8XSpUuxefNm5Obm4r777uNf12q1SE1NRWVlJXbs2IEff/wRCxcuxNtvv+3MTSfELUUFe0s+b7hg3jVnG4Z8ukWyfsUdaHUclu+/iJeXZOK5X/Zhxn9HsfKgMNgpr9KiRKJeRar2KLFliNXPPHW5FN9tO4OP1mSbXaayWoe7v9yGKcsPAdCPHTR/8ynkFZUDAC5cv4lvt5xGaUXdbq7YANTa/1BFtRbpp66istpzs3SEOJtTx8lZvXq14PHChQsRFhaGjIwMDBgwAEVFRfj++++xaNEi3HrrrQCABQsWoH379ti5cyd69+6NtWvX4siRI1i3bh3Cw8PRtWtXTJ8+HW+88QamTZsGtVrtzD+BkAZBx3E4WVCCQxf1PZkOXihEt5hgF2+VqWX7LmDiHwf5x/M36+tsUjsbJ+7sMT0NNyq1OPJuCnzUlk9R8RH+2HXGck1SmZmJQ1mbsgtw8EIRDl4ownvDE/Da0oNYdzQfS/acx4bXBuGuOdtw/WYVzl67gfeGJ1hdnzn2ZHLe/fcIft11DqN6RWPGfZ1r/ZmENGb1WpNTVKQ/AYeE6O++MjIyUFVVheTkZH6Z+Ph4xMTEID09HQCQnp6OhIQEhIeH88ukpKSguLgYhw8flvyciooKFBcXC34I8QTmej5duF4mGDPGUiKH4zjMWHUUP6fnOHjrrF+41x7Jt7wAgBs1QcmpAsd3sTbXlCfeX5uPFwAATtd0875eM5LzthNX6vT5OjvaqH7ddQ4A8Nvu83X6TEIas3oLcnQ6HV5++WX07dsXnTp1AgDk5eVBrVYjKChIsGx4eDjy8vL4ZdgAx/C64TUpM2bMQGBgIP8THR3t4L+GENcwV/NypbQCCjn7ovmL6aGLRZi/5TTe+vswzjh4rJZp/0jfeBjIzWy/eDoIwPRvvSmRkbElZGDXYy74q6/JPKkOh5D6VW9Bzrhx45CVlYXFixc7/bMmT56MoqIi/uf8eboTIp4txFctCHI4Tp+1+HhNNtYcFt4MlJYb60pe+G2fQ7fDWqbG3EV+2r+WgyNAP+bOrNXHUK3V4WRBKTiOsztoMJfJkYuiHHMZM5MpNeoQtVDvKkKcr17mrho/fjxWrFiBLVu2ICoqin8+IiIClZWVKCwsFGRz8vPzERERwS+ze/duwfoMva8My4hpNBpoNBoH/xWEuK+9OdeQX1zOP9bqOKw7WoAva3oU5cxMlXxf1kXTptztp+rWJGOJucv64VzTUZHXHy1A8yBhofVXm05h+f6LyC0qx0f3d7ap+YcNWMxlcsxlmMTYj5v2z2GsO5qPlS/2R6C3CquzLsFbrcTAtqE2vZ+DPtDR6jgoFTSaByHO4NRvFsdxGD9+PP766y9s2LABLVu2FLzeo0cPqFQqrF9vHKY+Ozsb586dQ1JSEgAgKSkJhw4dQkFBAb9MWloaAgIC0KFDB2duPiFuR2amXWXHqasoqzT2wtFynCDoEa7E8mfstlLIWxfWshfs65+uO47hX203WSa3psfTxD8OWqw9kvJTeg6GfrYFBcXl+DvzImb+d0x6fi0bgp6FO3Jw4XoZ/si4gILicjz3yz48/sNui38j20OM44CXFmei27tpuH6j0r4/hBBiE6dmcsaNG4dFixbh77//hr+/P19DExgYCG9vbwQGBmLs2LGYMGECQkJCEBAQgBdeeAFJSUno3bs3AGDIkCHo0KEDRo8ejVmzZiEvLw9TpkzBuHHjKFtDGh1z114dxwmyEVqd+Ykx6zJtQ12Zu/4bnn9lSabg+bNXb1pZn31RznsrjwLQB1CGgt4+rZqYLGduD0l1Y9fpOL4wGbA887pOkMnh8M+BXADAX/sv4sl+LaXfRAipNadmcubNm4eioiIMGjQIkZGR/M+SJUv4ZT799FPceeedGDFiBAYMGICIiAgsW7aMf12hUGDFihVQKBRISkrCo48+isceewzvvvuuMzedELejs5C24DggwFvFP7Y0ZYCtTTPOUGlmZGbD1i7PzLVvfbUcQ6aCed+V0gqTwM+eQmQOwgBTPEu8YFkmKGPHMrp6o8L2D3QQfU2T+9UFlVZUY+fpqxaPd0Js5dRMji1fIC8vL8ydOxdz5841u0xsbCxWrVrlyE0jpEFZuP0MPll7HLd3CJd8XaNSoFWoH/9Yq+XMZmzMNXnVh61mumDbMpaNFEs1OYa/U+rPVTJRiUwGQeqG44T77vw1y9mkAC+VYJ9WazlozJxZ2a1lgxy2GLw+cByHB+fvRKVWh2XP94HclZGvyOM/7EbG2euYelcHPNGXslukbqjajZAGYNq/R1BSUY1l+y9Kvn5X50j4qBX84/AAr3rrFu0IFwvLajVKc21v9tmeaDqdsHmK44SB0WM/CDs+SGGX11qqyWFeYgO0+k5aFJdVY3fONWSeL8Tl0vrPIlmScfY6AGDJHuoVS+qOghxC3FBBcTme+WmvzQWpKoVcUC+i4zizdSXuGvx0n55m93us9a4qLq/C8fwSk+fZLuOvLj0gyMQculgk2HfsWEJSHyfenxabWZiX2JY7S4GRPQ7nFtk0J1gV06SmcKMsDkvcrd9RNmUX4NTlUqesm7ifeulCTgixT68P9D0O1x5Jw+F3UqwurxONGaOV6jFUwz0vaUBRmemAgNZYC3Ju+2QzLpeYZiosXdjvmbsd/ubamyTIRA2DbJ0Qx3F45ucMaJRyzBnVTbC9WsGkqjZ/nEWpX2wDAPz1vz4Wp/Vga5JqO8/ZjYpqlFdp0cTPcR1ALhaW8b/LnXALfiyvGGMW7IGvWoHD7w51/AcQt0OZHELcnC13neev3cRTP+3lH4vrSljumsmpDQs1vgAgGeAApkGOeJeUmJmI0xCMiLM17MNeH6znm1wKSiqQdiQfKw5eQnFZtaAmhw14EluGYNfpq/hw9TGzxdQbswtMJjM1562/syy+XlFlrIGqMlMMbk336Wno8d46h3Z/f+m3/fzv7PH7444cPDA/XXLiVnt8s0U/V9qNWtaA1acqrQ7d3l2Llxbvt74wMYuCHELcnC032tP+PSJ4XMvrVoNjz1xQLC+VQvDY3sBP0LwkM82GvPmXfjbzlxdn8s/drKo227tKrZTjwW92Yt6mU2g75T/MTjsuWB/HcXhiwR6MW7QPeUXleH/lETy5cI/ZLIzUII+s8qq6Z3IM2SDDpLCOwDYtsv8nU/85jN1nruH7bWfqtH6NUmF9ITdx95fbcf1mFf7OzK11IEooyHG4ymodsi4WuWXXzMbkzJUbmL7iCK66WVGlvSIDvWr1Ph3HmfQYMvKcVE5tg5y6lqGIA4PVWdLz6KWfvip4j7neVWJfrD9h9vMul1Tg261nsOFYAdJPXcW6I/nIZZp5bFFezWZybNuH+89dxzWJrI29/wfnr93Er7vOoqJmG37bfQ5ra6YeYXt5Sf0X3TCTYbOVn6Z2Qc6lojJ8telkvQ7aePSSMVA1zNtWVqmtVbNuY0ZBjoNN+D0Td87Zhp93nrXrfbUd74NIe/2PA/h+2xk8zTThNBTsBU0uq93QfZXVOqw9bJxH6gW2GcBzYhy+Wai+sRd2GfSDC1qj1Qnrpti787etNC+x4wux4/Asz7yIp37ai9Qvttqw1UYVdmZytp64jHu/2sEXh7NBs7kC4cKblXwgw7r1k014868szN98Gmeu3MDkZYfwzM8ZAAAFuy6J9bL7r6isyupYOmWVWrzz72HsrAk2pbb1w9XHcOvHmywGD0kzNmDW6mxM+D3T4ufVVW5hGVZn5UGn45AUZxyk0nCs9Hp/Hbq8s7bOzXZ1wXEcVmflWR1awV1QkONgK2razL/edMrm92w5fhltp/yHGf8dddZmNTp7cvQXv33nCh263opqLc5edezM3WLlTL2ESiGrVS+Tb7acxrqjxiBnxcFLWHEwF8//klHrMWnc0XWJ2cttYTL4n42h5PWblRj00UZ0eHuN8b0S/z/H8kpMLvD6IMd4Uf56s/EccaXUcoZg+gpjcyR7POw4eaVmu+zbD8JMjvkbrGqtDmWVWny0JlvwPJv9kTo+r5ZWoOu7aRj00SbjZ9Yc14b3bj95RZAZ0k+vYTmTY3CyoBRd3lmLxxdY7t7/3dbTWLA9Bw99s9NkpYb/i3mbTuH0lRtYtOucxXUBwMbsy1aXsdfZqzeQU9ODb9BHm/DcLxn492AuWocZx70y/B8ZasWy80x7DNaXNYfz8dwvGeg/a6PLtsEeFOQ4iY8dvTMM43DM33zaWZtj0fzNp3DnnK0orfkCFd2skkxLOwrHcQ4ZzdSeJkFHNR8+uXAPBn60CTvqMInl7LTjuH/eDkEww2J7vlRpzfeSsmTbSdPtG79oP/7LyhNcXD2Zpd1WJapYPnvNtsD1ZqUWOVammjCY9o+4TkrYXHXqsu3BsmEKCkA4oKClgR0tHfNs4fH1m+a/60M/34qOU1fjlhYh/HNVWp0gmyTV9LfztH7+s0s184z9nXkR8W+txp8ZF/hldKL51Sq1OrDzlEr9aYbnluzRByTmBpc0OCfKNrDnWHEzndZaFbsTXCmtwMCPNmHQx5tQVqnlM3YvLc4UjFpeWa0TZNzsOX2WlFeh2oE1PTuZJtiGgIIcJ/FWuVeBW0l5FeZuPMkHMqwZ/x1D1sVivLBoH66UVqDLu2vRfXqa2YswAOzJuYY7Pt+KPTn2T+b4+II9GPr5ljoV07229AD6fbgRxTakbW9WVqPl5FVoMWllnYOd7Sf1X/C6DFT2xfoT2Hv2OpabGdiPzQDcrKx2WPdig8zzhY5dYQPEXmwB4M2/LDcX1cZvu4WZAa2om39tsRe7ixZqcfbkXJf8vgPCQHrqP4fNruNkQSl0HHAk11gfUlGtQ1U101wnEY2I63ReqinAfnXpAf65o5dK8L9f9/GPq7XC/WMpg2nrRZ7tRSf+7oun3zA36rPUDdmNimqz+9YeH/53jP9dPLUHG3RVVusE50tbJ60vKC5HwrS1eGB+OgD9+aSuKhpYaQUFOU7iyiDn/LWbJjNJJ0xbi4/WZKPT1DWC59lAZmP2ZfR8bx3/2HAXJmXk1+k4cqkYI79Ot2vbKqq12HL8Mo7nlwpOnPb6I+MCLhaW2dSllq1NMdd19HJJhdVmHPYkI3WSvVFRjc/WHbd44WEDGLOZHKZe4kaF1uJcSLVRUs9TCLgja81D9rC1PqJay0lO8GkvWy9U4xftQ6epa7DluGkTC9tL6fTlGzh4oRBZol5SbJ1gfonxXFBepUWaoCnUdL4xW4qRxUFCtZYTZC8MIUchk2kyvGxr8S0bgIkzN+LHCjNBlTjrV63VoePUNeg0dU2dez3tPCMsTBd8DpvJ0eoEwYWtU7OsqSno3neuEL/vPY8Ob6/B73UcSdoVGa+6oCDHSewZRfSWFuYH7aqN/rM24oH56SZ3klLYE82zA+MErzkyxTn6+11oMWkl/mYmYLxwXRgM2HqxYAMFpQ37uVmQN//77jOmqdaC4nLc8v46tH97tclrh3OL0GLSSqw9nIdCpu6hZRMfk2XHL9qHz9adQPInm81uS3GZcX+rlNJfP7bItFKro6J0N/eOqPu+OVodh6sOCK5Gfy9dhyLOahbUjBH0/krTWr+DF4QBzd1fbsedc7bh+21n8NrSA9DqOJQxQTjb9bq8SotrTNbh113nTHp3/W3nRKuA/lhnm8kN13E2oKnW6nDvV9vxhygTZw57ehDXsVRrdTh4oZBZVvpcUi0KhgqZ7bHU1AcA/xzIxey12WYzyOevGfdboaiuiv3eizM5JwtsG7GZDYxe/+Og/t8/D6KiWmtXL+Dfdp/DeyuOgOM4i5P/uiMKchyIDQrsqaNgjxlHjoewdK/liJ3jOMHEgOIvuSPSsYA+KDG0nbMnXPZvjZu8EgnT1qLFpJVW18eeDJQK4zbnFpZh7eE8ky/uDebOt7ImzX6zshpDP9uCt5Zn8aMLSzGMIPvMzxnC7qsS/8GGosQyC8187AlbfFIzYDM5AHDTwvqIee52Kq7U6nDnnG1OW/+zP2VIPm9tMEn2UJ6+4gj+yLiADccK8N8hY5aU7c7c78ONCPMXDm1guNgX3awCx3HYcKzA3s03maHecDHdddqYlf4x/Sz2S3QmuH6jUjITy95s3v/1DpP1s+dejUr6cigOcljHLlkuAH7xt/34YsNJk8y6FHHAxJ5vxEGOIWCxxlzT0vQVR3DnnG02B4uTlx3Cd9vOIL0Bzg5PQU4dcByHjccK+FFVy228496YXSDoocM2bd20s+dLbmEZJi87yH/B2QOwJ1MsKNZi0kq0nLwKK5kTmbj5pLSiGp+mHefTm8XlVRj88Sb8c8C+u7Stx43FgexFvrYTFLJ3ezcqjNt8++zNeObnDJPtu8ksU1Smf+/GY5dxLK/EpKu/pSCTvVMuq2XbNrsO8/USwv+H9UzTAGm4ip08vkm6mYLQah2Hu2qyNAZ9Whm7J9/dpZnJe57+aS/WHjF/3P2wXTgonwwyLN9/EV3eXYuWk1cJXrM1WzDqm514sGc0/9iQyfC10oljwfYz6DY9De3fXo0zV27gSmkFMs5eQ4tJK7F0r/EiLr7gV+uE87tVazks2nUOQz/bgktFxgwL21zVt3UTQdBjy+StAHD22k38eyAXBSXCEoCWTX3538VNcOz5oVIrrIMC9L3zDOf9nCs3UFBTxJ2dV4J3/j2M6zcqzTaJ/7JTn+W3d3DFK6WVfIYQ0HfRn7HqqNn6QndAc1fVwbJ9F/lCupyZqYI7CXPf612nr+KJBXv49wDCC2t5lRaB3irBe45eKsbSvRcw+Y54qEQVZ31mbgCg732RMzNVkLWwpesx2zW0XJRB2HL8Mr7dqv8SjOwZhed+zsCZKzfw4m/7oVbI+TsvrY6z2DynZpplesYGY2/N2CYTfj+AA+cLMe3ujoLlC0rK+TvFrItFiA7xEewTdlRU9kRgqLcRdwVl94mhq625CREvFZYjKtgb+85dR/vIAMFrbOZFXNtja5MSe4yYqwESn4wNJyRin9pkE5xpnQuD1UMXi3DoYhHG9msJAOjYLAA7TumDInM1Wpb2n7i5SyGX4eUlmZLL2noDc7mkAiql8Txi+E75WBnAj20uHPzxJsFrlrKq1VphdqRSq8O7NV31P1qdjdkPdq1ZTjhKtbkboa0nLuNwbjGeHRAHmUwmqLF5a3kW/702nPcB4eSv4oEG2f+X05dvmJzPZ/53DJdLKvC/Qa0wqObvzpmZipcW78exvBLkXLlhtcu7vYXw1Vodf9wA+mNkfs1UGcO7NbdvZfWEMjl1IO6VwEbNFdVaTFx6AC8t3o+F289g3K/7UFpRzX+JWGxtys1KLQqKy7E66xL/JR/2+Vb8sP0MH9BYwl70pQrE/L2UZi/IFaITAltAe7KgVJD2ZgMXQ6+A9Ufz0XLySpNxZNjMS5XojPdj+lkUi06ys1brA68VB3Nx55xt6PLOWsHrp5j2aKkRUMVf3GKJNnT2ObYm6vz1m5j4x0Hc/3U6Hv52p2A9bPBUVqlFfnE5rtSMqHy5VNwzQvrscdOmIIeapzwR2w3c1di6CkcEg5YKqm2dNqKsSivIuhqCAq2NIzLbq0rLCZrI2GbizPOF2FvTc7Ra1MtJ3Kw2O+04/jt0CS8vzsTM/47hxcWZeO7nDBzONQaC7I2LYXTn06JmxEJRJocNcqb+cxifrxOOgg3ou3OfENXnHKupPdprw0CZthTCs/9/X4nGf3tz+SHjutx0lH/K5NSBOHnB3jUcvVTCD0RnKMJLatUEh5keRYbBr9hg4nJJBQbPN/ZYYqN+c5MNstgvhuGCyh58vmqloHdGkI+Kz1CIm09WHTIOVf/NltPo06op/qpJS7LL/rgjBxNT4jH2R/3owgM/2iTY7itMAFAoUajH3s0A+p5TH4/sgvGLjKP0sgOFhQUY6wFKyqtRUa0V1BZFiKZCYDMwhTf0v09ZbuwyfOaKcSwNmQz4c58+xX1AdLe6+rBxf/y1/yK/L5aP62vyBS8pr8KV0grENvEVZN/YY0Rca5Nx9hqW789Fj1jHFqITYlBQUo55m05hwfYch67XUo89e6Z9YM+Phpsxc1nXuqrW6bCZyXSwzUWnr9zA/V+nY+WL/QQ3TZWi7A9gOgXHvzXN5ez5grX6cB5u7xBu0sNPXKMn7oixW2K4Do4TlijsO2cMbGwZ9NOW+NNSwTO7zZVanVvODUZBTh2wGQidjhM2RUikSUvKq3FXl2b8l+BmpdakvXnVIWGXaDYgiQv1hTk+agX/GQaGQIS9i5DLhMuwB+kVC/M8LbVQoBYbYrpd/xzIRadmAYgL9cO3W40DcEkNMig17843W4R3DCUV1Qjw0jdZsSej0opqtJsi7BUlbvdmC/qkekOwf7e5E4OfRmm2hmb43O34+tEegufGLNjDj0fz7MA4vJ4SD4VcJqjlMWSDvJQKBPqoMGKePri1d0oQQmzV633zRfZ18ZeFmgx7ghxx8xFQ+wlEranWcnxTCyB9bkj9Yhs6NjM2W+cVlZvUxthr2b6L2HL8Mt65u5PgeXGgWGJDx48jl4oxZ8NJ/vF9XxmLq23pBcX+36QdyceqQ5fw/r2d4KPWX5cKSspxqdD8UCKsguIKRIeY9jp1NQpyHOS1pQdw1coowQq5McoH9BdrX40SAV5KPmASz47M9iSIDfFBZbVO0FRkYGg3Zps6DAHM3hxjdC+TyZB+SrpAUdzk4qtWmB1XhnXtZqVJD44Xa+ZKypmZivxi43ql7vikMlQfrDpmskx5pRZLMy4I6ojWStwtnRaNJMuevMQpYTFzXXyt9TQTF5WyA+7N33waBcUV+PTBroLmqvziciTW9Oya/UAXi+snxJ1ZmhLBnhiFvegams+d1WVZfL42F6ix2aUrpZWCMYJq60ppJcYt2id4rqxKeI6xNTY0V3BuiyCm1tEwz1+QjwrPD2yFsAAvu4LiK6UU5Hi0ZTZUl4sLe6/frESzIG/BoFR5RcKsxiPf7eJ/35h9GQnT1uDZga0w4fa2guUMKV02a1N4sxJbjl8W9ACo1OoQ5CMsbDYQBxu2BDiAPjtjLl1dUGz9LuBYnvVBATdlXxbM32Mg1aOtQPR3sNmjm5XVFgcSe/1P27pm2vu+v/ZfxIz7EgTt5IeYwdcm/H5A6m2ENHgH7BhhWydqGnr8h93YLDGYoSP8uCPHpuVaNPERTOUhbp5ylHVH679Qft+5QmTnlaBdhD//3ILtObVqzhSfd91Fgwly5s6di48++gh5eXno0qUL5syZg169erl6s+ySnS8cU+GfA7l49fcDgqYtNushpaJahy/Wn4BOxyGPCSAMKV22eO76zSqTLo7lVVqzk1aKgzBbXS2tNDsolnjuGClbarqYy2TGu5eE5oGCIGCNmfZtW3o1Xb/B9Iqq0JoUMtcHX7UC8W+ZDjZIiKdjb9SsYc8XVVrOaQEOYHvBda+WITbPV1YXjhj009bsOyvlsy11/lxAP0aSO2oQvauWLFmCCRMmYOrUqdi3bx+6dOmClJQUFBS4tovo3Ie727W8eAqC+ZtP85XwBvvPW6+IB4AvN540Gcjp4IVCPPeLcUAwqSLfiiqdwydovHajAk/XFB2L3W/DtA8xNSnOgW1D+ecM41RoaprmbBlMyxw2AJPqjWWPYDNZMGvsPfEQQtyDuKm6ia/aRVtina1jtTmDePoLd9EggpzZs2fj6aefxhNPPIEOHTrg66+/ho+PD3744QeXbldq50hBUZojGLIpKR3D7Z59+u4vtwseF9cUOrPE3R8dYWP25Tq1mxtqgQK8VFDVjGBs6HngiMng2OYqa3VT1gR41y7IIYQ0TIdEc3pZm8rBlZxVpG2LKjedfsbtg5zKykpkZGQgOTmZf04ulyM5ORnp6dJZgoqKChQXFwt+nMXfyzktfr4aJcL8NXVeT74NNTGuZqjb8VYpbJ7Y1J4JUNkArK4nAfH/933d3XMALEKIY7DzSwH2FVI3JuvdbPBNA7cPcq5cuQKtVovw8HDB8+Hh4cjLk67TmDFjBgIDA/mf6OhoyeUcwd/LOXf2XioFgn3qnhYVj0FjC1smvQT03aodwdCzzFut4LsuGrQL95d6C0JclDIOEP1/d48JRniA5WDUXKG32Mz7Emq9XYQQ4kqG+QndjdsHObUxefJkFBUV8T/nzztvpNHaZHLY+UrM8VIq4K2u+8BKhh5Ti55K5OtbrGnqZ1sGSTz9xP09oswu2y0myOr6/DRKfrwfg3vNZEqa+gmDnOcHtcLycX2tfkZdiYMcL5VpYMZaPq4v5j3Sw+zrrLhQvzptGyHEft1iggTT0vSOMz/nX114mZkAtLbLWTKAqW+sTzdrOaefM7l9kNO0aVMoFArk5wvHJsjPz0dERITkezQaDQICAgQ/zlKbbEaoDUGERiU3O0fK23d2sPszA7xVJmPwmNPU37YsiTibEmthjATxbL3jBrfC4HbCL6JGKTcJ7Mz1OAgVzYIcE+JjNaMiNiW1veCxLQFrgLdwGS+V3GzTWai/Bl2jgxBqZX8mNA/EugkDHJYZI4TYrmdsMNTMqOTNg3xsuhG1VxNf285PlloHpMZIk+KIUof2kQF2F1mXWhj52lXcPshRq9Xo0aMH1q83Dkqk0+mwfv16JCUluXDL9MQT1dkiUeJOQZwV8VIqkHXRtJYoxFeNJ/u1xJkZd9j1mYHeKpvvEMxlcpLbh2N4V2Mhc6joi2SpUFo8AnRsiC/OirplalRyQSbHV63AoHbSdyTiz/aRaOoS79M7EoxBce+4EEQFewtebxYofCxFnMnxVikE28z+bgiarJ3cnh/UCq3D/OErmojQ1wGZPEIaq+T2YTYtV63joFQYT15qpdzihMPW/PdSf8nnL0qM7C4lQHSzNX5wa/53W88JtTl3dIkKFDxe+UI/mzpasOdVd5zWwe2DHACYMGECvv32W/z44484evQonn/+edy4cQNPPPGEqzcN7SOla0YskYqyxUOfmwtIDD2FZDKZXXf/0SE+FsfgYQMUc0HOw4nRgsDBW60QPO7YLFDqbQCAF29rI3isVMiQ1KqJ4Dm1Qg5vJlDx91Ihobn0OsVBjlRTl7hJK66pn2B59rOUcplgnez+YE944rssL5WwWZENggzLioMtMUMwJA7SghxQk0WIu4sI8LK+kJ3iI/zNZj3aRwoz+1odJ8jkaJTyOmVyxOu3lziwiGeuMeJs/F//64PTH5je8JobC82cduH++PmpRH6yYm+VAnK5DEM7GQMYc9eFqGBjBl98o+YOGkSQ8+CDD+Ljjz/G22+/ja5duyIzMxOrV682KUZ2hZ6xlttvpQ4MySnpRRX7XioFnh/UyuK6W4f5I+udFHz9qHG8Hn8vJX58shem39PR4nvF7c5eTARuLtXZPMgHGuZLplHKBePGiAMPlnidSoUc93UX1vColQr4qISZEJlMhsPvpJisTzxejZ9GCZVCzndBB4QZFG+VQrB9PmqlSQaGDUbahBkDIrYAXNyk5SXKPrFFxoY7MrnorvDLh7vxYwMZtt2wDSxXFVcTUp/evsv+5nex0b1jBY8VchlKK6THpvIX3Rje07W5SSbnveGdxG+rN+IbKV/m5kcc5HSLCTY5vwD2TYXx/r2dsOaVAQjwUmHpc31w+J0UHJ0+FAAEwd/ElLaS7w/2UWPdhIHY9NogKBXuF1K43xaZMX78eJw9exYVFRXYtWsXEhMTXb1JAKzPkLvk2d7IYi7SbcL84O+lwtS7OuC1IcaDpqSiWhB4eKnkCJcIGqSKh4d2isTJ94dh7sPdsX7CQAxsGyp54LMmprQTPGYzR+Z6dUUGeYnueBSCg1o8+zcb2DzQM1oQgCjlMpP2XrVSLpiw1HBHI57EFDANCPxqAgq2PqYJk8kpq9IimPk8X42wu7qvRikIYJoHGZuuQnyZwEV0l6VRCpvJLAVE+tdVuLNzM1Gzlspk2wHbe2URUlfsd7M+sM0pvVrWvdB3qihQkstk2GJmtGQ22zCkQzh6xAZDKTeexzLOXkd4gBcmD4sXvO8+0c2peFJea8Q1iOaw5w2lXCY455vrPLL2lQGCx0cvGUsd2Kz0j0/2EjTTZ0xJxiOJwgCRPd+y2bBIM835GqUcrcP80MIJdUyO0GCCHHd18brldtYALxX8NEo8OyAOAPDP+H4AgCf6tsT4W4VNOGFMMa2XSiG4KMfXzC2y6Onekp+jVMiR2jkSYTWp3+FdjV9IqcJYqWYXg0DRxbVHbDCGdYpAgJfK5At3ssA4MWdTPw1+elI/1cY3o3tgx6Rb8dmDXXFw2hColXLBZyjkMpMvRfqpq4LCXkuFwEq5XBBw+WtUJn+HOBPCNl/5qpUmTW9sALOLGWGZ/X/RKOWC7dJ3e5f+TKmmRFnNGYcd6dqfyfiYWxchdTXjvgT8+bx0HaM94045AnvxVMnlglq/2lAq5Pw5EhDe4IixF/FONc3hbJCXUTO/3DM152yDV0U3hrVp0npuoOXsPCBs8lYp5FAx+4o9v7E3pm1FQ22w50b2PjzYR4UtEwfj28d64si7KWhipRMMe76PC5X+ezUO6A3mTO69dQ3A2P4tLb5uuNBNvqM9cmamWuwWzl7gNEoFrjOj875zd0ccfXcoesQG27Rd7Lq2vjHY5HVx0S375RFncv58vg/m1dy1sAe0OJsC6Lsu5sxMxZCOEVAq5BjerTn/pWWXl5pywlejEMy6zn7Zd0y6VbBseICXYFsMmRwV8+VmRybtERssyCz5aJSCYE4lFxYbvsAEoF6iJjrhiUYhOGmygYmlruUsNmgylxUi7umFW1tbX0ikVwvndFE2J8hHf6N1b7fm6GGmed1a5tdW5lbzwb0JgmwLW6CqVMjQk9kntS36ZbOv79+bICjYZbE3H4bzhUqimUUmkyEpzlg3KK6ts3VIDgOFXIbnBsbh9aHtEBlovg6JLTxWKmSCgIUNbMSzKCx6OhGh/hrMH91DMJyGWhQkyeUy3N4h3Kbzk4xJA5mr//Ryw2JjFgU5dRTgpcLJ94dJFn8B9o15wAZAGpUcTZmLsp+X0q5xc2QyGc7MuAMn3x/G1wU1Y75Y7AGe2DJE8IW1VJnPnpx8NUq75nJiP7NdhGlx3v09opDJzFjMXvybBXkjZ2YqPnuwK8YNboW+rZsIvvxSBW+nLxsHQgwP0Ai6nStkMkH9T3Z+iSDwGtuvJZ7o2wLfPtZTEExpVArhSUcprMlpxpxopYJAqeeEzWbS9T0AMKK7+XGIiGuoa1GDYK3WztHWTxiIjLeSLQ4hUeiAyRUHtQvFb0ymuW24vq7NT6PEQ7dEo0t0EP8ae+FVKmRo0cSYJVDIZfhiVDe7P5+dsqZ5kLfgnCbotclcrA3febaGhY2xqplIQnyRl8pg9Ky5CRX3kAL05+QgHzX+N6g1HrolxuzfwZ739HWG7LnZ+Jq4VKJPq6bY82YyUjpGoAMz3RAb/NkbmK08mGv8bDNBDmVyGgGlQg65XCYZHMisTECVPvlW9G3dBOtfHSj4InqrFEjpaKxst3WMBfFnszUzOybfhimp7fF3TZS/7H99kNo5Ej+PFQ4UaOmgFQdDq17qj8HtQvlmKku8JVKtqQmR/HOB3ipBu7fU3dXwbs0xMSUeMplM8LpU18XvHu/J/67VcYITT35JuUmRHNuEp1bKMfWujri9Q7hJEx37X+olmoqCrd9h/z+npLaHl0qOOTUnb/Y19hjxttDc5u4nk8ZIZedFA7B9rBNHUchl9dK19/17EwTfoe8fvwXZ7w3FoWlDIJfLwDEXZUGQI5fDW218XK3V4c6ESL6nj1iLJj58AMUS91BVME1Q7Evs9+q7bWcACEeGf+k2Y61kpVa4TrZDgkapwLS7OgimdjHciEplo9inxNvKYv+vVAqZYF+xARBnpR7UgA1ubc0uG7A3bVLnY8A9u42z6KzpQIffHYqcmal2vScy0Bu/PtUbrUL9BBc4pVx/Ef/y4W6YPryTSVFvbT3VP46/o+oeE4y5D3eHWikX9JqydMfHfuF8NEpEBnpjwRO9bBphk72wGwr9nh1obPeOCfHB1LuNvcIW7sixuD6V0vREwn7x2a6NJeXVgmCiTGJW8JE9oqCQy5DaOVLwPPsllgpy2LvAmBD2jtS4r57qH4esaSnoFqM/cU++Qz8QobgYkb1bEnchr03WwJIYC4M3NhZ1rXuqzf9JXcZgqQ1xBlhqvrVHEs1nFljiugz2xk4hkwl6KRmCK8P3zlDbplbKoZAJl5Mzj3Wcvvns92eTsG7CANwWHyYYLys+IgBrXxlosm1v3dkBQT4qfpDPk/nGekF2nC62w4PU/8WIHsb9UyEa3+v/mAFEvVRyjOnbErMf6Mo/ZxjAVTrIYYIuk1eN2HO9uPaQPT/Y2oGqvEqL5wa2wou3tbHYA1bKqF7C4+LV29vCV61wyGCD9YWGWHWCw++kYN6mU3jYxhOHwRZm7g/DAXxn57oV5NlKmNpVIqF5IA5dLML80T1Ey7GD9dl3+LBfUMP8WAnNAzGqVzTC/L0gk8kE7d63xlsezEvqzsJc10nDndrkYfGYv+U0Xh8qLCKUyYA24f7Y+2ayRA8q4+eI78IVcplgpvQEZkCtU5dLBcuymaNHesWgQ6Q/OkQKxwESFB4zQY5cZtucYl2iAnHAxgEq6/la65a6RQfVaWLBrFz7BwO1dBdfWzkzU9Fi0kqT55c809vkTntIh3As23dRtE22fc7AtqF8M3D/Nk2RV1SOEzWdD+Ry4cVdfKGPDvHBT0/2gq9GiTf/OiR4je3dZCCTydA6zB/fj7kFAPi/z9BM0yEyAEeYXkTxEQHYN+V2vr5oP9P0zWID27u7mJ5f2e/4I71j8dbyLPRv0xQAoGG+w1IB7s7T+g4L4QFeuFIqrDsUBDlmjoHJw+IFY6+JMzm5Ng4oyDpz5QYmiXqK2apPqyYYP7g12tRkzl64rQ2eH9QK3d5N45cpr5Luqu8uKJPjBL4aJV5LaSdI9dmCzS50Zdqv6wN7IvRSyfHvC/1wZsYdgiYz/XKW62AsuVJqHIzQcKKXyWSYcV9nvHK7MUX8x3NJeCwpFt891tNkHayKKtMpHwpKhAMeGqbAMBTiPTuwFfa9dTvfHdIwKmrmW0MAAMG+apOTs7nB/gzYu1m23d5SkbhcLkOP2BCTu2w2cGRrcsylim1hbtRoRxWbNmTWmpOtEU9XYosqUROIrSPz1kYiUzhrIJWp/W33OUyvGRvm45FdBK+xQQEbaHdsFmiSnVVaCHIAfccEqe+FRIxjliFAWD6uL966swM2vjaIWY/xM1szTUu3dzCOqcb2vBrbT99xZO7DxrHGNArj3/RwrxgsfqY3f7PXjunBJTUmjOH8+PlDXdEjNhg/ss34zO7QMlHlt8x57un+caKibOHYX9dumHbYsKYumUOZTIbXUtrhHqa3rlIhx/+Yom5bJ3R2FQpy3MjnD3Xlf6/vdns10/RjqJaXugCwdSHmCtHMYduTLb23Z4sQvHtPJ6sXYVvmqnqyX0vkzExFuJlRVb97/BbkzEw16TbPYpsR/b2UJr0a7u8RhfaRAXx2aN4j3fFo7xiLE5aaw1402CEEbA1y2K7pBub2oqKOF3hPUNfzc23mG2OPpw6RAXigZ3St3ltbUlkTQD+g3uF3UnB/jyi+gBYQXljF2Ygn+xl7lypkMsiYo82ei5+5bZJiCBDUSjnG9mtptis3mxVmMyCB3qbDPHRmMrDsuVchl6F3XBO+lqWJnwabJw5CxpRkyc/sUzOKe+swf/z5fB8MZJrx2X3Hjg2U3D4MD/SMwuRh8ZDLhZmbkwWlgrqvEcw5xdrX99MHu6BtuB8y377d8oK18Hgf49g6Peu5t6C9qLnKjcSF+uHf8f2sTgPgbJYKXNVmKv1tcZnJspgLOuwxf3RP3PrJJsy8rzP/XL/WTbHt5BWHjjGTfuoq/7ufRmmSam7ipxHMVzMsIRLDEoR1PbZiaxr8RV1JbVEhMaGpuWxFfdeGuKO67gN7ejwasAW1Uk1Xj/aOwS87zwHQD0C3IbuA7/2UGBeCTdn6Qe5eG9IWH689bveM01LxRGLNRddw8/FAz2jsrRkvhsVmgaJDfJDY0jRTZGDPvrVlWZlMX0DcOSrIpnW+ntIOJeVVeCQxFr/vPY/DufqmLfb/zHBTERXsjf5tmtZMZGw54IptYhpULXziFvywPQfv3Ztg9n3sxMwD24ZiwRO3oG24P2QyGWbdb8yeiXtAsYHtVaYJzNoeu7dbFO7t5pwemT5qJQ5OG4KbFVq763zqG2Vy3ExCVCBimtR/QSibnbBULX/umnFSTXubq04xXbodIdRfg0PTUgSFwp880AWjekXjl7GOGxGbTa8rFXKTJjFHkpvpaWXPna6t7Bn63RHqOqePM7DZwtG9Y+1qJlYpZPz4TJaIi8vZoFOr4wT1MGzvHQCY/WBXwcWMHXV2/K1tcOL9YXzPxmdFg9eZIzVyrVwUCI/sKX1xZJ9/6JZoQZOqRiUXBG2WstHijKMtQc7qlwZgYko7m7vghwV4Yf7onhjQNlTQCSHMX4OmfmoE+6j44n6ZTIafxyZi7sPda9WEOahdGH56spegu7bBxyO7oG/rJniJmb9PJpNhcLswyeXZ82+wjwoqhZyvCXqkd4xgHa4U4KVyWIcYZ6JMDgEA7DpjzFZYOuH0bd2U/91SLyxXCQ/wwgwms+MIj/dpgS83nkSHmou0M4MD9qLLNlEp5TKLPTIMAr1VKCozHfPk2QFxmL/ltOA5R4yNYg83nNZG1KvH9v/XYB8Vdky6DW8uP2R9WQtZRS3HIeOscXTtF25rg93MdxEQ9sQRfzXZY2TyHe3x6pB2uHqjAo98twuPieZzMmjZ1BefPtgFTXw1eOqnvais1qFfm6aCZdgLaIivmm+y4jgIepAqFcCKF/SjuHupFIgO8UGHyAD4aZR29Txjm7YeS5Le7nYR/oKaGHs82bcF9p29jpE9o6BSyLFp4mAoZLJ6yWbe3yPKrqZrNjg03GD9PDYRHMcJ/l8oEWsbNzztEFew9eQR5q+BXKYPhOwdkXf9q/pun0NFxczuLtRfg8y3b8df4/oAcG5GwlydjFIhgy3XYKmpMGQAJg2Lx67/u03QFFrfQYfMaoLdMawNeCaskzA+z8F6nYNBgLcK3mqF1c+Sy4TbY6jZMNDpOJSUV/OP7+ocabKfDNmllk19TTIuYmqlHJGB3tjw6iCM6Wt+NPZ7u0VhQNtQrJ8wEDPvS8DT/c1ngT57sCv/u9R8ap2aB/LTIyjkMqx4oR+WPNvbrkwDG9xb61VZG2EBXvj9uSSMrKl/8tPYN7hqfWIDL3YfiventWOB6FEmhwCwPSujVMhxYOoQ6Dj7i6NbhfrZPY6Qu2DHrFn1Yj8cuFAkmCvHUczdWV64XgbOhlyO1HlPJtOfIMV1UBGB3sgvdl7Tm9R2OIOXSo5ypqedeC91ah6ArIvCCQsHtQvFpuzLuKdrM/ydqR/VleM4sxcOb5VCMNaKIRi1lq2Qy2R4Obkt/snMhUal4AeDNNBxwLv3dMLiPefRsqkvZDKZyX769IGu+GzdcTzUKwZKuQw/7zyL5PbhcIToEB881Et6qItFTyfi3NWbGNA2FOsmDEBlNWcy552U2vTaYzM5zmiabags7UkKcmxDQQ4BALyS3Babsi+bDP4kxZYTnSeTyWRO6+JfUV23MSeksyXSTTLNg7xw4HydPs4tqBTCIEfHNCcq5TIEeQszjjIAX4zqhvPXygTD33Oc+SYAfy+lIMgZV9OF1lqgL5frg8vD7w6VfF2r46BWygXBv/jiFeyrxjv3dOIfH5w2BP616NVlrz6tmqJPTflL6zDHB/QsS+PrNGZScczElHb4aE02Zt3v2GZ5T0UhMwEAdIkOwr63bsf7wztZX5g4TQtRzw1DtmjW/Z1taq46f/2mxdfZddRX8xH/ebW485w1wvqJvJmokJbtbdRFIhgd0jEC/l4qQYAD6PeNYSRqS3q1DOG78loNcqz8yVJ1QGP6tAAgPVAdoC/4dHXRaV0NqRm3xtCVWiloonHJJrklqWzNuMGtcWjaENxl5vggQhTkEF6Ir5oGiHOxJ/q2xF1dmuHDEfquqKtfHoCT7w+zeSwVa3OcsYOQ1ffFpDYfd0830xO5uHeNuLbiTWbo/Ykp7dCiqbC34oNm9qWO49A9Jthql+xopqeOWmG5mddak4JUkNOiqS+OTR8qGDfL03w0sgvev7cT5j+qH2SPsjfSdp6+Kvl8Y8+m24OCHELciFqpn8TzQWaWYsPIqrZMyNdEoicPe50Vz1xcn2oTVGmUChx+J0Xw3Iu3thE8Fq+XHb8pwEuFiSnxGNE9CnGhvpg/uofZQN6wZ+zpFWQ9k2P5j9aaDmsEQF8j19CzNZYEeqvwSGIs3/OMrcPx3L/afs2D7Rs1n5iiIIeQBuIpCz1gLGEvxM/Uch2OMKxT7XrVsaNjR4d4Wx2sjY0NlAr9fGifPNAFG14dZDJNCcsQ/1kLJtn1W+td1cFKTzxnzGPVEFEmR+jHJ3shKa4J/ny+j6s3pcGjwmNCGohoO2cNf6BnFH7fewHjmXlmXh3SFgcuFKJnbAiOF5ifd6l3XAg/2aCjJDQPsmt5qe7wQd5qqxkO9lV7Lp6G4MaewMNcJufrR7tjU/ZlvJTcRvJ1A209D8jortx9/qP6NrBtqGCoA1J7Tsnk5OTkYOzYsWjZsiW8vb3RqlUrTJ06FZWVwsnFDh48iP79+8PLywvR0dGYNWuWybqWLl2K+Ph4eHl5ISEhAatWrXLGJhPicT4c0RlZ76QIxvUxjOxq7eLrDPa2vmyZONjkOZsCEDaTY8fFc9eZazWfYfNbzAY5sU18MXNEZ8nRhVk6CnIACLudU1aHOJJTgpxjx45Bp9Nh/vz5OHz4MD799FN8/fXX+L//+z9+meLiYgwZMgSxsbHIyMjARx99hGnTpuGbb77hl9mxYwdGjRqFsWPHYv/+/Rg+fDiGDx+OrKwsZ2w2IQ0eG0jIZLJaTSDpLqRGCrYtxqndBfNizSSOUoEU+wy7RnPNVTE2Zt1cWSPlbtqF+8NPo+QHFiTEEZxyBhw6dCiGDjWOCxEXF4fs7GzMmzcPH3/8MQDg119/RWVlJX744Qeo1Wp07NgRmZmZmD17Np555hkAwOeff46hQ4di4sSJAIDp06cjLS0NX375Jb7++mtnbDohxI1JhQTiMEYuyOTYfh93e023ZnviDnGRct/W+hmrfWwcTZeaq4xWvtgP1TrOLaeLIQ1XvRUeFxUVISTEOCV7eno6BgwYALXaeLeWkpKC7OxsXL9+nV8mOVk4pX1KSgrS09MtflZFRQWKi4sFP4QQEQvXV2eMOmsueOgWE2Ty3OqX+5suaIa4Rqe28/uEB+i739tTk6MSBTm/PtUb3z7W0+aeUVR4bKRUyCnAIQ5XL0HOyZMnMWfOHDz77LP8c3l5eQgPFw5Nbnicl5dncRnD6+bMmDEDgYGB/E90tG1jjBDS0LUSzWJtrxZNfLBgzC14556OJq9NTGlXp3Wb/0xfk+fiI6R7JUn1fBrbTzhHU10rOtiPeCRR35X/tSFtJZe1d2oTMcrkEOJcdn1DJ02aVDO3ivmfY8eOCd5z8eJFDB06FCNHjsTTTz/t0I03Z/LkySgqKuJ/zp/3gLHrCbGBuZmnbbVp4mAMjg9Dq1A/bKiZUNVgHNNLqzZsmXvLXm/f2YEfPdegrsPLsNv53vBO2DHpVsG4RWytkFJRtw+jGIcQ57KrJufVV1/FmDFjLC4TF2cchyM3NxeDBw9Gnz59BAXFABAREYH8/HzBc4bHERERFpcxvG6ORqOBRmN55FdCPJHSgVOLO3wwOokLerDErNZS+rRqgh2nruJRURAXFextMuAeZ65K2EaCqS9kMjQL0veQ+uzBrvg78yLG32oM9jRKal4hxJ3ZFeSEhoYiNNS2vvsXL17E4MGD0aNHDyxYsAByURt/UlIS3nzzTVRVVUGl0p/o0tLS0K5dOwQHB/PLrF+/Hi+//DL/vrS0NCQlJdmz2YQ0GrYGDYBzMivWxEf441iecXwemUxm00jOP4y5BcfzS5Ag0fPGUbEYPxigmdeHd2uO4d2aC55ryL3XCGkMnFKTc/HiRQwaNAgxMTH4+OOPcfnyZeTl5QlqaR5++GGo1WqMHTsWhw8fxpIlS/D5559jwoQJ/DIvvfQSVq9ejU8++QTHjh3DtGnTsHfvXowfP94Zm01IgxcX6oepd3VwyLxHzhitZMETt1j8jHu7NceipxNN3uelUqBzVJBkdqk+Mk7m1LYm5+ORXQAAsx/oUqv3E0Js45TbkLS0NJw8eRInT55EVFSU4DXDXVtgYCDWrl2LcePGoUePHmjatCnefvttvvs4APTp0weLFi3ClClT8H//939o06YNli9fjk6daKZsQli3xYfxM2M/0bellaVdgwMQGeiN2zuEI+2IvhlaJhMGKZ8+2NWudRqakhxBqgDa+nuM4+GESIzrY879PaKQmhBpMrkoIcSxnBLkjBkzxmrtDgB07twZW7dutbjMyJEjMXLkSAdtGSGe6fsxt1hfSMSRvZe3vj4Y/WdtrMU7a5eF+WVsIs5fv+mQgeN+GZuITdkFGNO3BQD7mvFkMhn+e6k/Zqcdx6tmemCZQwEOIc5HDcqEEKu+frS72dfUSrlN82pJBVVv39UBG48V2L09/do0tfs9ltbFrs/e4K99ZAC+faynw7aHEOI4NAs5IY2UPRfzpDjzQYUhF7N90q0I9La98BkA7u7SzK7lbcEOKOertv8+rkMzyzOHE0IaDsrkENLARQZ61ep91ppl7K3nbR7kjTZhfth79nqtPs9RvNUKLHziFnAAfGvR+2liSjv4apS4s3Ok4zeOEFKvKMghpJGKCrZtEkkAFktn7A1d6mMmg0Htwmr9Xn8vFd4YGu/ArSGEuAoFOYQ0Ui8nt0FJeRXWHM5HUVmVqzeHEEIcjmpyCGngxHM32crfS4VZ93dBcvtwq8taarqytVWL5qIkhNQ3yuQQ0sDVNsgxeH1oOxy9VIxHescInpfZGL7YGrxQjEMIqW8U5BDSwNV1xN/wAC+seql/rd+f0snyXHKEEOIq1FxFCLHKUhj1wb3GEcgtZWtsmaOKEEIciYIcQogkWxNE/l62jY3DSfwGAEF2TCpKCCH2oOYqQohVDp8EkzFucGtkXSzCAz2j67yux5JiHbBFhBBPQUEOIcRhLDVJaWpm7BYv0tRPg6XP9XHI58c1tX+STUKI56LmKkIakJdua+OSz61rHqdzVCB6t2zikG2xxJkZJ0JIw0NBDiENyMvJbbD19cGu3gyzzOVx5ozqBrlcH4CE+Kqd9vn2zp1FCPFs1FxFSAMik8lsmvHbnb0xLB55xeUY1SvG+sI2mjWiM/afL6T5pgghAhTkEEKsMrQChQdokF9cUad1NfXT4OexiQ7YKqMHbonGA7fUvXCZEOJZqLmKEGKzX8YmCqaBEBf6mqs7piFyCCGuQJkcQojN2oT747vHe+L6jUocyytB77gQV28SIYSYRUEOIUSSpY5Kwb5qJLUy7S1FnZsIIe7E6c1VFRUV6Nq1K2QyGTIzMwWvHTx4EP3794eXlxeio6Mxa9Ysk/cvXboU8fHx8PLyQkJCAlatWuXsTSaEiNg8Cae55irHbQohhNjM6UHO66+/jmbNmpk8X1xcjCFDhiA2NhYZGRn46KOPMG3aNHzzzTf8Mjt27MCoUaMwduxY7N+/H8OHD8fw4cORlZXl7M0mxK2ldNTXxdwaH+a0z6AxZwghDZ1Tg5z//vsPa9euxccff2zy2q+//orKykr88MMP6NixIx566CG8+OKLmD17Nr/M559/jqFDh2LixIlo3749pk+fju7du+PLL7905mYT4vY+eaArPn9I/0MIIUSa04Kc/Px8PP300/j555/h42M6rkd6ejoGDBgAtdo4MFhKSgqys7Nx/fp1fpnk5GTB+1JSUpCenm7xsysqKlBcXCz4IcST+GmUuKdrc5snx6wrW5ubzC1HM5ATQlzBKUEOx3EYM2YMnnvuOfTs2VNymby8PISHhwueMzzOy8uzuIzhdXNmzJiBwMBA/ic6msbPIMSVKMQhhLiCXUHOpEmTIJPJLP4cO3YMc+bMQUlJCSZPnuys7bZo8uTJKCoq4n/Onz/vku0ghOip5DQkFyGk/tnVhfzVV1/FmDFjLC4TFxeHDRs2ID09HRqNRvBaz5498cgjj+DHH39EREQE8vPzBa8bHkdERPD/Si1jeN0cjUZj8tmEEPsYZg0HAIWNRchSzVIP9IxCTJOGPRUFIaRhsivICQ0NRWhoqNXlvvjiC7z33nv849zcXKSkpGDJkiVITNQP556UlIQ333wTVVVVUKn0dQVpaWlo164dgoOD+WXWr1+Pl19+mV9XWloakpKS7NlsQkgtNPXT4OXkNlAp5PBWK2x6T0rHCBy8UCR4btb9XZyxeYQQYpVTBgOMiRFOvOfn5wcAaNWqFaKiogAADz/8MN555x2MHTsWb7zxBrKysvD555/j008/5d/30ksvYeDAgfjkk0+QmpqKxYsXY+/evYJu5oQQ53k5ua1dyz8zIA4tm/piw7EC/JFxwUlbRQghtnFZQ3lgYCDWrl2LM2fOoEePHnj11Vfx9ttv45lnnuGX6dOnDxYtWoRvvvkGXbp0wR9//IHly5ejU6dOrtpsQogFKoUcdyREoomf2vrChBDiZPUyrUOLFi0k2+o7d+6MrVu3WnzvyJEjMXLkSGdtGiHEGag7FSHEDVCXB0IIIYR4JApyCCGEEOKRKMghhBBCiEeiIIcQQgghHomCHEIIIYR4JApyCCEOR52rCCHugIIcQgghhHgkCnIIIYQQ4pEoyCGEOFxUsLerN4EQQupnxGNCSOMyqlcMzl29if5trU/oSwghzkJBDiHE4VQKOabc2cHVm0EIaeSouYoQQgghHomCHEIIIYR4JApyCCGEEOKRKMghhBBCiEdqFIXHHKcff7W4uNjFW0IIIYQQWxmu24bruL0aRZBTUlICAIiOjnbxlhBCCCHEXiUlJQgMDLT7fTKutuFRA6LT6ZCbmwt/f3/IZDKHrbe4uBjR0dE4f/48AgICHLZeT0T7yja0n2xD+8l2tK9sQ/vJdvW5rziOQ0lJCZo1awa53P4Km0aRyZHL5YiKinLa+gMCAuhLYSPaV7ah/WQb2k+2o31lG9pPtquvfVWbDI4BFR4TQgghxCNRkEMIIYQQj0RBTh1oNBpMnToVGo3G1Zvi9mhf2Yb2k21oP9mO9pVtaD/ZriHtq0ZReEwIIYSQxocyOYQQQgjxSBTkEEIIIcQjUZBDCCGEEI9EQQ4hhBBCPJJHBTlbtmzBXXfdhWbNmkEmk2H58uWC1/Pz8zFmzBg0a9YMPj4+GDp0KE6cOGF1vQcPHkT//v3h5eWF6OhozJo1S/D6smXL0LNnTwQFBcHX1xddu3bFzz//bHW9586dQ2pqKnx8fBAWFoaJEyeiurpasMymTZvQvXt3aDQatG7dGgsXLrS6Xmsa0n66evUqhg4dimbNmkGj0SA6Ohrjx483mYfMGfsJcN2+Yi1evBgymQzDhw+3ut7GdkyxbN1PrjymXLWfFi5cCJlMJvjx8vKyul5XHU9Aw9pXjfGYAoDCwkKMGzcOkZGR0Gg0aNu2LVatWmVxvdeuXcMjjzyCgIAABAUFYezYsSgtLbX7s23GeZBVq1Zxb775Jrds2TIOAPfXX3/xr+l0Oq53795c//79ud27d3PHjh3jnnnmGS4mJoYrLS01u86ioiIuPDyce+SRR7isrCzut99+47y9vbn58+fzy2zcuJFbtmwZd+TIEe7kyZPcZ599xikUCm716tVm11tdXc116tSJS05O5vbv38+tWrWKa9q0KTd58mR+mdOnT3M+Pj7chAkTuCNHjnBz5syxul5P20/Xrl3jvvrqK27Pnj1cTk4Ot27dOq5du3bcqFGjnL6fOM51+8rgzJkzXPPmzbn+/ftz99xzj8VtbYzHlIE9+8mVx5Sr9tOCBQu4gIAA7tKlS/xPXl6exW115fHEcQ1rXzXGY6qiooLr2bMnd8cdd3Dbtm3jzpw5w23atInLzMy0uL1Dhw7lunTpwu3cuZPbunUr17p1a8F+sud7bwuPCnJY4v/s7OxsDgCXlZXFP6fVarnQ0FDu22+/Nbuer776igsODuYqKir459544w2uXbt2Fj+/W7du3JQpU8y+vmrVKk4ulwu+PPPmzeMCAgL4z3r99de5jh07Ct734IMPcikpKRY/2x7uvp+kfP7551xUVBT/uD72E8fV/76qrq7m+vTpw3333Xfc448/bvXi3ViPKXv3kxRXHFP1uZ8WLFjABQYG2rV97nI8cZz77yspnn5MzZs3j4uLi+MqKytt3r4jR45wALg9e/bwz/3333+cTCbjLl68aPNn28OjmqssqaioAABB2lEul0Oj0WDbtm38c2PGjMGgQYP4x+np6RgwYADUajX/XEpKCrKzs3H9+nWTz+E4DuvXr0d2djYGDBjAPz9t2jS0aNFCsN6EhASEh4cL1ltcXIzDhw/zyyQnJwvWn5KSgvT0dDv/etu5234Sy83NxbJlyzBw4EDBZ9f3fgKcv6/effddhIWFYezYsZKfT8eUnr37Scxdjiln76fS0lLExsYiOjoa99xzD39MGDSU4wlwv30l1hiOqX/++QdJSUkYN24cwsPD0alTJ3zwwQfQarX8ewxNf+x6g4KC0LNnT/655ORkyOVy7Nq1y+bPtkejCXLi4+MRExODyZMn4/r166isrMSHH36ICxcu4NKlS/xykZGRiImJ4R/n5eUJvuQA+Md5eXn8c0VFRfDz84NarUZqairmzJmD22+/nX+9adOmaNWqlV3rNbdMcXExysrKarUfrHG3/WQwatQo+Pj4oHnz5ggICMB3331n9bOduZ8A5+6rbdu24fvvv8e3335r9vPpmKrdfjJwt2PKmfupXbt2+OGHH/D333/jl19+gU6nQ58+fXDhwgX+PQ3leALcb18ZNKZj6vTp0/jjjz+g1WqxatUqvPXWW/jkk0/w3nvv8e8JDAxEu3btBOsNCwsTrFepVCIkJMTqMcV+tj0aTZCjUqmwbNkyHD9+HCEhIfDx8cHGjRsxbNgwwfTtM2bMwE8//WT3+v39/ZGZmYk9e/bg/fffx4QJE7Bp0yb+9fHjx2P9+vWO+FOcyl3306effop9+/bh77//xqlTpzBhwoRa/X2O5Kx9VVJSgtGjR+Pbb79F06ZNzS7X2I+puu4ndzumnPndS0pKwmOPPYauXbti4MCBWLZsGUJDQzF//nx+mYZyPAHuu68a0zGl0+kQFhaGb775Bj169MCDDz6IN998E19//TW/zL333otjx4457O+pDaVLP72e9ejRA5mZmSgqKkJlZSVCQ0ORmJgoSJ2JRUREID8/X/Cc4XFERAT/nFwuR+vWrQEAXbt2xdGjRzFjxgxBClC83t27d1tcr7nPDggIgLe3tw1/ce24035i1x8REYH4+HiEhISgf//+eOuttxAZGemy/QQ4Z1+dOnUKOTk5uOuuu/jXdTodAP1dT3Z2tuRdZGM7pmq7n9j1u9sx5czvHkulUqFbt244efKkxfW66/EEuNe+YtffWI6pyMhIqFQqKBQKfpn27dsjLy8PlZWVguYmdr0FBQWC56qrq3Ht2jWrxxT72fZoNJkcVmBgIEJDQ3HixAns3bsX99xzj9llk5KSsGXLFlRVVfHPpaWloV27dggODjb7Pp1Ox7eHmlvvoUOHBP/haWlpCAgIQIcOHfhlxHcLaWlpSEpKsvo3OoI77Cdz7wGM7c2u3k+AY/dVfHw8Dh06hMzMTP7n7rvvxuDBg5GZmYno6Giz621Mx1Rt95MUdzumnP3d02q1OHToECIjIy2u192PJ8A99pUUTz+m+vbti5MnT/J/JwAcP34ckZGRkgGOYb2FhYXIyMjgn9uwYQN0Oh0SExNt/my71Kpc2U2VlJRw+/fv5/bv388B4GbPns3t37+fO3v2LMdxHPf7779zGzdu5E6dOsUtX76ci42N5e677z7BOiZNmsSNHj2af1xYWMiFh4dzo0eP5rKysrjFixdzPj4+gu5sH3zwAbd27Vru1KlT3JEjR7iPP/6YUyqVgur1OXPmcLfeeiv/2NA9c8iQIVxmZia3evVqLjQ0VLJ75sSJE7mjR49yc+fOdUiXw4a0n1auXMn98MMP3KFDh7gzZ85wK1as4Nq3b8/17dvX6fvJlftKTKrXEB1TpmzZT648ply1n9555x1uzZo13KlTp7iMjAzuoYce4ry8vLjDhw+b3U+uPJ4a2r5qjMfUuXPnOH9/f278+PFcdnY2t2LFCi4sLIx77733+GWWLVtm0itq6NChXLdu3bhdu3Zx27Zt49q0aSPoQl6b770lHhXkbNy4kQNg8vP4449zHGfs0qdSqbiYmBhuypQpgm5qHKc/SQ4cOFDw3IEDB7h+/fpxGo2Ga968OTdz5kzB62+++SbXunVrzsvLiwsODuaSkpK4xYsXC5aZOnUqFxsbK3guJyeHGzZsGOft7c01bdqUe/XVV7mqqiqTv6lr166cWq3m4uLiuAULFtR6/7DrbCj7acOGDVxSUhIXGBjIeXl5cW3atOHeeOMN7vr16yZ/k6P3k2G9rthXYlIXbzqmTNmyn1x5TLlqP7388stcTEwMp1arufDwcO6OO+7g9u3bJ1jGnY4nw3obyr5qjMcUx3Hcjh07uMTERE6j0XBxcXHc+++/z1VXV/OvL1iwgBPnUq5evcqNGjWK8/Pz4wICArgnnniCKykpsfuzbSXjOI6zP/9DCCGEEOLeGmVNDiGEEEI8HwU5hBBCCPFIFOQQQgghxCNRkEMIIYQQj0RBDiGEEEI8EgU5hBBCCPFIFOQQQgghxCNRkEMIIYQQj0RBDiGEEEI8EgU5hBBCCPFIFOQQQgghxCNRkEMIIYQQj/T/5cZiMWyeHz0AAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 640x480 with 2 Axes>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"def normalise_dataset(data):\n",
" # combine the three channels into one by normalisation\n",
" # not sure if this is \"OK\" from a scientific standpoint, \n",
" # but it seems to work\n",
" data = np.linalg.norm(data, axis=1)\n",
" \n",
" # run a high pass filter, dont care about frequencies,\n",
" # just enough to remove the low frequency noise\n",
" b, a = butter(5, 0.01, btype='high', analog=False)\n",
" data = filtfilt(b, a, data)\n",
" return data\n",
"\n",
"plt.subplot(2,1,1)\n",
"plt.plot(x1, normalise_dataset(event1))\n",
"\n",
"plt.subplot(2,1,2)\n",
"plt.plot(x2, normalise_dataset(event2))"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.13"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment