{ "cells": [ { "cell_type": "markdown", "id": "290cde5c", "metadata": {}, "source": [ "## Single path absolute tropospheric delay time series" ] }, { "cell_type": "code", "execution_count": 12, "id": "2a6bd0b6", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Go to directory: /Users/yunjunz/Papers/2021_Geolocation/figs_src/tropo\n" ] } ], "source": [ "%matplotlib inline\n", "import os\n", "import datetime as dt\n", "import numpy as np\n", "from matplotlib import pyplot as plt, colorbar, ticker, colors\n", "from mintpy.utils import ptime, readfile, utils as ut, plot as pp\n", "plt.rcParams.update({'font.size': 12})\n", "\n", "data_dir = os.path.expanduser('~/data/geolocation/ChileSenDT156/mintpy_offset/inputs')\n", "work_dir = os.path.expanduser('~/Papers/2022_Geolocation/figs_src/tropo')\n", "os.chdir(work_dir)\n", "print('Go to directory:', work_dir)\n", "\n", "# file info\n", "lookup_file = os.path.join(data_dir, 'geometryRadar.h5')\n", "tropo_files = [os.path.join(data_dir, 'ERA5.h5'),\n", " os.path.join(data_dir, 'ERA5wet.h5'),\n", " os.path.join(data_dir, 'ERA5dry.h5')]\n", "lat, lon = -21.29, -69.22" ] }, { "cell_type": "code", "execution_count": 13, "id": "e814293e", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "input lat / lon: -21.29 / -69.22\n", "corresponding y / x: 109 / 220\n", "input lat / lon: -21.29 / -69.22\n", "corresponding y / x: 109 / 220\n", "input lat / lon: -21.29 / -69.22\n", "corresponding y / x: 109 / 220\n", "Done.\n" ] } ], "source": [ "dates, tsT = ut.read_timeseries_lalo(lat, lon, tropo_files[0], lookup_file=lookup_file, zero_first=False)[:2]\n", "dates, tsW = ut.read_timeseries_lalo(lat, lon, tropo_files[1], lookup_file=lookup_file, zero_first=False)[:2]\n", "dates, tsD = ut.read_timeseries_lalo(lat, lon, tropo_files[2], lookup_file=lookup_file, zero_first=False)[:2]\n", "print('Done.')" ] }, { "cell_type": "code", "execution_count": 16, "id": "668fbaa8", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "save figure to file: /Users/yunjunz/Papers/2021_Geolocation/figs_src/tropo/tropo_TS.pdf\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAewAAADMCAYAAACx8ZDiAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAACRSUlEQVR4nO2dd3gU1d6A37PpvRESQgm9NykiijSx4AVBwS5ybSg2xHZV9Bp7+ezliih2saDYxYIigqAISO+9BgjpPZuc74+zsztbs0l2k00y7/Psk91pe2YyO7/z60JKiYGBgYGBgUFgY2roARgYGBgYGBhUjyGwDQwMDAwMGgGGwDYwMDAwMGgEGALbwMDAwMCgEWAIbAMDAwMDg0aAIbANDAwMDAwaAcENPYD6wGQyyYiIiIYehoGBgYFBgFNcXCyllAGpzDYLgR0REUFRUVFDD8PAwMDAIMARQpQ09BjcEZCzCAMDAwMDAwN7AkLDFkKEAf8DxgCJwE7gPinlQi/2/RUYBYRIKc1+HaiBgYGBQfMjIy4RmAucBWQB95KRN8/NtjOB/wARwOfAdDLyynwxjEDRsIOBA8AIIA54APhUCNHe005CiMsJkEmHgYGBgUGT5VWgHEgBLgdeIyOul9NWGXFnA/cAZwDtgY7AQ74aREAIbCllkZQyQ0q5V0pZJaX8FtgDDHS3jxAiDngQuLu+xmlgYGBg0MzIiIsCJgEPkJFXSEbeMuBrYIqLracCc8nI20RGXg7wCPBvXw0lIAS2I0KIFKArsMnDZo8DrwGZ1R3PbDYzaNAg62vOnDk+GqmBgYFBM+GfD+CL6VCU1dAj8T2r3obXR8DrI2gRKRyttl2BSjLytuuWrQOcNWy1bJ3DdilkxCX5YpgBZ04WQoQAHwLvSim3utlmEHAaMANoU90xg4ODWbVqlU/HaWBgYNBsyD8MX92k3se1gdGzGnY8vmbQVeoFZN0gHGOhooE8h2V5QIyLIzluq72PAU7UdZj1omELIX4TQkg3r2W67UzA+yhfwc1ujmVCBajNMILMDAwM6pXSfNj2A5Q3szTR9Z/Y3m/+EnL3w6+Pwp6lDTakeqQQiHVYFgsUeLGt9t7VtjWmXgS2lHKklFK4eQ0DEEIIVBReCjBJSlnh5nCxwCDgEyFEJvC3ZflBIcTp/j4XA4Mmi7kclr0Auxa7Xl94HFa/q/42V767HT66GH66X302l8EP98LKNxp2XLVlxyJY/ITnCYiUsFYXEB0WCz89AL//H7w7Dj6/Tt07TZftQDAZcV10y/rh2mW7ybJOv91RMvLqrF1DYJnEXwN6AGOklJ4S1/OANN3ntsBKVIBaM36SGBh4yZKn4cBKOH82RLWwLf/zf7DoQQiJgrt3QYhDdcD5U2HfH7D1W7h8fv2OOVDYYDnvVW/BuOfhn/fVdQNIPxVSXLk1AxRzOXw4Sb2vKIazHnG9Xf4hyNK5b3P2QLHOj73hUwgKgYn/c95XShACCo7Cn69Ct3Oh3Sm+O4f6ICOviIy4BcDDZMRdC/QHJgCnutj6PeAdMuI+BI4A9wPv+GooARF0JoRIB65HXYhMIUSh5XW5ZX07y+d2UpGpvbAJ6aNSyiY9zfOKSjPkH4Ej69Tr+HZlvio8DmUFar2e0nw4vBY2LoClz8Lqd6A4uyFGblAfFGXB4sdg589K6GhICWs/VO8riiDvoPO++/5Qf3f8ZL+8vEjdR4GOuUydZ22pqrT/XFkB/3xo+/zaqfDRZVBRWvvvqE8KdfG6y1+CHT/Da8Ng05cO2x21/1x8AnL22i9b9xFUValrol2nxY/DM11h9xJYcB388SK8dbavz6K+uBGVV30M+AiVW72JjLh2ZMQVkhHXDoCMvB+Ap4HFwD7L60FfDSIgNGwp5T5AeFi/H+XMd7Vur6d9mwUbPlM/uPzDlgjOah5KIkhpTyIIyhxjKYDv7oSuZ0O/S6HLWRAc6pdhGzQAWTts7/f/aXt/aI29FpW7H07sguNblVaU2NH+OFVVypd54C9Y857SsG5ZY6+xBxLHt8ObZ0CLLnD1TxBUw0ffj7PsJzighFbBEftl276DnYugx7g6DbdeKDxm//nDyerv/KnQy/JcKMl1Fs56UnrD0Y0gq6AkG94coyZw1/0CS55S27w/Ua1vzGTkZQMTXSx3lk0Zec8Bz/ljGAEhsA3qgJSw8G416/V6n0ooL3S/vqpCmT23fguRLeDSj6Ht4LqPtTEiJZzYqR5aHUY0/slL9i7b+4N/K9P4ilch74D9drt+VaZeWQWLMmCcw/Nn0wL4/Brb5wrg4Crodo6/Rl43Ni2Asnw4tFpZntq4KPFwdBN8fzf0vQgGTrUtz9kHK15x3n77j84CG+zNxYFGVZWypEXEQ0wr99sVnYDSXHjtNDB78FC26q8ENigLTM4e9f7Di2zbOArrqiowBYRxt9FhCOzGTuFRm7BOHwbpQyE6Rb2ECcylUFHi+m9lOcSmKe0psSMktIcj62Hdx0p7Ki9UD59vZ8INS5UvqjlQnA17liihtWuxTZj1vwImvtqwY6srJ3QCuywfPrrE9WTPTkBJFWSkZ88S532KAjiE5PBa2/uDK10L7AXTlPDZtwwGXGm73zd/5fqYq992vTzQArAqSmH3YmhzMhxaBYsfVcsTOrjf55nO3mnFrfrCWsv7Ap2J/fgW9/uUF0K4Y9C1gTcYAruxc0z3wxh+B3QaXbfjdThdvc79P/jlIfhrNhzdoH7wro6942dlNh0yveZmxkChskJph7t+Va/Da1w/rNZ+CENvgpSe9T/GmmAug+Aw1+v0GjZ4b5lxtMjsW6H+xqdD7j7LsQJYszyy1vb+wEo4ZbrzNpqmCEr4xFo00M1fuj7miZ3qb1gsXP2D8mGDZ+tVQ7D0Wfj9aWh/OnQaZVuuacOucHX/h0RCdEt7E3lUsu19/iHvxmMI7Fpj2CUaO8d1tWWSe/juuKGRMOI/6kcKsPxl520O/wMfXqhSXJa/5Lvvri8qSuC3p+DpTvD2OeqhdmiV7WEVHAGdz4RRs5S1AqkCtgKZ7++Cx9Pg29tdp+pk7/a8f2of+8/t3WRKnrD4wjuOsN0j7ipgFWXBB5NV+lBDUJBpb7o+6FBEyVUgmubPL8lRZnRPtB4ILXtiDaWpKHa9nZSw6QsVO1BVqdLB3hnn/zS5359Wf/cutY9hqClRLeCsR22fO46y/e9BxdB4Q1mATWgaEc1CYFdWVjJt2jS++eabhh6K79E07PA4iEn17bEjE+EkS7ncXb9C5gbbuqoqFZymBbj99XrgmQLdIaWKhH3lZPjtcfvAu9Q+cNoMuPIr+M9euOIzGHE39LtMrd/6bfUP8IaivFhF+VeZYdVcmH26vQVGSjjhQWBPeBVadLVfNvZpiG/nfp+WPVWcAyjBXJqnzM96Ibj4cRWVvuRJ19Hn/kZvDgfI228z3/5wHzzaEj690n4bTWB7I0zbDFbmc014lbsR2Ht+h/n/VpHSb4xSMQJ7l8L6j709k7qz/Yfa7xuVDD3Gw7+eVf/3025VE3uNfBf+fFcsvEu53QKNbQvh61uJCyOooYfijmYhsIOCgpgzZw7jx4/335fk7FOaZl1msLVB07CTe/jHxzz0Rot2CSzX+TXXfqC0UY3CTNj4me+/v64cWQdbv1d+yA2fqTScd8erSNi8/WqblD5w/hy4cwfcsAzOfBg6joSQcNtxRtwNphD1/tdHnb4mINi/QsUlaGTvUvEHGoVHVcoWOAcc3bkTTrpCxTRopA1Q5v8xGe6/s2UPW2R4wRF4518wZ4QKhFz1lsq/3bfctv3xbSqI7ZMrVJphfaA3h2scWKm03FVvqWvm6KfWfsclOdUfv40lIFMTXnuXqkA+x/PTT/SO6MpNe4rCriuOY/DmfNyhmb8HXws3rlAuMjsN20uT+O7f4IvrIfdAtZvWK93GwnkvkVdGZfUbNwyN1OkYgMyfqkzEh/6Bq76rn++UEo5ZBHbL7v75joT20HOCMuVt/AzOeED9SBdlqPWxrZVGV3hUCfR+lwZOcNqGz+wjmR2JSFTnM2AqmKqZVCekw8B/w99vKGvDkqchc706xpgMZY1oaHb/pv4Kk0rH2/6Duie1qFx9wFnfi1ReLKhz0ISu/gHcc4L62+sCWP+pa+2sZS/bvvpAtJWWBju7Fivrj8baebaJXat+MPyuWp1qjdA07NjWSrOWlSrwLLm7+whoq0lcV5OgVX9lZZr0Bnx2tW15m0Hqb2iUCrzLXK9e6z6C636zxXaUukihBLXPpi/g18fgjP9Cz/NqeaI6CjJVAKkv70tXKXu1MYlrHN8K8W3rNqZmRrPQsP3OoTXqwQgqL7XCU6E2H1JwxGbO9aX/2pFTb1F/q8wqCO3XR23BSmc/BidPU++PbVLBaYFAaT78eJ/rdaZgGHID3LoGBl1dvbDWGH6n8muD8mVv+QbWvAv/G6pybxsaTWCnDYDOY9R7cynkW8zQ+oCzPheqPHxQWrI2yep7sXoIx7eDwZbJjhBwyTy49lf774tsAdHJ9oFHrsZUpSvWo7fCHPMQSexLtN9mu1Mgtbd6f3CVEqru0AS2vojQhe/AvQeh9yQ4x5Jj3G6oTSiGRNkfI3MDzLvQ5gbQ5z2fcpOy7IAyJc//t4oL+NRVx0Yv2LtMpaWBmsi/f7767u/vdL19at/qj5ncHdJPs32OdCWwddXwKsu8Hy+4Tokz8IghsH3Bqrm291UV9efj1D/w/KVhgwqqSR+m3v8911ZAosMI6DlRCT1rcJqLfNWG4Pf/s1VoOusxmL4cbvobbl0Ld++BsU9BRELNjhmTCqfcYPusCe/CTPhgkgr4cue/9DdFJ2wCqONIVSBEQ4tm1jRsU7Ca4HU+w7K9LnI4qRPcsRVuXgVhumZEpiBI649djSItWt5R84pqaXtfUezebx2d4sWJ1ZGCTFtFr1b9bebrw//YBLkr8g8p/7XehByRYDN7nzwNrlsMl+smIHp/rsauX9WEbtevtvsxbQCc87jt+tVVcB34W7ki3hyj7oPj2+DYZrXO3Tl2HFH9cU++zl5gCxfiIjTKeZm36C0+Bl5hmMTrSkkubPjcftn+FdB+mP+/218R4q447VaVn6pFwJpC4NxnlPYVmQj9L7eYi3+Bo5sbNvUpawf8+Zp63+ZkOOVG3xVqGHU/JHVRQqrjKNXF6Id7VKrKyjkqza3tEPVwj0hQ16bLWcqk7k/2/GZ733EkJOpybE/sUv5GTcOOT1dm2gvfUVpZa4ecZL0JW48pSK0rzVWfW1r+x46a15Qv4MCf8N0dFneJm5b15hpqZLVBH3CW1l9NEv5+U1ke1rzned+9S20mcRFkf11MJmg9wH77EAeBHRSq/ONl+aogS7AlJkKbqGhBogVuro+3LHlS/a0oVmP2JrBv8LUq8LLgCHQ9RwVTOhLV0t6M7+r/5VhvHiA02rvUtuoyFgyccCuwhRC/e3mMUinlWT4aT+Nj3cc2P5j2A9WXfPQnmoYdkaDyI/1J5zOhRTfI2qY+D70RknURxUNvVA9CpAq4aagCI1LCwv8oSwdC5ZP7sqpSUDCcdLnt84ApanL2xQ1KSOXscc5vjWsLM9b7t7qTZg4PjoC2J6sJVUikeohrAVRahHhSJ/U3NEptWxMiEnQC2zJJdDSJt+zhndZYlwAoPRUlyqRcmq+akoTpKkXqNcxW/SCps7o2VRVKkAK0OxX2L8eJPbpHYERC9bEZoQ7Vky/+QFX/+uNFZe7Wgha132qMJcCvyl1jQi/RBxqWZLt2z3Qdq3LPFz+mrAwJ7eGmlUqwbv7KtcCObqkmfr88rD67KrfqOEkBVa70gBfPwGwPeeAGLvH0BBmManfp6fUWMMTPYwwM1rzv/ICR0mYeTugAvS21eLUIVH9jjRDv7v9AL5MJRll8womdYPjd9usTO9p+0Bs+VdHBlWY12z/wd/2l82xbqLR8UOUl0/r7/zsTO8BV36sc1bSTlP83VGdOzjvgX22iOBt2Ws45/VRVNMVkUv8nUCbxqirbGLTltUHvRmhp6UylN4kHhSlN3FMqmIavBPZvT6iAuP3Lbf97DS1QLrm70pBj01Sd6+7jsJr3+11sv4+We753qW2M3rhPHE3i4XHQfrjtsyaYNQ1bK8xSV/RumKyd9lH5GoOvVWbwa35ScSfaeKNbup/sR7dUaY4XfwgXvKnuLUdMQep/rkeLE6iO7N21a8Yipe/unUaGJ5P4cinlu9UdQAhxmQ/HE5gse0G1HVw5R+XnakEm+5bbNM5BV6lo23Xz1Mz96CZVts9fSKl8VaAeRvVBr4mQamnwEBbtvH7oLSoQq7IcXhkM5QX2RUiu+VFpObWhNA9Wva0etC06u96mohR+vFe9D4+D0Q+43s4fmIJUcJ4WoAdwcDW8aakOd2St+3HXhU1fqHx4rcqY5pcG9X1HNyjtruCIzRKU5COBndxN/dX7uhPaq79xbZz3TexkH/imaepga8PoDfmHVYxCr/NV0JcW7Q72aUwndtlSurSId1D34CUfqvW5+1UsxjczbOs7jlTC+sROW2CeN9HWjtpmWKxynzhi1bB9ILC1Wvcameudg79Co+3vC0ei3AhsbXl1jUxCI6FE952xraHPRaqrm/5/7EhFkQrEi6lhLMNXN6uqgxfMUdkOzQi3GraU0sN/2G67pm0Ol9L2o89cb1+ZSAs2CwpVPlz9DHT/Cv+OK/+QzaTX0s/+az1Jndz7ONsNsQX1lOXZlzc0l8AX02tfXOX7u9Sk6eNL3c/KV7xiy2kddX/Dd45K6aUCvMA+79ZXZG5QpmBNWHc5Gwbp0tiSLBOE3AP2gZCOnbdqgiZkEjrYyku27GETVmMt0dOhUc6+7cvnK+GocfgfeP8C+OUReKq9xaXiBW+OUZatb2+3tfzU0FfR2qiLLek9yfk4SZ1UqU6TSQl+7bz0bgJtQh7hhcB2DMAKj4OoJBUzoMfqw/aBwC4+YS8U9eVVe4yHbv9SWrWnyZA7DdvVpNwVjhOViHiV+nb3btv9746aWt4qzaoGBFK17GxmGFHi1SGEMgf1uVB9PrYJ3h2nNOjNX6tlPSco4ZDY0ebP87fAPqYPOKsnDdsbzntF5e32u0zl2I57HgZepdYd22Qrk1gTygpt1zpru6r17Yi5XFWOAmWqHXS18zb1TUi4LRjQHwJbnyUw7gW47BP7Yi9W7U7aSscGhdryhmvDqbcozfbcZ2zLwuPg2kXK+qSvVa3XvMPi1O9j6te2NEBQJuylzyih890d1X//0U22Ah0ndjhXItMmsVKqPHxQ6VOaNcAdk99WFpnLP3Od8uSVSdyFwAbn4DSrhu2DyoT6lqhgbyruMQEunacmjp5wGa1fAxebY+CZdq1MQfYpfa4oL1RWue/vUm606tBbE5ohXglsIUQ/IcSvQohsIUS55VUhhGgUtSjrXJo0KBjOf91WnvL4VpgzyuaT0rQaIWwz9X0rauef8RZ9N5z61LCro2V3uPBtOP81GH2/Epxjn7bVqF76XM3T3rYttC9wsXGB8zY7F9lyw0+bETiNSDQXQOZ6398P+kIVvS9w1qL0JviDf6u/nUa7t5B4Q8vuKrq8yxj75Sm9lClZjz5+YNR9tvGFx7s/vpQqR9pdEY6/Xrf/XO5QyUuLaj66yaYd977A/fdpxLZSefYtOisN0VErro1JXBPgenM82AR2UIjna+EN+kwRp/G4iOB2haMmPeYhlQbpLY7nrT+n/leovz3OswXZ6ck/BJ9do9yNH0yqvvqdvjyyLywUeppQadKPgD+A4UAPy6u75W/A45PSpKYgVWt5gKXmsOYnSu6hCjJoaAK7MNO/JQc1DTsi0XPhikAgOBQmvqaiZGUlfHmj8jd7i2PJ040LVBCVnnUfqb8hUdX73OoTLY6hJMe553Rd0SKxQ6OVv9SRJBc+854TfTsGTwyZrrTbMx+2z1/3pK3uXARvngGvj3DtPnEsdqIvbAI287CdOdwLge2IY6xFbTRsbYLSc6KtlaUpBKJ1mrWr/1tNAlb3LPUwHhcR3N4w7LaapWU6mcR11+pfz8K1v8CF79pHs2us/9T2LD26AeZfpczeGuZy+8/6/78WL+ErGkFpUm8FdirwXynlRinlLv3Ln4MLOEwmGPeivZ9w0NX2mo1eePszvUvTsPVVqgKZ1D6q+xcoreA3Lzs36SOgI5PU34LD9mkjJTm2aOCeE+pWzMHX6B/8vjaLa1poTCvX90B4nH1AkSlEPZTqi3ZDYPoyZfHQ40n4fXSJ+lt0zHVt6tz99p/11cNApXaBKlQCKse8Ng/22ghsveASOiVNCLjmZ1Wrfdzz9oLUVZtJb/PTqyptlQVdReW7Srlyx/C7lL95wv+830fDcWKgDzgMCVcuGJNJtbF1xLEy4s6fVTnh8mJlJXm6A8wZaRPadg2IAlau+g1vBfa7QNOPBvcGk0nNGsf+n0ptcvSVpva1lSh0ldvpCxoiQtwXDJupqk2B8qke+Lv6fbZ8Y3M9nPuM7UGo16A2fWGbvTum6DQ0Kb2x+gN9LbA1DdtTepBey+40Wpl7GxpPY9D7PM0OVpjyIuf+3Y5FWTSTuFaMpLYFhbTgSQ1vAhj1gisoxH5ddLLFQudQelTv59fwtsTn4X9sPmtHszt4bxIH5b6696B9jQFv0U8MQqJsE2tHPJ1X78k2t9nmL1XzmA8vUj5uLdMB7AW2uzrwTRhvBfaTwCNCiE0WX7b15c/BBSxCwJBpMHqWs680KBjaWn7s/tKw8w7YKgkFkv+6OoKC4fzZKvBJVqm65NWhCeboFPVQ0oKaNn1pm3Wv+0T9jW3tvn9zQxEWbSsT6nMN2yKwXfkGNfR+7F4Tffv9tcXbkrCO/kxX3Z0cA5XK8pW7pMgSjBZdS3dRh+HKlN+qH7Q9xT663R36winVRUdruBLY3mrYO362ve/hSmDX0NJUEwFvt59OYCeku7f4uTKJa3QaDVd+bZvQ//O+rQY+QN4hpahoGRFQM7daE8Fbgf0ZsAd4DfjQ4WXgSDtLelfWdtUj2NcEaoS4N7TsYQtOqk6AFRxV+bCgIpNNQbb0nOIs2Pu7Kr6gmcf7XuR9I4/6RDOv+lJgV1V6p2Fr3x0cUb/mcE94LbDz7T87msPBtYZdkqNiJcB9jnF1CKFM+df/ruoHeGOZ0Asub+9DVz5sbwS2lLbYjtQ+9lUHreOppQCuKfrvcQzW03O+JWAwJs213zsyUf3OXZF3wPm6OFpgmgHehtL2B5KklI0iKrzBcfRj+zoIKlAjxL0ltY8qqpC9S/mq3AXHbP7SlsutVZHr/i9bCdiNn0Oszl/W9xK/DrvWtOoHG+ar5g8Fmb5J5yk6bhNKnjTs/lcoTaRVv5o3O/EX3kZGlzoI7DxXAtuFD7tIt8zfJXv16O9jU4j77fTUVsM+/I8txanPRc5lUR3H408cNWx39J6srGBJnWH2MFtfArBlLrhTQPIOOgvoZiiwvdWwlwIN2M2hkdFmkM0k5o98bE3DjmzR8MVBakOKpXShrPLcYlHLo41vZ8sdDo9TzTQANn9jiw5v1d+/Hcvqgj8Cz/RpT7EeBHZIOJx6M3QIIFdBrU3iLgS2/qEPSsPWC/H6zKDQm6AdfdjucCmwS1Wgpacsk01fWN4I6DNZafSOWmtNgs7qgv5/4EnDNpmg/WmqspljKpkmsN39hvMPOU9kDJO4W/YAPwkhXhdCPKx/+XNwjZbQKNtD2h8CWx8h3hjRgktABZS4ImcfHFyp3veeZO8X09J0yvIgd5963+9S34/TV+gLcfhKYOuba/iqJnV9ERyqCpTorSOu8EZgO1JVAR/rAqfqU8MODrW999aH7SpKfO2H8MEF8OaZ9ilNerQGMy262iZsei1bBClLVH2gd/t5az1yV2Qmzk0NepcathF05o5I4DsgFGire1Xzi2vGaPnYR9ap6FZfUVXVOCPE9SR2tPWSztzoeptNuuIojmUlu57jnELjqvRkoBARb0stqq3AdkyJ0WvYnkzigUqXM+GcalL73AWduarPrUdfTKW2PuzaoDf19/Cy5oMrH7YWjFl0zNbe0xHt2uh963qtNSSy/tI99dXVvJ0gOZrwtfNw7GinldHNO+isYVeZ3U9omiheCWwp5VVuXgFQ/zFA0eqKV5lV1TNfkbffZoIKVBNwdZiCbIUZjroR2Fp0eItuNhO6RmiUfQBVlzNrHw1cX9Qk8GzfCjikK7/60wPwRFvY+p1tmaZhi6D61SJ9iSvtUo+jwNZ80y1cCGxXEdEiqH799pGJKnXr5GkwapZ3+7gyiev5dib8OMu5Sp52bfT767XW+vJfA4yyNNwJi4O0AZ631bCzBpjsP7fUeV+7nK3+5h9ydn9As9Oy3QpsIYRXIYbebteQ1Lk0aW1oP8xmFnPVn7a2rHzD9r62OaaBgCaEMzc6Vy0rzbPlW/ac4FpT0ILQAPoFaLCZHk1g5x2AohPutzu6Cd4+B+aeadMol7+kHkwf60ohaCld0SmBGRnvDa60Sz2OUeJaVbO4ts7buvLjR7Xwbw9yV5x0herB7m3jjOquwdZvVVMbx37VLgW27n19RYiDSqW89heY/of3EwX95CI8zv43PvE1FZNyzpM2y1RluetGId6mwHlDIy9N6kUldgBclCMKLHxSmrSmhMdBW0urcF8J7BX/Uz9eUL4rfVehxobmxy4vsPmhNfQ+MXedpbqNhZH3qepprnJQA41UXeBZpgctW8vdrzK7bnKioVUBa2z+az2ONc3HPATTlqh7G+w17IpSW+2B6GSbS0XDlcD2pSvKX1SnYWtobjANLYLersGK3iRej9X+hFBBofEuJlLu0GvUjlkDaf3h+iVwynSIa21b7qqnfIUPNexGUJrUU2REuBDiPS+O4WU4ZDOk8xmq/d+JHSrisy61bzcugB/vU++jU1TQjreRqIGIXeDZRkjsYPus94m5a7ogBIz8j3/G5g/0vdGPrFeFIlxxQlftN3u3+4Yh1hzsRui/1nAU2MNuU381IaQX2HpfbmSS0tD05tBY3YNdQxPwgYy3Ajs4zP6zVcPWaeh6IVifGnZt0E8uPDWj0QtzrRiOnmaW2uVJw34M2OXF60k/j7Hx0lnX0Uirh10b9i6DL64HpPpRXj7fc75jY0Df8s8x8ExfftKbPsSNgeiWtuAwT35srQQjKOHtToPwpspZoOPOHKwt15vE9feEJrD16Cd8GvUZcFZbvO2cFqQT2FWVUGGxHrjTsOvTh10b7EziHtwC+vNwFYDnSw27EeBWw5ZSPlSfA2mSpPRRD42iY0pgD76m+n0cObYFPrpM+XBMwXDx+86NCRojYTHK4pCz1znwTN+ByZu2ho2FVv1U4xKPAlvX7zd7j2uzblmBLRK6MZvEg92kHbnSsJ0EtoOPOD4dzntZmY5lFaydBxc4tOIMRLzWsHXXSn9d3PqwG5HAdnRv2G2n+z87dmYDpWHvXgLbvlfHHH5X4FsX6kCANA1uophMyiy+7iPYs0S1inP3kHLHTw+ofGNQnXTcmVIbIym9lcDWF/QHh4dzExPY2xeqCm+l+c6ahblc5Z9rZO+yaVJ68nU52I1Zw3aHVcP2JLAdNOywaOiua/xy9uONo4udtwJbjzuBrdewg8NrP6b6QD+58OTa00883AnsQ6ttqXDD7/LN+AKUeg6hbIZoZvHyQjjwV832Lc2H3b+p9yddEXidqOqKVlAkd5+tyxLYTF8iSKWKNBX0lhHHSQqoYhhauVFQfupCF367An2Vs0asYbvDpYbtwoetx1HjbgzCGrwXrPpoaLcatv6auIl9CBT0YxUexJB+O8dObWAJRiyyHSfQJyp1JCAEthAiTAgxVwixTwhRIIT4RwjhsVOBEKKjEOJby/ZZQoin62u8NaLjKKztFWsaLb7rV1tryZ4TfTmqwCBVl199dJPtvfbDjEio/7Qcf1JdiVK9OVzDsRKcuazxF03Rc97Lqrb0Fbp2qXqBraX82cU1JLjWsBsj3k4svBLYumsgHVIlAw29Vu2pKpz+/+zKh20usQns0OiGnahlxCWSEfcFGXFFZMTtIyPOfUvqjLipZMStJiMun4y4g2TEPU1GXLUWb6+ehkIINw1OfUYwcAAYAcQBDwCfCiHauxlPKPAz8CuQiqq49oGfx1g7opKgtaWYQE0Dz7YtVH9DogKvbaQv0BdE0QeeadqUu766jZXYNFX/HVwL7KwdzsuOrLf/XFZor503dg17wJVwy2r7AE2rEJI2l4AmsMPj1MPeUaN21fyiKeFWYOvcKnrh7VjbINCo0lmSPAnsoFDbercatiUbwHESV/+8CpQDKcDlwGtkxPVys20kcBvQAhgCnAHcWd0XeKu+HBBCfCWEmGwRlj5FSlkkpcyQUu6VUlZJKb9F1S8f6GaXfwOHpZTPWfYtlVKud7Ntw6M9jI5usPc/eqLSDDt+tOw/WjVxaGrEt7OZvPWapFVgNyH/NajZv6eKZ5qGrffbOZrOC47AP5autu1PD4SHlO/RCx5NOGkPa20SV51JvDFx8Qe2hjbuqNQLbF30vF5g64OtAl3D1iyH4FlgC+H5HrfTsBvwt5ARFwVMAh4gI6+QjLxlwNfAFNfb571GRt5SMvLKycg7hGpVfVp1X+OtwE4HfgH+A2QKIeYIIYZ5uW+NEUKkAF2BTW42OQXYK4RYaDGH/yaE6ONm24ZHrz3s+tW7fQ6utOUjdw2QPsa+Rghbepdewy5poho22AR21jbVWlSPJrBT+9jS2Y46/ARWvm4LQjxluv/G2ZDoBbZWIMRRYOsf8im9G3c+eo/xKlXTU2cxfb6xO5O40BXoCnSB3XGU7f3Af3ve1lMRGHOZrWRpw05euwKVZORt1y1bB7jTsB0Zjnt5Z8XbWuLHpZQvSSkHA0OBY8D7Qojdlq5dPksKFkKEoGYb70opt7rZrA1wCfASkIZqTPKVO+3fbDYzaNAg62vOnDm+Gq53tB5oKwDgrR972/eWNwK6nu2PUQUGWgGVY1tsZjK9D7upoQlsWeUsjDWTeFIXSOqk3jvWSl5jqWUUn66aoDRF9LnJmnAqchDYWs/5qJZKQ20sQWae8KRpmstt7+00bJ3ANjUigR3fFq5ZBFO+gHZDPG/rKIj1E5MKBx92bVn1Nrw+Al4fQYtIUZvsqWggz2FZHlB9GkBG3FXAIOCZ6jatzcBSLa9YYA3QGvhHCPG0lNJlERUhxG8o/7Qr/pBSDrNsZwLeR/kBbvYwhhJgmZRyoWW/Z4D7gR6oWY0dwcHBrFq1qvoz8xemIJWOtWmBJZCssvr6z5r/uu3JjbPntbdogWfmElUopEWXpuvDBoeKZ2uh7WD1viQHii0lWZM6KZPhwb/dH2fIDY23hnh12JnEHTVsy2+h1/mqeFBiJ/uOVY0ZR4uLnko3Pmy9kNJPcBtDYSXt3q8OxyIwEfG2+8HsIx/2oKvUC8i6QTi3AMuI+w0PMgy4BSUT9cQCBc6b2x13Iqr42Bgy8rI8bouXAlsI0Qu4AuVILwTeBfpKKQ9Z1j8CrMdN1TMp5UgvvkMAc1EO+3OllBUeNl+PF/b+gKLzGCWwS3NVJyZPN2vWDpt5tFsTNYdr2AWerVdBVJp/q6n5sAESOii/fVmevR9bX5K0RRfndpp6QqPhpMvdr2/sOPqwpdQJbMs9IYSyXDUlHJud6HEVdBYaY59F0WYw9LpAlbT1tltYY8BRcw6JhKACVUxKr2H7s1hMRt5Iz+vjooBgMuK6kJGnRY/2w5OZOyPuHOAN4F9k5LnI83TGWx/27yjVfrKUsqeU8ilNWANIKfcCL3h5LHe8htKQx0spq6s39wFwihBijBAiCBVtlwVsqeMY/EfnM2zvqzOLa9o1QLdz/TOeQKFlD1se5tGNzgUymhpC2LTsTF2cpD6lK6mz+6YnoILNvC1p2RhxFNjlRTYNsyneE1Z0udOO7iCzi6Azx6IrQsCFb6vGGU1psusoiIPDbNXRzKW+MYnXlYy8ImAB8DAZcVFkxJ0GTEBZjF1sHzca5fqdREbeSm+/xluB3UpKebOU0u2BpZT/9fZLHbH4wK8H+qOC2gotr8st69tZPrezfNc2lMY/G8hBXZjzpJTlLr8gEIhJVaVKwXuBndDB1rmoqRISoXy2oALP9AUymkodcUc0P/bRzTbfpOa/FiYlrJM8COzG2gfdWxwFdlOfxLnC8Tw1gV14DP6xZLDWpkpaY8TR1B0cbmuGUlESSGldNwIRqBivj4DpZOQpDTsjrh0ZcYVkxLWzbPsAKoX5e8vyQjLiFro6qB6vTOJSynJL5PbJqLwxoVv3Vg1OyN3x9+mP6WL9fpRTX79sAWpG03jofIZKXzq0WgXRRLl4+BRnwwFLi8Vu5zaNYJrqSO2toqaPbnSuaNUU0QR2VQUc36I+a00/4tuph5EnDbsx90H3hlBDYBOZZG910SwMn11tW9ZsBXaYLc3VTsNuYIGdkZcNTHSzzl6GZeSNcrldNXhbOGUiqjPXw8DrKAf767jLMTNwjTW9S8Luxa632fGTLcKzqfuvNbRI8YIj9t2qmpJZT4+rimeaDzups/obkeDewtDUNeygYJsZtCy/eUziHHHsEW0uU+V79y6zLWtKfQU84VLDtpjES/NV73hX2zVBvDWJPwpcJaU8CSiy/J0GrPbbyJoibYfY/CzuzOJaOld4nC11pamjuQoA9vxue99UH85JnW0C6cg6VZXKKrC76Lbr5Hr/pu4mAfsWm8W64Nmmek844mhZM5fB/r+w+rnPfQZG3Vfvw2oQXPmwNQ1b3yO7MRfP8RJvBXY7KeV8h2XvAlf6eDxNm+BQ6GDJDNj5i3P5QHMZ7LQUVulylucuNk0JfU1xqwYhmm5glSnIZlU4sg7yD9nyrfVC2p1ZvAm3D7SimXvXvGfpBW+hqVpdAPpcqP5GtnDOozaXwr4/1Hthgr4XNw93GTgL4uBwiGur3h/SpesaGraVYxYfNqgKY0OBTkATTQT1I1q0eNEx58YOm76w9Tlu6tHheqJTbPm1WgRsRHzTzTMGm1k8c6Pq4azRQqdhJ7rQsJtqsRRHXPlnhcnZVNyUOPcZGPcCTFusuvPpqSyHfcvV+9S+zq1ZmzKOedjBYa5/B4bAtvIGoJUifR5YjCpQ8j9/DMrXVFZWMm3aNL755puGHorn9K6/Xld/o5Kh+7/qb0wNjRD2WjY0fdOnJrDNJbDtO9tydybxLmfB+a/D+bPrZ3wNjSuBHZHYtLq3ORIRr4p3xLeDHufBBW+oySyoaGgt3qG5uMo0XPmwu411bstZV5P4toXw9a3EhQWuIuptadKnpJSfW96/h6qbOlBK+YA/B+crgoKCmDNnDuPHj2/ooUBCe9tDWd+96+BqOLxGvR94lS1tobmQ2sf+c1NN6dLQB55t/kr9DYmEGF33rcQOtvdxbaHfJU2zXKsrXAnspj6J0yME9L0IupypPufus0WKN4cYBj2OtcSDw1T1x7YOJU3rqmF3GwvnvUReGZXVb9ww1Gq6KqXcL6UM3CIlgY4WLX7gLxX5CaqpA6h6woOudr1fUybFQWA39YdzcnfVOhBsaUtJnew1SL0PuxmY++wIc2Hybcr+a3cEuZi4e0r5a4q40rAB2gxy2M6Plc4CBLcCWwhxQAixv7pXfQ62yaAJ7CqzioouPAYbLSnlPcY3/h7HtcHJJN7EH85BIdCyp/0yLaVLIyIB2p0KCGjvt+Z4gYkrDVtL32lOaMJJj97y0hxw5cMGNem1267pR4l7KpxyhYd1BnWh/Wnqh2guVX7sY1tt9bNPvt7zvk2VFl2VxllpqfzV1AU2KLP4kbW2z3r/tcbUb1TqSnObxLkKqtKn8DQXgh0aEJpCILZNw4yloXAVJQ4uBHbTt0K5FdhSyiX1OZBmRUgEpJ8Gu36BHYtAWlwmqX2aX0CJRlAIJHeDTEvkfFP3YYO9HxvsI8Q1goKbn7AG1xq2K/NwU8dRw05IV/dEcyK2tf1nTcN2/L00A4HtbaWzMCHEY5b+13mWZWcJITy1wDTwhGYWzz+oKnwBnDyt+eRWuiJV13qyqfuwAVr1t//srlBKc8SpsYUJzn26YcbSkAQ5aNgJzcwcDmrCqk/jKrPUDnes0+AYnNYE8Tbo7HmgN6q9ptZSZhMw3R+DahZYy5RaiEiwFU5oruhbbTYHk3hKTxC6DBJHH3ZzRm8G7XI23L4VOo5ssOE0GI4adnMLONM461Hbe0fLlEYzsDx4K7DPBy6TUq4AqgAs7TVbe9zLwD0tuoC1cQsw4MrmUcHKE/qoz/j0hhtHfRESYfPDRbVsupXdaoO51PY+Ih5iUtxu2qRx9GEntG+QYTQ4LbrANT/DxNnQfZxtea8L1F9T86gK6e2UpNxxWyFEMnDC9eYG1SIEdBkDq95S5r5B1zT0iBqetierak+yytYzuqnT7Rw4tgk61ap5T9NFr1nGt3O/XVPH0W/fnK9F25PVS8+451XsTzNphOKtwJ4PvCuEmAkghGgFvAB87KdxNQ9G/EeldHU+QwWTGMDJ1zX0COqXUbNUVauUXg09ksCi5wT440WoKIZTb2no0TQcjibx+LYNM45AJSIeTr+9oUdRb3grsO8DngY2AJHADlS50of8NC6fopUmHT9+fGBUO9OISYVLPmzoURg0JKYgSOvf0KMIPILD4AZLI5jmHIjpaBKPMwS239i2ELYtDOjSpEJKWf1W+h2UKTxL1nTHBiQqKkoWFRU19DAMDAwMasb2H2HeRbbPD+Y27wlMPSCEKJZSBmTIuVsNWwjhKRwxRlhuGinlbl8PysDAwMAA57QuQ1g3azyZxHeiUrgEtlQu7W7Ra9cBaz4wMDAwaNTofdiOwtug2eE2rUtKaZJSBkkpTcC1qACzbkA40B2YBxihzQYGBgb+Qu/DjmrZcOMwCAi8DTp7BOgipSyxfN4hhLge2A6844+BGRgYGDR7zOW299GGwG7ueFs4xQS0d1iWjmEONzAwMPAfcbpGHycZ/ZiaO95q2M8Dvwoh3gYOAG2Bf1uWGwCLFy/m888/56yzzmLcuHGYTLVqNW5gYGBgI74tTJoL+Ydh4FUNPRqDBsYrqSKl/D/gKiAFOA9IBa6WUjbDavz27NixgwkTJjB69GheffVVJkyYQL9+/Zg3bx5ms/f9eysqKvjpp5+YPn06N9xwA5s2bfLjqP1HRUUFy5cvZ/ny5Wzbto2srCwqKyur3W/btm3cc889pKen06NHD9asWVMPozXwNyUlJdQ0A/Szzz7joosu4vbbb6esrMzlNvn5+VRVVfliiAFDdna262vVZzKcdisEmBJQVFTEK6+8wvbt22u034EDB8jJyfG4jf46HDlyhM2bNzs9RyoqKvj555/Jy8ur0fc3aqSUTf4VGRkpfUlVVZU8fPiwvO2222RwcLBERc07vTp27CjvvvtumZGRIZ988kn5/PPPy9dee02+9dZbct68efLzzz+X8+fPl9dee61MTEy02zcoKEjecsst8sSJE3bfffz4cblgwQK5ePFiWVVV5dPzqisrVqyQffr0cboOQgjZpUsX+Z///EeuWrXKOu68vDw5Z84cOXToUKd9wsPD5XvvvRdw51gbysrK5N69e+Xy5cvl/Pnz5dy5c+W7774r582bJ+fPny+/+OIL+cUXX8jPP/9cfvrpp/Ljjz+Wv/32mywuLpZSqvutqKjI5bHz8vLksWPHajSe3NxcuWPHDpfXdtu2bTI/P9/lfocOHZLz58+XBQUFUkopCwsL5fTp0+Wtt94qv//+e7lz507rtlVVVXLWrFnSZDLJCy+8UL766qty0aJFMi8vTx4/flzecMMN8qSTTpILFy60ux9mz55tdx+MGTNGZmdnW4+7ZcsWOWzYMAnIyZMn++T+KCkpkStWrJBlZWW12v/EiRPysssuk1OnTpVZWVlut6uoqJCvv/66nDdvnqyoqLBbd/PNN0tAnnzyyXL58uVO+/7999/yzTfflIWFhXbLi4qK5HPPPSffe+89t/fI4cOH5dNPPy1///33as9l3bp1cvz48bJr167yvvvu87jttGnTJCDj4+Pl+vXrnb5z6dKlsrKyUkqp7ofVq1fLH3/8UQYHB8v09HSZk5MjzWazzMnJsdt39uzZsmXLlnLq1Kly48aNMjIyUgKybdu28sCBA3Lfvn3ylVdekWlpaRKQ55xzTrXnVROAIhkAcsvVq8EHUB+vugrs7777Tt59991y0qRJsn///jI2NtZJIP373/+WW7ZskU8++aRs2bKlWyFe3Ss0NFSaTCbr56SkJPnMM8/I+++/Xw4ePFgKIazrevXqJd95550aP2iKi4vljh075K+//irfe+89+dhjj8np06fL8ePHy9NOO03269dPduzYUbZs2VImJibKk046SV544YXynnvukW+88Yb8448/rA9tKZUAuPHGG+3G5unVsWNHef7558uIiAindQMHDrQ7/759+8rZs2fbfV9FRYVcuXKlfOaZZ+Sjjz4qP//8c7l582ZZXl7u9pzLy8vlp59+KkePHi27dOkiZ8+e7fTQrAtVVVWytLRU7t27Vy5YsEDef//9cuzYsbJVq1Z1uhdOOeUU2alTJymEkJdffrk8evSolFLKvXv3ymnTplknjD169JBTpkyRV111lbzlllvkqlWrXI4zOztbtm/fXgKyf//+cs6cOdb755FHHpGATExMlF9//bXMzc2VUko5c+ZMuwnlwIED5f79++XMmTOdxjx16lS5fPlyOWvWLK/PMzw8XLZt29bjNu3bt5fPP/+87N69u93y22+/Xd5yyy3ynnvukTk5ObK4uFh+/PHH8oUXXpB79+61O/eNGzc6LZNSysmTJ1t/T3fddZe86aabZKdOneR5550nf//9d/nWW2/JhQsXyhUrVsjPP/9cPvTQQ/LJJ5+UK1askHfeeaeMj4+3jqdr165y586dMi8vT37wwQfyiSeekPfee69ctGiRvPjii63bde/eXa5atUp+9dVXsm/fvnbnZDKZ5LPPPivXrFkjJ06cKHv37m13/O3bt8vt27fLmTNn2t1fHTp0kP/73//kK6+8Ih955BHZsWNHmZCQ4PT7ev3112VZWZn88ssv5cyZM+X9998v58yZI88991yn6/7BBx/I//u//5MvvPCCzMnJkZs3b5ZSSpmVlSXDw8Pttm3Xrp0cPny4vOqqq6y/7dGjR8sPPvjAem/pX2PHjpVt2rSRkZGR8ttvv5WZmZly+vTpHu+Da665Rnbo0MFp+c6dO302uTcEdgO/wsLC5HXXXSe//vprL/5dztxwww1ub6ARI0bI1atX221fXFwsX3nlFdmlSxcZEhLi1YN5/Pjx8r333pO5ubly3bp1cuTIkV4/8Fq3bi0feeQR+e6778qPP/5YfvHFF/K7776zE8bjxo2T/fv3l0lJSbUWIPqXEEJ27txZTpo0ye6hERYWJh966CG5cOFC+cEHH8gXX3xRzpo1S55yyikex3/vvffKbdu2SSmlXLRokdM4Y2Ji5LXXXivHjRvnNGHSXsHBwbJHjx5y0qRJ8v7775fz5s2TK1askP/9739dCs7evXvLn376yev7oKSkRG7evFkuWLBA3nbbbfKkk06SLVu2lNHR0XaTDH++evXqJSdNmuTRsqO9LrjgAnnHHXfIq6++Wn7yySdu7+XBgwfLm266yWl5SEiIHDhwYL2cl+Prrbfekueff36djhEeHi4zMjJkdna2vPHGG63LTz/9dPn222/L8vJyuWHDhgY5v8b8mjFjhpwwYYJPj+nN/Vzda8yYMbV6vlvZ+r2UX90i48IolQEgt1y9alyatDFS19KkzzzzDPfddx/t27enY8eOdOrUiY4dOzJgwABGjhyJqKb6UFVVFWVlZZSWljr9LS8vp2vXrsTFxdntI6VkwYIF3HHHHezbtw+A/v37c/bZZzNmzBhWr17Niy++yJEjR2p9Xnri4+Np06YNycnJxMTEEB0dTXR0NEII9uzZw65du9i7d69Hf/SoUaOYPXs2Xbt2dbn+wIEDfP7558yfP5/t27czevRorrrqKs4880yCguwTDo4fP84bb7zB7NmzOXDggE/OESAhIYGoqCgOHjxoXTZu3DiGDBlCYWGh9VVQUGD3Pjs7m8OHD+Pt70UIQbdu3Rg4cCBdunShdevWpKWl0bp1a+Lj4zGbzVRUVFBRUYHZbEYIQVBQkDVYccuWLSxZsoQlS5aQk5NDUVGRS7/fuHHjGDRoEIsXL2b3blV08MiRIy7jJy6++GI+/fRTpJT06tWL/Pz8Gl3bK664gg8++MBp+dy5cwkKCmLevHn89NNPduuefPJJHnnkEZKSkjCZTISGhpKWlkaLFi249dZbWbFiBXv27GHz5s2cOHGCXbt2ce655zJ//nyEEHz00UcsWbKERYsWWc+vU6dOPPfcc1x88cWUlpY6jceRkJAQKioqnJa3bt2aY8eO2a0LCwtz6zd3R3JyMlOmTMFkMvHMM8943LZ3796ceeaZvPTSS3a/pf79+/Piiy8SGRnJ4MGDvf7uuLg4Lr74YjIzM/n666+d1p1//vnk5+czatQojh8/zjvvvMP+/fvdHq9Tp07cfffdHD9+nPvvv99uXdK5MwhP78+xzx6i4vheRHAoXc+ewsjuqbz+7GN227Zr146ysjKOHj3q9bloDBs2jBdffJGBAwdal51zzjn88MMPtrEkJfH+++9zzz33sH79ekD95o4dO0aLFi1q/J16Ark0ac3UcRWk1qqhZxk1fdXVJF5SUiLNZnOdjlGX7/7rr7/kkSNHnNaVlpbKt956S/bo0cPjzDM+Pl726dNHjh07Vk6bNk0+/PDD8q233pI///yz3LJli5252RPl5eVy+/bt8vPPP5f//e9/5YQJE2R6errs2rWrfOedd/zibzabzfKrr76S55xzjhRCyDZt2sgpU6bIt956S+7Zs0fm5eXJP//8U7799tvyrrvukuPGjZMdO3Z0Ms8PHjxYvv3227K4uFgWFxfLxx9/XEZHR9dJIzj11FPltGnT5MyZM+WsWbPkY489Jl955RX5+++/u/UD1+U6aBpnRESE/Pe//+3W7L1r1y55xRVXyJCQEBkeHu7keggLC5Pbt2+XZWVlcsaMGdZrFRsbK1evXi3nz58v7733Xjsz9bvvviulVKbQJ554Qj799NNy+vTp8sMPP7T77p07d8oBAwZIQN59991SSlmj+6KystLl9mVlZfLll1+WF154oVy7dq2UUsrMzEz5+uuvy0WLFsk33nhDhoWFSUDecccdctGiRfKkk06yO+8hQ4bIO++8U6ampjr9PydPnmz3fV999ZV13RtvvCGXLVsmv/vuO7l8+XJZVFQkf/zxR3n22WfL5557zm68H3zwgdV9MHPmTFlQUCB3794t77vvPvnJJ59Y3TBr166VEyZMkAMGDHD6P86YMUOCsmJ99tln8r///a+cPn26LCgokNddd50EZdrWWwwLCwvl1KlT5eTJk+W0adPkRRddJHfv3u10HbVnxt133y1ffvllmZ+fLz/55BM5ffp0q4VLSnW/vfbaa/KZZ56Rn376qewxeLhM/8+3Mv0/38rk81VsQr8bnpfp//lWnvrEL/L7H36UGzdulPv27ZObN2+2XpOsrCzZpUsX67VMS0uTU6dOla+99pp88skn5a5du+Rrr70mJ0yYIB9++GG5ZMkSq9/7iSeekKDieXbt2iXvvPNO2alTJ9mnTx/rPfDNN99Ik8kkk5OT5ZQpU1y6PGoKAWwS90rDFkLEA/8DJgMVUsooIcR5wMlSyvs97hwANPXmH1JKMjMzKS4uttPgExMTad26NdHR0Q09RJ9QWVmJyWSq1qIBUFxczPbt29mxYwedOnViwIABTttkZmZy//338+6772I2mwkNDbWzLuhfsbGxdOjQgc6dO9OpUyf69u1b79e1qqqKlStX0r17d+Lj46vdXrte+/bt45xzzmHbtm0APPhgBg/8978EmdR13LhxIwsWLGDU2Am8vq6Uwe0TmTGmC2azmc2bN2MymejduzeLtx7jxV92cOsZnRndPcXt95rNZrZu3UqvXr28+l/5ikOHDnH8+HH69+8PqPOfO3cuc+bM4eyzz+bBBx8kNDSUsrIy5s6dy2OPPcbhw4eZOHEib775JklJSXbH++mnnwgNDWXkyJHWZT9uyuR/v+3ijjO7MrxrMgDrD+YSHxFKu6RIAAoKCjh48CA9evSo1XmUlpby7LPP0rdvX5fdBSsrK50sUrWluNzMsfwy2rfwrFC+tWw3D3+7BYCWoRU8PrEX135qiw5/enJfLhrkupPYkSNHeOedd7jooovo1KmT12OrrKzkf//7H+3bt/fYZdFsNmMymXyWShvIGra3AvtjIAd4GNgspUywdO1aLqXs4ucx1pmmLrAN6kZ5uaomFRradGs1nzhxgltuuYXQsDDksOtZezCf+TcMpWtKjHWbp3/Yyv9+24UQsO7Bs4gND7E7Rvt7vrO+3/vkv+pt7P6ioqKCnJwcWrb0voKY4zX4a/cJLp7zJzHhwfw8cwSpceEe9g4sqqok57z4O9uPFvL2VYMZ1c39dfj32yv5bdtxAAamJ5AWH8E36w5b13dNieanmSPc7r/+YC73fL6Biwe3Zeqp7X12Dv4gkAW2t1OSM4BbpZRHUKYNpJTHAaNWXhOhqkqy81gBB7KLG3oo9U5oaGiTFtagfH7z5s3jyRdeY8mOE+SVVLBgzSG7bf7YmQWAlJCZV71vGNR2S3cc52BO47tvQkJCaiSsXTFvpfIHF5SaOeWJX3jsu81UVDaO/PDs4nK2Hy0EYMZH/5BVWMb7f+7jaL79/760opI/d5+wfj6UU8LWI/l222w/Wkhllb3ytzeriI9X7qeozMwlc/5k85F8Hvy6cdaXCBS8rXSWB7QArBFOQoh2+s8GjYdycxXZReUczith9d4c/tqTzd97s8krUcE33VJiOLt3Kmf3SqFnq9h6NWsa+JdDuSXW9+sP5lrf5xVXsP6QrQBFZl6pnfZdWuEcbPjP/hwuen0FFZWS0CATv989KmA1TCklK3adIDUunI7JNXdlSCk5XmgfjJarE3gabyzdw6mdW3jUVgOF4wW288kvNXPLvH9YsfsEH/65j+9vPR2TxWVyMKeY0grbJCQzv5SjBc7HK62opKDUTJm5kvSkKMa/vIyCMjNbjuRTXF598SSD6vFWYL8JfC6EmAWYhBBDgceB2b4YhBAiDOUjHwMkolp73ielXOhme4FqSHIVEA38A9wkpWyW07eqKsmWzHxyiyvIL6kgz/LKLi4np6ic7KJyTlj+ZheWU1DmuQLbtqMFbDtawEu/7KBjiyhmTxlo9/A2aLwc1gnsDQfzqKqSmEyCFbtPoPeOZTpoWXpBD+qee/DrTVRUqp3KK6tYfzCX1LhU/w2+DizdkcWVb60kLiKEZf8ZRYyDuR/AXFnFn7uz6ZkWS2KUzeIipeTKt1aydEeW3fb/7M9la2a+42E4UVjutCwQOVZgPwFZYdGit2YWcNIjP3NmzxRuHNmJ9QedK4lp90r31Bi2ZirpfSSvhImvLqfcXMV715xsfc68u2Kf3b6VVdIaP9GkyIhLBOYCZwFZwL1k5M3zYr9fgVFACBl5Hh/O3grsp4BS4FUgBHgLeB140cv9qyMYVaN8BLAfOBf4VAjRR0q518X2FwJXA8OAfcCjwPuAc2RRM+CmeWtYuDGz1vunxYUzpGMSg9snkl1Uxg+bMtl4SD2IdmcVcd+CDXw2/VRfDdegAdEL7IIyMzuOFRIWbGLpjuN22x3MKeGG91ez90QRU09tTysHzfn7jUecHuQnigJXUC3fpYRRXkkFGw7mcWpn59Sfd5bv5dHvtjCgXTwLbjzNunzLkQInYQ3w8d/7rYJr0oA2fL5GpQqWlHtfkrg+WLknm1lfbODKoelcPiSdl37dQXxECFFh7h//eSUVfLb6IJ+tPuh2G4CeabFWgb1mfy6FFiF9yZw/3e5TXG52OWFqArwKlKNKePcHviMjbh0Zee4VyYy4y/FeDnu3oSXU/QXLy+dIKYuADN2ib4UQe4CBwF4Xu3QAlkkpdwMIIT4AZvpjbN5irqxi+9FCuqXG1OvsscxcyU+bXec6hgQJEqNCSYwKIzEqhMSoMJKiQi3LQmkRHUrv1nG0SYi02+/m0V04kF3MMz9t46u1h1m1L4cNB/Po0ybO5fc0dcrMlRzLL6NNQkSjdw8cdtCUp8z9y0nTAnh72R6rhnTvgg3Ehts/Kn7a5HzPnSh0nb8speSvPdm0jo+gbWKky238zXadDXedG4H96HcqCnrN/lxKKyoJD1GR2N+sP+y0LcCPlmtgEnD3Od2sAtuT+VdKWe/30DXv/k1BqZkHvtpEalwELyzaAUCL6LrHbfRsFcsCVCyEt3EPxeWVTU9gZ8RFAZOA3mTkFQLLyIj7GpgC3ONmnzjgQeBKYIU3X+OVwBZCjHazqgw4KKXc52Z9rRBCpABdAXczk4+Bi4UQXYE9wFTgBzfb1guPfreFd5bv5YpT2vHoxD719r07dMEeM8d05cyeKcRFhhAbHkx0WHCtHw5tEyO5+5zufLv+CJVVkneW7+XZi/o5bVdZJamSkpCgwGpMUFeKysws2X6cHzZm8uvWYxSWmblldGfuOKtbQw/NI0u2H+erfw4xfWQnurhwYxzKtX+ouhLWgJPbJL/U/rNmPu3XJo7dx4soKDO71bAXbzvG1e+sIjkmjGX/GUVYcP135d2WaRPYGw7lVrv9nqwierSKRUrJt24EtkaHFlG0iA6zfi5x4e8HNVma9Npy2iRE8P41Q6wTAn9ToPvfLdNZUrLqaLqPCAkiPckWTH3ES4F97otLOad3Ko+dX3/PyXqgK1BJRp6+E8o6lNXYHY8DrwFem0e9fcrOBRZaXh/o3n8M7BRCrBZC+CS9SwgRAnwIvCul3OpmsyPAUmAbUIIykbvVsM1mM4MGDbK+5syZ44uhWimtqLSajn7adNTrali+YKvuQXRO71R6psXSOj6CmPCQOs/kW8dHcHYvlW/7zbrDZDloUMXlZi743x/0f+gnVu/LrtN3BQJ5xRV8vvog1723igGP/MyNH67h63WHrWa+137bxe7jhdUcpeGQUnLX/HUs+OcQV761kmMFzg9QRw3bkaQoe61rQv804iKctSEtYKlnWiyJFk3Nne/28e+3WvfZeaz+r19BaYVDsJ29Kf+XLUe51MGEq43zeEEZB7I9X7MerWIJMglCg9XjtMSNhv39hiMcySvl7705zF22hw0H8/h+wxHrhNvsh+hyx2fRz26scd7QOj6CLi1tAXstY8OICrVNOjLzPF8njRNF5Xz41367SVSDs+pteH0EvD6CFpHCaxO1jmhUcLaePMB18E9G3CDgNODlmnyJtwObC8QB/5VSlgghIoCHLAN6AXgWFTR2pqudhRC/4X6m8YeUcphlOxPKF10O3OxhPA8Cg1F9uTOBK1D9untJKZ3yS4KDg1m1alU1p1h7ft9+3PpQP1ZQRmZ+Ka3iIvz2fXq2WNIrQoNMdEz2ferg1KHt+X5DJuWVVXy8cj83j7bNy57+YRvrLA+/uz5bz8IZpzeI9lQTDueWsPFQHhWVEnNVFeZKSV5JBYu3HWPFrhOYHVJTokKDOK1zC37echRzleSJhVt548pBDTR6z2w7WmDVmI/klTL9gzXMu26I3f/kcDUP1RHdkq3pXiYB947twdCOSdyzYIPL7XumxbEts4B9J4o5UVTG9qMFLN+ZxfCuyXy97jDnn9SaUJ315UB2CSXllRzOK2Vcn1bWSGR/st0hpPlgTgknCstIig7jWEEp17zr/GzYZZmYZRdXr4X2aBULQGRoEOXmKrcmcb2m+38/buP/flSFbF657CROFJbz+PdbeGBcT644Jd27E/OCnGL7kqyHvdSCXZEcE8YD43rw77f+pqDMzNCOSYTrBLa3GrbGsYJSuqUGSDDroKvUC8i6QTgHIWTE/YYHGQbcAsQ6LI8FnGclGXEmlLycQUaemQzvXY3eCuwZqJKkZgCL0J4FHJZSPiaEuANwG50gpRxZ3RdYIr/nohz250opnYv/2ugHfCKl1L7zHSHEC0BPwH+S2Q2OAV/rDuTWm8DWolQ7t4z2i1n65A6J9GgVy5Yj+bz/5z6uH9GJkCATK3ad4J3le63b7T5exOzfdjNjTODV0ZFSsnJPNm//sZefNmdSVY0BJD4yhDN7pHBO71RO69yC8JAg7vh0HZ+vOcjPm4+yYtcJhnZK8nyQBmCZQ2DU6n05PPPjNmb9qyegzPy5lgd4l5bR7NBpu6HBJu45p7td6tLJHRJJjQvnwkFt+WTVAf7Zn+v0nT1bxfK7xRy8LbOQya8ttzOfrzuQS6Tuof7XnhO8/cdeAGLDgxlZD+lP2zKdtfoNh/IY2a0lz/7oupfzruOq0FJusafHkKKnJrBDgsilwq3A1tImHdl4KJ/ZS3YBcP+XG30isMvNVRzJKyHbh4GAyTFhDExP5PsZp7N42zHG9U2zy9l2zCyoDsfJRECTkTfS8/q4KCCYjLguZOTtsCzth2u3biwwCPjEIqy1H8hBMuIuJCNvqbuv8fYJX4TSaPUMBDRt1he2nNeAHsB4KWV1tpW/gQuFEClCCJMQYgoqen2nD8ZRI8rMlSxyMDP9cyC3Xr5bSsmWI2oCp83yfY0QgqsslYmO5pfxw8ZMCsvM3PXZOgDCQ2ya/auLdwaUyVhKyVdrDzHu5WVcPOdPftjkXlinxIZx5dB05l07hFWzxvB/F/bjjB4pVj/jXWd3IzxE/Vwe/W4zWYVlPP3DVl5dvNNljnJDoBU+aRUXTp/Watb+7XpbqYQjOu16bJ9Wdvtuf3QsVw/rYDfpG9NDuUOCTIJ5157Ci5f0t9tHCJXWowUvZRWWOfm6F287TqnZdn3e16X4fLTSfRMKX7It02aF0lh/MA8pJd9vcF1KYpdlMpOr07CnnJJOvzZxfHnTaXbb9kxTv70Iy8Tk8zUHmTL3Lx7/fotdISK9j79Xmu33esRLU7I7tPtcP2H799srGfF/v/H8oh0e9qwZLWPUxKxtYiRXDm1PYlQoETo/vDeTGz3emtAbBRl5RcAC4GEy4qLIiDsNmICyGDuSB6ShIsn7o7KiQMnUvzx9jbca9n+Bn4QQX6PSr9oA41FmAFCV0D7z8lhOCCHSgetRQWyZOt/r9VLKDy1FWjYDPaWU+1FpZi2BtUAUSlBPklLm1nYMtWXp9ixrgI5JQJVUWkV9cLygzDqD7tHKf6al8/qn8cTCLeQUV/DO8r38ufsEB3PUj+2ec7rTt208k15bTnllFbO+2Mi864YERDT1W3/s5ZFvN1s/hwaZOK9/GhcObEN8ZChBJkFIkCAkyERqbLhH82xqXDjXD+/Ei7/sYNPhfE578lfKzGqe+vnqg/zfhX0ZmJ7o93NyR7m5ir/2qDiCYZ1b0CYhkg2H8jiSV0pRmZmosGC7gLNTOyXx0i/qYa7/V00a0Jq5S3fTIiaMS09uZ10eERrE2b3sc6zbJkQSFRZsl7Psiv0nbEJL73KoLyvUNotJvHurGHKLK9ifXcz6g3kczit1W5Ngd1Yh5soqOy1w+shOpMWrMU/sn8aXaw8TGx5sFWSRobbH6dIdWSzdkcU7f+zl6mEduHl0Z2sUvZY2dtHsFazcm11tXEF1/Lr1GDM+XgvAn/eeQViwyZrG9vv24y73GdUtmcXbXK/Tb1NSUcmfu9V95SqWISK09i6wmprQGwE3olKejwEngOnWlK6MOKsMIyNvP/pAswxrzuRRn+RhSynfE0KsQoWtpwHbgaFSys2W9d8C33p9Ws7H3we4fVpahHS07nMpcJPl1aBoM/TwEBPn9Erly7WH2XAwr16KA2zWlQf0l4YNEB4SxKUnt+N/v+1i9b4cVu9TbR6HdkziyqHtMZkEVwxJ5/0/97Fi9wkWrDnEpIFt/DYeb9ibVcT//aiCnRIiQ7jqtA5cNqSdXTRvTbl+REc+WrmfYwVlVmENKld98uwVXHVqB+48uyuRocEUl5s5lFPCsYIyeqXFEh/p39Kn6w7mWk2xw7q0sLv39mQV0bt1nJ1gaJMQwd3ndOPtP/by1CRbtG56UhR/zRpDeLCJYAcXS1iwidAgE+WW4KiuKeonmRTl+Zo6at0aVfUQnCmltAY3dUuJoaSikv3ZxWw4lOvk29ZTWlHFmv255Og07ATd//DB8b1Ii49gdPeW1slphIuo7/LKKmYv2cXOY4XWyXWS5R5Mi1fP6cO5dRNcc37fbX2/7mCunSXBHZMHtq1WYE8b3om/92ZbBbZjfAfgMtI9IiTIbaS8Hm/TwBoNGXnZwEQ36+xkmMO6vXiQf3q8joazCOfN1W7YjCgzV1qjLkd3b8nQTkl8ufYwReWV7DxW6PeACn2EeHc/f9cVp6Tz+u+7rRGtUaFBPD25r1UrveucbvywKZPjBWU8+t1mBqYnEBwkKK2opLSiitbxESRUo4n5iqoqyT0L1lvLKc6+YiBDOtbd5xwZGszDE3px87x/SI4J47/jenIgp5hnftpOubmKt/7Yw9frVMCWPmWma0o0P9423K9WB7059NROLezKTu46XmgnsIWAlNhwbhzZmRtHdnY6VrSbghpCCCJCgygvUde1c0t1zyU55PO+e/XJPPTNJnYf99xwpz58mMcLy6zf0y01hiop+Xb9EY7ml7F0u3MxFD2/bTtmFVJhwSY7bTIhKpS7z+lut72jtnnRoDasO5DHtqMFrNiVZRVumgtB09Zr6vt1RO8zzy+pcNkP4OQOiZzRvSVPLNxKu8RIzuqVwiMTerEnq5jwEBP/+22X0z7JMaG0TbRZQVJjncvOupqk9G4dy997nXu3O1LX826OeJuHnQjcibK3280SpJTDfT+swKO43Gxn8gL1kNRMamN7t6JLiu3SrDuQ63eBrUWIt4wJs87a/UWaJcXr+w3KknPfv3rYFcGIDQ8hY3wvbpq3hpziCkY+85vd/jHhwXx7yzC7vE1/8dHf+61awRWntPOJsNY4p3crVt2fREx4iFWLPaNHCnd/tp7V+3Jc5rZuP1rIgewSa/tFX1NZJfndkl/bPTWG5JgwO6GrCU4ttSklJrzWAYr6wCktxcfRajGiazKf3XAqAx752eOxcmoZELXjaAGPfb+FyQPbMK5vGvmlFdz04RpKKyqd8pv1qUPdUmPszvv9P/cCKshQ73/Vym0u3nac3hZfc4IXFpJIB4E9vGsyfdvEc/+XGykqr6TIIlg1i0Qri8B2bJpRUwp1Zv3DuaUuY2guHNiGCwe15exeqSRFhxISZGLK0PaA+1iCFtFhjOubxoI1hygur7RzkWiEBAmCTMLuHHq08lJg11LD3nQ4j0Wbj3H5KXWzmDVGvNWw5wFhwKfYAs2aDQeyi5n02nJuPaOLXQTndxZzeFiwidHdWxIeEkRkaBDF5ZWsPZjLRYNd94f1FVuPaL45/5nD9dx5Vjf2ZhUzID2ey1z8eM/tk8ro7i35desxp3UFpWae+mEr/7t8oF/HeCSvhCcseb9pceH8x0EL8gWO5u1OydF8ev1Q5v21jyXbs0iKCqVNQgQSeO5nFYW88XCezwV2VZXk9d938/6KvdZ0ndMsFbwiQoNoHR/BodwSdmcpga1p2Joptq50tghsfRvOIR2UHz8hMoSY8GC7VCZHcrxImXLF2S/8TpVU/tlxfdO4Zd4/1tKhS3dkcWZPW69uO4GdEkNSdBiD2yfw994cax30nq1irT5fgLN6pbI1s4AtR/KJsAQaxkdWX5nLUcOOiwhx6afXfP6tffB/MFdWcSjH5uo4lFvMWhfR/Oef1BrAZd9rzQfvSFyEquXw/jVD3H6/EILIkCC7WIDk6DBeuvQklmw7bq3+5opjBWWYK6ucXC/V8a+XlgGw41gBr1zWvKpReyuwTwWSpZSuyyI1YaSU3PjhGo4VlHH/lxvZfbyIWf/qQWWVtJrDR3Vraa3L26d1HH/tyXb5o/ElZeZKa66oPwPO9HRMjub7Gae7XS+E4MVL+vP1usNUVknCg4MICzHx9drD/LL1GN9vyGT1vhwGpif4ZXxSSmZ9sdGqcTx2QZ96K4EYZBJMGdreqrWAKqDxwqLtVEmVRnSuQ2R2XVmzP4enfrDVFooOC+aCAa2tnzsmRymBbblPNF+pZoqtK50sArt7qxg6JkdxJLeUJy5Q/nAhBO2Totig6wD208zhvLl0N5+uUg/xTYfzueLNv7hsSDuW7sji8iHt6N3ac07qwZxia6R/lVQ51kt0gVWOqVNa5HxiVCjJMWEIIXj7qpO5+u2/WblXWWG6psTYCeyT2sZb36+x/I69EdiOGnZseAidW0YjBHaNVZIcTOJ1YX92sTWmAJQ1RxOeA9rFExJk4j9ju3sUiikuTN2A1y6c8FB7gR0fGcJ5/dI4r18a32844tafXVklySosr1GHtzJdxsG364/wymVe79ok8HZqsx4VGd7sEELw0IReVr/TW3/s4br3VvHjpkyr9nBuX9uDuL/lx77taIHbike+YOexQqt/rUdq/WjY3hATHsLlQ9K5cmh7Lhrclgn9W5NxXi9rFajHv99Sq0pwi7cdY8jji6z5qq74au1hq3Z/wUmtG7zFYURoEF0sft6NOsHlK/bqoq/vHdudP+4ZTa80m8DraNGm9mQVUVklrelDresgKDrpivNoZveQIBM/zBjOmgfOtGtd2cFBm0tPiuTpyf24ZlgH67JlO7O48cM1fLRyP+NeXkZBaQV3zl/H23/scfn9js0oHAPHCkttAntPVpE1sGps71SrAIoOC+adqwfzrz6tSI0N58JBbXjt8gG0iA7jofN6udRCvTOJ2+s/cZYGG+0d3EAtrEFndRfYju09d+ly6ycNbMMn1w9lQDvPE2RHDfuUjok8f7FzGWJ3OPqx43TXqrqJjqtqfJ6orupcU8dbgf0r8IMQ4j4hxNX6lz8H5ysqKyuZNm0a33zzTa32H9AugS9vOs0a2KXSKP4BlDn8jO42waAJ7MoqyabDvn9Ia2j51+DfCHFf0DYx0prLvXpfDj9uqllnMSklj3y7maP5ZTz383byXAQrVVZJa+WoFtGhPDCuZ53H7Qt6tVb/m02H831eslZftOKKU9Kd0m404VlcXsmmw3lWE3BdBMUrlw1gTI+WzJli79oIdQjKArhmWAdr0FL/tvHWimsJHh7iT/+wjc9WH+Shbza7bCby1277EriOub96E/y7usI+/7bcfxqRocG8evkAVtyrJjlj+7Ti71lnMPXU9rRJiHDK8PAmyt9RcMVa/h+OFjBNw44ND3EZtFUTn7Y+UwTsa8A7avzucIx/+XjaUM4/yXv9zPEc4nX3YXUpf4WlZvZkFfH9hiNUeFGadW+WLZDRXXBkrdm2EL6+lbgwap+r5me8FdinoyqZnYnqPqK9rvDTuHxKUFAQc+bMYfz48bU+RpuESObfMJRR3ZIBrGa5kd2S7drU9dOZ09b6MR97q59LkvqaG0d1ts62n/phm1c/To0/d2dbA6fKzVV87aIZwx87s6xBVTPO6FJvEenV0dui8WYXldepLKQrtEjwmLBgl60S9ffF4q02s3FdBHaPVrG8OXUwZznkZLuiX9t4ltw5ih9uO50Pr7X5QT0JP33VQFf+76MOGpljfXtNYOWXVjB/1QEATu/SwmUjFLA3+2rvQ4JMtE2wv0aeJhkajhMWrcPZWT3tr5U+DS42wvn/Vm72/rfxz373wV2uJgOuqGv6abgL373Gdad3JC0u3C7aXM9BS0OUGz9cwzXvrqK4mtake0/YBHZKrI8DzrqNhfNeIq+MwKiE5AKvBLaUcpSbl7suXk2SmPAQ3pw6mKtPs5n0JvRvbbdNq7hwki0mJn8K7C1+Lknqa+IiQrjFUod8T1ZRjapczXPY9jPLg1jPJ3+rZeEhJiac1NppfUOhb0nqa7O4pmG3dPPg0pun31yqcnWFqL+YBwCTSdA9NdZuQuHJvKwXwEUuHt7H8u0F9HGHbmMFFpP4/FUHrVHZV53WvsbjdjSL1yZKXPMbT+ifZnUDtE2MsNM6XWmJ3grsqippjZWJcXGciFDvNdCTLcGCjpYIb9AC8zT0PumJJ7Vm+b1nsPTu0U4tWgEe/XazNT/99+3Huej1Ffy69ajbCb1eYPtcw24E1PhJLxQm7eWPQQUyQSbBf8f35K1/D+LpSX0Z29t+9iyEsJrF1x3M9csY9CVJu9fjw7euTDklnXaWVLAXF+2wPlw9kVVYxg8bbdH4oPoZ632XJwrL+Gmz0szO7dPKLmq5oenRKtZaSWyTvwR2jOugnVax4dZyqrb0w1Sn/uf1jTfaKihzqd3nMrNdChM4a9hakRatJWb7pEhGdq15LIOj/92rKHE3Gq0Qgvv/1YPPp5/K5zecaqfRugqKLKv0TsHbccwWYHZ2b2eLh7cmcYA5Uwby+pSBtcqq0J93SJAg2U2qlStTv/b/0mJcNh7K5+p3VnH3Z+tdHmOfLm6jrAaWiKaCVwJXCNFaCPGFEOIEYAYqdK9myejuKVw0uK3LSEpNYB/ILnHph6sr+pKkPQPcf60nNNjE3eeoftInisqtTSA88dnqg1bf66MTe9st1/jin0PWbS4Z7Jxu1pBEhwVbH/4bfC6w1b3lzjRoMgk6tLAvruSqUEp94627wlE4H3NRaMNRw9aEvJbCNqh9Yq06gjkKbL21wh36oDPHx4IQgoHpCbR0iMiOcaF1llV4J4jW6Mzhw7smO6331iQOyk1xdq/UWpUa1e+TGue+xK+rSmkaj0zoxaQBNr/5DxszeXHRDk55/BeufGslVZZ99+h82IFSw78+8VZDno1qeXkGUAgMAL4GbvDTuBo1/drEW9/7Q8veYlfhrPEIbIB/9WllLWnp2F3KkaoqaTWdd24ZzeSBbaymuwVrDlFRWYWU0moO79giisHt/ZMyVhe0RhwbD+dXs6X3SCmtwspdWg7Y+7GHd02uNm2qPvDGvAzOAvtovvPk97ijD7u0gipLuhC4zzGuDn20/SWD2zKgXXy1++g12hCTd49WlyZxL+M7ftmiMiJiw4PtUtFcjcef6AvVpHmoD+8pmC4pKoxnL+rHzaPUhLKkopLnF20nM7+U37cfZ2tmAVJKu/rjpV5ObJoS3grsU4GrpZRrASmlXAdcA9zhr4E1Zvq2tf3Y1x7wfaT4Frsa4o3HJA6apqGE7qbDedaZsyuW7zphNYFddnI7hBBcaKlRnlVYxu/bj7Nmf661TaQ7i0dDowWeHS8oc6kl1obc4grrg91Ra9PTWacZ3jiyk0++u64kRHlnEncMOnOVApRVYF98paDUTHZxuVU4JNdSYA9MT+DpSX157qJ+PHFBH6/uK72mGRLk3X3oSsP2xoedXVTOb9uUwB7b27UbqC6NOWqCXpP3lDKoTy9zrHeuXYdR3V27L7YdzcdcJe2Evr4LXHPBW4FdiTKFA+QKIZJRLTcDJ7ongIgND7Hmq/oj8EyLEE+uh5Kk/kDTOIvKK9lzwn296Q//Uq0Yw4JNVnPZuX1aWTWH+asO8snfSgMPNgm7oiGBhJbaBb4zi+ujpT1pkZcNacfA9ASuH97RWoWsoQkLDuKFi/tXu52zhu0ssB2LchSUmu3M5LUV2KAmgBcMaOP1JNDOlxvsrYbtwofthcD+dv1hq4n5/AGtiQpzFs6RIfUTlKUPevWUgfDMhf0YlJ7AHWd2taa8abhLgdPYeqTA6bp46zpoSngrsP/C1rPzR+ATVO/PVf4YVFNAS+9adyDX5/m3/u6B7W/6tK4+cvpYQam1kty4vmnEWYJ+osKCrRXDftl6lG/WqYC0M3q0dBt81dDozasbD/nGLK6PlvZkEk+JDefz6ady77k9Asr6MPGk1jw9qa/HbRyDzlyZxB0pKK2wF9j1OKG1M4l7mbkRXUsNe5WlVndyTBgnt08kOMhkDTDUqC8NO18XPNrKQ7nVdkmRfDb9VG45owvRDhMMTcN2LD6jsTWzgDKHyVmpudLnz9ZAx1uBPQVYYnl/G6qQykagmRWG8x7Np5RXUsGuaroW1QS7kqR+bi7iL7qmRltNhhsOuhbY81cdtGoQlw2xDySbbDGLV1RKq4YVaMFmeuIiQki31BHfWItiOlJKp5Kbem3T5/mo9YQrc7Aedxp2lAdBVFReyZVvrbR+9uQu8DVaYRiAEC8D3VylOpWbq9h5rNBjFoW2Tt/HXe8PDzYJa+S1v8nWNXHxthmHo2DWa9x9damQZ1nqwm/NzHfSsKX03t/fVKj2PyqECAJeRJnAkVKWSCkflVL+R0p5xN8DbKyc3MHWIWrZDs99Z2vCjqO2kqSNKaVLT1hwEF0thSxcCTApbcFm3VNjnAJ+hnRItKaHgXpouYqSDSQ0P3Z1udhVVZLHvtvMcz9vt2oP932xgf4P/8RXaw9Ztzum0yID1bJQHY5mUUecfNgWDdtV6VB31MUkXlP0ZulhXVp4tY+roLNFW44y5rklnPfKH3y8cj+f/u1cd0CbzOj31+e61yRCvK7oSzB7G1CoH7cQEK0T4I9N7EPv1rE8MK4n/S2//aP5ZXYBZxrNLfCsWoEtpawEzgIa7ZWpa2nS2tA1JdraP/b3aqKhvaWqSvLEwi3Wz3110eiNDc0svulQvlPg2Ymicg5aOhBN6N/ayZQrhLBq2QAXDmpT52pN/kaLzj6SV+qUO6xn6c4s3li6h5d+2WH1d3+08gBSwoyP11q307TN2PDgejN9+prqNOwiBw1biwjXT9Y0XAV5RYQEedTGfU1SdBj3ndudf/VtxaxzvSuN6yoP+x1LSdU9WUXcs2ADd3++npV77EuyFpYpIakX0nohWJ/3hJaq2SounH5tvctC0E9uosOC7VLB+rSJ49tbTueaYR2stfjBuW482DcDqTNNqDTp88BDQojAqUhRA3xRmrSmCCEY3lXNslfsOuGTG2vusj38sVN1FZo8sA2dvMgNDVQ0AVZQZmZftn3HVn2f5NYJroNYLhnclpYxYSTHhDmZzAOR3rrAM09a9rZMm497p66RgyOatlmfJl9f4yisnr2wHy9depI17c/RJK7VNGgRHeakQbrq+BQbEVzvfvtpwzvx6mUDrDEX1eHKh+2KlXtO2H0uLFMmcf2kRy+86yulC2BgeiK/3DGChTNOt3MLeCJSN1ZPhY70ZU6zXfRP92ngWVMpTQrcAtwFFAghDggh9msvP46t0aOZaUsqKq1BIrVl0+E8nv5RtVJMT4ok47xedR5fQ9LbQ+DZCd0PM9GNia1lbDhL7hrF73eNctlzONDorQs82+QhH1tfGGLfiWK3QTValHhj9V+Ds/920sA2nNcvjVTL/1PfyKKisspaFSsxKtRJIKW6mLh4E6TW0HhbXtMxiK3IqmHbrkOMnYZdv2U7OyVHe9UgRUNvAvdkadH/n3Nd9E9vbsVTvP2vNoomH4HGsM4tMAnVKGTJ9uOc1tk7v5YjJeWV3PrRP1RUSoJMghcu7t/o6+h2T40h2CQwV0k2HspjfL806zr9TNpTt5/GZApOiAqldXwEh3JLPGrY+gDF/dnFbn10moad0kj91+DaHAw2waNvlam3uiRFh6r/vS6W01OkfCDjKujMFY4BZFoEvT4trKE07NoQqZtoeNKw9eeRXeQchGf4sF0gpVzi7uXvATZm4iNDreldv2+vfeDZY99vtj7IZ47pwknV9LdtDISHBFk7KDnmJp9weDg3FTSzuKdcbL2GvfdEkcvuRVJKaxGRxmwSdxfFrE1G9SbxEw6TOEeB5Co6eeaYrr4Ypl/x1iSuv1Zl5kprdLQ+Paqhgs5qg17h8Kxh29a50rB96sNuBDS75h31zQiLWXxrZgGZtWiv+MuWo3zwp/I8nNw+kekBUAvaV/ROUwJs46E8O9OvXpvypulCY0Ezix/MKXH58HHMId5/opjicucHUk5xhbV2em1LbwYymhDT52E7Wl0cTb5dde0zbx3dmWcv7Me04R39PNK6462lLEjni9fM4Y776wVfoFuf9ILYU866XhPPcWkSr+JAdjHLdmSxZn8O5iae5mUIbD+jTzf6vRbpXc/9vB1QP8bnLu4X8NHQNUFrPZlfauZAdol1ufZwjgkL9jqIpTGg99u78mPvdsjXP1FU7rIcp30OduPVsN2hCaGi8kprKUo7q0tUmFP0d8uYMD67YSjPX9yPmWd2ZdLANgEvtACivPQ169tN6qPn9Vq1/lj1lYNdW/SWAU/PtMgQvQ/blUm8koUbj3DF3L+44H/LKW3iHbwC+7/aBOjXJt4a6bikhmbxrZn51gf7tcM6NnhbRF+jF2B6M7H2cE5sQuZwcH++GnpzuMbmI/apLObKqiZRNEXjkQm9aBUXztypg6zL9Jqi1hM7W5cK58okHhkWxKD2iZx/kvelRAMBbzuJ6YuG6PPT9ddKb14P9ApgeguJp2sQHGSy1h13qWGbK+0sDoHuCqgr3rbXvNPN8tt9O5ymR5BJWIsoLNuR5bFjjSOf61pIBmqd7LrQIzUW7beqF2CaSdxTwFljJDkmzCpgXQWe7T7unMa19Yi9Jl5UVmlXNKWxa9hThrZnxb1ncEaPFOsyvZlXM4vrTeIJkSFOJvHGHoRZHXqBXVTuWsPWa601ec40BHoRXV2fFM0snuNCwy6rqLLGeUSEBDUpC6QrvNWw/+tm+f2+GkhTRvNj55VUeN1u01xZxZdrDwOqsldbF8UiGjsRoUHWwgibDrvQsGuQJtJY0PzYLk3iFg1b77ff4iCwC8oq2GrRuoWo30pe9YVeCGnmX+2eiI8MITjIZGcqBfc1qBsD3tQR0NcX1/v23VU6C3B5TaXOAlCdlcHxf62n1FxpjfNw1QClqeFRYAshRgshRgNBQohR2mfL61rAufSMgRPDu+j82F6axZfuzLIGIE3SVfVqamhm4g26wLPsInXeTU3DBtv57skqsmuaADYfdp/WcSRYhPbWTPufWG5xBV+vUxO5k9sn2vUibiroTbtaLna2g9XF0T/dmB/Wj03szT8PnOkyl1xDXzNbHz2vF9j6lpWe2tYGAsE6IV2dHz/Sg/WktKLKKrAbQ8xCXalOw55reYUDb+k+vwlcjSqoEvA0RGlSPalx4XSzRLF668fWzOHhISbG9k7129gaGi3VKbe4goM5JUgpybHkWzY1HzbY+7E367RsKaXVh92xRRTtklS9bMco8YUbj1hLmzbViVyMC5O4pmEnWQS2Pn+5e2pMo85HF0KQEBXq0Zyrr+hlJ7B118Gk891XBbgP+5zeqSTHhBEZGsRNozxnvnjKKS+tqLRaYbwN4HNLIyhN6vEMpZQdAIQQ70kpr6yfIfkerTRpQzKiWzLbjhaw7kAuucXlHqsC5ZVU8JOlteQ5vVLdFphoCvSxi5zOIz4yxKpNNEmTuEOJ0lM6qiYxmfml1s5jHZOjyS2pYJ2LXurvLlc9wiNCgqxtRpsaeiFU6EbDPq9/a37ddoxereK4f1wPr4O3AhlP3bnKK20TN3dR4iad+lUZ2PKayNBgfr9rFBVVVR4Lp4BzIFlESJD1t1JmtmnYdS4W020sdBtLXtnLAZvc7W3hFKuwFkKY9C//Da1poZnFqyQs2+m5Gcj3G45YfVZNVYvS6JkWi9AFnnlb5ayxkhobTguL5UAfeLZHl9LVMTmKdDcxC5oAO6d3apMNtPIUdJYYpXz2nVtG8+0tp/PU5L5NZkJbUOZcJEej3E2UuF6rHJieaH1/czVaayAQERpUrbAG+0kJKKtjmCVtrayi0hqE57hdU8TbKPEBQogVQogioMLyMlv+GnjBoPYJ1plidX5szRyeGhvOqZ1qV860sRAZGmxtYrLhUH6TrXKmIYSgl9ZqU2cS36VL6eqgM4m7Y9KApjuRi9GV2ywoM1NZJa0pPUlNcBKnobdiO5Ys1QtsTcOODLWPio6LCGHxnSP5fPpQTu6QSFPB0TcdFhxkFdilFZXW9p5NPaULvI8SfxdYDAwCOlpeHSx/DbwgPCSIUzqqH9GS7cfd5knuO1HEqn2qUcjEk1o3+TQF0LfazCO7UK9hN70IaLCZxXcdL7SmpGgadliwibS4CNonuc8KSI0NZ2inJLfrGzv6ALLCUjO5xeVWYdYUrS6ucCw5q0/rctULW6NDiyg7Tbsp4FgkJyzEZA22LK2oMjRsF6QDs6SUW6SU+/Qvfw6uqaGldx3NL2P7UdetEz9fc8j6flITzL12hRaIdaKonM26NKam6MMG2wRFSlvg2e4sdT90aBGFySRo50FgD+2U1KQncsFBJqu2VFhWYecmaYpWF1c4WhLKvRTYTRHHlL2wYJ3ANldSXOYjH3YjwFuB/QVwlj8H0hzQlyldsv2Y0/qqKsmCNcoc3rdNnLU5RlNHH3imj6JvilHigNUkDjY/tjVCPFmZwpOjw9w+gPq1iXO5vClhrSdeZnZq/NEccNQWtUDMBWsO8u36I4D3jUMaO65M4uEhmg+7eWnY3p5hOPCFEGIZkKlf4avocSHEB8AZQJTlO56WUr7pYfuZwH+ACOBzYLqUMqAb4HZoEUXbxAgOZJfw+/Yspg3vZLf+773ZHMxRNbWbso/SES3wTEr4Z79yB4QGm5xMYU2FNgkRxEWEkFdSwcbD+ZSZKzmQXQyoewSUr7tdYqRTHjZAX0sHuKZMdFgwxwvKKCyrbPKBiN5QZq6i3FzFQ99sti5LrybOoangZBLX1UnPLiq3tthscB92RlwiKu35LCALuJeMvHketu8IvASMAMqAt8jIu9vTV3irYW8GngL+AHY5vHzFE0B7KWUscB7wqBBioKsNhRBnA/egBHx7lC/9IR+OxS8IIazR4iv3ZDu1Tvzcol2HBAm7/tBNneiwYKug0uo9JEWFNqqa0DVBCGG1Kmw8lMf+E8XW8+7YItq6Xbobs3jPVrEulzclNHPvN+sOc+OHa6zLk5poXAPAh9cOoWNyFI+d39spxqXMXMUfO7PIK1FxvgPaxfPQeb0aYpj1jmMZ2rAQE70snf5W7s22Lg+A4jmvAuVACnA58BoZca7/SRlxocDPwK9AKtAG+KC6L/BKw5ZS+l0YSik36T9aXp2A1S42nwrM1fYRQjwCfIgS4gHNiK7JfPjXfsorq/hrdzajurcEoKS8ku83KOPFqG4tm50m0ad1nF23qoQm6r/W6NU6lmU7s9hxrNDOb6+ZxMG9BtUUq5s54s4/mxDVNFK4XHFa5xb8esdIABZZ6jBolJur+Ga9qnAXFmzivWuGNBsftrOGHcQ1wzrw8d8H7JY3aHnajLgoYBLQm4y8QmAZGXFfA1NwLZf+DRwmI+853bL11X2N13nUQohQIUQfxxKl3u7v5Xf8TwhRDGwFjgDfu9m0F7BO93kdkCKEcBk6azabGTRokPXVkEVUhnZKspbl0/trf9qcaQ0maeq5167onWbvl23qwUXa+VZWSRZusHmZPGnYCZEhPHdRv/oZYAPjyj/b1NqteuLuc7rbfS43V/LLFhX3Mqpby2YjrMGVD9tEl5QYu5ggqKOGvepteH0EvD6CFpGiNhe3K1BJRt523bJ1KFnlilOAvWTELSQjLouMuN/IiOtT3Zd4NTAhxDBgPhAGxAL5QAxwAB+mdkkpbxRC3AIMBUai7PquiAb07Y609zHACceNg4ODWbVqla+GWSdiwkMYmJ7AX3uy7fKxP7PkXidEhjCqW8uGGl6DoS/ZCU3fV6kPtPt1m3oQJ0WFEqdr/JGeaNOwLxnclicu6NNk3QSOxLgQSE01CNEVPVrFsuj2EWR8vYllO7M4mFNiTe0a3IRyrL3BVZQ4qJK0+mdoREgdJjGDrlIvIOsG4b6CjXscZRKWz+4ih9sAo1Du31+AGcBXZMR1JyPPuY+oBW/P8HlUENjzQogcKWWiEOK/QLE3OwshfkM51l3xh5RymPZBSlkJLBNCXAFMRznlHSlETRw0tPeNohnJ8K7J/LUnm91ZRRzILiYkyMQflupnE/q3Dvjm8/6gV2t7v2xTN4m3S4wkJiyYgjKzNWVH8+Nr6DXsyNDgZiOswbWG3dQncY50bhlN6/gIwD4Pu1Ny8wg203BlEgdoFWefq+5XH3ZG3G94kGGovhqOwSWxuJdJJcAyMvIWWo7/DKr7ZQ/srcd2eCuwuwIvOix7EtgDPFPdzlLKkV5+j55glA/bFZuAfsCnls/9gKNSSiftOhAZ0TWZ//txG6DM4gWlZmvQUVPse+0NseEhtE+KZO8JNQdsyhWtQLUU7JkWy197bEEzHR0exGnxEUSFBlFUXkmLmKZ9PRxxZfIN8H4WfsHV5F2rDNhccDKJW1K6WsVF2C33qw87I2+k5/VxUUAwGXFdyMjbYVnaDyWrXLEeOK2mw/BWlcvDNns4IoToCSSgzAB1RgjRUghxiRAiWggRZIkCvxQVQeeK94BrhBA9hRAJqJnJO74YS33Qs1WstZ70ku3HrbnXXVpG25lKmxt6s3hzMH86/q87tLD/OQWZBM9d3J8pp6Rz8aC29Tm0BseVhq1P72ouhDkI7LBgk1Xrbi60iLbPDNCuSVp8PWrY1ZGRVwQsAB4mIy6KjLjTgAnA+272+AA4hYy4MWTEBQG3oVLBtnj6Gm8F9gLgXMv7uagypatRfm1fIFHm74NADkprv01K+RWAEKKdEKJQCNEOQEr5A/C0ZRz7LK8HfTQWv2My2dK7Fm89xo5jqsrVpIFtmpXZ0xG9AGvqGjY4++0dNWyAs3ul8sjE3iRFN910Jle48mE34eJubnHUsLVKeM2JNgkRdr3CbSZxBw27Lj5s33Ajqi7IMeAjYDoZeUrDzohrR0ZcIRlx7dTnvG3AFcBslMybAJznyX8N3qd13aZ7/6wQYiVKu/6xZufj9vjHce8fQEq5HwdtXkr5HPCc6z0Cn+Fdk1nwzyHMFlu4ScDE/s3THK4xMD3B+r5NgvvSnE2F3g5+++bmm/SEq6pVD45vHnnHehwFdqeWzcscDqpuwaD2CdYKbyFB6po4TuojGzoPOyMvG5joZp2TDCMjbwFKGfaaWk1JpJRLa7OfgY3Tu7SwVvcClYOZ6hBE0dwYmJ7Ag+N7IiXWwghNmQ4tookMDaK4vBKTgLZuWmo2R/SGptHdW/LfcT1p36L5TWicBHYzvAYA/dvGWwX23hOqXoOjpSGqIfOw6wm3ZyiEWIoyVXtESjncpyNqJiRFh9E7LY4NllrSk5th7rUjQgiuOq1DQw+j3ggyCXq2imXVvhzaJkY2mxxjbyiyNHQAaBkT1iyFNeB0TzQHy5Mr9G2G2yS49uFr9cWbMp6mJG7reBv4hrN7pbDhUB7xkSGc1TO1oYdj0ABcPawD+7KLuWZY85moeMPZvVJ5+JvNVFRVcf0Id8kiTR9HDbtlbPOKZdDomRbLTaM6sfFQPlNPbe9ym+YQ/+NWYEsp363PgTRHrh/RibiIEE5ql+CUumDQPDi3TyvO7dOqoYcRcCTHhPHTzOFUVFY55ac3J8KC7AV2SmzzdZvddXZ3p2UvX3oSt3z0D2f3SmmAEdU/Ho3+luYbZVLKjZbPycALQG9gBXCnlNJ1Y+cAorKykmnTpjF+/HjGjx/f0MOxEhJkYsrQ9g09DAODgKS5msH1OGrYzVlgu2J8vzSGdEykhS+awmxbCNsWEhdGwGpPwrErjN1K5cd+SEq5yPL5KyANlfN8KbBeSnljPYyzTkRFRcmioqLqNzQwMDAIIBZuOMJ0XbeyPU+c2yxMvw2JEKJYShmQs8Xqwup6AEsBhBDxwFigt5RyuxDia2A5KvfMwMDAwMDHOGrYhrBu3lQXVheM6u8JqrtIppRyO4CU8gAQ77+hGRgYGDRvmmNfAQP3VHc3bAIutLy/BFikrRBCtMa5O4mBgYGBgY8I1QWdxbgo12rQvKjuDvgP8I0QYjZQCQzTrbsY1aXEwMDAwMAPhIXY4p+aegc7g+rxqGFLKZcB7YAzgY5Sym261d8BM/04NgMDn5GRkcEVV1zR0MNwydixY3n3XSOLsiZU9//s1asXv/32m1+++4YbbuCRRx7xy7E3b97MoEGDrJ/PHdqXkr1rAUjQ9Up3xRdffEHbtm2Jjo7mn3/+8cv4asrjjz/OtddeW+12F1xwAT/88EM9jKhxU62DREpZIKVcLaUscFi+TUp52H9DM3jiiSc499xz7ZZ16dLF5bKPP/7Y7XHeeecdhg0b5nZ9oNK+fXsWLVpktyzQzmXkyJG8+ab3NYZcCZqFCxcydepUXw+tQZg3bx6DBg0iOjqaVq1aMXbsWJYtW1bv49i0aRMjR46s83Fc3W+zZ8/mgQceqPOxXfHAAw9w5513Wj/rY8wGt0/0uO+dd97JK6+8QmFhIQkJCQghMJvNLrf9+eefSUlJISsry7qsrKyMHj168Prrr9ftJHTcd999Xv0+7rnnHmbNmuWz722qGBENAczw4cP5448/qKxUZRozMzOpqKhgzZo1dst27tzJ8OFGhdja4u6hZlAznnvuOW677Tbuu+8+jh49yv79+7nxxhv56quvGnpojYIjR46wePFiJk6caF0WEmRiTI+WjOqWzMwzu3rcf9++ffTq5V2DlDPPPJNx48YxY8YM67JHH32UVq1aMW3atFqNvy6cfPLJ5Ofns2rVqnr/7kaFlLLJvyIjI2VjpKysTEZERMhVq1ZJKaX85JNP5L///W85fPhwu2WdOnWSubm58uqrr5apqakyLS1Nzpo1S5rNZrl582YZFhYmTSaTjIqKknFxcQ14RjUjPT1d/vzzz3bL3n77bXnaaafJp59+Wl5wwQV2626++WY5Y8YMKaWUu3fvlsOHD5fR0dFyzJgx8qabbpKXX365lFLKPXv2SEC++eabsm3btvL000+XlZWV8pFHHpHt2rWTycnJcsqUKTI3N1dKKWVJSYm8/PLLZWJiooyLi5ODBg2SmZmZ8r777pMmk0mGhYXJqKgoedNNN0kppbz11ltlmzZtZExMjBwwYID8/fffpZRSLly4UIaEhMjg4GAZFRUl+/btK6WUcsSIEfKNN96wnsecOXNk9+7dZXR0tOzRo4dcvXq17y+uj8nNzZVRUVHy008/dbtNaWmpnDFjhmzVqpVs1aqVnDFjhiwtLZVSSrl48WLZunVr+dRTT8nk5GSZmpoqv/jiC/ndd9/JLl26yISEBPnYY49Zj/Xggw/KSZMmyYsuukhGR0fLk046Sa5du9a6Xn/vPPjgg/LCCy+UU6ZMkdHR0bJnz57y77//tm77xBNPyI4dO1qv94IFC6SU0u1vZ+rUqXLWrFnW/efMmSM7deokExIS5Pjx4+WhQ4es6wD52muvyc6dO8v4+Hh54403yqqqKpfX591335VnnHGG3TL9eVRWVlrHmpiYKC+88EJ54sQJWVpaKqOioiQgIyMjZceOHWXbtm0lIKOiomRUVJRcvny5y/9Z69at5bfffis3bNgg4+Pj5c6dO12ObfLkyTIlJUXGxsbK008/XW7cuFFKqZ5R/fr1ky+99JKUUkqz2SxPPfVU+dBDD1mvvfa7c/c70rj22mtlRkaGy++vT4AiGQByy9WrwQdQH6/GKrCllHLkyJHyueeek1JKedNNN8m5c+fK++67z27ZVVddJSdMmCCnTZsmCwsL5dGjR+XgwYPl7NmzpZQ2IdfY8CSwDx8+LCMjI2VOTo6UUsqKigqZnJxsnciccsopcubMmbK0tFQuWbJERkdHOwnsKVOmyMLCQllcXCznzp0rO3XqJHft2iULCgrk+eefL6+44goppZSzZ8+W48aNk0VFRdJsNstVq1bJvLw8KaWzsJVSyvfff19mZWXJiooK+cwzz8iUlBRZUlIipbR/gGnoj/Hpp5/KtLQ0uXLlSllVVSV37Ngh9+7d68Or6h8WLlwog4KCZEVFhdttHnjgATlkyBB59OhReezYMTl06FB5//33SymVwA4KCpIPPfSQLC8vl3PmzJEtWrSQl156qczPz5cbN26UYWFhcteuXVJKdR2Dg4Pl/PnzZXl5ufy///s/2b59e1leXi6ldBbYYWFh8rvvvpNms1nec889csiQIdZxffrpp/LQoUOysrJSfvzxxzIyMlIePnxYSun6t6MX2L/88otMSkqSq1evlqWlpfLmm2+Wp59+unVbQP7rX/+SOTk5ct++fbJFixZy4cKFLq/PnXfeKW+88Ua7ZfrzeP755+WQIUPkgQMHZGlpqZw2bZq85JJL7L5rx44dUkrbPe7p/yGllF9//bVs06aNHDx4sHz++efdbjd37lyZn59vnXT169fPuk4T9ps3b5aPPvqoHDJkiDSbzVJK+/vd0+9ISimfffZZef7553scb30QyAK7WeQJ+KQ06cJ7IHOD7waV2gfGPlntZiNGjOD3339n5syZLF26lBkzZpCWlsbrr79uXXb77bczbdo0cnNziYiIICoqipkzZzJnzhyuv/76Og/1tttuY+3atXU+jkb//v154YUXvNp24sSJBAfbbtPy8nIGDBhAq1atGD58OPPnz+e6667jhx9+oEWLFgwcOJD9+/fz999/s2jRIsLCwhg+fLjL/3tGRgZRUaqg0Ycffsjtt99Ox44dARU/0Lt3b95++21CQkI4ceIEO3fupG/fvgwcONDjmPU+6jvuuINHH32Ubdu20a9fv2rP98033+Tuu+9m8ODBAHTu3Ln6i+SCh77ZxObD+bXa1xU902I99qM+ceIELVq0sPtfOfLhhx/y8ssv07JlSwAefPBBrr/+emsAV0hICLNmzSIoKIhLLrmEadOmMWPGDGJiYujVqxe9evVi/fr11v/RwIEDmTx5MgC33347zz77LH/++Senn36603cPGzbMGvsxZcoUu/vvwgsvtL6/+OKLeeKJJ1i5ciUTJkyo9rp8+OGHXH311QwYMABQ901CQgJ79+6lffv2gPLPxsfHEx8fz6hRo1i7di3nnHOO07Fyc3NJSkpy+12vv/46r7zyCm3aqM5+GRkZtGvXjvfff9/jdffE+PHjee+999izZw+33nqr2+2uvvpq6/uMjAwSEhLIy8sjLi6O3r17c//993P++edz9OhRVq5cSVCQc3XP6n5HMTEx5Obm1uo8fEIjKE3aLAR2UFAQc+bMqdtBMjfAvvoPnhk+fDivvvoqOTk5HD9+nC5dupCSksLUqVPJyclh48aNdO/enYqKClq1sjWRqKqqom3btj4Zw9q1a1myZIlPjlVTvvzyS8aMGWP9/M4771iDWKZOncprr73GddddxwcffMCUKVMAOHz4MAkJCVZhDJCens6BAwfsjq2/PocPHyY9Pd1ue7PZzNGjR5kyZQoHDhzgkksuITc3lyuuuILHHnuMkBDXUbvPPvssb775JocPH0YIQX5+vl1wjycOHDhAp0517061+XA+f+3JrvNxvCUpKYmsrCzMZrNb4eHqGh8+bItbTUpKsj7oIyJUC8WUFFtTh4iICAoLba0L9P8/k8lEmzZt7I6nJzXV1g0vMjKS0tJS61jfe+89nnvuOfbu3QtAYWGh1/+vw4cPW4U1QHR0NElJSRw6dMgqsB2/W38OehISEigoKHC5DpSP+vzzz8dksoUeBQUFcfToUVq3bu3VeF3Rq1cvwsLCrMeNjo62rtu8eTOtW7dm1qxZzJ8/n+PHj1u3y8rKIi4uDlC/xVmzZjFp0iS6dOni8nuq+x0VFBQQHx9f6/OoM93GQrex5JW9XFn9xg1DsxDYPiG1T4Mcb+jQoeTl5TFnzhxOO+00AGJjY0lLS2POnDmkpaXRrl07wsLCyMrKcvmwrGs5w/79+9dpf38db+LEiUyfPp2NGzfy7bff8vTTTwPQqlUrcnJyKCoqsgrt/fv3O10H/ee0tDT27dtn/bx//36Cg4NJSUkhODiYBx98kAcffJC9e/dy7rnn0q1bN6655hqnYy5dupSnnnqKX375hV69emEymUhISFD+J6r/X7Rt25Zdu3bV/qJY6JkWW+dj1OR4Q4cOJTw8nC+//NKq9TqiXWMtMGr//v2kpaXVekz6CVhVVRUHDx6s8fH27dvHddddxy+//MLQoUMJCgqif//+Xv+/HO+boqIiTpw4USsB2rdvX4/pfW3btuWtt96yPgc8UZffvOOE4v333+err75i0aJFtG/fnry8PLt7GuDGG29k3Lhx/PjjjyxbtsxlJkdISIjb3xHAli1bvLJCNWcMge0tXpiv/UFERASDBg3iueees0t7GDZsGM899xxjxoyhVatWnHXWWdxxxx088sgjREdHs2fPHg4ePMiIESNISUnh4MGDlJeXExpa8+IL3pqv65vw8HAmT57MZZddxsknn0y7du0ApbkNGjSIBx98kMcff5yVK1fyzTffcN5557k91qWXXspTTz3F2LFjSU5O5r777uPiiy8mODiYxYsX06JFC3r27ElsbCwhISFWTTAlJYXdu3dbj1NQUEBwcDDJycmYzWaefPJJ8vNtpumUlBR+/vlnqqqq7DQljWuvvZbbb7+dYcOGMWDAAHbt2kVISIidZuoNnszX/iAuLo6HH36Ym266ieDgYM466yxCQkJYtGgRixcv5umnn+bSSy/l0UcfZfDgwQghePjhh+uUG7969WoWLFjAeeedx0svvURYWBinnHJKjY5RVFSEEILk5GQA3n77bTZu3GhdX91v57LLLuOSSy7hsssuo0ePHtx3330MGTLEql3XhDPPPJMZM2ZQWlpKeLhzV64bbriBWbNm8e6775Kens7x48dZvny5S9N9cnIyJpOJ3bt307Wr5+jy6igoKCAsLIykpCSKi4u577777Na///77rF69mnXr1vH1118zdepU1q1bZ6epAx5/RwBLlizhgw8+qNNYmzpGWlcjYMSIERw7dsxu1nr66adz7NgxazrXe++9R3l5OT179iQhIYHJkydz5MgRAEaPHk2vXr1ITU2lRYsWDXIO/mLq1Kls2LDBag7XmDdvHn/99ReJiYk89NBDXHnllR6Pc/XVVzNlyhSGDx9Ohw4dCA8P5+WXXwZU6tzkyZOJjY2lR48ejBgxwipoZsyYwWeffUZCQgK33norZ599NmPHjqVr166kp6cTHh5uZ7rV/KVJSUl2plT9+lmzZnHZZZcRExPDxIkTyc6uP9N2Xbj99tt57rnnePTRR0lOTqZt27a88sor1jSl+++/n0GDBtG3b1/69OnDgAEDuP/++2v9fRMmTOCTTz4hISGB999/nwULFrh1U7ijZ8+e3HHHHQwdOpSUlBQ2bNhgp8FW99s544wzeOSRR5g0aRKtWrVi165dHmsieCIlJYXRo0e7TYObMWMG5513HmeddRYxMTGccsop/PXXXy63jYyMZNasWZx22mnEx8fz559/1mpMAFdeeSXp6em0bt2anj172k2K9u/fz2233cZ7771HdHQ0l112GYMGDWLmTOeaWp5+R3///TdRUVGcfPLJtR5nc8Bje82mgtFes+myf/9+unfvTmZmJrGxvjUDGxjUN5s3b2bq1KmsXLmyWXXmmjRpEtdcc41TUaiGIJDbaxoC26DRUlVVxe23305+fj5vvfVWQw/HwMCgCRDIAtvwYRs0SoqKikhJSSE9Pd2oQWxgYNAsMAS2QaMkKirKbXqMgYGBQVPECDozMDAwMDBoBBgC28DAwMDAoBHQLAS2Vpr0m2++aeihGBgYGBgEItsWwte3BnRp0mYhsLXSpLWuI15H6lwWtQlhXAsbxrWwYVwLG8a1sFGv16LbWDjvJfLKCNjSpM1CYDviStN2p33XZLm7bd3ddN4ew59ja2zXor639ef3GdfCRmP7jfjiGMa1qH7bQPiNBBKGwPawrKbLa/rPDtQfYCBfi/re1p/f5w7jWtT8GE1JSLnDuBb+/75AF9jNonCKEKIKKNEtCgIns4erZTVd7m7bYMBch2P4c2yN7VrU97b+/D7jWthobL8RXxzDuBbVb9sQv5FQKWVAKrPNQmAbGBgYGBg0dgJyFmFgYGBgYGBgjyGwDQwMDAwMGgGGwDYwMDAwMGgEGALbS4QQYUKIuUKIfUKIAiHEP0KIsbr1ZwghtgohioUQi4UQ6bp1oyzL8oQQe10ce68QokQIUWh5/VRPp1Ur/HktLNvMEELsEUIUCSG2CCG61sNp1Rp/XQ8hRDvdPaG9pBDijno8vRrh599JfyHEUsv6g0KI/9bTadUKP1+LU4UQKy3HXS+EGFZPp1Ur6ngt7hJCbLTst0cIcZfDsdtb9im2HGNMfZ5bvSKlNF5evIAoIANoj5rojAMKLJ9bAHnAhUA48H/An7p9TwamANOAvS6OvRcY09DnGCDX4lpgPdATEEAnILGhz7mhrofD93RARba2b+hzbqB7YzPwGCqStxNwBDivoc+5vq8FkAhkWfYNAq4AcoCEhj5nP12Lu4EBqIjxbsA+4BLd+hXAc0AEMAnIBZIb+pz9ch0begCN+WURLJMsP6rluuVRqDSy7g7bj3HzINpLIxLY/roWlh/yAeCMhj6fQLgeLo75ILC4oc+toa4FUAz01H2eD9zb0OdX39fCIuw2OSzbDlzT0Ofnz2uhW/8S8LLlfVegDIjRrV8K3NDQ5+ePl2ESryVCiBTUzbIJ6AWs09ZJKYuAXZbl3vKhEOK4EOInIUQ/nw7Wz/jwWrSxvHoLIQ5YzF8PCSEa1X3qh3tD40rgXV+Msb7w8bV4AbhSCBEihOgGDAUW+XTAfsSH10JYXo7LevtmpP6nttdCCCGA0y37Ydlmt5SyQLfZOlf7NgUa1YMwUBBChAAfAu9KKbcC0SiTjp48IMbLQ16OMg2lA4uBH4UQ8T4ZrJ/x8bVoY/l7FtAHGAVcClzjm9H6Hz/cG9pxTwdSgM98Mc76wA/X4ltgMkr72grMlVL+7aPh+hUfX4vlQJoQ4lLL5GUqykUQ6csx+4s6XosMlNx62/LZJ7+vxoIhsGuIRdt7HygHbrYsLgRiHTaNRfloqkVK+YeUskRKWSylfALlgzndNyP2H364Flo1uqellLlSyr3A68C5dR+t//HHvaFjKvC5lLKwToOsJ3x9LYQQicAPwMMoP2db4GwhxI2+GrO/8PW1kFKeACYAtwNHgXNQloaDPhqy36jLtRBC3IyyMv1LSllWk32bCobArgEWc8xclKYzSUpZYVm1Cein2y4KNePd5HQQ75A4m7wCCj9di22oH3KjK7/nz3tDCBGBCshpFOZwP12LjkCllPI9KaVZSnkQ+JgAn8z5676QUi6RUg6WUiaigtO6ASt9OXZfU5drIYS4GrgHFd+in5hsAjoKIfQadT9q/+wNbBraid6YXsBs4E8g2mF5MsoMMwk1+38K+yhHk2X5WFSEYziqXi1AO+A0INSy/C7gOJDU0Odb39fCsv49lOkzBmUi30ojCKbx1/WwbHOZZZ1o6PNsqGuB0ppyLdfCBKSiooMfa+jzbYj7AjgJCLFclxeAPxr6XP14LS4HMoEebo77J/CMZd/zMaLEjRfKvyyBUpQZRntdblk/xiJcSoDf0KXeACMt++pfv1nW9UJFSxYBJ4BfgEENfb4NcS0s62NRmlMBKmL8vwS4oPLn9bBs8yPwSEOfZ0NfC2A08Dfq4Z4JvAFENvQ5N9C1+MhyHfKAT4CWDX2+frwWe4AKh/1m69a3t+xTgrLSNeqMG08vo/mHgYGBgYFBI8DwYRsYGBgYGDQCDIFtYGBgYGDQCDAEtoGBgYGBQSPAENgGBgYGBgaNAENgGxgYGBgYNAIMgW1gYGBgYNAIMAS2gYGBgYFBI8AQ2AYGBgYGBo0AQ2AbGBgYGBg0Av4flsY+Puzo/RoAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "fig, ax1 = plt.subplots(figsize=[7, 3])\n", "ax2 = ax1.twinx()\n", "h1, = ax2.plot(dates, tsW, '-', color='C1', lw=2.5, label='Wet')\n", "h2, = ax1.plot(dates, tsD, '-', color='k', lw=2.5, label='Hydrostatic')\n", "h3, = ax1.plot(dates, tsT, '-', color='C0', lw=2.5, label='Combination (left Y-axis)')\n", "\n", "# axis format\n", "ax1.tick_params(which='both', direction='in')\n", "pp.auto_adjust_xaxis_date(ax1, dates, every_year=1, buffer_year=None)\n", "ax2.spines['right'].set_position(('axes', 1.0))\n", "ax2.tick_params(which='both', axis='y', direction='in', colors='C1')\n", "for ax in [ax1, ax2]:\n", " ax.yaxis.set_major_locator(ticker.MultipleLocator(0.2))\n", " ax.yaxis.set_minor_locator(ticker.AutoMinorLocator())\n", "ax1.set_ylabel('Slant range delay [m]')\n", "# ylim\n", "y1min, y2max, yspan = -3.17, 0.03, 0.8\n", "ax1.set_ylim(y1min, y1min + yspan)\n", "ax2.set_ylim(y2max - yspan, y2max)\n", "fig.tight_layout()\n", "ax1.legend(handles=(h1,h2,h3), ncol=3, loc='upper left', bbox_to_anchor=(0.02, 0.18),\n", " handlelength=1.5, columnspacing=1.0, frameon=False)\n", "\n", "# output\n", "out_fig = os.path.join(work_dir, 'tropo_TS.pdf')\n", "print('save figure to file:', out_fig)\n", "fig.savefig(out_fig, bbox_inches='tight', transparent=True, dpi=300)\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "id": "cf935eaf", "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 }