diff --git a/applications/plot_dis_ts.ipynb b/applications/plot_dis_ts.ipynb index 7ee4584..175db96 100644 --- a/applications/plot_dis_ts.ipynb +++ b/applications/plot_dis_ts.ipynb @@ -16,7 +16,9 @@ "name": "stdout", "output_type": "stream", "text": [ - "input lat / lon: -1.0352 / -91.1917\n" + "read from file: /Users/yunjunz/data/archives/Galapagos/GalapagosSenDT128/mintpy/geo/geo_timeseries_ERA5_ramp_demErr.h5\n", + "input lat / lon: -1.0352 / -91.1917\n", + "corresponding y / x: 1803 / 906\n" ] } ], @@ -27,21 +29,22 @@ "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from mintpy.utils import utils as ut, plot as pp\n", - "from mintpy.defaults.plot import *\n", + "plt.rcParams.update({'font.size': 12})\n", "\n", "# work directory\n", - "proj_dir = os.path.expanduser('~/insarlab/Galapagos/GalapagosSenDT128/mintpy')\n", + "proj_dir = os.path.expanduser('~/data/archives/Galapagos/GalapagosSenDT128/mintpy')\n", "\n", "## file in radar coordinates\n", "#geom_file = os.path.join(proj_dir, './inputs/geometryRadar.h5')\n", - "#ts_file = os.path.join(proj_dir, 'timeseries_ECMWF_ramp_demErr.h5')\n", + "#ts_file = os.path.join(proj_dir, 'timeseries_ERA5_ramp_demErr.h5')\n", "\n", "# file in geo coordinates\n", "geom_file = None\n", - "ts_file = os.path.join(proj_dir, 'geo/geo_timeseries_ECMWF_ramp_demErr.h5')\n", + "ts_file = os.path.join(proj_dir, 'geo/geo_timeseries_ERA5_ramp_demErr.h5')\n", "\n", "# read dates and displacement\n", - "dates, dis = ut.read_timeseries_lalo(lat=-1.0352, lon=-91.1917, ts_file=ts_file, lookup_file=geom_file)" + "print(f'read from file: {ts_file}')\n", + "dates, dis, dis_std = ut.read_timeseries_lalo(lat=-1.0352, lon=-91.1917, ts_file=ts_file, lookup_file=geom_file)" ] }, { @@ -53,12 +56,12 @@ "name": "stdout", "output_type": "stream", "text": [ - "save to file: test.png\n" + "save to file: dis_ts.png\n" ] }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZ0AAADSCAYAAACRph2vAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydeVhVVdfAf+cCggMqCKg44ACKoOKUmpWZ85hDWimar6amVlam5piYQ5aZOad+pmKaiXNmTq+aYc44JpjgLDiAE4ggyvr+uNz7glzkAhcZ3L/n2Y/37Omsc7KzXHuvvZYmIigUCoVC8TzQ5bQACoVCoXhxUEpHoVAoFM8NpXQUCoVC8dxQSkehUCgUzw2ldBQKhULx3LDOaQFyO0WLFpUqVaqYbLt37x7FihWzWNuzxty6dQtnZ+fncq+80JbZ+Sz9HrNDxtzSpt6jZdos/R5zy3Ol13b06NH7IpK6UURUeUZxcnKStOjfv79F2541pm7dus/tXnmhLbPzWfo9Pqstt7yrzLap92iZNku/x9zyXOm1AbfExDfVys/Pz6SWUuhZtGiR35AhQ9Jsr1q1qkXb0qpfuHAhAwYMeC73yittmRmTHe/xWW255V1ltk29R8u0Wfo95pbnelbbhAkTwv38/GY/Xa+JOhz6TOrVqydHjhzJaTGoV68euUGOvI56j5ZBvUfLkJ/fo6ZpR0Wk3tP1ypEgj/Csfw0pzEe9R8ug3qNleBHfo7J00iG3WDoKhUKRl1CWjkKhUChyHKV0FAqFQvHcUEpHoVAoFM8NpXQUCoVC8dxQSkehUCgUzw2ldBQKhULx3FBKR6FQKBTPDaV0FAqFQvHcyHVKR9M0R03T1mua9kDTtEuapvV4Rt/PNE27rmnaPU3TftI0zTZZ2x5N0+I0TYtJKmefGttM07QQTdNiNU3brWmaW3Y+l0KhUChyodIB5gKPgJKALzBf0zTvpztpmtYKGAk0AyoAlYAJT3X7SESKJJWqycY6AeuAcYAjcAT41fKPolAoFIrk5Cqlo2laYeAtYJyIxIhIILAJ6GWie29gsYj8IyJ3gInAf8y8VRfgHxEJEJE4wA/w0TTNM6vPoFAoFIq0yVVKB6gCPBGRf5PVnQBSWTpJdSee6ldS07QSyeq+1jQtUtO0fZqmNUlrrIg8AMJM3efWrVvUq1fPWBYuXJjhh1IoFIoXgYULFxq/lYCTqT5pZg7VNG2vmfeJE5GWmZDPFEWAe0/V3QPszehr+G0PRAFfAGfQL9W9C/ymaVotEQlLGnvLnPs4Ozvn29DjCoVCYUkGDBhgjJytaVqkqT7PSlf9EjAwnXtowMxMSWeaGKDoU3VFgWgz+hp+RwOIyMFkbcs0TesOtAVmZ/A+CoVCobAQz1I6f4vIsvQmeJZ3WSb4F7DWNM1DRM4l1fkA/5jo+09S2+pk/W6ISFQacwt6JWkY29vQkLSXVDmN+ygUCoXCQqS5pyMizcyZwIJLa4a9lXXAV5qmFdY07RWgI7DcRHd/4H1N07w0TXMAxgJLATRNK65pWitN0+w0TbPWNM0XaAxsSxq7HqiuadpbmqbZAV8CJ0UkxFLPolAoFIrU5DZHAoDBQEHgJvALMEhE/tE0rXzSeZvyACKyFfgW2A1cSirjk+awASah37eJBD4GOonI2aSxt9B7yU0G7gAN0O/7KBQKhSIbMStzqKZpPsAMoBb6TXjQL1WJiBTIPvFyHpU5VKFQKDJOWplDn7Wnk5xfgLXAEOChJQVTKBQKxYuDuUqnFPClmGMWKRQKhUKRBubu6SwDLOmlplAoFIoXEHMtnanAfk3TRgM3kjeISFOLS6VQKBSKfIm5SmcNcAG9q7Ha01EoFIocIiYmBltbW2xsbHJalExhrtKpBZQQkUfZKYxCoVAonk2HDh3w8fHhhx9+yGlRMoW5ezp/AV7ZKYhCoVAons2ePXu4cOEC/v7+RERE5LQ4mcLcczpzgW7ol9ee3tP5MntEyx2oczoKhSK30KRJE/r06cOJEydITEzM1dZOVs/pFAJ+BwoA5ZLVKxdqhUKheA7s2bOHq1ev4uvrS6tWrfDy8uKLL76gdOnSOS1ahjDL0nmRUZaOQqHIDRisnN699bGKhw4dmqutnbQsnTT3dDRNK2jmxGb1UygUCkXm2LNnD6dPn6ZatWoEBQURFBRE8+bNWbRoUZ7b23nW8toNUuecMcU1wNEy4igUCoXiaS5dukT58uUZODBlirOqVaty/vz5PLXElubymqZpj4BVZszRWURMZfbMF6jlNYVCocg4mXEkmGzm3FMzJ5JCoVCkZNeuXQQEBDB//vycFiXPEhUVhU6nw8HBIadFMYlyJEgHZekoFM8HEaF+/fqcOnWKv/76i5deeimnRcqTtG7dGltbWzZu3JijcmTYkUChUCieJ1u2bCEuLo5p06YxYcIEi84dGxtr0flyK/v37yc4OJgjR44QFBSU0+KYRCkdhUKR44gIfn5+jB8/nv79+3P8+HEOHz5skbkTEhKoWbMmy5Yts8h8uZkJEyYwZswYRo4caXHFbSmU0lEoFDmOwcrp0qULdnZ2jBo1ymIfTX9/fwC++uorEhISLDJnbmHt2rXEx8cD/7Ny/vOf/9C/f/9ca+2YpXQ0TeuWRn1Xy4qjUCheNJJbOTqd/pP0/vvvW8TaSUhIYNKkSSxdupQKFSrw888/W0LkXMGhQ4fo2rUrP/30E/A/K6dAgQLY2dnlWmvH3Nhr90Uk1ZkdTdNui0i+PqOjHAkUiuzlxIkT1KpVi3LlyhmVDsCtW7fo1asXP/74Y6bnXrx4MatWrWLHjh3s3buXPn36EBISkmfTAiSnXbt2lCxZkh07dvDzzz/Tpk0bVqxYQYECBQCIi4ujZ8+e7Nu3jzp16jx3+TIVe03TtEpJP3WaplUEtGTNlYA4y4moUCheRGrWrMnly5d58uRJqjZnZ+dMz2uwcpYvXw5A48aNjdZOnz59Mj1vbuDQoUOcOnWKdevW0blzZ9atW8cbb7zB//3f/6Xo17RpU6Kjo3NIStOkF/AzFH1QTw0Ie6rtOuCXDTIpFIoXCE3TKFeuXKr6jz/+mNq1a9O3b99Mzevv74+DgwMuLi78+++/APTo0YNJkybRs2dPi1o7Dx48QKfTUbDg84kKNmHCBEaNGoWtrS3jx4+na9euhIaGYmtr+1zunxWeuacjIjoRsQL+SvqdvLiKyMLnJKdCocjHnD59OsV1SEgIy5cvZ9y4cTx8aH6y4qioKKZMmYKI8O+//xITE0P79u2NZcqUKURFRXHx4kUSExN58OCBWfM+evSIZ21F+Pr68v7775stZ1YwWDkGZdygQQNq1Khh3NvJ7ajDoemg9nQUiuzl4MGDNGzYkL179/Laa68B0LNnT6pVq8bhw4dp2rQpQ4YMMWuukSNH8s0337Br1y7eeOONVO1jx47l66+/ZsGCBVy+fJmdO3eyb98+NE0zMRs8efIEKysr2rdvT+3atZk4cWKqPseOHaN9+/YkJCTw559/Uq1atQw8fcZ58803iYqKomXLlsa606dPc+DAAS5cuIC1tbkZa7KXtPZ0zHUkqIg+LE4toEjyNhEpbykhk+7lCCwGWgKRwCgRWZlG38+AL4CCwFpgkIjEa5pmC8wDmqMPRhoKjBaRP5LGVQAuAMn/mfONiKT6G6WUjkKRvbRt25bY2Fisra3ZuXMnISEhvPbaa4SFhREWFkb79u0JDQ1Nd+nq1q1beHp68vnnn7Nt2zb27NmTQplERUVRpUoVZs+ezciRI4mNjaVo0aLMnz+fVq1apZpv6dKlzJs3jx9++IGuXbsSHx/P2bNncXBwwMrKCtB73nXu3JmmTZsSGxvLyZMnWbnS5OfKYvj7+xMWlnK3IzIykmPHjrF3795cr3QQkXQLsB/4GWgDvJ68mDM+IwX4BfgVvXJ7FbgHeJvo1wp9JGxvwAHYA0xNaiuMfr+pAvolxPZANFAhqb0C+r0q6/TkqVu3rigUiuzhwIEDUq5cOYmJiZGKFSvK3r17xdfXVyZNmmTs07FjR5k5c2a6c33xxRcyaNAgSUhIEA8PD9m1a1eK9jFjxki/fv1ERKRSpUry2muvyapVq6Rhw4aSmJiYom98fLxUqFBBXF1dpXbt2rJgwQIZMGCAfPrpp+Lq6ipHjx6V7du3i6enp7i6ukpsbKzcv39fnJ2d5cyZMxZ4Mym5c+eO/P3332m2d+vWTaysrFI9c04CHBFT33hTlak6wX1AZ07frJQkZfEIqJKsbrlBmTzVdyUwJdl1M+D6M+Y+CbwlSukoFLmGNm3ayPz580VEZPHixdKwYUNxcnKSe/fuGfsEBQUZP+xpERERIY6OjnL58mUREfH395fGjRsblUlkZKQ4OjrK+fPnJSoqSooVKyZlypSR2NhY8fLykq1bt6aYb9GiRdKiRQv55ptvxMbGRuLi4uTixYtiZ2cnpUqVkg4dOkiDBg2kQIEC0rNnT+O4r7/+Wrp3757h93Dt2rVnKpVPPvlEihUrJnfu3EnVdurUKSlZsqTMnz8/xTPnNFlVOpuBuub0zUoBagMPn6obBvxmou8J4J1k105JiqSEib4l0bt3e0pKpXMNuAosAZxMyVS+fHmpW7eusSxYsMBi/1EUihcZg5UTFxcnIiKPHj2S4sWLi06nkyJFiqQoQCrFYCA2NlaKFSsmHTp0MNY9be0kt3LGjh0r/fv3l5YtW8qiRYtSWTsRERFSoUIFCQwMlNatW0vZsmVl/fr1cvv2bbG1tZV+/fpJiRIlpGzZsmJnZydlypSRv/76S/bv3y/bt2+XggULZtjaefvtt6VUqVImFWt4eLg4ODhIu3btxM/PL1V7t27dZNq0aWlaeM+TBQsWGL+VwEUx8U01d09nDvAusA69q7QREfky3QnMRNO014AAESmVrK4/4CsiTZ7qGwZ8KCJbk65t0FtJFUXkYrJ+NsAfQJiIfJBUVwTwBI4DJYC5gL2IpFrYVXs6CkX20LZtW4oWLUrXrv8LbLJ+/XquXr3K5s2bU/W3tzedtuu7775jxIgR2Nvbp9hcP3XqFBUrVmTLli0UK1aMkiVLYm9vz8mTJ6lUqRLR0dEUKVKEgIAA2rZty5IlS3Bzc8Pb25v69eszY8YMOnTowFdffcWcOXNo3rw5586dIzAwkOLFi2Nra0uJEiU4c+YMLi4uODk5ARAdHc0333xD8+bN2bFjBx06dDAp9+PHjwkICKB69eo0b96c2rVr06ZNGz755JMU/T799FM0TePDDz+kYcOGhIaGUrx4cUDvQNC8eXPCwsIoXLgw/v7+LF68ONV+Vk6QqcOhySgM/AbYAKkd6i1HDKmzlRZFvx+TXl/Db2NfTdN06JfnHgEfGepFJAYwaJIbmqZ9BERomlZURO5n6QkUCoVZuLu7Ex4ezqpVKXNFent7p6lgnubhw4d8//33fP/990ybNg1PT09q1qwJQL169ShSpAiapnHw4EGio6MJDg5m0qRJJCYmUrhwYXQ6Hb169aJEiRKEhYWxZMkSrKysEBH27t2LnZ0dU6ZMISIiguDgYEqVKoWdnR06nY7Y2Fh+/fVX7t69y8CBAwkMDCQ8PBx3d3dCQ0MJCwvj008/5eTJk9SoUSOV7CtXrqR37968/vrrfP7557Ro0YJ27doxYMAAo9NEREQE/v7+/PPPP5QuXZr27dszc+ZMxo8fD+jjyQ0bNozChQsD+nNIEydOZM+ePSa993IDucplWtO0wsAd9I4D55Lq/IFwERn5VN+VwAURGZN03RRYabCSNL2a/wn9UlpbEUnT2V/TtJLoLbjiInIveZuydBSK3MusWbPYvXs369evZ+XKlcydO5fAwEBAr3Ru3brFuXPnzDo0GRISQqNGjbC2tiYqKoqiRYsaPcEePHhA+fLl+euvv2jTpg3vv/8+169f5+jRo6xZs4ZmzZoxcOBA/v77b65du2aMGdeiRQsePnxIQEBAins9fvyYatWqUbNmTX777TeioqKwt7enc+fONGnSxGjtGKycGTNmABAaGmq0dh4+fIirqyteXl7Y2dkZ575y5Qqvvvoq69aty/oLzgJZ8l5LUkzVgHHAnKTrqkBNc8dn4D6r0HuwFQZeIW3vtdboFYUXeu+1XSRzOAB+BA4ARUyMbZAkvw798tqvwG5T8ihHAoUidxIbGyuurq4SFBQkIiKPHz+WqlWryvbt22XTpk1Ss2ZNadWqldFRIT169OghU6ZMERGRGTNmSOfOnY1trVu3lhIlSoiDg4Og3w82lgIFCkiBAgWkY8eO4ujoKLdu3ZIaNWpIjRo1JCYmRkqWLCknT55Mca9ly5bJ66+/Lt26dRMnJyfZuXOniOidJkqXLi2xsbESEREhdnZ2MmLECJkxY4axeHh4yIQJEyQxMVFOnjwpR44cSVXCw8Mz9C4TExPl+PHjGRqTHmTRkaAbcCvpQ34/qa4esNOc8Rkp6M/VbEB/huYy0COpvjz6JbXyyfoORe82fR+9M4BtUr1b0l+IuKQxhuKb1N6d/53TiQD8gVKm5FFKR6HIncycOVOaN28uN2/eNJZZs2bJyy+/LHXq1JF169alclZIi+DgYHF2dpb79++LiF6hlS5dWo4dO2a2PAMHDpSRI0dKbGysODk5iZOTk8TGxsq0adOka9euxn4JCQni7u4uS5cuFU3TpEOHDuLm5iaTJk2SyZMnS+nSpWXWrFly8+ZNGTp0qHzyySepyvr16zP30tJg69atAsixY8ckISFBbt++neU501I65joSBAPdReS4pml3RMQhaYM+XEQyH5EvD6CW1xSK3Imvry/bt29PVZ+YmIirqysnTpxAp9PRtm1b3nzzTQYOHPjMuWxsbPjwww+Ndf/3f//HrVu3zFqmunz5MrVr1+bs2bOsWLHCuJHfuHFj+vfvT8WKFQkICOD111/H39+fn376iZ9//pn58+fz5MkTFi1aRKtWrahQoQIArVq14vXXX8/4S8kEIkKjRo0oWLAgxYsXx9PTkz/++IOjR4+miPqdUbIakSAKvUuxGNIZaJpmjV7puGRaqjyAUjoKRd5BRKhXrx5jx46lc+fOgD7MTrdu3Z65t/PWW29x+fLlVPUuLi78/vvv6d538ODBHD16lJ49ezJ27FiqVq2KTqfjn3/+YdiwYcyePRtN07h58yZVqlShWrVq1KpVyzh+7969AOzZsycTT501tm3bxtChQzl48CDu7u7ExcVRsmRJvvnmGzp16pTpebPqvXYU6IV+GcrAu8ChTEukUCgUWWDevHkULlyY3r17G+u2b9/OiRMn2Lp1awor6O7du6xYsSLNiNVr167Nkiy1atVCRDh16hRVq1bl8ePHAFSrVo1Lly7x8OFD7O3t2bVrFx999BH37qXwV6Jp06aULVs2SzJkBhF9Ar0vv/ySIkWKUKNGDUJDQ/n222/x8/OjY8eOFne9NtfS8QS2o98HaYg+5EwVoKUkeZnlV5Slo1DkPu7cuYO7uzs2NjbGMyoA165dY9OmTSbHNGvWjCpVqphs8/f356WXXsqWYJ09e/bEy8uL8uXLs2DBAvbu3ZvjZ2gMGKyckydPcufOHapWrYqVlRXbtm2jb9++jB8/PtPWTpaW15ImKIQ+hpkbcAXYLPrzLvkapXQUitzH+PHjuXLlCtHR0TRs2JDPP/8803NFRERQqVIlmjZtatZSWkYICQmhcePGhIaGUqhQIby9vZk3bx7NmjVLc8z69eupWbMmlStXtqgspnj55Zfx9PSkdevW/PLLLzx48MB4fql37974+fkRFBSUKSWZltIxe5dIRGJFZLWITBORVS+CwlEoFLmPO3fuMHfuXMaMGcOXX37JtGnTzM6LY4pvv/2WPn36cOrUKQ4dsuyOwaRJk/j000+NZ37GjRvHsGHD0pQ3MjKS9957j+HDh1tUjrSoWbMmsbGx/PLLL/z+++/ExMQQExND4cKFcXFx4datW2zcuNGi9zR3ea08MB59bLSnUxuYtlfzCcrSUShyFwYrx5C0rFu3bpm2diIiIvD29uaff/5hw4YNbN682WLWTkhICD4+PixfvpyiRfUBU+7evUv37t3p1asX/v7+qcaMHj2aq1evsnPnTrZs2ZLC2SA7OXHiBAMHDsSUPujTpw8ffPBBhufMqvfaQSAECABSnOwXkf9mWJo8hFI6CkXu4c6dO3h4eHDw4EHj8tOpU6do0aJFir0dc/nss88AmDFjBvHx8Xh4eLBmzRrq16+fZVkPHDiAn59firrQ0FDu3r3Lo0ePuHz5sjGGGuitnKpVqxIUFMSGDRv4888/LRZVICIigsmTJxu96J4HWc2nc4/nkNogNxZ1OFShyD3MnTtXChQoIF5eXimKtbW1+Pv7Z2guQ/Tm5Kf3582bJ23btrW02CKiz4lTokQJOXfunPTt21fGjx+fon3UqFEyYMAAEcnc4dRnMWTIENE0zeKHSp8FWTwc+jOwWER2W1gZ5nqUpaNQ5B7i4uJSZc004OHhQYECBcyea/To0SxZsiRFYMz4+HjWrVuXZpDOrODn58elS5dYsmQJ58+fp379+saI0cmtHDc3NwBmzpxpEWsnPDyc6tWrM3XqVObPn59px4CMktXlNUfgbyAMfdgZIyJi2vE9n6CUjkKRPzlz5gzHjx9PVa/T6Wjfvj1FihQxMSpz3L17F3d3dw4cOIC7uzsA77//PuXKlcPPzw8/Pz8WLVpEx44djWOio6P5+eefOXjwYJaW+z755BOsra357rvvqFu3Ll9++WWWDn2aS1YPhy4BngDBPLWno1AoFNnNo0ePMmTFmIOXlxdeXl4WnTMtZs6cSalSpdi7d68x+kDp0qWZMWMGn376KR06dMDFJWVwl127dqHT6diwYUOmlU54eDjLly/nzJkzaJqGn58f48ePz5ZDn+ZirqUTDbiKiKm8NvkaZekoFDlLZGQknp6ebNmyxSIb/DnBsmXL+PPPP1PVW1tbM3nyZJydU4awjI2NpXLlysyePZtBgwYRGhpKsWLFMnxfg5Uzffp0AMM+NePGjTOGCcousrq8tg99hOaL2SBbrkYpHYUiZxk9ejQrVqygevXqjBs3jqlTp7J+/fpcc6o/O5gxYwaBgYGsXbuW3r174+7uzrhx4zI0R2xsLI6OjlSpUiWFwjp16hSFChUiPDycMWPG8M477xgT31mSrCqdicA76JfZnt7T+clSQuZGlNJRKHIOwwb7gQMHaNq0KaVLl+aff/5h1apVaaaBzusYrJytW7fi4+PDuXPnaNSoUYatHREhKCiI2NhYY110dDQ9evRA0zSWLl3KW2+9RatWrSweiQGyHpHgVeAa0BJ94E9D6WkxCRXZwueff57lYIYKRU7x/fff8/bbb+Ph4UG3bt04ffo0S5Yswc/Pz+RBxvzAggULaNSoET4+PoDeK69t27bMmjUrw3MtXryY0qVL89prr/Haa69x8OBBunTpwtChQ/n444+ZOnUqJ0+etHgkhmdiyo9alfxxTic4OFgKFy4slStXlvj4eJkwYYI8ePAgp8VSKMzi1q1b4ujoKJcuXRIRkTfeeEMcHBxk//79UrNmTdm0aVMOS5g9VKpUScqWLSs+Pj7G4ubmJi4uLhmaZ+vWraLT6aRnz54iInL79m0pUaKEhIaGyp9//ik6nU5OnTolc+fOzZazSWTlnA6ApmklgLboM2xO0zTNFf2B0avZqBNznLy8vGaIbrtjxw68vLyYN28e06dPZ+jQoTktmkKRLqNHjyY0NJTp06dz8OBBhg4dSq9evTh+/Djvv/8+kydP5siRI/lub+fSpUvcuXMnVX3RokWpVKmSWXOI6BOz9enThzFjxhAYGMjKlSuN4YM6depEQkICTk5OLFy4EHd3d9auXWtRR42sRiR4HYgEtgLRyep+M2d8Xi551dIxpN+9d++e/Pe//5UCBQrIpEmTpFSpUsraUeQJOnToIGXLlpWyZcuKra2tODg4SNmyZcXT01Pi4+PztbWTVbZu3SpeXl7y+PFjmThxonTr1k0cHR3l7NmzcvjwYXF1dZVr165JiRIl5N9//80Wa4c0LB1zlc4xoFnS7ztJf9oBN8wZn5dLXlU6vr6+MnnyZBERCQgIEHt7e/npp5/krbfekunTp+ewdAqF+Vy5ckUAcXBwECcnJ2Oxs7OTrl275rR4uY7ExERp2LChrFq1SkRE7t69K/b29qLT6USn0wmQolSoUEFq164tgJw7d85icqSldMz1XrsjIg5Jvw3pqnXALREpkWn7Kw+QF5fXkufwKFKkCD4+PvTo0YPFixezevVq2rVrR1hYGIUKFcppURUKs7hz544xG2dy7O3tsbOzywGJci/JE7NZWVkB+hQLZ8+eZfny5YSHh3Pp0qVU46ytralbty46ndkZb55JViMSnNE0rZWIbEtW1xw4ZRHpFBbF39+f+/fvU61aNR4+fEhMTAyzZ88mIiKCu3fv8sorr/Djjz+qvR1FnsHBwSGnRcgzzJ49m8jISBo0aGCsi42N5ezZs8yYMQNXV1dcXV1zTD5zLZ2GwGbgd+BtwB/oAHQUkcPZKmEOkxctnfj4eCIjIwFo164djx49wsPDA4CCBQty+fJlwsLCuH79er7bhFUoXnSuXbtGREREqnpbW1uqV6+e46kNzLJ0ROSApmk10Z/L+Ql9uur6ks8913KKuLg4EhMTM738ZWtrS5kyZQD44YcfuH79eqo+9vb2SuEoFPmQMmXKGP//z42YpXQ0TbNFv3/zbbI6G03TbEUk3pICJUW0Xoz+IGokMEpEVqbR9zPgC6AgsBYYZJAnvXk0TWsGzAXKAweB/4hI6oXOHKBfv348fPjQIoc6mzRp8sz2OXPm4OLiwttvv53leykUCkV6mLtjtAOo+1RdXWCbib5ZZS7wCCgJ+ALzNU3zfrqTpmmtgJFAM6ACUAmYYM48mqY5AeuAcYAjcAT4NRueJcOEhISwbds29u/fz4kTJ7L1XpGRkYwbN47PP/+c+HiL/ttBoVAoTGK29xrgKMk6J3mvRRm82iwijKYVBu4A1UXk36S65cA1ERn5VN+VwEURGZ103QxYISKl0ptH07QB6C2bRsnuGxIopbMAACAASURBVAnUFpGQ5Pd53ns6vr6+eHt7U6hQIf76669sDWEzevRobt++zdWrV2nbti2DBw8G4MGDB+zcuTNFbg+FQqHICFk9HHoRfSSC5HWlgavmjDe3ALWBh0/VDcPEIVTgBPBOsmsn9H7nJdKbB5gJzH+q/TTw1tP3KV++vNStW9dYFixYkBFX9QwRHBwsTk5Ocu/ePWO62u+//17u3Llj8XsZQoxcvHhRDh06JGXLlpW4uDgREZk8ebLodDoJDg62+H2zQlBQkOzbty+nxVAoFGmwYMEC47cSvVGQ6cOh04FdQHWgEFAD/ZLb9+aMN7cArwHXn6rrD+wx0TcMaJ3s2iZJ6VRIbx70ez1Tn2rfh976ybHDoT169DAe6BQRGTt2rGiaJl988YXF7zVq1Cj54IMPjNft2rWTuXPnyr1798TZ2Vl69+4tvr6+Fr9vZklMTJRatWqJu7u7JCQk5LQ4CoUiHUjjcKi5ezpj0GcNPQREAweAs8BoM8ebSwxQ9Km6okn3TK+v4Xe0GfNk5D7PhZCQENavX0/Dhg05cuQIR44c4fjx41hZWTF//nxu3bplsXtFRkayYMECRo0aZawbP348kydP5ocffqBly5bMmjWL7du3ExISkmq8iJg8XJadbNq0CQBXV1dWrjTpV6JQKPICpjRRWgXQAGeS9oIsXYDC6Df/PZLV+fOUVZJUvxKYnOy6KUnWTXrzAAOAfU/dNxbwfPo+z8vS+eOPP1Is41WvXl2srKzEx8dHWrVqZVFr5/vvvxc7OzupUaOGsXh5eQkghQoVMi6rTZ482aS1s2XLFrG2tpbz589bTKZnYbBy1q9fL7t371bWjkKRByArsdfkfx9ne6Aiek+xSkCljIw38x6rgF+SFMErwD3A20S/1sB1wAtwQL/8N9WceZIU5z3gLfQx5L4BDpiSJ6dir/Xr10/GjRsnIiKXLl0SBwcHuXnzpkXmfvDggZw4cSJFGTRokBQrVkycnJwkMTFRRMS41BYcHCyJiYkyaNAgOXnypNSvX1/q168v/fr1s4g86bFhwwapVauWUa7XX39dli1b9lzurVAoMkeWlE7Sh/0YkAg8SfbnE3PGZ6Sgd2HeADwALgM9kurLo18WK5+s71D0mUzvo89qapvePMnamwMhwENgD1DBlDw5oXTOnz8vJUqUkKioKGPdwIEDs2VvR0Tkzp074ujoKMWKFZNixYpJz549Zfz48TJ+/HipU6eOvPfee7Jr1y6xtraWl19+Wby9vY2OCNlt7SS3cgwoaycl8+fPT7EXqFDkBtJSOua6TO8BgoCvgAvoN+u/Bv4WkZ/TnSAPkxNhcAYOHEhgYGAKl+Vr166xZs0arly5YvE4VBMmTCA0NBQvLy+OHj3KwYMH6du3rzFigZeXF3PnzqV79+4MGTKESZMmMWLECMaNG8f169dZtGiRReVJzpYtW+jUqRMffvhhiggKCxcuZNGiRXTv3j3b7p0XiI6OpnLlyiQkJBAcHEypUqVyWiSFAkjbZToj53RcRCRB07S7IlI86WzLaRGpmA3y5hpyQuls3rzZ5MFQGxsbPvzwQwoXLmyxe929exd3d3cOHDiAu7s7T548oUaNGsycOZMWLVoAsHv3bgYMGMCMGTPo378/DRo0YMOGDdy+fRsPDw+OHDlCxYrZ89fgwoULrF+/3mRb27Zt8fT0zJb75hUM6YadnZ2xtrZm+vTpOS2SQgFkXelEAJVFJFbTtFD0m/Z30B+2fNoLLF+RFwN+ZoQJEyZw4MABZs6caaxbs2YNv//+O4GBgYA+lE7fvn2ZN28e9erVY9WqVUycOBF3d3fmzJlDyZIls9XaUZjGYOX8+eefFCtWjOrVq3PmzBll7ShyBVlNbfAX+ujSS4E1wB9APPrNe0Ue5sKFC4SFhdG+ffsU9QULFuTBgwccPnyY8PBw3nnnHTZu3Mi5c+dwdnbGz8+PWrVqAVjU8lKYz9y5c2nevDnVqlUDoFevXkybNk1ZO4pcjVmWTooB+vA3vkARwF9EHmSHYLmF/G7ppEfXrl3Zs2dPiqi18fHxnD17litXrlC2bNkclC5tunXrhq+vL506dcppUbKF5FaOQemEh4cra0eRa8iqpWNERBKB5RaRSpHrmT9/PuHh4anqbWxscm349L1797J9+3b++ecfOnToYMyemJ/YvHkzUVFRNG3aNEV9dHQ0q1evZsiQITkkmULxbNK0dJICZKZrBonIe5YWKjfxols6eZFmzZrh6+vLokWL+OSTT3j33XdzWiSLk5iYyI0bN0y2GZwKFIqcJDOWTmg2yqNQZAt79+7l4sWL9OrVizJlyvDZZ5/RrVu3fGft6HQ6SpcundNiKBQZJk2lIyIT0mpTKHIrEyZMYOzYsdjY2NCyZUuKFStGQEBAvrR2FIq8iLkBP9E0rammaYs0Tfs96c9m2SmY4sVjz549BAUFmdX30aNHqeoMVk7Pnj0B0DQNPz8/vvrqK548eWJRWRUKReYwS+lomjYUfSyz28DvQBSwUtO0z7NRNsULxKNHj+jduzf9+/cnPY/KR48e4eXlRUBAQIr61atXc/HiRYoVK0ahQoUoVKgQnTp1Ijg4mFOnTmWb7I8fP862uZ/mhx9+ICoq6rndT6GwNOZaOp8DTUXkCxGZJ/osnk2T6hWKLLN06VKqVq3KkydP2Lx58zP7Llu2jMePH+Pn58fVq1fp1q0bjx49YubMmURHRxMZGWksUVFRxMbGGs8UWZrHjx9Ts2ZN1q1bly3zJ2f//v189tlnfPvtt9l+L0sSFRVF7969TVqnihcQUwHZni7ANcDuqbqC6CMSWDTgZ24rORVlOj+TmJgos2bNkocPH4qISHx8vJQvX17+/vtvWbdundSpU8cYUfpp4uPjxc3NTfbt2ycNGzaUNm3aiI2NjSxcuDBTsjx69Ej69Okjt27dytT4JUuWSJkyZcTb21uePHmSqTnMpVWrVjJu3DiLRhx/HowePVp0Op0sWrQop0VRPEfIYpTp/sAKwCNJ2VRBf1anH3prSQfozJkrrxWldCzPli1bBJDZs2eLiD7FbatWrURE5MmTJ+Lj4yObNm0y9k+ugBYuXCgtW7YUEZEVK1aITqeTNWvWiJubm8THx2dYlp9++kl0Op2MGjUqw2MTEhKkcuXKsmfPHqlfv76sXr06w3OYy99//218xkGDBsmIESOy7V6WJDIyUhwdHWXFihVSoUKFTP03UuRNsqp0EpOVJ2lcWzzNQW4oL7rSOX78uBw8eNBi8yUmJkr9+vXliy++kDJlysi9e/eMVo6B5NbO3r17pWrVqvLw4cMUVo6IyJAhQ6RUqVKyatUqadmypUlr586dO9K+fXu5d+9eqrZHjx5JpUqVZPny5eLo6JjC2jl79qxcunRJRET+/fdfk1bMkiVL5I033hARvSLNTmunVatWxue7fPlynrF2Ro8ebUyL3rx5c2XtvEBkVem4mVPMmSuvlRdZ6SQmJkqdOnWkatWq8vjxY4vMuWXLFqlevbo8efJE2rdvL76+vmJvby+DBw82loEDB4q1tbVs3bpVmjRpIi4uLjJnzhxZuHChNGjQQC5duiSHDh2SYsWKyXfffSdeXl6yd+9ek9bOhAkTpECBAjJp0qRUsvz000/StGlTERH54IMPjNaOIYdPixYtJDw8XOzs7GT58uUpxia3cgxjXnrpJfnyyy8t8p6Sk9zKMZCetfPkyRPZv3+/xWVJj927dxtzHxmsnIsXL4qIyF9//aWsnReILCmdVIP0S2wFMjM2r5X8qHSuXLkiixcvTrffpk2bpFatWtKoUSNZsWKF2fMnJibK77//nmpfxmDlBAQEiIjI4cOHxcXFRWbMmCFz5sxJVVavXi2VK1eW/fv3S5kyZWTgwIFSrlw5KVeunBQpUkSKFCki5cqVk0qVKsnVq1dTWTt37twRJycn2bx5szg5OaWwdgxWzp9//ikiIhcvXpTixYvLzJkzZf369VKrVi1xc3OTbt26ySuvvCJVqlRJkTRuyZIl8tJLL8n169eNZeDAgQLI4cOHzX5X5tCxY0fx9PSUnj17GkuLFi2kUKFCEh0dbXLM2rVrBZDjx4+nO/+BAweMWWqzwuPHj8XT01NKliwpDx48SGHlGFDWzotDVi2d74D6Sb/boc+2GQt0MGd8Xi75Uen07dtXrKys5Ny5c2n2MVg569evl+3bt2fI2tm8ebMAsn379hT1ya0cA+3btzfu7TxNkyZNZMmSJcZ+c+bMERH9x61gwYJSrlw58fT0NBYXFxepXbu2cfyECROkd+/eIiLSs2fPFNbO4sWLjVaOgWrVqgkg1apVk40bN8q3334r1tbWEh4eLo0bN05h7QwaNEhcXFyMxdnZWaysrMTOzk5atGhh1nsyl8DAQPH3909VVq9enSp76tq1a+XGjRtSs2ZNadOmjXTp0kXi4uLSdJRITEyU1157TWxsbOT06dMm+5w6dUo+/PDDdOVcuXKlvPLKK9K5c2f57rvvxN7eXtzd3aVevXrGUq5cOalRo4aIiOzbty9PLBEqMkdWlU4EUCjp90HgLfTpnk+ZMz4vl/ymdMLCwqREiRLyySefyH/+8580+23atEl8fHwkMTFREhMTzbZ2EhMTpV69evL2229Lo0aNUlg7r776qrz00ksycOBAY2nWrJm4urqmsoqmTJkilSpVMn5UDx8+LGXKlDF6vF26dEnOnDmTqly/fl1E/mflGBRrSEhICmunePHi4uDgIPXr15f69euLj4+P6HQ6sbKyEkdHR0lMTJSPPvpI7O3tJTAwUHbt2pXK2knOypUrpVGjRhIbGyuurq5y9OjRbPdme5qLFy+KjY2NtGvXTurWrSsPHjyQUqVKSY8ePaR27domPQL/+9//SpUqVWTSpEnyzjvvmJy3Y8eOYmVlZVxKNIXBytmxY4ccP35cSpUqJUFBQXLo0KFU5fz583L//n1xdHSUfv36Wez5FbmLrCqde0l/lgBuJau/b874vFzym9Lp27evfPnll3Lnzh0pUaKESWvHYOWsW7fOWJeWtZOYmJhijX7z5s1So0YNSUhIEE9PzxTWzsaNG2XevHmpytP7Jf/8848A0qVLF4mOjjaWFi1aGK0dU/z8888ybNgwEUlp5RgwWDtHjx6V4sWLi6Ojo+zevVsOHDggnTt3lvfee088PDyMisbBwUG+/fZbadGihSQmJqaydgwYPriGZ501a5Z4e3vLW2+9laas2cEHH3wgffr0ESsrK1m2bJmIiPj5+YmNjY14enrKxo0bJS4uzqh8DFbO8uXLJTo6WlxcXFJZO0FBQeLq6irz5s0zOk2YwqB0DXN36dJFpk+fnmb/KVOmSNu2baVEiRJy/vz5rD66IheSVaVzGH0OnfHAyqQ6J+CGOePzcslPSsdg5dy+fVtE9B8kU9bOli1bxMbGRoYNGyYjRoyQESNGyPDhw6VQoULy66+/pug7btw4ad26tYj8z8pZs2aNiKT+EJlLly5dxNbWVjRNk8KFC6coXbp0MTkmLi5OypYtK/b29nLu3DlxdXWVkiVLiru7u7GULFlSypUrJ2+++abMmjVLunbtKtOmTZOLFy+Ko6Oj+Pv7S5EiRaRt27ZSsmRJ8fDwkL59+4pOp5OpU6fKiBEjxMPDI5W18/Rz3rx5U3Q6nRQvXtzi+ztpYXiGJUuWiIuLiwwZMkRERAYPHiyFChWS77//XmrVqiX16tUzLmcarBzD80ydOjWVtdOxY0eZOXNmKseJ5DytdEXEaO08ePAgVf/79++Ls7OznDlzRsaOHausnXxKWkrH3HTVLwEzgUfA+yISpmmaL9BaRHqlO0EeJj+lNnj//fcpW7YsEyboY7nevXsXd3d3Dhw4gLu7u7FfaGgoa9euNTnHm2++aUwaFhUVRZUqVbCzs2P16tXcvXuXUaNGcfz4cXQ6HU+ePKF69erMmjWLFi1amJzv/PnzrF+/ns8/1we3OHPmDG+88QZhYWG89957vPbaa3z22WfpPtv8+fP57bffaNCgARcvXuTrr78mOjo6Vb8LFy7Qp08fwsLCOHfuHC1atKBdu3aUKlWKdu3asWLFCuLj41m6dCmlSpXCxsYGADc3N4oWLYqdnR1LliwxZksVEby8vKhbty4NGzYEYPv27Rw+fBhnZ2fc3Nz47bffjH137txJ8+bN0TQtlWwiwtKlS+nRowe2trbpPnNyBg4ciIODA1u2bMHT05PNmzczYsQIpk2bRv369SlevDjHjx8nJiYGa2trQkNDad26NQMGDDDGqouJicHNzY2NGzfy6quvcuzYMdq3b09oaCgFCxZk6dKl+Pv7s2tXyoTBGzdupGvXrnTo0CHFc23bto3p06fzwQcfGOvi4uKYOHEiFy5cYOXKldy+fRsPDw+OHDlCxYoVM/TMitxNWqkNMpw59EUjvyid8PBwypcvz5tvvknx4sWN9Xv27KFly5b8+OOPGZ5zzJgxREVFUbduXQICArh37x6lS5emcePGxj67du3i7t27BAYGmpzD19eXX375hZMnT1K9enW6d+9OrVq1+OKLLzh+/Dht2rQhLCyMQoUKpSlHfHw87u7urFmzhqpVq5pUpAY6duxI8+bN+fjjjwHo0KEDf/zxBx07dsTR0dHYb/fu3TRt2pSFCxc+8x0kJiYyevRoYmJiAH1cuJ9//plOnTrh4+PD7Nmz2bBhA/Xq1WPHjh20bNmSbdu20bJly1RzGe45b948Bg0a9Mz7JufSpUvUqVOHkJAQJk6cSHR0NIcOHSIiIoLSpUtTv359KlWqxJIlSwDw8fHBw8ODadOmUatWLWPuHRHh2LFjVK5cmX///ZdOnTrRtGlTY0K4x48f4+zsTLdu3VK8l6ioKPbs2WNStoYNG6ZI9jdo0CAWLVrEiRMn8Pb2BmDcuHFcv36dRYsWmf3MitxPhpWOpmmNRWRv0u+mJjsBIrIrrbb8QH5ROg8fPiQgIMBkcMrq1atTv379DM1nsHKCgoIoXbo0VapUoXXr1qmUQ0JCAoUKFeKbb75JNUdISAiNGzdm0KBBBAcH4+fnZ7RyihQpAkCXLl3StXYMVs6WLVsAfXqDixcvGj+yBoKCgmjXrh1nz57Fzs4OgMOHD9OmTRumTJmSyrqoXr06L7/8cobey9SpUzl58iQrV64EYM6cOWzbto1Nmzbx6quvUqpUKSIiIti3b18Kq0BEaNKkCT4+PmzYsIFz584909qJjY1l69atdOnShY8++oidO3fy5ptvGtsvXLjAmjVrWLBgAVWqVCEwMJBly5YRHR3NyJEj+eabb1izZg0FChQwjvnjjz8ICAjg+vXrrF69mqZNm+Lp6WmUIyEhgeDgYGxsbAgPD6dEiRIZejfh4eG4u7vz+PFj3n33XapXrw7A7du3mT59OqGhobi5uWVoTkXuJTNK57SIVE/6fSGNeUVEKllOzNxHflE6lmbMmDFERkayYMECABYtWkRAQADbt29P0W/48OH8+uuvhIaGpvjAgd7KqV69OkOGDKFy5crUq1ePyMhIOnbsaOwTHBzM9u3buXTpksmPcHIrp0GDBkDay4aTJ082Li0mR9M0tm3bRpMmTTL9PkCfKrpSpUr8/vvvxiXIuLg4atSowbhx45g7dy4nTpygZs2azJw5M4W1s3v3bj744APOnDlDx44dad++/TOtnW+++YZRo0Zx8uRJIiIiOHbsWIr28PBwNm/eTOnSpdHpdAQFBeHm5oaTkxN9+/Zl/fr1NGnShE8++QTQWzHe3t78+OOP7N+/n9OnTzNmzBgSEhKMc06dOhU7OzseP35MhQoVmDRpUobezyeffMKVK1cACAwMpHfv3kbFq9Pp+PTTTylVqlSG5lTkXtJSOjm+UW8ogCOwHngAXAJ6pNP/M+A6cA/4CbBNqrcFFifNEQ0cA9okG1cBfRrumGRlXFr3yU+OBJYiMjJSihYtKrt375bz58/L+fPnJSQkRJydnSUwMNDY78aNG+Lg4CANGjSQH3/8McUcwcHB4uzsLPfv3xcRkW+//Vbq168vI0eOTFUmTJiQ5in2P/74QwBxdnZOUWxsbMTPzy/7XoIJ9uzZYzy0mrwULlxY3Nzc5JdffhERvePByy+/nMKLrHHjxuLv7y8iIgcPHpRy5cpJXFycyfsYPM369Okjb7/9drpybdiwQSpXriznz5+XCxcuyIULF2TdunVSqlQpiY2NFRGR5cuXS+PGjSUxMVHu3bsnTk5OEhwcbJzj6tWr4uDgIDdu3JDz58+Lo6OjREZGmv1url27Jg4ODhIRESGJiYlSu3btFN6RivwHloxIkB0F+AX4FSgCvJqkTLzT6NsKuAF4Aw7AHmBqUlthwC9JueiA9knKp4KkVDrW5sillE5q9u7dKxUqVDBZ5s6da+w3bNgw+eijj+Tvv/+W8uXLp1AcPXr0kClTphivY2JixMXFRU6dOpUhWRITE+XGjRsmy6NHj7L+sBagdevWUq5cOaO7ucHba9u2bSIismvXrlRecW3btpV58+aZnM/gZRYTEyMlS5ZM95199dVX4ubmlqpUrFhRTp06JQkJCVKlShXZtWuXcczkyZNTeLJ99NFHRnd0EZF+/frJmDFjzH4HQ4YMkc8++8x4vXHjRvHx8XnuZ5kUz4+0lM6zlteuJH2cn4mIlE+vT3pomlYYuANUF5F/k+qWo0+dMNJE/5XARREZnXTdDFghIiZtc03TTgITRGStpmkVgAuAjYikm31LLa9ljps3b+Lp6cnJkycpW7YsrVu3pnPnznzwwQdcvXqV8uXL88Ybbxj3bkC/59KiRQt++umnHJTcshw6dIiXX36ZIkWKMHHiRONy0n//+19u3brFvn37aN68ObGxsTRq1Mg4LiQkhBMnTnDp0iUWL15M0aJFeeedd4iJiaFy5crs2rULb29vpk2bxpEjR/j1118zLePPP//M559/zuTJk43yXbt2jfHjx7N69WoaNWpEjRo1CAkJwcXFBdDvGdWrV49///033b2d8PBwKleuzP/93/8Zl89EhF69ejFv3jw6d+6cadkVuZfM7Om8nuzyJaA3MAv9spUb8BHgLyLTLSBcbeBvESmYrG4Y8LqIdDDR/wQwRUR+Tbp2Am4BTiIS9VTfkkky1xKRkGRKJxy9Ut0BDBeRSFOyubm5ibOzs/F6wIABDBgwIAtP+2IwfPhwHj58yJw5cwB9ArJ3332Xc+fOoWka27dvN+nU4OHhgZeX1/MWN9to164dOp2OgwcPUq9ePSpXrmxsK1myJGPHjmX16tVcvnw51dgiRYrwzjvv4O7uTsGCBQkNDWXmzJkcO3aMVatWAfDgwQMqV67Mzp07jRvzGWXr1q2psrAeOnSIixcvUrduXd59910GDRqUak/u0aNHrFy5kjfffJNNmzbRvXt3k/OfPHmSoUOHYupb85///IdevfL1qYsXioULFxo9G48ePXpJRCo83cfcczqngVYici1ZXVlgqyQ5G2QFTdNeAwKSWyqapvUHfEWkiYn+YcCHIrI16doG/RmiiiJyMVk/G+APIExEPkiqKwJ4AsfRR1iYC9iLSCtTsilLJ+PcvHmTChUq8PXXX1O2bFlj/fDhwxkxYgQDBw7MQemeH4cOHeKtt94iNDSUHTt2MHbsWIKCgtDpzE3Y+z9PvKioKJo0acKUKVNo2bJlCgeJLVu2UKVKFaPHXFYxWKlHjhyhSZMmBAQEpJl5tUCBAkanhqCgIGrXrm0RGRR5nyyd09E07Tb6D/q9ZHXFgQsi4mDG+D3A62k07wM+BvaJSKFkYz4HmjzD0pksIquTrksAkSSzdDRN0wErgaJARxFJeHqepH6l0MeWKyYi959uV0on41y4cIERI0aQmJiYqq1Tp0459i9bESEhISHVv9izi3bt2tGuXTsGDx6MiPDSSy8xZswYs5aTHj16xPHjx2nbti0HDhzg3r17vPnmmwwePJgnT56k6u/p6cnbb79tEbmTW6k//vgjv/32G7///rvJvtHR0VSuXJkuXbpw48YN1q9fb2z78ccfadu2LeXLZ3kFXpEHyZL3GrAU+BNoAVQDWgK7gWXmjDdj/sLoLRWPZHX+JDkHmOi/Er3SMVw3Ba4nu9aAJUkyFkzn3iXRL7MVM9WuHAnyD3PmzJGXXnopw2F5MsPBgwelbNmyKTzQfvvtN7M3z6dNmya2trbSo0cPY12HDh3SjMhtKQweh1euXBERfXihcuXKpZnI7+uvv5bu3btLbGyslC5dWoKCgkRE752o0+meGVRWkb8hi7HX7ICpQBj6tAZhSdfP/KBnpACr0HuwFQZe4dnea63Ru0t7ofde25VcQQE/AgeAIibGNgCqovdsK4HeY253WnIppZM/MER/LleuXIpU2NnF4MGDxdbWVkqWLGksLi4uAsjRo0efOTYmJkacnZ3F2tpaPv74Y2P9kSNHUkTazgo3btyQ6tWrS0RERIr6YcOGSbdu3SQsLMxYxowZI23btk01R/IYaiIiM2bMkE6dOomIiK+vrwwbNizNoLKK/E+WlM7zKOjP6WxAf07nMsnO6QDl0Z+nKZ+sbih6t+n7SVaN4ZyOW5LlEkfKszi+Se3d0TsSPEC/rOYPlEpLLqV08gczZ86Ujh07pkiFnZ3ExcVJREREqnLjxo10x06bNk1q164ttra2YmVlJT4+PlKnTh2pU6eOWFtby8qVK7Ms37Bhw6RYsWIp3JhFRDp16iQVK1ZMVWrUqJEqwrjByjFgsHbWrl1rTCORVlBZRf4nLaWjYq+lg9rTyfs8fPgQd3d3Nm/ejI+PD3Xq1GHixIl06JBquzDHMXijbdq0CSsrK/z8/HB1dU3hMent7W0M45MZDI4C27dvp2XLlpw5cybDkQDi4uIoW7YsTZs2pVKl/wUl2bFjB5GRkQwYMIAxY8akGR1Ckf9RAT8zdKMzLAAADeBJREFUiVI6eZ9Zs2axa9cuNmzYAMD69euZNGkSR44cMRnt+Xlw5coVjh07liJeGsB3333HoUOHWL16NaCP+N2wYUNCQ0NTBGrNCsOHDycuLo7Zs2fz6aefotPp+P777zM0R0JCAvPmzePhw4cp6m/evMnixYu5cuUKRYsWBdKOhafI3yilk0mU0snbPHz4kMqVKzN//nyjO29iYiLNmzdnxowZOWbt+Pr6snHjRi5cuIDhHNiDBw+oWLEigwcPNkZgBvj666/p2LEj48ePz/J9DVbOqVOnKFOmDBEREXh7e2fK2jFF//792blzJ6+88oqx7t69e/zxxx9cvnwZV1fXLN9DkTdQSieTKKWTtzl9+jTt27c36b79/vvvW+RDnlEM0bVbtWpFmTJlmDp1KgCRkZEMGTIkRZBNA6+//jofffRRlu+d3MoxkFlr58KFC7i5uaU4d3T06FGCg4NT9bW2tqZTp05ZWhZU5C0ypXQ0TasLxIvI6aRrZ+AHoDqwHxgmIjHZI3LuQCkdhaUxRNfu2bMntWrVIiQkhORRL7KTSpUqcenSpRSKIjExkTJlypiMipAWhgPACxYsUBEFFCbJrNL5C33Msp1J1xsBV/TndroDJ0VkcLZInEtQSkdhSQxWTlhYGPb29gwaNIhixYoZrZ3s5smTJyatPp1Oh5WVldnzDB8+nMDAQKKiojhz5owxEZxCYSCzSicSKCMi8UkRCG6SFJRT07Ry6OOllcs2qXMBSukoLInByhk1ahQAly9fplatWpw9e/a5WTtZJXkwV19fX/r166esHUUqMqt07gIOIiKaprUGFkqyqNKapkWLiH22SJxLUEpHYSnCwsLw8PCgb9++2Nv/73+btWvX0rt3byZOnJiD0plP8jA5yZPPKWtHkZy0lE56f0v+AboBq4F3gZ3JJiyDPmqAQqEwA0dHR2bMmJFqeeuzzz4zZj3N7Rhcok+ePAlAkyZNKF26NL/88ouydhRmkZ6l8yrwG/oT/k+AV0XkbFLbUKCBiLzzPATNKZSlo1D8j+HDh3Po0CFGjvxfmqvAwEACAgKUtaNIQaYsHREJ1DStPFAF+FdEopM1/44+XppCoXhBKFy4MIUKFWLWrFkp6r29vXn48GGKZUOFwhTmpjYoArwMOKFPIXDgKQWUb1GWjkKhUGSczO7poGnap8BE9JGmI9ErnjhN08aLSMZOkykUCoXiheaZKQw1TfsPMBJ4H7ATkdLolU8/YISmaX2yXUKFQqFQ5BvSs3Q+A3qLyDZDhYg8AX5Ncqf+Dn1aAYVCoVAo0iW9ZO2VSeYm/RT/hf9v795j5CrrMI5/H22h0FIQqcVQaQVbkUIoyRJUhJZLQPCGVETuUXARgjHGIHIrlwpE/AcTr5UiBSsglxqqkRJMG6EmBG80VkoNodBKCltg1162BcrPP96zcLrd2Vl2zzkzs/t8kpOdOe/7nvO+b2b2N+f2vhxQI83MzGwn9YLORmC/Gmn7ZelmZmYDUi/o/A74qaQdhoaVtBvwE2BRWRUzM7Php941ne+RTqOtkfRH0vTOHwROBtYB55RbPTMzG076PdKJiC7S8zlXke5aOyL7exXwyYjoLL2GZmY2bNR9Tici3gDmZ8vbJL1X0g0RMaesypmZ2fBS75pOf0aRjnjMzMwGZChBB0CF1MLMzEaEoQad+gO3mZmZZeoNg3NcrQU4tsiKSNpb0iJJmyU9L+msOvm/LWm9pC5Jt0vaNZe2TNJWSZuy5ZleZY+XtErSFklLJU0usi1mZta3ejcSzK+T/kJRFSE99/M6MBGYAfxB0lMRsbJ3RkknkW7nPg54kfS80PXZuh6XRsRtfZTdB3iQNH7cYtJgpvcCHy+wLWZm1od68+l8uIpKSBoLzAYOiYhNwOOSHgLOZcdA0uN8YH5PQJI0F1hYI29vpwErI+K+rOx1wAZJB0XEqiE3xszMahrqNZ2iTAO2R8Tq3LqngOk18k/P0vN5J0p6f27dzZI2SFouaVatshGxGXi21r46Ojpoa2t7e5k3b96AG2VmNpLMmzfv7f+VpGlwdtIsc8uOA7p6resCak1D2Dt/z+s9gFeAy4F/k07XfQVYLGlGRDyble0Y6L4mTJiAJ3EzM6uvvb2d9vZ2ACRt6CtPJUc62YX9qLE8DmwCxvcqNp7aA4r2zt/zeiNARDwRERsjYltELACWA6fUKFtvX2ZmVpBKgk5EzIoI1Vg+BawGRkmamit2GLDTTQSZlVl6Pu9LEfFKrSrwzjNFO5TNricd2M++zMysIE1xTSe7rvIgcIOksZKOAr4A3FWjyJ3ABZIOlvQ+4GrgDgBJe0k6SdIYSaMknQ0cA/RMRLcIOETS7Gz07DnACt9EYGZWvqYIOplLgN2Al4G7gYtzd6ftnz1vsz9ARDwM3AIsBZ7Plmuz7YwGvk+6brMB+CZwakQ8k5XtIN0pdyPwGnAk6bqPmZmVTBEeVKA/bW1t4RsJzMzeHUl/i4i23uub6UjHzMyGOQcdMzOrjINOi/BDqcVwPxbD/ViMkdiPDjp1dHX1fmb1HYsXLy40rb8y/X04i95XK6QNdntF92N/ac3SV4NNcz8Wk1Z0PzZLu+qlAXv2tdJBp47Oztozclf54eyPv8gDT+tPs9SxWdLcj8WkFd2PzdKuemnAXn2t9N1rdUjqAv5TI3lPdh6+Zyhp/ZXZh3QLeBX7aoW0wW6v6H4so47NkuZ+LCat6H5slnbVS5saETsd7TjomJlZZXx6zczMKuOgY2ZmlXHQMTOzyjjomJlZZRx0KiZpV0nzJT0vaaOkf0g6OZd+vKRVkrZIWippci7t2Gxdl6Q1fWx7jaTubHDUTZIeqahZlSuzH7M835L0nKTNkp6WNK2CZlWurH7MDdKbX0LSdypsXmVK/l7PkPRYlr5O0pyKmlUKB53qjQLWAjNJtxteA/xW0hRJ+5CmeLgG2Bv4K3Bvruxm4Hbgsn62/7mIGJctJ5bRgCZRWj9KuhC4APgMaabZz1L7ttZWV0o/RsQLuc/hOOBQ4C3ggTIb00Blfq9/A/w5KzsTuFjS58toRCUiwkuDF2AFabqFduAvufVjgW7goF75TwDW9LGdNcAJjW5PK/cj6YfYWuD4Rrenlfuxj21eCyxtdNtasR+BLcDBuff3AVc0un2DXXyk02CSJgLTSDOXTgee6kmLNLnds9n6gVooqUPSI5IOq599eCiwHydlyyGS1man2K6XNCK+KyV8HnucBywooo6toOB+vBU4T9JoSR8FPgE8WmyNqzMivkjNStJoYCGwINLMpePY+eneLmCPAW7ybGAKMJk0wd0SSX0ORTGcFNyPk7K/J5JOCR0LnEk63TaslfB57Nnu0cBE4P4i6tnsSujH3wNfIh0drQLmR8STBVW3cg46DZL9cr4LeB24NFu9CRjfK+t4YONAthkRyyOiOyK2RMTNQCdwdEFVbkol9GN39veWiOiMiDXAL4BThl7b5lXG5zHnfOCBiNg0pEq2gKL7UdLewMPADcAY4EPASZIuKarOVXPQaQBJAuaTfv3Njog3sqSVwGG5fGOBA7P1gxGAhlDVplZSPz5D+ocxYsaHKvPzKGk34HRGwKm1kvrxAGB7RNwZEW9GxDrgHlr4R5CDTmP8DPgY6U6z7tz6RaRrCbMljQHmACuyQ3QkvSdbPzq91RhJu2Rp+0s6StIu2frLSIMJLq+yYRUrvB8jYgvpzqLvStpD0iTg66RTHMNV4f2Y80XSEffS0lvReGX04+ps3VlZvn2BM8hdI2o5jb6TYaQtpOstAWwlHXb3LGdn6SeQztt2A8uAKbmys7Ky+WVZljaddLfMZuAV4E9AW6Pb22r9mKWPJ/2a3Ei6k20O2eC4w20psx+zPEuAuY1uZyv3I3Ac8CTpOtB64JfA7o1u82AXjzJtZmaV8ek1MzOrjIOOmZlVxkHHzMwq46BjZmaVcdAxM7PKOOiYmVllHHTMKiTpSkm3VbCfWZLeyuax+fQgtzEtK789m+7BbMhGNboCZsOJpPz4YrsD24Dt2fuLIuKmCqvzYkRMqp+tbxGxGhgnaVlxVbKRzkHHrECRJiwD0kyuwIUR0XLD0EsaFRFvNroeNvz49JpZhSRdJ+nX2espSlM4fzWbu+c1Sd+QdISkFZI6Jf24V/mvKU2f/ZqkJcpNe1xnv7tKelXSobl1H1Ca3nxCdjpunaTLJa0HflVow80yDjpmjXckMJU0kOOtwFWksbqmA1+WNBNA0qnAlcBpwATgMeDugewgIraRxpM7J7f6TODRiOjI3u9LmhJ5Mmm2S7PCOeiYNd7ciNgaEY+QBmy9OyJejoj/kgLL4Vm+i4CbI+Lp7NTXTcCMgR7tkKYXOCs3C+q5pLlferwFXBsR22LHUZLNCuOgY9Z4L+Ved/fxvuc60WTgR9lpt07gVdJ8SfsNZCcR8QQpqM2UdBDwEeChXJaOiNg6uCaYDYxvJDBrHWuBGyNi4RC2sYB0im09cH+vIOMh5610PtIxax0/B66QNB1A0p6STn+X27iLNLHaOcCdBdfPrC4HHbMWERGLgB8A90j6H/Av4OR3uY11wN9JRzWPFV5Jszo8iZvZMCTpGNKsnduAMyJiSS7tdtKDo1fX2cZU0oyVuwCXRMQd5dXYRgoHHbMRRNIU4J/A4RHxXGNrYyORT6+ZjRCS5pJOyf3QAccaxUc6ZmZWGR/pmJlZZRx0zMysMg46ZmZWGQcdMzOrjIOOmZlV5v/5tB3IJ/TTtQAAAABJRU5ErkJggg==\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaQAAADMCAYAAAA8nNe2AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABLEUlEQVR4nO3dd1zV1f/A8ddhuAAH4syVKJIjF+YsZ2pmZmqWZqY5Gt/SSi23KDlzlSvJleXKUZqaM81R+nPlVhBRzJGKIiiIcnn//rjcG8i6wL1cxnk+HuchfOb7fiLevD+f8zlHiQiapmmaZm8O9g5A0zRN00AnJE3TNC2L0AlJ0zRNyxJ0QtI0TdOyBJ2QNE3TtCxBJyRN0zQtS3CydwDpUbBgQfHy8kp2/b179yhUqFCa12V031u3blGsWDGrH9uWMWfHfW11nfW+CenrnDn75rbrfOTIkXARSXqliGS75uHhISnp169futZldN+6deva5Ni2jDk77mur66z3TUhf58zZN7ddZ+CWJPO7PVvesitcuHCK61955ZV0rcvovqlJ77FtGXN23Dc12fEzZcV9U5MVY86O+6YmK8acwc8bltwKJdlwpAYfHx85fPiwvcNIxMfHh6wYV06jr3Pm0Nc5c+S266yUOiIiPkmty5YVUlbVv39/e4eQK+jrnDn0dc4c+jr/R1dImqZpWqbRFZKmaZqW5emEpGmapmUJOiFpmqZpWYJOSJqmaVqWkOxIDUqpEAuPESUiVawUj6ZpmpZLpTR0kAfwUir7K2C99cLRNE3TcquUEtIqEfkjtQMopVZbMZ74x60MnATWiEgPW5xD0zRNyzqSfYYkIr0tOYCI9LVeOAnMAQ7Z6Niapmk5xunTp2nbti0xMTH2DiVD0tSpQSlVUClVOn6zRVBKqTcxjne00xbH1zRNy0nGjh3L3r17Wb58ub1DyRCLEpJSqpVS6iJwF/gnXrti7YCUUgWBccCg5La5desWPj4+5ubv72/tMDRN07KFU6dO8ccff7B69Wr8/PyyZJXk7+9v/n2NsX9CkiwaOkgpdRnwA1YCUfHXiYghY6EmOtfXwDURmayU8gUqPfkMSQ8dpGk5j8FgwNHR0d5hZDtdu3alXr16DBkyhObNm9O7d2969uxp77CSZY2hg/IBi0XkvogY4jfrhQlKqVpAK2CGNY+raVrWdvbsWUqXLs2NGzfsHUq2YqqOPvzwQwDGjBmTZaskS1iakGYAnyullC2DAZoBFYAQpdQNYDDQWSl11Mbn1TTNjvz8/HBycuKrr76ydyjZyrhx4xg8eDAuLi4ANGvWjDJlymTbZ0mW3rKrDGzFeO/vdvx1IlLRasEoVQAoGG/RYIwJ6gMRuWVaqG/ZaVrOcfbsWZo1a8aePXto2LAhZ86coWTJkvYOK8uLiIigSJEiuLm54ezsbF7+4MED6tevz++//27H6JKX0i27lN5Dim8NsBdYzRPPkKxJRCKBSNP3Sqn7wMP4yUjTtJzFz8+Pzz77jCpVqtCjRw+++uorpk2bZu+wsjwXFxd27dpFlSqJB8pxdXW1Q0QZZ2mFFA4UFpFY24eUOl0haVrOYKqOgoKCcHV15erVq9SoUUNXSRb46aef6NatG+fOnaNy5cr2Dsdi1ujUsB5oYb2QNE3TjNVR7969iYqK4tatW+TJk4cOHTroZ0mpiI2NZezYsTRr1owvv/zS3uFYjaW37PICG5RSe4F/468Qkazbv1DTtCztyJEjbN++nYULFyZYXrVqVTtFlD2sWbMGNzc31q5dS6VKlQgMDMxWVVJyLE1Ip+Oapmma1Zw/f97eIWQ7pupo2rRpFC5cmAEDBvDll1/y/fff2zu0DLMoIYnIWFsHomla7vLw4UNGjBjBhAkTyJs3r73DyTZM1VGbNm0AGDBgQLJVkohw584dihYtao9Q08zSoYOGKqXqPbHsOaXU57YJS9O0nO67775j+vTpLFmyxN6hZBum6qhOnTqsWrWKlStXsmXLFurUqZPks6QFCxZQrVo1IiMjkzha1mNpL7vrGIfweRBvmSsQICI2GWA1JbqXnaZlbw8fPsTT05Phw4czZcoUAgMDyZMnj73DyvIeP37M+++/z4MHDxKtq1q1KqNHjzZ//+jRI7y8vHBzc+Pdd9/l008/zcxQk2WN95DyAI+fWPYI45BCmqZpafLdd99Rt25d/ve//7Fx40YWL17Me++9Z++wsjxnZ+dEHUCS8/3331OlShUmT57MSy+9xHvvvUeBAgVsHGHGWNrt+wjw4RPL3gf0kD6apqXJw4cPmTRpEmPGjAGM469NmDCBR48e2TmynOPRo0eMHz+eMWPGUKtWLRo2bMj8+fPtHVaqLE1In2Icy+6IUuqnuLHlvgAG2C40TdNyIlN1VLduXQAaNGhA1apVWbx4sZ0jyzlM1VGjRo0AGD16NFOmTMnyz5IseoYE5mdG7YGyGOdB2igi920YW7L0MyRNy74qVKhAoUKFKFOmjHlZSEgIDx8+JDAw0I6R5QymZ0effvopDRo0MC//6KOP6N69u92fJVnjGRJxyWel1aLSNC1XWr58OXfv3k20PLt0Tc7qQkNDzSN+xx/128HBIcnrnpUkm5CUUt+LyDupHUAptVhEels3LE3TciKDwcD27dv54osvyJdP94myhVKlSrFv3z57h5EuKVVIXZRSS4DU5kB6DdAJSdO0VK1cuRJfX188PDz43//+Z+9wtCwm2WdISqlLgCUPmKJFxNuaQaVGP0PStOzHYDBQtWpV3n33XWbNmsWFCxeybJV07949nJyczBPf2VJoaKhVb1dm9dEZ0jXat4hUEJGnLWiZmow0TcueVq5cSfHixfn888+pXbu2xe/T2MMbb7xB//79bX6ea9euUaZMGXbs2GG1Y65YsQIvLy/Cw8OtdszMYmm3b03TtHQzGAyMGzcOX19flFKMGTOGiRMn8vDhQ3uHlshff/3F2bNn2b59O+fOnbPpuSZPnkylSpXw9fXF0h7PKTFd59KlSzNr1iwrRJi5dELSNM3mTNVRixbGadV8fHxsXiVdv36dP/74I837jR07lhEjRvDpp5/i5+dng8iMrl27xg8//MCWLVu4desWO3fuzPAxV61ahYeHB6tXr2bmzJnZrkqy+D2krEQ/Q9K07KVmzZo4Ojri5eVlXnbp0iWuXbtGSEiITc7Zo0cPtmzZQnBwMHv27OHIkSOMHj2a9evXExwczCeffJJon7/++os333yTwMBAoqOj8fT0ZM+ePXh7W//JxMCBA3FycmLatGksW7aMefPmsXfvXpRKrR9Z0gwGA9WqVWP27Nm0atWKHj168MwzzzBixAgrR54xKT1D0glJ0zSb+/3337l582ai5YULF6Zt27ZWP9+5c+d4/vnnadiwIQ0aNGDp0qVcuXKFv//+m7Zt23L79m0CAgIoUaJEgv3atm1Lp06dzM+PJk6cyKlTp1i2bJlV47t27RrVq1c3T9Vu6vAxZ84cWrVqla5jLl++nLlz55qTmukaBAUFUbBgQavGnxEZTkhKqTsi4p7E8psiUtwKMaaJTkiapqXEVB106tSJ+vXrU61aNZo3b87+/ftxcnKiWrVq5MmTh6lTp5r3+euvv2jdujXr1q0zjzweERFBx44dGTlyJL6+vlaLb+DAgVy+fJnhw4ebl61du5b9+/enq0oyVUe9e/emefPm5uWDBg2ibdu2WapKskZCihARtyeWOQM3RCTT+xbqhKRpWnLiVwYuLi4UKVKErl278uWXX/LUU0+xfPlymjRpQo0aNTh79qy5StqwYQPTpk1LcKyoqCiOHTuGm5sbV65csVo38L59+3LixIlEy11cXNi2bRvOzs5pOt6dO3d49dVXiY6OTrSuWbNmTJkyJd2xWlu6E5JSai/Gd5EaAn89sboMcFpEXrFioHmBuUArwB24AAwXkd/ib6cTkqZpJvfv3+fUqVPmcdviPztZvnw506ZNIyQkhHHjxjF58mRat26Nv78/AwYMSFQlPalv37489dRTnD17lnr16jFkyJDM+lg5VkYS0jsYR2qYh3G6CRMB/gV+F5En50nKSKAuwBBgCRACtANWADVE5JJpO52QNE0zGT58OLNnzyY4OJi7d+9SpUoVPv30U1xdXZkzZw7t2rVj586dREVFsWTJEnr16sWRI0dwdnZOVCXFd/HiRZ577jkCAwO5evUqLVu2ZOHChXh7e1OpUiU7fNKcwRq37LxFxLYd8pM/9wlgrIisNS3TCUnTcjfTaAQiQpUqVWjYsCG1atXi448/Zt68ecTGxhIcHMyPP/6Ih4cH0dHRREREkC9fPqKjoylevDilSpUiICCAUaNGMXToUDZs2ICI8OqrrwL/VUdjx44FoFOnTmzfvp3nn3+ezZs38+uvv1KtWjUqVqxoUcx///031apVS/PtuJzGKr3slFKtgVqAa/zlIjI6yR2sQClVArgM1IqfEHVC0rTcbfHixQwdOpSePXsSERHB0KFDqVu3LgEBAeYhcwwGA+fPn0dEePz4MRcuXDDvX6pUKfPsqZUqVSJPnjxUrlyZ2NhYLly4wLVr18zVUZEiRQAYMmQI33zzDcWLF2fRokV06dKFFi1a8PPPPycb50cffcT7779P0aJFqVixIlOnTk12DL+YmBg++eQTxo8fT6FChax1qbIca1RIs4GuwC4g/gxPIiLvWiXKxOd0Bn4DgkQkwdzG5cuXl2LFipm/79+/f6YM86Fpmv09fvwYb29vnJ2duXLlCmfOnKF8+fL079+f4sWL8+WXX6b5mPPmzWPjxo04ODjQpk0boqKiGDp0qDlpiQiRkZGICO+//z47duygXr167N69m02bNrFr1y5atmxJzZo1zcfcvXs3rVq14pVXXqF8+fKcPn2as2fPJjuG35IlS3j33XcZN24cI0eOTDLO69ev06dPH9auXUv+/PnT/Dntxd/fH39/fwCOHDlyWUQqJLmhiKTagFCgrCXbWqNhHEFiJbAZcH5yfd26dUXTtNxp0aJF0rx5c+nTp4/kz59f7t+/LyIiwcHBUqRIEbl06VKajvfw4UMpW7asHDx4UA4dOiRPPfWUPHjwQMLDw81t/Pjx0rFjRwkPD5dr166Jg4ODrF27VmbMmCEvvviiODk5SevWrRMct1mzZjJv3jwpXry4FCxYUK5duybt27eX2bNnJ4rh8ePH4unpKd999514eHhIWFhYkrEOHDhQ8ubNK19//XWaPmNWAhyW5H73J7dCEiaIAMDNkm0z2jB2oliMsRrLn9Q2OiFpWu4UFRUlFStWlA0bNoirq6t4eXlJx44dZdasWTJr1iwpUaKEeHh4SExMjMXHnDt3rrRr1878ffv27WXWrFnm7+/fvy8lSpSQEydOiIjI6NGjpVGjRvLSSy9JZGSk5M+fX959910pV66c7N+/X0REdu3aJZ6envL48WN54YUXxNPTU0TEnPBu3LiRIIbFixdLs2bNRESkZ8+e4ufnlyjOa9euSZEiRWTz5s1SunRpiYyMtPgzZiXWSEjvARswdv+uGL9Zsn9aGvAtcABwTW4bnZA0LffZsWOHFCtWTJo2bSq3b9+WgQMHyptvvin58+eX/v37y3vvvScFChSQkiVLyvLlyy06ZvzqyMSUNKKiokREZObMmVKsWDEZPHiwfPTRR5I/f37p2bOnAPLLL7+Ii4uLtGvXTubPn2+ukpo1ayZLliwxJ5HixYvL0aNHRUSkdevWki9fPtm0aZOI/Fcd7dq1S0REAgICkqySBg4cKJ988omIiHTs2DFDVVJYWJgcO3Ys3ftnhDUSUmwyzWDJ/pY2oDzGLuUPgfvx2lvxt9MJSdNyl9jYWGnQoIE4OjpKx44d5eDBg+ZWv359mTp1qsyYMUNee+012bZtm3h7e1tUJf3www+SN29eeemllxK0fPnyyZIlS0TEmKCmTJkiU6ZMkTZt2oibm5s0btxYGjduLJUrV5Y6deoIIMeOHZNy5crJrFmzzNXRwIEDpUGDBvLWW29J3bp15YcffpD69euLk5OT+Pj4SGxsbILqyOTJKsmU2K5duyYiIkePHs1QlTRgwAApXry4+XZnas6dOyfHjx9P17melFJC0mPZaZqW5W3fvp0PP/yQAgUKcPbsWZ599lkcHP6brKBZs2b8+OOPbN68mZo1a9KkSRM++ugjunXrluJxQ0NDOXjwYJLr6tevn2iSu8OHD/Pnn38CcPv2baZPn87IkSNxcXGhe/furF27ltGjRxMTE0PZsmW5desWkZGRGAwGwsPDadeuHTt27KBDhw6cOXOGr776imHDhhEZGUnZsmXN57l16xZXr14lNDQUpRSffPIJSilmzJhh3ua1116jefPmDBgwIE3X0jSOXu3atWnXrh2DBg1KcXsRoXHjxkRERHD8+PEE1z09rDa4qlKqLPCUiBzIUEQZpBOSpuUeIpIgwXTo0IEXX3yRjz/+2LzNzJkz2bNnD+vWrQOMCWzAgAGcOnUKR0dHm8Q1cOBAVqxYQe3atc3LHj16xO7du/npp58SvTzr7OzMihUruHnzJt999x1r1qxhypQpLFiwgNu3byc6fsGCBfHx8cFgMODq6krBggUT9KyLiIigdOnSnDx5Ms1xOzk58c4779C6dWvzEEtJERG2b9/OwIEDcXNz4/PPP6dLly5pOt+TUkpIlt5KKwfsBx4A9+OWdQEWWLK/tZu+ZadpuceTt+AOHz4spUuXNj/jiYyMlJIlS8qPP/4op06dklOnTsnJkyelSpUqqT5LunLligQEBKQrruDgYNmyZUui1q1bNzl8+HCi7W/fvi3u7u5y8eJFERExGAxSvXp187Ok5Kxbt04AWbp0qQQHBydooaGhaYr56tWrUqRIEbl+/bqIiHTu3FmmTp2a5LaTJk2SPn36SMOGDWX58uWyadMmqV69uhgMhjSd80lY4RnSb8BwjN2x78YtKwRctmR/azedkDQtd4iNjZVGjRrJ7Nmz5datW+bWokUL+eabb0TE+HyjRo0aUrVq1UQtqd5q8bVp00Zq164tsbGxVon36NGj4ujomKDXnsmIESOkTZs2cuzYMXMbO3as1KtXL9nzGwwGqVmzpvTo0UMaNWqU4TgHDBggn332mfn748ePS4kSJRI9S7p7964ULVpU3NzcpGLFihITEyOxsbFSr149Wb16dYZisEZCCgUc4r6+E295mCX7W7vphKRpuUNoaKiULl1aihYtmqh169YtQ8f+888/pVy5clKzZk1Zv369VeJ99dVXZfLkyYl67omIvPfee1KzZs1ErUGDBokSws2bN6VatWqyYMECqVu3rjx+/Fi8vb1l27Zt6Y7t2rVrkjdvXhk7dqzMnz/f3CpWrJioShozZoz06tVLypUrJ40bNzYvt0aVlFJCsnSkhjNARxEJMM2NpJSqCqwUkWfTcRsxQ/QzJE3TMso0GV/x4sUZN24cR44cSfdsrQDHjh3j5ZdfJigoiCVLlrBx40Y2bdqUrmN98cUXzJ8/H2dnZxYvXkz79u1ZsWIFs2fPZt++femK8+rVq/j5+REbG5toXevWrc3PhsLCwqhUqRIzZsygT58+ODs7069fP4oVK4aIMGnSJJYuXUqnTp3S9dms8QzpXYwvx/YGwoFuwEme6I6dWU1XSJqmZYSpOoqOjpbY2FipVatWhqukV199VWbOnCkiSb/fZKmbN2+Ku7u7TJkyRRwdHeWff/4REZGYmJgMV0mWGDNmjPTu3VvOnDkjw4cPl6ZNm0q1atVk+PDh5mZ6pyo9sEa3b6VUR6A/xneFrgDfisgv6UqRGaQrJE3LuUJDQxk7diwzZ87McBfj5Dw5Vfkvv/ySoSopfnVk6glnGh8vrVXSF198wf379/nzzz8JDw+nRIkSvPKKcdq5P/74gwcPHrB37940x2iJsLAwPD092bdvH56engCEh4fj7e3Nvn378Pb2zvA5UqqQnCw9SFzy+SXD0WiapqVg2rRpzJ49myZNmtC1a1erH//gwYPs3LmTdu3a8e233wLGO0UXLlzg119/pUOHDmk+5jfffMPjx49p2bKleVl0dDRHjx4lMDCQypUrW3ScW7dusWDBAv7++29WrVrF5cuXWbBgAXXr1sXV1ZXatWtTqlSpNMdnqf379xMREZFgkFiTbdu2WSUhpSQtFdLzQG0STz8xwQZxpUhXSJqWM92+fZsqVaowceJEvv76a06ePGn1KunEiRPMnTs3yXU9evSgSZMmaT7mtWvXuHTpUqLlTk5O+Pj4WPwZTNXRnDlzzMsGDRpETEwMX3/9dZrjssS1a9cIDg6mcePGNjn+k6wx/cQsjNNP7AWi4q0SEelplSjTQCckTcuZhg8fTmhoKN9++y2NGjXi008/tUmVlFVVrlyZoKCgBAksNjaW0qVL888//9jknN27d2fnzp1cvHgx2RdkrckaCekOUF1Erlk7uPTIiQlp27ZtXLp0Sc/rpOVapuro6NGjlC9fnq1bt/Lhhx9SpkwZtm/fTp48eewdos3FxsYm2QtOKZXhESd27drFggULWLZsmXnZ2bNnadq0KT4+PrRs2TLVYYSswRoJ6TjQQkRCrR1ceuS0hGQwGKhWrRrXr18nKCgIDw8Pe4ekaZnOVB3Nnz8fMD7X8fDwIDIykq+//lr/sZYBIkKDBg04efIkmzdvplmzZoCxOqpRowbt27fnxRdfTHEYIWuxRkLywThSwwrg3/jrRGSPNYJMi5yWkJYvX87cuXOpVq0aRYsWZcKETH8sp2l2FRMTQ5EiRahYsaJ5QNP79+9z9OhRKleuTGRkJIGBgbmiSrKF3377jcGDBzN48GCWLl3Krl27zNVRUFAQbm5uvP766zRo0MDmVZI1EtJ7wEyMY9k9+QypnDWCTIuclJBM1dHs2bOpXLkyderU4fz587pK0nIVEeHQoUPcv3/fvMzPzw8vLy8GDBjAZ599RufOnXWVlA6m6mjQoEF06tQJb29vFixYgL+/PzVq1GDYsGEAnDx5MlOqJGt0+54AvCIiO6wXlgawatUqPDw8aNmyJUopXn/9daZPn66rJC1XUUrx3HPPmb8/deoUZ8+eZePGjbi4uDBmzBi6detGr169dJWURlu2bOH+/ft06dIFBwcHRo0axRdffMHhw4cpX758gt81Tk5O+Pv78+mnn9olVksrpBCgkog8sn1IqcspFVL86qhVq1YAXL58WVdJWq73xhtv4OPjw5AhQ8zL2rRpo6ukNIpfHZl6K8bExFCpUiWaNGlCuXKJb3C9+OKLNG/e3GYxWeOWXS/gOWAccDP+OhFJ3CXExnJKQtq8eTPt27encePGCd4OP3LkCB999BGTJ0+2Y3SaZh+RkZEUKVIEZ2fnBD3LoqOjadCgAbt377ZfcNnMuXPneOaZZyhatGiCruTh4eF069aNxYsXJ9pHRDI0pl9qrJGQTEkn/sYK4zMk28x+lYKckpCioqI4dOhQouVDhw4lMjKSY8eO2fQHQ9OyqgcPHhATE5Noeb58+cibN68dIsq+QkNDMRgMiZYXKlQo0bWcPn06Bw4c4KeffrJZPNZ4hvS0FePR4uTPn58XXnghwbLDhw8TEhJCkSJF2LRpE+3bt7dTdJpmP5nxgmZu8eQ07MmJiIhg0qRJiAhHjx6lTp06No4ssbROYe4AlBCR67YLKXU5pUJKyiuvvELbtm0pWbIkEydO5NChQ7pK0jTN5iZOnMjJkydp2LAhO3bsYP369TY5T0oVkkUDLCmlCiullgMPgQtxyzoopb60Xpja4cOHOXbsGH369OG1117j0aNH6Z5PRdM0zVIRERHMmDGDUaNG0a9fPw4fPszRo0czPQ5LRy38FriHceoJU0+7v4A3rB2QUspdKfWzUuqBUuqyUqq7tc+RUQ8fPrTJXw9jx45l2LBh5MuXDwcHB8aMGYOvry+pVbHXr18nMjLS6vFompY7mHr6PvPMM+TLl4+hQ4cyduzYTI/D0oTUEhgQd6tOAETkFlDcBjHNwZj0SgBvAfOUUtVscJ50mz9/Ph07duT48eNWO+aRI0fYunUrjo6OLFmyhCVLlnDv3j0CAgJSrJJiY2Np27YtgwcPtlosmqblHhEREUyfPp1PP/2U+/fvc//+fbp168bBgwczvUqytFPDPcADMD87UkqVi/+9NSilXIDOGAdyvQ/sU0ptAN4GhlrzXOkVFRXF5MmTeeeddxg3bhxr1661ynEdHBzo3r07Bw4cSLC8U6dOuLq6JrOXcWKx2NhYVq1axbBhwyhbtqxV4tE0LXc4cuQIUVFRSb57tGPHjkzt3GBpt++hQAdgBPAz8BLG0RvWi8hMqwWjVG3gTxHJH2/ZYKCpiLxiWmbPTg1ff/01u3fvZtmyZXh6erJly5YkJ7PKDLGxsdSuXZsvv/yS/fv3Ex4enuw8L5qmaVlBhjs1AJOBnzDeTnMGFgHrAWvPGOWKsRqL7x7gFn/BrVu38PHxMTd/f38rh5E0U3U0evRoChQowJAhQxg3blymnDspv/zyC05OTrRv355BgwaxatUqrly5Yrd4NE3TkuLv72/+fY3xbluS0tTt29biKqT9IlIg3rJBQLOsUCGZqqOff/4ZML5Rbq8qKX519MorxkszdOjQLFsliQihoaF6OCRNy+Ws0e17qFKq3hPLnlNKfW6NAOMJAJyUUvEnoK8JnLbyedLs0aNHTJo0ifz58+Pn54efnx/Tpk2jfPnyfPll5vd+/+WXX7hx4wYODg5s3ryZzZs388wzz7Bo0SJCQkIyPZ7UrFixAi8vL+7de7IA1jRNM7K0U8NAYNYTy84AvwBTrBWMiDxQSq0Dximl+gK1gFeBRtY6R3oppRgwYACRkZE8evTfGLMtW7ZMsdOBrTx+/Ji6desyZ86cBMtbtGhBWFhYkoMm2ovBYGDcuHGUKVOGWbNmMXLkSHuHpGlaFmRpp4ZQoFT80b6VUnmAGyLibtWAlHLH+IzqRSAUGCoiy+Nvk5VGavD392fgwIFcvHiRUqVK2TucZO3Zs4dnn32WwoULZ/q5TRMQLlq0iMaNG3PhwgUKFSqU6XFommZ/1ujUcAT48Ill7wNW76QuIndEpKOIuIhIuSeTUVby6NEjxo8fT8OGDZkyxWqFokUiIiISVGopuX37Ni+//DJ+fn42jioxU3Xk6+uLl5cXL730ErNmPVlsa5qmWZ6QPgU+V0odUUr9pJQ6CnwBDLBdaFnfkiVLeOaZZ1i2bBnff/89169n3hB/r732WqqTaA0cOJDdu3czffp0XnzxRRYvXsyNGzcyKUKj+BMQAowcOZKvv/5aP0vSNC0Ri3vZKaVcgfZAWeAKsDHu5dVMlxVu2T169IjKlSuzcuVKGjZsaE4OM2bMsPm59+zZw9tvv839+/f5+++/k3wZ9tixYzRs2JBnnnmGkJAQjh07xtSpU3F2dmbatGk2jxH+m4DwnXfeSfDS3eeff07r1q31syRNy4UyPB9SVpMVEpK/vz/r1q1jy5YtgHE8uWrVqnH69GmbP0tq2bIlb731FgEBAURERCTq2ADQsWNHmjZtyoQJE6hXrx6bN2/m6tWr1KhRgzNnzlCyZEmbxghw584dXnvtNaKjoxOte+GFFzL9NqemafZnjQn6nDA+Q2qK8aUm83wIIvJCcvvZir0Tkqk66tq1q+lFLwBmzZpFvXr1bFol7dmzh169enH+/HnCwsLw9vZOVCUdO3aMl19+mYMHD1K1alUqVqzI33//be4pmJlVkqZpWnzWSEizgBaAPzAe4xBCHwArRcTXeqFaxt4JKSIigv/9739J/uXv4+PDkCFDbHZuU3X07rvvAsaXYZ+skjp27Ejz5s25desWa9asISwsjJdffhlvb2+uXbvG/PnzuXz5MsWKFbNZnJktODgYpRQVKlSwdyiapqXAGgnpKtBQREKUUmEiUlgp5Q3MF5GmVo43VfZOSPbi7+/PkCFD2LJlC87OzoBxGKXOnTtz/vx5ypYta66OgoKC2Lp1K/v37+f8+fMcOnSIt956C6UUTk5ODBkyBHd3q/bYtxsRoUmTJjg5ObF79+4cM6Hh2bNnWbZsmV1evNY0W7HGFOYFMHZkAIhSShUQkXNxQ/1omcBgMODr64vBYODjjz9OsK5q1apcunSJsmXL8t133xEREcGzzz5rXh8bG8uNGzfo3r27XaYltrUdO3Zw584dDAYDu3fvTnLU4uxo5MiR/PLLL3Ts2DHBrWFNy6ksrZD+BD4Rkf9TSv0KnAXCgbdE5Bkbx5hIbqyQli9fzpw5cwgLC2PatGm0bds2ye0iIiKS7Nrt4OBAxYoVc0z1YGKqjj766CNiYmJYsGBBjqiSjh8/Ttu2bfnss8/Ys2cPv/76q71D0jSrSKlCQkRSbUA9oE7c15WBHcBB4HlL9rd2q1u3ruQmMTExUqVKFdm+fbusXLlS6tevL7GxsfYOK0vYtm2beHt7S0xMjDx+/FgqV64sv//+u73DyrBOnTrJ9OnTJSoqSp566ik5dOiQvUPSNKsADksyv9stejFWRA6JyNG4rwNFpJWI1BeRvRnPl1pq4r9c2qVLFyIiIti6dau9w7KpXbt2pTqyhIjg6+vL6NGjcXR0xMnJiVGjRlk07XtWdvz4cf7880/ee+898uXLx7Bhw+wynbSmZbZkb9kppVpYcgAR+d2qEVkgN92yM71caprzHowJasaMGfz111/Z/tZUUkSEOnXqcO7cOU6dOoWnp2eS2+3YsYMOHTrw1Vdf4ejoCEBMTAyDBw/mt99+y7bPkjp37kyTJk3ML1s/fPiQSpUq8csvv+hnSVq2l95ODQstOLYAFdMVlWaRtWvXcvXqVXbt2sXu3bsBY5I6evQoO3bs4MUXX7RvgDawYcMGwDiiw/jx41m0aFGS2xUsWJCePXty8uTJBMt79eqFm5tbkvukJiIiAldXV7sl+suXL7Nu3TpCQkLM826BMdHOnDmTH3/80S5xaVpm0CM1ZHEnT55k/fr1Sa574403qFy5cpLrsitTdeTr60vTpk2pVKkSBw8eTLZKsqZHjx5RpUoVRowYQd++fW1+vqQYDAYOHjxITExMonUVKlTIUtOKaFp6WGXoIKWUI9AAKA1cBQ6KiMFqUaZBbkpIuc369evx9fXl6NGjKKUYM2YMV65cSbZKsiZ/f3+mTZtGdHQ0AQEB5MmTx+bn1LTcxhozxj4LBAKrgSHAGiBQKVXLWkFqmqmTgq+vr/mW2SeffMKGDRsICgqy6blNU4ksWbKEKlWq8P3339v0fJqmJWbp9BOLgDnAUyLyHPAUMBvLnjNpWrIePnzI2bNnAfjzzz/5+++/GTNmDLVr16Z27dq0aNGCR48eMX/+fIuPOW/ePFavXp2mOExTiTRs2JAxY8Ywfvx4i+eb0jTNOix9MTYcKBL/Fl3cLby7IlLQhvElSd+yyz4eP35sHuYoKWPHjmXWrFkEBwdToEABTp48mWSX7aefftqi2W7v3LlDpUqVcHFx4cKFC+TNmzfVfZ6cSgSgTZs2dOnShX79+qW6v6ZplrPGjLGbgQ5PLHsF2JSRwLScbdeuXZQvX56IiIgk14eFhTF79myqVavGnDlzcHR0pFatWubqKH6zdOr1GTNm0KlTJ2rWrMnChZYV8EuWLMHd3R1XV1dOnjzJyZMn6dy5M4MHD+bUqVOWftx0u3Xrlp6KQ9PA4pEaVgPRwJ/Aqrh/o4GfgKWmZsmxrNFy20gNabFz504JCQmxdxgSGxsrzz//vJQtW1YmTpyY5Da+vr7Su3dvOX36tBQrVkzCw8MzdM7Q0FBxd3eXoKAg+b//+z8pU6aMPHz4MNX9RowYIdWrV0/Qnn76aQGkQ4cOGYrJEoMGDRJA9u7da/Nz2dr27dstuuZa7kUKIzVYmpDGWNIsOZY1mk5ISQsPDxd3d3fp2rWrvUORnTt3SuXKleXkyZNJJpu7d++Kh4eHXLhwQURE3nzzzWQTl6VGjhwpffr0MX//8ssvy5w5c9J1rJYtW8rXX38txYoVk9OnT2corpTcuHFDihQpImPHjpVWrVrZ7DyZISAgQBwcHGTmzJn2DkXLwjKckLJa0wkpaRMmTJCOHTtKiRIl5OTJk3aLw1QdLV26VEREunXrlijZmKojk/RUSYGBgebPGb86MklLlRTfnj17pGLFivLo0SOZOHGivPnmm2naPy0GDRokH3/8sTx69EgqVKhgUZW0bt06eeutt2wWU3r17NlTXn/9dSlVqpRERkbaOxwti7JGhdQceDru65LA9xh73pW0ZH9rN52QEgsPD5dixYrJ2bNnZcqUKZlaJUVFRcnUqVPNA76aqqPHjx+LiMiZM2cSJJu7d+9KoUKFxM/PT5YvX25uVatWtbhKMiW9GjVqiMFgkJkzZ0q+fPmkdu3aCZqjo6OsXLkyTZ+nZcuWsmjRIhH577raokoyVUf//POPiIgsWLAg1SrJNNCum5ub/PXXX1aPKb0CAgLEw8NDwsLCpGPHjrpK0pJljYR0FigX9/XyuLYQ2GDJ/hYcP2/c8S4DEcAx4KXkttcJKbEJEyZI9+7dRUTk/v37mVolffPNNwLI+vXrzYli8uTJcuXKFXNr27atOdlcv35d3nrrLXnzzTcTtQULFlh0zp07d4qXl5fUrVtX1qxZIxEREXLkyJEkW1oqpPjVkYmtqiRTdWRiSZW0bNkyadKkicydO1fatm1r9ZhERDZv3my+lWqpnj17yrhx40RE5OjRo1KyZElp1apVmo+j5XzWSEjhcf86AaGAK5AHuG3J/hYc3wXwBSpg7PnXPi4xVUhqe52QEopfHZkkVyUZDAZp3ry51aZoiIqKktKlS8vIkSOldu3acu/ePalUqZI89dRTiZrpl/q9e/dk+/bt6T6nKen98MMP8uuvv5qrJGvo3LmzlCxZUho1amRutWrVEqWUXLt2zSrnEDH+N8ufP780b95cunbtam5eXl7JdqQwVUc7duyQhw8fStmyZS2qkv755x+5fPmyRXGFhYVJkSJF5JVXXrH4s8Svjkxq1aolzs7O0rt3b4mNjZXFixdLdHS0xcfUcq6UEpKlM8aGK6VKANWBMyJyXymVB0j+BZM0EJEHcQnJZKNSKhioC1yyxjlysm+//ZY8efIkGHjzwYMHrF69mrFjx+Lt7W1e/vPPP3P06FFGjhzJvn37MjyI6HfffYePjw9jx47l119/Zffu3QQGBqa4z+TJk5k2bRpBQUE89dRTaT7nrl27+Pfff3nzzTdxdHTE19eXn3/+mc6dO6f3Y5hNnz6dK1euJFqeJ08eSpYsmeHjmxQoUIAVK1bw8OHDBMs7depEhQoVktxn1apVFCtWjBYtWqCUYvjw4YwdO5bNmzen+N+xV69eREVFsXfv3lT/e8+aNYtWrVqxf/9+jhw5Qt26dVP9LBMmTMDT09PczT4mJobAwEBcXV3ZsGEDy5Yto3fv3kRHR/Pee++lejwt97L0xdgvgP9hrIo+EZGVSqnmwCQRqW/1oIzJ7zJQS0TOPbm+fPnyUqxYMfP3/fv3p3///tYOI9vYt2+feSTw+JRS9O3blxIlSgDGqcxr1arF+PHj+fzzz/nmm28yNFr4w4cP8fT05Ndff6VOnTr8/PPP+Pn5ceTIkWR/8YWGhuLl5UWrVq0oUaIE33zzTarniYiI4O7du5QrVw4RoWnTpvTv358ePXoAsHHjRoYPH87ff/+Ng4Olr9ZlL6ZpSObMmUPLli0BiI6OplSpUtSqVYvff096Fph9+/bx9ttvkzdv3gRTmCTl3r17VKpUif3797N161a2b99uHnk9JYsWLUrwvtaJEye4cOEC3bt3RynFwoULef3119mwYQOBgYF6jMBcyN/fH39/fwCOHDlyWUQqJLlhcqXTkw3wAjyf+L6Gpfun4TzOGGeknZ/cNvqWXfqsWbNG6tWrJ7GxsbJ8+XJp1KhRhmaenTFjRoLbSwaDQWrWrCnr169Pdp/hw4dLv379Ej3QFxE5efJkkt20+/btK1WrVhWDwSC7du0SZ2dnGT9+vEyePFkmT54skyZNEhcXF1m3bl26P0tWd+jQIQHExcVFXF1dxdXVVVxcXASQfPnyycGDB5Pcr1WrVvLdd9/JsmXLpHHjxin+9/bz85O3335bRP67FXv48OE0xRkdHS3lypWT/fv3i4jI+vXrxcHBQc6fPy9t2rSRb7/9Nk3H03IeUrhllynTTyildgNNk1m9X0SaxG3ngLHDREHgVRF5nNQOeuigtDNVR5MmTaJdu3YYDAaqV6+eoEpKbZif+AICAvD29qZBgwZUrVrVvPz48eM4Ozvz559/JtrHVB0dOXKEChUqMHjwYB49emSuktq1a8eOHTsIDAykfPnyAAQHB1OvXj3Kli3LsGHDqF69OkuWLEkypjfeeMOiW0zZ1f3794n//+u8efP4448/aNeuHZs3b2bTpoQDp5iqo/Pnz+Po6Jhoosf44ldHXl5egPH2XXJV0sGDB6levTouLi4Jlvv7+7N27Vq2bt1qrmY9PDwoXLgw/fr1480339RVUi6X0tBBKfasi/f1FSAkqZbc/mltgAIWA7uA/CltqyuktItfHZnEr5LCw8OlVKlSsmvXrhSPExUVJSLGXlXu7u7StGlT+e677xK0TZs2JbmvqToyiV8lHThwQMqWLSuDBg2S/v37m7fp27evjBo1SjZv3myukrKT0NBQ6dKli9Xfy4mMjJTSpUubexGWLVtWpkyZkqDbvKk6MkmpSpo0aZKULFlSBg4caG7vvfeeAHL8+PEE2968eVNcXFxkzJgxiY7ToEEDcXNzk9KlS0vRokXF0dFR3N3dJU+ePBIVFaWrJC19FZJSqomI7Iv7OrnqBhH5Ix1JMqnzfQvUAlqJyP2UttUVUtp169aNVatWmaf6Bsx/bQcHB/Pjjz/y7bffUqFCBXbv3p3kM6Dg4GBq1qzJqlWr6NmzJwcOHOC5557j77//pmzZsimePywsjLJly/L2228nmGRu2bJltGzZkoCAADp06MDrr7+Ol5cXR48eJTY2lnr16hEQEECRIkVo0KABgwYNomvXrla6KrY3evRoJk6cyNSpUxk4cKDVjvvNN9+wc+dO8+SNs2fP5osvvsDBwYFz585x6dIlWrRowaxZs8zViMFgYMCAAaxfvz5RlXT48GH27duX6DxKKbp37078Z7ZDhw7l1KlTHDhwgAsXLiQYZzA8PNw8dmHXrl2Jjo7Gx8cHJycn8ufPz5kzZzh58iSXL1+226y8mn1ZZYI+W1JKlcfYmy4aiD9V5nsisuzJ7XVCSrvY2FgMhsTzKSqliIqKwtPTk99//51OnToxf/58mjdvnmjbfv36sXPnTgB69+7NqFGjGDp0KOHh4cydO9e83b179yhUqFCCfcPDw/nqq694/DjxXdhChQoxb948AgMDyZs3LyNGjOD27dvExsZSqlQpxo0bB8Bvv/3G4MGDOXnyZLbovHDnzh28vLzw9/fno48+IigoiPz582f4uKb/Xr/88gt16tQBYOHChQwYMIAuXbrg7u5Oz549mTNnTpL79+vXj8aNG6fr3Ldu3cLb25u///6b0aNHU758eXx9fZPcduXKlUn2WHRzc+O9997TCSmXSu8tu3GWtOT2t2XLTbfsYmJizCMe2Er8l2qXLl0qL7zwQqLbOhcvXhR3d3fZtWuXKKXk//7v/0TEePvG3d3dPKDrzp07xcXFRW7cuGHx+V966SWZN2+e+fvbt29LoUKFxNnZWfz9/WXNmjWyZs0aWb16tRQvXlxWrVpl8bH/97//ydy5cy3e3ppGjRplHlvPNHqBNQa+3bdvnzg5OYmjo6M4OjqKg4ODAOLg4CCfffZZos4iyUlPh5YvvvhCPvjgAxExDt1UtGhRuXv3bpqPo+VepPOW3eJ43+YDOgOHMHbHLgc8B6wVkW4ZzZhplZsqpP79+2MwGCyeSiGtIiIi8PT0ZM+ePXh7exMTE0PVqlUTVUn9+vWjZMmShISEcOHCBUqUKGH+C3zMmDE4OTkxZ84cmjZtyq1bt3j55ZeZOnVqquc/ffo01atXp02bNgnmLtq/fz9FihTh2WefTbTPW2+9RadOnVI9dkBAAPXr1ydv3rwEBQUlegBvS6bq6P/+7/+oWLEix44do1WrVty7d4+AgAAqVqxotXMtWbKEpUuXmrt+Dx48mOjoaGbNmpXifh988AFPP/00n3/+OZGRkeTNmzfBLd0nxa+OTLdoe/funWKVpGlPyvAtO6XUSmC1iKyNt6wT8LpOSLZz8eJFnnvuOQAOHDhApUqVrH6OiRMnsnfvXvM7AgA//vgjs2bN4tixYxQvXpzg4GB8fHwIDAykVatWXLt2jVu3buHu7o6Tk/Hd6ooVK+Ln58f777/Pzp07qVWrFmfOnDG/A5WcqKgoduzYQWxsbKJ1VatWpXLlyun+bD179sTLy4vjx49Tv3597t27x3PPPccrr7yS7mPGJyJMmjSJjz76CDc3twTrRo8ezbVr11iwYIF5Ww8PD1xcXHjxxReT/QPj4sWL+Pv7M2nSJItiiImJwdvbm4ULF9K0qfFR77///kuVKlXYsGEDL7zwQpL7BQYGUr9+fRwcHLhw4QKdOnWiXr16TJ48OdlzDRs2jE2bNvHBBx8kijckJCTRbVpNS4o1EtI9wF0Szxh7R0Qy/acwtySkPn36UKZMGRwdHbl48WKy3Z0zolu3buzduzfBssjISO7evcsnn3zCjBkzzLOmDhs2zLzNjBkzuHPnDsuWGR/xSVwX3379+vH2228zcOBAnJ2dLaqSbCEgIIDGjRtz4cIFQkJCaN68OVFRUZQpU4YzZ86kWAk8KTY2lrfffpvRo0dTpUoV8/JNmzbRvn17Jk6cyNChQxPs4+HhQYECBcy/pO/fv88///xDuXLlCAsLY+fOnWzatInhw4cneJbSs2dPfvzxR/744w+ef/75VGPbvHkz7du355lnnkmw/Ny5cxQrVoxr164l+bztnXfewdPTk+DgYESE33//nfv373Pu3DmKFy+e5Ll++OGHJLvzOzs7M3bsWIoUKZJqvJqWrmdI8RtwBBjwxLKPgaOW7G/tlhueIQUFBUnRokXlzp075rmDAgMDbX5eg8EgVatWFX9/f3F3d5ebN29K8+bN5emnn07UGjRoYN7vyRG+r169Ku7u7ml6lmRNb7/9tvj5+Zm/r1ChgrRq1UqaNGkiP/74Y5qOtW7dOnF2dpZu3bqJiMjDhw8lNjZWfHx8xM/PL8lpM4KDg+XEiRNy4sQJOX78uNSuXVsmTpwoV65ckVGjRkmtWrUEkK1bt5r3OX/+vHh4eMj06dOlZcuWCY5nMBjk9u3biWKLiYmRU6dOJWgrVqyQkiVLSq1atWTNmjWJ9gkICDA/+wkMDBQnJyeZM2eOfPjhhzJkyJA0XRtNSyusMLhqbYy94P4BDsb9ewmoY8n+1m65ISG9++67Mnr0aPP3Y8eOlXfeecfm5121apXUr19fYmNj5cMPP5TPP/881X1Mg50OGzZMjh07Zm6vvvqqDBo0yOoxpvYw3vSL3TTY54ULF6Rw4cJSrFgx+fXXX8XLy0tiYmIsOpdp9Illy5ZJsWLF5MCBA+Lh4SETJkyQZ599VgwGQ6qTC27fvl2qVKliPmdQUJAopeSzzz6Thg0bmj+PKYk+evRInn76admzZ4/5GFOnTpUKFSokGIU8Oe3atZO5c+fKxo0bkxx4tmfPnjJ27FgREdm1a5e4ubnJqFGjJCQkRIoUKSL//vuvRddG09IjwwnJeAycgeeBN4AXAGdL97V2y+kJKSgoSFxdXeXYsWNy6dIluXTpkhw/flzc3NxsWiWZqqPffvtNRERCQkLMVVJKoqOjpVGjRlKzZs1EzdpJdNKkSQlmhU3KxIkTRSklefPmlbx584qDg4M4OjoKIBs3bkxTlbRu3TqpU6eOxMbGysSJE6V69epSvnx5cXNzk7Vr14pIypMLxsbGSuPGjeXjjz+WHTt2yI4dO+Ttt9+WcuXKSa9evcTb21u2bt2aKIkuXLjQXCWZphOpVq1aqtNzHDx4UMqWLZugiotfJa1evVoKFy4soaGhYjAYpFmzZjJp0iRzxaSrJM3WrJKQslLL6Qlp7dq1Uq5cuSRbWro8p1X86sjE0iopM9y9e1eKFi0q7u7uCabaeFJsbKxERUVJVFSUnD59WooWLSrXr183jzKxY8cOi6okU3W0YcMGETEmaKWUjBkzRvLkySM7d+40b5tclWQwGKRDhw7SokULadGihTz//PPi7OwsDRs2lA8++ECWL18uDRs2THSLMX6VNHXqVOncubPs27cv1SrJVB2ZxK+SYmJipHjx4gKIUkqABG3WrFmyZcsWcXV1tbhKio2NlV69esmpU6cs2l7TUkpIWeLF2LTKLZ0aMlvDhg25ceNGgpEU7t69S1BQEHfv3rX7+GO+vr5cvnwZLy8vTp8+nWC6jeSMGDGCGTNmJJjmQkQICgpi586dtGjRItl9f/75Z7788ksOHz6MUgpfX182btxISEgI+fLlIzo62jziwZUrVzhz5gz//vtvih0mxowZw9WrV8297wwGA97e3ly8eJHatWvj6upq3jYwMBAfHx8OHjzItm3bePbZZ2ndujVvvPEGffr0SXTsw4cP06BBA7766itzF3oRYcSIESxatIiHDx8yZ84cHjx4wLhx47hx4wYzZ85MMs4FCxbQqFGj5C9snN9++41XX32VDh06sGbNmlS317QsP1JDWumEZBvnz5/nxo0biZa7uLhQt25du75ZHxYWRqVKlThw4AAlSpRI8O5USiIjI7l69Wqi5UopKlasmOKID23atOHQoUMUK1YMg8FAcHAwJUqU4Pr168ydOxdfX1/69u1rHlzW1dWVDh06JHudoqOjKV68OLVq1Uowt9KJEydwcXFJskfitm3bCAgIMP+y379/Pz169CAgICDRQLinT59OdnSGrl278v777zN79mwiIiJSnSbEEiJCgwYN+PDDDxk6dChbt25N8r0xTYtPJyQt2zNVR4sXG9/XnjhxosVVUnqFhoZy+/ZtwDjy9enTpxk5ciR58uShTJkyLF++nO3btyc5F1VSYmNj2bRpE5GRkYnWlS1bNlFF8uDBA/P7XfF/0ffp04fPPvssySopOcuXL2fu3Lns3bsXEaFOnTqMGzeODh06WHyMJ/32228MGTKEEydOMHPmTP78809dJWmpsmpCUkq9iHHm2IMikvilhEygE1LuEhYWhqenJ3v37jW/HBwREYG3tzd79+5NtUqyhu7du/PXX38lWl6gQAEOHTpEgQIFrH7O4OBgevTokeQYhJ06deLzzz+36Dimyf3iTz1hyWSKKTFVR4MHD+b1118nMjIST09PXSVpqUopIaU4hblSagWwU0QWxH3/BcYx7E4AXyql3heRH6wdsKbFt3//fiIiIqhVq1aidVu3bs2UhLR8+fIU19+5c4cmTZqwefPmZKcgT6unn36a/fv3Z/g4a9eu5cqVK2zbto3t27cDxmrt9OnT/Pbbb7Rr1y7Nx9yyZQsPHjwwTxtfoEABhgwZwrhx43SVpKVbihWSUioE8BGRm3GT5/0LvC8ia5VSL2GcwrxmJsVqpiskLasZNWoUc+bMoUuXLgmGYcoKzp49m+xU5J07d07XkFQdOnRg3759CYaGio6OJjg4mBs3bqQ6ZJSWe6X7lp1SKlxECsZ9XRfYDRQWEYMy1vl3RaSw9UNOmU5IWlZy584dKleuzLZt22jdurV5RtycLDQ0lJs3byZanidPHjw9Pe0QkZZdpPuWHXBbKVVBRC4BzYG/5L/x7FyAxDe3NS2XmTFjBp06daJu3bp88MEHTJgwIctVSdZWtGhRihYtau8wtBwmtVnOFgCblFLTgaEYpxg3eQE4a6vANC07uHPnDnPnzmX48OEAfPrpp6xdu5ZLly7ZNzArGTlyZJIzyWqaLaSYkERkAjAF47BBA0VkRbzVxYBpNoxN07I8U3X09NNPA8bKwVQlZXdnzpzhq6++YsiQIWTH10O07Ee/h6RpGVCjRg1OnTqVaLmHhwe3bt2yQ0TW061bN6pXr87y5cuZPn06bdq0sXdIWg6QkWdIKKUqAL7Ai4AHcBvYCfiKyEXrhalp2c+JEyfsHYJNnDlzhp07d+Lv70+lSpXw9fWldevWdh2tQ8v5Urxlp5R6BjgKFAdGAB3i/vUADset17RcSymVbMvO/Pz8+Oyzz3Bzc6NLly6Eh4ezbds2e4el5XCpdfteD5wQkVFJrPsSqCki1pkPOg30LTtNs50zZ87QrFkzgoKCzFOzr1q1imnTpnHgwIEUx//TtNRk5JbdC8A7yaybBgRnJLCkKKUqAyeBNSLSw9rH1zQtZfPmzSM6OjrBFOoGg4FTp04xdOhQpkyZYsfotJws1RdjgVIi8iCJda7AdRFxs2pASm0D8gOXk0tIukLSNNu5c+cOISEhCZZt2LCBhQsXYjAYuHDhAvny5bNTdFp2l1KFlFrtfQjoncy6XoBVs4JS6k0gDGOnCU3T7MDd3Z1atWqZW/Xq1Vm6dCnff/89tWvXNs/lpGnWltotu1HAVqVUFWANcB0oBbyO8Vae1fqBKqUKYhy4tSVg+bj6mqbZ1LJlyyhbtizNmjXD1dWVjh070rdvX10laVaX2ouxfwKtgZoYq5Zzcf/WBNpaefoJP2ChiFxJbcNbt27h4+Njbjl9mBZNs5eYmBj8/PwYM2YMAD4+PrpK0tLM39/f/PsaYy/tJFn8YqxSKj/gjnFA1UillCMwRkRGW7DvbqBpMqv3Ax8By4DaIvJIKeULVNLPkDTNvr7//nu++uorVqz4b5CWw4cPM3LkSIKCgnSVpKVZhl6MNRGRKCD+XNBOGN9JSjUhiUizVAL8BKgAhMS9v+EKOCqlqopIHUtj1DTNui5fvoxSiu7duydYXqxYMa5evapH9tasKt1DByml8gJRIpLhlxKUUgWAgvEWDcaYoD4QkUTjr+gKSdM0LXuySoWUDKsMhCcikUCk6Xul1H3gYVLJSNM0TcuZUpvCvEUKq/NYORYzEfG11bE1TdO0rCm1CmlhKutDUlmvaZqmaRZJMSGJyNOZFYimaZqWu+lREjVN07QsQSckTdM0LUvQCUnTNE3LEnRC0jRN07IEnZCsSI+plzn0dc4c+jpnDn2d/5MtE9K9e/dSXP/rr7+ma11G903tByu9x7ZlzNlxX1tdZ71vQvo6Z86+ue06A4WSW5EtE1JYWFiK6+31HyE1OiFZZ9/UZMfPlBX3TU1WjDk77puarBhzBj9v4eRWpHssO3tSSt0DAlPYpBCQXBmV0rqM7usB3LbBsW0Zc3bc11bXWe+bkL7OmbNvbrvOlUUkySopWyYkTdM0LefJlrfsNE3TtJxHJyRN0zQtS9AJSdM0TcsSdEJKglIqr1JqoVLqslIqQil1TCn1Urz1LZVS55RSkUqpXUqp8vHWNY9bdk8pdSmJY19SSkUppe7HtW2Z9LGyHFte57htBiqlgpVSD5RSZ5VSXpnwsbIcW11npVS5eD/HpiZKqUGZ+PGyDBv/3qillNobt/4fpVSqM3VnRzohJc0JuAI0xdhbZBTwk1KqglLKA1gXt8wdOAysirfvA2ARMCSF478iIq5xrbUtPkA2YbPrrJTqC/QBXgZcgfak3JMpJ7PJdRaRkHg/x65ADSAWWGvLD5OF2fL3xnJgT9y+TYEPlFIdbPEh7EpEdLOgASeAzkB/4M94y12AKMD7ie1bAZeSOM4loJW9P09Wbda4zhj/0LoCtLT358mqzVo/z09sMwbYZe/PlpWaFX9vRAJV432/Ghhm789n7aYrJAsopUoAXsBpoBpw3LRORB4AQXHLLbVMKXVLKbVNKVXTqsFmY1a8zmXiWnWl1JW423ZjlVL65x2b/Dyb9AS+t0aMOYGVr/NMoKdSylkpVQVoCOywasBZgP4fNBVKKWdgGfC9iJzDePvnyRe+7gFuFh7yLaACUB7YBWxVShW2SrDZmJWvc5m4f1tjvI3UHOiG8RZermaDn2fTcZ8HSgBrrBFndmeD67wR6IKxqjoHLBSRQ1YKN8vQCSkFcX9R/wA8Aj6KW3wfKPjEpgWBCEuOKSL7RSRKRCJFZCIQBjxvnYizJxtc56i4f6eISJiIXALmA+0yHm32ZYuf53jeAdaKyP0MBZkDWPs6K6XcgS3AOCAfUBZoo5T60FoxZxU6ISVDKaWAhRj/6ussIo/jVp0GasbbzgXwjFueHgKoDISardnoOp/H+MtAD0MSx5Y/z0qp/MDr6Nt1trrOFQGDiCwVkRgR+QdYSQ78A0snpOTNA57B2CMuKt7ynzE+m+islMoHjAZOxJXlKKUc4pY7G79V+ZRSeeLWlVNKNVZK5YlbPgTjOFb7M/ODZTFWv84iEomxB9PnSik3pVQZoB/G2x65ldWvczyvYaz0d9n6Q2QDtrjOAXHLusdtVxJ4g3jPpHIMe/eqyIoN4/MdAR5iLLVN7S35ryfMOYy3hnYDFeLt2yxu3/htd9y6ahh73TwAQoGdgI+9P29Ou85x6wti/CsyAmOPu9HEjd2Y25otr3PcNlsBP3t/Tns3G/88twAOYXzudAP4Dihg789s7aYHV9U0TdOyBH3LTtM0TcsSdELSNE3TsgSdkDRN07QsQSckTdM0LUvQCUnTNE3LEnRC0jRN07IEnZA0zQJKqdNKqWaZcJ4lSqlHSc2Jkxnizh+llPrHHufXcjedkDQNeGKSuViVcBLFt0SkmojszqRwpohIhUw6VwIi0gt4KbXtNM0WnOwdgKZlBWKcYA4wzuoL9BWRHDe8f3LixmBTIhJr71i03EtXSJpmAWWcer5V3Ne+SqnVSqkflXGq6pNKKS+l1DCl1M24OZhax9u3kDJObX1dKXVVKfWlUsrRwvO+rpQ68sSyQUqpX+K+zquUmqqUClFK/auU+jZusFOUUkWUUhvj5t66G/d1mXjH2a2UGq+U2o9xAriKGb5QmpYBOiFpWvq8gnGKgSLAMYzjuTkAT2GcJmB+vG2/B2KASkBtjPM09bXwPBuAp5VSz8Rb1iPu3ACTMU4CVyvu+E9hHLePuHgWYxxjrRzGMdRmP3H8tzHOZuoGXLYwJk2zCZ2QNC199orIVhGJwTiddDFgkhinG1gJVFBKFY6bNfQl4BMReSAiN4EZwJuWnEREojGOXN4DQClVDeMEjxvjbrP1Az4VkTsiEgFMMB1bREJFZK0Y596KAMYDTZ84xRIROS3GaQ0eo2l2pJ8haVr6/Bvv6yjgtogY4n0PxllCS2OcUuC6MX8Axj8Er6ThXN8DK5RSIzFWND+JSLRSqjhQADgS79gKcARQShXAmPzaYqzkANyUUo7xYk1LHJpmUzohaZptXQGiAY+4airNROSAUuoRxpmFu8c1gNsYk181EbmaxK6DgCpAfRG5oZSqhfH2YvwJIfVw/1qWoW/ZaZoNich1YBswTSlVMG6CNU+l1JO3zlKzFOPznxgR2Rd37FiM8+LMiKuWUEo9pZRqE7ePG8aEFaaM02CPscJH0jSb0QlJ02yvJ5AHOAPcBdYApdJ4jB+A6vzXmcHkC+ACcEApFQ7swFgVAcwE8mOspA4AW9IRu6ZlGj1Bn6ZlIUqp74BuwL8i4hlveX7gJlBHRAJteP6FwOvATRGpZKvzaFpSdELStGxAKfUZ0F5EWtg7Fk2zFd2pQdOyuLiRIxTQ0b6RaJpt6QpJ0zRNyxJ0pwZN0zQtS9AJSdM0TcsSdELSNE3TsgSdkDRN07QsQSckTdM0LUvQCUnTNE3LEv4fdY/qtIP6EREAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] @@ -72,20 +75,29 @@ "source": [ "# plot\n", "fig, ax = plt.subplots(nrows=1, ncols=1, figsize=[6, 3])\n", - "ax.scatter(dates, dis, marker='^', s=6**2, facecolors='none', edgecolors='k', linewidth=1.)\n", + "ax.scatter(dates, dis*100, marker='^', s=6**2, facecolors='none', edgecolors='k', linewidth=1.)\n", "\n", "# axis format\n", - "pp.auto_adjust_xaxis_date(ax, dates, fontsize=font_size)\n", - "ax.tick_params(which='both', direction='in', labelsize=font_size, bottom=True, top=True, left=True, right=True)\n", - "ax.set_xlabel('Time [yr]', fontsize=font_size)\n", - "ax.set_ylabel('LOS displacement [m]', fontsize=font_size)\n", + "ax.tick_params(which='both', direction='in', bottom=True, top=True, left=True, right=True)\n", + "pp.auto_adjust_xaxis_date(ax, dates)\n", + "ax.set_ylim(-5, 5)\n", + "ax.set_xlabel('Time [year]')\n", + "ax.set_ylabel('LOS displacement [cm]')\n", + "fig.tight_layout()\n", "\n", "# output\n", - "out_file = 'test.png'\n", - "#plt.savefig(out_file, bbox_inches='tight', transparent=True, dpi=fig_dpi)\n", - "print('save to file: '+out_file)\n", + "out_file = 'dis_ts.png'\n", + "#plt.savefig(out_file, bbox_inches='tight', transparent=True, dpi=300)\n", + "print(f'save to file: {out_file}')\n", "plt.show()" ] + }, + { + "cell_type": "code", + "execution_count": null, + "metadata": {}, + "outputs": [], + "source": [] } ], "metadata": { @@ -104,9 +116,9 @@ "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", - "version": "3.6.7" + "version": "3.7.9" } }, "nbformat": 4, - "nbformat_minor": 2 + "nbformat_minor": 4 }