{ "cells": [ { "cell_type": "markdown", "id": "5e58fd47", "metadata": {}, "source": [ "## CDF of mis-registration" ] }, { "cell_type": "code", "execution_count": 1, "id": "97b46fd3", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Go to directory: /Users/yunjunz/Papers/2021_Geolocation/figs_src/stats\n" ] } ], "source": [ "%matplotlib inline\n", "import os\n", "import string\n", "import numpy as np\n", "import seaborn as sns\n", "import pandas as pd\n", "from matplotlib import pyplot as plt, ticker, patches\n", "from mintpy.objects import timeseries, sensor\n", "from mintpy.utils import readfile, utils as ut\n", "from ipynb.fs.full import utils\n", "plt.rcParams.update({'font.size': 12})\n", "\n", "work_dir = os.path.expanduser('~/Papers/2022_Geolocation/figs_src/stats')\n", "os.chdir(work_dir)\n", "print('Go to directory:', work_dir)\n", "\n", "proj_dirs = [os.path.expanduser('~/data/geolocation/ChileSenAT149/mintpy_offset'),\n", " os.path.expanduser('~/data/geolocation/ChileSenDT156/mintpy_offset'),\n", " os.path.expanduser('~/data/geolocation/KyushuAlos2DT23/mintpy_offset')]" ] }, { "cell_type": "markdown", "id": "9b2e9a7f", "metadata": {}, "source": [ "### Calculate" ] }, { "cell_type": "code", "execution_count": 2, "id": "429381f2", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "max residual for S1: 0.25 m\n", "max residual for A2: 1.64 m\n" ] } ], "source": [ "# read from HDF5 file\n", "ts_S1_list = [[], [], [], []]\n", "ts_A2_list = [[], [], [], []]\n", "for proj_dir in proj_dirs:\n", " suffix = '' if 'Alos2' in proj_dir else '_S1Bias'\n", " fnames = [os.path.join(proj_dir, f'timeseriesRg{suffix}.h5'),\n", " os.path.join(proj_dir, f'timeseriesRg{suffix}_TECjhr.h5'),\n", " os.path.join(proj_dir, f'timeseriesRg{suffix}_TECjhr_SET.h5'),\n", " os.path.join(proj_dir, f'timeseriesRg{suffix}_TECjhr_SET_ERA5.h5')]\n", " tsDict = utils.read_ts_files(fnames, print_msg=False, print_max=False)[1]\n", " for i, ts in enumerate(tsDict.values()):\n", " if 'Sen' in proj_dir:\n", " ts_S1_list[i] += ts.tolist()\n", " else:\n", " ts_A2_list[i] += ts.tolist()\n", "\n", "# calculate CDF\n", "# link: https://stackoverflow.com/questions/24788200/calculate-the-cumulative-distribution-function-cdf-in-python\n", "num_S1 = len(ts_S1_list[0])\n", "num_A2 = len(ts_A2_list[0])\n", "p_S1 = 100. * np.arange(num_S1) / (num_S1 - 1)\n", "p_A2 = 100. * np.arange(num_A2) / (num_A2 - 1)\n", "\n", "print('max residual for S1: {:.2f} m'.format(np.nanmax(np.abs(ts_S1_list[-1]))))\n", "print('max residual for A2: {:.2f} m'.format(np.nanmax(np.abs(ts_A2_list[-1]))))" ] }, { "cell_type": "markdown", "id": "6a339c34", "metadata": {}, "source": [ "### Plot" ] }, { "cell_type": "code", "execution_count": 11, "id": "c2995b74", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "save figure to file /Users/yunjunz/Papers/2021_Geolocation/figs_src/stats/stats_cdf.pdf\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcgAAADMCAYAAAD3RzNIAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABdGklEQVR4nO2deXgUVdaH30NIQkjCFgzIvgiCMCDLyCooKDgsiqDoiCOggrgDOi6oI+qo+I3oMOIGOoO74ojIMoKo7AKCCGJk35RVEAIECGQ53x/VXXQ33Ul10mtZ7/PUk1ruvXV+fevk1nLvuaKqODg4ODg4OHhTJtoGODg4ODg4xCJOA+ng4ODg4OAHp4F0cHBwcHDwg9NAOjg4ODg4+MFpIB0cHBwcHPzgNJAODg4ODg5+cBpIBwcHBwcHP4SsgRSRu0RklYicEpEpPse6i8gGETkhIvNFpK7HMRGR50TkN9fyfyIirmNlReRDEckWkc9FJN0j3yMiMipU9js4ODg4OHgSyifIPcDfgX977hSRqsA04DGgCrAK+MgjyXCgH9ASaAH0AW5zHesPKFAVOOreLyL1gb7ASyG038HBwcHBwSRkDaSqTlPV6cBvPof6A1mq+rGq5gJjgZYi0sR1fDAwXlV3qepuYDwwxHWsPrBAVfOB+UAD1/5/Afe79js4ODg4OIScSHyDbAasdW+o6nFgq2v/Wcdd6+5jPwLdRCQJuBTIEpGrgYOquiTchjs4ODg4/H4pG4FzpAEHfPYdAdI9jh/xOZbm+g75P+BijNeyy4EPga+AHiLytOvYj8BIVT3teYLU1FQtKCgAoEyZMmRkZFC+fHkAypYtS0pKCseOHSMvL4+kpCTS0tI4ceIE7jypqank5eVx+rRRbEFBAWlpaZw8eRKAxMREkpOTycnJMc+RmppKdnY2CQkJfssoV64cIsLJkycpKCigXLlyfss4fvw4hYWF5nnLlStHXl4eACkpKagqubm5ACQlJZGYmEh2djaJiYkkJCRQvnx5rzLS0tI4deoUeXl5pg5/ZRw/fhyAhIQE8vLySEhIwB2rNz09nZMnT5Kfbzy0ly9fnoKCAk6dOkVBQQHly5cnISGBEydOnPUbuyksLCQpKclvGQDJyckkJCRw5MgREhMTzypDRLzqqaCggAoVKgT8jd31lJuba9aJ+zf2V9cnT54kISHBbxme9VRQUEClSpUC/sbuejp06BApKSl+f+MyZcpQvnx5cnJyyMvLo0yZMqSkpJCbm0t+fj6qSlJSEoWFhZw6dQpVJT8/n6SkJHJzcyksLERVKVOmDKdOnfI6t6rmqKr5vd4KVv0FIC8vjypVqhTpL1avdd8yjh49atZVqOpBVTl+/DgJCQl+r3V//vLbb7+RnJzsVUYgf3HbkZOTU6y/nDhxgoSEBPNaD4W/uMtw2xHIXzzr6fjx4yQmJhb7G5cpU4b8/HwSExNLXdf+6sntm76/cdmyZRERjh07Rn5+PoWFhZw+fdr0A7ePuP+nlZLA/qKqIV0wvkNO8dieALzik2YdMMC1fgS4yONYG+BYgLL/AYwAegFzAQHeAEb4pm3Tpo1aoXHjxpbSDRs2LKbTRUtHOMqM5TrJz8/XwYMH6549e3Tz5s26Zs0aXbp0qc6ZM0eff/55HThwoF500UV64YUXalJSkjZs2FBr166t1apV08qVK2tqaqomJiYqxrf1cC2rNEi/teovqtbrJ5jfNZbSqganMVZsDiZ9NPXl5ubq+vXrdcaMGdquXTu9/fbb9aqrrtIOHTpogwYNNC0tLaz+kZqaqg0aNND27dvrVVddVaS/ROIJMgvjOyMAIpIKNHTtdx9vCXzr2m7pcQyPfM2BjsCDwF+B71RVRWSlK09Y6du3b0yns0o4zmsXLe50GzZsYMKECaxcuZLjx497Le6797feestSmVu3brWUzs6E41oKZ9pgiRWbw6WxJDafPn2abdu2sWXLFjZv3szmzZvN9Z9//tl8kgRYsWJFqewTETIyMsjMzKRatWpkZmYyd+5c7r//fnPb/TczM5PU1NSz8gcsW0M03ZWIlMV4Zfs4UAsYBuQDlYEtwM3AbOAJoKuqtnflGwHcC1yG0cLPA15S1dc8yhZgAXCfqq4SkYHAXa48U4DVqvq8pz1t27bVVatWFWt3mzZt+O6770qsO1awiw6IrJb8/HwWL17M999/z7p161i2bBkbN26MyLkTExNJSkryuyQnJ5+1r0KFCqSnp5Oenh5wPS0tjbZt236nqm2DscWqv4C9rrVA2F1jKPSdPn2a7du3ezV+7vWdO3d6NYLBkpiYSLVq1c5q4Pz9rVq1KmXLej/rBaNPRAL6SyifIB/FaBzd3Ag8oapjRWQAMBF4F1gBXO+R7nWM3qnrXNtvuPZ5MhT4UVXdHjwNo3fsAYxvk77pLXPVVVeVNGtMYRcdEBktqsrEiRN57rnn2L17t+V8IkL58uVJTU09a6lTpw4XXXQRf/jDH0hNTeXNN9/k3nvv9dvoJSYmFnnnGsvY6VoLhN01WtWXm5vL9u3b2bp1K1u2bPFqCEvaCIoIderUoVGjRpx33nk0atSIOnXqeDV6FStWLJV/hKr+QtZAqupYjCEc/o59CTQJcEyBB1xLoLL/jcf4SjWGd1wfKH0wdOnSJRTFRB276IDSaSksLDQd+dChQ2RnZ3P48GEOHz7stZ6VlcXBgwf9lpGSkkLPnj0ZOnQoDRs29GoEU1JSLDvuNddcQ6NGjUqsJVax07UWCLtr9NR39OhRtm7dajaCnn937dpFSd4yuhtBdwPo/tuoUSPq169PuXLlQinnLEJVf5H4BungEBKOHz/O6tWr2bt3L7/99huHDh3y+rt//35++ukns3dhMNxyyy20aNGCFi1a0K5dO7P3qYNDvKOqHDx40KvhW7p0KWPGjGHLli0cOOA7yMAaIkLt2rXPagDPO+88GjRoEPZGMBLYtoE8cuQIw4cPp2/fvkV+ZL7gggsiaFX4sIsO8NaSn5/P22+/zauvvsr3338fii7dJhUqVKBdu3b87W9/o3PnziEr102k62TmzJnMnDkToGKwea36C9jrWgtEvGrctm0bX3/99VlPg57DR4KhTJky5pNgw4YNvRrCWG4ErdSfFX8JWSedWMNqp4O9e/dy7rnnRsCi8GIXHXBGy3//+19GjhwZ1DdCgGrVqtG8eXPOOeccKleubC6VKlUy12vVqsV5551HmTLhi5URrTopqtNBIILppGOnay0Q8ajxo48+4i9/+Ys5BtQqSUlJNGjQwKsRdP+tW7cuSUlJYbI4fARTf5HqpBOXbNy4Me4cwR920QHGMIslS5YwcOBAr/0iQrNmzWjcuDEZGRlkZGRQpUoVr/VGjRpRvXr1KFnujZ3qxBO76vIk3jSuXLmSIUOGBGwc09LSvBq+/Px8evfuTcOGDalZs6YZnMEuhKr+fvcNpJuHH36YatWqMXLkyGLT9u/fn+HDh3PFFVeE37DfCcePH2fhwoVkZWXxwANn99d64oknuOuuu6hSpUoUrHPwxdNfFixYwI033siuXbv8ph09ejSNGzdmxIgREbby98GePXvo16+fGYWmYcOGDBo0yKtBPOecc7w6ly1YsIBLLrkkShbHEYEiCMT7YjUyyKZNm/TXX3/VGjVq6IkTJyzlWbFihbZu3dpS2kixadOmaJtQIrZs2aL33nuvVqhQIWDki88++yzaZpaIaNUJYY6ks3z5ci9/mT9/vtasWTNg+j179mitWrX01KlTJdITDeLJn7p37276SuXKlS3ZHk/6SkIw+oryF8tPkCLSA7gQI3aqZwP7t5C01FGiVq1aTJw4kV69elnuuXjRRRdx9OhRVq1aRdu2QX3qCRu1atWKtglBUVhYyL///W/uuusuMzqNL+XKlWPLli3UrFkzwtaFhnirE6t89dVXQfnLueeeS5MmTZgxYwbXXHNNmK0LDfFSd9nZ2Xz11VeA0aFm6tSploYWxYu+khIqfZZ6KIiIe5B/G6C2xxL3v/KKFSv4/PPP6dq1q7nv8OHD9OnTx+zk0adPn7NeH11yySXMnj070uYGpLThmiLJ6dOn6d27N8OGDfNqHBs1asTtt9/O3XffzZdffsnRo0fjtnGE+KqTYPjvf//r5S9unnnmGapWrUq9evV47733vI7Fmr8UR7zUnTtgOEBGRgaXXXaZpXzxoq+khEqf1S58fwbaqOp1qjrUY7k5JFZEmXXr1nH++eeb24WFhQwdOpSdO3fy888/k5KSwl133eWVp2nTpqxdu9a3KIciOHHiBB9//DHdunVjzpw55v66devy+eefs2HDBl555RX69+9P9+7dSUxMjKK1DoHYtm2bl78A7Nu3j4MHD7J7927eeusthg8f7hWyz/GX8OA57CmcPbJ/r1j9RX8DssNoR8hxj+tyjXMJSFpaGtnZ2aSnn5ntJCMjgwEDBlC+fHnS09N55JFHWLhwoVe+9PR0srOzw2F6iUhLSys+UZQ4fPgwEydOpGHDhgwcOJClS5eax6677jp++OEHrrjiCtPBY1lLMERax8yZMxk+fDiUYhxkcf4CkJOT4+Uvbp566imSk5Pp2rUrvXv3ZurUqeaxWPOX4oiHa7CgoMCrU6G/OglEPOgrDVb0WfKXQB8nPRfgNmAG0AEjbqq5WMkfjSWYTgfnnHOOfvvtt+b28ePHdfjw4VqnTh1NT0/X9PR0BTQ/P99MM378eL366qstn+P3wunTp/W///2vDh8+XFu1aqUVK1YM2Pnm6quv1ry8vGibbDsIcycdX3+ZP3++Vq1a1SvN/fffryNGjDC3P/nkE23VqlXJBDmcRUFBgQ4ePNjLn55//vlomxWXFOUvVp8gXwX6AEsxZuZwL5st5o9Zli1bRosWLdi0aZO5b/z48WzcuJEVK1Zw9OhRFi1aBOC+WQBg/fr1tGwZ9lm2LLNs2bKonr+wsJDFixfTpk0brrnmGiZNmsT333/PkSNHvNJVrVqVxx57jKVLl/LJJ5+cFYUfoq8lVNhFhy9169b18hcw3hJ4hvj7+eefqVGjhrkda/5SHLFcd3l5edx2221eU67dfffdjB492nIZsawvFIRKn6UGUlXLBFjifnTpqVOn6NWrl9cr1GPHjpGSkkKlSpU4dOgQTzzxxFn5Fi5cyJ/+9KdImlokgXqCRoLnnnuOmjVr0qVLF9atW3fW8ZSUFJo1a8bQoUNZs2YNTz75JB07dgwY9DuaWkKJXXT48sc//vGsTw4Ajz/+OKdPn2bx4sXMmjWLa6+91jwWa/5SHLFad7t37+bSSy/ljTfeMPfdeuut/POf/wxq9otY1RcqQqUvqEABIlIHqAnsUtVfQmJBDHDTTTdx4YUXcvLkSVJSUhg5ciQ33HADVatWpUaNGtx3331Mnz7dTL9y5UpSU1O56KKLomd0jPDcc8/x0EMPnbX/z3/+M8OGDeOCCy4gMzMzbqd2cjibnj17cuedd5r+AlC9enUqV65MjRo1KF++PK+99hpNmhgT+Ozdu5effvqJfv36RdHq+Gf+/Plcf/31/Prrr+a+v/zlL7z22mtOB51wEejdq+cCnAssBE4De11/FwE1rOSPxmL1m4r7G9jDDz+sL774oqU8/fv319mzZ1tKGymi8S3vxIkTmpiY6PUdZNiwYbp69epSlWuX75LR0kGYv0Hm5eUF5S+jR4/Wl19+OQgF0SfWrsHnn39ey5QpY/pZmTJl9P/+7/+0sLCwROXFmr5QE4y+ovzFagM5HfgXkOraTgUmADOs5I/GYtXhs7KyLP+QsUykdfzyyy86YsQIr8Zx7969ISnbqZPSEe4G0i71UxSxpHHy5MlefpaZmanz588vVZmxpC8cBKOvKH+x+lzeGbhPVY+7njqPY0xw3LGkT67hxmq3dc/XFfFMJHVs27aNpk2b8tprr5n7hg4dGrIg4U6dlIxIDfOwS/0URaxoXLRoEXfccYe53bFjR77//vtSx1GNFX3hwoo+K/5i9RvkYeACwHOk7/nE8NjIihUrMmnSpGibYRvWrFnD6tWr2b59O3//+9+9jjVq1Ignn3wySpY5uHHP5Th58uQjxaf2xvGX2GPNmjUMGDDAnKHjwgsv5IsvviA1NTXKltkDK/5itYH8P+BLEXkT2AnUBYYCj5XezOjSvHnzaJsQEsKh4/jx43z22We8/vrr5lAXX0aOHMnTTz9N+fLlQ3Zep05iG7vq8iRaGk+cOMHUqVOZNGmS11CFzMxMPvvss5A1jnavw1Dps9RAqupkEdkK3AC0APYAf1bVr0NiRRSxS3fnUOkoLCxk3rx5LF68mIkTJ541jtGTTp06MX78+JD3oHPqJLaxqy5PIq1xzZo1TJo0iffee4+jR496HUtKSuLTTz+lTp06ITuf3esw4sM8XI1h3DeIvmzevDmuA2K7CZWOwYMH8+677/o91qBBA/785z+bs4937NgxLN3LnTqJbeyqy5NIaMzLy2PKlClMmjSJVatWnXW8bNmy9OvXj4ceeog2bdqE9Nx2r8NQ6QvYQIrII6r6tGs94AcmtTjdlYjUA17BCFd3CvgvMFJV80WkO/AyUAdYAQxR1Z2ufDcA44FcYKiqLnDtbwi8A1ysqgU4lJpZs2ad1TiWK1eO2267jWuvvbbIwf0ODg7Bcdddd/n97nveeecxbNgwBg8eTLVq1aJgmYObop4gPaeyqh2Cc70C/IoxprISMA+4Q0TeB6YBtwIzgaeAj4D2IlIWGAe0xphqayLgfrn8L2B0aRvH2rVDIS36lFbH+vXrGThwoNe+d955h8suuyxkvVOt4tRJbGNXXZ6EW+OpU6d4//33ze2kpCQGDBjA8OHD6dq1a9hvRO1eh6HSF7CBVNXbPdaHhuBc9YGJqpoL7BOROUAzoD+QpaofA4jIWOCgiDTB6D27W1X3isiXGAHSEZFrXPuXl9You9yhlUbHtm3b6NWrFydPngSMqChz5syJWuxMp05iG7vq8iTcGhctWkROTg4A9erVY9WqVWRkZIT1nJ7YvQ5Dpc/SN0gROaSqVfzs/1VVMy2eawJwvYgsACoDf8LoBXsJHsNHVPW4q0NQM+BTIENEagGtgCwRSQMeBboXdbK9e/eac9YlJyczdOhQWrVqBRjTWTVr1sy8SCtVqkTnzp1ZvXq1+YG8bdu27N+/n19+MSLqNWrUiOTkZH788UfA6FXWuHFjlixZYp6jQ4cOrFq1yrzw27Vrx65du9i9ezcA559/PgkJCfz000+A0RDVr1/f7K2WkpJCu3btWLFihdlYdejQge3bt7Nv3z4ALrjgAgoKCsy59mrWrEmtWrVYsGABaWlppKWl0bZtW5YtW2Z+qO7cuTObNm0yxwY1b96cU6dOsXnzZjZt2sTo0aO9Ak0/9thjtGzZkiVLlpCfnw9Aly5dyMrK4rfffgOgZcuWHDt2jG3btgGGk1epUoXVq1cDULlyZVq2bMnChQtRVUSErl27snbtWg4fPgxA69atOXToEDt27ACMb5zp6eksXbqUtLQ0r3oC45tMPNXT2rVrzTkt3fXknsg12HoC4664WrVq5veqChUq0Lp1a5YsWcKnn37KrFmz3FMeVSVIrPoLGNNd9enTJ27qwZ+/FFcPS5YsoVy5ckHXQ3H+smXLFr755huvCaXbtGnDrl27yMjIKJG/uOfZDMZf1q9fT1paWtzXUyB/ycnJoWnTpqX3l0ARBDwX4JiffYnAb1byu9I3Bb4D8jEiQkwBBHgTGOeTdinGd0gwGsLlGKHuLgReAG7BaFjnA3OB5r7nsxoZpLQRKWKFkug4ffq0tmzZ0itKx6233hp644Lk91wnoYAwR9KxS/0URag15ubm6ptvvqlNmjQ5a9q3r776KqTnsoLd6zAYfUX5S5FPkCKy2FWJ5UTEdyBcLeCbovJ7lFPG1ZC9jhF9Jw34N/AckANU8MlSATgGoKpfAV+5ymkBtAX+CuzAiPBTG3gDaG/FFl8qVPA9dXwSrI7Dhw8zaNAg8+6zbNmyfPTRR1x99dXhMC8ofq91Ei/YVZcnodJ49OhRXn/9df75z3+yZ88er2NJSUncddddXHrppSE5VzDYvQ5DpU+MBjTAQZHBGE95rwIjPA4psB/4WlXzij2JSFXgAFBJVY+49vUD/o7R2WawqnZy7U91pW2tqhs8yhCMp8h7MAKmf6+qNUQkGTisql4j1du2bav+uk47GN8cr7jiCvOVEcC4ceN48MEHo2iVQ6gQke9UtW0weRx/CQ2qyubNm5k7dy5z585l/vz5nDhxwitNhQoVGDFiBPfee6/XnJkO0aFIfwn0aOm5AE2spCumjG3AQxjfPSthfF98DzgHOAIMAMphPFUu95N/GDDBtV6WM+HvrgB+9E1v9ZXR4sWLLT+KxzJWdZw8eVKbNWvm9Yrn9ttvj6no/r+3Ogk1hPkVq13qpyiC0Zidna3Tpk3T2267TevVq3fWK1T3Ur16dR03bpxmZ2eH0XJr2L0Og9FXlL9YDRTQUUT8BiZX1X9bLKM/8E/gQaAA4/vhKFU9ICIDMIZwvIsxDvJ6z4yuJ9B7cQVHV2Ps5F0YgQtyMcLelQj3R/V4x4qO3NxcevbsSVZWFmB8gH/nnXe8JraNBX5PdRKP2FWXJ0VpLCws5LvvvjOfEpctW0ZBQeDRZk2bNmX06NH85S9/ITk5ORzmBo3d6zBU+qw2kH/x2a4ONMToTGOpgVTVNRgda/wd+xJoUkTeg5wZ/+je9x7GE6iDRSZMmOAVU/WJJ56IucbRwSFW2bBhA6+++irvv/8+Bw8eDJguLS2NSy+9lJ49e9KzZ0/OO++8CFrpEEqK/AZZZEaRm4GmqvrX0JoUGqx+UyksLLTFbNxF6VBVxowZw7hx48x9t9xyC5MnT47JyDi/hzoJJ+H+BmmX+ikKt8a8vDxmzJjBK6+8wtdfB4602apVK7NB7NixI0lJSRG0NnjsXofB6Cv1N0h/C1AGOFTS/OFezjvvPB02bJjOmDGjyPfPP/zwQ7HvqOOBQDoKCwv13nvv9foWkpGRoQcPHoywhdaxe52EixkzZuiwYcMU2Kxh8hdV+9RPUXz55Zc6duxYrVGjht/viZmZmXrjjTfqO++8o/v27Yu2uUFj9zq0os+KvwTTGHouaRi9WrdZyR+NxRkHaTBu3LiznHvt2rWRNS5I7F4n4QZnHGSJ2bVrl1533XVapkyZs/ymTJkyevXVV+u8efO0oKAg2qaWCjvXoWqExkF64B7c78lujJ6lDjHKypUrefzxx83tZs2asXDhwoiGtHJwiBfmzZvHoEGDOHDggNf+atWqMXz4cIYNG2b7GKYO3lhtIOv7bB9Xo+NM3BOteKOhxlfH2rVr6du3rxmWqX79+nz77bchndg4XNi1TuyC3XQVFBTw5JNP8tRTT7nfmAFGuLg777yTfv36xfw3xWCxWx36Eip9lr5iqupONaafygcyAdtcLceOHYu2CSHBU0dOTg49evRg//79gBEXde7cuXHROII968RO2ElXdnY2PXv25MknnzQbx+rVq/PBBx+wcOFCBg4caLvGEexVh/4IlT5LDaSI1HGFndsJzAZ+FpElIlI3JFZEEXew7XjHU8crr7xiBvAFmDlzJo0aNYqGWSXCjnViJ+yka+zYsXz11Vfmdrdu3VizZk3Ep3iLNHaqQ3+ESp/Vfr5vYQQar6jG7B2VgJWu/Q4xxLRp07xCxt1888106tQpihY5OMQu8+fPN9cfeOABvvjiC9tPBeVgHasNZBvgr6p6HEBVczAi4rQJl2Gl5ciRIwwfPpyZM2cWma5evXqRMSjM1KtXjyVLlpw16XE8xle1U51EkpkzZzJ8+HCAisHmteovYJ/6UVWvJ42//vWvJCQkAPbRGAhHnzV/sRQoQES+AJ5Q1aUe+zoCY1W1hwV7I47Vgc9Hjx6N+8j2qsrKlStp166d1/7PP/+cK664IkpWlRw71AlET0e4AwXYpX5+/fVX82kxPT2dI0eOmIEz7KIxEI6+MxTlLwGfIEXkSfcCbAX+JyLvi8hzIvI+8D9gSwlsjyncE/zGIwcOHGDo0KGce+65ZzWOU6ZMicvGEeK7Tjyxiw5f7KJrzZo15nrjxo29okrZRWMgHH3WKGqYh++An2muv5nAKYzZOMqFxAqHoCksLOSyyy7jhx9+OOvY008/zeDBg6NglYND/PDtt9+a6xdddFEULXGIVQI2kKpa4hky4onKlStH24QSMWTIEK/GMTU1lQ4dOnDHHXfQr1+/6BkWAuK1Tnyxiw5f7KJrxYoV5rrvGxi7aAyEo88aARtIEamnqjtc6w0CpVPVuO4vHI8DZhcsWMA777xjbjdo0IDNmzfbJvhwPNaJP+yiwxe76Nq+fbu53qJFC69jdtEYCEefNYr6j7rOY30LsNn113PZ7CdfXLFw4cJomxA0//rXv7y2V65cyeLFi6NkTeiJxzrxh110+GIXXbm5ueZ6amqq1zG7aAyEo88aRb1iTfdYt8ejiR+s9OKNJfbu3evVFX/BggVUqVIl7nQUhV202EWHL3bQpaqcOHHC3PadyNgOGovC0WeNYhs+EUkQka0iEhtTYVvE6riuWJwPsSj+9re/mbNl169fny5dugDxp6Mo7KIl0joiNQ7SDvXzxBNPsHfvXgDKli1LpUqVvI7bQWNROPpCOw5yE3CRqmZbNzG6BDOuK14oLCykdu3a7NmzB4CXX36ZO+64I8pWOcQa4R4HGe+88cYbDBt2ZiKiESNG8Oqrr0bRIodoUqJxkD78E/hIRLqKSEMRaeBeQmZllFi7dm20TbDMqlWrzMYxKSnJy8njSUdx2EWLXXT4Es+6Pv74Y0aMGGFu9+zZ86xv+hDfGq3g6LOG1emuJrr+Xu6zX4GEkFgSJQ4fPhxtEyxx6NAhevbsaW5ffvnlJCYmmtvxosMKdtFiFx2+xKMuVWXcuHGMGTPG3Ne6dWs+/vhjLz9yE48ag8HRZw2r012VCbAE1TiKyPUisl5Ejru+a17s2t9dRDaIyAkRme85S4iI3CAie0Vku4hc4rG/oYh8IyJx3UBbIT8/n27dupGdnW3uu+CCC6JnkINDHHHq1CmGDh3q1Tg2atSI2bNnk56eXkROh989qlrsAvwrwP5/WsnvSns5xnRZ7TEa5pqupSpwBLgWIzLPP4DlrjxlgZ+Bc4E+wI8e5c0G2gc6X5s2bdQKR44csZQuWhQWFuq1116rGE/rCmivXr103759XuliXUcw2EVLtHQAq9SiX2qQ/qIaX/Vz5MgRvfjii738p2vXrnrw4MFi89kZR98ZivIXq98ghwTY/xeL+QGeAJ5U1eWqWqiqu1V1N9AfyFLVj1U1FxgLtBSRJkAGsFtV9wJfAg0AROQa1/7lQZzfL4cOHSptEWFl6tSpfPzxx+Z2//79mT179llT8sS6jmCwixa76PAlXnQVFhZy0003eY0Rvvnmm/niiy/IyMgoMm+8aCwpjj5rFPkNUkRudqfzWHfTADho5SSu16BtgRkisgXjSXE68FegGWB+UVXV4yKy1bX/UyBDRGoBrYAsEUkDHgW6F3XOvXv3cv755wPGGKehQ4fSqlUrADIyMmjWrBmLFi0iJyeHXbt20blzZ1avXs3Ro0cBaNu2Lfv37+eXX34BjFcyycnJ/PjjjwBkZmbSuHFjlixZYp6jQ4cOrFq1ipycHMAIX7Vr1y52794NwPnnn09CQgI//fQTYMxcXr9+fZYtWwZASkoK7dq1Y8WKFZw8eZL8/HxGjhxpamrbti0TJkxg7969bNy4EYCaNWtSq1YtfvzxR3bs2EFaWhpt27Zl2bJlnDp1CoDOnTuzadMmcxLl5s2bc+rUKTZvNuI81K5dm2rVquHuxVihQgVat27NkiVLzCElXbp0ISsri99++w0wIlUcO3bMnC6oXr16VKlSxQwSXLlyZVq2bMnChQtRVUSErl27snbtWvP7QOvWrTl06BA7duwAjIhA6enpphbPegKjO34s1hNAhw4d2L59O/v27QOMV+CbN282tbnryR3eLNT19OmnnzJr1iz3K8OqBIlVfwHIycmhXr16MV8P7777Lp999pmpcfjw4dx3330UFBSwYMGCIuthw4YNZt3Fur+4O6QE4y/r169nx44dMVFPYPhLQUHBWf/XSuovOTk5FBQUlN5fAj1aGk+ezHct+R7r84GvgQ8o4hWnTzk1MF5vrMJ4XVoVWAo8DbwJjPNJvxQY4lrvDiwHFgIXAi8AtwCXuGyZCzT3PafVV0bz58+3/CgeaYYMGeL1amjjxo0B08ayjmCxi5Zo6SDMr1jjoX7mzJmjImL6zqhRo4LKHw8aS4Oj7wxF+UuRT5CqeimAiPxdVR8tKm0xnHT9fUmN16WIyAsYT4KLAN+JuyoAx1w2fAV85crTAuNJ9K/ADqAzxqwjb2B82wyaBg1ic6TK22+/zZQpU8ztUaNG0bhx44DpY1VHSbCLFrvo8CUedN13331mNJUuXbrw3HPPBZU/HjSWBkefNax+g5zgerXpjqwzVERuEhGrvWAPA7sw7uZ8yQLMyLIikgo0dO3HY79gDDe5B+MJNEFVdwIrAe9Iw0EQq73YJkyYYK537tyZZ555psj0saqjJNhFi110+BLruvLz883XfQAfffSR36EcRRHrGkuLo88aVhvIWUAj1/ozwP3AaGB8EOf6D3C3iGSKSGVgpKvcT4HmIjJARMoBfwN+UNUNPvlvBb5X1TXAb0CKiFwAXAqUeEaRWBwwm5eXx4YNZ+S///77lCtX9NSbsaijpNhFi110+BLrug4ePGg+PVatWpXq1asHXUasaywtjj5rWA0U0BhY41ofBHQEcjCe8kZZLOMpjCe/TUAuMBV4WlVzRWQAxtPhu8AK4HrPjCJSFbjXdV5UNV9E7sL4FpoL2Gruyh9++MEMpFynTh1q1/adu9rBwSEQ+/fvN9d9e3s7OASD1QayAEgSkcbAEVX92fV6Nc3qiVQ1D7jDtfge+xJoUkTeg0Bzn33vAe9ZPX8giuvuHQ3mzp1rrnfo0MFSnljUUVLsosUuOnyJdV27du0y188999wSlRHrGkuLo88aVhvIzzGe+DKAD137LgB2h8SKKNKsWbNom+BFdnY248aNM7fbt7fW9yjWdJQGu2ixiw5fYl2XexgEQN26dQMnLIJY11haHH3WsPoN8laMyDVvAs+69lXFGNQfk1idvsc9ZihWmD17NseOHQOMKVuuvfZaS/liTUdpsIuWSOuI1HRXsVo/mzZt4pFHHuHpp58295W0gYxVjaHC0WfNXyw9QarqKWCSz74FVvJGi4oVKzJp0qTiE8YYH374obk+bNgwatasGUVrHOKJvn370rdvXyZPnnwk2Lzx6i/Z2dl89NFHvPXWW+bAdE9atChxB3cHm2PFXwI2kCIySVWHu9bfwf8QDVT1plJbGkXKlrX6ljn8/Prrr17fHy+++GLLeWNJR2mxixa76PAl2roKCgqYN28eU6ZMYfr06WZkFU8yMjK4/fbb6d27d4nOEW2N4cbRZ42AEyaLyMOq+qxr/fFABajqEyGxJMTE4wSwc+fO5YorrgCMsEgHDx4MevyWg4OdJ0xWVXr27Mm8efPOOla2bFl69+7N4MGD6d27N0lJSVGw0CHeKNGEye7G0bX+RKAlHAZHEncsxFjAc3Dz9ddfH1TjGEs6SotdtNhFhy/R1LV8+fKzGsdWrVoxYcIE9uzZw/Tp07n66qtL3Tjate7cOPqsUdQr1m5WClDVr0NiSZRwB/GNBdwBgyH4+R5jSUdpsYsWu+jwJZq6PL/Rd+vWjRdffDEs3xntWnduHH3WKOpF7Zs+2zUxvkP+hjHcQzDCx9k7qF+EOHz4MNOmTTO3mzdvXkRqB4ffH3l5eUydOtXcHjNmjNMJxyGsBGwgVbW+e11ExmA0io+p6gkRKQ88idFYxjVt2wb1qSZsjBkzhuzsbMCYcsZqgAA3saIjFNhFi110+BItXU888YQ5PVK1atW45JJLwnYuu9adG0efNayOgxwFPKSqJwBcfx/GiMcak1gd1+UZlipa7Ny5k9dee83cfuGFFyhfvnxQZcSCjlBhFy2R1hGpcZDRqJ+FCxd6Bez/61//SkJCQtjOZ5drMBCOPov+EmgeLM8FY2qpTj77OgI7reSPxhJP80F26tTJnLeuZs2amp+fH3QZsaAjVNhFizMfZGg4dOiQ1qpVy/SRbt26aUFBQVjPaZdrMBCOvjMU5S9WB4s8BswRkZnALxhzMPYB7rSY3yEAixcvZunSpeb2PffcE9Y7YweHeOPll18246tmZGTw9ttvU6aM1ZdfDg4lx+p8ju8A7YD1GJMZbwDau/bHNY0aNSo+URiZPXu2uV63bl3uuOOsWO6WiLaOUGIXLXbR4UskdRUWFvKf//zH3B4/fnxEokvZte7cOPqsYTncgKr+BPxUbMI4Izk5OWrnPnjwoNdM54MGDSItzfIEKV5EU0eosYsWu+jwJZK6Fi1axLZtxnSvlSpV4rrrrovIee1ad24cfdb43b+n8Bx7GGk8x3QB9OjRo8RlRVNHqLGLFrvo8CWSujyHdQwaNKjYicNDhV3rzo2jzxq/+wYyWhw4cIAXXnjB3L7nnnvo2rVrFC1ycIg9fv31V3O9S5cuUbTE4feIbRtIq93WMzMzI2TRGbKzs7nqqqvYvn07AAkJCYwaNapUZUZDR7iwi5ZI64jUMI9I6srLyzPXI/la0C7XYCAcfdb8JWCwcr+JRWoDNVV1ueVMUcJq8OX8/PyIRrb/8ccfGThwIOvXrzf3vfrqq4wYMaJU5UZaRzixi5Zo6Qh3sPJI6vrTn/7EnDlzAKNDW69evSJyXrtcg4Fw9J2hRMHKfQqoIyJLMXqvfunad42IvGHR3phlyZIlETvX4sWLadWqlVfj+Oyzz5a6cYTI6gg3dtFiFx2+RFLXzz//bK5XrVo1Yue1a925cfRZw+or1teB2UA64H7nMQ+4PCRW/A5Ys2YN/fr1Iz8/39z3xhtv8NBDD0XRKgeH2OXUqVNs2rTJ3A42gL+DQ2mx2kBeBIxT1UJcEyer6hGC/NYhIo1EJFdE3vXY111ENojICRGZLyJ1PY7dICJ7RWS7iFzisb+hiHwjIqUeUR+J7xorVqygffv2HDp0yNz3j3/8g1tuuSVk57BTt227aLGLDl8ipWvTpk3mDWX9+vVLPASqJNi17tw4+qxhtYHcD5znuUNELgB+9p88IC8DKz3KqApMw4jUUwVYBXzkOlYWGAe0Bu4GJnqU8y9gtKoWBHn+swg2KHhJeOihh8xZz5OSknjmmWe4//77Q3qOSOiIFHbRYhcdvkRK1+bNm831pk2bRuScbuxad24cfdaw2kA+D8wSkaFAWRH5M0ZD9lzR2c4gItcD2cBXHrv7A1mq+rGq5gJjgZYi0gRj9pDdqroX47tnA1c517j2h6SjULhnUf/6669ZsGCBuT1nzhwefvjhkJ8nHmaDt4pdtNhFhy+R0uXu5Q3QoEFkZ9Wza925cfRZw1I3H1X9t4gcAoZjxGK9CWPqq+lW8otIBYzpsboDnu8VmwFrPc5zXES2uvZ/CmSISC2gFZAlImnAo65yQkJOTk6oivLLu++ab5Pp168fl156aVjOE24dkcQuWuyiw5dI6XJH0AHjFWsksWvduXH0WcNSAykiCa7GcHoJz/MU8Kaq/iIinvvTgAM+aY8A6apaKCK3A/8FTgHDMBrZl4A/iMjjwGngPlU9K2zC3r17Of/88wHjffTQoUNp1aoVYAQ8btasGYsWLSInJ4clS5bQuXNnVq9ebc5E3bZtW/bv388vv/wCGLH9kpOTzQgNmZmZNG7c2OwtlZycTIcOHVi1ahU5OTmoKseOHeOdd86Eqx0yZAi//vorP/1kROyrXr069evXZ9myZQCkpKTQrl07VqxYwcmTJwHjVcH27dvNefAuuOACCgoK2LhxIwA1a9akVq1a5OTksGDBAtLS0mjbti3Lli0zX+t27tyZTZs2mYOumzdvzqlTp8xXWLVr16ZatWrmXVeFChVo3bo1S5YsMb8BdenShaysLH77zZgCtGXLlhw7dsz8J1avXj2qVKnC6tWrAahcuTItW7Zk4cKFqCoiQteuXVm7di2HDx8GoHXr1hw6dIgdO3YAxlNCenq6qcWzngDKli0b8noCaNeuHbt27WL37t0AnH/++SQkJJS6nvLy8sy3B+56WrFiBUDI6+nTTz9l1qxZpKenAwTd3dOqv8CZfz7hrocNGzaY9h0/fpwNGzaEzF+Kq4fc3Fyz7mLdX9auXeu3noryF7ePxZK/lKSeAvlLTk4OW7duLb2/BJrmw3PBaMRewWfKK4t5LwSygCTX9ljgXdf6BOAVn/TrgAF+ymkBLAISMJ5i6wKdgeX+zmt1+p4TJ05YShcMhYWFOmzYMHN6HkCrVKmip0+fDvm53IRDR7Swi5Zo6SDM011FSleTJk1M/1mzZk1EzunGLtdgIBx9ZyjKX6x+g+wB5AAfiMgOEXlWRP5gMe8lQD3gZxHZB9wPDBCR1a6Gs6U7oYikAg1d+/HYLxiddO7BaO0TVHUnRoefFhbt8It7Gp1Q8s033zB58mRzOykpicmTJ5OYmBjyc7kJh45oYRctdtHhSyR0FRYWen2DjPQrVrvWnRtHnzWsTnf1vao+oKp1gMFAZeArEfnBQvZJGI3eha7lNYwxlT0xvjM2F5EBIlIO+Bvwg6pu8CnjVuB7VV0D/AakuHrRXgpsoxS4XxOEEs8PxG3atOGnn36if//+IT+PJ+HQES3sosUuOnyJhK5p06aZr9IyMzOpUKFC2M/piV3rzo2jzxoliTW0EWNeyF+AYifdUtUTwAn3tojkALmqesC1PQDj6fBdYAVwvWd+11CQe4GOrvLyReQu4GsgFxhaAg1h5cCBM59Ve/XqRcOGDaNojYNDfJGbm8sDDzxgbv/lL3+JojUOv2esdtKpBAwAbgDaA19gDPGYEewJVXWsz/aXQJMi0h8Emvvsew94L9hz+8PdMSGUuGNHgvExPhKEQ0e0sIsWu+jwJdy6XnrpJfP1akZGBo888khYz+cPu9adG0efNaw+Qe4BvgHeB/qrEUXHFiQklDoYjxfr1q3ju+++A4weYFdddVVIyw9EqHVEE7tosYsOX8Kt6+WXXzbXx44dS+XKlcN6Pn/Yte7cOPqsYbWTTkNVvUxV34yXxtHq9D3ursmhYtq0aeb6ZZddRkZGRkjLD0SodUQTu2iJtI5ITXcVTl07d+5k586dAKSmprr1RBy7XIOBcPRZ85eAT5Ai0kVVF7k2m4qI31hPqvp1sZZEgYoVKzJp0qSIn9c9rgmgd+/eET+/w++Xvn370rdvXyZPnhz0TWy0/MUX9xg+gE6dOpGUlBRFaxzsjBV/KeoV6yuc+fb3ZoA0iisEXLxSvXr1kJW1d+9eZsw481k2kjOgh1JHtIlFLSNGjKBmzZo89thjlvP407FgwQJuvPHGuO5mH876Wbx4sbkeSf/xJRavwVDi6LNGwFesqtrcY71+gCWuG0cI3fiq3NxcunXrZm5XqVKFxo0bh6RsK0R6nFg4CUbLkiVL6NixIxUrVqRKlSp06tSJlStXFp+xCKZMmULnzp299r322mtBNY5Qsjr58ccf6dmzJ1WrVsUn6lTMEM5r7dtvvzXXO3bsGLbzFIed/MkfodB3ySWXULlyZXM4DhjRwh599FG/6Xft2sWgQYPIyMggNTWViy66iFmzZnml+eyzz7jwwgupUKECVatWpXv37mbkIF82bdrEVVddxTnnnEOVKlXo2bOnGYknVPVndcLkzwLsn+ZvfzzhDodUWiZNmuQVGmvChAlhDQzgS6h0xAJWtRw9epQ+ffpw9913c+jQIXbv3s3jjz8eM1P5lKROEhMTGThwIG++GeilTfQJ17V24sQJM+SZiNCmTZuwnMcKdvInf5RW344dO1i8eDEi4vXWLBCHDh2ic+fOJCUlkZWVxcGDBxk1ahQ33HAD//3vfwHYsmULN910E+PHj+fIkSNs376dO+64gzJl/DdT2dnZXHnllWzcuJH9+/dz0UUXmZ0iQ1Z/gULseC7A0QD7D1nJH43Fauis+fPnW0pXFPn5+VqxYkUzLNajjz5a6jKDJRQ6YgWrWlauXKkVK1YMePzNN9/UJk2aaKVKlbRHjx66Y8cO8xigr776qp533nlaqVIlveOOO7SwsFB/+uknTU5O1jJlymhqaqpZ/uDBg/WRRx4x7atZs6Y+//zzes4552j16tX13//+t1l2bm6u3nfffZqZmamZmZl62223maGv3HmLY/PmzWq4Z/AQ5lBz4brWli5davpQ06ZNw3IOq9jJn/xRWn1PPPGEduzYUUeNGqW9e/c293v6iSePPvqoNmvWTAsKCrz2jxs3TuvUqaOFhYX68ccfa8uWLUts02+//aaAHjx4MCh9RflLkU+QIvKkiDwJJLnXPZZ3gZ2haaajR0pKSqnLeOCBBzhyxPjOW6VKlYCvGMJJKHTECla1NG7cmISEBAYPHsznn39uBnUGmD59Os888wzTpk3jwIEDXHzxxfz5z3/2yj9r1ixWrlzJ2rVrmTp1KnPnzqVp06a89tprdOjQgZycHLKzs/2ee9++fRw5coTdu3fz5ptvcuedd5rnf/DBB9m0aRPvvPMOW7ZsYffu3Tz55JMl+zFikHBda57zP7ZoUaoIkqXGTv7kj9Lqe/vttxk0aBCDBg1i7ty57N+/v8j08+bNY8CAAWc9DQ4cOJCff/6ZTZs20bp1azZs2MCoUaOYP39+0DNyLFq0iOrVq5ORkRGy+ivuFWtt11LGY702UAsjks61IbEiirRr165U+V955RVeeOEFc/uaa66Jyiu+0uqIJaxqqVChAkuWLEFEGDZsGOeccw5XXnkl+/fv5/XXX+fhhx+madOmlC1bljFjxrBmzRpzCAEYE1lXqlSJOnXqcOmll7JmzRrLNiYmJvK3v/2NxMREevXqRVpaGhs3bkRVmTx5Mi+++CI9evQgPT2dMWPG8OGHHwb7M8Qs4brWohl71Rc7+ZM/SqNvyZIl7Ny5k4EDB9KmTRsaNmzI+++/X2SegwcPcu655561373v4MGDNGjQgAULFrB7924GDhxI1apVGTJkiKWGcteuXdx5553m/+JQ1V+RDaSqDlXVocCd7nXXcrOqPqyqW0JiRRiwOq7LPZ1KSTh69Cj33HOP177HH3+8xOWVhtLoiDWC0dK0aVOmTJnCrl27+PHHH9mzZw8jR45k586d3HvvvVSqVIlKlSpRpUoVVNUrRqNnT7fy5csHdceakZFB2bJnOoG78x84cIATJ07Qpk0b0tPTqVSpEldccYVX+EE37733HmlpaaSlpfGnP/3J8rkDEalxkOG61jyfICMVgSoQdvInf5RG31tvvUWPHj2oWtWYJeqGG27grbfeKjJP1apV2bt371n73fvcZbVv356pU6dy4MABFi9ezKJFi3j66acBTF9JS0vj559/Nss4cOAAPXr04I477jDfElnRV6pxkJ6o6mQAEUnHmE1DPI6VKlh4uLA6rss9P1lJmDNnDgUFBeZ2VlYWNWrUKHF5paE0OmKNkmpp0qQJQ4YM4fXXX6d27do88sgjDBo0KOhyStN7tGrVqqSkpJCVlcXmzZu55JJLAqZ1v6IKFZEaBxmOa+3gwYN8+umn5vYf/mB1sqDwYCd/8kdJ9Z08eZKpU6dSUFBg3mCeOnWK7Oxsc15Kf1x22WV88sknPP74416vWadOnUrt2rX99vj/4x//SP/+/c2OW/5uYA8fPkyPHj248sorvUISWtFnxV+s9mJtKiLfY0xmvMW1bHYtv1u2bDnzAH3LLbdwwQUXRNGa3x8bNmxg/Pjx5pjCX375hQ8++ID27dszYsQInn32WbKyjJnTjhw5wscff2yp3GrVqrFr1y5Onz4dtE1lypRh2LBhjBo1yvwmuXv3bubOnWspv6qSm5trnjs3N9erG71dmThxovlP7cILL6RDhw5RtsjBH9OnTzcnR16zZg1r1qxh/fr1XHzxxbz99tsAFBQUkJubay6nT59m1KhRHD16lFtuuYV9+/aRm5vLBx98wNNPP80//vEPRIQlS5YwefJkcwLkDRs2MGPGDNq3b+/XlqNHj9KzZ086derEuHHjwiM4UO8dzwVYALwAVAIOYzySvgLcaCV/NBarvfJyc3MtpfNHv379zF53b7zxRonLCQWl0RFrWNWya9cuvfbaa7VGjRpavnx5rVGjhg4fPlyPHDmiqqpvv/22Nm/eXNPT07VWrVo6dOhQMy+gmzdvNrc9e9+dOnVKe/XqpZUrV9aMjIyzjvvriVq3bl2dN2+eqqqePHlSH374Ya1Xr56mp6drkyZNdMKECQHzerJ9+3avibYBrVu3rqXfw0NbWHuxhvJa++GHH3TYsGGanJxs6v3ggw9CVn5JsZM/+aOk+nr27KmjR48+a/9HH32k1apV00GDBp11/Xbq1ElVVXfu3KnXX3+9Vq5cWcuXL69t27bV6dOnm2WsW7dO+/Tpo5mZmZqamqp169bVBx54IOBE81OmTFFAy5cvr6mpqeayc+fOoPQV5S9WG8jDQKJrPdv1NxXYbiV/NBarDr9+/XqLP+PZNGjQwLwIli9fXuJyQkFpdMQadtESLR3hbiBLqysvL08/+eQT7dq161n/TBs0aKB5eXmlKj8U2OUaDISj7wxF+YvVYOW5gHvU+0ERqYPxejYykbjDyL59+0qUb8eOHWzbZnx+LV++PK1atQqlWUFTUh2xiF202EWHLyXVdfDgQcaNG0eDBg0YMGAACxcu9Dp+4YUXMnPmTK/OT9HCrnXnxtFnDatX4mJgIDAF+C/wOXAKY9Li3yVfffWVuX7xxRc7QZUdHIpg5cqVdO/enWPHjnntT0hIoH///txzzz106tQpZsPrOfw+sdqLdaDH5hggC0gD3g6HUaHA3W3d3VMpECXtWPPuu++a6927dy9RGaHETh2E7KIl0jpmzpzpHqZR4mEexfkLlEzX5MmTvRrHqlWrMnz4cG6//XZq1aoVdHnhxi7XYCAcfdb8RYxXsPajbdu2umrVqmLT7d271+8A1qLIz8/3irO6efNmzjvvvKBtDCUl0RGr2EVLtHSIyHeq2jaYPFb9BUqmq23btuZE4k8//TSjR4+mXLlyQZURSexyDQbC0XeGovylqPkg38H4cF4kqnqTJStilI0bNwZ9oXz99Zk3y+ecc07UG0comY5YxS5a7KLDl2B15eXlsW7dOnP7tttui+nGEexbd24cfdYo6hVrzEbJiTaer1evu+66KFri4BD7bNmyxRzXWatWLTIy4r5vn8PvhIANpKo+EaqTiEgyxrjJy4AqGI3vGFX93HW8O/AyUAdYAQxR1Z2uYzcA4zF60g5V1QWu/Q2Bd4CLVbWAElKzZs2g0mdnZ/POO++Y26GMhFIagtURy9hFi110+BKsLs/JoRs2bBhqc8KCXevOjaPPGpY66YhIt0DHVNVKT9ayGMHNuwI/A72AqSLyByAHmAbcCswEngI+AtqLSFlgHNAaaANMBNwTOf8LGF2axhEIuoOAe74xMALtBoryEGlisaNDSbGLFrvo8CVYXZ4NZLz8Y7Zr3blx9FnD6jjIN32WGcAc4A0rmVX1uKqOVdUdqlqoqrOA7RiNXn8gS1U/VtVcYCzQUkSaYIyz3K2qe4EvgQYAInKNa/9yi/YHJJigvUeOHGHRokXm9l133VXa04cMOwVXtosWu+jwJVhd7nB/ED8NpF3rzo2jzxpWh3l4zT0jIgnAo8Ax/zmKRkSqAY0xhovcDphRblX1uIhsBZoBnwIZIlILaAVkiUia69wRH1vheScsIjz00EORNsHBIa44ffq01yeJiy++OIrWODgER4lCVqhqgYg8DezCiNFqGRFJBN4D3lLVDa4Gz3cuoCNAuqoWisjtGMEJTgHDgCeBl4A/iMjjwGngPlX90bOAvXv3cv755wOQnJzM0KFDzWg3GRkZNGvWjEWLFnHixAmWLFlC586dWb16NUePHgWMbun79+/nl19+AaBRo0Z4doNv27YthYWF5hNlcnIyHTp0YNWqVWbU+Xbt2rFr1y5ziqXzzz/fDPQLxnRL9evXZ9myZYAxiWm7du1YsWKFGbi5Q4cObN++3YwMccEFF1BQUMDGjRsB4468Vq1anDhxggULFpCWlkbbtm1ZtmyZGeS6c+fObNq0yQwC3Lx5c06dOmVOL1S7dm2qVatm6qtQoQKtW7dmyZIl5OfnA9ClSxeysrL47bffAGjZsiXHjh0zownVq1ePKlWqsHr1agAqV65My5YtWbhwIaqKiNC1a1fWrl1rBvFu3bo1hw4dYseOHQA0aNCA9PR0U4tnPQGULVvWUj0lJyebMwBkZmbSuHFjlixZEvF6SkhIYMGCBV715L6zDXU9ffrpp8yaNYv09HQwZtwJCqv+AmdmSrBSDy+99JKpp2bNmpQvX54FCxZE3V+Kq4eCggKz7mLdX9yzaATjL24fiyV/KUk9BfKXEydOsHXr1tL7S6AYdMUtwJ+APUHmKQN8CPyPM7FdJwCv+KRbBwzwk78FsAhIwPimWRfoDCz3TRtMbEmrtGvXzowZ+cADD4S8fAeHUECYY7EGQ9++fU2fGTt2bFjO4eBQGoryF6vTXf0iIj97LAeBjwHL7xjFiCH1JlDN1fjluQ5lAS090qUCDV37ffNPBO7BaPET1OjputLVcJYI911Ocdxzzz1e77W7du1a0lOGBas64gG7aLGLDl+C0bV161Zz/corrwyHOWHBrnXnxtFnDauvWG/02T4ObFLVo0Gc61WgKXCZqnrOZvkp8A8RGQDMBv4G/KCqG3zy3wp8r6prXL1bU0TkAoyhISWetNnKXHtff/01L730krmdkZFBu3btSnrKsGCnOQPtosUuOnwJRtfBgwfN9WrVqoXDnLBg17pz4+izhtVOOguLTxUYEakL3IbxHXGfR0Di21T1PVfjOBF4F2Mc5PU++asC9wIdXfbki8hdGMHSc4GhpbGvKPbt28fll19ubtepU4fFixc7g50dHIqhsLDQ/AYHOD7jEHdYisUqIhUxXm22wghSbqKqPcJjWumwGlsyPz+/yOl1brvtNiZNmmRuv/vuuzETHMCT4nTEE3bREi0d4Y7FalVXdnY2lStXBoxOFr4zecQydrkGA+HoO0NR/mJ1HOTHwCUYT2wf+SxxzaZNmwIeKywsZNasWeb2pZdeyoABAyJhVtAUpSPesIsWu+jwxaoud+9LgCpVqoTLnLBg17pz4+izhtUGsj1whapOVNU3PZeQWBFF3F2E/TF27Fj27NkDGF2E586dG7NBlovSEW/YRYtddPhiVdehQ4fMdfeTZLxg17pz4+izhtUGcglGB5u4wT2/nWu+r6D55ptveOqpp8zta6+91muKKweHWGPmzJkMHz4cSjEfZEn9xR+ePVirV68esnIdHEKBFX+x+g0yE2Ps4gpgv+cxVX2ydGaGB6vfVA4ePEjVqmePEx09ejQvvviiub13796YdvJAOuIRu2iJlo5wf4O0quvBBx/k//7v/wB46KGHePbZZ4MxKarY5RoMhKPvDKH4Bvk0UBtjDGMjjyX6EyGWkkDdgT3fYT///PMx3ThCdLpt16tXj5SUFNLS0qhevTpDhgwxo20ADBkyBBFhxowZXvlGjhyJiDBlyhTACEd23333UatWLdLS0mjdujWjRo3yex73Euo4uKrKgw8+SEZGBhkZGTzwwAMUdfP41Vdf0aRJE8qXL8+ll17Kzp07vY6vXr2aPn36kJaWRrVq1ZgwYUJI7Y0mVq81z3HDbdq0CZc5YSEc/hQuf6lfv37Q/lJafcH4y44dOxARL3s8386B4S9dunQJmb+ErP4CRRDwXDBirp5rJW2sLFYjg8yfP9/v/saNG5sRQNasWWOprGgSSEc4qVu3rs6bN09VVffu3astWrTQMWPGmMcHDx6sjRs31v79+5v78vLytEaNGtqwYUP9z3/+o6qqY8eO1S5duuju3bu1sLBQP/jgA33rrbf8nidcvPbaa9q4cWP95ZdfdNeuXdq0aVN99dVX/aY9cOCAVqhQQadOnaonT57U+++/X9u1a+d1/JxzztExY8Zobm6uHj16VH/66aew2u8JYY6kY+Va++yzz0z/AXTHjh1Bqogu4fCncPnL9u3bg/aX0uoLxl+2b9+ugObl5fk97vaXd999N2T+Eoy+ovzF6hPkNiCv2FQ2IS8vz4yZCHDeeXH/oBx2qlevTs+ePVmzZo3X/r59+7J06VKzR+OcOXNo0aKF1xP5ypUrufrqq6lRowYiQvXq1bnpppsiaT5vvfWWeVdes2ZN7rvvPvOO3Zdp06bRrFkzrr32WsqVK8fYsWNZu3YtGzYYsS1eeOEFevbsyeWXX05ycjLp6ek0ber/E7777vo///kPtWvXpnLlyrz22musXLmSFi1aUKlSpZiaNcYKv/76K7feequ5fd1111G3bt0oWhR7hNJf6tWrF9P+Uhxufxk0aFDM+YvVBvIdYIaI/FlEunkuIbMkStSuXfusfStWrDCDDtepU4fU1NRImxU0/nREkl27dvH555+fdTNRrlw5rrzySj788EMA3n777bOcuX379rzwwgu88sorrFu3rlRzub3//vtUqlQp4PLzzz/7zZeVlUXLlmbEQ1q2bOk1TVNRaVNTU2nYsKGZfvny5VSpUoXRo0eTmZlJ3759A57XzYoVK9i8eTMfffQRI0eO5Omnn+bLL78kKyuLqVOnsnBhqWJ1hJTirrURI0Zw4IAx/8C5557LK6+8EgmzQkq4/SmU/qIW+pH44tYXCX9xU7duXWrVqsXQoUO9Iiy5/aVjx44h8xf36INSE+jR0nPBmLvR37LNSv5oLFZfGR07duysfffff7/5amjYsGGWyok2/nSEm7p162pqaqqmpaUpoN26ddPDhw+bxwcPHqyPPPKILl68WNu3b6/Z2dmamZmpJ06c0E6dOpmvjPLz83XixInasWNHTUpK0urVq+uUKVPOOk/FihXNZdKkSSHVUqZMGV2/fr25vWnTJgW0sLDwrLQ333yzPvjgg177OnbsaOpp1KiRVqxYURcsWKAnT57Uu+++Wzt27Oj3vO7XT7t27TL3ValSRT/88ENzu3///vriiy9a1kKYX7EWda3Nnz/f69XqnDlzLJcbS4TDn8LlL+eee27Q/lJafcH4y7Fjx3TlypWal5en+/bt0wEDBmiPHj3M425/+fbbb0PmL+PGjbOspSh/sfQEqar1AywNQtNMhx6r3dZ9e+6pKp999pm5fdVVV4XFvlBjtQdiqJk+fTrHjh1jwYIFbNiwwevO0E3nzp05cOAAf//73+nTpw8pKSlexxMSErjzzjtZunQp2dnZDBw4kJtvvpn169d7nSc7O9tchg0bFlIdaWlp5pRAAEePHiUtLQ2PsIgB07rTu6bOISUlhauvvhpVpVy5cjz++ON88803HDlyJOD5PeOUpqSknLXt2ZkjEJEa5hHoWlNVHnnkEXN78ODB9OzZM1hTYoJw+VM4/OWRRx4J2l9Kqy9Yf2nbti1ly5alWrVqTJw4kS+++MLM7/aXP/7xjyHzFyuBAqz4i9VXrHFHxYoVmTRpEn379g0q37p168x531JTU+nePeLzMsclXbt2ZciQIdx///1+j994442MHz++2G8lbmepXLmyOb9cMLz33nteveV8l0Cvbpo1a2bOqwewdu1amjVrZint8ePH2bp1q5m+RYsWXv8o3OtagldhwdC3b193WMTA/1kCUFJ/8eR///sf33zzDQBJSUk88cQTJS7L7oTSX+68886Y9hdffP0hlv2lpNNdmUuojI0WFSpU8Np+/vnnzfWePXvGbOQcX3x1RIORI0cyb968szoegDFd2Lx58+jSpctZx/75z3+yYMECTp48SX5+PgsXLuTYsWPmhL3BMGjQIHJycgIuderU8Zvvpptu4oUXXmD37t3s2bOH8ePHM2TIEL9pr776an788Uc++eQTcnNzefLJJ2nRogVNmjQBYOjQoXz66afs3buXvLw8nnrqKTp37kylSpWC1hOLBLrWnnnmGXN9xIgRcd0xJxL+FCp/eeutt4L2F7e+SPjLihUr2Lhxoxm8/p577uGSSy6hYkXjwc3tL2vWrAmZvyQnJ5c4rydWnyBvBP7isTwA7AbGh8SKKNK6dWtzfevWrbz//vvm9g033BANk0qEp45occ4553DTTTedNcYJjFic3bt39/sKJiUlhfvuu4/q1atTtWpVZs+ezSeffEKDBmfe4Pft29frzvbqq68Oqe233XYbffv25Q9/+APNmzend+/e3HbbbebxZs2a8d5775k6P/nkEx555BEqV67MihUrzE4VAN26deOZZ57h/vvvJzMzky1btnhdV/GOv2tt9erV5tNjYmIiDz/8cKTNCimR8KdQ+cvLL78ctL+UVl8w/rJt2zauuOIK0tPTad68OcnJyXzwwQdmWre/9O7dO2T+UqNGjVLlNwn0cbK4BagOrClp/nAvVjsdLF682Fx/7rnnzM4FLVq00NOnT1sqIxbw1BHv2EVLtHQQ5k46/nQNGTLE9J0bb7yxtBKijl2uwUA4+s5QlL+U5hvkKaB+KdvnqOMezgEwb948c/3uu++Oq9irnjriHbtosYsOX3x1ffzxx7z99tvmdryN2/SHXevOjaPPGpYmzBIR33ir5YFewOchsSIG2Lx5M19++aW53a1b3A/xdHAIO7Nnz+aGG26gsLAQMDqftGvXLspWOTiEBqvByv/js+s4sAZ4R1UjHwTUAlaDLxcWFqKq9OnThzlz5gDG+/nvvvsu3CaGlMLCQsqUsUenZLtoiZaOcAcrd+v6+uuv6dWrlxn3skmTJixcuJDMzMzgjY4x7HINBsLRd4ZSBytX1aE+y12q+kasNo5gfVxXVlYWzz//vNk4AowZMybc5oWc4qJYxBN20RJpHZEaB5mVlcXWrVvp16+f2Tg2aNCAL7/80haNI9jnGgyEo8+ivwT6OOl6suwEPBfg2DigfVH5o7lY7XQwefJkTU5ONjsYXHjhhX6jQcQ60QhWHi7soiVaOghzJ525c+dqmzZtTJ+pVauWbt++PbQiooxdrsFAOPrOUJS/FPcEOQZYFODYAuCRAMfihldeecW8C65ZsyazZ8/227XawcHBYNKkSeYniMTERKZPn069evWia5SDQxgoroG8EJgT4NiXQHxN8uZDQUEBP/zwg7k9d+7c0I2fiTCegYPjHbtosYsOT7755hs++eQTc/v555+Pu7kerWDHuvPE0WeN4hrICkBSgGOJQHoojBCRKiLyqYgcF5GdInKDa39tEVkuIodEZLxPnjkiElRHBF+WL19OQUEBAJmZmZZDJcUir7/+erRNCBl20WIXHZ5Mnz7dXO/Tpw9333139IwJI3asO08cfdYoroHcAPQIcKyH63goeBk4DVQDBgGvikgz4GHgLYzxlv3cDaKIXIcxk0iJI+7m5eUxevRoc7tDhw5FprfSeSGa6f7zH9+OxpE5bzjKjJYWu+gIJ4sXLzbXb731VkufI8JxLYUzLVivu3DaEazNwaSPBX3hLDsYfUVRXAP5IvC6iPQXkTIAIlJGRPoDrwEvlNYAEUkFBgCPqWqOqi4BZmCEtKsPfK2qR4CVQAMRqQA8hPF9tMR88cUXfPvtt4ARHX/s2LFFpo/1f8ZWiYcGMlrntYuOcHHixAmv4U+dOnWylC8WGpBw/naxYnO4NMaKzdG4/osdBykio4EngGTgIFAVyAUeV9UXS22ASCvgG1VN8dh3P9AV4wl1J/AusAroBwwDVqvqW8WUewzvG4ADLvt9qRpgvy8VsTZLQrTSRUtHOMp06iT4dFWBc1zrhaoa1OePIPzFfS4rutyE41oKZ1oITmOs2BxM+ljQF86yi9NnyV+KjaSjqi+IyBtAByAD+A1YpqpHi85pmTTOFn0E4/vms8CrGI3iyxjfPVsAT4jI+0AtYKqqTvRjd0i+jzo4/B5w/MXB4WwshZpzNYZzw2RDDkZnIE8qAMdU9RBwHRivdjGGnIzAeMX6IzAEWC0iX6tq8JOhOTg4ODg4BCAWYg1tAsqKSCOPfS0B31AIw4Hlqvoj8AeMwZ2ngXVA84hY6uDg4ODwuyHqDaSqHgemAU+KSKqIdAKuAt5xpxGRTOBOYKxr13bgUhFJA9oC2yJqtIODg4OD7Yl6A+niDiAF+BX4ALhdVT2fIJ8HnlTVHNf2s0A34Bdghudwj0BjKv0hIqNEZJ+IHBGRf4tIcknKCRch1LJARHJFJMe1bIyMAvP8lnSISHMRmSsiB0XkrN5j0a6TEOqIan1YIVTXXqwSRF0OEZECj7rKEZFLImtt8IjIXSKySkROiciUYtLGY/1Z0lfq+gsUgy5eF4wG9iOMzj+dMTr8NPOTriewH2gGVMYInTcu2HLiRMsC4NY4qJPzgVsw3iBoScuJAx1RrY9IXnuxugShbwiwJNr2lkBff4xe/68CU4pIF6/1Z1Vfqeov6kJD/KOlYgQcaOyx7x1/FQ68Dzzjsd0d2BdsObGuxbUdtX/IJfktgfN8G5Zo10modES7PkKttbhrLxaXIPXFZQPpYf/fi2lA4q7+gtRXqvqLlVesoaIxUKCqmzz2rcW4O/KlmeuYZ7pqIpIRZDnhIlRa3DzreuW3NMKviEL1W0a7TkJ9/mjVhxVCfe3FGsHWZStXXW0SkcdExFLv/zghHusvWEpcf3ZrIIsaU1lcWvd6epDlhItQaQF4EGgA1AQmATNFpGHoTC2SUP2W0a6TUJ4/mvVhhVBee7FIMPoWYfSSz8SI+PVn4K9htS6yxGP9BUOp6s9uDWTAMZUW0rrXjwVZTrgIlRZUdYWqHlPVU2pEIFoK9AqxvYEI1W8Z7ToJ2fmjXB9WCNm1F6NY1qeq21R1u6oWquo64EngmgjYGCnisf4sU9r6s1sDaXVMJa59LX3S7VfV34IsJ1yESos/FIjUpJeh+i2jXSfhPH8k68MK4bz2YoHS1GWs1VVpicf6Kw3B1V+0P7KG4aPthxg91FKBTgTunXYFsA+4AKP31td49/y0VE6sawEqYfRUK4cROWkQcBw4PwZ1iMvOC1wXcjkgOVbqJBQ6YqE+InXtxfIShL4/AdVc600wIng9Hm37Legr67rGnsXogFQOKGuj+rOqr1T1F3WhYfjhqgDTXf90fgZucO2vg/E6oY5H2tEYXZyPAv/x+Wfst5x404IRkHclxiuTbGA5cHks6gDqYTQonsuOWKmTUOiIhfqI1LUXy0sQdfm8S9txjIAkTwKJ0bbfgr6xfq7BsTaqP0v6Slt/xc7m4eDg4ODg8HvEbt8gHRwcHBwcQoLTQDo4ODg4OPjBaSAdHBwcHBz84DSQDg4ODg4OfnAaSAcHBwcHBz84DaSDg4ODg4MfnAbyd4KI1HHNhZYQZTsujub8hyIyRkTeCEO5Y0Ukz/Ubp5awjK9d80QuCbV9DsHh+It5/t+1vzgNpB9EZIeInHRV3j4RmSIiadG2qzSo6s+qmqaqBVG2Y7Gqnl9cOpcDvVuac4nIJSKyy+f8z6jqraUptwg+cv3Gx0uSWVW7ASNCbFPYcfwlrHY4/hKASPiL00AGpq+qpgEXAq2Ah6NrTmwQK1P9iIFz/cYOjr/4wfGX+Mb5wYpBVfcBczEcHwAReUhEtorIMRH5SUSu9jg2RESWiMjzInJYRLaLyJ88jtcXkUWuvF+KyMued34i0l5EvhGRbBFZW9Rcga4797+KyA8iclxE3hSRaiLyuUf5lV1p64mIuh3WZec2V7rtIjKoiPOoiNwpIpuBza59fURkjcvOb0SkhUf61iLyvavsj0XkIxH5u+uY112qiDwoIrtdaTeKSHcRuQIYA1zneipZ60q7QESeFpGlwAmggYgMFZH1rvzbROQ2V9pU4HOghquMHBGp4XunLSJXikiWS8cCEWnq8/ve7/p9j7h0lAv0O/n53aaIyCuu+sgRY+7H6iLyT9e1sUFEWlktLx5w/MXxF1v5S7Rj6sXiAuwALnOt1wLWARM8jl8L1MC4wbgOI87fua5jQ4A8YBiQANwO7AEzrN8yjPiASUBnjPiH77qO1QR+w5j6qAxwuWv7nCLsXA5Uc+X9FViNcQefjBF4+HFX2noY8QrLYgRoPoorQDZwLkUE/Xblm4cRvzIFaO06VzuXxsEuW5JdunYC9wKJQH+M2dv/7irrEmCXa/184BeghoeNDfVMrMV3fexYgBE3s5lLRyLQG2iIESC8K8Y/gta+5/IowywXY+Lc467fORF4ANgCJHn8vt+66roKsB4YEeA38mfvFOAg0AYjmPLXwHbgJtfv9ndgvk+eIcTZDPY4/uL4i039xXmCDMx0ETmGcUH+CjzuPqCqH6vqHjXmGPsI4y7xIo+8O1V1shrfL97CcKhqIlIH+CPwN1U9rapLgBke+W4E/qeq/3OVPQ9YRdFzBb6kqvtVdTewGFihqt+r6ingUwzn90ch0FxEUlR1r6oWN9XPs6p6SFVPYvwze12NeQ0L1JjT8BTQ3rWUBf6lqnmqOg3DafxRgPFP4gIRSVTVHaq6tRg7pqhqlqrmu8qfrapb1WAh8AVwcTFluLkOmK2q81Q1D+MfcQrQ0SPNv1x1fQiYiceTkUU+VdXvVDUXoz5yVfVt17XxEYHrJ95w/MUbx19s4C9OAxmYfqqajnFX1QSo6j4gIjd5vC7JxpixuqpH3n3uFVU94VpNw7izOuSxD4x/KG7qAte6y3WV3RnjH0Yg9nusn/SzfVZnCTU+il+H8YF7r4jMFpEmLm1ZHq9YPB3H1877fOys7dJXA9itrts7P3k97dgCjMS4m/xVRD4UkRpFaD2rLBH5k4gsF5FDLjt64V0XRVED4+7dbU+hq/yaHmn2eayfwM/vWQxB10+c4viL4y9gM39xGshicN1lTcG4W0JE6gKTgbuADFWthDHHmJVJOPcCVUSkvMe+2h7rvwDvqGoljyVVVceVXok3qjpXVS/H+GeyAUMTqtpMjZ5laaq62DOLj51P+9hZXlU/cGmsKSKev4enRl873lfVzhj/RBR4zs/5vLK4V0QkGfgEo26querif5ypi+KmqtnjOq+7PHHZuruYfA4BcPzlTBYfOx1/iUOcBtIa/wQuF5ELMb5HKHAAQESGYtwRF4uq7sR4BTRWRJJEpAPQ1yPJu0BfEekpIgkiUs71kb5W6KSAGB0TrnR9mD+FMX9aMN3ZJwMjRKSdGKSKSG8RScf4ZlQA3CUiZUXkKrxfp3nacb6IdHM5bi7GHaLbjv1APSm6510SxiunA0C+GJ07engc3w9kiEjFAPmnAr3F6OiQCNyH8Xt8U/xP4FAE/8TxF08cf4lTnAbSAqp6AHgbeExVfwLGY1zY+4E/AEuDKG4Q0AGjM8HfMd6rn3Kd5xfgKoweaQcw7jz/SujrqQzGxb0HOITxsf4Oq5lVdRXGd5WJwGGMD/VDXMdOY3Q0uAVjQuAbgVm4NPqQDIzD+DC/D8jE0A7wsevvbyKyOoAdx4B7MBz3MHADHt+oVHUDxqzx21yvtmr45N/osu8llw19MYYrnLb0Qzj4xfEXbxx/iV+cCZOjjIh8BGxQ1ceLTRyniMgK4DVV/U+0bQkXIvIoxti/PKCmlmDws4jMw+i08a2qdg+xibbA8Rd7EC/+4jSQEUZE/ohxF7od4/XGdKCDqn4fTbtCiYh0BTZi3GUOAl4DGqjq3qga5hB3OP7iEE1iIsrD74zqwDQgA9gF3G4nZ3dxPsZrnDRgK3CN4+wOJcTxF4eo4TxBOjg4ODg4+MHppOPg4ODg4OAHp4F0cHBwcHDwg9NAOjg4ODg4+MFpIB0cHBwcHPzgNJAODg4ODg5+cBpIBwcHBwcHP/w/zoFoKz1QgLEAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, axs = plt.subplots(nrows=1, ncols=2, figsize=[6.5, 3], sharey=True)\n", "labels = ['w/o ancillary refine', 'w/ ancillary refine']\n", "#axs[0].plot(ts_pre_S1, p_S1, label=labels[0], lw=2)\n", "#axs[1].plot(ts_pre_A2, p_A2, label=labels[0], lw=2)\n", "#axs[0].hist(ts_S1_list[-1], bins=20, range=(0, 0.2), density=True)\n", "#axs[1].hist(ts_A2_list[-1], bins=10, range=(0, 2.1), density=True)\n", "axs[0].plot(np.sort(np.abs(ts_S1_list[-1])), p_S1, label=labels[1], lw=3, color='k')\n", "axs[1].plot(np.sort(np.abs(ts_A2_list[-1])), p_A2, label=labels[1], lw=3, color='k')\n", "\n", "# axis format\n", "axs[0].set_xlim(-0.01, 0.18)\n", "axs[1].set_xlim(-0.10, 1.6)\n", "for ax in axs:\n", " ax.tick_params(which='both', direction='in', top=True, bottom=True, left=True, right=True)\n", " ax.yaxis.set_major_formatter(ticker.PercentFormatter(xmax=100))\n", " ax.set_ylim(0, 100)\n", " ax.xaxis.set_minor_locator(ticker.AutoMinorLocator())\n", " ax.yaxis.set_minor_locator(ticker.AutoMinorLocator())\n", " ax.grid('on', linestyle='--')\n", " ax.set_xlabel('Range mis-registration [m]')\n", "axs[0].set_ylabel('Cumulative distribution', labelpad=-4, loc='bottom')\n", "for ax, num in zip(axs, ['(a)', '(b)']): ax.annotate(num, xy=(0.05, 0.87), xycoords='axes fraction', ha='left')\n", "std_S1 = ut.root_mean_sq_error(ts_S1_list[-1]) #np.std(ts_S1_list[-1])\n", "std_A2 = ut.root_mean_sq_error(ts_A2_list[-1]) #np.std(ts_A2_list[-1])\n", "for ax, title in zip(axs, ['Sentinel-1\\n\\nRMSE = {:.2f} m'.format(std_S1), 'ALOS-2\\n\\nRMSE = {:.2f} m'.format(std_A2)]):\n", " ax.annotate(title, xy=(0.95, 0.08), xycoords='axes fraction', ha='right')\n", "\n", "fig.tight_layout()\n", "# output\n", "out_fig = os.path.join(work_dir, 'stats_cdf.pdf')\n", "print('save figure to file', out_fig)\n", "plt.savefig(out_fig, bbox_inches='tight', transparent=True, dpi=300)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "031f438f", "metadata": {}, "source": [ "### Implication for Stack Coregistration" ] }, { "cell_type": "code", "execution_count": 6, "id": "193ce687", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "S1 (304 images) - rg: 0.06 m (1-sigma) / 0.05 pixel (2-sigma)\n", "S1 (304 images) - az: 0.36 m (1-sigma) / 0.05 pixel (2-sigma)\n" ] } ], "source": [ "# Conclusion for stack coregistration - S1\n", "s1_ts = ts_S1_list[-1]\n", "s1_num = len(s1_ts)\n", "s1_rg_std = np.std(s1_ts) # Yunjun et al. [in prep]\n", "s1_az_std = 0.36 # Gisinger et al. (2021)\n", "s1_rg_std_pix = s1_rg_std * 2 / sensor.SENSOR_DICT['sen']['range_pixel_size']\n", "s1_az_std_pix = s1_az_std * 2 / sensor.SENSOR_DICT['sen']['azimuth_pixel_size']\n", "print(f'S1 ({s1_num} images) - rg: {s1_rg_std:.2f} m (1-sigma) / {s1_rg_std_pix:.2f} pixel (2-sigma)')\n", "print(f'S1 ({s1_num} images) - az: {s1_az_std:.2f} m (1-sigma) / {s1_az_std_pix:.2f} pixel (2-sigma)')" ] }, { "cell_type": "code", "execution_count": 7, "id": "717bebe1", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ALOS-2 (49 images) - rg: 0.55 m (1-sigma)\n", "ALOS - rg: 0.06 pixel (1-sigma)\n", "ALOS2 - rg: 0.38 pixel (1-sigma)\n", "NISAR - rg: 0.09 pixel (1-sigma)\n", "\n", "ALOS-2 (49 images) - rg: 0.31 m (1-sigma)\n", "ALOS - rg: 0.03 pixel (1-sigma)\n", "ALOS2 - rg: 0.22 pixel (1-sigma)\n", "NISAR - rg: 0.05 pixel (1-sigma)\n", "\n" ] } ], "source": [ "# Conclusion for stack coregistration - ALOS-2\n", "a2_sub = 0.311 # from Table III\n", "\n", "a2_ts = ts_A2_list[-1]\n", "a2_num = len(a2_ts)\n", "a2_rg_std = np.std(a2_ts) # Yunjun et al. [in prep]\n", "for rg_std in [a2_rg_std, a2_sub]:\n", " rg_std_pix_a1 = rg_std / sensor.SENSOR_DICT['alos']['range_pixel_size']['stripmap_FBD']\n", " rg_std_pix_a2 = rg_std / sensor.SENSOR_DICT['alos2']['range_pixel_size']['stripmap_ultrafine']\n", " rg_std_pix_ni = rg_std / sensor.SENSOR_DICT['ni']['range_pixel_size']['24MHz']\n", " print(f'ALOS-2 ({a2_num} images) - rg: {rg_std:.2f} m (1-sigma)')\n", " print('ALOS - rg: {:.2f} pixel (1-sigma)'.format(rg_std_pix_a1))\n", " print('ALOS2 - rg: {:.2f} pixel (1-sigma)'.format(rg_std_pix_a2))\n", " print('NISAR - rg: {:.2f} pixel (1-sigma)'.format(rg_std_pix_ni))\n", " print('')" ] }, { "cell_type": "markdown", "id": "8281722b", "metadata": {}, "source": [ "### Backup - an alternative to the boxenplot to summarize geolocation bias correction" ] }, { "cell_type": "code", "execution_count": 17, "id": "fcf86723", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAA1QAAAFcCAYAAADPrf8ZAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAACTBUlEQVR4nOzdeXxU1fn48c+ZJfu+QAhbWATCElAjq0jFryibUrUiIhCrIq1FaQW1thb52ZaqYGvrvgJWUQSLooUKxVAQZFGCJhL2LWHNnsk+M+f3x0wmmSyQZbKQPO/Xa165Z7nnPjfX4H3m3nuu0lojhBBCCCGEEKL+DC0dgBBCCCGEEEJcriShEkIIIYQQQogGkoRKCCGEEEIIIRpIEiohhBBCCCGEaCBJqIQQQgghhBCigSShEkIIIYQQQogGMrV0ADWJiIjQMTExLR2GEEIIIYQQQvDtt99maK0ja2prlQlVTEwMe/bsaekwhBBCCCGEEAKl1Ina2uSWPyGEEEIIIYRooDolVEqpXyml9iilSpRSy6q03aCUSlVKFSqlvlJKda/UppRSzyqlMp2f55RSysP7IIQQQgghhBAtoq63/J0G/gjcBPiWVyqlIoBPgPuBdcAzwEfAcGeX2cAUYDCggY3AUeC1xocuhBCXOa3R1lIoLUKXFqGLC53LxdiLC9ClhdhKitGlRdhLi7CVFjvaSkvQ1hLsJcXoTldhN/tjLytD26zOnza0s6ytVoxTxqMNCjt20GDXduzY0Vq7LWut3eu1HY37cnmbW73WFeUa1q+6jrfRm8m9Jrf0b18IIerM8W+k42TW8W8e2J11jnac/96599Ea57+7ulq9xtFZ17C+3a4vPa6z3l65T6XlmtaviLnymDXHZa9l21Tb3kXGrbJ+eUzl/Wsd19nvzvgumIyt/4a6OiVUWutPAJRS8UCXSk23ASla64+d7U8DGUqpflrrVGAWsFRrneZsXwo8gCRUQrR7lU+6rdqKzW7Dpm1Y7VZs2kaxtZhdZ3dxLPcYdrsNZbVBWRmUWjGUOctWK9jtYLWBtQxVXAglBaiSIlRpEZTZUHY7ymYHu0bZNMr1015pGVebobxs1xjsjuXyOqO9vK1i2WADg3b+tGsMdjDawVDpY6xUV/Wnyd7Y3+Rnl+wx3fYKZebWc3NAuE+4JFSiQex2jU1r7Fpjt1NpWWOza+zOEzW7dpbtjrLN2ceucfar1EdXqi/v4+pfQx/XulRs29leEccl+lSK217jSbCudkJqr3yyWsOJZ+WT80uerNZ0olzDCXvVk2CN43daed3qJ8GV1q9hP2odV1eNtfb1cdtH998DNcRit9dhXGd7TYmPaDk/vbIzJmNLR3FpjZ2UYgCwr7ygtS5QSh1x1qdWbXcuD6hpIKXUbBxXtOjWrVsjwxJC1JXVbqXYWkyRtYiikgKKSgsoKS6gqNRCYVE+WZYLZBdmkFOYSW5hFrbSUuw2qyOJsVnRVhvaZkXZ7GirI9kxldgwllkxltowldowldoxldkwl9kxlWnMpXbMVo1XGZhsYLZpzDYwWx0fRx30sEEfq2NZNJzRDmUtHUQlzlOby46u8aS6lpP2i5yQ2+yOEzibLl92TwYafNJeeZzybdeaSFQZt8a43cfV5fHriphtzpPgWvfLOU5F/0rjVOrjlmjUmvi09H8BQojmdrn8/6KxCVUAcKFKXS4QWKk9t0pbgFJK6fKvJpy01m8AbwDEx8dfHr89IZqJttvRJSXokhLsJaXostJK5RJ0aRm61L1sKy7Ekp9FwcmjlB4/AWfOYS0uwlZWCjY7Rpt2u6JS+UqJ2fkJAqJaaJ/bE7sCqwHsBrBV/hjBppz1xvJ25epnNyjsRgM2owG7EexGhd2gsBkU2rlsNyqig8KxeZtQKJRSGDA4fioDBmVwqzeoijaFqihXWsetb/m6lcYq71O+XLneoAz4GH0pLrNhtWtsNsdJtNVux2bXWG2OE2ir84Tbaqs4EbfZ7W7l8vWt9srr2Kus4xgjr7iMs7nFnMkt5mxuMYVlVscVjhpO2mu7CiLfVAshlAIFGJRyLjsqlKtNufo4/h0EVOX+FfXK2aiqrG9Qjj6u7V1q3Crbpso2Ko+LW/8q/Srtiyvei43rVq6+H5XXp1J/Qw374b5/FWXDZTL1QmMTKguOc67KgoD8WtqDAEvVZEoI4c5eWEjpqVMUff89hd/sxPLVV9gLCxs0lgK8ncveF+vYytkNGm0E7frpWMbgWMYA2mgEkwnMXo6PyYQ2GlAGA5gMaIMBZTKC0QBGIxiNKKMBTOXLjvWV0YgymVAmx0+M5ctmZ73jp8Hs5Wo3eHlhMJlRXl4YjGaU2cvZbsbg5Y0ye2Mwe2E0ezv6mH0wePuizN4YjabqiYoH/ydyv8dGqp/iMpsrkflPylm+T8vlwNl8zueX8Nq/NrRQVOJyVn6CZXSevBkNVZYNjgTfqMqXq/cxKMfHsQwGZ7tBKQyGWvq4lhUGZ52x8nK1dSrVV+tfKTbnOhc9WVXVT3qrnnhWPllVgMFwsZPgKutXPiF21dUybqWT5uonwXUY1y0m95PuyokBVcetoX9t+1Q1rmqJz6XGrW3bl8mJvWgZjU2oUnA8JwWAUsof6OWsL28fDOxylgdXahOiXbParZTaSimzl2HZuZOijZuxpuxHnzkH2bktEpNdORIXx8+KREUp7fgfnsHxUxm0639MGDSqUh9l1I6PSWM0aAwmjcGoMRjtGI0ao1Fjcv40mhTKbHQmK45kR5lNGMrrzCbHx8uMMprAUP4xg8EIXgHQ+wbHJ6gzGM0t8ntrbbTWFJXZsBRbyS+xkl9sxVJsxVJSRom14kqQ1e64MlRmc1zZKbOVXxGyU+a8QlRms1f0tdmd61Radv10bDOvuIy8Iit5xWWUWhv9gFirU/PJc6WT8lpOniuflCulHDm9Kl+uLSmg1gShYvu4TvKNqqZYKvWpEm95YlJzUlBlnUr75YrZQPXY3Pax+n6V/x5qGtc9pipJTKXfhRBCtDZ1SqiUUiZnXyNgVEr5AFbgX8DzSqnbgS+APwDfOyekAFgB/EYp9W8cz/s9CvzDs7sgROtitVs5df4wJ0/v53TGUc5lniAj8xS5OedRxSWYiq2YSqz4lGoi8jQ37q3bBdsyI5SYodQEZSZHufxnqUlhdZYrtxu9fSjzV5QFFqMDSgmnlAGlJXS1WfFVdnyVHS/AoLQzcWrEjisDhHSDDgMgegj0HQ++oZWSIKMzEapUlpOjOrPbNdmFpZzPL3F88oq5YCnhfF4JF/JLOJ9fzIX8ErIKSrGUWFv18yYGBWajAZNBuU60jQb3sntbRdlkMOC44Giosd29bMDo7OvnZaRTiC+dgnyICvYhyMdc+9WIqlcsXMvy36sQQojq6nqF6vfAwkrle4BFWuunncnUS8A/gZ3AXZX6vQ70BH5wlt9y1glx2dBaYz1/gbL0NOwWC/aCAtfHVlBAcW4WmVnp5GSkYzt9hoALBQQVajoAHRq4TasBMoLgRAfF4WhFUk/FiQ7UmIAEeQXRJ7QP0QHRdPLvRGffDvTWRnplpeG38Q8VHcsndjA6P7VSjsQovLfjE9YTfEPA7Of4ePnVvGzylgSpEcpsdudVJCtZBaVsTj1PyulcZ/JUQoalBGtrzpJqYDIogn3NBPqY6NMxkMmDo+kfHURMuD9GSU6EEEK0Eao1Ps4UHx+v9+zZ09JhiHbMeuECeRs3kvfFvyk5cAC7xdLk27QZYMXdHTnXNYDiEF/MJi+8lAEvFGYUPijCDF6EYyIcA+E2TRet6Flmw1CSB8W5UHABMg+D3Vr7hsz+4B0A3oGOW+bKf0b0hoF3QGQ/MPs0+f62N1prMgtKSc8uYsfRTDb+eI5MSwkW5y15JR6+Pc7HbCDA25HMBPqYCPA24e9twttkcF0dMhmV60qP2agwldcbDM62SnVGhdnZ12RUmI0V6xkNBswGhbfZQJCPmSBfM0E+ZnzMnn0WTAghhGgpSqlvtdbxNbU19hkqIdoUbbNx6sE5FGzb1qhxrEYo9jeDrw8GPz/M/kH4BIZg9g/A6OePKSAAk58/xogIvGJi8OreHa9u3Xjuh4/gx88g9RsoyvbQXgE3/AHi7wPvIMfTyqJRtNakZRdRVGajuMxGUanNuWynuLzO+SkutXHwnIWthy5QUNq4+d8DfUx0CPSmQ6APHYK86RDoTWR5OdCbDkHehPt7E+BjwnwZvAhRCCGEaAskoRLtXkFxPudPppLx178R8NV3NfYp9IbToWDxVRR7QZEXFHk7fpZ6GwkM6UCHyBi69hpC3wHXEd1joGPGuLrIOASntsOXs+DMPs/sVGgPiBoEUXEQcy10Gy634zVQidVGWnYRJ7MKScsq5GhGARuSz3Imt9ij2zEaFAHeJtenQ5A34wd2om9UAB0CfYgM9MbHfBm83VAIIYRoZyShEu2ONTOTvWvf4sSGT+h4Mp/QPI1RO16aVtXbNxrY1VeRHQAohb/Zn76hfYkNj6VfWD9iw2LpGdwTc31ml8s7DQf+DZlHIPkTsJytpaMCnyDwDgafYMeyT7DjKpPPRerCezvqhRu7XXMmr5ijFyycyCwkr7iMolIbha6PlcJSm7POsZxTWMa5/GKPvX8o0NtE51BfuoT6MahzMOMGdCQ8wItAb7k9TgghhLhcSUIl2qR0SzoHsw5yIu8EaZY0ckpyUKfOctW/jxC3N4cADQMusv66oYrPbwiiY1hXrg7oQpfALsSGxXJj9xvrlzwBWEvh8EY4uAHSv4PzP4K+yPMyBhP8/EuIvrLd355nt2syLCVkFZZSXGanpMxGsbXmnyVWx+125T9zCsvIKigls6CUTEsJ2YWllNk8+8xozwh//LyN+JiM+HoZ8TE7Pr5mA77minKgj4nRV0TSp2OAJE1CCCFEGyMJlWhz/v7d33nzhzddZaNNM+ffdsYk134yne9vwBLqg7V7Jzo++Vse7DKQBV5BjTv5Lc6D7z+C7f+AnBOX7j90NnS5BnrdAP7hDd/uZcjmnBI801JKyulcNu0/R+rZfNKzizw+WUN9GRR0Cvala5gv3cL86BbmR9cwP667IpJQf68WjU0IIYQQLU8SKnHZyy3JJd2Szom8E7y+73WO5B7Bt1gz/ICmT7rm6kOakEL3dQ5Gw7f9zMROuJs7b3gEo6+v5wLKOwOpn8OG34K9rIYOyvFcU6/rHc86hXaHjoPA1DZPzivPbpeeU8TpnCLSnMvp2UWcyysmq7DUY7fVXUyon5mekQH0iPAn3N8LXy8jfl5GfL1M+Fda9nMuB3qbiQr2wcvUvq8UCiGEEKJ2klCJy9ri7X/kf9s/JDxfE2KBgQUw4YJm6AGNdw0zh+uBffCdncD1P7mByeYADMqDJ8pnvodVMyD7ePU231C4+l7ofYNjoojL/BmnUqud/OIy8oqt5BWVkVdcRm5RGXlFVvKKy1x16dlFfJ+WS2ZBaYO3FeJnJjLAG18vI94mAz5mx09vs+NWO2+zodpPb5Nj+u7wAC/C/b0JC/Ai3N9LJnUQQgghhMdJQiUuO/bSUtLXrebEu68x4dgFptRlJmqzmfD77yPy4Yc9+wyL1nBoI+x52/GMVE1+8iQMm+1Iqi5Dh89beGHjAQ6dsziTJStFZY2b/rsmIX5mwv0dCVC3cD8mDIoiPiaMIJ96PrMmhBBCCNGMJKESlwV7YSHZH3xAzv8SKU7ah6HUyqWeMjJFdyJ02jR8Bw3Cp39/jEEevCqUfw7OJcOuN+Hg+pr7XH0vDPoZxIzy3HabgdVmp6jMRkGJjf+mnuMv/04lv+QiLwquA38vI13D/Ogc4kvnUF/Xz+gQx3KYv5e8N0kIIYQQlyVJqESrt//A1xTf/xt8LuQBUPW0OyvYSHCvfoRG98AUEYEpIhxz584EjB2Lwdvbs8GcS4Ev5sPJ7TW3950IcT+D2FvA0LpuLzuVVcgriYc5k1tMUaljdrzi0kovoHV+6jITnkFBkK+ZYF8zQT5mgnxNjp+Vl33NhPp7MSA6iB7h/hgMMrudEEIIIdoeSahEq6W15oVvXyBg8dtcd8H9JP90KKT0MhE5M4EpY3+Jr8mDk0rUxFYG3y2HLx6tud0/EqavhughTRtHA5Va7cz557eknM5r0PqdQ3x5/o44YiL8CfI14+9llOm/hRBCCCGQhEq0UkdzjrLs7w8Qu+MMQ45VJFPnQuDDhwcS2+9abut9G12DujZPQOsegaT33es69HdMc94lHuKmgsnDV8M84MgFC+9sO8ZnSafrfNueQeF6f1L3cD/G9u3APcO7yxThQgghhBA1kIRKtCr2ggK27PuM1H8sZtpe9ynHS0ID6LthLa8Fd26eYEoLIfuYY7KJqsnU+OcdE020YptTz/Hge99Wu4Wva5gvz94e53rxrK/Z+VJakxEfLwNeRoNcfRJCCCGEqCNJqESL0TYbpSdOUPDNN+R++iklBw+hi4qIAqKq9LVf0Z1+S/6OT1MlUwWZjpfwXkiFrKOOT156zX3n/QAh3ZomjkYos9k5k1NMWnYhJ7IK+e0nP7i1Rwf7MPWabvzy+l4yAYQQQgghhIdIQiWaVcmRIxSnplKyfz85n36K7ULGxVfw86Hnqo/x6tWr6a6aFGXD66NrT6DKdRwIt73Z7MmU1ppD5y1cyC9xvfsp3/n+p/xiK9mFpaSezefw+fxaJ5R48a4hTI6LlokhhBBCCCE8TBIq0WzO/e1vZL32+kX7lJogxx8sgSZ69R1GzK8exbt376YLynIBVtxaczKljI7kKbw39B0PQ6aD2afpYnHSWpNbVEZadhHHMwv49UdJdZp5rzaP3dyXW4c0022SQgghhBDtjCRUoslpu509i35NwEdfVmsrNUJKd8XRKPjfQANnwuDGmHH8+do/42NqouTlyFfw46dwei+c/R60vaKt2wgY/SiE9XQkU8bmeams1pojFwr4y/pUdh/PIreo7NIr1aBjkDddQ/3oGuZHzwh/xg/qRO8OAR6OVgghhBBClJOESjQprTXJTzxMwGf/datfN1SRFqHY2VdR6KPwNfkyqeckbul1C0M6DPF0EHAmCXa+Ace3Qe7Jmvv5hMAd70JQJ89uvwa5RWV8dzKb0zlFpGUXsSH5LMcyCi653o39OxLoU/7OJxOBzvc+9YgIoF+nQIJ8micBFEIIIYQQDpJQiSaTt2ED5597HtPp0666I1Hw0bRohg6dwjCfMG72DiXUJ5T+4f0J9g72fBBbl8KuNyH/TC0dlOOq1ICfwsDbwT/c8zE45RaWsea7NPal5bDxx3MUltou2r9PxwC6hPrRLyqQa3tHMLxnuDwDJYQQQgjRykhCJZpE3saNpP/6N46rQ06ZgZD30u94b+A9Tbvx4lzH7XwH1sPO16q3G72h781w1SyIvhL8wpokjANn81mblM7xjAJOZReSnH7pl+qOviKCa2LCuGd4d8LkvU9CCCGEEK2eJFTCo6xZWWS8/ApZH3yAciZTpUbYOlDxv4ld+WjA9Kbb+Jl98M2rkLwGbKXV23vfCENnQ8+fgKnpkhWrzc4fPkvhg5213Fro1CnYh1G9I4gO8aVziA8je0XQNcyvyeISQgghhBCeJwmV8AitNWlrPyLvT89jsBRSfmPamVD4wz1GCoLMvHnjn5pu6vNjW2HFLe4TTJQL6QY//w8ERTfJptOyC1m37wxnc4s4l1fChpSztfbtHOLL7Vd1ZkzfDgzuEoxJ3gclhBBCCHFZa3RCpZSyVKnyBV7RWs9VSsUAx4DKT9s/q7V+prHbFa3HZwf+xeFXlzJ+QyaV04PDUfC3KUZGDprA9NjpDI4c3DQBaA07XnJPpjoMgC7x0OUax7NRXp6/8lNUamPZ9uP8bdNBSqw1JHKAr9nIn346kO7hfnQJ9SMywFuegxJCCCGEaEManVBprV1zMiul/IFzwMdVuoVora2N3ZZoXQ5mH+QXm37B3R+cYfz+imelzobA8v8zcHRAGG/e9BZ9w/o2TQDnUuB/zzumQS/Oqaif/CJcndA028RxNe6/+8/z/z7/kZNZhbX2G9krnFfvuZpgX5l5TwghhBCirfL0LX93AOeBrR4eV7QyZfYy/vrPh5i58QxDD1YkUxeifNn9xM3c1fc6ru18LQFeTfAOJK1hx8uwaSHYq+Tp/SY1STJlt2tOZBXy3Yls3vjfUQ6cy3dr7xDozbSh3ejdIYCOQT50CvaR56GEEEIIIdoBTydUs4AVWlea2s3hhFJKAxuBBVrrjKorKqVmA7MBunXr5uGwhCedzTvNB3+awcOfnnar1wP7MPqDj7nOq4lnp9v9Fnz5O/c6r0CIuRYm/c3jm/v6cAa/+9cPHM+sfjUqwNvE4+P7cdc1XTHL81BCCCGEEO2Oqp77NHAgpbrheF6qt9b6mLMuAOgHJAHhwMtAoNb6pouNFR8fr/fs2eORuIRn5RRl83nC/3H1Pvfkwu+aa+j89xcxhYY23cbtdtjyLGz5S0Vd9FUwYYlj+nOD5xKa3MIy9pzI4tkNqRw8V/UxQYeEkTH88vpedAj08dh2hRBCCCFE66OU+lZrHV9TmyevUM0EtpUnUwBaawtQnhmdU0r9CjijlArSWl/6pTyiVbFkn+d/c253S6bK/Lzo/OhjhN19d9PN4AdwYAN89Sc4+31FXXhvmLUOvD13W+HXhzN4am0yRzMKamwf1TucKzoEcsfVXRjYuQleRCyEEEIIIS4rnk6o/nKJPuWXw2Sas8tM8ldryH3iD1yRWzGbXUZcV0a+8wnGgCZ4Tqpc3mn45hXY/g/3+rBecM+aRiVTJVYbu49lc+h8PkcvFJB0Kocf0nNr7BsV5MOr91zFld2a8AqcEEIIIYS47HgkoVJKjQQ6U2V2P6XUMCAHOASEAn8HErXWNZ+1ilZpxb536Pub5wkrqqg7NK4fk174GIOpCV9ltu8j+Nfs6vVDZ8PY34NPw64QldnsvPLVEd7dfoycwrJa+/XtGEiPCH+u7xfJ1GvkuT4hhBBCCFGdp86GZwGfaK3zq9T3BP4MdADycExKMc1D2xTNIPPIj8TMXkJQpWQq5Zf/xx0P/6P2lRqrxAKfzIYDX7jXR/SFqf+EyD6NGv63n/zA6m/Tam3v3SGAZ28fxNXdwxq1HSGEEEII0fZ5JKHSWj9YS/1KYKUntiGaX9769Zx7dD6R9oqJS8puHNnEyVQ+vHcbpO1yr//Jk3DtPDB5N3jorw6c5+nPUjhRaba+ziG+XNs7gp6R/vSKDCCua7BMMiGEEEIIIeqsCe/XEpcbbbNRduoUJcePk/fFv8lbt87tYbfc2M4MXfJK02z8h9Xw3/8HOSeqt81OdMzi1wgbks8y55/futUN6hzM2odGYTTII31CCCGEEKJhJKESABQlp5A+bx5laTXfCvfxdSZ+8beVGLwbfoWoRpbzsOsN+N/zNbf/7iyYfRs0tN2u2XsqhyPnLTy25nu3trguwSy+bZAkU0IIIYQQolEkoRJYtmzh1INzam7zgcV3Guk7ehKRfpGe3fD3H8NnvwJrsXt9RF+IuAJGPtzgZKrEamPGW7vYdTzLrd6gYONvxtArsglnJhRCCCGEEO2GJFTtXN6XX5L+8CNudWVxfdjpncYPkcXs7Kso9TXz6lXzPLfR/HOw/jH4cW31tp8thwFTGjW83a751Qd7qyVTAOMHdpJkSgghhBBCeIwkVO2YvaSEM0/+zq3u4Jwb+X3oV86SAZPBxJIxS+jo39FzG15zHxzfWlEO7ARDpsMV46DbsEYP/9JXh9n44zm3uodvuIK+HQP5v/4dGj2+EEIIIYQQ5SShascyXn0Vu8XiKr+REMUmVzIF/mZ//vqTvzIieoTnNmqzwqmdFeXIWJi5FgKjGjVsqdXOJ9+lseXgBdYnn3XVzxrRnYWTB2CQZ6WEEEIIIUQTkISqnSpLTyfztddd5Y1XGdnUKcNV7h3Smxevf5FuQR58oW3OKdj7T7CVVtT94mswGBs8ZHGZjU37z/GrD/ZWa4vvHspTk/pLMiWEEEIIIZqMJFTtVNY/33ctl3oZWHWtY9mkTDwa/yhT+03FbDB7boN73oHPf+1e1+uGBidT2QWl/Df1PPM/3ldj+6DOwbx091WYjIYGjS+EEEIIIURdSELVDhVs307Wu++6ysvGQq6/wqiM/HPiPxkQPsBzGyuxON4vtet193qDGUb/pt7D5RWXMfeDvfzv0AW0rt7uYzbw8YMjGdg5CKXkypQQQgghhGhaklC1I9pm4+wf/0jOyg/d6rf1dyQe/cP7ezaZKsqGN2+ArCPu9VfNghEPQWTfOg+ltWb59uMs/fIg+SXWau29Iv15Y2a8zOAnhBBCCCGalSRU7YS2Wjk1ezYF23e41f93sKLY25FQ3X7F7Z7boOU8rP65ezIVGeu4KhV3Z72GKiq18fx/DvDO18eqtT18wxXc0K8Dg7uGNDJgIYQQon0qKysjLS2N4uLiS3cWog3z8fGhS5cumM31e+xFEqp24tTsB92SKUNEOAtvKyUlrBCA8THjub2PhxKqve/DvxdAWUFFXc/rYca/oJ634b2w8SCvbTlCqdXuVv/zUT34/cRYmXBCCCGEaKS0tDQCAwOJiYmR2+VFu6W1JjMzk7S0NHr06FGvdSWhagfOL1lCwfbtrnLgzTfz5c9iSDn8FqCI9o/m6ZFPe2ZjFw7Ap790r4saBLe9Ua9kymqz887Xx/j7fw9Va3vvvqGMviKysZEKIYQQAiguLpZkSrR7SinCw8O5cOFCvdeVhKoN01pz7pk/kv3BB646r+7dOf6bn/L3xIdddbPjZuNn9mv8BnPTYW2VZOqOdyD2FjDW7dJpXnEZ/2/dj2xIPoulyrNSD17Xk7uHdaN7uH/jYxVCCCGEiyRTQjT870ASqjbM8tVXbsmUITCQ4Lf+wV1b7sambQB0DujM5F6TG7ehgkzYugS+WwGlFS8K5oaFMLDutxGezyvmvuV7+CE9163ey2Rg/SOjZcIJIYQQQgjR6khC1YblfPKJa9lnwACC//YXZu35NcW2iodOX7z+RbyMXg3fSHEuLJ8M51Pc64dMd8zkV0fJ6bnc+foOCkttrrqoIB/Gxnbg4bFXEBXs0/AYhRBCCCGEaCLy1tM2StvtFH37nasc8dAvWfDjYo7nHXfVPTTkIfqG1X3q8hp99Wf3ZCqwE8z8DKa8AibvS66eXVDKun2nmfH2Trdk6vcTY9nx27H8+aeDJJkSQggh2rFt27YxcuRIgoODCQsLY9SoUezevdvVXlBQQEBAABMmTKi2bkxMDL6+vgQEBBAVFUVCQgIWi6VaPyEaQxKqNirzzbewZWcDYPP15hcFb7Hr7C5X+y29buHBuAcbt5HzqbDztYpy3F0wLxl6jqnT6n/+936u+uNG5q7cS3Zhmat+6c8Gc//onnI/txBCCNHO5eXlMWnSJObOnUtWVhbp6eksXLgQb++KL21Xr16Nt7c3X375JWfOnKk2xrp167BYLCQlJbF3714WL17cnLsg2gFJqNqg0rR0Lvztb67y/64oZV/WD65yoFcgi0YualzCknkEXhlWUfYOhkl/BeOl7yI9n1/M/1v3I2/87yhaV9T7mA2sfGA4t1/dpeFxCSGEEKLNOHjwIADTpk3DaDTi6+vLuHHjiIuLc/VZvnw5c+bMIS4ujvfff7/WsaKiorjppptISkpq6rBFOyPPULVBme++Q+VM5aPRFXlzR7+OLBmzBJOhgYc+5yTsehO2/929fsLz4HXpmQLf/foYi/+dSqnN/b1S9wzvxozhMfSNCmxYXEIIIYTwiJgnvmi2bR3/y8SLtvfp0wej0cisWbO46667GD58OKGhoa72kydPkpiYyEsvvURYWBjLly9n/vz5NY6VlpbG+vXrGTt2rEf3QQi5QtXGlBw/Ts77FTP7rR2uyApSRPpG8s5N77Dxjo0M6TCkfoNaS2H3W/D2TfDikOrJ1OhHYfDUSw6TnJ7LonU/uiVTfl5GNv76Ov44ZZAkU0IIIYRwExQUxLZt21BK8cADDxAZGcktt9zCuXPnAFixYgVxcXH079+fadOmkZKSwt69e93GmDJlCoGBgXTt2pUOHTqwaNGiltgV0YZ5JKFSSiUqpYqVUhbn50ClthuUUqlKqUKl1FdKqe6e2Kao2d41r7uWi82wf1gn/nb939j0s01cE3VN/W7z0xp+/NRxa98Xj8Kpb0BXTByBTwhMXAo3/KEOQ2me/qxi8opAbxN/nDKQxAU/4YqOkkgJIYQQomaxsbEsW7aMtLQ0kpOTOX36NPPmzQMcCdX06dMBiI6OZsyYMSxfvtxt/bVr15Kfn09iYiKpqalkZGQ09y6INs6Tt/z9Smv9VuUKpVQE8AlwP7AOeAb4CBjuwe0Kp3RLOqnfrKf8yaZPrjPz1uz1DZsW/cJBWPcInNxeva3bSBg5F3qNBfPFZ+DTWrP62zRe3XKEoxcKXPWv3nM1114RUf+4hBBCCNGkLnUbXkvq168fCQkJvP7662zfvp1Dhw6xePFili5dCkB+fj4pKSksWbIEk8n9NHfMmDEkJCQwf/581q5d2wLRi7aqqZ+hug1I0Vp/DKCUehrIUEr101qnNvG225XzheeZ8cmdvPxDiavunp8907Bk6tQuWDYJbBVj4RPsuLVvyHTwv3gipLUmr8jKmbwiHl/9PfvS3F/UOymukyRTQgghhLik1NRUvvjiC6ZOnUqXLl04deoUK1euZPjw4Sxfvpwbb7yRFStWuPoXFRURFxfH+vXrmTx5crXx5s2bR0xMDElJSQwZMqQZ90S0ZZ5MqBYrpf4CHAB+p7VOBAYA+8o7aK0LlFJHnPWSUHlAXmke7//4Pq/se4X5aypux9N+PvS5toHfMO183T2ZGv4QXDcf/MIuutoHO0/y9rajpGUXUWK119hnXP+OPHdHXI1tQgghhBCVBQYGsnPnTl544QVycnIICQlh0qRJPP/883Tv3p0VK1YQFRXlts6MGTNYvnx5jQlVZGQkM2fO5JlnnmHNmjXNtRuijfNUQvU48CNQCtwFrFNKDQECgAtV+uYC1R6aUUrNBmYDdOvWzUNhtX3zE+ez48wO7v7KxtCDFTP7Rc37DcrUgMN7NhmSV1eUb30Zrrznoquczili+fbjvP6/o7X2GRoTxl/vGkLnEN/6xySEEEKIdqlz586sWrWqxrZs5/s2q3rllVdcy8ePH6/W/uqrr3okNiHKeSSh0lrvrFRcrpSaBkwALEBQle5BQH4NY7wBvAEQHx+vq7aL6pLOJ7HjzA76ndJM+abiV+Z37UjCZs6o32Baw4H18O9KU412ueaiyVRxmY1nPv+RlbtOYq9yxPy9jHQI8iEy0JvhPcOZd8MVGAzyol4hhBBCCNG2NNUzVBpQQAowq7xSKeUP9HLWi0Z678f3MNo0C1bb3Oo7P7+k/oNt/iNsrbSe0Qv+7+lauy/5zwFe+upwtfreHQJY8fOhRMuVKCGEEEII0Q40OqFSSoUAw4AtgBWYClwHzAOygOeVUrcDXwB/AL6XCSka70TeCb488SXjv9MEFlfUd37xRUyVXnh3SWf2wYYn4cQ29/qfLYeYa2tcZc/xrGrJ1MDOQUwb2o3bruyCr5ex7tsXQgghhBDiMuaJK1Rm4I9AP8CGY7KJKVrrAwDOZOol4J/AThzPWIlG+ujAR3TI1ty7qWLyB9+rriLopnF1HyT137BqJtjL3Osf2g2Rfap1t9s1L2w8yCuJ7snUgpv68uB1PTEZ5T3RQgghhBCifWl0QqW1vgBcc5H2TTiSLeEhP2b+yOc/fMwjX7jf6hfxizl1G0Br+PpF2LSwok4ZHFOij/09BEbVuNrq79LcrkwZDYr1j4ymj7yYVwghhBBCtFNN/R4q4WGnLaf5w1t387t1JXSr9KLv6Gf/QsDo0XUb5IvfwJ533Ovu2wRdrq51lcQD53ls9feucpi/F3+5bZAkU0IIIYQQol2ThOoysnfX5+z66+9YlFSKodKsekGTJxN0yy11GyT92+rJ1N2rak2mkk7l8PGeU7y/86SrzqBg46+vIzzAu767IIQQQgghRJsiCdVloiglBZ+ZC7iuUp02GYlasIDQmTNRqg5TkmsNm56uKId0h5mfQlgPt252u2b38SxeSTzCloPurxEzGRS/nxgryZQQQgghhBCAzCLQytksBaQ/9hjHb7/Drb44tju9v/iCsFmz6pZMAex4CY79r6I8/eNqyVSJ1UbCst1MfeObaslU1zBfNswbTcIo93WEEEIIIUTjbN26lb59+7rKMTExbNq0qQUjEnUlCVUrd+ap35P32Tq3uq+v8mHw6i/w6t697gN9uxy+/H1Fud8kiOzr1mXdvtPE/3ET/6uSSN12ZWfeSYjny3lj6N1BnpkSQgghRPPYtm0bI0eOJDg4mLCwMEaNGsXu3bvd+hQUFBAQEMCECROqrR8TE4Ovry8BAQFERUWRkJCAxWJpUCxz5swhICCAgIAAvLy8MJvNrvL48eM5fvw4SilXXfnno48+co2xa9cuJkyYQEhICGFhYQwdOpR3330XgNGjR3PgwIEGxSZaliRUrVjuZ5+Rv36Dq/xjV1j8MwPnHr4Dg7Ee73o6tBHWPVxRDuoMt/zDrcuWgxeYu3Iv+cVWV12nYB/W/GIkL0wdwth+HeX9UkIIIYRoNnl5eUyaNIm5c+eSlZVFeno6CxcuxNvb/bGD1atX4+3tzZdffsmZM2eqjbNu3TosFgtJSUns3buXxYsXNyie1157DYvFgsVi4cknn2Tq1Kmu8vr16139cnJyXPUWi4WpU6cCsGPHDsaOHcuYMWM4fPgwmZmZvPrqq27rNpTVar10J9FkJKFqpbTWXPjHS65yUg/F0/eY2NvbwN2xd9d1EPhqMbxf6XbBiL4wYy34hbl1feWr6u+W2vrY9VzdvR4vCRZCCCGE8JCDBw8CMG3aNIxGI76+vowbN464uDi3fsuXL2fOnDnExcXx/vvv1zpeVFQUN910E0lJSU0Zdq0WLFjArFmzePzxx4mIiEApxdVXX82qVasASExMpEuXLm7rJCUlERcXR3BwMFOnTqW4uNit77PPPktUVBT33ntvs++PqCCTUrQy2molZ/UactasoezUKQCKzfDirY7cd0b/GfQIruMzTN8thy1/ca/72bJqL+398XQeO49lucqfz72WgZ2DG7wPQgghhLiMPd2M5wBP59ba1KdPH4xGI7NmzeKuu+5i+PDhhIa6f9F78uRJEhMTeemllwgLC2P58uXMnz+/xvHS0tJYv349Y8eO9egu1EVhYSE7duzgmWeeqdd6q1atYsOGDfj4+DBq1CiWLVvGnDmO946ePXuWrKwsTpw4gd1ub4qwRR3JFapWJn3BAs4+/TTFP/zgqvsqTlHgq3j4yod57JrH6jZQcR589Wf3uolLoWN/t6rTOUXc/dY3rnLPCH9JpoQQQgjR4oKCgti2bRtKKR544AEiIyO55ZZbOHfunKvPihUriIuLo3///kybNo2UlBT27t3rNs6UKVMIDAyka9eudOjQgUWLFjVp3BEREYSEhLg++/fvJzs7G7vdTqdOneo11sMPP0x0dDRhYWFMnjzZ7eqawWBg0aJFeHt74+vr6+G9EPUhCVUrUnL0qNszUwBfxyrev97A8E7DuW/QfXUf7NOHwOL8B8c3DOYlwzX3V+s278MkcgrLAPAyGnhqcv9qfYQQQgghWkJsbCzLli0jLS2N5ORkTp8+zbx581ztK1asYPr06QBER0czZswYli9f7jbG2rVryc/PJzExkdTUVDIyMmrc1tatW10TSQwYMKDBMWdkZJCTk+P6xMbGEhoaisFgqPEZr4uJiopyLfv5+blNqBEZGYmPj0+D4xSeI7f8tSKlJ064ls3du7F4hh+7bI5nm+4bdB8GVcf8N+Mw7P+sojzuGQjp6tZFa83v1iaz63jFrX6v3nMV1/ft0PAdEEIIIcTl7yK34bWkfv36kZCQwOuvvw7A9u3bOXToEIsXL2bp0qUA5Ofnk5KSwpIlSzCZ3E9zx4wZQ0JCAvPnz2ft2rXVxh89enSDZwC8FD8/P0aMGMGaNWu4/vrrPTJmnV+bI5qcXKFqRfLWfe5aNsdf6UqmAK4IuaJug+SmwYpbKsodBsCV97iKJVYbz/8nlUFPf8kHO0+66q/uHsoNsR0bHrwQQgghhAelpqaydOlS0tLSADh16hQrV65k+PDhgGMyihtvvJEff/yRpKQkkpKSSE5OprCwsNaZ8+bNm8fGjRtbZGKK5557jmXLlvH888+TmZkJwL59+7jrrruaPRbhWZJQtRLFqank/fvfrvICn4rkKso/inDf8LoNtPaXkJdeUb7O/cHMpz/7kZe/OoKlpGJ6zRE9w3np7isbFrgQQgghRBMIDAxk586dDBs2DH9/f4YPH87AgQNZunQpxcXFrFq1irlz5xIVFeX69OjRgxkzZlS77a9cZGQkM2fOrPfkEPUREhLi9h6qF154AYCRI0eyefNmNm/eTM+ePQkLC2P27Nk1vj9LXF6U1rqlY6gmPj5e79mzp6XDaFZn//Rnst97D4ATVwSx4I5CV9vd/e7mt8N+e+lBinPhL90qytf/DsZUTGKxfPtxFn6W4rbKzQOiWHLnYAK85e5PIYQQoj3av38/sbGxLR1Gu7d582buv/9+jh492tKhtGu1/T0opb7VWsfXtI6cRbcC9tJSVzIFsDbWQvnFw2FRwy49s5+1BDb+Ab5bUVFn8oHrFmC3a5JP5/LY6u9JPZvvau4XFcinvxqFt0le1iuEEEII0dKSk5Pp0aOOr8YRrYokVC3MeuECJ2bOqigbYHcfx0OGN3a/kaVjll76ocOv/w47X3OvGz2fD3efYvH6VHKLytyaekb689GDIySZEkIIIYRoBR555BE+++yzWm9VFK2bPEPVgnRpKem/eZTSY8dcdV/3V5SaFWaDmXlXzbt0MlViga/+WFH2CYbb3iQxahZPfPJDtWQK4G9ThxDsa/bUbgghhBBCiEZ48cUXOXbsGNddd11LhyIaQK5QtaALL79C4e7drvLHo02sGel4pu2vP/kr3YK61byi3QbJa+DEdvj2Xbcm/dAu/vpNHn//7263+pG9whl9RSQPjO6BySh5tBBCCCGEEJ4gCVUL0DYbGa++RqbzPQoAW0eH8vG1+YCid0hvrutSyzcUeWfglWGOCSiqip3MqtRS/v7fQ64qpeCrR39CTIS/h/dCCCGEEEIIIQlVC0h/dD75Gza4yjZvE/8cnAc4bu97aMhDNd/ql/oFfPIglOa715t84KpZ7Oj9ax5/5zu3pmdvi5NkSgghhBBCiCYiCVUzKzl82C2ZAnj7ng5kB54HwNfky9huY6uveHgTfHQPaLt7/fQ10PUatp0q4563d7qqw/y9SFzwE4J85FkpIYQQQgghmookVM2sqMqbuZPemsumI68CYFAG1t+2HoOq4Rmn/y2tSKbM/vDTVyH2FlAKm13z0Acb3bq/fPdVkkwJIYQQQgjRxBo9O4FSylsp9bZS6oRSKl8ptVcpNd7ZFqOU0kopS6XPU40P+/KV/9/NruUDU4fyZ2cyBTBrwCzCfcOrr7R/HZzcXlGe/jH0v9XxgBTw2pYjbrP5ffDAMEb0qmEcIYQQQgjRKm3dupW+ffu6yjExMWzatKkFIxJ15Ynp3kzAKWAMEAw8BaxSSsVU6hOitQ5wfp7xwDYvS/aSEgq2VyRGb/l/61pWKO4fdL/7CgUZsGqW41a/coGdIGaUq5hVUMqriUdc5RnDuzOyV4TngxdCCCGEaGbbtm1j5MiRBAcHExYWxqhRo9i9230m44KCAgICApgwYUK19WNiYvD19SUgIICoqCgSEhKwWCwNimXOnDkEBAQQEBCAl5cXZrPZVR4/fjzHjx9HKeWqK/989NFHrjF27drFhAkTCAkJISwsjKFDh/Luu44Zm0ePHs2BAwcaFFtDvf322/Tr14/AwEA6duzIxIkTyc93PKufkJCAl5eX274MHjyYrVu3usr+/v7V9vnkyZPNug+tQaMTKq11gdb6aa31ca21XWv9OXAMuLrx4bUtli1b0CUlAOSEeXGig6M+0BzIexPeI8grqKJz1jF4ZTj8uNZ9kHEV75w6m1vM9Ld2YimxAtA93I8nxvdryl0QQgghhGgWeXl5TJo0iblz55KVlUV6ejoLFy7E29vbrd/q1avx9vbmyy+/5MyZM9XGWbduHRaLhaSkJPbu3cvixYsbFM9rr72GxWLBYrHw5JNPMnXqVFd5/fr1rn45OTmueovFwtSpUwHYsWMHY8eOZcyYMRw+fJjMzExeffVVt3Ubymq11nudLVu28OSTT7Jy5Ury8/PZv38/d955p1ufxx57zG1f9u3bx+jRo13llJSUavvcrVstr/1pwzz+QiKlVEegD5BSqfqEUipNKfWuUqrdXj7J/uf7ruX/9isDpTAZTHw06SMGRw5277znbSi4UFHuONAxAcWgOwCw2zU3/nUL+8/kubrMH9cXf295LE4IIYQQl7+DBw8CMG3aNIxGI76+vowbN464uDi3fsuXL2fOnDnExcXx/vvv1zQUAFFRUdx0000kVXmevbksWLCAWbNm8fjjjxMREYFSiquvvppVq1YBkJiYSJcuXdzWSUpKIi4ujuDgYKZOnUpxcbFb32effZaoqCjuvffeeseze/duRowYwZVXXglAWFgYs2bNIjAwsJF72v549OxbKWUG3geWa61TlVIBwDVAEhAOvOxsv6mGdWcDs4E2mdkWp6ZSuGuXq7x1oCOXvf2K2+ka1LWiY246bP8H7Kx4tor+U+DO5W7jrfv+NPnFFd9GPDimJ5MHRzdJ7EIIIYRoPwYtH9Rs2/ph1g+1tvXp0wej0cisWbO46667GD58OKGhoW59Tp48SWJiIi+99BJhYWEsX76c+fPn1zheWloa69evZ+zYGmZTbmKFhYXs2LGDZ56p35Mvq1atYsOGDfj4+DBq1CiWLVvGnDlzADh79ixZWVmcOHECu91+iZGqGzZsGE899RQLFy5k3LhxxMfHV7v6J+rGY1eolFIG4D2gFPgVgNbaorXeo7W2aq3POevHKaWCqq6vtX5Dax2vtY6PjIz0VFitRu6/1rqWj0bB6XDFkMghPBr/aEWnwix483r3ZEoZYMLzbmNlWEpY/O9UV3lEz3CeuFlu9RNCCCFE2xEUFMS2bdtQSvHAAw8QGRnJLbfcwrlz51x9VqxYQVxcHP3792fatGmkpKSwd+9et3GmTJlCYGAgXbt2pUOHDixatKhJ446IiCAkJMT12b9/P9nZ2djtdjp16lSvsR5++GGio6MJCwtj8uTJblfXDAYDixYtwtvbG19f33rHOXr0aD755BO+++47Jk6cSHh4OL/5zW+w2WyuPkuWLHHbl1mzZtV7O+2BRxIq5XgL7dtAR+B2rXVZLV11+Sqe2O7loiw9nexKDySuG+r4tT8x7Al8TZX+AD6bC5aKfyTw7wA/WwYBHbDbNf/df47F6/cT/8dNnM1zXPL19zLy8vSran4RsBBCCCHEZSw2NpZly5aRlpZGcnIyp0+fZt68ea72FStWMH36dACio6MZM2YMy5e739Wzdu1a8vPzSUxMJDU1lYyMjBq3VXmyhQEDBjQ45oyMDHJyclyf2NhYQkNDMRgMNT7jdTFRUVGuZT8/P7cJNSIjI/Hx8al13QEDBrj2Z+vWrTX2GT9+POvWrSMrK4tPP/2UZcuW8dZbb7na58+f77YvVX+3wsFTt/y9CsQC/6e1LiqvVEoNA3KAQ0Ao8HcgUWud66Httnq2/HxO/PznaOc9r6fDYHt/xejOoxkQXumP9dvlkPp5RXnog3DTn8BoprjMRsK7u/jmaFa18eff1Jcwf6+m3g0hhBBCtBMXuw2vJfXr14+EhARef/11ALZv386hQ4dYvHgxS5cuBSA/P5+UlBSWLFmCyeR+mjtmzBgSEhKYP38+a9eurTZ++WQLTcHPz48RI0awZs0arr/+eo+Meakv08snjKgLg8HADTfcwNixY0lOTm5saO2OJ95D1R14EBgCnK30vqnpQE9gA5APJAMlwLTGbvNycuEf/6DsRMX0keuGGtBKce/ASg8P7vsI1j1cUY6MhZsXg9HxYt5Pk9JrTKZ+0jeShJExTRW6EEIIIUSLSU1NZenSpaSlpQFw6tQpVq5cyfDhwwHHZBQ33ngjP/74I0lJSSQlJZGcnExhYWGtM+fNmzePjRs3tsjEFM899xzLli3j+eefJzMzE4B9+/Zx1113NXssAJ9++ikffvgh2dnZaK3ZtWsXW7Zscf1+Rd15Ytr0E1prpbX2qfSuqQCt9fta65Va6x5aa3+tdSet9Uyt9VlPBH65KN73vWs5PQwS4xSjOo8ivmO8o3L32/Cv2RUrhPaAe1aDwQhAXnEZf/x8v6t5UOdg/jHtSr789XW8m3CN3OonhBBCiDYpMDCQnTt3MmzYMPz9/Rk+fDgDBw5k6dKlFBcXs2rVKubOnUtUVJTr06NHD2bMmFHrrWmRkZHMnDmz3pND1EdISIjbe5leeOEFAEaOHMnmzZvZvHkzPXv2JCwsjNmzZ9f4/qzmEBoayptvvskVV1xBUFAQ99xzDwsWLHDdQgmOJLDyvkREtNvJui9Kaa0v3auZxcfH6z179rR0GI1mLynh0IiR2AsLAZjzkJGsIMV3M77DfP4AfJwAmYcqVgjtAQ9sBr8wACwlVu56YwfJ6RVTo380ezjDeoY3524IIYQQog3bv38/sbGxLR1Gu7d582buv/9+jh492tKhtGu1/T0opb7VWsfXtI68tKiJaKuV0/MXuJKprADICnRMk27GCB/cCXnp7ivNWudKpgD+sn6/WzJ1TUwo8TFhCCGEEEKItiU5OZkePXq0dBiiASShagLWjAxO3vtzSg5VXH3aNkDha/bjl7Ez4PXR7slUSHe44x0IqXgf1cnMQlbuOuUqzxnTiwU39cVokFv8hBBCCCHakkceeYTPPvtMZtG7TElC5WG6tJS0uQ+7JVPFZtjW38Diob+nwzsTIb/SlJmxk+HO96DKs1CvbjmCze64HXNw1xBJpoQQQggh2qgXX3yRF198saXDEA0kCZWHnfrFLymq9EK5/V3gxVuNmKI6ct2GRe7JVHA3uONdVzJVYrXxadJpvj6cwadJp13dnri5nyRTQgghhBBCtEKSUHlQ0Q8/UPD1165ygTcsvc1Inr/iFa8YzJmrKzoP+hlMeQ2MjkNw+Hw+v1m1j+/T3F/RFd89lOE95bkpIYQQQgghWiNJqDyk5MgRjv/sTlc5zxcW3OdIpqYG9WP03krJVKch8NM3wGCgxGrj9/9KZvV3aVSdcLFXpD9LfjZYpkYXQgghhBCilZKEygNslgKOTpzkVvfHu4xkByr6BHRj7g+bKhq6jYB71pByNp//JJ/lrW3HKCy1ua17dfdQnpzQj8FdQjAZG/2qMCGEEEIIIUQTkYTKAwr37HYrfzZMcTxKoVD89WgKwXa7oyG0B9uvXcZzbyaRdCqn2jjXxITy+4n9Gdw1pOmDFkIIIYQQQjSaXP7wgOLvf3AtH+xu5p9jjQDcGtSXbgXZrras+Hnc917NydRP+kby0ewRkkwJIYQQQrRDW7dupW/fvq5yTEwMmzZtusgaorWQhMoDir7/3rW8Ic5x+57JYGL+kSRXfenox7kpsQtFZRW3990yOJoX7xrC5kfH8G7CNRhkJj8hhBBCCJdt27YxcuRIgoODCQsLY9SoUeze7X5nUEFBAQEBAUyYMKHa+jExMfj6+hIQEEBUVBQJCQlYLJYGxTJnzhwCAgIICAjAy8sLs9nsKo8fP57jx4+jlHLVlX8++ugj1xi7du1iwoQJhISEEBYWxtChQ3n33XcBGD16NAcOHGhQbA319ttv069fPwIDA+nYsSMTJ04kPz8fgISEBLy8vNz2ZfDgwWzdutVV9vf3r7bPJ0+erHcctW0LqPZ7jYmJ4S9/+Uu1MRITE1FK8dxzz1VrU0rh7+/vGuP++++vd4wXIwlVI9ny8ynctctVPhztSIqmEkyw5byrfvR/Y7iQX+Iqv5twDX+fdiW3DulMz8gAmXhCCCGEEKKSvLw8Jk2axNy5c8nKyiI9PZ2FCxfi7e3t1m/16tV4e3vz5ZdfcubMmWrjrFu3DovFQlJSEnv37mXx4sUNiue1117DYrFgsVh48sknmTp1qqu8fv16V7+cnBxXvcViYerUqQDs2LGDsWPHMmbMGA4fPkxmZiavvvqq27oNZbVa673Oli1bePLJJ1m5ciX5+fns37+fO++8063PY4895rYv+/btY/To0a5ySkpKtX3u1q1bg/ahpm1VVr6N1atX88wzz7Bx40a39uXLlxMWFlbry5H37dvnGvutt95qUIy1kYSqkXI/+wxdWgrAyUg4GwpPqg4sOFLxLqq1tpGcswe7yhPjOnF9vw7NHqsQQgghxOXi4MGDAEybNg2j0Yivry/jxo0jLi7Ord/y5cuZM2cOcXFxvP/++7WOFxUVxU033URSUlJThl2rBQsWMGvWLB5//HEiIiJQSnH11VezatUqwHGFpUuXLm7rJCUlERcXR3BwMFOnTqW4uNit77PPPktUVBT33ntvvePZvXs3I0aM4MorrwQgLCyMWbNmERgY2Mg9bVrx8fEMGDDA7TgWFhayevVqXn75ZQ4dOsSePXuaNSaZlKIRtNVK5rJlrvKmIQZutBqYllZxEPfae7OwLAGAED8z943qwX2jezRzpEIIIYQQdbO/X2yzbSs2dX+tbX369MFoNDJr1izuuusuhg8fTmhoqFufkydPkpiYyEsvveS6OjF//vwax0tLS2P9+vWMHTvWo/tQF4WFhezYsYNnnnmmXuutWrWKDRs24OPjw6hRo1i2bBlz5swB4OzZs2RlZXHixAns5ROg1cOwYcN46qmnWLhwIePGjSM+Pr7a1b/W6JtvviE5OZnf/va3rro1a9YQEBDAz372M95//31WrFhBfHy823rXXXcddrudkSNH8sILLxATE+OxmOQKVSOcWvYG1lNpAFh84Ks4xSNn013th+3R3F76NL27d+X3E2NJnP8T5t5wBX5ekscKIYQQQlxMUFAQ27ZtQynFAw88QGRkJLfccgvnzp1z9VmxYgVxcXH079+fadOmkZKSwt69e93GmTJlCoGBgXTt2pUOHTqwaNGiJo07IiKCkJAQ12f//v1kZ2djt9vp1KlTvcZ6+OGHiY6OJiwsjMmTJ7tdlTEYDCxatAhvb298fX3rHefo0aP55JNP+O6775g4cSLh4eH85je/wWareN5/yZIlbvsya9asem+nri61rYiICHx9fRkxYgS//OUvmTJliqtt+fLlTJ06FaPRyN13383KlSspKytztW/ZsoXjx4+TmppKdHQ0kyZNatBtkrWRhKqBzqfupWDJP1zl/w1UTC3Kp7vz4GyzDeCX5mf4+BfXsuYXI7l/dE9C/LxaKlwhhBBCiMtObGwsy5YtIy0tjeTkZE6fPs28efNc7StWrGD69OkAREdHM2bMmGrP0Kxdu5b8/HwSExNJTU0lIyOjxm1VnmxhwIABDY45IyODnJwc1yc2NpbQ0FAMBkONz3hdTFRUlGvZz8/PbUKNyMhIfHx8al13wIABrv3ZunVrjX3Gjx/PunXryMrK4tNPP2XZsmVuzxfNnz/fbV9qez7pYt5//323yTtqc6ltZWRkYLFYWLJkCYmJia6E6dSpU3z11Veu/w5uvfVWiouL+eKLL1zrXnfddXh5eRESEsKLL77IsWPH2L+/9quj9SWXShqgzFrMxqdmcZWznOcLfoMLsObdwsSSfpzQHbDgx6ezR8k06EIIIYS4rFzsNryW1K9fPxISEnj99dcB2L59O4cOHWLx4sUsXboUgPz8fFJSUliyZAkmk/tp7pgxY0hISGD+/PmsXbu22vjlky00BT8/P0aMGMGaNWu4/vrrPTLmpSY0K58woi4MBgM33HADY8eOJTk5ubGhuZk+fbor2Wkso9HIo48+yr/+9S9eeeUV5s2bx3vvvYfdbmfy5MmufsXFxaxYscLtKlZlSim01h6JCeQKVYO8/f5dxB6ouIx4anQZm7MeZqnlZlJ0DBb8+O34fpJMCSGEEEI0UGpqKkuXLiUtzfF4xalTp1i5ciXDhw8HHLd53Xjjjfz4448kJSWRlJREcnIyhYWFtc6cN2/ePDZu3NgiE1M899xzLFu2jOeff57MzEzAMfPcXXfd1eyxAHz66ad8+OGHZGdno7Vm165dbNmyxfX7bc2eeOIJnnvuOVfitHDhQtd/A0lJSaxZs4YvvviCzMxMUlJSSEpKwmazYbFYePTRR+ncuTOxsZ57VlASqnr6OukdAlcdxNcxsR+5wfBH32fYrfsBYDIoEkbGcN+1MvGEEEIIIURDBQYGsnPnToYNG4a/vz/Dhw9n4MCBLF26lOLiYlatWsXcuXOJiopyfXr06MGMGTNqvTUtMjKSmTNn1ntyiPoICQlxe5/SCy+8AMDIkSPZvHkzmzdvpmfPnoSFhTF79uwa35/VHEJDQ3nzzTe54oorCAoK4p577mHBggVuV5Oee+45t32JiIhosnjqs62JEycSGhrK888/z/Hjx3nooYfc/ju45ZZb6N27NytXruTcuXNMnTqVoKAgevbsyfHjx/n8888xm80ei1158nKXp8THx+vmnu6wrj6fMoheqRUPsb0yZBLrYn5C346B/Cy+C1Ov6Uqgj+cOkBBCCCFEU9q/f79Hv60XDbN582buv/9+jh492tKhtGu1/T0opb7VWsfXsIo8Q1VX2mbj2Jxb3ZKpH6I7cOCaCST9YpRMOCGEEEIIIRosOTmZHj3kDqfLkSRUdXR+yTOUbD3iKluN8NRVc/l46pWSTAkhhBBCiAZ75JFH+Oyzzxo0i55oeU3+DJVSKkwp9S+lVIFS6oRS6u6m3qan2UuLOP7xR67yzj6K39wxgTuv7Udcl5CWC0wIIYQQQlz2yqfyvu6661o6FNEAzXGF6mWgFOgIDAG+UErt01rXfS7HllRayOrfjmKQcxZNqwG+uL4j8d0T+H+3NvwdBUIIIYQQQojLX5NeoVJK+QO3A09prS1a623AZ8CMptyuJ+1/+Xf0/U9RRbmrorT0IR67qe8l5/8XQgghhLgctMZJyoRobg39O2jqW/76ADat9cFKdfuAy+bSTtnWb/CyOhKnQi/IfvAJPvz1Twn1l+emhBBCCHH58/HxITMzU5Iq0a5prcnMzMTHx6fe6zb1LX8BQG6VulwgsGpHpdRsYDZAt27dmjisujNYK2b1Wx+veHjSFHy9jC0YkRBCCCGE53Tp0oW0tDQuXLjQ0qEI0aJ8fHzo0qVLvddr6oTKAgRVqQsC8qt21Fq/AbwBjvdQNXFcdXfr3ew/sJUybSN66CD8Tf4tHZEQQgghhMeYzWaZrluIRmjqhOogYFJKXaG1PuSsGwxcHhNSAAPv+zUD+XVLhyGEEEIIIYRohZr0GSqtdQHwCfD/lFL+SqlRwK3Ae025XSGEEEIIIYRoDk3+Hirgl4AvcB5YCfzispkyXQghhBBCCCEuosnfQ6W1zgKmNPV2hBBCCCGEEKK5qdY4RaZS6gJwoqXjqCICyGjpIESTkePbtsnxbbvk2LZtcnzbNjm+bVtbO77dtdaRNTW0yoSqNVJK7dFax7d0HKJpyPFt2+T4tl1ybNs2Ob5tmxzftq09Hd/meIZKCCGEEEIIIdokSaiEEEIIIYQQooEkoaq7N1o6ANGk5Pi2bXJ82y45tm2bHN+2TY5v29Zujq88QyWEEEIIIYQQDSRXqIQQQgghhBCigSShEkIIIYQQQogGkoRKCCGEEEIIIRpIEionpVSYUupfSqkCpdQJpdTdF+n7a6XUWaVUrlLqHaWUd3PGKuqvrsdXKTVQKfUfpVSGUkoeMLxM1OP4zlJKfauUylNKpSmlnlNKmZo7XlE/9Ti+dymlDjj/bT6vlFqulApq7nhF/dTn/7+V1tmslNLy99u61eNvN0EpZVNKWSp9ftK80Yr6que5c0+l1OdKqXznOdZzzRlrU5OEqsLLQCnQEZgOvKqUGlC1k1LqJuAJ4AYgBugJLGq+MEUD1en4AmXAKuC+ZoxNNF5dj68fMA/H29uH4fg7nt9MMYqGq+vx/RoYpbUOxvFvswn4Y7NFKRqqrscXAKXUdBzHVrR+9Tm2O7TWAZU+ic0VpGiwup47ewEbgc1AFNAF+GczxtnkZJY/QCnlD2QDA7XWB5117wHpWusnqvT9ADiutX7SWb4BeF9rHdXMYYs6qs/xrbROb+CQ1lo1X6SiIRpyfCut+xvgeq315KaPVDREQ4+vUioAeAWI0FpPaJZgRb3V9/gqpYKB3cBMYAdg1lpbmzFkUUf1PLdKAO7XWl/b7IGKBqnn8Z0NzNBaj27+SJuHXKFy6APYyv+DcNoH1PQtygBnW+V+HZVS4U0Yn2ic+hxfcflpzPG9DkhpkqiEp9Tr+CqlrlVK5QL5wO3A35o8QtEY9f37/TPwKnC2qQMTjVbfY3ul81awg0qpp+R2zlavPsd3OHBcKbXeeYwTlVKDmiXKZiIJlUMAkFulLhcIrEPf8uWa+orWoT7HV1x+GnR8lVL3AvHAkiaKS3hGvY6v1nqb85a/LsDzwPEmjU40Vp2Pr1IqHhgF/KMZ4hKNV5+/3f8BA4EOOL4ImQYsaNLoRGPV5/h2Ae4C/g5EA18AnzpvBWwTJKFysABVH1wOwvEN56X6li/X1Fe0DvU5vuLyU+/jq5SaAvwFGK+1zmi60IQHNOjvV2udDmwAPmyiuIRn1On4KqUMOG7hfERu8bts1PlvV2t9VGt9TGtt11r/APw/4I5miFE0XH3+bS4Ctmmt12utS3F8kRkOxDZtiM1HEiqHg4BJKXVFpbrB1HwrUIqzrXK/c1rrzCaMTzROfY6vuPzU6/gqpW4G3gQmO//HLVq3xvz9moBeTRKV8JS6Ht8gHFeUP1JKncXxHBVAmlKqzT6XcZlrzN+uBuQZ5tatPsf3exzHtM2SSSmclFIf4jjY9wNDgH8DI7XWKVX63QwsA8YCZ4A1wK5LPfwuWlY9jq8CvHHMEJYC+AJaa13SrAGLeqnH8R0LfAz8VGv9v+aOUzRMPY7vdGArcAroBqwAMrXWtzVrwKJe6nJ8nf82d6y0WldgF45biS44v/UWrUw9/nbHA99prc8ppfoBq4GPtdYyi3IrVo/j2xfYC9wCfAU8DPwKiG0rf7tyharCL3GcPJ8HVgK/0FqnKKW6Od+H0A1Aa70BeA7HfxAnnJ+FLRSzqLs6HV+gO45L0+X/GBQBB5o9WlFfdT2+TwHBwL8rvetkfQvFLOqurse3P7Adx60oX+P4232gJQIW9XLJ46sdzpZ/gAvOdc+1lROyNqquf7s3AN8rpQpwnJR/gmMCEtG61fXc+QBwD/AajpkBbwVuaUt/u3KFSgghhBBCCCEaSK5QCSGEEEIIIUQDSUIlhBBCCCGEEA0kCZUQQgghhBBCNJAkVEIIIYQQQgjRQJJQCSGEEEIIIUQDSUIlhBBCCCGEEA0kCZUQQgghhBBCNJAkVEIIIYQQQgjRQJJQCSGEEEIIIUQDSUIlhBBCCCGEEA0kCZUQQgghhBBCNJAkVEIIIYQQQgjRQJJQCSGEEEIIIUQDSUIlhBBCCCGEEA1kaukAahIREaFjYmJaOgwhhBBCCCGE4Ntvv83QWkfW1NYqE6qYmBj27NnT0mEIIYQQQgghBEqpE7W1yS1/QgghhBBCCNFAdUqolFK/UkrtUUqVKKWWVWm7QSmVqpQqVEp9pZTqXqlNKaWeVUplOj/PKaWUh/dBCCGEEEIIIVpEXa9QnQb+CLxTuVIpFQF8AjwFhAF7gI8qdZkNTAEGA3HAJODBRkUsxKVo7fhhs7VwIEIIIYQQoq2r0zNUWutPAJRS8UCXSk23ASla64+d7U8DGUqpflrrVGAWsFRrneZsXwo8ALzmsT0Q7UtJPuSfReemk5d5jKwT+8lLO07R2fOUXcjFnl2CMc+OT5E/xpJSlvy/OLKKs8gpyWHdT9cR4RvR0nsghBBCCNGmlFrtFJXaKCyzUlhqo6jURlGZzbnsqCuvL3T2Kyovl5XXWyvaS22UWG3s/t3/cTnc3NbYSSkGAPvKC1rrAqXUEWd9atV25/KARm5TtEFlJQVkZx0kO/sI2bknyM4/TbblDIXnL2A9n4POLsKcY8UnDwLzFCF5EJ4HRg1eOD4VFFAIwJFT+yjwdfwhZhVnSUIlhBBCiHbHbtcUW6skNZUSmKJKSU1FguOeEBVVqnckSxXJk9WumyTu4jI7vl7GJhnbkxqbUAUAF6rU5QKBldpzq7QFKKWU1trtN6+Umo3jFkG6devWyLBEa1BiK+FA1gGyCjPIzthPdt5JsgvOkV2cRXZpLjllhRQWlhCYYSMgFzrkQmSuJtL5c1AemOyVR6z/NxTh+VDg61jOKc7xxG4JIYQQQlSjtcZm19g12LXGXl62O8o2Z52rbHeWNdjs2rF+lfZSm73GqzzVrv7UcpWnuKwiAbocFZZa20VCZQGCqtQFAfm1tAcBlqrJFIDW+g3gDYD4+PimSXNFsziUfYg1qR+x7shn5NmKXPXepZoe56DXGc2QM5qeZzTR2Y3bVl6AIi/ERGGYL6WRgdijO2Ds0g2f7j3x69KdPwV1JMQ7hFCfUAK9Ai89oBBCCHEZsTtPyt1OxLVG28HmPCl3najrSv2rnMzbq6xvt1c/0XcsV+pfdX2tsTnHcK3vWq60viuxqNS/hvVtuvr2XAmLa78r9a8hmSlf3zGWe5/qY1Xqf9HfG9X2x9ZEV2guFyaDwtfLiJ+XEV+zEV8vE36Vyn5e7nU+zjpXvauPET9nPx+zkVA/r0tvvBVobEKVguM5KQCUUv5AL2d9eftgYJezPLhSm2hDCssK+c+Bj1n94wd8X3Qas1XT/RyMOKvp5UyeumSCoZ7/3pQGmCmLCERHRWDq2hWf7r0J6N6TkJg++HeNweDj0zQ7JIQQokZa135CWfMJa5UTfedJc00n5eXr13QibbfjfqJf6QS5phP1ep3o1+lEuuYT9Won+jWc+F8smakxIapHYiBEXfmYDfh5mVwJTuUExpEEVaozmyq1V06CqtSbTfh6GfEyte83MdUpoVJKmZx9jYBRKeUDWIF/Ac8rpW4HvgD+AHzvnJACYAXwG6XUvwENPAr8w7O7IFpKYVE23+7/mC2H1/PD0cN0Omtn2FnN3Wc03S5UvV2vFkYDph7d8O3VB3OXLpg7d8arc2fHcnQ0Bl/fJt8PIYSojc2uKbPZKbXZsdqcy1Y7Vnv15TKrnTK7dvy01b5stWtKKy1bbZUTC+eJe00n+lW/Ya/8Df3FrgBUOlGvOQm61BUA98Sg+j0mQohyRoPCoMCglHPZWTYojEqhlMJooNKye7tBKQzOOqNBYTYaGn2Vp3x9g6H1T+5wuarrFarfAwsrle8BFmmtn3YmUy8B/wR2AndV6vc60BP4wVl+y1knLkO2wmxSf/iYH3Z+zvnDxzGds9L9rObWC3B7XZIngwHvXr3wGTgQn4ED8B00CO++fTF4ezd57EKIlqW1diUPpTZnMmGrIVmpZdnxqb5stdkprbRc5hy/puXq41TUVY2rfFmSB9HaqfKTd6VQzpPwysvlJ+jGSif25Sf7SjlO7B3LjhN9g6pIAioSgkr9qyQGFQmAwlh5fKUwGKomFhXj1jexUM4xjQZcy4ZK8VZOXi6+31Qaq8rvqXIcNaxfERe1/g5F+6RqeJypxcXHx+s9e/a0dBjtnj3jJGmb3+Pw7kTyjp3B75yN6AzHzHp14RUTg8/AgfgOGuhIovr1w+Dv37RBC9EOaK0ps2msdjtlVmfSUNtyrclE3RIUTyYrou0oP4ksP7m8+Alw1ZPWKifZVU6qa0oGDOXbu9j4NcRQcZJMLWPVklhUOgF3W7/WxKLKSXu1ZKKGxOKS+02lsWpOcAwKOYkXopkopb7VWsfX1NbYZ6hEG1T2w/848sgvsZ6xYtSKTkCnS6xj7xBI8OBr8B18pSN5GjAAY6BMAiFEOZtdk19cRm5RGXlFVvKKy8grcpaLK+oc7WXkFVvJLSqjoMTqurXMlUTZWt8XYW2VUuBlNGA2GjAblfOnY9nkXPZyLbu3X2zZVGk9U5Vv1uueWFT5tr5yclPTt/V1OPGvfsWi0rf1chIvhBA1koRKVHP6iXno0zaMtUxTnhtuwtAzmo7xY4gYdgM+/WMxBlWd7FGItkVrjaXESl6xtSIRciY+VROjiuUy8p2JkaXE2tK70OIczwPULdFwSzoMBswmZx+DAbOp/mOUL5uNCq86LJePaZRnDoQQQlyCJFTCjeWTdyg8UjHVeXo45Ef74tuvNz1G3UqvUbdgkitP4jJWWGrlQn4JGZaSaslPXrGV3EJnuYarSa1pRi1zeaJhVHiZDM6kw5kI1LJsMhjwqmXZbHImK8bKy8qRyNRxDNdyTQmKwSAPRAshhGiTJKESLlpr0v72d1f5v4MVwQt/y4z+M1owKiEurcRqI8NSyoX8Etcnw1Jl2VkuLG2ZlxsqBYHeJoJ8zQT5mAn2NRPkayLIx0yQr7PsY6pYdvYL8DHhZTQ4bjszVSRRcsuVEEII0TpIQiVcCvf9G32+BIAyI+ye0I23+01r4ahEe2W12ckscCZJVZOjKst5xc1zO52fl9GZAJmcCVD1ZCjIVW+qlDiZCfQ2yRUaIYQQog2ShEoAYLeV8detv+d2Z7nYC56Y9FdMBvlPRLhzzfBmtVNis7mWyydOKKm07KovL9vsbn1LneuWWu1kF5a6JUpZhaVNMmW1l8lAZIA3EQFeBPt5uSdDVa4cVVwpcrSbje37xYVCCCGEqE7OlgXplnQW/mcOZ4rLXAmVr8mX2PDYFo2rPbPbHdNMV01AymyOhKV8xje3OueLRUur9K05uXGOb61IiKolPLWM05qeIypnNCgiAryICPAmMtDbkTA5f0YGOj7lbUE+JrldTgghhBAeIwlVO6a1ZvWh1SzZvYSi0gIW/afi2ZKgAYNbMLLmUf4un6qJSGnVBKRKW01XV2rsd7HkxlbxnqCartxYW2PW0syUgjA/r2oJkSNZ8iIywMdZ70Won5fcTieEEEKIFiEJVTt1KPsQz+5+lp1ndgIweZemX5qz0Wik42MLWi64JvDN0Uz+sj6VE5kFFVdk5EWjDWJQjtvmzEYD3s6fXiaD61095cteJoOzn8LLZHTWKfd+lcYJ9jW7EqcOgd6E+XthklvshBBCCNHKSULVzmQUZfBy0st8cugT7NqRUAxLtTP9Kxs43zsV8Ytf4NO/fwtG6VlbDl5g9oo9lFgvrwSqelJS/hLRKolMlaSmfL2KOlVzv6pJUKXkpnJfs1HhbTS6liXJEUIIIYSoIAlVO1FsLea9H9/jrR/eotBa6KofcBLmfWrD4EymfPr1JOLB2S0VpkfZ7JqNP57j4ZV7a70aZTKo2hMN59WUGhMNo+NFoxdPbpyJjNHoSogqb6PmbZYnOTItthBCCCHE5UASqnbAZrcxZ9Mcvj33rVv9zV5D+Pnqb8HuOHH3CvOi67v/RJnNLRFmo5RYbRw8ayHldC4pp/NIOZ1L6tl8t3cOdQ7x5Z2Ea+gU4uO6kiPP3QghhBBCiMaQhKod+OLYF27JVO+Q3jw6aDY9fv078kockx8YfWx0/dufMYWGtlSYdZZfXMaPp/OciZMjeTp83nLRiRy6hvnywf3D6Rrm14yRCiGEEEKItk4Sqjau1FbKy3tfdpWn9ZvGYzFTKPnrPZw8UOSq7/LL8XgNndgSIV6S1prP9p3mPylnSTmdx4nMwkuv5BQR4M3IXuH8dkI/OgX7NmGUQgghhBCiPZKEqg2zlFqYv2U+pwtOAxDqFcRDB1LI+ONrZB/2o3wSiqChffCb/WILRlo7u12zaF0Ky3ecuGTf7uF+DIgOYkB0MP2jgxgQHUSHQJ9miFIIIYQQQrRXklC1UWcsZ/jlf3/J4ZzDABhtmj/85zRnv83AXubv6mfw9abD4ldaKsyLstrsPLbmez75Lt2t3mRQXNEx0Jk8BdG/UxCx0UEE+Vx+z34JIYQQQojLmyRUbdDZgrPc/e+7ySjKACA8V7PkgzL8c3yoPNed3+BYOj7zF8ydO7dMoJewYPX3/GtvRTI1Ma4TvxjTiys6BuBtMrZgZEIIIYQQQjhIQtUGfbD/A1cyZdKaP20rwD+n4tY3c3QHOj75FAE33NBqp+Y+esHilkxNG9qNP04ZiFFm5RNCCCGEEK2IJFRt0NHco67l32Vm0eWMD+XTOITffx+RDz+M8vJqmeDq6FhGgWv56u6h/PmnA1tt8ieEEEIIIdovQ0sHIDwvLf+Uazn2KBRlVSRPoTNmtPpkqtRq598/nHWVe0cGSDIlhBBCCCFaJblC1cZ8n76DI7lH8SvWzNpkx/BDMOVvZ/Lq3h1zx44tGt+lHD6fzyMfJpFyOs9V16uD/0XWEEIIIYQQouVIQtWG6Nx0Xv98DqMP25n+lZ0wS0WbMTSUTov/3HLBXYLWmn9+c4I/frGfEmvF1BnX9o5g2tBuLRiZEEIIIYQQtWt0QqWUslSp8gVe0VrPVUrFAMeAgkrtz2qtn2nsdoW74/u2sXXpr5j5QylBRe5tgeNvJuqppzCFhbVMcJdgt2v+3+c/smz7cVedl8nA4zf3496RMRhkIgohhBBCCNFKNTqh0loHlC8rpfyBc8DHVbqFaK2tjd2WcGcvLiZnw3oOLXuZoNR04qu0G8PDifrDHwi6aVyLxFcXJVYbv1m1jy++P+Oq6xcVyN/uGkK/qKAWjEwIIYQQQohL8/Qtf3cA54GtHh5XVGIvLSVr2XIy3nwDnW+hatpRFAjRM+cQcc9MTKGhLRJjXVhtdu5fvoethzJcdRMGRfHCnUPwMct7poQQQgghROvn6YRqFrBCa62r1J9QSmlgI7BAa51RfVVRF5atWzn3xz9ReuKEW73VAAd6avp3z+XKa8eipj3SQhHW3YaUs27J1KwR3fnD5AHyrikhhBBCCHHZ8FhCpZTqBowB7qtUnQFcAyQB4cDLwPvATTWsPxuYDdCtm0xCUFXpqVOcW/wXLJs3u9WfC4HNQ4z07prLtLI8zADxCS0QYf3985uKpHDWiO48fcsAmR5dCCGEEEJcVjx5hWomsE1rfay8QmttAfY4i+eUUr8CziilgrTWeZVX1lq/AbwBEB8fX/UKV7tWlJTEyZ/fh72w0FVX4A0fXWdg+7AAXriQybBc569zyHS44sYWirTujlyw8M3RLACMBsUvftJbkikhhBBCCHHZ8XRC9ZdL9ClPlOTMuY5KDh3i5INz3JKpzXGKD35iwCeiI+/E3kffNb9wNAR2gpsXt1Ck9fPN0UzX8vV9I4kK9mnBaIQQQgghhGgYjyRUSqmRQGeqzO6nlBoG5ACHgFDg70Ci1jrXE9tt68rS0zl53/3Ycx2/LmNoKP+YFkhi8GkA/jTiKfqeO1axQo8x4BPcEqHW26FzFbPtX9mt9U6cIYQQQgghxMUYPDTOLOATrXV+lfqewAYgH0gGSoBpHtpmm3f+xRexnj8PgMHPD79/LHYlU95Gb4Z1GgYXDlSsEN6rJcKsF5td89bWo3y4+6Sr7ooOARdZQwghhBBCiNbLI1eotNYP1lK/EljpiW20RyX797uWo597llM9IiHFUe4R3ANfow8c+HfFCp0GN3OE9XP4fD4LVn/P3pM5rjo/LyNXd5crVEIIIYQQ4vLk6WnThQdZL1RMKe47ZAhl9nRX2cvgBae/gxznTHnewdDzJ80c4cVprTl4zsLXhzPYfiST/x28QKnN7mrvFxXIkp8NJjzAuwWjFEIIIYQQouEkoWql7MXF2HJyHAWDAWNoKOdOfedq9zZ6w+Y/VazQbyKYWjYx0VpzIrOQ7Ucy2X4kgx1HMsksKK3Wz2xUzB17BXPG9MLL5Km7ToUQQgghhGh+klC1UqUnK54xMnfujDIa2XB8g6tuCN5w5L/OkoKhDzRzhA55xWVs+vEc249ksuNIJuk5RRftf2W3EP5yWxx9owKbKUIhhBBCCCGajiRUrVTp8eOuZa8eMRSUFbDl1BZX3aSULys6D3sQOl/VjNE5FJRYmfj3rZzKqj2JCvUzM7JXBCN6hTOyVzg9IvzlfVNCCCGEEKLNkISqlbJeuOBaNkdHsz9zP6V2x+1zvZUPPfOdV7CCusDY37dEiLy25Ui1ZMrfy8iwno7kaWSvCPpFBWIwSAIlhBBCCCHaJkmoWiur1bWovLw4kXfCVe6bX/FSXG55Ebyb//a59Jwi3vjfUVd51oju3HplZwZ1DsZslOeihBBCCCFE+yAJVStVduasa1kZTRzOOewqdy8rcyxcfS/0/r/mDo3iMhuPrNxLidUxY9/AzkEsnDxArkQJIYQQQoh2Ry4ltEKlx4+T/cEHrrLPgP58deorV3lASSn4BMO4Z5o9Nptd8+uPkthzIttV99TE/pJMCSGEEEKIdkkSqlZG2+2ceXoRutTxvJTPwIEcvboT6RbHO6gCbXZGFBVD5/hmv9Uvt6iMBR/vY31yxdWz30+MZVjP8GaNQwghhBBCiNZCbvlrZS789W8UfvONo2Aw4PXEwzyxdYGr/YbCQszQrLP6aa354oczLFr3IxfyS1z1Px/Vg/tH92y2OIQQQgghhGhtJKFqRXI/+4zMN990lQNmTuVXBx7ntDUfAB+7nZm5+eDfAYbc3Swxncoq5A+fJvPVgQtu9bdd2ZnfT4xtlhiEEEII0XTsdjtpaWkUFBS0dChCtCh/f3+6dOmCwVC/m/gkoWolSo4c4czvn3KVA0Zczd9DP+awc7I/k9a8cCGLK668F65/EnxDmywWrTXfnsjm3e3H2ZB8Fptdu9o6BHqz6JYB3DwwSt4nJYQQQrQBGRkZKKXo27dvvU8khWgr7HY76enpZGRk0KFDh3qtKwlVK3H++SWu56a8unSgU7fNfOUVSvljbs/oMEbPXAlRg5oshhKrjS++P8O7Xx/nh/Rctzal4J5h3Vlwc1+CfMxNFoMQQgghmldOTg4xMTGSTIl2zWAw0LFjR06cOCEJ1eWoYOcuLImJjoKCzgOSsRjtlBgckz34G7yYNHOLI6tpAlpr3vn6OK8mHiHDUlKtfXjPMB67uR9XdWu6q2JCCCGEaBk2mw2zWb4sFcJsNmOt9C7YupKEqhXIevdd13Jw90IKw+38IbJi5ryOgV2aLJkCWLb9OM98/qNbnbfJwJQhnUkYFUNsp6Am27YQQgghWp7cxi9Ew/8OJKFqBUrTTrmWTw4q44nOncg0GV1110Rd02TbPplZyHMbDrjKUUE+zBjRnWlDuxHm79Vk2xVCCCGEEKItkJtlWwF7bp5r+fFuYW7J1D2x97DgmgU1rdZoWmue/NcPFJXZAOjbMZAtj/2Eh67vLcmUEEIIIVqFbdu2MXLkSIKDgwkLC2PUqFHs3r3b1V5QUEBAQAATJkyotm5MTAy+vr4EBAQQFRVFQkICFoulOcMX7YAkVC1Ma40tryKhKvBx/Az3CefV/3uVx4c+jrfRu0m2vf9MPtsOZwBgUPDcHXF4V0rmhBBCCCFaUl5eHpMmTWLu3LlkZWWRnp7OwoUL8fauODdavXo13t7efPnll5w5c6baGOvWrcNisZCUlMTevXtZvHhxc+6CaAckoWphurAQXeKYCKLUBCVmiNVm1tyyhms7X9uk29744znX8vhBnRjcNaRJtyeEEEIIUR8HDx4EYNq0aRiNRnx9fRk3bhxxcXGuPsuXL2fOnDnExcXx/vvv1zpWVFQUN910E0lJSU0dtmhnJKFqQbacHM6/+KKrnOcHKMWjPj0I9w2vfUUPsJRYWZuU7irfPCCqSbcnhBBCCFFfffr0wWg0MmvWLNavX092drZb+8mTJ0lMTGT69OlMnz6dFStW1DpWWloa69evp3fv3k0dtmhnZFKKFmAvLCTrvX+S+dZb2PPzXfVnQhX9S0oY2i2+SbdfYrXx4Ht7OJbheCO6j9nAmL6RTbpNIYQQQlweYp74otm2dfwvEy/aHhQUxLZt23j22Wd54IEHOHv2LBMmTODNN9+kY8eOrFixgri4OPr3709ISAiPPfYYe/fu5corr3SNMWXKFJRSWCwWxo4dy6JFi5p6t0Q7I1eomlnu519w+KabuPDXv7olU0ei4J1xBnrihRo6u8m2n55TxNwP9vL14UxX3aJbBsjLeoUQQgjRKsXGxrJs2TLS0tJITk7m9OnTzJs3D4AVK1Ywffp0AKKjoxkzZgzLly93W3/t2rXk5+eTmJhIamoqGRkZzb0Loo3zSEKllEpUShUrpSzOz4FKbTcopVKVUoVKqa+UUt09sc3LUVl6OqcXLMB2oeIP2atjCKk3FvHbBCPpEYrQK24Gf8/e7meza746cJ77l+9m9LOb+bLSs1Pzx/Vh6jXdPLo9IYQQQoim0K9fPxISEkhOTmb79u0cOnSIxYsXExUVRVRUFDt37mTlypU1vpx1zJgxJCQkMH/+/BaIXLRlnrzl71da67cqVyilIoBPgPuBdcAzwEfAcA9u97JR8M1O0BoAY3g4kY88TEjeO7xYet714t6Y6KEe216GpYRVe07xwc6TpGUXVWtPGBnDQ9fLfcRCCCGEqHCp2/CaU2pqKl988QVTp06lS5cunDp1ipUrVzJ8+HCWL1/OjTfe6PbcVFFREXFxcaxfv57JkydXG2/evHnExMSQlJTEkCFDmnFPRFvW1M9Q3QakaK0/BlBKPQ1kKKX6aa1Tm3jbrU7hnj2u5bCEWYROvpGC53/Bd106uurjoxr//FRxmY2lXx5g2fbjlNl0tfZre0cwc0R3buzfUd6MLoQQQohWKzAwkJ07d/LCCy+Qk5NDSEgIkyZN4vnnn6d79+6sWLGCqCj3ibVmzJjB8uXLa0yoIiMjmTlzJs888wxr1qxprt0QbZwnE6rFSqm/AAeA32mtE4EBwL7yDlrrAqXUEWd9u0qobHl55G/a5Cr7xcdz5tAGfhUVQa7R8e6nMJ8wegT1aNR2vj2RzYKP93HUOeFEuRA/Mz+7ugt3D+tOjwj/Rm1DCCGEEKI5dO7cmVWrVtXYVnXGv3KvvPKKa/n48ePV2l999VWPxCZEOU8lVI8DPwKlwF3AOqXUECAAuFClby4QWHUApdRsYDZAt25t75merOUrXJNQeHXvztFoA7/a+Fcyvb1cfR4a8lCDrxgVl9n468aDvLn1KPZKF6UGdw1h5vDuTIzrhI9ZXtorhBBCCCGEJ3kkodJa76xUXK6UmgZMACxAUJXuQUB+lTq01m8AbwDEx8dXv0/tMmYvKCCr0v29frMTmLn5l+TZHM81mbRmYedxTOl7Z4O38eQnP/DJ3or3SgV4m/jdxFjuuqar3NYnhBBCCCFEE2mqZ6g0oIAUYFZ5pVLKH+jlrG83CnbudF2dMnfvxg9DQsjblgdAkM3G385ncM21Exq1jR1HK6ZBv7Z3BM/eEUfnEN9GjSmEEEIIIYS4uEZPm66UClFK3aSU8lFKmZRS04HrgP8A/wIGKqVuV0r5AH8Avm9vE1IUbN/hWg668UZOFqa5yrdYCrimuASMXjWtWmeFpTbX8j+mXSnJlBBCCCGEEM3AE++hMgN/xPGsVAYwF5iitT6gtb4A3A78CcgGhuF4xqpdKdy927XsN2IEp/JPucpdyt+ToBr+fNOWgxfIKy5zlX295FkpIYQQQgghmkOjb/lzJk3XXKR9E9Cvsdu5nNlyc13L3j16cPC7g65yz1JnIuQX1qCxtx/OYPaKPeWvt2JoTJhMPiGEEEIIIUQz8cQVKnEJ2lbxtu6zpRkcyjnkKvcrT6iCu9Z73JTTudy3fA8lVjsAnUN8eWHq4MYFK4QQQgghhKgzSaiagS4ucS3P2vhzSmyOcpTVSqjdDgEdwTug3uO+/NVhisocz05FBfmw8oHhdAn180zQQgghhBBCiEuShKqJ2SwFrhn+yoyQYS4GwKA1D+Q4ZvpjzGP1Hje3qIxN+8+7ym/NiqdbuCRTQgghhBCXo61bt9K3b19XOSYmhk2bNrVgRKKuJKFqYpnHKyY0zAgCrRS9S0t578w57iyywm1vwjX313vcL1POUuq81W9AdBADOwd7LGYhhBBCiNZg27ZtjBw5kuDgYMLCwhg1ahS7K032BVBQUEBAQAATJlR/BU1MTAy+vr4EBAQQFRVFQkICFoulQbHMmTOHgIAAAgIC8PLywmw2u8rjx4/n+PHjKKVcdeWfjz76yDXGrl27mDBhAiEhIYSFhTF06FDeffddAEaPHs2BAwcaFJtoWZJQNbG9Kf91LWcFKn6ZncOq9LPEeUfCz/8DcQ17me/h8xX/GNzYv2Oj4xRCCCGEaE3y8vKYNGkSc+fOJSsri/T0dBYuXIi3t7dbv9WrV+Pt7c2XX37JmTNnqo2zbt06LBYLSUlJ7N27l8WLFzcontdeew2LxYLFYuHJJ59k6tSprvL69etd/XJyclz1FouFqVOnArBjxw7Gjh3LmDFjOHz4MJmZmbz66qtu6zaU1Wq9dCfRZCShamJpJ5Ndy5Hh4fwiJw8zQL9JED2kweOWT0QBEOxrbniAQgghhBCt0MGDjlmRp02bhtFoxNfXl3HjxhEXF+fWb/ny5cyZM4e4uDjef//9WseLioripptuIikpqSnDrtWCBQuYNWsWjz/+OBERESiluPrqq1m1ahUAiYmJdOnSxW2dpKQk4uLiCA4OZurUqRQXF7v1ffbZZ4mKiuLee+9t9v0RFSShakJaay6kV8zoFxYSUtFoaNzU5oWlFd9EeJtkmnQhhBBCtC19+vTBaDQya9Ys1q9fT3Z2drU+J0+eJDExkenTpzN9+nRWrFhR63hpaWmsX7+e3r17N2XYNSosLGTHjh3ccccd9Vpv1apVbNiwgWPHjvH999+zbNkyV9vZs2fJysrixIkTvPHGGx6OWNRHo99DJWp3MPsghqyKd1B1CPSpaPQJadTYRy4UuJY7BftcpKcQQgghRD083YzPZT+dW2tTUFAQ27Zt49lnn+WBBx7g7NmzTJgwgTfffJOOHR2PO6xYsYK4uDj69+9PSEgIjz32GHv37uXKK690jTNlyhSUUlgsFsaOHcuiRYuadJciIiLcyjt27CAoKAi73U6nTp3qNdbDDz9MdHQ0AJMnT3a7umYwGFi0aFG1WyBF85MrVE1o86nNBFfkPXibSysKQfX7g6rMZtf8eDrPVR7QOajBYwkhhBBCtFaxsbEsW7aMtLQ0kpOTOX36NPPmzXO1r1ixgunTpwMQHR3NmDFjWL58udsYa9euJT8/n8TERFJTU8nIyKhxW1u3bnVNJDFgwIAGx5yRkUFOTo7rExsbS2hoKAaDocZnvC4mKirKtezn5+c2oUZkZCQ+PvKlemsgCVUT2ncuiZ5ntats0hXTnBMU3eBxz+QWud4/Fe7v5X7lSwghhBCiDerXrx8JCQkkJzueT9++fTuHDh1i8eLFREVFERUVxc6dO1m5cmWNkzSMGTOGhIQE5s+fX+P4o0ePdk0kkZKS4tHY/fz8GDFiBGvWrPHYmEopj40lGkdu+WtCoQfP0SXTsax9vfFT+x0FZYTOVzd43MoTUgTJhBRCCCGE8KSL3IbXnFJTU/niiy+YOnUqXbp04dSpU6xcuZLhw4cDjskobrzxRrfnpoqKioiLi2P9+vVMnjy52pjz5s0jJiaGpKQkhgwZ0ly7AsBzzz3HuHHj6N69Oz//+c8JDw9n3759LF68mA8//LBZYxGeJVeomtCAb85VFK6JwWB2Xq3qOgx8Qxs8bpmtIqEyGeTbCSGEEEK0PYGBgezcuZNhw4bh7+/P8OHDGThwIEuXLqW4uJhVq1Yxd+5c19WpqKgoevTowYwZM6rd9lcuMjKSmTNn8swzzzRZ3CEhIW7voXrhhRcAGDlyJJs3b2bz5s307NmTsLAwZs+eXeP7s8TlRWmtL92rmcXHx+s9e/a0dBiNctpymtRxN9Apy1EOmhFF57LvHIX/exqu/XWDx35vx3Ge+tRxKXp4zzA+nD2ikdEKIYQQor3av38/sbGxLR1Gu7d582buv/9+jh492tKhtGu1/T0opb7VWsfXtI5coWoiq374Jx2ds3vaFXQqciZTygCDGvYy33Krv01zLU8c1PDJLYQQQgghROuQnJxMjx49WjoM0QDyDFUTSDqfxNavV3KT8+KfPdwPQ/lvuveNENy5wWOnns1jX5rj3mYvo4HJgxs+uYUQQgghhGh5jzzyCJ999lmttyqK1k2uUHnYsdxj/Grzrwi/UDFFenDHwIoOPcc0avz3dpxwLY8b0JEQP69GjSeEEEIIIVrWiy++yLFjx7juuutaOhTRAJJQeVBhWSG/2PQLckty6VTpZd7eoZVm4vNp+Mvy8ovL+NfedFd5xvDuDR5LCCGEEEII0XiSUHnQ5lObSbc4Ep5OeUZXvVfZ4YpOfhFVV6uzxAMXKCx1vH+qT8cAhvYIa/BYQgghhBBCiMaThMqDUjIqXgI3gIrJIkzmQsdCSPdG3fL3zdFM1/LNAzvJC92EEEIIIYRoYZJQeVByRrJrObigot7kYwOTD0z9J5h9GzS21podlRKq4XJ1SgghhBBCiBYns/x5SLG1mJRMxxUq/yKN+VjFs05mfxvc/Bx0imvw+N+n5XL0giNL8zEbuLJbw18MLIQQQgghhPAMuULlId9f+J4yexkAE0+Gg9XxrJNPWCkmXzt0Hdao8VfuOulanjgoGl8v40V6CyGEEEKIy8nWrVvp27evqxwTE8OmTZtaMCJRV41OqJRS3kqpt5VSJ5RS+UqpvUqp8c62GKWUVkpZKn2eanzYrc+3577Ft0Qza5ON2/91wVUf2LXYseAV0OCxi8tsfLbvtKt897BuDR5LCCGEEOJysW3bNkaOHElwcDBhYWGMGjWK3bt3u/UpKCggICCACRMmVFs/JiYGX19fAgICiIqKIiEhAYvF0qBY5syZQ0BAAAEBAXh5eWE2m13l8ePHc/z4cZRSrrryz0cffeQaY9euXUyYMIGQkBDCwsIYOnQo7777LgCjR4/mwIEDDYqtod5++2369etHYGAgHTt2ZOLEieTn5wOQkJCAl5eX274MHjyYrVu3usr+/v7V9vnkyZOX2Grb44krVCbgFDAGCAaeAlYppWIq9QnRWgc4P894YJutitaasi8T+esbNibu1ijteKOv0dtGcEwh+HeA4C4NHn/38SzX7H4x4X5c1S3EE2ELIYQQQrRaeXl5TJo0iblz55KVlUV6ejoLFy7E29vbrd/q1avx9vbmyy+/5MyZM9XGWbduHRaLhaSkJPbu3cvixYsbFM9rr72GxWLBYrHw5JNPMnXqVFd5/fr1rn45OTmueovFwtSpUwHYsWMHY8eOZcyYMRw+fJjMzExeffVVt3Ubymq11nudLVu28OSTT7Jy5Ury8/PZv38/d955p1ufxx57zG1f9u3bx+jRo13llJSUavvcrVv7++K/0QmV1rpAa/201vq41tqutf4cOAZc3fjwLg9nnniCG9/5nrBKX3j4dVZ0vyEDs68dRv4KDA2/RS/xQMUVr5/07SCz+wkhhBCizTt48CAA06ZNw2g04uvry7hx44iLc38mffny5cyZM4e4uDjef//9WseLioripptuIikpqSnDrtWCBQuYNWsWjz/+OBERESiluPrqq1m1ahUAiYmJdOni/gV8UlIScXFxBAcHM3XqVIqLi936Pvvss0RFRXHvvffWO57du3czYsQIrrzySgDCwsKYNWsWgYGBjdzT9sfjz1AppToCfYCUStUnlFJpSql3lVI1vohJKTVbKbVHKbXnwoULNXVplawZGeR++pmrnO0PHZ6aR7dr0/EOsoF3MMT/vMHjn8oqdHt+akyfyEbFK4QQQghxOejTpw9Go5FZs2axfv16srOzq/U5efIkiYmJTJ8+nenTp7NixYpax0tLS2P9+vX07t27KcOuUWFhITt27OCOO+6o13qrVq1iw4YNHDt2jO+//55ly5a52s6ePUtWVhYnTpzgjTfeqHdMw4YN4z//+Q8LFy7k66+/pqSkpN5jCAePzvKnlDID7wPLtdapSqkA4BogCQgHXna231R1Xa31G8AbAPHx8dqTcTWl0kr3iZ4OhTcfuYIPR10Jh5yVEb3Bu2GZvtaaJ//1g+t2vys6BHDtFQ1/MbAQQgghxKUMWj6o2bb1w6wfam0LCgpi27ZtPPvsszzwwAOcPXuWCRMm8Oabb9KxY0cAVqxYQVxcHP379yckJITHHnuMvXv3uq66AEyZMgWlFBaLhbFjx7Jo0aIm3aeICPdztR07dhAUFITdbqdTp061rFWzhx9+mOjoaAAmT57sdnXNYDCwaNGiardA1tXo0aP55JNPeOWVV3jxxRexWq3Mnj2b559/HqPRcWfVkiVLeOmll1zr3HrrrSxfvrxB22vLPHaFSillAN4DSoFfAWitLVrrPVprq9b6nLN+nFIqyFPbbWllp065lk90VHSPjoWinIoODZyMQmvNc/85wNZDGQAYFDx3Rxxmo0zMKIQQQoj2ITY2lmXLlpGWlkZycjKnT59m3rx5rvYVK1Ywffp0AKKjoxkzZky1E/61a9eSn59PYmIiqampZGRk1LitypMtDBgwoMExZ2RkkJOT4/rExsYSGhqKwWCo8Rmvi4mKinIt+/n5uU2oERkZiY+PT63rDhgwwLU/W7durbHP+PHjWbduHVlZWXz66acsW7aMt956y9U+f/58t32RZKpmHjk7V46Het4GOgK3a63LaulafuWpzTwEVHbuvGs5IwhigmMg5ZOKDsFd6z2mza753dpkXk084qr7+age8u4pIYQQQrRb/fr1IyEhgeTkZAC2b9/OoUOHWLx4MVFRUURFRbFz505WrlxZ4yQNY8aMISEhgfnz59c4fk2TLXiKn58fI0aMYM2aNR4b81LP1KekpLj2Z/To0RftazAYuOGGGxg7dqzr9yvqzlO3/L0KxAL/p7UuKq9USg0DcnDcABcK/B1I1Frnemi7Lc6Wk+NazvdVdLbZ4YePKzpcnVCv8aw2O79etY91laZJv6FfB+bf1PciawkhhBBCeMbFbsNrTqmpqXzxxRdMnTqVLl26cOrUKVauXMnw4cMBx2QUN954o9tzU0VFRcTFxbF+/XomT55cbcx58+YRExNDUlISQ4YMaa5dAeC5555j3LhxdO/enZ///OeEh4ezb98+Fi9ezIcfftissQB8+umnFBUVcdNNNxESEsLu3bvZsmULf/vb35o9lsudJ95D1R14EBgCnK30vqnpQE9gA5APJAMlwLTGbrM1sWZUTKCR7wuBP64Du/Nbke6joOs1dR7Lbtc8tuZ7t2RqypBoXptxNT5meZGvEEIIIdqPwMBAdu7cybBhw/D392f48OEMHDiQpUuXUlxczKpVq5g7d67r6lRUVBQ9evRgxowZtd6aFhkZycyZM3nmmaZ7i09ISIjbe5leeOEFAEaOHMnmzZvZvHkzPXv2JCwsjNmzZ9f4/qzmEBoayptvvskVV1xBUFAQ99xzDwsWLHDdQgmOJLDyvlR9Pkw4KK1b3/wP8fHxes+ePS0dxiVprTn8k+uxnjsHwO9nGHnKcI4hJaWODjP+Bb3G1nmsZz7fzztfH3PVzRjenUW3DMBgaDN3SAohhBCildm/fz+xsbEtHUa7t3nzZu6//36OHj3a0qG0a7X9PSilvtVax9e0jsxw0AglBw+5kimLD5zspOlfnkwNf6jOyRTAP7854ZZM3XVNV/7frZJMCSGEEEK0B8nJyfTo0aOlwxAN4NFp09ub4uSKe4yTuytiy0rxAugyFG6s35Sc73593LU8fmAUf/rpIHmBrxBCCCFEO/DII4/w2WefySx6lym5QtUIttw813JmEAwpLgWfEPjZu2A013mcghIrxzILADAaFEvvHIxRrkwJIYQQQrQLL774IseOHeO6665r6VBEA0hC1Qj2Su8CKPSCQSUlcOvLENylXuMkHrhA+aNsPSP88fOSC4dCCCGEEEJcDiShaoTzGSdcyyVecGW/2yB2Ur3G+CEtlwWr97nK1/QI81h8QgghhBBCiKYlCVUj7D/xjWu5F1YibqjfFJxp2YXcu2w3haU2ADqH+DLvhis8GqMQQgghhBCi6UhC1UDnCs6Sk5flKo+O7Af+4fUa49kNB8iwlAAQ7Gtm+c+voUOQj0fjFEIIIYQQQjQdSagaKG3vciJzK97hFT1qZr3Wzy0s4z/JZ13lV++5it4dAj0WnxBCCCGEEKLpSULVECX5XPj2LaKyK6q8Bo2s1xCf7Uun1GYHIK5LMCN7yZunhRBCCCHaq61bt9K3b19XOSYmhk2bNrVgRKKuJKFqiJM7KUovJajIUbR6Gfn/7d15dJX1ncfx9zc3hCULSQzmCohAgULRAJWRbQItjFDAFEbEqJGCHqVUC2UKqGVoKTJzEESn6SBaKBVwGJQBXClqWwwFiywtwYayiaxKkLBIQhIM5Dd/3HglCQnkmuQJ8Hmdcw/32b9P8iP3+d7fFp6YeNmHF50vZv66rybxHfbtqo0KKCIiInItWL9+PT179qRx48bEx8fTq1cvNm/eXGqfM2fOEBUVxaBBg8od37JlSxo2bEhUVBR+v59Ro0aRd8EozVUxZswYoqKiiIqKIiIignr16gWXBw4cyP79+zGz4LovX6+88krwHJs2bWLQoEHExsYSHx/PbbfdxosvvghAcnIyu3btCim2UC1YsID27dsTHR1NYmIigwcPJjc3F4BRo0YRERFR6l46derEunXrgsuRkZHl7vngwYNVjqOiawHlfq4tW7bkqaeeKneOjIwMzIxZs2aV22ZmREZGBs/x0EMPVTnGyiihCsXZ0zTZGhFc/Kx7G8znu+zDl205xMET+QDENqrHv367WbWHKCIiInIlO336NHfccQdjx47lxIkTfPLJJ0ydOpX69euX2m/58uXUr1+fd999lyNHjpQ7z5tvvkleXh6ZmZls3bqVGTNmhBTPCy+8QF5eHnl5eUyePJnU1NTg8urVq4P7nTp1Krg+Ly+P1NRUADZs2EDfvn3p06cPH330EcePH+f5558vdWyozp07V+Vj1q5dy+TJk1m6dCm5ubns2LGDu+++u9Q+jz32WKl72bZtG8nJycHl7du3l7vnFi1ahHQPF7vWhb68xvLly5k+fTp/+MMfSm1ftGgR8fHxFU6OvG3btuC5f/vb34YUY0WUUIWg8KP9XHcoMFdUMXDsrktPwlZ0vpjf//0IIxZs5N9fzQqu/2HvbxDT4PInARYRERG5FuzevRuAe++9F5/PR8OGDenfvz9JSUml9lu0aBFjxowhKSmJJUuWVHg+v9/PgAEDyMzMrMmwKzRp0iRGjhzJ448/TkJCAmbGrbfeyrJly4BADUvz5qVbLWVmZpKUlETjxo1JTU2lsLCw1L4zZ87E7/fzwAMPVDmezZs306NHD7p06QJAfHw8I0eOJDq6bvfp79q1Kx07diz1e8zPz2f58uU899xz7Nmzhy1bttRqTEqoQnBq9XvB9x90MBq0bl3hvgeP5zPr7Z30mLGGR5b8jXV7coLbro+uz8ieN9VorCIiIiJXonbt2uHz+Rg5ciSrV6/m5MmT5fY5ePAgGRkZpKWlkZaWxuLFiys83+HDh1m9ejVt2rSpybAvKj8/nw0bNnDXXXdV6bhly5bx9ttvs2/fPj788EMWLlwY3Jadnc2JEyc4cOAA8+bNq3JM3bp145133mHq1Km8//77nD17tsrn8MIHH3xAVlZWqd/jihUriIqKYvjw4QwYMOCi5aB37974/X7uvPNO9u/fX60xhVfr2a4BrqiI05v2Bpf/2NmYFF2+avPTUwVMfvXvZOw6Vm6bGXz3m9czZXAHGkXoVyAiIiJ1x472HWrtWh127qhwW0xMDOvXr2fmzJk8/PDDZGdnM2jQIObPn09iSd/1xYsXk5SUxLe+9S1iY2N57LHH2Lp1a7DWBWDo0KGYGXl5efTt25dp06bV6D0lJJQeaGzDhg3ExMRQXFzMDTfcUKVzjRs3jqZNmwKQkpJSqlYmLCyMadOmlWsCebmSk5NZuXIlc+fOJT09nXPnzjF69GiefvppfCVdWWbPns2cOXOCxwwZMqTCJnVf16WulZCQwNmzZyksLGTChAkMHTo0uG3RokWkpqbi8/m47777GDduHM888wz16gVaga1du5bu3buTn5/PlClTuOOOO8jMzCQ8vHqew1VDVUVn/vwnzhcEhkvPiYbtLYwWMeUTqvQ/7imXTCXG1Gdcv7asf7wvvxv1T7RuElUrMYuIiIhciTp06MDChQs5fPgwWVlZfPrpp4wfPz64ffHixaSlpQHQtGlT+vTpU+6B/7XXXiM3N5eMjAx27txJTk4OF3PhYAsdO3YMOeacnBxOnToVfHXo0IG4uDjCwsIu2serMn6/P/i+UaNGpQbUaNKkCQ0aVDx/aceOHYP3s27duovuM3DgQN58801OnDjB66+/zsKFC0v1L5o4cWKpewklmVqyZEmpwTsqcqlr5eTkkJeXx+zZs8nIyKCoqAiAQ4cO8d577wXLwZAhQygsLGTVqlXBY3v37k1ERASxsbGkp6ezb98+duyoOJmvKiVUVZS3akXw/YYOxvVRfq5rUH5C379/8nnw/T+3SWDeiFt5//G+/PT2djSLbVgrsYqIiIhcLdq3b8+oUaPIygr0Rf/LX/7Cnj17mDFjBn6/H7/fz8aNG1m6dOlFB2no06cPo0aNYuLEiRc9/8UGW6gujRo1okePHqxYseLSO18mM6t0+/bt24P3k5ycXOm+YWFh9OvXj759+wZ/vtUlLS3tooN3hMLn8zFhwgQaNGjA3LlzAXjppZcoLi4mJSUFv99P69atKSwsrLT5p5nhnKtwe1WpvVkVuPPnyd3wt+DyX9sYvZv3Llegj3xewN5jX32DMOe+LsQ2ikBERESkrqusGV5t2rlzJ6tWrSI1NZXmzZtz6NAhli5dSvfu3YFAM6/bb7+91INzQUEBSUlJrF69mpSUlHLnHD9+PC1btiQzM5POnTvX1q0AMGvWLPr3789NN93Egw8+yHXXXce2bduYMWMGL7/8cq3GAvD6669TUFDAgAEDiI2NZfPmzaxdu5Zf/epXtR5LVT3xxBOMHj2aMWPGsHjxYqZOncqYMWOC2zdt2sTw4cM5fvw42dnZFBUVccstt1BQUMCUKVNo1qwZHTpUX9NW1VBVwYnf/JpzJwPDnZ+pD7uaG8nNSmf87+38jEHp6zh7LjBp7w2NGyiZEhEREami6OhoNm7cSLdu3YiMjKR79+7cfPPNPPPMMxQWFrJs2TLGjh0brJ3y+/20atWKESNGVNg0rUmTJvzgBz9g+vTpNRZ3bGxsqfmUnn32WQB69uzJmjVrWLNmDa1btyY+Pp7Ro0dfdP6s2hAXF8f8+fNp27YtMTEx3H///UyaNCnYdA4CSeCF91K2f1h1qsq1Bg8eTFxcHE8//TT79+/n0UcfLVUOvv/979OmTRuWLl3K0aNHSU1NJSYmhtatW7N//37eeuutYP+q6mDVWd1VXbp27epqe7jDSzm7dy/7hqTgzgV+Xq8kh/F+/xv4/Z2rqeerR9H5Yma/u4vfrP04eEyYQfo9XUjp1NSrsEVEREQqtWPHjmr9tl5Cs2bNGh566CE+/vjjS+8sNaai/w9m9lfnXNeLHaMmf5fBOUf2pEeDydS+RHith/GzWx6mnq8eO46cZsKybfzjyOngMf6YBvz63i7c1ireq7BFRERE5AqRlZVFq1atvA5DQqCE6jKcWfdn8v9xAIBzYTB3sI/rohNJaT2EOWv2kP6nPRSd/6qm7zvfbMKzd3cmPlJN/URERESkcj/5yU944403amxIcqlZSqguwTnHsae+mq9gTSfjQKKR5h9F6m+2lBrNr354GJMGfJMHe7UiLKzykVdERERERADS09NJT0/3OgwJUY0nVGYWDywA+gM5wM+cc/9b09etLoVZWRR+HJgz4ItwWNErjPCiG3lhVSzwVTLV+cZYZg/vRJvrNbeUiIiIiMi1ojZqqJ4DvgASgc7AKjPb5pyr3gH+a0jhlj8H329ua5yMCiP/wGC+HCAxwhfGT/u34+Hk1vhUKyUiIiJXIOfcJec1ErnahTpYX40mVGYWCQwDbnbO5QHrzewNYATwRE1euzoUfHGejHdfpn3J8sEmxtnPvsf5gpZE1w/nnttuZGTPljSPa+RpnCIiIiKh8vl8FBUVERGhvt9ybSsqKiI8vOrpUU3XULUDzjvndl+wbhvQp+yOZjYaGA3QokWLGg7r8uSdO0HRyePB5dzIOPz2PR5MacVdXW8kqr66oImIiMiVLTY2lqNHj9KsWTPCwjRFqVybiouLOXr0KI0bN67ysTWdEURxYUejgM+B6LI7OufmAfMgMA9VDcd1WZo0iKdZdBw5MSdIOA2D+/2C2YO+q6Z9IiIictVISEjg8OHD7Nq1y+tQRDwVGRkZ0uTFNZ1Q5QExZdbFALk1fN3qEebjtpfXkbE3i8jio/Rr1w9TMiUiIiJXkbCwsDrTOkjkSlTTCdVuINzM2jrn9pSs6wRcEQNSANTzhXF7uySvwxARERERkTqoRhvKOufOACuBJ80s0sx6AUOAl2ryuiIiIiIiIrWhNnoePgI0BD4DlgI/ulKGTBcREREREalMjQ9T55w7AQyt6euIiIiIiIjUNgt1AquaZGbHgANex1FGApDjdRBSZ6g8SFkqE1KWyoRcSOVBylKZuLLc5JxrcrENdTKhqovMbItzrqvXcUjdoPIgZalMSFkqE3IhlQcpS2Xi6qHZ20REREREREKkhEpERERERCRESqgu3zyvA5A6ReVBylKZkLJUJuRCKg9SlsrEVUJ9qEREREREREKkGioREREREZEQKaESEREREREJkRKqSzCzeDN71czOmNkBM7vP65jEO2ZW38wWlJSFXDPbamYDvY5LvGdmbc2s0Mz+x+tYxHtmdo+Z7Sj57NhrZslexyTeMLOWZvZ7MztpZtlmNsfMwr2OS2qPmf3YzLaY2VkzW1hmWz8z22lm+Wb2npnd5FGY8jUoobq054AvgEQgDXjezDp6G5J4KBw4BPQBGgM/B5aZWUsvg5I64Tlgs9dBiPfM7HZgJvAAEA30Bj72NCjx0lzgM+AGoDOBz49HvAxIat2nwH8Av7twpZklACsJPEvEA1uAV2o9OvnalFBVwswigWHAz51zec659cAbwAhvIxOvOOfOOOd+6Zzb75wrds69BewDbvU6NvGOmd0DnAL+5HEoUjdMA550zn1Q8nfiE+fcJ14HJZ5pBSxzzhU657KBtwF9MXsNcc6tdM69Bhwvs+lOYLtz7v+cc4XAL4FOZta+lkOUr0kJVeXaAeedc7svWLcN/SGUEmaWSKCcbPc6FvGGmcUATwITvI5FvGdmPqAr0MTMPjKzwyVNvBp6HZt4Jh24x8wamVkzYCCBpEqkI4HnSiDwpS2wFz1nXnGUUFUuCvi8zLrPCTThkGucmdUDlgCLnHM7vY5HPDMdWOCcO+R1IFInJAL1gLuAZAJNvLoAUzyMSby1lsAD8mngMIFmXa95GZDUGXrOvEoooapcHhBTZl0MkOtBLFKHmFkY8BKB/nU/9jgc8YiZdQb+Bfgvj0ORuqOg5N//ds4dcc7lAM8CgzyMSTxS8lnxDoF+MpFAAhBHoI+diJ4zrxJKqCq3Gwg3s7YXrOuEmndd08zMgAUEvoke5pwr8jgk8c53gJbAQTPLBiYCw8zsb14GJd5xzp0kUAvhvI5F6oR44EZgjnPurHPuOPAiSrAlYDuB50og2Hf/G+g584qjhKoSJW1ZVwJPmlmkmfUChhComZBr1/NAByDFOVdwqZ3lqjaPwIdf55LXC8AqYIB3IUkd8CIw1syuN7M4YDzwlrchiRdKaij3AT8ys3AziwVGckG/Gbn6lfzuGwA+wGdmDUqGzn8VuNnMhpVs/wXwoboRXHmUUF3aI0BDAkOeLgV+5JzTNwfXqJL5IX5I4OE528zySl5p3kYmXnDO5Tvnsr98EWi+UeicO+Z1bOKp6QSG0N8N7AC2Av/paUTipTuB7wHHgI+Ac8C/eRqR1LYpBJoDPwHcX/J+SslnxTACfx9OAt2Ae7wKUkJnzqlVgoiIiIiISChUQyUiIiIiIhIiJVQiIiIiIiIhUkIlIiIiIiISIiVUIiIiIiIiIVJCJSIiIiIiEiIlVCIiIiIiIiFSQiUiIiIiIhIiJVQiIiIiIiIhUkIlIiIiIiISov8Hsi9BZHYTB3AAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "fig, axs = plt.subplots(nrows=2, ncols=1, figsize=[12, 5], sharey=True)\n", "for i, (ts_S1, ts_A2) in enumerate(zip(ts_S1_list, ts_A2_list)):\n", " axs[0].plot(np.sort(np.abs(ts_S1)), p_S1, lw=3, label=list(tsDict.keys())[i])\n", " axs[1].plot(np.sort(np.abs(ts_A2)), p_A2, lw=3, label=list(tsDict.keys())[i])\n", "for ax in axs: ax.legend()\n", "\n", "fig.tight_layout()\n", "# output\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "id": "d3d6c893", "metadata": {}, "outputs": [], "source": [] } ], "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.8.12" } }, "nbformat": 4, "nbformat_minor": 5 }