{ "cells": [ { "cell_type": "markdown", "id": "cbce6e72-320a-4c2d-afe8-1dcbebbd079e", "metadata": {}, "source": [ "# Real case: Mannsworth - MASW2D \n", "\n", "This notebook implements a 2D MASW (Multichannel Analysis of Surface Waves) workflow based on a moving window approach. \n", "Data comes from [Barone et al. (2021)](https://doi.org/10.1007/s10712-021-09631-x)" ] }, { "cell_type": "code", "execution_count": 1, "id": "55253c08-ed34-47b6-bf3b-9ab77d623503", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/home/alberto/anaconda3/envs/pyswapp/lib/python3.9/site-packages/obspy/core/util/base.py:26: UserWarning: pkg_resources is deprecated as an API. See https://setuptools.pypa.io/en/latest/pkg_resources.html. The pkg_resources package is slated for removal as early as 2025-11-30. Refrain from using this package or pin to Setuptools<81.\n", " import pkg_resources\n" ] } ], "source": [ "from pyswapp import *" ] }, { "cell_type": "code", "execution_count": 2, "id": "76653f4f-2c4b-4f6b-857d-e3300682eed3", "metadata": {}, "outputs": [], "source": [ "# directories\n", "prj_dir = '../../data/real_data' # project directory\n", "path2raw = os.path.join(prj_dir,'raw') # path to raw data\n", "path2geom = f'{prj_dir}/Geometry_mannsworth.csv' # path to the geometry file" ] }, { "cell_type": "code", "execution_count": 3, "id": "0e03faf9-ca1d-45e1-aa99-481c5f586b9b", "metadata": {}, "outputs": [], "source": [ "# processing and plotting settings\n", "settings = create_settings(fmin=5, fmax=30, # frequency range\n", " vmin=10, vmax=500, velstep=1) # testing phase velocity range and step" ] }, { "cell_type": "code", "execution_count": 4, "id": "1ed7b696-c51e-4253-8fb6-f38b14ba7a37", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Warning: Ignoring XDG_SESSION_TYPE=wayland on Gnome. Use QT_QPA_PLATFORM=wayland to run on Wayland anyway.\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "Creating project:\n", "Copying raw data into project directory ..... 0.03 s\n", "Copying geometry into project directory ..... 0.0 s\n", "Reading seismic shot files ..... 14.85 s\n", "Read 82 files and applied geometry and settings.\n", "\n", "Existing procsets: raw proc1\n", "Active procset label: proc1\n", "\n" ] } ], "source": [ "swam = MASW2DManager(f'{prj_dir}/proc/masw2D', # project directory path\n", " path2raw=path2raw, # optional, copy & per default rename data from path (outside project directory)\n", " path2geom=path2geom, # optional, copy geom from path (outside project directory)\n", " settings=settings, # optional, define settings for visualisations and processing\n", " rename = False, # optional, rename copied raw data to preferred filename format (Shotfile_), default value is False!\n", " overwrite = False, # optional, overwrite database .db file if it already exists --> create new project data base\n", " )" ] }, { "cell_type": "code", "execution_count": 5, "id": "333672a3-8c99-4896-9c0d-a48981cc7d9d", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Moving window along (SIN,REP) = (43, 2) ..... 23.39 s\n", "Applying check_traces to (SIN,REP) = (43, 2) ..... 26.29 s\n" ] } ], "source": [ "# offset in meters\n", "dx = 2\n", "min_offset = 2*dx\n", "max_offset = 10*dx\n", "w_len = 22\n", "w_m = 2\n", "\n", "# Run moving window along data\n", "swam.moving_window(minoffset = min_offset, maxoffset = max_offset, wlen = w_len, wmove = w_m)\n", "# remove traces of zero amplitude\n", "swam.preprocess(type='check_traces')" ] }, { "cell_type": "code", "execution_count": 6, "id": "2f4ddea9-0ab0-4a10-82dc-12fada236b38", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Applying fdbf transformation to (SIN,REP) = (43, 2) ..... 251.79 s\n" ] } ], "source": [ "# transform\n", "tr_type = 'fdbf'\n", "swam.transform(method = tr_type)" ] }, { "cell_type": "markdown", "id": "90d8fa9d-cc71-46f1-8789-1e073c45234d", "metadata": {}, "source": [ "you can pick either in `fdbf` of `FK`" ] }, { "cell_type": "code", "execution_count": 8, "id": "ae60c595-a876-4cc2-a06c-85d4097447b3", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Loading the GUI ..... 6.01 s\n" ] } ], "source": [ "# manual dc picking\n", "swam.gui_interact('pick',\n", " use_windows=True,\n", " cmap = 'jet')" ] }, { "cell_type": "code", "execution_count": 25, "id": "477e7822-fa9f-4ad2-9530-9ec7cd740bc1", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Applying smooth to dispersion curve of (SIN,REP) = (43, 2) ..... 9.76 s\n" ] } ], "source": [ "# smooth the DCs\n", "swam.process_curves(type = 'smooth', method = 'FK')" ] }, { "cell_type": "code", "execution_count": 27, "id": "f6509a33-ef3e-416d-b331-46a9a1ccaf15", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAxUAAAEiCAYAAACP/f82AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAAC8oklEQVR4nOydeXwURfbAv13dcyQh3Ge470MuCYcCIsh9CqigEGURj3U55FTxxJ+u4gUoroquC6uI166KiiviKigrui54rBoBBWSR+zYhmZnu6t8fM9OZIQkzCQNDOv310x+Lnkr1e/PqvanquhTTNE0cHBwcHBwcHBwcHBxKiUi2AA4ODg4ODg4ODg4OZRunU+Hg4ODg4ODg4ODgcFo4nQoHBwcHBwcHBwcHh9PC6VQ4ODg4ODg4ODg4OJwWTqfCwcHBwcHBwcHBweG0cDoVDg4ODg4ODg4ODg6nhdOpcHBwcHBwcHBwcHA4LZxOhYODg4ODg4ODg4PDaaElW4CzwcGDB1m9ejWNGjUiJSUl2eI4ODg4ODg4ODgkkLy8PHbs2MHAgQOpXr36KfPu3LmTgwcPJuzZ1atXp0GDBgkrr6xSLjoVq1evJisrK9liODg4ODg4ODg4nEGWL1/O+PHji/18586dtG7ViBN5ZsKemZqaSnZ2drnvWJSLTkWjRo0AePHFF2nTpk1yhSln3JR5S7GfPb3x4dMuX0rJ3r17qV27NkI4s/nswrlg13jrbjz5TjdPScs6VzkX7FoUZfk7TQZZ9y4vdE/TBFde1IIBPTuVyLbj/vhSsZ+tuCPYMBy9pPg8b9xY0Hgc+toLxeZbNeYaAEYsLyx7mLcjXj4OfWdZ8WUN/10wz1t/LT7PyAlWesQHzxf/zAGTgmXFITvAyH8+V2y+t/peD8CVnz5VbJ5XLvoDANf/e2GxeQCe6zoDOHd9tiiys7PJysqy2nzFcfDgQU7kmfz1yZq0au4+7ef+uNXPhCn7OXjwoNOpSLYAZ4PwlKc2bdrQqVOnJEtTvqioVCn2M8cWDucy8dbdePKdbp6SluVQMpzvtGSkVl9b5P03/pvD7dM7l7CsdcV+Fv7uvXVi5wHwfP5RzHyetWvjK+u/a2Lmc3/1YXxlbf1HbLnikB3Auz0jZr7UPXVi5kk/UqvYPCc/s6wR7zT3Fs00OrRznfbzpClPuwy7cG53OxOMlI7h7YaUkoMHDzq2tRmOXe2JY1f7oqkKAzs1dGxrM+zssxIzYZdDkHLVqXCwJydOnEi2CA5nAMeu9sSxqz0RikKzOpWSLYbDGcCuPqsjCZjGaV869utwlZZy1ak41+cDOpQcIQQNGjRwbGszHLvaE8eu9sWvS/606lvHtjbDzj5rJmiUwnRGKizsV0tOgR2H78o7Ukr27dvn2NZmOHa1J45d7YumKoy8oKljW5thZ581MBN2OQQpFwu1HeyNruvJFsHhDODY1Z44drUnCgpVKniSLYbDGcCuPhseqUhEOQ5BytVIhR2H78o7Qgjq1q3r2NZmOHa1J45d7UvAkCz98AfHtjbDzj5rmGCYZgKuZGty7lCuRirsOHxnB/qLK4r9bI18/ZR/Gx6arVWrli2DXnnlnLCrUnbr0+n41JnknLCrw+mjFL7lUgWX9WiGlBIhBB1uLv4chG8en2Gl1UACW2Ra7LKkK77niaOxtxrVjsTXhFJ+SY0rX6IQSmwdVRFfe8jOPmtCQpZYO32KAspVp8LBwcHBwcHBwcEhUeshnDUVBSS12/nggw/SpUsX0tPTqVmzJiNHjmTz5s3F5r/xxhtRFIVFixaV6nl262U7BG1ap04dx7Y2w7GrPXHsal8ChuSVT7Y4trUZdvZZ3YRAAi69FH2KeNu/2dnZjBgxgkqVKpGens4FF1zAzp07rc99Ph9Tp06levXqpKWlMWLECHbt2nU6X8tpkdRasm7dOiZPnsznn3/OmjVr0HWdAQMGkJubWyjvW2+9xRdffEFGRvGnScbCmf5kP6SU/Prrr45tbYZjV3vi2NW+uFTBxH5tHNvaDDv7rETBSMAli5oPGIN42r8///wzPXv2pFWrVqxdu5ZvvvmGu+66C6/Xa+WZPn06b775Jq+88grr168nJyeHYcOGYRhGQr6jkpLU6U/vv/9+1L+XLl1KzZo12bhxI7169bLu//rrr0yZMoXVq1czdOjQsy2mwzmOpjmz+OyIY1d74tjVnpiYHMnxJVsMhzOAXX1WmsErEeWUlHjav3fccQdDhgzh4YcftvI1adLESh87doznn3+eF198kX79+gGwfPly6tevz4cffsjAgQNLoc3pcU6NZx07dgyAqlWrWveklFx99dXMmTOH8847L2YZe/bsYdOmTVFXdnZ2VB4ppdXrjjdtmsFaYxhG3GnTNAulgbjTJZXxXNRJqArulOCCN6EJK61qwpJFdam4vFpB2hNMa241ph4ANWrUQAjh2MlGOpmmSc2aNRFCJE0nza0C4PJoqK5Q2luQjszvTnEhNFGQVpVCOnnS3ChCKUiHXm6F8yhK8D6AIhQrLVQFT6rb0qM4fwqnpZTF+lNYp2TVPdM0qVWrFoqinFN1z5Pqtmx2sp3s4E+JjhGqUHBralRaN0ze/XIbYTSh4FKDPuFSRVQ6ys9C37tbU1FF4bRhGIiQs3g11UqnaKrVgInUL0UN1nclIi1CeSB48neKqlppbyitKkqUXB4RvK8pCu7Q1J/ItJQSVyjtFgItJFdkOtJOHlWghvVQC/SIlr1ApxRVtd5/R9opRRToVJBW8CqaZSdPKK2i4FZC+iGstJQSFwX3w2kNFY3oGCGEoEaNGpZdk1334vGneJGQoJGKINnZ2YXan3v27IlLlpPbv1JKVq1aRYsWLRg4cCA1a9akW7duvPXWW9bfbNy4kUAgwIABA6x7GRkZtG3bls8++6xE30WiOGc6FaZpMnPmTHr27Enbtm2t+w899BCapjFt2rS4ylmyZAmZmZlRV1ZWFgAHDhwAYN++fezbtw8IdkLC93ft2sXhw4cB2LlzJ0ePHgVg+/btHD9+HIBt27aRk5MDwNatW8nLywNg8+bN+HzBtzTZ2dnouo6UkuzsbKSU6LpudW58Pp81dy4vL4+tW7cCkJOTw7ZtwaB8/Phxtm/fDsDRo0etOXSHDx+25ssdOHDAqrDnqk7NuzVhxkvXAdCuTytuevZqADqP6GDp1PPKLoy7fyQA/a7ryei5gwEYOq1vTJ1++eUXtm3bhpTSsZPNdNqxY4elRzJ0GjrtEgBG3zaIfpN6ADDuvkvpeWXnQjrd9OzVtOvTCoAZL11H825NCuk0b81MajWuDsD8DXOpWDMdT5rb0qlizXTmb5gLQK3G1Zm3ZiYA9c+ry20rJ1s6FedPExeMsXQqzp+GTuubcDuVtO7t3LmT3Nzcc6ru3bZyMvXPq1ukneziT4mMEZ2a1uSWyzIB6NE6g6nDO+DWBHeP7cb//vc/AEZ0bMyE7kGfGNOlGWO6NANgQvdWUTpd0qEBAFNHdKBHm+D05jmXZ3J+05qWTu1qBhtajw/oTrMqFQFYMrQXGRXTonRKUTVe79afFFWjqtvL6936A1AvpYKlU/P0ijzVtScAHapU5bHMCwC4oHrNKDvNbd8BgKH1GjCjTbBNckWjJtzQopVlp6tbBHW6oXUrLm8a9Pnp7dsypGGDQna6t2sm3WoFG+cPd+9K+2pVC9npr936UC+1AgCv9+hPVbeHFFWNstNLzYcDUNedznNNBgHQzFuFxY37WXa6p1ZQ704p9bitZh8AeqY14ubqQb0PHz7M+IrBt+H9UztwWYXgdzAsLZNhaUG7huuelJItW7Zw6NChQjqda79PYb+Jl0R0KMIXQFZWVqH255IlS2LKUVT7d//+/eTk5DB//nwGDRrEBx98wKhRoxg9ejTr1q0DYO/evbjdbqpUqRJVXq1atdi7d2+JvotEoZjhLmSSmTx5MqtWrWL9+vXUq1cPCPbChg4dyqZNm6y1FI0aNWL69OlMnz69yHL27NlTqGeYnZ1NVlYWX3zxBV27drV6s+G32/GkFUWx3q4JIeJKA1ZPP5xWVdXqacdKm6ZZIhnPRZ0GusaguTX8eQGEJtBcKv68AKomeC//ZYQQDPJciVAVAvk6qktFCIWAT0dzq6w6seKUeui6zpEjR6hWrRpAUuw0rMJ4Syc1rJ9LRagK7514GSklw9LGWzoB6H4Dl0dDSpP3fa9YOg3UxuJOcWEEDAxd4k5xoQcMpC55P/BKual7gUCAo0ePUq1aNUuWs63T0NRxUXYyAgYur4Y0irZZ2E7uFBe6Xw/mC9lsgHYlnjQ3/rwApjSD6RN+TBPe9wfr+EBtDO5UN75cP4oIjvD5cv0IVcHlcfFuzktR9e1kf1JdKqtyVwTrW4XxRfoTwD/yo2UP62QEjKi6t1p/NeH+pOs6x44ds34Ez5W6N6zCeAK+ANIwC9npnePLy7w/JTpGdLthIaoQ+HUDVSioQiBNSf8ODbnnhhFomkan2YuCvmxIa5QinP7yoamW7BdNehzdMHG7VAxDYsiC9GfLZmAYBh3uWYw0Tbyait+QSNMkRVPJNwy+v3+GpVPzRxaRoqnk6cFmnjeUFkD2rGmoqkrTl+bjEYI8aSAIjjzkSwMVhS1X3WLJ1W7hYnyGRBMKAgW/LEj/MPtmpJS0XvwEASlxC4E0TXTTtNJbps+w7NTsiUV4hIpuSgzTxKuq+GVQj61Tb0YIQePnHiVF1fAZOpLgKEu+oWMCP0+aadnpivWz8BFAQcGNho8AAgUNjb9ftADTNBn56Wz86AgUVFQC6KgIBIK3ez2GlJIh624lgBG6rxDAsEYpdAzWXPwoYQ4dOkSVKlXQNC3pde9U/rRp0ya6dOnCxo0b6dSpE8WxadMmMjMzefndmrRu5y42X7xk/9fPVcP2s3z5clq3bh31WZ06dahTp84p/76o9u/u3bupW7cuV111FStWrLDyjhgxgrS0NF5++WVWrFjBxIkTrU5bmP79+9O0aVOeeeaZ09atpJwTE+WmTp3K22+/zSeffGJ9oQCffvop+/fvp0GDBtY9wzCYNWsWixYtYseOHYXKOpUBw/MCI3cxKGlaDQ2Vnk5aUZS40kpoePR05E22TtIw8ecFgmld4teDQcDQCwKKETAwglmC6VDZur8gABUnu6ZpUUOzZ0Onk9Nh/QxdYoT1i9BJCEHAp1s6hQnfO1mncHknp8tT3XO5XHHb9UzpFLZVpJ0C+aW3mS/XX2Q6nMc0C+6b0rTS0jDxnfBbehTnT+G6J4Sw5DzZn4qSPVKnSNnPRIxwuVxUr16dSM6Fuhf+fqGwncLylGV/SnSMMKSJIY1C6X9s2sF9od9ZXZqEd/APGAXTUgJGQewXQqCHTg7zBwrqZ2RaVVVk6N1nvl5wP0+PznPyfTMiLSPySEzypGGl88N6YEbJ5QvJHKlHZFoIQSDUqPVHTLuJTEfaxicL5M03ipHd0ItMR9VVAiH9TCstMfGH0oqi4Ee37stQ2kBiEBEjMCLuB9EpOkZExuJk171T+VPkv+MhcpThdAiX0bp161N2ZoqiuPZv9erV0TSNNm3aROVv3bo169evB6B27dr4/X6OHDkSNVqxf/9+unfvXlp1ToukTn8yTZMpU6bwxhtv8NFHH9G4ceOoz6+++mq+/fZbvv76a+vKyMhgzpw5rF69usTPK+l8O4dzHymlNU3GwT44drUnjl3ti1tTmTWqk2Nbm2Fnn9URBEz1tC+9FE3pWO1ft9tNly5dCm0zu2XLFho2bAhAZmYmLpeLNWvWWJ/v2bOH7777LmmdiqSOVEyePJkVK1awcuVK0tPTrTlglSpVIiUlhWrVqlnTWsK4XC5q165Ny5YtkyGywzlIxYoVky2CwxnAsas9cexqTwwp+XLLPsYOTrYkDonGrj4rEzRSUZotZWO1fwHmzJnD2LFj6dWrF3369OH999/nnXfeYe3atVbeSZMmMWvWLKpVq0bVqlWZPXs27dq1s3aDOtsktVPx9NNPA9C7d++o+0uXLuV3v/tdwp9X0qExh3MfIUTUbmEO9sCxqz1x7GpfDGnyyfe/Or+zNsPOPmuYAsM8/fpamjLiaf+OGjWKZ555hgcffJBp06bRsmVL/v73v9OzZ08r/8KFC9E0jTFjxpCXl0ffvn1ZtmxZ1PSys0lSOxWlWSNe1DqKeLHj8F15R8rgrk+NGzd2fsxshGNXe+LY1b64NZVbLsuMWoDrUPaxs8+aKMgErAIwSzFSEW/799prr+Xaa68t9nOv18vixYtZvHhxiWU4E9irhsQgvLjHwT4oikL16tUd29oMx672xLGrfTGkZPWmXxzb2gw7+2yit5R1OEd2fzpblFWn6C+uKPazNfL1Mv+8eDlX5XJIPIqiUKlSpaTKoNWqmbCylCQNRSeCRPrdAHVMXGXF88x45XLiRuIxi/gt1U3490/7rd9ZGWeVN1zxPLAEwsVAORTfFqKKETuPqcUnmEyNXZhayR8zD8AvR6vEzPPrsdix89CJtLiedy7E4jOFTND0J5mAMuxCufomDCOOKOFQpnCnuLj1jZsc29oMwzDYunWrY1ebEfbX8OnfDvbB41J5IOtCx2dthp1jsSS4yPr0L4cw5apTYbf5gA6g+3XeevQDx7Y2QwhB7dq1HbvajLC/6n49dmaHMkVAl6z4ZIvjszbDzrFYIjAScCViXYZdcKY/OZRppGGy+bOfHdvaDEVRSE9PT7YYDgkm7K8O9kOaJt/tPOTEYpth51ismyoB8/SbwbpZdqe3Jppy1b2y4/BdeceT6uaeNTMc29oMwzD48ccfHbvajLC/elLjm9fuUHbwulQWXnuR47M2w86xWJoKRgIuaTod6TDlqlNhx+G78k7AF2DZrNcd29oMIQQNGjRw7Gozwv4a8AWSLYpDgvHrkiff+9bxWZth51iciKlP4cshiDP9yaFMIw2TX77d5djWZiiKQmpqarLFcEgwYX91sB/SNPl57zEnFtsMO8diEyUhOzeV5pwKu1Kuuld2HL4r73jS3Dz42W2ObW2GYRj88MMPjl1tRthfPWnO9Ce74XWrPPP7Po7P2gw7x2JnpCLxlKuRCjsO35WURO3ZfqpySlrW6eDPC/D41c/zwpY/AYndu97Z3z55CCFo0qRJcn02zhNPUWLLaOr2n+4Tj7+E/dWfZ//vw84osrBv+H069736BSs7dQAg3vWv0h37La+veny+qNeOfdaDUtMXV1n+GrF3KKvS9EhcZcVznoUMJC7WmXHM8fdq8e3Adk7E4jNEeE1FIspxCFKuOhXOsKz9MKXJ3p8POLa1GYqi4PV6ky2GQ4IJ+6uD/TBN2H0o14nFNsPOsTh4xkQCDr9zpj9Z2K/reQrsOHxX3vGkuVn47T2ObW2GYRh89913jl1tRthfnelP9sPrVlk2Y4DjszbDzrE4uKXs6V/OlrIFlKtOhR2H78o7/hN+5vVf4NjWZgghaNmypWNXmxH2V/+J2NNUHMoWPr/BjOfWOT5rM+wci4PTn8RpX870pwLK1fQnB/thmpCfE98cWYeyhR1/xMo7YX+Nd7mKQ9nBBPKck9JtiV1jcaIWWTsLtQsoV9+ElDLZIjgkGE+am/kb5jq2tRlSSrKzsx272oywvzrTn+yH163yzOS+js/aDDvHYjN0cN3pXvEsjC8vlKtOhV172+UZX66f2y580LGtzRBC0Lp1a8euNiPsr75cZ/qT3cj3G/z+T/90fNZm2DkWGygJ2lLW6VSEcaY/2YTyuv2pooC3gifZYjicAaSUyf0hM+J7M6eoziK9eAn7q/+E35kCdQ5y4bjHiv1sw4pZVtpUCzeiFMCb4rL+LeNsXRxpE9vPhD++Rps4GvuhqhbfgmPPvthlHT5UIb6yKsWeoptSIb5pvK445O+W8UvMPC2r7I/reXAOxOIzhIlI0OF39vtuSku5+ibsOHxX3nGnupm3ZqZjW5shpWTz5s2OXW1G2F/dqc70J7vhcassuraX47M2w86xODhSkZjLIUi56lSozhtF2+HL9TOj/b2ObW2Gqqq0bdvWsavNCPurM/3JfuT7DSY8scbxWZth51ismyJBW8qWq6b0KUnqN/Hggw/SpUsX0tPTqVmzJiNHjmTz5s3W54FAgFtvvZV27dqRlpZGRkYG11xzDbt37y7V80xnvN12KEKhdtMajm1thmma5OfnO3a1GWF/VYTzZs9uKArUrZrm+KzNsHMsNs3g9KfTvcxSdCpitX9P5sYbb0RRFBYtWhR13+fzMXXqVKpXr05aWhojRoxg165dJZYnUSS1U7Fu3TomT57M559/zpo1a9B1nQEDBpCbmwvAiRMn2LRpE3fddRebNm3ijTfeYMuWLYwYMaJUz7Pj8F15x53i4uYXJzm2tRlSSrZt2+bY1WaE/dUdMffewR54XCp3jenq+KzNsHMsNhJ0ToVRit2fYrV/I3nrrbf44osvyMjIKPTZ9OnTefPNN3nllVdYv349OTk5DBs2LGmHFSZ1ofb7778f9e+lS5dSs2ZNNm7cSK9evahUqRJr1qyJyrN48WK6du3Kzp07adCgQYmeZ8fhu/KOL9fP3O7zbb0YvTyiqipt2rRJthgOCSbsrw72I99v8PtnPmbjsx2TLYpDArFzLDZRkAlYD2GWooxY7d8wv/76K1OmTGH16tUMHTo06m+OHTvG888/z4svvki/fv0AWL58OfXr1+fDDz9k4MCBpdDm9DinJoIdO3YMgKpVq54yj6IoVK5cucjP9+zZw6ZNm6Ku7OxsAKvnJqW0et3xpsNDf4ZhxJ02TbNQGog7HX6+qgnrzZ7qUnF5NSsdzqO5VTR3sNPk8miorlDaq0XpoWpBk7tTXIjIdGg3D8MwrLQnzW1NU/CkuVGUaBkVBWu/eUUoVlqoipVHqIolu4jUQxMF+p2kk8ujWTpFyR6hUzjtTXPT+PwG1vdVlE6eVHeUfmE7hXUKp8M6WfqdpJMntbB+pmlG6VSUnaSUUToVZafI+uZOcRVpp2TUvdL4SiL8Sdd1cnNzi9TjbOmkhWyjuVVUl4iwmShsM2/x/hRV307yp8g8p/KncN2LrG8n+1M4LaUs1p/CdS9S9kh/iqx7kXY6HX8Kp03TJKWih4bt66G5RZQ/hWWPJ+5Fyn46/uROcRXIHqHTyXaygz+VJEYoSvDMCQChKHhC6cg8qlBwh77HcFooCi0yKlt5NKHgVkN+IwSu0O5BblVE+1mo0niEGpVWlYjfpFDaq6pWOkVTrQZMpH4poZeHSkRaUBDXBQpeoRWklWAeFSVKLk9I9kg9NKHgFgX+4VJCOikCLZT2RKQj7eRRVNRQ49OrqAgK+1CKolnN08h0pJ28ihbST7HSIiJtmiYuQjECgVZEWkqJFnqnrKKioVpplWgfMk2TnJycYttPyWwbFecT8ZKYUYrgBZCdnV2o/blnz564ZCmq/Sul5Oqrr2bOnDmcd955hf5m48aNBAIBBgwYYN3LyMigbdu2fPbZZ6d83ubNm5k3bx59+/aladOm1KlTh/bt2zNhwgRWrFiBz1e6Q4XPmU6FaZrMnDmTnj170rZt2yLz5Ofnc9tttzFu3DgqVqxYZJ4lS5aQmZkZdWVlZQFw4MABAPbt28e+ffuAYCckfH/Xrl0cPnwYgJ07d3L06FEAtm/fzvHjxwHYtm0bOTk5AGzdupW8vDwgaKCwEbKzs9F1PerQGF3Xrc6Nz+ez5s7l5eWxdetWAHJycti2bRsAx48fZ/v27QB0HtGBiQvGANDzyi6Mu38kAP2u62lV2KHT+jJ0Wl8ARs8dTL/regIw7v6RUTp1HtEBgJuevZp2fVoBMOOl62jerYmlU/3z6gIwb81MajWuDsD8DXOpWDPd0smT5qZizXTmb5gLQK3G1Zm3ZiYA9c+ra+nUvFsTZrx0HQDt+rTipmevtnTauXNnkTqNnjvY0inSTpE69byyCwATF47lusVXIqVk+/btRep028rJlk6RdgrrFHkgV8Wa6ZadTtbptpWTLZ0i7RSpU1F2OnDgQJRORdkpsu5NXDCmSDslo+4dPXrUstPhw4etuZoHDhyw6t6Z8KctW7awc+dOS49k6DTod52RJ04w8g896HNFB+SJE4yd05sLh7SK0kmkpfL7Z7NoP7AdIi2V6cuvpUWvVoi0VEsn4Q7uelS7RR2EO1jfKtetRkqVCpZOlepUZf6GuSiai1rNajNvzUwUzUWDDg2tupeTk1OsP4Xr3uHDh4v1p3Ddi7RTlD9F1L1IO52OP4VjhM/n4+7VM/jdY1fQuGPDKH8K6xRP3Iu00+n404yXrrPqXqROJ8c9O/hTvDFC9yrUyKjAw1N7o3sV6jeqxP/d2BPdq0Tp1KlJDW4dlRm0TesMpg3tgFsTzBje0dJpVJvGTOzUClPAVR2acVWHZpgCJnZqFaXTsNoNEX6FO87rSN/qdRF+hfs7dObCyrUsndpmVMFIMVncpztNaqZjpJj8uV8valdPjdLJkyJ49eJ+eFIElSu4efXifki3SUblNMtOjZTqPNGoH/4TLtqodXiwQR/8J1x0ctePstOtfdqitP6NId1qMu3iliitf+Oy7nW5rldTy04T6rZEqJLfNziPsXWaIlTJzEYdGV6zQZSdjIDKPfV60jmlLkZAZX6D3rT11MIIqJad0jx+nmk8lGZpKaR5/Py16UjqpqpU80bbaUm9K2hW+RAXVJU8UfdSmlU+xEVVNR7JGGzVvfGeq6jqyqWDqz5jPaOp6sqls7sZoz1Drbp3qbcflV0nuNjTkcHeXlR2naC/pxv9Pd2o7Dph1T0pJb/88guHDh06o3UvEf4U9pt4MUnQ4Xeh7l9WVlah9ueSJUtiy1FM+/ehhx5C0zSmTZtW5N/t3bsXt9tNlSpVou7XqlWLvXv3Fvk3X331Ff3796dDhw588skndOnShenTp3PfffeRlZWFaZrccccdZGRk8NBDD5W4c6GY58jqm8mTJ7Nq1SrWr19PvXr1Cn0eCAS44oor2LlzJ2vXri22U7Fnz55CPcPs7GyysrL48ssv6dy5s9WbFULEnVYUBUUJvqEWQsSVhuj9naWUqKpq9bRjpU3TRAjBIPdYVJeKPy+A6lIRqkIgX0d1qbyXtwIhBIO9VwKg+w1cHg0pTYyAgcur8W7OS5YeQ7xXYegSd4oLPWAgw2m/zurAaxiGwRDvVUjDxJPmxp8XwJSh9Ak/q/XXkFIyyHUlihLcItKX60cRwbeNvlw/QlV4L/9lVFVloGsMmlvDnxdAaAItrIcmeC//5aB+niujdBJCIeDT0dwqq06sKJA9ZZylkzSC+rlTXBgBg/f9ryKlZLD3qiidpGHiSXUT8AUs/YQQDFDHWDqZZvCNpC/Xj6LAP/yvoKoqA7QxUTq5PC58J6L1M02TYRXGWzqdbKf3TryMlJJhaeMtnU620/u+V6z6NlAba+l0sp3eD7xy1uteaXzlXPeneGUcXm0Sut9A82iYYX8K2ey9o0ut/EOq3YDbq6EHJNKIrnurDi4J1vGUrCL9yTThvdwXgnm844v1J5fHxTtHl0XVt5P9SXWprMpdEaxvFcYX6U8A/8iPrm9F+ZOhS1brr6IoCgPUMUX6kzTMqDo5LH18kf7kTnXz7m8vYZomwytmFelPmltjVe4Ky04DtbFF+lMgX2e1/qole1H+ZAQMK48QgqFp44r0J3eKi3dzXgrKXmG8pdPJdnrn+PJy408XXP94cHTJpZLvD44QuDVBfsDg38/PsHTq+vtFqELBr0tUoVhpTVXY8KdpCCFoe+8iBAp+Ka1RioCUuIXg67umWrK3XvQ4umniUVUMKa20LiVbZs7AMAxa/ukJpGmSoqr4pCxIGwY/3TzT0qnxksdIUVXyjOAmn95QWgBbr5+Bqqo0eeEhPEKQJw0EwZGHfGmgorA1a44lV8eVD+IzDTRFBPWw0vDfUXcGZX/1YQKmxK0IJKCbEreiIjHZctVtlp2avzYfj6KimxIDE69Q8UsZzHfFrQgh6Pr+HXgVDZ9pYGKG0jom8PmA+yw7XfP5NPwEUFBwoUWlX7xwcbBR+O14dAIoCFRUdAIIBAKVBzoEY8TsbyZgoCNQUVAw0K1RCgODhR1eKnOxfNOmTXTp0oWNGzfSqVMnimPTpk1kZmYy4ZVLqN26SrH54mVv9hH+euVHLF++nNatW0d9VqdOHerUqXPKvy+q/btx40aGDh3Kpk2brLUUjRo1Yvr06UyfPh2AFStWMHHixEIN//79+9O0aVOeeeaZQs9q2LAhc+bMYdy4caecFbRhwwYWLlxIx44duf3222N+B2HOicPvpk6dyttvv80nn3xSbIdizJgxbN++nY8++qjYDgWc2oBKaNg08hCXkqYj12WUNq0oSlzpsLyGLjH0oPMYAQMjgJUOy6b7CxblBHx6QTpft/IIIaxy/HkBK09kWlVVpBHsZ0Zu+xhOR8pomgX3TWlaaWmYVh5pmFb5Upf4w3roBQHlZJ3Cmuh+I1r2gGHphJVHp3m3JgVBpgj9fCcK9Ii0R1H6mWZBnpN1CpcTqZ+iKNazirOTEMKySXF2iqxjp7LN6aZLWvdOx1dOx5+EEOTk5FChQoWk6RS2lR7pT8XZLKJOFmezoupbZJ5T+VO47kXWt5P9KVz3hBCWj5zsT0XJHihG9uLqZEn9KTJ2BPIDtOzelK1fbIvyp3D58cS9SLlOx5/8eQHrWZE6naxHOE9Z9qeSpE0zuD4CQJom+YHwtKgCnQxpYkgzKi0UhVZ1q1q66tIEgnkCEdNS/BGNSSEEeui9pi9iYWlkWlVVZChPXsT9vJPynHzfjEjLiDwSkzxpWOn8UNrAjJLLZ4b83yyQPTIthCAQ+rc/4r7fLNrPfBH3w888WfZ8Uy8yHZnHTyCkn1lkWlEUdOu+RCc0PSj0X1guAz10v0AWg8Kym6ZJbm4uFSpUKKRTsttGJ8tS0gP6zNBIw+lihspo3br1KTszRVFc+/fTTz9l//79UWuHDcNg1qxZLFq0iB07dlC7dm38fj9HjhyJGq3Yv38/3bt3L/J5W7duxe2OfU7QhRdeyIUXXojfX7Ltv5M6/ck0TaZMmcIbb7zBRx99ROPGjQvlCXcotm7dyocffki1atVK/bySzrdzOPfR3BojZw9wbGszpJTs3bvXsavNCPur5j4n3mc5JBCXJhh3UQvHZ22GnWNxMs+piNX+vfrqq/n222/5+uuvrSsjI4M5c+awevVqADIzM3G5XFEbGu3Zs4fvvvuu2E5FrA5FeGpbvPlPJqmRffLkyaxYsYKVK1eSnp5uzQGrVKkSKSkp6LrO5ZdfzqZNm3j33XcxDMPKU7Vq1RIr6+z+ZD/8eQEeGv00/eQlyRbFIYGoqkrz5s2TLYZDggn7q4P98AUMbn9pAwMvviDZojgkEDvHYpmgkYrSlBGr/VutWrVCL9FdLhe1a9emZcuWVt5JkyYxa9YsqlWrRtWqVZk9ezbt2rWzdoM6FQ899BCNGjVi7NixAIwZM4a///3v1K5dm/fee48OHTqUWK+kjlQ8/fTTHDt2jN69e1vTlurUqcOrrwbnw+7atYu3336bXbt20bFjx6g8sVa2F8U5snzEIYEITdChfxvHtjbDNE2OHTvm2NVmhP01vAOTg31QhUKXZjUdn7UZdo7FiTj4LnyVlFjt33hZuHAhI0eOZMyYMfTo0YPU1FTeeeeduF6iL1myhPr16wOwZs0a1qxZwz/+8Q8GDx7MnDlzSqwTJHmkIlYlbdSoUUIrsh2doryjuVR6X3NBaEtO55Reu2CaJgcPHiQ9Pd2xq40I+2v2+q3WehAHe6AKhUHnN3Risc2wcyyWKBgJOKeiNGddlKY9umPHjkL3vF4vixcvZvHixSUub8+ePVan4t1332XMmDEMGDCARo0a0a1btxKXB+fQlrJng5Iu4nE49/HnBXj86r84trUZQgiaNm3q2NVmhP01csG0gz3w65L7Xv/S8VmbYedYLKHILWJLfCVbkVJSpUoV/ve//wHBw/jCU6YizwUpKeVqtZwdFxqVd1RN0HlEh6jt6RzKPlJKjh49SuXKlZNn1wS+lTPjCNCKK3HhWDlH14+F/fU/b39j7ezkcO5geOKr84an8D1VKFzUIsOKxdqJ+J7pORL7mdITZ13xxvYz9WB8azGNrekx84hfvXGVpefG9u0D+yrFVVa6Kz+ufLGooMZ3/sA5EYvPEGYppy4VVU5ZZPTo0YwbN47mzZtz6NAhBg8OnnXy9ddf06xZs1KVWTa/CQeHEKpLpWP/NskWw+EMED5UycE+hP01fPK1g33QhELXJjWTLYbDGcCusViiJOwqiyxcuJApU6bQpk0b1qxZY20bvGfPHv7whz+UqsxyNVJht162Q3A6xZKbXuLyG0cmWxSHBCKEoFGjRskWwyHBhP3VwX74dMnD733FuMG9ky2KQwKxcyzWpSAgT/8Fhy7LVtvy9ttvZ+TIkXTt2pXZs2cX+jx8uF5pKFvfxGniTH+yH6pL5eKrL3BsazOklBw8eNCxq80I+6szUmE/NKEwuF0Dx2dthp1jcULWUyRoW9qzyZ49exg2bBh16tThhhtu4L333it0KndpKVedCgf7IVSFRh0Kn8LuUPY5cSLOSdkOZYawvwq1bP0IO8RGCIVmtSonWwyHM4BdY7GZoKlPZhmb/rR06VL27dvHa6+9RuXKlZk5cybVq1dn9OjRLFu2jIMHD5a67HLVqXCmP9mPQL7OX2f/zbGtzRBC0KBBA8euNiPsr4F8PdmiOCQYvy5Z/OG3js/aDDvHYkmCRirKWKcCQFEULrroIh5++GF+/PFH/v3vf3PBBRfw3HPPUbduXXr16sWjjz7Kr7/+WqJy7VdLToEdh+/KO6pLZeBNFzu2tRlSSvbt2+fY1WaE/dWZ/mQ/NKEwOrOJ47M2w86x2DSVhBx8Z5ax6U9F0bp1a2655Rb+9a9/8b///Y8JEybw6aef8vLLL5eonHK1UNvBfgihULlWxWSL4XAG0PXkvs1O5LasiuaKmccMxKevSE2NXVa8e4wrZ/e9UthfhVAo3S7oDmcSGWeLoMgdNFWFKhXi22I1krxasRur6on46qm2L/Z2sVqjnLjK0tNj19BA1QTGKBHfYWhbjtWImeebo3Vj5tl2onpcz4Pkx+IzRaLWQ5S1NRWxqFmzJpMmTWLSpEkl/tty1amw4/BdeSfg03l13jtcd/c1yRbFIYEIIahbN/YPo0PZIuyvDvYjYEj+/OkPTL18YLJFcUggdo7FidoOtixOfwLIz89n8eLFfPzxx+zfv7/QaNSmTZtKXGa56lTYcfiuvKO5VYZO6+scfmczwkPutWrVcuxqI8L+uuqJf6L7nbEKO+FSBWO6NHNisc2wcyw2pEjIdrBGGdtSNsy1117LmjVruPzyy+natStKAg58LXGnYseOHXz66afs2LGDEydOUKNGDc4//3wuvPBCvN6SD306ODg4ODg4ODg4nE1MMzFTl8z4Zq6dc6xatYr33nuPHj16JKzMuLtXK1as4IILLqBJkybMmTOHt956i08//ZQ///nPDBo0iFq1avGHP/yBX375JWHCJRq79bIdQPcbrHz0A8e2NkMIQZ06dRy72oywvzqjFPYjYEhe+nyL47M2w86xuDzv/gRQt25d0tPTE1pmXLWkU6dOLFiwgKysLHbs2MHevXvZuHEj69ev54cffuD48eOsXLkSKSWdO3fm9ddfT6iQicKZ/mQ/XB6NsfOGO7a1GVJKfv31V8euNiPsry5PuZp5Wy5wqYLrLmrj+KzNsHMsTsQZFYlal5EMHnvsMW699daEDgbEFdnvu+8+hg4dWuznHo+H3r1707t3b+6//362b9+eMAEdHE6FlCZH9x1PthgOZwBNcxqediPsr1KW0fkCDsVimiaHc/OTLYbDGcCusdhM0O5PZXVL2c6dO5Ofn0+TJk1ITU3F5YrepfDw4cMlLjOumnKqDsXJVK9enerV49+q7Gxix+G7ckUR218ausnqZz5l9p+mJEEghzOFEIJatWolW4gElhX7R8fUA3EVJRN5uq15dt8+GgGD1U+vO6vPdDgDFNEn1A2TNzZu454JQb8JxDmrwvVbbD8z4tjeFUCvEnvr08CutLjKqtb4SMw8Q8//Nq6yhDe2/BWrxOfXTSvGPu1YKLE77emaL67nnROx+AxR3reUveqqq/j111954IEHqFWrVnIWal9yySVcfPHF3HPPPVH3jxw5wmWXXcZHH3102kKdKew4fFfecXk1xt13qbPjiM2QUrJr1y7q1avn2NVGuLwa4+4fyYo733JO1bYZblXw+z5tnVhsM+wci8v7Qu3PPvuMDRs20KFDh4SVWeJOxdq1a/nvf//LV199xUsvvURaWrDX7/f7WbfOeQPlcHaRhsmOb3clWwyHM0BqHIe8OZQtpGGy45tdSKOM/go7FIs0TbbuO5psMRzOAHaNxeV9pKJVq1bk5eUltMxSdTs//PBD9u7dywUXXMCOHTsSKtCZxG69bIfgdIp1L37h2NZmCCGoXr26Y1ebEfTXzzECzu5PdkOXJv/4707HZ22GnWOxRGCYp3/J0jWlk878+fOZNWsWa9eu5dChQxw/fjzqKg2l+ibq1KnDunXraN++PV26dGHt2rWleviDDz5Ily5dSE9Pp2bNmowcOZLNmzdH5TFNk3nz5pGRkUFKSgq9e/fm+++/L9XznOlP9sOd4uLGp8c5trUZUkp27Njh2NVmBP11PO4UV+zMDmUKjya4dfD5js/aDDvH4oRsJ1vK0Y5Y7d9AIMCtt95Ku3btSEtLIyMjg2uuuYbdu3dHlePz+Zg6dSrVq1cnLS2NESNGsGtXfLM3Bg0axIYNG+jbty81a9akSpUqVKlShcqVK1OlSpUS6wSl6FSEF3J4PB5eeuklbr75ZgYNGsRTTz1V4oevW7eOyZMn8/nnn7NmzRp0XWfAgAHk5uZaeR5++GEWLFjAk08+yZdffknt2rXp378/v/32W4mf52A/jIDB1x/8kGwxHM4AFStWTLYIDgnGCBh8veYHZ6TChujS5Itt+5MthsMZwK6x2DSDOzed/lXyZ8dq/544cYJNmzZx1113sWnTJt544w22bNnCiBEjosqZPn06b775Jq+88grr168nJyeHYcOGYRixY+zHH3/Mxx9/zEcffRR1he+VhhKvqTBP+vbuvPNOWrduzYQJE0r88Pfffz/q30uXLqVmzZps3LiRXr16YZomixYt4o477mD06NEA/PWvf6VWrVqsWLGCG2+8sUTPs+PwXXnH0CVfvPm1Y1ubIYSgatWqyRbDIcEYuuSLN75KthgOZwBDmqzd/KsTi22GnWNxMreUjdX+rVSpEmvWrInKs3jxYrp27crOnTtp0KABx44d4/nnn+fFF1+kX79+ACxfvpz69evz4YcfMnDgwFPKcPHFF5dY7liU2Pu3b99eaMvYyy67jM8//5y//OUvpyXMsWPHAKwKvH37dvbu3cuAAQOsPB6Ph4svvpjPPvusyDL27NnDpk2boq7s7GwAdD2424iU0hrKizcd7kwZhhF32jTNQmkg7nT4+aomrOkCqkvF5dWsdDiP5lbR3CoQPGBKdYXSXi1KD1ULmtyd4kJEplXFkj2c9qS5UUREWomWUVGC9wEUoVhpoSpWHqEqluwiUg9NFOh3kk7hg7E0txote4RO4XRalVRufnGiZa+idPKkuqP0C9sprFM4HdbJ0u8knTyphfUzTTNKp6LsJKWM0qkoO0XWN3eKq0g7JaPulcZXEuFPgUCAn376CSll0nQK20lzqwV1rzibeTWEWrQ/RdW3k/wpMs+p/Clc9yLr28n+FE5LKYv1p7BOkbJH+lNk3Yu00+n4Uzhtmibp1dO4+cVr8VZwR/lTWPZ44l6k7KfjT+4UV4HsETqdbCc7+FNJYoSiQEroexSKgjeUjsyjKgqe0PeoimDaownuvbSL9TurCQV3qIPhEgJXKO0WItrPwt+1KtCUgrSqRPwmhQ4XSxFqQVrVrAZMpH4parCeKBFpQUFcF4pCihqhXyitKkqUXG4l5P+IqLRLKahjKqG/RUMQzlOQjrSTR6ioIdm9EXpEya5oKGH9FM06Ui3STm5CvoUSlXZR4ENuwnoruEJpFWGlpZRo1n0VNSodLbuUkp9++qnY9lMy20bF+US8JGaUQrE6FdnZ2YXan3v27IlLlpPbv8XlURSFypUrA7Bx40YCgUBUGzkjI4O2bdsW20beuXNnXPKE+fXXX0uUv8SdioYNGxb5JqJt27alGq0IY5omM2fOpGfPnrRt2xaAvXv3AhTaI7lWrVrWZyezZMkSMjMzo66srCwADh06BMC+ffvYt28fEOyEHDhwAIBdu3ZZh33s3LmTo0ePAsHOTXjRyrZt28jJyQFg69at1sr5zZs34/MF933Ozs5G13WklGRnZyOlRNd1q3Pj8/msuXN5eXls3boVgJycHLZt2wbA8ePHrUMEO4/owMQFYwDoeWUXxt0/EoB+1/W0KuzQaX0ZOq0vAKPnDqbfdT0BGHf/yCidOo8Ibh1207NX065PKwBmvHQdzbs1sXSqf15dAOatmUmtxsEO5PwNc6lYM93SyZPmpmLNdOZvmBu0SePqzFszE4D659W1dGrerQkzXroOgHZ9WnHTs1dbOoUr98k6jZ472NIp0k79b7wY4fUw/oFRXJR1AcLr4eqHRrN76wEURWH79u10HNEZrVo1Zrx6I636d0SrVo3b3plK4+5tCtkprJMnzc38DXMtncJ2Olmn21ZOtnSKtFOkTkXZ6cCBA1E6FWWnyLo3ccGYIu2UjLp39OhRy06HDx+25moeOHDAqntnwp+2bt1K5cqVURQlaToNvrEnqIJRMy+h7zVdQRVceedAelzWMUonM9/HDY+N4rwLGmDm+5j2zFiatauDme+zdJL5Pu5++yZq1ElH5vt44KPpVEh341JNS6dT+VO47uXk5BTrT+G6d/jw4WL9KVz3Iu007v6R9LyyS6G6F2mnyBhx28rJVoyIx5/COvl8Pm5/eyprX/icOi1qRflTWKd44l6knU7Hn2a8dJ1V9yJ1Ojnu2cGf4o0RUoPa1dJ4PKsXUoPGdSry0NjuSC1ap8wmNZg7LBNThZ4tM5g+oAMBTHYfy7XmfI9u1ZhJHVthChjXthnj2jbDFDCpY6sonQY2q4/hMbmlWwd6Na6D4TG5p0cmnevXsHTqWL0KeAyePv8iWlStAB6Dv2b2pl6llCidUivA6936k1oBqlXSeL1bf5RUnfpVvZadzmuu8VSP7ngb/EaXlik8dkEXvA1+46JW6VF2uqlibw4dSOdCpS0T0npy6EA6fbXzucJ7gWWnmocuZ9WmDjQ6Opb0/SNZtakDbX4bj7p7WJSdhCa5v+kFXFi1FkKTPNa8J+dXqorQpGWnnN+8LK4/nMr+6uT85uX5RqNx51fGyK0QZaepKVPJCXjxGnW4wXsjOQEvlWQDfue51qp7UyuORJoKbVyNuSF9CNJU6ORuwYQK/a26N8zTnxru3+jlOZ9Bnl7UcP9GX88F9PVcQA33b1bdUxQFRVE4cuTIGa17ifCnsN/EiyRBaypC3b+srKxC7c8lS5bElKOo9u/J5Ofnc9tttzFu3DhrOtrevXtxu92F1j+cqo3cpUsXrr/+ev79738XK8+xY8d47rnnaNu2LW+88UZM+SNRzJPnMxVDePpRLEoqQJjJkyezatUq1q9fT7169YDgHro9evRg9+7d1KlTx8p7/fXX87///a/Q8BEEg/DJPcPs7GyysrL48ssv6dy5s9WbFRFvTGKlw45lGAZCiLjSQNSe3VJKVFW1etqx0qZpIoRgkHssqkvFnxdAdakIVSGQr6O6VN7LW4EQgsHeKwHQ/QYuj4aUJkbAwOXVeDfnJUuPId6rMHSJO8WFHjCQ4bRfZ3XgNQzDYIj3KqRh4klz488LYMpQ+oSf1fprSCkZ5LoSRQF3qhtfrh9FBN82+nL9CFXhvfyXUVWVga4xaG4Nf14AoQm0sB6a4L38l4P6ea6M0kkIhYBPR3OrrDqxwpJ9WNVJQZ0i9HN7XRi6wXvHliGlZHjGFKQhcXuDOklp4klxEfDprNr3tGWbAeoYSyfTDL6R9OX6URT4h/8VVFVlgDYmSieXx4XvRLR+pmkyrMJ4S6eT7fTeiZeDsqeNt3Q62U7v+16x6ttAbSzuFBdGwChkp/cDr5z1ulcaXznX/SleGUfUm4zuN9DcKqZJVN1btedpK//gtAnB+hYwgnUv5E/SMHkvZ1mwvmlXFulPpgnv+4N1fKA2plh/cnlcvJvzUlR9O9mfVJfKqtwVwfpWYXyR/gTwj/zo+ubyakgj5E8RdW+1/iqKojBAHROlkyfVTcAXQBpmVJ0clj6+SH9yp7p597eg7MMrZhXpT5pbY1XuCstOA7WxRfpTIF9ntf6qJXtR/mQEDCuPEIKhaeOK9Cd3SvA7VRQl+H2FdDrZTu8cX15u/On82U+gKODVVPICBkJRcGuC/IDBNwumWzqdP3MRmqrg0yWqUNBEMK0Jhf88Mg0hBG0eWIRAwS+lNUoRkBK3EHx721RL9taLHkc3TTyqiiGlldalZMvMGRiGQfOXHkVikiJUfFKG0ho+qbPtmtssnZq8NJ8UoZEndRTAG0oLYMtVt6CqKm3fuhe3opJv6sE3+YrAZxqoKHwz4q4CuVY8hk8aaEpw3MFnSjRFBMu6+haklLRYupCAKXErAgnopsQjBIYJP0+aZdmp+esP4lFUdFNiYOIVKv6QHlsuvw0hBE1fftDSSQIpQiNf6pjA1rG3WnbK2nAzfgKh0QktKv1S9ycwTZMRn8zBT1A/DRU/OioCFcE7Fz+KlJJ5312FQQCBioKCgY4Ij2Sgc1/bv5W5WL5p0ya6dOnCxo0b6dSpE8WxadMmMjMzaf3ERNKa1S42X7zk/rSX7GlLWb58Oa1bt476rE6dOlHt16Ioqv0bSSAQ4IorrmDnzp2sXbvW6lSsWLGCiRMnWp22MP3796dp06Y888wzhco6fPgwDzzwAH/5y19wuVx07tyZjIwMvF4vR44c4YcffuD777+nc+fO3HnnnQwePLhE30XcayoqVaoU9e8VK1YwfPhw0tPjPDbzFEydOpW3336bTz75JOoLrV07aOy9e/dGGWX//v3FnvB4KgOG+0+RIy0lTauhodLTSSuKElc6vCje0CWGHnQeI2BghA7eNQIFDqr7CxblBHwFh0oF8nUrjxDCKsefV3B6b2RaVVVrD3lfrt+6H05HymiaBfdNaVppaZhWHmmYVvlSl/jDeugFAeVkncKa6H4jWvbQ4s5I/VBg5ou/wzCMkOwh/fILdPKdpN/JOkWmTbMgz8k6+U4U1k9RFEu/4uwkhLBkLs5OkXXsVLY53XRJ697p+Mrp+BPATz/9RJMmTZKmU9hWcdksP7bNiqpvkXlO5U/huhdZ3072p3DdE0JYB8ud7E9FyR55CF2k7MXVybAs8egXqZOiKJjS5NY3bmLh+D9H+VO4/HjiXqRcp+NP/ryA9axInU7WI5ynLPtTSdKmCXmhWCtNk/xAeFpUgU6GaWLooSkt0sSQJh5N5b4RXa3fWV2ahI/eDkRMS/FHNCaFEOih/L6IhaWRaVVVkaFy8mTB/TypR+U5+b4ZkZYReSQm+aZupX1msEwDM0ouX+hZuikJP0mPOIFeCEEg9G9/xH2fjM5j3TcLZM+X0foVpVNx+vkJhPQzi0wrioKfAv3CaQOJQUGMMEL5JQWyFGhaILthGGzbto0mTZoU0inZbaOTZSnpeh5pKhjm6a8BCq/LaN269Sk7M0VRXPs3TCAQYMyYMWzfvp2PPvooatF87dq18fv9HDlyJGq0Yv/+/XTv3r3I51WtWpVHH32U+++/n/fee49PP/2UHTt2kJeXR/Xq1Rk/fjwDBw4sdsQkFnF3KpYuXRr177/97W88/PDDVkUrDaZpMnXqVN58803Wrl1L48aNoz5v3LgxtWvXZs2aNZx//vlAwSF7Dz30UImf5ywgsx+6X+ftRR/Rd3jiFxw5JA8hBLVr13Z81mbofp23Hv0A3e+cpm03AoZk+RdbGNizW7JFcUggdo7FyTz8Llb7Fwo6FFu3buXjjz+mWrVqUZ9nZmbicrlYs2YNY8YEp4ru2bOH7777jocffviUz/d6vYwePTruWUjxUuLdnxLJ5MmTWbFiBStXriQ9Pd2aA1apUiVSUlJQFIXp06fzwAMP0Lx5c5o3b84DDzxAamoq48aNK/Hzwr1bB/sgDZPNX+xwbGszFEVJyCiow7mFNEw2f/ZzssVwOANI0+TbXw85sdhm2DkWm2bwSkQ5JSVW+1fXdS6//HI2bdrEu+++i2EYVp6qVavidrupVKkSkyZNYtasWVSrVo2qVasye/Zs2rVrZ+0GdbZJaqfi6aefBqB3795R95cuXcrvfvc7AG655Rby8vL4wx/+wJEjR+jWrRsffPBBqSp5PPv2OpQtPKlubnl1kjX9ycEeGIbB1q1bad68uWNXG+FJdXPbysnMv/RPUVOOHMo+XpfKI5d1d2KxzbB3LC7Yuel0yykpsdq/u3bt4u233wagY8eOUXk+/vhj6+8WLlyIpmmMGTOGvLw8+vbty7Jly5Jmq6R2KuJZI64oCvPmzWPevHmn/Tw7Dt+VdwK+AC/MXUmPzy9ItigOCUQIQYMGDRyftRkBX4Bls14n4AvEzuxQpvDrkic++pZeXTOTLYpDArFzLI7cDvZ0yyn535y6/duoUaO42sher5fFixezePHiEstwJoi7UxHuMYWRUvLPf/6T7777Lur+yaf9nUs4w7JnHyWBvWXF4yl0zwR2bi0YctdD2wYn5Hnxyq7YL9gmG0VRSE1NTa4M8dpfjW1/rWqV2HlOmi97OqgVKiSsrEQiDZNfvt2VbDEcisEf58HJRbWhDEy2HDhmxWLPkfjKMlJjny3gSolvDY55zB07T5wNQM2buHU/9avH/jIS88Y8SMO02M9r6o3v9PNzIRafKZK5psKuxN2pGDlyZKF7J59oHd4y7FzlXJbNoXR4Ut3c9fdrnSF3m2EYBps3b6Zly5aOXW2EJ83NvDUzmdd/QdQuSw5lnxSXypNjezmx2GbYORYnc03FuUBubi5paWkJLTPuV6yRpysWd53rjXY7Dt+Vd/z5ARbf9JpjW5shhKBJkyaOXW2GPy/A41c/H7W1q4M9yNcN7nn3347P2gxbx2IzMadqU0Y7FbVq1eLaa69l/fr1CSvThrWkeJzpT/bDlCb7dhx2bGszFEXB6/U6drUZpjTZ+/MBTFlGf4UdisU0YdfRXMdnbYadY3GiT9Qua7z88sscO3aMvn370qJFC+bPn8/u3btPq8y4OhUbNmyIu8Dc3Fy+//77Ugt0JjnXR1IcSo4n1c2jn9zs2NZmGIbBd99959jVZnjS3Cz89h48abHnvjuULVJcKi9f29/xWZth51iciFGKRC32TgbDhw/n73//O7t37+amm27i5ZdfpmHDhgwbNow33ngDXS/5uqK4OhXXXHMN/fv357XXXiMnJ6fIPD/88AO33347zZo1Y9OmTSUW5Gxgy+G7co4/z899lz3v2NZmCCFo2bKlY1eb4T/hZ17/Bfid7WRtR37AYPIrnzg+azNsHYvNBF5lmGrVqjFjxgy++eYbFixYwIcffsjll19ORkYGd999NydOnIi7rLgWav/www8sWbKEu+++m/Hjx9OiRQsyMjLwer0cOXKEH3/8kdzcXEaPHs2aNWtKfby3g0NJMU3IdxZ82hJb/oiVc0wT8nN8ZXZho0PxmEBewDkp3Y7YNRYnc0vZc4m9e/fywgsvsHTpUnbu3Mnll1/OpEmT2L17N/Pnz+fzzz/ngw8+iKusuDoVLpeLKVOmMGXKFDZt2sSnn37Kjh07yMvLo0OHDsyYMYM+ffpQtWrV01LsTCNl7K3rHBKLWrNG7Dxxbn+puApXV0+qi/vfuQEpZcJ3phBxyqVWqZSwZyqaK2FllWWklGRnZ9O6devk7TgS564YIi32dotG44zYBcWxNS3Et/WseY7GOk8FL/M/u5Xbuj901nZ/cnwqfgIV4+vtafmF73ldKs9de0nJY7EW+5mBo4W3Ey+S1NhTdBpVOxxXUW537A5Si2bxzT/fcyzOvXoTRMCMHUvWHW15ys+vDf3/nIjFZ4oE7f5UVkcq3njjDZYuXcrq1atp06YNkydPJisri8qVK1t5OnbsyPnnnx93mSU+/K5Tp0506tSppH92TmDX3nZ5xnciwJ3Dn+XNbY8nWxSHBCKEoHXr1o7P2gxfrv+sdigczh75AYPrn/uIzxadl2xRHBKInWOxmaATtc0yulB74sSJXHnllfzrX/+iS5cuReZp0qQJd9xxR9xlJvVEbQeH00VRwJvqLPq0I1JKW/6QlWcUBbxpHvwn/M4UKJuhAClup0lhR2wbi02KPsmxNOWUQfbs2RPzYMOUlBTuueeeuMu0YS0pHmf6k/1wp7i489XfOba1GVJKNm/e7NjVZrhT3cxbMx238yLAdnhcKk9M6OX4rM2wcyw2ZeKuskh6ejr79xc+Wf3QoUOlnupWrjoVtpsP6IDvRIA5ff/k2NZmqKpK27ZtHbvaDF+unxkd7nOmP9mQ/IDB1U+tcXzWZtg5Fpf3LWXNYoaLfT4fbnfpXvyUq7HK4r5Ah7KLIhRq1q+CaZq2PJynvGKaJj6fD4/H49jVRihCoVbj6uzbftA5AM9mKApkVE5zYrHNsH0sLodh6IknngCCBxv++c9/pkLEpjSGYfDJJ5/QqlWrUpVd4k7F9u3bady4cakelmzsOHxX3nF7NaYsvuyM7P7kkDyklGzbto2WLVs6drUR7hQXN78wkXkDFjmjFTbDo6ncc1lXJxbbDDvH4vK6pezChQuBYIfxmWeeibKr2+2mUaNGPPPMM6Uqu8SdimbNmtGrVy8mTZrE5ZdfjtfrLdWDk4HdHMIhOP3prhHP8Y99TydbFIcEoqoqbdq0SbYYDgnGl+tnbo+Hky2GwxkgP2Bww58/5qunOiZbFIcEYutYnKiD68rYaMf27dsB6NOnD2+88QZVqlRJWNkl7lR88803/OUvf2HWrFlMmTKFsWPHMmnSJLp27Zowoc4UzvSns4++Z2/sTFqc1bCIfEIo1G9Z0xpyV7q1L6GExWMcOxZXPqVyHOdUKPEtX1JE2XrjcaYwTZO8vDxSUlKSN+TujXNOaUrsFyvqzjj8IM6XHvqhQzHzKK74zmaI54wYkZISV1nCE/t7EKpC/fMy+N/3u5HG2YnHph44K8+xBXGaRBbhGkJRaFKjohWLzTjf4WmHEzcLW02JfbbEvpz0uMqqX+VozDwHcuM7y0hTEzdLIiPlaMw8tTy/xczTqeIvcT3vnIjFZwwldCWinLLHxx9/nPAyS+zNbdu2ZcGCBTz88MO88847LFu2jJ49e9K8eXMmTZrE1VdfTY0asQ88SwbO9Cf74fJoZN3RzxlytxlSSnbu3Enz5s0du9oIl8fF7x69nPkjn8Z3wpn+ZCfcmuDmvu2dWGwzbB2Ly+FIxcyZM7nvvvtIS0tj5syZp8y7YMGCEpdf6lcEmqYxatQohgwZwlNPPcXcuXOZPXs2c+fOZezYsTz00EPUqVOntMWfEWznEA748gL8Mesl/vHrBckWxSGBqKpa6oViDucuvhN+7h3gHFRpR/IDBlNf/pRvunVOtigOCcTWsdhUQCbinIqyM1Lx1VdfEQgErHRxlHZUqtSdiv/85z/85S9/4ZVXXiEtLY3Zs2czadIkdu/ezd13382ll17Kv//979IWf0Zwpj/ZDyEUmp1f19lxxGaYpklOTg4VKlRw7GojhKrQvFtjtn6x/axNf3I4OwhF4byMqk4sthm2jsUmiTmEswyFssgpT2di+lOJz6lYsGAB7dq1o3v37uzevZsXXniBX375hfvvv5/GjRvTo0cPlixZwqZNm2KW9cknnzB8+HAyMjJQFIW33nor6vOcnBymTJlCvXr1SElJoXXr1jz9dOkX5DrTn+yH5lYZfuOFjm1thpSSvXv3Ona1GZpbY+TsAWjOycu2w6UKsi5o4fiszbB1LDYTeJWQBx98kC5dupCenk7NmjUZOXIkmzdvjhbPNJk3bx4ZGRmkpKTQu3dvvv/++6g8Pp+PqVOnUr16ddLS0hgxYgS7du2KS4Zjx45x+PDhQvcPHz7M8ePHS64UpehUPP3004wbN46dO3fy1ltvMWzYsELHtzdo0IDnn38+Zlm5ubl06NCBJ598ssjPZ8yYwfvvv8/y5cvJzs5mxowZTJ06lZUrV5ZUbMCZ/mRH/Pk6j93wumNbm6Gqqj3n8JZz/HkBHhr9DP48Z/G03fDpBrf+fYPjszbD1rHYJDh16bSvkj963bp1TJ48mc8//5w1a9ag6zoDBgwgNzfXyvPwww+zYMECnnzySb788ktq165N//79+e23goX406dP58033+SVV15h/fr15OTkMGzYMAzDiCnDlVdeySuvvFLo/muvvcaVV15ZcqUoxfSnrVu3xszjdruZMGFCzHyDBw9m8ODBxX6+YcMGJkyYQO/evQG44YYbWLJkCf/5z3+49NJL45Y5jDP9yX4IVdC2eyNnyN1mmKbJ8ePHqVixomNXGyE0Qbs+Lfnvx5uRug3ffJZjVEUhs1ENJxbbDFvHYhOUJE1/ev/996P+vXTpUmrWrMnGjRvp1asXpmmyaNEi7rjjDkaPHg3AX//6V2rVqsWKFSu48cYbOXbsGM8//zwvvvgi/fr1A2D58uXUr1+fDz/8kIEDB55Shi+++KLIxdi9e/fmjjvuKLlSlKJTsXTpUipUqMAVV1wRdf/111/nxIkTcXUm4qVnz568/fbbXHvttWRkZLB27Vq2bNnC448Xv9Bvz5497NmzJ+pednY2gNVzCw/jCSHiTiuKgqIoGIaBECKudPhZkWlVVTFNM660aZoIIVA1gepS8ecFUF0qQlUI5OuoLtUqX/MEt5DU/QYuj4aUJkbAwOXVrDxSSlRNYOgSd4oLPWAgw2m/bn1HQlWQhoknzY0/L4ApQ+kTfktGCJ6g6k5148v1owgFd4oLX64foQa/A1VVcbdpgeZS8ft0hKqgaoKAz0BVRYHslSsghBK87wp+h7rfQIvQT0qJWjkdQ5e43GpQP13ireim1xUdre9L2fhdlE7SMPGkugn4AgX6hexRoFMw7cv1oyhYsqspKVE6uTwufCei9TNNE7Hn16BtirBT2O7u9FR0n47mVi07aR7NOlXY+k49HtxeDUOXQTt5NfSARBoyKXWvNL6SCH/SdZ0DBw6Qnp5eSI+zpZOo6EEPSDSXwDSJqnuRNiM1BZdHxdAl0jBxe7Sgb0mzwB4pXjwpLvw+PehPKS78+YFgueE8eXm4U9z4ToT8yVtQ31yh6UKmaZJatzb+/ABCFaguEYwFmkBoBd+LOz2FgE8v5E9hLNkVJSpeRNY9y05ud7Q/pbgIhNKRNnN7lCL9yZ3qtmRPr1mZ3hMuZOumXZgmln5ahH7hlz/Fxb3IaRiR/hSpR6Qt3SkujIBRZNwLN4DDMaKouBeWpyz7U7wxAoLTFzyqSp5hIBQFjxDkGUaUTkYlBU0o+AyJJhRURUECQzo0xDAMNE1DuIMyBQyJK/SMgJS4ImK/lMG/16WJRxUY0kQ3g2k95GeGYaD6VKRpkqKq+KQsSId+08M6GUfdpIRkVwBvWA8K4npOTgpuoZIvdQQKbkWQbxqoKFFyNa9wHJ0AAhWBgo6OQEWhIK7f3uJDTEVHMV2AGUq7AQncZ32nDaocwYWGgYHExI2GHkqHZa9SJRevouEzDUzMUFrHJPp3a8vRhuSbARQUPIpGvhlAoOBSVMuH/nu4MX6C+mmo+NFRQ5pYtpYepBJAMVUUBFIJIMygH0ol2s8OHDhAWloamqadU22jk2Up1RStBL5rDrc1I6lTp05cmxYdC21hX7VqVSB4lsTevXsZMGCAlcfj8XDxxRfz2WefceONN7Jx40YCgUBUnoyMDNq2bctnn30Ws1Ph8/nQ9cLbMAcCAfLy8mLKXBQlnv40f/58qlevXuh+zZo1eeCBB0olRHE88cQTtGnThnr16uF2uxk0aBBPPfUUPXv2LPZvlixZQmZmZtSVlZUFYM0d27dvH/v27QOCnZADBw4AsGvXLivPzp07OXr0KBA0bnh+2bZt28jJyQGCozbhL37z5s34fD4gWLF0PeiU2dnZSCnRdd2qcD6fz5o7l5eXZ43+5OTksG3bNgCOHz9uHVDSeUQHJi4YA0DPK7sw7v6RAPS7rqfVgRo67RKGTrsEgNG3DaLfpB4AjLvv0iidOo/oAMBNz15Nuz7BHR1mvHQdzbs1sXSqf15dAOatmUmtxkFbz98wl4o10y2dPGluKtZMZ/6GuQDUalydeWuC25PVP6+upVOzdrWZ+kCwwrftUo8b7+oLQObFjdm5cycAPYadx9hZvQG4ZExHRt7UHYDBv+sSZadLRrcFYOzU7vQY3DKo3809+fc/f0YIwfbt24vU6baVky2dIu0U1smT5mb+hrmWTmE71WxUjbvf/UNQpzZ1uOW16wBo3qVRlJ1uevbqU9rpwIEDjJ4dfIsw5A8XM+QPFwftNLsffX8X3LUqsu5NeGAEmYODhw3d8PjltO3VFEhO3Tt69Khlp8OHD1tzNQ8cOGDVvTPhT1u3bqVevXoIIZKm05DL2wXtdPX59B0WrFdXXd+FHn2bFtLpxrv60rZLPQCmPjCAZu1qF9LpzmVjqVkveKbJ/a9dTcWqqXhSXJZOFatX4I+rg/WtZoMq3PXGJADqt6rNnOXXWDpNffpyANr2bMyNj44EIHNAKyb83xBLp7G3Bv3skqsyGTm1FwCDr7uAwdddUMhOV945iB6jg3Ehsu5F2mnakitpltkAgDnLr6F+q6B+kXXvwfW3ULFG0J8eXH9L0J9qpPPg+lssO932t+tYfN1L1GpcnVteDerXvHNDbv7L1YXsdCp/Koh7fRk6Lajr6LmD6Xdd8Lchsu5NXDCm2LgXrnuRMeLkuGcHfypJjKhbIY2/9AnGqOaVKvLkRT0K6dSlbg3uvjgTgF4NM5jdowM+Q/LZrn3s3r0bgFHnNebazJDfdGzGVR2bAXBtZqsonYY0rB+0QWYHetfLAODebpl0q1XD0qlDlWBD66muPWmeXhGAZd17Uy81LUqnFFXl9Qv7k6KqVHV7eP3C/gDUS61g2al5SmWebBqsM+3TavBw46Cu3dLrRNlpiGtUMI/aiX7aUAA6qxfQS+tn2ammHvTFOnoW1fVgnrqB66hqXFLITjdWHEJbVyMAplYcSTNX3UJ2erzupWS4gudoPFv/ciqrKXgVLcpOj9QItmdqqZX4v2pB/2igVeeOqqMtO82oPAKAdu6G3FQpOBuks6cZEysG9T58+DBt5VUANDb70coM6trMHEIzMxhHwnVPCIHb7bb0OJfaRif7U9hv4iYhU58Ua/enrKysQu3PJUuWxBbDNJk5cyY9e/akbdtgO2fv3uDZRrVq1YrKW6tWLeuzvXv34na7Cx1eF5nnVHTp0oVnn3220P1nnnmGzMzMmH9fFIpZwjlBXq+XH3/8kUaNGkXd37FjB61bty5170ZRFN58801Gjhxp3Xv00Ud57rnnePTRR2nYsCGffPIJc+fO5c0337SGek6muJGKrKwsvvjiC7p27VrmRioGuccWO1LxXt4KhBAMThkHFD1S8e5vL1p6DPFeVexIxerAaxiGwRDvVcWOVKzWX0NKySDXlaccqXgv/2VUVWVo2zuKHal4+5v/QwjBsCYzih2pWPnTo5bsI9rdWXikIs1Fh+6NuPuJiQAM9l5V7EhFWD8hBAPUMcWOVPzD/wqqqjIo7ZpiRyrePbbUstOwCuOLHal478TLSCkZXnVSsSMV7x1fZtW3IVWuK3akYtXh55L+JihR6Vj+FAgEOH78OFWqVLFkOds6De9+X7EjFe99cY+Vf2jbO4sdqXj3v/chhGBIk1nFjlS8+9MjwTzVrj/lSMXK3U9jmiYjG0wrdqTi7f89iZSSkY1uLnak4p1dT0bVt+JGKt478mcURWFw1euLHalYdXCJZbPhla8pdqTinaMvYJomozNupF3vFmxa/T1CVaNGKt45+Lxlp4Ha2GJHKlbrrwIwUBtb7EhFOI8QgqFp44odqXg35yUURWFYhfHFjlS8c3x5mfeneGNEi4WPFztS8dOsGZZObe9ZVGikwjBNejeqy8KsUWiaRvs7FxU7UvHVvVMt2ds89kSxIxWbb5mOYRi0ePrxYkcqfp46y9Kp8ZLHih2p2Hr9DFRVpcVrfyx2pCL78rmWXNO+uq7YkYo/ZS5FSsnffu5R7EjFmOYbCmLE+hnFjlS822MBQggy/3FXsSMVXw6817LTkLVzix2peL/PQ8HfpHW3FDtSsar3I0gpWfzj8FOOVNzcsmB6zuHDh6lcufI5P1KxadMmunTpwsaNG+nUqRPFsWnTJjIzM6k992Y8DeoVmy9efDt3sffBx1m+fDmtW7eO+iyekYrJkyezatUq1q9fT716QXk+++wzevTowe7du6P+/vrrr+d///sf77//PitWrGDixIlWpy1M//79adq0Kc8888wpn/uvf/2Lfv360aVLF/r2DXY4//nPf/Lll1/ywQcfcNFFF8X9HYQp8fSnmjVr8u233xbqVHzzzTdUq1atxAIUR15eHrfffjtvvvkmQ4cG3wC0b9+er7/+mkcffbTYTsWpDBiuwOH/lyYduViptGlFUeJKh+cvhn/kgeCPY2iNoxEocFDdX7AoJ+ArGM4K5OtReofLiVwoGZlWVdXa6tGXW3A4VTgdKaNpFtw3pWmlpWFaeaQ08YfkkYaJtIarCwKKEZCEpTcCBcOXeoR+kbIHInQ1dEm7C+pbeWQR+kUeshVpj6L0M82CPCfrFC4nUj9FUaxnFWcnIQR66DuItJMeYafIOubP14tMJ6PunY6vnK4//fbbb1SpUiVpOumhuqhH1MnIuheZP+AruO/3FW0zX2SdPMnnAGs6EITqXkR9C+dXFAV/fiB0P9jZhOi6J4SwYsDJ/lS07EXXt6g6GUN2KN6fImOH7jfo0LclX32QHaWfP0K/WHEvUq7i4l5xsp+cDj8rMkacrEc4T1n3p3jTEsgLTxU2TSsdqZNhmhih3wldmugEOwLd6tUo+E2SBdviBCKmpQQiYr8QBdOcfEZBnsi0qganPhEh18npSD3C981IPaL0M8mXekHaDP0mYUbLRcjPMAhLU/BLFcxjKsFyTKWgXpmKPyqPpTcRfkbRMSLf1ItMR+cJhPQzrbTExBfKryiKVb7EtNJGhCZCBDsRQXkNzJBeUinsQ1JKcnJyrGk551Lb6GRZIv8dFzJ0nS6hMlq3bn3KzkxRTJ06lbfffptPPvnE6lAA1K4dHA3eu3dvVJt2//791uhF7dq18fv9HDlyJGq0Yv/+/XTv3j3ms3v06MGGDRt4+OGHee2110hJSaF9+/Y8//zzNG/evER6hCnx9Kcrr7ySadOm8fHHH2MYBoZh8NFHH3HzzTeXerV4UQQCAQKBQKFKoqpq6ebNQckrnMM5T8Bn8Of7PnJsazOEEDRq1Mixq83w5wd49ubXrY6Rg33wGZI/fvKV47M2w96xOFFTn0q+gN00TaZMmcIbb7zBRx99ROPGjaM+b9y4MbVr12bNmjXWPb/fz7p166wOQ2ZmJi6XKyrPnj17+O677+LqVAB07NiRFStW8P3331vnz5W2QwGlGKm4//77+eWXX+jbty+aFhoqk5JrrrmmxGsqcnJy+Omnn6x/b9++na+//pqqVavSoEEDLr74YubMmUNKSgoNGzZk3bp1vPDCC6U6Ojwsp4O9UDVBj8Eto4ZSHco+UkoOHz5M1apVHbvaCNWl0uPy8/nX377CCMTe8tCh7KAJhUHN6jux2GbYORYrCdr9qTRlTJ48mRUrVrBy5UrS09OtNRCVKlUiJSUFRVGYPn06DzzwAM2bN6d58+Y88MADpKamMm7cOCvvpEmTmDVrFtWqVaNq1arMnj2bdu3aFTub52QMw+Ctt94iOzsbRVFo06YNI0aMKPUWwiXuVLjdbl599VXuu+8+vvnmG1JSUmjXrh0NGzYs8cP/85//0KdPH+vfM2cGF/pOmDCBZcuW8corrzB37lzGjx/P4cOHadiwIX/84x/5/e9/X+JnOdgTIRQatCi8cYBD2efEiRPWkLuDPRBCoVG7DDa88TVOl8JeCEWhRbXKyRbD4Qxg21hcyoPriiynhIQPcg4fmRBm6dKl/O53vwPglltuIS8vjz/84Q8cOXKEbt268cEHH5Cenm7lX7hwIZqmMWbMGPLy8ujbty/Lli2Lq1Pw008/MXToUHbt2kXLli0xTZMtW7ZQv359Vq1aRdOmTUusV6mPNW3RogUtWrQo7Z8DwS/zVOvEa9euzdKlS0/rGZHYrZcdiZqWelafp2iuuPLJbb/EzKP/ujuusvStPxe+Byz7/RbG3zA0rjJKRJsmcWUTHm/MPIo7vu9LSYldVnlACEGDBg2SKkOgckqcGQtvyVcItzt2HjXO+CRiD7Wb/vimF8VT38x8X8w8gLU98qnwHf2NZTNejau8mCgJjOeJLKssUy8/rmyBtCLuIXno22+5fnBwwaevanwtLemOnc+sGIePQVyNO48nvrLW/nKKKSChjXHu/Kr487LGNitI1009FvN5dSqU7gTjM825EIvtSDx7JCmKwrx585g3b16xebxeL4sXL2bx4sUllmHatGk0adKEDRs2WJ3GQ4cOkZWVxbRp01i1alWJyyxxp8IwDJYtW8Y///lP9u/fX2hK0UcffVRiIc4WzvQn+6G6VPpd19MZcrcZUkoOHDhAjRo1HLvaiLC/fvjn9c70J5uhCYXLmzZ2YrHNsHMsTub0p3OBdevW8fnnn0eNQlWrVo358+fTo0ePUpVZ4k7FzTffzLJlyxg6dCht27a13wmLDmUKIRQq16qYbDEczgBFHcrjULYJ+6sQijP9yWYIFKp5nZFWO2LbWBxxxsRpl1MG8Xg8/Pbbb4Xu5+Tk4I5ndL0IStypeOWVV3jttdcYMmRIqR6YTOzWy3YIbiP56rx3uO7ua5ItikMCEUJQt27dZIvhkGDC/upgP/xS8qf//sCMIac+xdehbGHrWGySmC1ly+hIxbBhw7jhhht4/vnn6dq1KwBffPEFv//97xkxYkSpyixxK9vtdtOsWbPYGc9BnOlP9kNzq1w6e4BjW5shpWTPnj2OXW1G2F/Dh9Y52AeXEExs3cLxWZth51gcnv6UiKss8sQTT9C0aVMuvPBCvF4vXq+XHj160KxZMx5//PFSlVnikYpZs2bx+OOP8+STTzpTnxwcHBwcHBwcHMoeSdz96VygcuXKrFy5kq1bt/Ljjz9imiZt2rQ5rYGDEncq1q9fz8cff8w//vEPzjvvPFyu6F1t3njjjVILc6Zxpj/ZD91vsPLRD5jy8PXJFsUhgQghok4RdbAHYX91sB8BKVmavYU7RpyBnfgckoatY3E571SECZ+DkQhK3KmoXLkyo0aNSsjDzzZ2HL4r77g8GqPnDnZ2HLEZ4SH3OnXqOHa1EWF/fePBfxDw2XTxZznFLQTXn9fKicU2w86xuDzu/hQ+Dy4eSnPQdIk7FYk8N8IhgcR5DkKiEBWK2Ki8CJSWsc96EOe3KbUcUpoc3Vewv7co5Y4FRZadEl9Z/ovaxs4U7+mUcX6v5QFNK/UxOglBzYuv0avXjGP3sXjOsqhTI67n4fHEzhPvmRdxfMdKWpzndZixX9qE/VXGONNCrVAhIc+LF9XxOwCM/PjiVNqewvc0YXK8VsE5F9ITX0tLBGJPoxbu+PYKq9vsQMw8Pl98cWXaeR+f4tO7AUhPi+9cj3jI3lU7rnzVvSdi5ulaNfb5UCUh2bH4zJGg3Z8oO0sBvvrqq7jylXZ5Q6lqiq7rrF27lp9//plx48aRnp7O7t27qVixIhXi+TFIEnbrZTuAETBY/fQ6Zv9pSrJFcUggQghq1aqVbDEcEkzYXx3shy5N3ti0jXuucX5n7YStY3E5nP708cen6iyfPiX2/l9++YV27dpx6aWXMnnyZA4cCL4ZePjhh5k9e3bCBUwkzvQn++Hyakx49HLHtjZDSsnOnTsdu9qMsL+6vHZ981l+cauCaZe0d3zWZtg6FktQEnAlZFvaJPLTTz+xevVq8vLygPhO+y6OEncqbr75Zjp37syRI0dISSkYFh81ahT//Oc/Sy2Ig0NpkIbJjm92JVsMhzNAampqskVwSDBhf5VGGXq15xAX0jTZuv9ossVwOAPYNRaX9y1lDx06RN++fWnRogVDhgxhz57gvMbrrruOWbNmlarMEncq1q9fz5133lnotL2GDRvy66+/lkqIs4Uz/cl+GAGDdS9+7tjWZgghqF69umNXmxH2VyPgnKdtN3Rp8o/vdjo+azNsH4vNBFxllBkzZuByudi5c2dUx3Hs2LG8//77pSqzxLVESolhFP5B2LVrF+np6aUS4mxhy+G7co47xcWNT493bGszpJTs2LHDsavNCPurO+XsbizhcObxaIJbB53v+KzNsHUsTkSHogx3LD744AMeeugh6tWrF3W/efPm/PJL6Rb7l7hT0b9/fxYtWmT9W1EUcnJyuOeeexgyZEiphHBwKC1GwODrNT8kWwyHM0DFinHsquRQpgj7qzNSYT90afLF9v3JFsPhDGDXWFzepz/l5uYWObXt4MGDeOLZZbAISrxabuHChfTp04c2bdqQn5/PuHHj2Lp1K9WrV+fll18ulRBni7I6fKdosd/qKSlxbvuYIJQ4t7D11Yq9G5ieHl/lFUXoaAJf/uNHy7bS74+rrLj47Ou4splq7K3XFFecruZLoPxlGCEEVatWTaoM0hvf9prSHTufXrtyzDymJz6fMqvG8QOfUTO+snJjb0+JPxBXWXHhTeU/H2wFbyqn/NYSuEW2WrlyzDwyL3Fbg5YH8qoXddfkH4d+5fEz8Dtr5MUXP4/neWPmGdwsO66yFnzTt9jP/tAy+P8aaTlxlfXz8Wox81zfcX1cZfWrFvslWquU3THzDEv/Lq7nnQux2OHM0KtXL1544QXuu+8+IDhIIKXkkUceoU+fPqUqs8SdioyMDL7++mtefvllNm3ahJSSSZMmMX78+KiF2+cithy+K+e4vS5ufHKMc+CSzZBSsn37dho3buzY1Ua4vRo3PHEFz057HX++c/idnfCognk9Mp1YbDNsHYvL4ZaykTzyyCP07t2b//znP/j9fm655Ra+//57Dh8+zL/+9a9SlVmqff1SUlK49tprufbaa0v10GRR2sM8HM5d9IDBuhX/YeiVA5ItikMCURSF6tWrOz5rM/SA5JOXN6IHnBc8dkOXJm//9AujL+qRbFEcEoidY3F5PFE7kjZt2vDtt9/y9NNPo6oqubm5jB49msmTJ1OnTp1SlVniTsULL7xwys+vueaaUglyNrCjU5R3pCH59qPNjm1thqIoVKpUKdliOCQYaUi+/XhLssVwOAMYpsmG3fudWGwzbB2LTRJzxkQZ7VQA1K5dm3vvvTdh5ZW4U3HzzTdH/TsQCHDixAncbjepqanndKeiqF2rHMo27hQXN//lagzDQFXjmwPvcO5jGAbbtm2jSZMmjl1thDvFxbTnxvHE9Svw5yVwrYZD0vGqKvMv7urEYpth51hc3kcqGjduTFZWFllZWbRs2TIhZZZ4gtyRI0eirpycHDZv3kzPnj1LvFD7k08+Yfjw4WRkZKAoCm+99VahPNnZ2YwYMYJKlSqRnp7OBRdcwM6dO0sqNlB2F2o7FI/u13l70UeObW2GEILatWs7drUZul/n7cVr0f3Oegq74ZeSZd9tcXzWZtg6FidxS9lY7d+cnBymTJlCvXr1SElJoXXr1jz99NNReXw+H1OnTqV69eqkpaUxYsQIdu2K/zDgqVOn8v7779O6dWsyMzNZtGiRdQBeaUlILWnevDnz588vNIoRi9zcXDp06MCTTz5Z5Oc///wzPXv2pFWrVqxdu5ZvvvmGu+66C6839g4PReEMy9oPaZhs/mKHY1uboSgK6enpjl1thjRMtnzxi3Oitg2RpsnX+w85PmszbB2LE7WdbCnCWaz274wZM3j//fdZvnw52dnZzJgxg6lTp7Jy5Uorz/Tp03nzzTd55ZVXWL9+PTk5OQwbNizuWTkzZ87kyy+/5Mcff2TYsGE8/fTTNGjQgAEDBsRc6lAcpVqoXRSqqrJ7d+xtzCIZPHgwgwcPLvbzO+64gyFDhvDwww9b95o0aVJqGcvq9CeRVngf4UJUiCNPnGiNGsbMk9cxdh4A7/+Oxszj+nZ7XGWJCmmF7nlSXMxeOs4actfqZsRVVjxozZvGlS/lv3G8GTjpBHqHU2MYBlu3bqV58+ZJG3LPrxafzaQ79rsZ7UhuzDxGWpzP88YO28rh43GVpaTG3rHPDMQ3qiA6tImZR5N+bls5mfmX/gnfieK3TzYOH4ldVo0i9zUtFXHFWAcLX53CU9dShMozHXpZsdh0xdfSMkXsfPUbHIyrrN92xz6A970t5xX72eJOcT3mjPDq9sxiP5sb4Vo3tlwXs6xL02K3xZpqsbd7h3MjFp9RkvR+I1b7d8OGDUyYMIHevXsDcMMNN7BkyRL+85//cOmll3Ls2DGef/55XnzxRfr16wfA8uXLqV+/Ph9++CEDBw6MW5YWLVpw7733cu+99/L5559z0003MXHixFItZyjxSMXbb78dda1cuZJnnnmGq6++mh49Erfrg5SSVatW0aJFCwYOHEjNmjXp1q1bkVOkItmzZw+bNm2KurKzo/elllJa28vGmzbNYM0zDCPutGmahdJA3Onw81VN4A41JFSXistTkA7n0VwCzRU0p+ZWUbVg2uVWo/QI33enuBCR6dBZC4ZhIEQw7UlxoUSmlWgZFQW8IVki00IoVh4hFEt2oSq4PKqlk6XfSTppoX3/tZNlD+nn8mhWGgVW/HE1QgTLC+vh9mpReoTTkXbypLkJv3zxhBp0ilLQ+VSEgid08q8QCh5vYf1M0yyQXRMReghLVyllgU4uFc1VoF9Yj8j65vJoBfbzaFG2Odt1rzS+kgh/Mk2T+vXrI4RImk6aGrKBJtBCNohMR+Z3u1WrjkWmo+pbihbhT5pV98J5FAWrjkWmhVDwhOqSaZq4Q2mhKrhD9UqNSEspC/zMJYqse4XqW6RvaQV1Mix7ZIzwpLqLrJMeb4FORelhmiZCU1g263UMXceT6i7QI+RnheJe6L7qUnGFyomUPR5/cnsLdHJ7NYRaEPcKbBOtn2WnVLct/KkkMUIAKSJU91BIEQV1L5xHVRQ8oftaKO2Tkod/+powmqLgDk2ZcSkClxL63oWI9rNQpfGoalRaVSJ+kwjdVzQr7VVcKBSuhymhhq8SkRYQoV+BTgIFbyitokTJ5VHC+gncUemCOqaFTlvRUFFDzSlXRDqqHlIguztCjyjZFc3SKZimSJudbI9C9c30htIq0vSE0pqVllIizaD/SVOLSLuQpitKdiFE1InL51LbqDifiJsET3/Kzs4u1P4s7XSinj178vbbb/Prr79imiYff/wxW7ZssToLGzduJBAIMGBAwc6XGRkZtG3bls8++6zEz/v3v//N9OnTGTVqFJs3b+byyy8vldwl7lSMHDky6ho9ejTz5s2jffv2/OUvfymVEEWxf/9+cnJymD9/PoMGDeKDDz5g1KhRjB49mnXriu+pL1myhMzMzKgrKysLgEOHDgGwb98+9u3bBwQ7IQcOHABg165dHD58GICdO3dy9OhRALZv387x48E3f9u2bSMnJ3jgzdatW8nLywNg8+bN+Hw+IFixdF1HSkl2djZSSnRdtzo3Pp+PzZs3A5CXl8fWrVuB4By6bdu2AXD8+HG2bw++wc8c1JoJfxwOQI/R7bnyjmAl6pvVxaqwg8d1ZPC4jgCMvLYzl4wKvo0ZO+XCKJ06j+gAwE3PXk27Pq0AmPHSdTTv1sTSqX7z4BvAO5deQc26wV0f7n91PBWrplo6eT0alSp6efSPwQpeq2YF7rvrEgAa1q9k6dSsXW2mPhDM07ZLfW68O3igUGavxtbamB6j2nPl7WGdOjNqeu+gTtd3j7LTJVcF3+SMvbUvPS5tB8DVdw+iRv0qKIrC9u3baXthIwCmLhpJs451AZj97OXUb1GjkJ3mb5hLxZrpeNLczN8wF0+am4o10y071axbkTufHQVA/WbVmP1E0AbN2teOstON84cGderXnAl39w/qNPw8xs4OHh5z4MABRt7UPajT77ow+Hddgna6qTuXjAnaLLLuTbirP5l9mwNw44NDLJ2SUfeOHj1q2enw4cPWfM0DBw5Yde9M+NOWLVtQVRVFUZKm06iQDa4c3JpBPYP+MXFkO3p3aVBIpyk3dqN921oAzL65By2aVyuk053PjKZmRvDguvtfGEvFKil4UjRLp4qVvDz02BAAatVO5977g3WpQcPKzA35Vk5ODtNvD9ardh0z+P2sngB07t6Q3/3hAkuncN27ZOz5jPxD8GXP4IldGDyxSyE7jZ3dmx7Dg/Fiwt0DyOzXopCdbl52Dc27NALglteuo36b4HaDkXXvgedGUbFyCh6vxgPPjcLj1ahYOYUHnhtl2enu92fwy7e7qNsqg9tWTgagebcmzHjpukJ26jyiAxMXjAGg55VdGHf/yEJ2GnzdBQy+Lqj3yKm9rBgR5U/3DyNzYGsAblgwmrYXBW057ZmxVt2bs/wa6reqDcBdb0yiZoMqAPxx9R9s4U8liRH1vBX4a4dgfWuRVomn215cSKcLq9Tij626AnBJ9brc0bwTEpPmaZX59ddfARhTvyk3Ngl+7xMatWBCo2C9urFJ6yidhjQM+tPc8zvQJzTS/H9dMulWs6alUyt3sL7Nq34pDV1B33qwxmXU1ipG6ZSiqrx+YX9SVJWqbg+vXxj0oXqpFSw7NU+pzJPNg/p1qFCdR5pcBEC3irWj7HRb/aB+Q6o05ua6wSGMy6u34Pra7S07DUntDMDotO70TekIwFUVLqaHt00hO02pMpD2nuAI/+wqw2gR0inSTk/UG0GGKzja8lyDy6iippCiaCWue/uOLwzeD1zI/uMPAJDr68vB3+4EgnXv19+mAXDwxGXszQn63/7c8ezPHQ8U1D1FUTh8+DBHjhwppFOy20Yn+1PYb+Il0SdqZ2VlFWp/LlmypEQyhXniiSdo06YN9erVw+12M2jQIJ566il69gzG/L179+J2u6lSpUrU39WqVYu9e/fG9YwtW7Zwzz330Lx5c3r06MEPP/zA/Pnz2bdvH6+++mqp5FbMcBcyySiKwptvvsnIkSMB2L17N3Xr1uWqq65ixYoVVr4RI0aQlpZW7KLwPXv2FOoZZmdnk5WVxRdffEHXrl2t3qyIeGMSK60oCoqihN7ki7jSQNRBQFLK4LBwqKcdK22aJkIIhta4AVUT+PN1VFfwLWjAF0y/s/dphBAMb38HENwHXnOrmNLE0CUut8pbX/2fpccQ71UYusSd4kIPGMhw2q+zOvAahmEwouVtSGniSXHh9+mY4XR+gFU/PYKUkr5DHw2+hXRr5Pv0qLQQCmvemYWqqgw973Y0t4o/X0eoCqomCPgMVE3w9jf3BfWr9fsonRQFdL+B5lZ5e/dTluwjGkzFCEhcHg0pJUZAUqFKCrf+NYvzO3dAURSGN5mJNEzcXg3db1h6BHw6q35ZZNlmgDoGT5ob/wk/phkcqfDl+lEU+If/FVRVZUjrubg9Gr68AEIouNwqvvygfm9/d79lp0ubzQrKrgmEKkJ6iGC+nxcipeTS5rODOoXequqBoH6mafLutoVWfRvacEZQP0MG7efRMHQDaZi8u2PBWa97pfGVRPiT3+9n69attGrVyvq7s61Tn3GPoxsSlyYwTRPdMK30uldmWvn7Dn4Et1tF1yVSmlHpD9+bjRCCIW1ux5Oi4fcZIX/S8OfrmCa8+939wTrZ58FgfcsP+lM4LYSCyyVY9cGtmKbJ8N4P4PcZCFVBUwV+v4GqKqiq4N1P7kBKycjms4J+5gp+n5F1753tEfWt/s1R/hRZ91b973EURWFQ6jVWjJCGiSfVTcAXQBom7+Uss2w24sL7gvHCDI5OROrx9oa7ME2Ty2tfy9yVU/i/QQuRuonvhD+oh1tjVe4Ky04DtbGomkB1qfjzAsG4pyoE8nVW68Efu6F1JhfpT0ZAsmrPnyxbjsi4CUMP6uT2augBiTSCcW/lrj+hKAqXZtxEIEI/f34gaKdUN2/t+lOZ96d4Y0TTFx9FAB6hkSd1BAoeIciTBtsn3GLp1OQvj6ApAp800BQFVREIReGF8/twfpu2uFwumi15DKEo+KW0RikCpsQtBD9eN92SvfXCx9FNE4+qYkhppXUp2TJrBoZh0Ofj25GYeBSNgGkgMfEqLnymzqf9HrJ0arzkMVJUlTzDQCG4I1WeYSCArdfPQFVVmr/yoKWTIDiaki8NVBR+HHOrJVe7vz+EzzTQFIFAwW+l4YcrbkdKyfBP56BjoKFiYmIgcaEiMXnv4gUFMWLtHNxo6ARldysaekiPDy9+GCEE5793NymKRr5pYGKG0jomsHHQvBLVvWO7WyKUfExTxURDKD5MU8NEpXLdrcG3+vvaIxQ/0tQAEUqHRuaVANT8EQiODPz444+0aNECl8t1TrWNTpZl06ZNdOnShY0bN9KpU/Fz2TZt2kRmZiaNr5lJSu16xeaLl7y9u9j+wgKWL19O69atoz6rU6dOzDMfTm7/Ajz66KM899xzPProozRs2JBPPvmEuXPn8uabb9KvXz9WrFjBxIkTrQ5bmP79+9O0aVOeeeaZmHILIejcuTPjxo3jyiuvpHbt2vErXQwJW1ORaKpXr46mabRpEz1Pt3Xr1qxfX/xx9qcyYHg+YOQuBiVNR84pLG1aUZS40uGFUeEfRAAjYBBeGWIEChw08jAp3V+wdiTgL8gjhLDKidzOMTKtqipSBvuZvoj74XSkjKYJ+T69UFpK08ojpWmdnCsNExkawjT0goBysk6RekTJHtIx4CuY4517LJ8np/6dpV+cHzxiPrQANPK0Xt9J+ln3c/2F0qZZkMeUpvW3Upr48gvrpyiKJU+0naSlkxDCsol+kn5E5AkTqV9kOhl173R85XT8yeVy0bRpU+vHKBk66UaovukFvhWZjszvj7BlZDqqvuXpRaYj/SlcxyLTUpr4fOEphwr+UFoaJv6wPxkR01dEsOMOWD4D0XUvnvoWpV9kLIhYCxGlX6TPFaGHoij8diiXx69+nvwcP2YozkjDtMo/ZdwLFJYrHn+KjAVR6byA9SxfMfr5TvitPGXZn0qSlkCeDNU9TPJkQd0L5zFME8MMxTTTRDeDDfdZ32/gk/YdrfuE3lkGzIJ66I9oTAohgvkAX8Sax8i0qgYb6QA+s8B++WbRcT0vPE0mIi2j9CvQSWKSH0obmFFy+Sz9InzIjPZ/PRTldSJ+cymmHhJR98yi43pexP28YvLEU9+Ekh9KGyiE7aejhGQQQoASrOdCifB5peA7DctumiZNmzZF07So+yenz4Xfp5LuUKWQoC1lQ/9v3br1KTsz8ZKXl8ftt9/Om2++ydChwZkQ7du35+uvv+bRRx+lX79+1K5dG7/fz5EjR6JGK/bv30/37t3jek64s5hIStypmDlzZtx5FyxYUNLiLdxuN126dLGGwsJs2bKFhg3jWyR8MrbcvaCcY0qTfb8cdmxrMxRFKfUubw7nLqY02fvzgWSL4XAGkMDOvBwnFtsMW8fiUu7cVGQ5CSQQCBAIBAp1koIvfYMd28zMTFwuF2vWrGHMmOA00T179vDdd99FbW50KhLdoYBSdCq++uorNm3ahK7r1mEZ4fnPkT20eAJLTk4OP/30k/Xv7du38/XXX1O1alUaNGjAnDlzGDt2LL169aJPnz68//77vPPOO6xdu7akYgNld/cnh+LxpLq4/50bnAOXbIZhGGRnZ9O6dWvHrjYivH7ptgsfjBopdCj7pKgaf+s8wInFNsPOsTiZh9/Fav9efPHFzJkzh5SUFBo2bMi6det44YUXrJf1lSpVYtKkScyaNYtq1apRtWpVZs+eTbt27azdoJJBiTsVw4cPJz09nb/+9a/WkMuRI0eYOHEiF110EbNmzYq7rP/85z/06dPH+nd4FGTChAksW7aMUaNG8cwzz/Dggw8ybdo0WrZsyd///ndroUpJseXhLeUcf16A+8cu4+XvHkm2KA4JRAhBy5YtHZ+1Gf4Tfub1X4D/FNvJOpRN8g2da776iI3nFb9lq0PZw9axOIkjFbHav6+88gpz585l/PjxHD58mIYNG/LHP/6R3//+99bfLFy4EE3TGDNmDHl5efTt25dly5YltfNX4k7FY489xgcffBA1h6tKlSrcf//9DBgwoESdit69e1tbjRXHtddey7XXXltSMW2F4nbFkSlxQ86ySuHzIE4m5etf4ipLj6OsvB5xDsGlFB6CNYF8WRDszNCOE7EQqbH3pddrxN7zHEAeir2nPl5PXGXhisPW5YRk/4j5K8TnU1pO4T37T0avGtsPXHuOxfU89VjsOm5Wia/uEsforagY3372YlfsnVdMBPm5AUxEwUTk0hKnryiVYn8XSnp8OtqdHb+7Na584kTRjZZ8X8SWnmqc51TUzY+Z59e9VWLmAfDUPBEzjyLi23a0LM/iqpixM2YeUXtL3OUlOxafMZLYqYjV/q1duzZLly49ZRler5fFixezePHikgtwhihxTTl+/HiR23bt37+f3377LSFCnSlKvIexwzmPJ8XF/a9f49jWZkRuOehgHzxpbuZ/dqt1LoyDfUhRVV7r1c/xWZth51isJPAqy/j9fjZv3oyux3fQ6akocadi1KhRTJw4kb/97W/s2rWLXbt28be//Y1JkyYxevTo0xboTGLb3nY5xpcX4M4rXnBsazOEELRu3dqxq83w5fq5rftDznoKG5JnGIz55EPHZ22GrWOxCYo8/StZp3KfLidOnGDSpEmkpqZy3nnnWWd+TJs2jfnz55eqzBLXkmeeeYahQ4eSlZVFw4YNadiwIePHj2fw4ME89dRTpRLCwaG0KAp4U50pQ3bEjm/GyjuKAt40T5meWuJQNAqQqp2zu9Q7nAa2jcUJPlG7rDF37ly++eYb1q5dG7XDV79+/Up9+F2JOxWpqak89dRTHDp0yNoJ6vDhwzz11FOkpcWeN5xMbOsY5Ri318Wdf73Ksa3NkFKyefNmx642w53qZt6a6bhTnelPdsOrqizr3tvxWZth61hczjsVb731Fk8++SQ9e/aM2rG1TZs2/Pzzz6Uqs9TjWeGTq1u0aEFaWlrMBdfnAnbbDs0hOP1pzrDnHdvaDFVVadu2rWNXm+HL9TOjw33O9CcbkmcYDP94teOzNsPOsTi8pWwirrLIgQMHqFmzZqH7ubm5pT5vpsSdikOHDtG3b19atGjBkCFD2LNnDwDXXXddiXZ+SgZloePjUDIUoVCrQWXHtjbDNE3y8/Mdu9oMRSjUbloDRTjzn+yGABqklo0XjA7xY/tYXE5HKQC6dOnCqlWrrH+HOxLPPfccF154YanKLPEEyBkzZuByudi5cyetW7e27o8dO5YZM2bw2GOPlUqQs0GZHb6LY56q6U7cXFY9PY7tTytVjKssf7XYZWm/xd6SE8CsVHjbR3eKxpTHRiClDJ422axeXGWJqrG3KDQ88b2ZiWfLTVkhJa6y8DrTQiDoq9u2baNly5ZJe0OWejC+wzJPZMRh2zje+hzpUviNUVEYlWI/L1Atvqmosn6t2JmK2Mq5SCrE3qY5tUl9bn7xKu6/+mV8ecX7vVYrju/CFWfMi2fr2ThkdyjAWy+n0L0UofFIo0usWFyn4aG4ypJHEhfzDD32e9IKab64ympUI7b8l9b+Jq6yfld/Q8w84xr/J66yzjbnQiw+UyTz8LtzgQcffJBBgwbxww8/oOs6jz/+ON9//z0bNmxg3bp1pSqzxCMVH3zwAQ899BD16kU33po3b84vv8R3dkGysJtDOIAvT+fO695wbGszVFWlTZs2jl1thi8vwF2Xv3DKDoVD2SRP6kzY9rbjszbD1rG4nK+p6N69O//61784ceIETZs25YMPPqBWrVps2LCBzMzMUpVZ4tfbubm5pBZxcNjBgwfxeOI84CtJ2Hb4rhwjhEL9JlUxTbPUcwAdzj1M0yQvL4+UlBTHrjZCCIX6LWrwvy0HkNKJx3ZCoNDMW8WJxTbDzrG4vI9UALRr146//vWvCSuvxCMVvXr14oUXXrD+rSgKUkoeeeSRqCPHz0XK7PQnh2JxuVWuvrm7Y1ubIaVk586djl1thsujkXV7X1weZ+tRu+FWVGbW7ub4rM2wdSw2AZmAq4x2KjZt2sR///tf698rV65k5MiR3H777fj9pdtMo8SdikceeYQlS5YwePBg/H4/t9xyC23btuWTTz7hoYceKpUQZwtbDt+Vc3z5OvdPfcexrc1QVZVWrVo5drUZvrwAf7zm1OspHMom+abO73f8w/FZm2HnWFzed3+68cYb2bJlCwDbtm1j7NixpKam8vrrr3PLLbeUqswSdyratGnDt99+S9euXenfvz+5ubmMHj2ar776iqZNm5ZKiLOFM/3Jfgih0KJdbce2NsM0TX777TfHrjZDCIUWneoinN2fbIdAoUNqLcdnbYatY3E5X1OxZcsWOnbsCMDrr7/OxRdfzIoVK1i2bBl///vfS1VmicagA4EAAwYMYMmSJdx7772lemAyseXwXTlHc6mMyOpo7TjiYA+klOzdu5cmTZo4drURmltl+PUXsHjGSvz5erLFcUggbkUwoXp7JxbbDDvHYsU0URLQWUpEGcnANE2rXfzhhx8ybNgwAOrXr8/BgwdLVWaJRipcLhffffddmV2sYzeHcAC/T+fRW993bGszVFWlefPmjl1thj9f57Gb/u50KGxIvmkwc+cax2dthq1jcTkfqejcuTP3338/L774IuvWrWPo0KEAbN++nVq14thuvAhKvFrummuu4fnnn2f+/PmlemAyKbPDd3Gcg2CkxbfXd/6l3WLmcWf/L2ae4+2qx/W8nDqxq1hO/fj2wc+vU/h7EEKhQ9ta1o4jgUpx7kCWFnuvf1OLr/N8omuTmHnEifj2RvfXSo8rn90xTZPjx49TsWLFpL3E8FeI752LGYd4/iqx/bPi1ty4nifdsX/c86vHFw+U/NiL8U60y4irLL1utZh5BCZtuzfiu892II1TjBx7Y/uxmRrf2S/+epVj56nt+F1JqFflaKF7AoV27oZWLN7za9W4ylKrxY6NTeoeKKmIxXJe9b1x5buo+s8x8/yh5cenK45Ft9TYz0sG50IsPmMkaj1EGW1aLlq0iPHjx/PWW29xxx130KxZMwD+9re/0b1791KVWeJOhd/v589//jNr1qyhc+fOpKVFH7K0YMGCUglyNiiznQqHYtE0wSW9mjjbGNoM0zQ5ePAg6enpjl1thKoJeo1uT/a/d566U+FQ5tBQ6Z3SzonFNsP2sbgcNwvbt28ftftTmEceeaTUI1Ml7lR89913dOrUCcBaNR7mXK9wQpR4XbrDOY7fb/DYk58xfNhFyRbFIYEIIc75jR8cSk7Ap/PkjLeSLYbDGcCPzuPH3mGYuCTZojgkEDvHYkUGr0SUYye83vhmjxRFXJ2Kb7/9lrZt2yKE4OOPEzfcd7ZxFmrbD1VV6JpZDyml02m0EVJKjh49SuXKlR272ghVE2T2a8HGD7dg6E48thMqgs6eZk4sthl2jsXl/fA7wzBYuHAhr732Gjt37ix0NsXhw4dLXGZcNeT888+3VoI3adKEQ4cOlfhBDg5nAlUVdOpQJ9liOJwBjh8/nmwRHBKMUAXtL2qCUO3VOHEIdio6ehonWwyHM4BtY3E5X6h97733smDBAsaMGcOxY8eYOXMmo0ePRgjBvHnzSlVmXJG9cuXKbN++HYAdO3Yk7I3/J598wvDhw8nIyEBRFN56661i8954440oisKiRYtK/Ty79bIdgtOf/vTcvx3b2gwhBI0aNXLsajMCPp0/3/EeAZ+z+5Pd8KOz5Phqx2dthp1jsUKCDr8rxbPjaf9mZ2czYsQIKlWqRHp6OhdccAE7d+60Pvf5fEydOpXq1auTlpbGiBEj2LVrV9wyvPTSSzz33HPMnj0bTdO46qqr+POf/8zdd9/N559/Xgqt4uxUXHbZZVx88cU0btwYRVHo3LkzTZo0KfIqCbm5uXTo0IEnn3zylPneeustvvjiCzIy4tuFpDic6U/2Q1MFl/Rq7NjWZkgpOXjwoGNXm6G6BL1Gt0N12a+BUt5REVzsbev4rM2wdSw2zcRdJSRW+/fnn3+mZ8+etGrVirVr1/LNN99w1113Ra13mD59Om+++SavvPIK69evJycnh2HDhmEYRlwy7N27l3bt2gFQoUIFjh07BsCwYcNYtWpViXWCONdUPPvss4wePZqffvqJadOmcf3115Oefvpb8A0ePJjBgwefMs+vv/7KlClTWL16tbWHbrnjeE7MLIGK8W0hebRx7BX9P97VLGYeocfnRLl1Y+f5rX58/fx9XQrr6FYFdc4v2MoyUCG+HQv29a4ZM4+eEl/Dx30k9racptcVV1mug7FtXV44ceIEVavGty3lmSCQFuf7pziyyTiqZW6D1LgeF0iPXZcMT3yy59WvFDNP6k/xzav9rWHsxX1mw1rUz2yI+e1BDH/xP3xmpdjbaMu0+LaPVuJYu2Gq5/YmI+caP/1aOH56FJWRDepb/1ZTAwl73s6DiYsDL3V7Lq58d7d9O2HPnND8s5h5Lm60JWaeZJHsWHymSOaailjt3zvuuIMhQ4bw8MMPW/ciX9wfO3aM559/nhdffJF+/foBsHz5curXr8+HH37IwIEDY8pQr1499uzZQ4MGDWjWrBkffPABnTp14ssvv8TjiXN7/pOIe/enQYMGAbBx40ZuvvnmhHQqYiGl5Oqrr2bOnDmcd955cf3Nnj172LNnT9S97OzsQuVCcFgv3rSiKCiKgmEYCCHiSoefFZlWVdU6xTBW2jRNhBComkCogoBPR3UJhChIh8vXtOAzdF3i0gSmCbohcUXkkVKiCgVDmng0gW6YGKaJR1MJhLZ3NAwDoShI0yRFVfEZBhJIUVXyDSPqBEZFAY9bI9+nR6WFCH4HqqoiFAW3EOQbBqqioAkFnyHRFKVAdqEgFAW/IdFE8HsOGBKXiJZdEwq6NHGrAmma6NJEAZ74/L9cP7hvKK+ClCYet0pAl1FpS7+QPbwuFV/AwAyl8wMGSiiPqqqFdHJpAp/fiNLPNE3cbhW/30BVFVRVBNOaQIR2Q5NSorlU9ICBFnpDqweC98LbHEe+BXJ7VAzdxDAkbo+KrkukYSal7pXGVxLhT6ZpUr9+/aTq5FJFsB6qAhMT3TBxa8G6V8hmLhXDkEHfcoXqnmlG1zePhs+vY5rR6QK9weNWyfcZUelw3YPg9o4et2rVQ00V+AOhuhfxvbg1gV+XaGrIn3SJFrGOIVJ2l0tgymC8cLtUDCkxDNOyE4Dbq6H7jaA/eTUCoXSkzbwhuUzA61bJ9wf9yRM6VyP4/Su8sOTfCKHg8Wj4Qr4V9ouwnSC4XkrVFPy+kD8JJfTcAtmj4p5LIE0wdBllS7c7ZBsj6Kt6KC64Paq1/anHoxIIBO+fbKdIHy2r/pToGKGioAmBTxpoioKqBNOP7fyKG7oHGziaIhCA35S4lOAzAqbErZwU1xWBbko8ioqBRDdNPIqKbkbEbBQkJl6h4ZcGEpMUoeGTekGekB4pQiNP6iiAN5QWEBWzy4udEqVTvXr1LJ87l3Uq8WhKotZDhMo4ua0JUKdOHerUKdm6Tyklq1at4pZbbmHgwIF89dVXNG7cmLlz5zJy5Egg2BYPBAIMGDDA+ruMjAzatm3LZ599FlenYtSoUfzzn/+kW7du3HzzzVx11VU8//zz7Ny5kxkzZpRI5jAlHoNeunTpWelQADz00ENomsa0adPi/pslS5aQmZkZdWVlZQFw4EDwAJ19+/axb98+INgJCd/ftWuXtdp9586dHD16FAieLhheqLRt2zZycoJvk7du3UpeXh4AmzdvxucLHuKTnZ2NrutIKcnOzkZKia7rVoXz+Xxs3rwZgLy8PLZu3QpATk4O27ZtA4ILo8LrWDL7tWDC3cGK02P4eYyd3RuAS8aeb3WgRgxuyYjBLQG4YuR5DLgkuAXc1WM7ROnUq1lwCtntgzLp3LAGAPcP70rbjKqWTs3TKwKw7MLe1EsNnkPy2kX9qOrxWDp5PRqV0r0svKs/ALVrVOCBOX0AaFi3kqVTu5pVefCSrgB0zajB3RdlAnBxwwxrbuDAFvWZfGFbAEad15hrM1sBcFXHZlF2GtU6uAhwSre2DGoWfCM2p2dHZnZvj5SS7du307FV8A3arddfQKvGwRGMeVMuomFGpUJ2evaGS6ic5sHrUnn2hkvwulQqp3ksO52s0z039wKgZZNqUXa6aXrwkJjOF9Rn4o1BXXte3JhxvzsfCNa7UROC2zAPHtOewWPaB3Wd0Im+I1oD0XVvwoyLyOzVCIAbbr+Etp2DAT0Zde/o0aOWnQ4fPmzN1zxw4IBV986UP+3evdvSIxk6XdYjOGKX1acVw7oE3xBdP6Atfds3KKTTjKzOdGwZrHtzJ11Aq8YF/hTW6YE5fahdI/gGfuFd/amU7sXr0SydKqV7eHxuyJ+qV2D+zGDda5RRiXunXGTpNGdKsL51OK8WU68P1rdunepy/TWdLJ2uHRYc0h5yYRPG9Q/WsdG9mzO6d/NCdrpmbAd6dW8Y1O+aTnTLrFvITlMfGEizdrUBmL1oOPWbVbPsFK57T87qR+V0D163ypOz+uF1q1RO9/DkrH6Wne55eBADh7eiUdOq3Hpf8H7z1jWYfnvvQnbKvKghE27uAUCP/s24MuRbkXYactl5DLks+MJp1PiO9BsajIGRde/aG7rSpVswXvxhWnfahTZ2mDGnl1X37ritNw1Cozb/N68vtWoF7fTw/EG28adExogLKtXm/iYXAHBJlfrc3qgzmiK4q1EX/ve/4OGpl1dvzg11grHu6pqtubpmsB7eUKd9lE5Dqgbj+twGXelTKWin/2vUnW4V61g6tU8L/lYtbtKX5imVAfhzs4HUdadH6ZQiNF5uNYwUoVFV8/Jyq2EA1PWkl0s7JUInKSU//fSTtUHPuaxTWPa4kaAYp38R6stkZWUVan8uWbKkZDIB+/fvJycnh/nz5zNo0CA++OADRo0axejRo1m3bh0QnLrkdrupUqVK1N/WqlWLvXvjO+Bx/vz53H777QBcfvnlfPrpp9x00028/vrrpT7gWjHPkRPhFEXhzTffjOqFDR06lE2bNllrKRo1asT06dOZPn16seUUN1KRlZXFF198QdeuXcvcG4ZhDW8udqTi7W0LEULQb/ijQNEjFavfmmnp0emWJ4odqfjmsekYhkGLZx4vdqTi58kzkVLS+7IFpxyp+Oj1GaiqSpv/W1TsSMW3d05DCEH7OxcVO1Lx1f9NtWTvcO8ThUYq0t0aWR1aMGVYsDHWe+zCYkcqPnl9lmWbzD8sKnak4t9PTkNVVXqNeqTYkYqwfqZpMqzvQ8WOVLz38W1IKRnR9f+KHalYtfFeq74N6XBXsSMVq77+v6S/CUpUOpY/BQIB9u3bR0ZGhiXL2dapx01PFDtS8fmSGVb+nlctLHak4tMV0xFC0GvUo8WOVKz9e9A/Lxrz2ClHKj5aMQPTNBk4ekGxIxVr3poV9M+Jjxc7UrH+hQLZ+wx7tNiRio/fnY2iKAxpc3uxIxXvfne/ZbPeExcVO1Lx/+3de3gTVd4H8O9MLm0ppfRG0wBFQIRCa5EKIt5ABUSh66oURRCRdUEX2EpB9EWkvi/ibRdYcRdcRWFBFvUVXHV9FHihXJZFoYXFQoFiWxBoKUIp9JbLzHn/SDMktGmmzWCSk9/nefIwnUwn5/Q7Z5iTM5e8D7PBGMPDI9/E/b/uh398fBCA4DZS8VXefyk53T8g1+NIxdcFuQCA+wYv9DhS8c2el5Qs77/3dY8jFf/cNBeCIGDU6Lc8jlT888ucoG9PWu8jrv/o9SYjFYwxTOucihm3jYRer0fvT17zOFJx6JG5StlT/vcNjyMVx8b9FyRJQt8Nr3scqTg69iWlTr0+fdXjSEXRwy+GXE5a1AkAzpw5A5PJBL1eH9B1KigowMCBA5Gfn688T605BQUFyMjIwI33ZKN9TBePy6lVU3UKB/9vKdauXYuUlBS399SMVFx9/HvmzBl07twZjz32GNatW6csl5mZicjISPz973/HunXrMHnyZKXD5jR8+HD07NkTK1as8LlebdHqh9/9Unbu3InKykokJycr8yRJQk5ODpYuXYqysrJmf6+lAPV6R3WdG3Jbpl2fMtjWaUEQVE07Tz2Q7LJyT3fJJkPClWln2ewu5w3bXKddlhFFEZLs6ENaXJax2K+c26zT6ZTTO+pdLvZxTruWkTGgofEuLq7TssyUZWTG0ND4uxJjkCTHuu0uB4p2+coYpOu0TXYvu72x7FaXJ/FettqxfO9hzMwcqXw2AFhcztd2nXbNo8EmNZlmLstcXSfnelzrJwgCrI3zpcZTlIDGzBrXLYoi7I3rt9uulN3u8vmu25jVIjU77Y9tz5e24kt7MhgMbkPu/qiT85RAm8v2ZnVpN26ZuWRpsXnY3lzueOQ67b69SU2mXbc9QRDctkOr3HTbE0VRKaddcmlbUvNlt7lsk1ZP22TDlfJaGpqWHQAaXNqZc5q5TAuCgJrLVnzyt/1wJctM2c7d9nuSDOcu6Or25GS/al/XbNldyuU2bZGUz7K4tLOrc3IuE8ztSet9hAQGqXHbszMGO3NMv33qIJ7TP9A43yUbl2kru2q/3viehbm0IebehuTGbbhBvpJNvdz8duicz1ymZZX14y0nrerkui8O5Dq5/qxKGy+ybnY9AFJSUlrszKgVHx8PvV6Pvn37us1PSUnBrl27AAAmkwlWqxVVVVVuoxWVlZUYMmSI6s86duwY8vLyUFlZ2eT0sZdffrnVZQ/YTsXEiROVi0+cRo4ciYkTJ2Ly5MltWmerz7cjAc8gihh/Iz1wiTeyLOPs2bNITEykXDmi14u4/6F++HrDIbcOAQl+BkHEE0l9aF/MGZ73xYH68Duj0YiBAwcqp4I5HTt2DN26OU5VzcjIgMFgwObNm5GVlQXAcaZOYWGh28XdLXnvvffwzDPPID4+HiaTSemsAY7OW9B1KmpqanD8+HHl59LSUhw4cACxsbFITk5GXFyc2/IGgwEmkwm9e/f+pYtKCCGEEEJ44qcLALwd/86ZMwfjxo3DnXfeiWHDhuGbb77Bl19+iby8PABAdHQ0pkyZgpycHMTFxSE2NhazZ89GWlpaky/kPVm4cCFeffVVzJ07V7N6+bVTsW/fPgwbNkz5edasWQCASZMmYdWqVZp/XrD2sq29vT+fo9akLsoufz/u+c3Fjn/6LDvneZnpjn8uJav7PFuk9xYrRahaFTqcaG5dEjaWHsULDzuG3GvM6m4pG3XC++0O9fXqvklVdYvP9upuz8bC1d0amHeiKLb6jhlas3ZQeVvWeO/7FWuU9+0yrFrdA+GM1d63XWt7dfs6S6z3dmzrpO7GHFFl9d7XJTP8438LHT+Inv++lk6RXtdV31ndLXitMd7bVG0StbvWMIQ3v63+7cJBvCT+CgAQ26FW1bquN7Xw/02jmKg6Ves6nvWSquWIeoGwL75W/DlS4e3499e//jVWrFiB1157DTNnzkTv3r3x2Wef4fbbb1d+Z8mSJdDr9cjKykJ9fT3uuecerFq1yu3UspZUVVVh7NixrS98C/zaqRg6dKhyuz41PF1HoRad/sQfg17E+HtpyJ03siyjvLwcSUlJlCtHDAYRDz16IzasP+h2DQQJfkZBxG86pdO+mDNc74s1vqaiNdQc/z711FN46qmnPL4fHh6OZcuWYdmyZa3+fAAYO3YsNm3ahGnTprXp95sTsNdUEKIGYwxVly3eFyRBx3ljBcIPmQEXqxogB8Q9B4mWZADn7d5Hq0jw4XVfHKjXVFxLb7/9tjJ9/fXXY/78+dizZw/S0tJgMLifddGaxzk48bmleMBdL5vALjF88a8fMW8qZcsTURSRmJjo72IQjUl2Gd9+dcTfxSDXgJ3J+Pj8EfwP/T/LFa73xQwQtBgwDaJOxZIlS9x+bt++PbZv3648/8JJEATqVHhDpz/xx6gXMeX+VBpy54wsyzh16hS6dOlCuXLEYNBh/OQBWPdhAWwut68lwc8o6DAzKYP2xZzhel8sA5oMmwbRoaXzwYHXSkh1Kgh/ZMbw45lqfxeDXAPt2qm7EJcED5kxlJVcUJ6FQ/ghg+Fo/QV/F4NcA9zui688xsf39QQ55/UdrreVbQvOup0t466XTWCXGDbtO0HZckYURcTHx1OunJHsMrZv+VF5mCfhh53J+LLqOLVZzvC8L3ZeU6HFK1itXLkSqampCA8PR3h4OFJTU/H++++3eX38bSUtoNOf+GM06PDc2AGULWdkWUZZWRnlyhmjUYepM2+F0ajulockeIQJOszvMoTaLGf43hezK3eA8uUVpEMV8+fPx+9//3uMGTMGn376KT799FOMGTMGzz33HF56qW23Z6bTn4LApeu8P+OgtrPKIauo9l4XabguxusyPca18LyLVrJ0Undu9eXkpnXUiQy7zlUiy7nMdeoatz3S+0FNXSfvz58AAFt77+uqSVb3MA57FN0v36lDhw5+/fyDi59TtZysYi/aEOu9fap9VoK1g7rtUo2oUu/3/790XbiqdTUkeN9PXTKF4fuSc7hkCoPUwrnMst7730tW2S/Z8Y85Xpe53NW3If9QE9Wuock8HUTst5Th8caf945apGpd39611Osydyb92IrSEa35e198rYTi3Z9cLV++HO+99x4ee+wxZV5mZiZuvPFGzJgxAwsXLmz1OkOqU8Hj8F2ok2SGrSWnKVvOiKKI2NhYfxeDaEySGXYdOO3vYpBrQIKM7bUltC/mDNf74hC/pkKSJNx8881N5mdkZMBuV/cw1quFVOvnc/gutIXpRCy8dyBlyxlZlvHjjz9SrpwxGnSYO/kWGA10+hNvjIIO8zsNpzbLGZ73xYLMIEgavIL0wTsTJkzA8uXLm8z/61//iscff7yZ3/AupEYqfL2qnQQeu8zw1dET+PVdt/m7KERDgiAgPj6e2ixnJEnGlj1lkCT+DlBCncRkfHP5CH4l3OXvohAN8bwvFhiDoMGd6LRYh7+sXLkSmzZtwuDBgwEAe/bswU8//YQnnngCs2bNUpZbvHixqvVRp4IENYkx7PmpkrLljCAIiI6O9ncxiMYkmSG/6Ky/i0GuAQkMe+t/on0xZ7jeF4f46U+FhYUYMGAAAODHHx3XLSUkJCAhIQGFhYXKcq1p0yHVqZAketgSb8L0OiwaPgiSJEGno1MqeCFJEkpKStCjRw/KlSNhBh1enDIYr63cAws9/I4rYYIeCxJH0L6YM1zvi5W7N2mwniC0bds2zdcZUtdU0AVk/LFJMv62/xhlyxlRFGEymShXztjsMj7ZdAQ2ek4Fd2xMwrqLBdRmOcP1vlirZ1QEZ5/imgipkYpgHZatTfJe7jCVDzL98UmT12VOTvb+DeLG2/6s7vNmz/K+kA9kxvCfivNKtlGl6jK2Rqm5xae6dVmivS/377/nqFrX/219UdVyvBMEAVFRUf4uhiq1nb0vY/F+l2bUqNzeLnfzfkvZmm6qVoU6s/dbHdvbqWwHHb1/i3kpScD3DVVAkgDA83ov3uC9jg0x2n1revhVdbcPJg5rUld5fE8QZmj+eX9I/1jzdRJ1gmlf3CZBOsoQqDjsenpGpz/xJ1yvw/LMOyhbzkiShCNHjlCunAk36LDs0TsQTnd/4g5jEZAvraU2yxme98WCrN2LOIRUp4LL4bsQZ5VkLP7XQcqWM6IoIjk5mXLljNUu4+2tB2Gl0584ZIHQ7lVqs5zhel8sM+1eBACd/kSCnMwYis9XU7acEQQB7dq183cxiMZkxlBcWe3vYpBrQBBkQF9E+2LO8LwvFqDRLWXpogoFh11Pz3gcvgt1EXodVj08jLLljCRJOHz4MOXKmQiDDu9PHIYIOv2JO4y1g1y9gdosZ7jeFzvv/qTFiwAIsU4Fl8N3Ia5BkvDS5u8pW86IoogePXpQrpxpsEtY8OX3aLBzeIAS8hogtH+O2ixnuN4Xyxq+CAA/dyp27NiBMWPGwGw2QxAEfP7558p7NpsNc+fORVpaGiIjI2E2m/HEE0/gzJkzbf48GpblD2PAqUu1lC1nBEFAeHg45coZxoDTF2vpiz0OCYIMQXeC2ixneN4XO5+orcWrtVo6/r3a1KlTIQgCli5d6jbfYrFgxowZiI+PR2RkJDIzM3Hq1KlWl0VLfr2mora2Funp6Zg8eTIefvhht/fq6upQUFCA+fPnIz09HVVVVcjOzkZmZib27dvXps8L1uE7Q633ZdTc1hIAWC/vKxvfb6+6lWkktutFVcvVdW2aX4ROh08HD1ceuNSQoO4z936g3a1uG2L429n6myRJKCoqQkpKSsA/cMkWb/O6TENPi9dlaq9Tt3+qTVK1mCoX+nj/XqlgubrbrdaYvbcDWz/g08HDMXbPZtS3sD+uSfb+n3T+e3QbWH/p2/V0k3nONiuZ6eF3PAmmfXGb+OkbjpaOf119/vnn+O6772A2m5u8l52djS+//BLr169HXFwccnJyMHr0aOTn5/stK792KkaNGoVRo0Y1+150dDQ2b97sNm/ZsmUYNGgQTp48ieTk5FZ/HpfDdyGuQZIwae827OvXz99FIRoSRRG9e/emNssZZ3ttCNIveIhn1Gb5xHWufnyidkvHv06nT5/G9OnT8e233+KBBx5we6+6uhorV67EmjVrcO+99wIA1q5di65du2LLli0YOXJkq8ukhaC6+1N1teMuPx07dvS4THl5OcrLy93mFRUVAQBkWXb7VxRF1dOCIEAQBEiSBFEUVU07P8t1WqfTgTGmapoxBlEUoRMF6EUBFrsMvShAFARYJce0c/2Gxs+wyTKMoggZDHaZwagTlWVkWYYeAuxgCBN0sDMZEhjCBR2szFFfSZIgQASDDD3CIMEKBgY9wmCHVSmjr3Vy/fvqITrqxCTo4fgb2pgEA64quyDAzhjCRBESY7AzBqOog1WWrpRFECAxhnCdDlZZhswYwvU6WKUr9VObkyA4noNRb5MgCgKMOhENdse0c2SEMYYwvQiLXW42p7Zub4Gy7fladl/q5BxuD/Q6AUCY2NieGEO42LjtgbmVPULUwyLbITdON8h2MNdtEkC4Tod6SYIIIEyZFhCmc6yDNW7bDZIEndC4vUmOtqETr2xvRp2o7CNECLDKsrKPcN0mATj2EcyxvwjTNbYtmbV624vQ69Bgl8DguIFCvV1y1Emvcyt7nWSHAAHhOlGpn1G8Uj/W+B+0o04iLJIEvdC435Nlak8BWifBZX+n0+m4qBOPObV22luugVQn1/2aKlpdD9G4DuexpqukpCQkJbV+eFmWZUycOBFz5sxBv2a+NM3Pz4fNZsOIESOUeWazGampqdi9e7ffOhVB0/VsaGjACy+8gPHjx6NDhw4el3v33XeRkZHh9powYQIA4Ny5cwCAs2fP4uzZswAcnRDn/FOnTuHCBcejqU+ePImLFy8CAEpLS3Hp0iUAQElJCWpqagAAxcXFqK+vBwAcPXoUFovjFIeioiLY7XbIsoyioiLIsgy73a5scBaLBUePHgUA1NfXo7i4GABQU1ODkpISAMClS5dQWloKALijlxnZ96YDAIb37Yppd6UCADL7d1c6UI+nXI/HU64HAPzmxj54uFd3AMDMm1Ld6jS0o2OE55Xk23BLlGNDf6P7XbgxMkGpUywc51JlYjai4Jj/sPASIhClWZ0uXryIkydPAgDujeqF38bdAgAYE90XE2MGAADGdrzRLaesLj0BAM/1uhGjkxyPDJ6fchP+fsu9kGUZpaWlGGR2lPe1uwchrVMsAOBPw4fg+pgOrc6pc4dI/PnhOwEAPeM64I+/GgIASE2KdavTvOEZAIA7e5oxa6gjpxF9uuLZ2x05nTt3Tskp2LY915wuXLignK95ret07NgxHDlyRKlHoNfp1T6DcGvHRADA4r5D0D86rkmdVqfegy7h7QEAn6bfh1hDOCJEvVKnWGMYPh0yHADQpV17rB40DABwQ1Q0lmfcodTp9WGDAACDzAlYcLtj27urmxnPD05X6jTzJse293Cv7vhNWh8A7vsI1zpNvyUV913fFQAw+7Z03HmdWXVOrtve6qy7ERMRhgi9Dquz7kaEXoeYiDCszrpbyWn1zcPw6eDhSIuOwfKbHHXq3zEOf0wf3CSnYWYzXkx31On+5GRkp1J7CuQ6ybKMo0ePclUnHnNqbZ2cuZ4/fz7g6+Qsu1qCLGv2AoAJEyY0Of589913W1UmpzfeeAN6vR4zZ85s9v2KigoYjUbExMS4zU9MTERFRUWbPlMLAnN+LeRngiBg48aNePDBB5u8Z7PZMHbsWJw8eRJ5eXktdio8jVRMmDAB33//PQYOHBh03zAMmLXU40jF3j/MhCiK6Pfanxx/K7npSMXB52co9ej32WseRyqOjp0HSZLwyuEsjyMV/536mebfMAz6+mWPIxX/vi9XWf6G95c0GakIE3UwigIOPvV7MMaQtvBtjyMVh1/OblVOaXOXeByp2L9ohlKnm//rbY8jFQWv/z4kv93ytU52ux0AlG/GArlO3Vf9weNIRcmkHIiiiO5r3/A4UvHj+DkQRRE9VvyxxZGKot8+B8YY+i9a5nGk4uA8x/Z241vLPI5UHHrhyjbZb8GfPI5UHHrl963a9jJeftvjSEX+f88EYwx9Vy4FAFgkGWFXjVQcefo5JadeS5d6HKkozs6m9hSAdXK2W1EUaaSCozp5yzWQ6lRQUICBAwciPz8fAwY4vpxsTkFBATIyMnDr9U+hQ4TvF6ldqi/Hv49/gLVr1yIlJcXtPTUjFVcf/+bn5+OBBx5AQUGBci3Fddddh+zsbGQ37v/WrVuHyZMnKx02p+HDh6Nnz55YsWKFz/Vqi4A//clmsyErKwulpaXYunVrix0KoOUAnRuw89+2TLte/NLWaUEQVE0rp3/IDFLjExvtMgNwZdpZNpt8ZQzP6jotyW71tjf+roVdOae5wWVap9OBNY7l2XFlY3VOu5bRlzq5/k3tkJ1Vcpu24aqyN/Z/LW51ldBeH6YsIzUu43rOtuvtK1uTE2NAva3x1CrGlPXIjLnVydL4dGBPOfmyvfl729Oq7G2pk2vHItDrZJFdtje5+e2tXrY3O61sb4By8bLsNs2UaUEQlG1bYgyS1Li9MQa7dGV7c57u57oduu4jXMvuXBZwHOw3t4yabazepZ05p5nLtLPsscYwNEgWt/o5/2auOTnqJCnTzvOWqT0FZp3YVfnwUKdrNR1MdVKbayDUyfVnVRg0uqbC8U9KSkqLnRm1du7cicrKSrdrhyVJQk5ODpYuXYqysjKYTCZYrVZUVVW5jVZUVlZiyJAhPpehrQL69Cdnh6K4uBhbtmxBXFycT+uTXf5TJXwI1+mweuAwypYzziF3ypUvzvYa7nJAQfhAbZZPXOcqa/jS0MSJE3Hw4EEcOHBAeZnNZsyZMwfffvstACAjIwMGg8Hthkbl5eUoLCz0a6fCryMVNTU1OH78uPJzaWkpDhw4gNjYWJjNZjzyyCMoKCjAV199BUmSlPPEYmNjYTQaW/15OvqPjDv1koTR//oGZY3nWxM+6HQ6pFKm3HG2V8IfarN84jpX1rZnTDS3ntZq6fg3OTm5yZfoBoMBJpMJvXv3BuC4Q+qUKVOQk5ODuLg4xMbGYvbs2UhLS1PuBuUPfu1U7Nu3D8OGDVN+njVrFgBg0qRJyM3NxRdffAEA6N+/v9vvbdu2DUOHDm315wXI5SOtdjHV7nUZdoOKh1kAGNKt1Osy/5O2UdW6tFJw/0JVy5VNnd1kHmMMFosFjDEIgoD6zt7/VmoVvqnuPviFb9H98rXmzDUsLEwZ6g5UJ5563usyw/oc9bpM2TNNt+/mHH3Z+/ZmSfT+7AwAaOik3Vdshxd5L1fpb3NU5fpj4/8FJHgEU5sl6vGdq0a3lEXr19HS8e+qVatUrWPJkiXQ6/XIyspCfX097rnnHqxatcqvX6D7tVMxdOjQFg/0te4EcDl8F+JkWUZJSQl69+5NI1EcoVz5RLnyi7LlE9e5yszx0mI9reTt+PdqZWVlTeaFh4dj2bJlWLZsWas//1oJ+Au1tcRdgyDQ6XTo27evv4tBNEa58oly5Rdlyyeuc5UZoMWXzVp0TDgR0Bdqay1YT38injHGUFdXR9lyhnLlE+XKL8qWT3znyq48VduXVxtOf+JVSHUq6PQn/siyrDx4ifCDcuUT5covypZPXOfqPP1JixcBQKc/kSCn0+nQp08ffxeDaIxy5RPlyi/Klk9c58pkx0uL9RAAITZSwefwXWhjjOHy5cuULWcoVz5RrvyibPnEda7Oh9/5/PJ3RQJHSI1UBOvw3YnfzPG6TPHYl1Sta/Wglb4WJ6DIsoyKigr06NEDOp1O1d+KBL6rcw12Hw788Bf9vBNTvN/mFgBKZ+Rc45K44y1XcgVlyyeuc2UanbrEY4erjUKqU8FdgyDQ6XTo1auXv4tBNEa58oly5Rdlyyeuc1UutNZgPQQAnf5EghxjDNXV1ZQtZyhXPlGu/KJs+cR1rrKs3YsAoE4FCXKMMfz888+ULWcoVz5RrvyibPnEda6MadOh4PFv00YhdfqTKIZUHyokiKKInj17+rsYRGOUK58oV35RtnziOlc6/UlzIXWUHawXahPPZFnGhQsXKFvOUK58olz5RdnyietcNbnzk0YdE06EVKeC8OnSpUv+LgK5BihXPlGu/KJs+cRtrs67P/n6ok6Fgk5/IkFNFEVcd911/i4G0RjlyifKlV+ULZ+4zpUxME0efkedCqeQOsrmcvguxMmyjJ9//pmy5QzlyifKlV+ULZ+4zlWLUQrniwAIsU4F4VNdXZ2/i0CuAcqVT5QrvyhbPnGbK11ToTk6/YkENVEUkZyc7O9iEI1RrnyiXPlF2fKJ61xlGZAkbdZDAITYSAWXw3chTpZlnD17lrLlDOXKJ8qVX5Qtn3jOlcmyZi/iEFIjFYRPdrvd30Ug1wDlyifKlV+ULZ/4zVWrU5fo9CenkOpU0OlP/BFFEZ07d/Z3MYjGKFc+Ua78omz5xHWuMoMmHQK6UFvh16PsHTt2YMyYMTCbzRAEAZ9//rnb+4wx5Obmwmw2IyIiAkOHDsWhQ4fa/Hk8Dt+FOlmWUV5eTtlyhnLlE+XKL8qWT1znyhjAZA1ere9UtHT8a7PZMHfuXKSlpSEyMhJmsxlPPPEEzpw547YOi8WCGTNmID4+HpGRkcjMzMSpU6d8/av4xK+ditraWqSnp+Odd95p9v0333wTixcvxjvvvIO9e/fCZDJh+PDhuHz58i9cUkIIIYQQwgvGGJiswasNnYqWjn/r6upQUFCA+fPno6CgABs2bMCxY8eQmZnptlx2djY2btyI9evXY9euXaipqcHo0aMhaXHxeRv59fSnUaNGYdSoUc2+xxjD0qVLMW/ePDz00EMAgNWrVyMxMRHr1q3D1KlTW/15dPoTf0RRRFJSkr+LQTRGufKJcuUXZcsnrnNlMgAtHn7X+nW0dPwbHR2NzZs3u81btmwZBg0ahJMnTyI5ORnV1dVYuXIl1qxZg3vvvRcAsHbtWnTt2hVbtmzByJEjW18PDQTsNRWlpaWoqKjAiBEjlHlhYWG46667sHv3bo+divLycpSXl7vNO3DgAACgsLAQwJXToERRVD0tCAIEQYAkSRBFUdW087Ncp3U6HRhjqqYZY60qYyjWyW634/z580hISAAALurEY06trZNrrs6yBHudeMyptXWSJAkXLlxAXFwcAHBRJx5zakudAODcuXOIi4uDXq/nok485tTaaW+5BlKdjhw5AgCor6+HGjVSNZgGnYpaOM6eKSoqavJeUlKSJp2y6upqCIKAjh07AgDy8/Nhs9ncjpHNZjNSU1Oxe/du6lRcraKiAgCQmJjoNj8xMREnTpzw+HvvvvsuXnnllWbfmzRpknYFJIQQQgghAaWsrAy33Xabx/fj4+PRrl07HKr7XrPPNBgMmDBhQpP5CxYsQG5urk/rbmhowAsvvIDx48ejQ4cOABzHyEajETExMW7LJiYmKsfP/hCwnQonQRDcfmaMNZnnaurUqU3OO/vuu+/w7LPPYuXKlejfv/+1KCbxk6KiIkyYMAFr165FSkqKv4tDNEK58oly5Rdly6dgyrW+vh5lZWVev6VPTk5GUVERfv75Z80+W5Ik6HS6JvN9HaWw2Wx49NFHIcsy/vKXv3hd3tsx8rUWsJ0Kk8kEwNEbcw2lsrKyyeiFq5aGmvr3748BAwZoW1ASEFJSUihbDlGufKJc+UXZ8ilYcm1phMJVcnJywD8p3GazISsrC6Wlpdi6dasySgE4jpGtViuqqqrcRisqKysxZMgQfxQXQAA/Ubt79+4wmUxuF6tYrVZs377dr38wQgghhBBCrhVnh6K4uBhbtmxRrkNzysjIgMFgcDtGLi8vR2FhoV+Pkf06UlFTU4Pjx48rP5eWluLAgQOIjY1FcnIysrOzsWjRIvTq1Qu9evXCokWL0K5dO4wfP96PpSaEEEIIIaRtWjr+NZvNeOSRR1BQUICvvvoKkiQp10nExsbCaDQiOjoaU6ZMQU5ODuLi4hAbG4vZs2cjLS1NuRuUP/i1U7Fv3z4MGzZM+XnWrFkAHBdUr1q1Cs8//zzq6+vx7LPPoqqqCrfccgs2bdqEqKgofxWZEEIIIYSQNmvp+Dc3NxdffPEFADS5Dnjbtm0YOnQoAGDJkiXQ6/XIyspCfX097rnnHqxatarZazt+KQJry1M7gkx5eTneffddTJ06ld/7LYcoypZPlCufKFd+UbZ8olxJa4REp4IQQgghhBBy7QTshdqEEEIIIYSQ4ECdCkIIIYQQQohPqFNBCCGEEEII8Ql1KgghhBBCCCE+CYlOxV/+8hd0794d4eHhyMjIwM6dO/1dJNIKubm5EATB7eV84jrgeCx9bm4uzGYzIiIiMHToUBw6dMiPJSbN2bFjB8aMGQOz2QxBEPD555+7va8mR4vFghkzZiA+Ph6RkZHIzMzEqVOnfsFakOZ4y/bJJ59s0oYHDx7stgxlG1hee+01DBw4EFFRUejUqRMefPBBHD161G0ZarPBSU221GZJW3Dfqfj444+RnZ2NefPmYf/+/bjjjjswatQonDx50t9FI63Qr18/lJeXK68ffvhBee/NN9/E4sWL8c4772Dv3r0wmUwYPnw4Ll++7McSk6vV1tYiPT0d77zzTrPvq8kxOzsbGzduxPr167Fr1y7U1NRg9OjRkCTpl6oGaYa3bAHgvvvuc2vDX3/9tdv7lG1g2b59O373u99hz5492Lx5M+x2O0aMGIHa2lplGWqzwUlNtgC1WdIGjHODBg1i06ZNc5vXp08f9sILL/ipRKS1FixYwNLT05t9T5ZlZjKZ2Ouvv67Ma2hoYNHR0WzFihW/UAlJawFgGzduVH5Wk+PFixeZwWBg69evV5Y5ffo0E0WRffPNN79Y2UnLrs6WMcYmTZrEfvWrX3n8Hco28FVWVjIAbPv27YwxarM8uTpbxqjNkrbheqTCarUiPz8fI0aMcJs/YsQI7N6920+lIm1RXFwMs9mM7t2749FHH0VJSQkAx6PtKyoq3DIOCwvDXXfdRRkHETU55ufnw2azuS1jNpuRmppKWQeBvLw8dOrUCTfccAOefvppVFZWKu9RtoGvuroaABAbGwuA2ixPrs7WidosaS2uOxU///wzJElCYmKi2/zExERUVFT4qVSktW655Rb87W9/w7fffov33nsPFRUVGDJkCM6fP6/kSBkHNzU5VlRUwGg0IiYmxuMyJDCNGjUKH330EbZu3Yo//vGP2Lt3L+6++25YLBYAlG2gY4xh1qxZuP3225GamgqA2iwvmssWoDZL2kbv7wL8EgRBcPuZMdZkHglco0aNUqbT0tJw6623omfPnli9erVy4RhlzIe25EhZB75x48Yp06mpqbj55pvRrVs3/POf/8RDDz3k8fco28Awffp0HDx4ELt27WryHrXZ4OYpW2qzpC24HqmIj4+HTqdr0muurKxs8u0KCR6RkZFIS0tDcXGxchcoyji4qcnRZDLBarWiqqrK4zIkOCQlJaFbt24oLi4GQNkGshkzZuCLL77Atm3b0KVLF2U+tdng5ynb5lCbJWpw3akwGo3IyMjA5s2b3eZv3rwZQ4YM8VOpiK8sFguKioqQlJSE7t27w2QyuWVstVqxfft2yjiIqMkxIyMDBoPBbZny8nIUFhZS1kHm/Pnz+Omnn5CUlASAsg1EjDFMnz4dGzZswNatW9G9e3e396nNBi9v2TaH2ixRxT/Xh/9y1q9fzwwGA1u5ciU7fPgwy87OZpGRkaysrMzfRSMq5eTksLy8PFZSUsL27NnDRo8ezaKiopQMX3/9dRYdHc02bNjAfvjhB/bYY4+xpKQkdunSJT+XnLi6fPky279/P9u/fz8DwBYvXsz279/PTpw4wRhTl+O0adNYly5d2JYtW1hBQQG7++67WXp6OrPb7f6qFmEtZ3v58mWWk5PDdu/ezUpLS9m2bdvYrbfeyjp37kzZBrBnnnmGRUdHs7y8PFZeXq686urqlGWozQYnb9lSmyVtxX2ngjHG/vznP7Nu3boxo9HIBgwY4HbbNBL4xo0bx5KSkpjBYGBms5k99NBD7NChQ8r7siyzBQsWMJPJxMLCwtidd97JfvjhBz+WmDRn27ZtDECT16RJkxhj6nKsr69n06dPZ7GxsSwiIoKNHj2anTx50g+1Ia5ayrauro6NGDGCJSQkMIPBwJKTk9mkSZOa5EbZBpbm8gTAPvzwQ2UZarPByVu21GZJWwmMMfbLjYsQQgghhBBCeMP1NRWEEEIIIYSQa486FYQQQgghhBCfUKeCEEIIIYQQ4hPqVBBCCCGEEEJ8Qp0KQgghhBBCiE+oU0EIIYQQQgjxCXUqCCGEEEIIIT6hTgUhhBBCCCHEJ9SpIIQQTqxcuRIjRozwaR2VlZVISEjA6dOnNSoVIYSQUEBP1CaEEA5YLBb06NED69evxx133OHTumbNmoVLly7h/fff16h0hBBCeEcjFYQQwoHPPvsM7du397lDAQCTJ0/GRx99hKqqKg1KRgghJBRQp4IQQgLIuXPnYDKZsGjRImXed999B6PRiE2bNnn8vfXr1yMzM9Nt3pNPPokHH3wQixYtQmJiIjp27IhXXnkFdrsdc+bMQWxsLLp06YIPPvjA7ffS0tJgMpmwceNGbStHCCGEW9SpIISQAJKQkIAPPvgAubm52LdvH2pqajBhwgQ8++yzLV4vsXPnTtx8881N5m/duhVnzpzBjh07sHjxYuTm5mL06NGIiYnBd999h2nTpmHatGn46aef3H5v0KBB2Llzp+b1I4QQwie6poIQQgLQ7373O2zZsgUDBw7Ef/7zH+zduxfh4eHNLnvx4kXExMRgx44dbqc/Pfnkk8jLy0NJSQlE0fEdUp8+fdCpUyfs2LEDACBJEqKjo/H+++/j0UcfVX531qxZ2L9/P7Zt23YNa0kIIYQXNFJBCCEB6A9/+APsdjs++eQTfPTRRx47FABQX18PAM0u069fP6VDAQCJiYlIS0tTftbpdIiLi0NlZaXb70VERKCurs7XahBCCAkR1KkghJAAVFJSgjNnzkCWZZw4caLFZePi4iAIQrMXVhsMBrefBUFodp4sy27zLly4gISEhDaWnhBCSKihTgUhhAQYq9WKxx9/HOPGjcPChQsxZcoUnD171uPyRqMRffv2xeHDhzUrQ2FhIW666SbN1kcIIYRv1KkghJAAM2/ePFRXV+Ptt9/G888/j5SUFEyZMqXF3xk5ciR27dqlyefX1dUhPz/f5wfpEUIICR3UqSCEkACSl5eHpUuXYs2aNejQoQNEUcSaNWuwa9cuLF++3OPvPf300/j6669RXV3tcxn+8Y9/IDk5WZNnXhBCCAkNdPcnQgjhRFZWFm666Sa8+OKLPq1n0KBByM7Oxvjx4zUqGSGEEN7RSAUhhHDirbfeQvv27X1aR2VlJR555BE89thjGpWKEEJIKKCRCkIIIYQQQohPaKSCEEIIIYQQ4hPqVBBCCCGEEEJ8Qp0KQgghhBBCiE+oU0EIIYQQQgjxCXUqCCGEEEIIIT6hTgUhhBBCCCHEJ9SpIIQQQgghhPiEOhWEEEIIIYQQn1CnghBCCCGEEOIT6lQQQgghhBBCfPL/7amLotbzZH8AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# plot pseudosection\n", "c_map = 'viridis'\n", "\n", "fig,ax = plt.subplots(figsize=(8,3))\n", "swam.plot('pseudosection', method = 'FK',\n", " cmap = c_map, axes = ax,\n", " vmin = 120, vmax = 270,\n", " fmin = 10.0, fmax = 25,\n", " width=3\n", " )\n", "fig.suptitle(f'min_offset = {min_offset}, max_offset = {max_offset}, method={tr_type}, wlen = {w_len}, wmove = {w_m}', fontsize=10)\n", "fig.tight_layout(rect=[0, 0, 1, 0.95])\n", "fig.savefig(f'{prj_dir}/proc/masw2D/04_figs/pseudo_masw2D_smooth.png', dpi=150)" ] }, { "cell_type": "code", "execution_count": 30, "id": "76fcc423-46af-493d-a748-5e6587d2e6d9", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Save dispersion curves corresponding to (SIN,REP) = (43, 2) to \"../data/real_data/proc/masw2D/03_proc/proc1_FK/Mode0\" ..... 1.1 s\n" ] } ], "source": [ "# save DCs\n", "swam.save(method = 'FK', # tr_type\n", " dc_mode = 0, apply_to = 'all', use_windows=True)" ] }, { "cell_type": "markdown", "id": "7ba3f22c-9933-467a-8334-72bdb3e172a2", "metadata": {}, "source": [ "## Comparing raw and smoothed dispersion curves\n", "\n", "Here's the pseudosection after picking \n", "\n", "\n", "and this one after smoothing \n", "" ] }, { "cell_type": "markdown", "id": "2e45924a-2c68-4601-bfd0-55ebeb5f45c0", "metadata": {}, "source": [ "```{image} ../../data/real_data/proc/masw2D/04_figs/pseudo_masw2D.png" ] } ], "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.23" } }, "nbformat": 4, "nbformat_minor": 5 }