Skip to content

Instantly share code, notes, and snippets.

@SusanEAllen
Last active June 9, 2020 00:26
Show Gist options
  • Save SusanEAllen/4e82270c2f65d12f8320d0b6ddbcfd67 to your computer and use it in GitHub Desktop.
Save SusanEAllen/4e82270c2f65d12f8320d0b6ddbcfd67 to your computer and use it in GitHub Desktop.
Test notebook for lat lon selector
Display the source blob
Display the rendered blob
Raw
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"import datetime\n",
"import importlib\n",
"import matplotlib.pyplot as plt\n",
"import numpy\n",
"from pathlib import Path\n",
"import rasterio\n",
"from rasterio.plot import show as rio_show\n",
"import xarray as xr\n",
"\n",
"import moad_tools.midoss.random_oil_spills as ros\n",
"\n",
"%matplotlib inline"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<module 'moad_tools.midoss.random_oil_spills' from '/Users/sallen/Documents/MIDOSS/moad_tools/moad_tools/midoss/random_oil_spills.py'>"
]
},
"execution_count": 2,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"importlib.reload(ros)"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [],
"source": [
"random_seed = 101\n",
"random_generator = numpy.random.default_rng(random_seed)"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [],
"source": [
"geotiffs_dir = Path('/Users/sallen/Documents/MIDOSS/ShipTrackDensityGeoTIFFs/')\n",
"spill_date_hour = datetime.datetime(2015, 9, 11, 13)\n",
"dataset = rasterio.open(geotiffs_dir / f\"all_2018_{spill_date_hour.month:02.0f}.tif\")"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [],
"source": [
"mesh = xr.open_dataset('/Users/sallen/Documents/MEOPAR/grid/mesh_mask201702.nc')"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [],
"source": [
"watermask = numpy.load('gifmesh.npy')"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"48.192264556884766 -123.41216278076172 281 473 7.7146745\n",
"48.85030746459961 -123.05535888671875 237 497 8.250972\n",
"48.22787094116211 -123.74916076660156 279 451 9.837744\n",
"48.12992477416992 -123.43057250976562 285 472 50.48799\n",
"48.44931411743164 -122.96807098388672 264 503 1.6171553\n",
"50.3343620300293 -125.4255142211914 139 339 13.6485195\n",
"48.88280487060547 -123.30116271972656 235 481 26.81064\n",
"48.7255744934082 -122.51237487792969 246 533 119.27002\n",
"49.04998016357422 -123.25282287597656 224 484 0.40582165\n",
"49.309959411621094 -123.16914367675781 207 489 36.488113\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAABIcAAAFlCAYAAABxxYi1AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nOzde3zb5X3//dclyfLZOdomcRKSkIRTIBQCgdJTSiEt64l09EjbbXSsHazrGHeb7NDTvS482o7f2l/ZOu6ylQ3o1rKkZS1rOJRCC6UkoQTCMZCE4DixncSJHZ8lXfcf1/cbybJkS7Zsnd7Px8PIkr6HS5KJrbc+1+cy1lpERERERERERKQ8BfI9ABERERERERERyR+FQyIiIiIiIiIiZUzhkIiIiIiIiIhIGVM4JCIiIiIiIiJSxhQOiYiIiIiIiIiUMYVDIiIiIiIiIiJlLJTvAaQyd+5cu3jx4nwPQ/Lo5R17Rt224oKleRiJSGHZsWPHYWttY77HISIiUqjCoRpbHZ6R1zHY/sG8nl9EJJ0eulK+nyjIcGjx4sVs374938OQaXZ54OqT368xS0Zv8BQ8EPvRNI5IpPAYY17L9xhEREQKWXV4BpcsuzavY4jueimv5xcRSedBe0/K9xMFGQ5JeUkMhURERERERERkeikckryYaCB0eeBqVQ+JiIiIiIiI5FBGDamNMfuMMc8aY542xmz3brvaGPOcMSZmjFk9xr7vNMa8ZIx5xRizIVcDl+J0eeDqSVcKqdJIREREREREJHeyqRxaa609nHB9F7Ae+Jd0OxhjgsCtwOVAK7DNGHOvtfb5iQxWileuAx1VEImIiIiIiIjkxoSnlVlrXwAwxoy12UXAK9baPd62/wm8D1A4VAZU4SMiIiIiIiJS+DKaVgZY4H5jzA5jzHVZHL8FeD3heqt3m5SwXEwdy/Q8IiIiIiIiIjI5mVYOXWqtbTPGNAEPGGNetNY+msF+qcqKbMoNXeh0HcCiRYsyHJYUinwFNZpeJiIiIiIiIjI5GVUOWWvbvMsOYAtuulgmWoGFCdcXAG1pznGbtXa1tXZ1Y2NjhoeXfJuuKqHxxiAiIiIiIiIiEzNu5ZAxphYIWGt7vO+vAL6a4fG3AcuNMUuAA8CHgY9OdLBSGBTGiIiIiIiIiJSOTCqHmoFfG2N2Ak8CP7PW/twYc5UxphW4BPiZMWYrgDFmvjHmPgBrbQS4AdgKvAD80Fr73FQ8EJl6hVAllE6hjktERERERESk0I1bOeStNLYqxe1bcFPMkm9vA65MuH4fcN/khin5Eju0gnXzR738BUn9h0RERERERESyl+lqZVJmCrlKaCzFOGYRERERERGRfMp0tTIpE8nhyrr5q9jatnPEdREREREREREpHQqHZMxqm8RgyL9e6AGRppeJiIiIiIiIZE7TysrYRKaOFXow5NP0MhEREREREZHMqHKozGQamiRXDEHxBEM+VRCJiIiIiIiIjE/hUJmYTCVNsYVCIiIiIiIiIpI5hUMlLFdTq4qhz1A6qh4SERERERERGZvCoRI0mVCoFKaTJVNAJCIiIiIiIpKewqESMtFQKFUgVGoUEImIiIiIiIikpnCoyE3VqlzFXi0kIiIiIiIiIpnRUvZFaiLL0KdSitPI0tHy9iIiIiIiIiKjqXKoiExFuOEHQeUwtQw0vUxEREREREQkmSqHikCuqoTGklgtVOpBkSqIREREREREROJUOVSgpjvASA6E/OulOsVMRERERERERBxVDhWY6agSSmXd/FUpg6BSrSJS9ZCIiIiIiIiIo3CoQOQrFEqWKiRSQCQiIiIiIiJSuhQO5dHlgauJHVqhkCKP9NyLiIiIiIhIuVPPoTxIDCTWzV9VUP19SrVKSERERERERERSUzg0TcaqUCmEUMiXGFb510udlrcXERERERGRcqZpZVNsvF5CW9t2qlqnAGh6mYiIiIiIiJQrVQ5NkfECIV+hVeaUc1ClCiIRERGR4mf7B4nueimvYwiuPD2v5y8E+X4NRCQ7CodyKJvqk0ILhdLZ2razaMYqIiIiIiIiItnTtLIcyGYZ+kIPW5LHVshjnQqaXiYiIiIiIiLlRpVDEzSREKFYpmwlN6UuN5peJiIiIiIiIuVElUNZyqZKqFiVczDkK/XXWERERERERMSnyqEM5CooKPeKHBEREREREREpPKocGkOuq4SKIRhKHmMxjHmqqHpIREREREREyoEqh1Iot1AgXQBUbs2oU1H/IRERERERESl1Coc85RYIwfhVQYW+stp0UUAkIiIiIiIipazsw6FyDIV8yT2Q/OsKhERERERERETKR0bhkDFmH9ADRIGItXa1MWY28F/AYmAf8EFrbVcm++Zi4JNRLIGQH9xMZ1ijYCg1VQ+JiIiIiIhIqcqmIfVaa+15CeHOBuAha+1y4CHveqb75kU+l6HPNuiZ6mBoa9vOsm42PRHFEiqKiIiIFBpjzD5jzLPGmKeNMdu922YbYx4wxuz2Lmel2fedxpiXjDGvGGPGes8hIiITNJnVyt4H3OF9fwfw/skPZ2rkMxSCia345YdCuQ5wFApNjgIiERERkQnL+sNmY0wQuBV4F3AW8BFjzFnTNWARkXKRaThkgfuNMTuMMdd5tzVbaw8CeJdNWew7ijHmOmPMdmPM9s7OzkzHn1bs0Apih1bk5M18PsOUsSqHMh2XHwiNtSqZppOJiIiIyDTL5MPmi4BXrLV7rLVDwH96+4mISA5l2pD6UmttmzGmCXjAGPNiFucYta+19tHkjay1twG3AaxevdpmcfwRYodWjLg+0QbLyUFK8vWJHDNXAUw2U86KdZn6fPRbypT6D4mIiIhkzf/A2AL/4v3tP+LDZu/9QrIW4PWE663AmuSNvA+hrwOooibXYxcRKXkZhUPW2jbvssMYswWX4LcbY+Z5/5DPAzqy2HdUODQZidVBW9tG359NQJJJNU42gcVEQw5/v3Th1kRDk0IMW4qRAiIRERGRrEz0w2aT4rZRHyQnftDcYGZP+INmEZFyNW44ZIypBQLW2h7v+yuArwL3Ap8EbvYuf5LFvjmRaspY8vLsYxlru/Gmc01lyJIYKCVOB8vmsRWrxMdX6EGWAiIRERGRzEziw+ZWYGHC9QVAio+DRURkMjKpHGoGthhj/O3vttb+3BizDfihMeZaYD9wNYAxZj7wPWvtlen2neygx+sjlBgqTKYZ9Fj3ZxIQJZ471TjGGqd/X2Jj6lTHK/QAJRPjPTciIiIiUrwm82EzsA1YboxZAhwAPgx8dFoGLiJSRsYNh6y1e4BR79SttUeAy1Lc3gZcOda+EzHRxtKZTB2bTP+gdCFRJuHNRMaRvIrZZKqYpiJgynQ8pRIIqXpIREREZFwT/rDZWhsxxtwAbAWCwL9aa5/Ly6MQESlhmTakzpupWDo8eXpWqpAkm9BlMkHHZHsHpauMyrR3kl+RlOuAKPEc2d5fbBQQiYiIiKQ3mQ+bvev3AfdN5RhFRMpdQYZDL+/YMyWhUKJU4cp4K5SlOsZ0hBtjVfgkBl3J3yfum27/qRjneLdNx1hEREREREREJDMFGQ5Np/EqcMbaJ1vZTCPLZtvkx5Bu26lspD2Zx1YKVD0kIiIiIiIixarsw6FEyYFGulXCxuszlM05xtou2xAlVR+kVNVEqc6Ty+Aok4bepUgBkcj0mDt3rl28eHG+hyFSHnq9RaFq5+d3HJKxHTt2HLbWNuZ7HCIiUlwUDnnSBSepvk9eOSwb6YKatXY/17KLRvropIbGW+uBhhH7TSRUGS/UykUVTylWAk2UAiKRqbd48WK2b9+e72GIlK4nvpT69ou/Mr3jkAkxxryW7zGIiEjxKchwaMUFS+GpfI8ivVxMK0sMatba/dzIDqqIAtBMH/amfizA+oasz5mq6miqqnlKrbm0iIhIWUoXCCVvo4BIRESkJBVkOATwQOxHU96UerqlWoYewFy4F9MaHbGt6bd0XN/DNTe8ecS2kw1gEkOpyVT8qFooPVUPiYhIwcskDEq3nwIiERGRklOw4VApS56ixoFIyu0a6Uu5z3jGCm4m0stovOOrYmg0BUQiIlKQJhoKiYiISEkr6HCoFKuHkq2bv4o77QGaE4IgXyc1KfcZq4ooXV8h//ZU08DS9VDKtOl2rqqaSo0CIhERKQi5DoRUPSQiIlJyCjocgvIIiG5n5YieQwADBLmdlcDYlTrJYdBYDacz7Q+UuILZWMZq0q2gSEREJE+mozpIAZGIiEhJCeR7AJko9eqLh80ibuEC7IIQ1oBdEOIWLuBhs2jUtsmhy3jTxBKXtE93jHTbZxLwqGpobKUebIqISIF44kvxLxEREZEsFXzlULl42Cxiw7au+PX58WBossGLHyAlVgSlakqd6jypKpGSg6bEgCix2bYCI0fTy0REJOcKIQTyx6AKIhERkaJXFJVDUPrVQzCyWscPcTJtHj3Wdqn6DeVqGth4U97EUQWRiIhMWqFWBxXaeERERCRrRVU5VA79h1LJNMQZr1on3RS0TIOhsY6f7nZVEImIiEyCgheRohTd9VJez//qP16c1/MDrPje6fkeQkHI98+CSKaKKhwqR5kEK+kaVKfqB5SLpewTz5VqjLk8RynR9DIRMMbsA3qAKBCx1q42xnwDeA8wBLwK/KG19lj+RimSR8UaBqlBtYiISFErmmllPr25HlvySmNjrV7mT2ObTJCjqqDslGPlm0gKa62151lrV3vXHwBWWmvPBV4GNuZvaCJ5UKjTxbJV7OMXEREpY0UXDkHpB0SZrhSWyXF8uVpVLJNeSFrSfmwKiERGstbeb62NeFefABbkczwi00qBioiIiBSAogyHoPQDosRVv7KVPK0sXd+fbAOjTEOrYmtKXQxjFCkhFrjfGLPDGHNdivv/CPjfaR6TyPRJrBIqxWCoFB+TiIhIGSjacEjG54c0iVPNkoOi5GAkuSooVZVQNvsXS+VQNivD5YKqh6SMXWqtPR94F3C9MeYt/h3GmL8GIsBdyTsZY64zxmw3xmzv7OycvtGK5EIph0GplMvjFBERKSFF3ZC6lFcvm8qmzqkCnXS9icZanSzVcYrRVEy/y4QaVEs5sta2eZcdxpgtwEXAo8aYTwLvBi6z1toU+90G3AawevXqUfeLFBSFI2pQLSIiUmSKOhyC0g6IJiJVQJMutEkVgCSGPtkscZ/p8YvFdFU/KSCScmKMqQUC1toe7/srgK8aY94JfAF4q7W2L6+DFJkoBUIiIiJSxIo+HILSDYgmWj002Z4/idPQ0gUj4x23WIOhVM/ddFYTiZS4ZmCLMQbc75+7rbU/N8a8AlQCD3j3PWGt/XT+himSIQVCY1P1kIiISNEoiXCoFCUGM2OFNP796fZNZbyQY7zgJ9X9pRicJIdz4z2vE6XqISkX1to9wKj/iay1y/IwHJGJUSCUHf/5UkgkIiJS0EqmIXWpvblOFUxkI12D5WzDjcTjjHXMYu03NJ7kFdqm6nGWYuWbiEjJKLeG0iIiIlJ2SqpyqNSmlyWHEpOd3jSZqpexQpFSDYbSmaoKIhERKSAKgnJLU8xEREQKWslUDpWTVEvLZxJWpFqmPhUFH6OlqiDKdShWSsGmiEjRUoWQiIiIlKGSqhyC0qse8iX2H0q8TAwsxpvelXyMyYZA5R4iqf+QiEiJUBg0PVQ9JCIiUrBKsnKolN9c57IHTiZBUrrbyy0YSq4Umurl7UVEZIqpj1B+6PkWEREpSCVXOeQr1QqiZLkMKTKtOiq3YEhEREqAQonCoQoiERGRgpNR5ZAxZp8x5lljzNPGmO3ebbONMQ8YY3Z7l7PS7PtOY8xLxphXjDEbcjn4cpY8RSwTY/XJyfQ4EzlvsfKfr1Q9nqZaOQSbIiJTTtVBhUuviYiISEHJpnJorbX2cML1DcBD1tqbvdBnA/CFxB2MMUHgVuByoBXYZoy511r7/CTHnZFyqR6aTDVP8lSpbEOiUpTqOcjH41X/IRGRCVDoICIiIpK1yfQceh9wh/f9HcD7U2xzEfCKtXaPtXYI+E9vv2lTym+uUwUW2YYYmfYvmorVuQpJqiohv79TPoOwcgg3RUQmTRVCxUmvl4iISMHINByywP3GmB3GmOu825qttQcBvMumFPu1AK8nXG/1bptWpRwQ+SYa3JRy4JOJdKFXKVdGiYiUBAVCpUGvn4iISEHIdFrZpdbaNmNME/CAMebFDPczKW6zKTd0odN1AIsWLcrw8JJoskFPJv2E/PtKITzJRy+hidD0MhGRBAoTRERERHIuo8oha22bd9kBbMFNF2s3xswD8C47UuzaCixMuL4AaEtzjtustauttasbGxszfwQZKtU311MRaGQylaqUKo7yPXUsE5peJiJlTVVCpU2vrYiISN6NGw4ZY2qNMfX+98AVwC7gXuCT3mafBH6SYvdtwHJjzBJjTBj4sLdfXiggGlviFKtSCn/GUyyPVQGRiJQVBUIiIiIi0yaTaWXNwBZjjL/93dbanxtjtgE/NMZcC+wHrgYwxswHvmetvdJaGzHG3ABsBYLAv1prn5uKByK5k2lYUuxTzLJZnU1ERKaBgqDy9sSX4OKv5HsUIiIiZWnccMhauwcY9e7fWnsEuCzF7W3AlQnX7wPum9wwc6dUl7fPpF/QZI+bbmWzrW07YXM3ndf30EgfndRwOyt52BR+76hiC4jUf0hESo4CIUmkgEhERCQvMm1IXVJKNSCaipBjVAiUyuZuzE0dNHu9xpvp40Z2gKUgA6LkoMu/DYqjCkoBkYgUPQVCMhYFRCIiItMu06XsS04pvrkeN9jY3I25cC9m/m7MhXthc3dOzms2HcH0j1yErooo17IrJ8efCumCLr/vUrol7gtFKYabIlLi1ENIREREpGCVZeVQKUs7Tcqr7jkZ4rRG4KYOV+uzvmFyJz0QSXlzk+mb3HGnQKoKoXTT5YqhikhEpKApCJKJUvWQiIjItCrbyiEozeohSL00e6rqHtNvMZuOTP6ELWkyxpZQwVbgpBtTuqCoEKl6SEQKkiqEJFf0MyQiIjJtyjocgtINiCCp8iVNdU/a27NgN87BVpsRtw0QxG6cc/J6IQQtmVYFFVPFkAIiESkYCoREREREilbZh0NQ2gHRSWNU90za+gbsN5uwC0LEgHZquIUL3HQ1v8/RvN3cH7iH++093GnvY63dP/nzZqmYqoJERIqCqoRkqulnS0REZFooHCpxfiCSqrrHVpsR1T2Tsr4Bu20JHFzONeZKt0rZhnbMDe2Y1ggGMFEwxFczy0dAlGisptTFRNVDIjKtFAjJdNPPmoiIyJRTQ2pPqS5vD/GAaK3dz8YFT7mpZC0hFwxNthl1ClvbdrqKoRu6MTb1Nv5qZg8zPUvdZxL4JG9TbNPLyqICTkTyQ2/OJd/UoFpERGRKKRxKUMoBEcDDZhEbtnVNy7nMpiNpgyFfI4WzmlmxVQulooBIRHJKgZAUGgVEIhk77XNP5HsInLhqTb6HQN3uY/keAv0F8DxUb/ltvocgRUDhUJnxq2EmHIZs7nYrnI1XfZRBo+tOaiY2hjTW2v1cyy4a6aOTGm5npZveNgHFVDUkIpIzCoSk0CkgEhERmRLqOZSkXCovJhR+bO7G3NTheghZ3OVNHbC5e9SmHXbs4MdWGxpvrc9+DGmstfu5kR0000eAkX2NxgrCtrbtLOrpZMlKufJNRKaQegiJiIiIlDWFQykoIErNbDqC6R85V8z0W1dJlKTx1noGCI64zfpfC0LYbzbB+oaT4cxkp3Vdyy6qiI64rYqo67E0jmIOg1JRQCQi40psKq1QSIqNfmZFRERyTtPK0ij1/kO+dfNXZR7MtKaZKpZqCtn6Bm65/tQR07wab60fswG2P45US84nBzjJYzbz+1zyNMbYUoVAicFUKYVE6j8khcoYsw/oAaJAxFq72hhzNfBl4EzgImvt9vyNsMTpTbWUCk0vExERySmFQ5JZH6LN3W4d+lQBTEvqH6MNB7tYN//K+A03wNb12a8alnxbyhCnJZQ6vPLGNla4lPz4SykkEilQa621hxOu7wLWA/+Sp/GUvkxCoZgd+W984r/5ATPy9nSiMYjEIBwaezsRERERKSgKh8ZQLtVDmUi3+pg1uKbUSdIFLJNuiO3tm1zxtKn1fG5kx4ipZQMEuaX1fB6evyjl/mONOd02xUTVQ1IsrLUvABijNCGnfv23YK0Ldiwu/BmOuAAHA6GAuy9mIRKFwQhUVkDQuy0ag2Hv39RgAGIWe+QEQ9EolTWV0FANM6ohaqEiCLEYvNrB1+79Hc0zqvnUunOgeUb8HJVl/ieH/ztUP+a5o+ohERGRnFHPoXGU05vrdfNXnfwaJd3qY5aUU8XGC38SzzWRECb5+A+bRYRvnev6GRlop4ZbuGDUamWZhj7FHgz5FG5KAbLA/caYHcaY6zLdyRhznTFmuzFme2dn5xQOr8g98SX4zZfg8S+6UMYPZqyXTISCLugZisDRE9DeDT0DgHFhT20YqsJQUwn1VTCr1t1eVQGVIczMGjq7B7jnsd388L6d8NoRGBx2x7fAojn89YfXEIlZzvzcXbzrs//BPT98kmMvtOXzWSkM/UMuQEtVgSsTp6mSIiIiOVHmH+NlphwriEb1Iko3bWtB+h+hbKZppesHlCm3bQNX3PBmd4M+mRUpVJdaa9uMMU3AA8aYF621j463k7X2NuA2gNWrV+vtdbLEN8gGMF61kDGcLPu0FiJeiBMKuNAnZl11UE8/9AehJuy+jOHkvLKKoPuyQHWYBQ3V/P7SRp577gA3/3gHZ7XM4r3vPNftZ4FZtXz645fy6cvO4oPffoCrv/MAAN/8yMX85Z+83VUp1VZO7vFGvYqo5H/ro9ZVPhWi6rALh6baUNS9XgX6NEwJVRCJiIhMmiqHJK3EwMZunIOtHvmXpq02KaeUJZvoamTZVhWVSrVPrpVbsCmFzVrb5l12AFuAi/I7oiKVyWpj1psa5vcBGoq6YCZgXFBRFXbTyKoq3OXgMAeeb6N1216OPr0f9nbypj/8/6h/+yY+++Ut7ngx6yqPZtRw9sXL2PDxS2moDvPl7/6CB3/6tAuajIFoDBu1fOvjl9LxnU+wcHYtN/3gCc78yD/x4b++B473uyBnrJjPn9YWte7Ln+oW9S6HIm4qXCShGqdQgyGI93SayiFGY2x/+AU4dNw9NyIiIiIZUjiUoXKaXpbS+gbsN5tOTttKXI4+U7latn4ySj1A8p/b5OdYAZEUAmNMrTGm3v8euALXjFoylcnS8xYXoAx7YdBw1AUpw1EXqvQPQWcP7D9MbP9hju5uZ6ijG2KWlqZ6FiyYxeyZNbxy8BhvOHUuJwYjdHT3c/ddj7PnyT3umOEQVFdAYz1vu3wlX776ImLWctsPfsMjDz5Hz95OAn9wG39x1+M0ntXC/jv+hJe//mHWr15CT/8Qf/WdB3j9yT1wuCfe1yhZIBCfGuc3xB4YduOPRKH1qLscirjKp74hF4AlN9bOlP+85VrUugALXMXWVLHAcJTV5y6AoOHwcwem7lyFKJP/N0RERCQtY+0U/CE0SatXr7bbtxfmSsbl9CZ7qkOcbIKaTMYy3vFKocn0WMZd0Q2FnKXAGLPDWrs63+OYCGPMUly1ELhpzXdba79mjLkK+L9AI3AMeNpauy7dcQr5d8SUyPYNr18p5FfeDEWgb4jB3gEGhmNUVQSorA676qH6KqgIxZtTpzIcZeDlQ7x86DhP7z/C7NpKojHLaU0NnDq3jvpTZrhz1Fby9I59NFSFeeFgF/+9bS+ffvtZXLTiFHcua6GuEjp7uHf7Xp59/Si1lSHOO3Uuly5vpuLUufGG2Iks0DfoKmECxvVDCgbi20VjcOSEuz0UiFcvxbwpceGgdyAzciraYAS6+/n0//k5393wHrddIOD28TafNL/fkzHpn99c8BuOHz3hnuu6SU7bK2aaXlbUvydypcHMtmvMZfkeRtnrv2pNvodA3e5j+R4CJ5bPzPcQqN7y23wPQQrIg/aelL8n1HNI8lbJk0mYkXj/eOMcq8dRqQZD+azCEsmWtXYPMOp/RGvtFuKhkcDEKyCiFrr7oX8Y2zPAcDTGYCRKXWUFlXPrqayrcquGBbOoYKkIUnV2C+ee3cK5EJ/ete8wB4728ugjL7L70HHOX9xIRTDAcHSAHz25h0Vz6vinh57j0PE+5tZVccGSuVTGYtBQzXvftYr3dvfz1PMHONo7yA+eeJXw9r2sbJnNyqWNcMoML6QxruF11LrQww9wEgUD0JRQxVoRckFMJOpCsr4hdz0UdJfDUVd9NByFuiq++8X3u2Na6yqTrHUhUzbPUTp+SDWV091iXmVST7+bIjjZfk7FTv2HREREJkThUJZKsTl1JsHLlNrczZ32Phrpo5MabmflqFXGIPNxpgudSi0gSn4uMnlsWt5epMBNdlrMwBCR9m5CFUFMOEi4uprwnDoXdOSqciVgIBCE5c20AC1rlkLfEH17Otndfpx7n9pH/1CEIycGqK0MsbJlNsGgoWcgQmVlhGOHTzAcjTGrtpLz33AqHOvjtcM9BDC0dvXybw8+x8zaStaeOZ+Zi+a4Rte1lS5g8YMpY0YHLns6Oe1zd/LCzR8ifPopbpwVXiDUN+QClMFhF6DMrEkd/lQEoXcQbDB+v8Vr5h2DQ8dg3sx4ddF4Evp6T1mvoYCB3gH33NRWllcjahEREckZhUMTUKoBkW9ag6LN3ZibOmj2GkQ008eN7ADLpAIiXzbVSeVCAZFIAcpVr5SBYUIB45ag96tkMg0yJqMmTM3KFlatbGHVW8+AIyf40QO7+OCtD3LkxCAfv3Q5BsPMmjDnnTqXmbNrYSiCPdoLwNy6Kg4e76OmMsRVq5fQ2TPAT57aR/iZ11m3cgGzT53jpo35K7EBDFtX8RM00DPADx/fzZ7OHh558SCXn9YUH1vMuv2O97lQqKE6/ePwp60NDLlAKeZVERnjwqKW2ZlXAVm81clSrKqWS1HrqqIGht30vqqKKTxZkVD1kIiISNbUkHqC9OY6N8ymI5j+kX2vqohy7Rg9aica8hRCQ+xcmkzYVWrhpkhRGm+1sYnwp3EIMzsAACAASURBVE+FQy4EaazP3bEzFQpAcwNXr3dT2f/rt6+yeG49/cMR/vO3r/KN+57mkaf303a0l6i1mNpKamfXsuy0Zs6YN4O2Y33sPnScJY0NhEMB7tm2hzv/53d0PHfABTbRmKvkATd1bGAYTgzwwfe+AfvLv+Lyy8929wUDrgpoKAIxy+9eOuSCn/H4PYIiUfdYwiE3FS8UyC4YGo66y6leQS1oXO+oOfVwYnBqmmoXIzWnFhERyYoqhyS/DqReareRvik5naqHRCSvpvoNa1UYwl4/HTvV66aPozJewXL2G5dzNvCBEwM8/2wr//vMftq7+3nXuQup94OeaIxgIMCi2XWcMW8G3f3DPH+gi/6hCK92dPMXdz/ORUub+PP3vMH1JAoFXGVP/zDMrHUVUtFYfDpYV6+bkuZVHL3hsrPGnlr3+lFaD3Sx4Kz5buxVFam397MXk/C99RpPRxMaUBvGPl/EG2uuXqKpDqGKkSqICooxJghsBw5Ya99tjFkFfBeoA/YBH7PWdqfYbx/QA0SBSLk32xYRmSoKhyahFKeXQTxAmZYKm5YQtI4OiMyCEIyxCu9ExzhW0+piNpHHo+llItNkOisYAsYFIsMRV1FTHZ7aVbLGEjTY/73JTeXyVdZx1ptP56zzFvH4Y7vZsmMfl53VQsuZ86C2ilA0Rp11fYVm9g/xxhWn8MbufoaO9rLz9SPc8auXWfSn32dJYwOPfPH9rv+PwU3f6h5yoc6w9zulMuTCMv/hH+8DDMyoHv2cWNj3+lH6hyIMH+jinm17+cgHL4o3xbY2vvJYJBqf2uZXMUVtvFm2P/VtrD5P/lS3XIlZFwjm67UWycyfAy8Afgf57wE3WWsfMcb8EfD/AH+bZt+11trD0zBGEZGypWllk1Qub67XzV81JYGK3TgHWz3yj1lbbbAb56ScArbW7udOex9b7T3cae/j5nmzJnTeUphalvgYJvp4SjHcFCkYuZ4ylqlgwIUUJwbdqlz5nGWUGAz5ggZm1vDGd53LJ97zBh7ffYgf3fs7eLXdBVrGuJAnZjGX3czvfWkz4TPnc+HlK/nOV65i//ev40NrTsN87J/5gy/+N12vdDB8oIuBjm66Xu2AngFXTdTeDfsPw77D0NpFrKObY68dpve5A/DyIdjdDq92MPDcAV557GXaj/cRDBh27j/Kxac1uSlsvUPeamkxb+wBV1FU6U0185tlN1S7IK6qwo09HEwf1ERtPDzKVZZj8Rp1o4bUyTS9rCAYYxYAv4cLhHynA4963z8AfGC6xyUiInGqHJKs5Hxls/UN7n3LpiNuillLCLtxDqyPL0vsn+/mebO4kR1UEQXizatvnncBGw525W5MBSrd8z7Z0E4VRCI5lO83on7o4E9z6u53t9cU4CpWAQNz67j6IxdDaxd/f8+T3PzTp/mTtWeyaE49j+0+BMB9z7zu+gaBe0y1lXzi0hWc3TKLf3roOWb/6ff54EVL+X8/cCHH+oY4PWCY4VcHGa8fT8AQGBiiobqCwNx6d/vhHp7ee5iKYIAXDx7jA285HWbXudDHbyLtVwglTiHzTfT5nIoAJ2hc6BQIxMcscZpeVgj+Efg8kNgIbRfwXuAnwNXAwjT7WuB+Y4wF/sVae9tUDlREpFypcigHSvWNdbrQIecVROsbsNuWYNuWY7ctGREMJdq44KmTwZDPb16dtrJpczfmwr2Y+bsxF+6FzfGp7MVUPTTWWIvpcYiUpKloLD1RsZj7qgi6r2DAVeMMDk9PBVEkFq+yyVQwAKfO4a8+ewXPfO1qdu4/wmfvfIz/+u2rJzcxb/t797V2ExztpW7BLN58+inc+sk38/ymD9I7GGHV39zDj558laqKgKtYamyAOXWuoqe2Ck6dS2DRXDcFKxxiT0cPvYMRmhqqOXfhHGzvEJwYgMHI6PDGb9+U+DURU/UaDEbcwfuHXBNuGa0Q/v8sU8aYdwMd1todSXf9EXC9MWYHLjQaSnOIS6215wPv8rZ/S5rzXGeM2W6M2T6M/j8QEclWxpVDaiI3tlLsPzRWf55UFUQ5rypKlkHz6hFj2NyNuakjvhpaawRu6nB/m3sB1Na2nQXXfyjb5zAX41f1kEiWCvWN5sCwm042FHGVMzOqXSXJYMRdhkNTW0EUMBMPQCqCLH7jMu5/4zJ3jKEIHO7h/t++ym2/fIH/3rYXgM/880NUhgKsOGUmK06ZwfmnzuWnX76Kl146xBkb/osrVy1i7WnN8cbUMesCs2jMfQ9QGWLpJctYGjAQszS2H+el3e2cZi0n2o4xc9EcmFVbeNVWqQxGXKhVW+mtrKal7KXgXAq81xhzJVAFNBhj7rTWXgNcAWCMWYGbdjaKtbbNu+wwxmwBLiI+HS1xu9uA2wAazGwt2ycikqVsppWpidw4Si0g8oOWTJo4+/clbpPzoChN8+pORva08MdtNh2JB0Me029h0xFsmuqkfBrv+ZrqRuEKiEQyVKjBELgqmYFhiAboOtrLrGjMhQa1la6RMriKolw2Lo4l/Dvrr9Q1WQY3vatlFlesX80V6xM+V+rsYc9Lhzh0vI/ewQj/8/Rr7HtwFw1VYQAu//rPuLtngA+uOc1VEPkrj4WCbmx1VfGx+iulDUdZOKeOV9u7OX1Zk3u+snkYgxEXPtWEx94uZlM2wybqhVeZhncWF3gNeo3H66py/7qWIk0vywtr7UZgI4Ax5m249w/XGGOavMAnAPwN7kPnEYwxtUDAWtvjfX8F8NXpG72ISPnIKBxKaCL3NeBG7+bkJnJbSR8OSZFKDCSSQ6JMKoVyWU1kN86BxEogXPPqxm/Ws3X9yAqgdfNXcT+7Ux8oqQIp3yuYpavASg6D/NsSn9NcjbkQK6hECkohB0KJ/KXcG+uZFQ5xvKObGTFvha3ZdW6bgeH4lLPkMCFq3WpfVWNUn0RjbvrYgS6efPkQF71phQs0wqGxV+jKlcZ6ls6tZ6l/Gj9c6eln/YVL+Og/PcSHbn2QD9364Khd582s4Z3nLOTchbNZOKeO5oZqKiuCPPlqB2cvmM25C2a7x+Ef1+81NByFUJrHFrWucqd/CGpmj77f4u4D73kPJt3vhVTJ0/GGoq6xtS9mof24q2hKnPo2OAwzanL3vA9HvSAtN4crOAqICslHjDHXe99vBv4NwBgzH/ietfZKoBnYYlzoHALuttb+PB+DFREpdZlWDqmJXIZKrXpoPInVRamqhiYbDCUHFlu/86sRzas3tZ7PhvVdJ881Yvs0lUa0pP+xn+6QJFUwlHg53va54p9P1UMiCYolEEoUibk39gAzqplhwB7vp/9YPzXBoKukCYdcEDE87LZNDCAMLnioCMUDi2hsZB8ba91+c+q46NLlLhDqH3KXFUnBx1QxSd+HAjCrllMvWcZja05zj6H9OPvajnHwWB8PPtfK//zuNWbWhGk71ssvX2xjb2fPuKf51FvP4MMXL+OSZc3UnDLDPcaasHsOhqPucR/vh7pK95x19cVXMfPFYm47Y1KHbt4UuuH9R6homeWqgGIWjvdBY8KfXf6UvYFhd5xwyB0zlMOKoaEoO3/1EqveuByqS3h6mgKivLHW/hL4pff9t4BvpdimDbjS+34PoE+vRESmwbjhUGITOa8U1PdHwLeNMV8E7mXsJnJtxpgm4AFjzIvW2lHzhI0x1wHXASxatCjLh1FYSjEgSqxmSVe1Mh3Byrob3hy/cgC2HhwZliSOYVPr+Wys3jaq0shunBM/XorKqOmUSRXQVE8nS6aASMpeMYZCvj4vxAkHXVhRE8b0DtJ1rI+acNBVEPkBkV8dMhx1IYPxymSGo+44JwbcNLXKCrdUux86GUaHEYUUJASMG8/iuSw+dS6LDVxy5bn8rd9zKLEqJmZdoDYwDANDcGKQA509fOXH23np4HG+98iLfO+RF08e+kNrTuMfPnIJLafMcAFQfbWr5PGfj8EIdPW6EMmfIhaJ8YOf7aSpoZrLzlvktq2rcuPo6aevq4+a+ioqaisZPtDFwHCUwUiMuY317rVIDNxm1rhKocRAqCY8+UbXUes1LbesWnNaYb2eIiIiMi2MtWP/RWGM2QR8HIjgNZEDNntN5PxtVgB3WmsvGudYXwZOWGu/OdZ2q1evttu3b8/oARSqUguHfOkqXdLdnwuZBE5pz7u5G5NQaWQ3zhm1Glqq8GW6p1iNFxCN9bzmcqyJ4ZoCosJkjNlRqo39MzUlvyOKORBKtO8w9zy2m99fsxRCQY4d66OuqoJQRZDBwWECxlARdFU2BAMuZAgF4pUw4IIScIFQVUV597GJxNzz0tHNC3s7OWvjD0/e9cPr38H7zl9MuLbSVfgEAy58O9bH8cMnePnQMZY0NnBiYJi6qgrmzqpxvYwsnGjvpm7eDFdlVBE6+RwP7+kgZi2V82ex97kDBAKGU09rcuFNyKtaisbc6mvBLF+XaEKPJXDfx7zXPRhwxy/l6WSplGj1kH5PuIbUa8xl+R5G2eu/ak2+h0Dd7mP5HgInls/M9xCo3vLbfA9BCsiD9p6UvyfGDYdGbBxvIvfupCZy3wd+aa3916Ttk5vIPQB8dby5wqUQDkH5BETTZVIhUQbHziT4OhkkfedXmE1HsK0ROqnhdlbysJlYxdtYlVjJY0n1+KYqyFI4VJj0R38Of0eUSiCUqKuPf/iPXzN/Zi3vOLuFoycGmVkTJgYc7Opl3qxa5s2ohoogsaEIgeqwW/I9YFxY0O8FQzXhqV/ZrBh19bJ71wFWfP4/U979p5edxdc/dDG1s2oZ7umnorHBPZfBgKvM8kO44QjMjK+G9uTPn6WyIkh9ZQVLz5rvVpkbijKwp4P/feZ17t/Vyp9dvpKz3nw6dPfFG1BD5s2w/e18/piCgeyDplJSggGRfk8oHCoUCocchUNSaNKFQ9msVpZMTeTGUYrTy2B0UDHlS9h70jVuThxTro495m2buzGfa8cMu7/rm+njJraDZUIBUXLfpvFCqumi6WVSskoxFPJVhrjsrBZ27DvMno5ujDG8dOgYP3t6P9XhIF96/2o3VayqgoBfNRIw7itq471syjksGMusWpa/eQX28S/Cq+08+UoH/UMRdrUeJRwKct2/Pco/PfT8mIe4YuUC/vq95/OWgSG46zfQcZzzGqppv+IcFl51AdRXQSAAVQGqzpzPVcuaeeOyZtq7+6F3wB0kEHCX/ud7wUD6E0ZiMOT13yv3SrBU1H9IREQEyDIcUhO57JVqQAQjQ4vp7ouTfK5UK3lNFfM3nZjhkbeFiXE9T/Mw2YdDmY53up/bdfNXKSCS0lHKgVCiqgrOu3Apy5pm8NqRHvYfOcEjLx7k3t/tY6G/WplfFRQMxBtOm6CbUmStqoUyYYBlzVy0rBmAt3o3//Efen3x/KXmAwlT9o718YV/e5TvPvw83/37e7kAqPX2Cx/vZ+G9T8GyJrfP3Dr3GgVck+/mcxdy7IlXefjhF1h70Wlugv9gxE0LA/daDgzHpwn6YxiOul5CwYALBfXaioiISBpZTSubLqUyrcxXquFQOvmadjZdzLzdKf++tsAV5vcndMx0z1kmYdd0VBYpICocmi6Q5e+IcgmFfBavyXLUhT6DEejs5h+3Pstf3P0bwE19+qv3nE/LimZvtTJ/1auA29evJJLcsYxsgn3Vt6Cje/R2c+t45nPvZN/hHt771jNcM3Br3eszMMwrr7bTNxjh3DPmueCovsrtF/FWTgsFXQiU+PpZ66rCasLx4EhGK6HqIf2e0LSyQqFpZY6mlUmhmYppZZKhUq4eSiUfVUSFYiLT3MZ6nsrxORTJiXILhXyJH/hUBF24MKOazy2YzUcvWc637n+W2375At/9xfN85u1n82eXr+T0M+a5VbD8XRMvlRHlhv88Wlxw05kiGALs4ROce/FpnFtV4a0aN+Rex5owNFSxbE4dtB+n89BxGv0m2NGYa3RdX+3O46/KFrNu1TqMqx4aHIZIwE0tExEREUmicGialFtAVNJmGehKUXE3K/4uKtO+QZmEP4n7TndYpOllUlTKNRBKJbFyxBiorKDpvEV87cx5fO1Tb2Xns61sfeZ1vvfIC1za1sWVqxYRbpkFdZWABb9vsX+pHkS5YXABUdMMaD8++u6mBlf5Ay4Uqq9yrx/G7RsEGutprKqArl7ad7XSfPq8kcvdB4Kuz5CNxgO+ypD7pm/Q+3nQn3+jqPeQiIiUOf11MI3KLSAq1Qoi+3dN4DWkPnlbhXd7GpmuPJZsvGBoqqeUJR5fAZEUJAVCKXhTw2KWeNDjLWNeVQFVYVZdupxVFy6Fnn6e393OfTv38/5ZNe7+YCAeLlmvCoVAfEUz0LSzyTDAZ94Om37qqnk8vUDtH78tvl0g4MIgv3F4NOZCnYqQqxSKRIkd76dv/xFqzpo/OhAMeYFRJOrCI39lskg0P+FQMVSiKSASyalCmMoUzfcAgBu3vJLvIXALH8v3ELjx63fl9fw/OXJ+Xs8P0HZx6srhQqHJ5zLl/EbRJWN9A/Yfm7ELQliDu/zHZljfkNHu/spkY8nkOSup51RkIhQMpZY8PWwEr59QZYULGMIhDLC0qQFqKuPhj9/oGC9kMuQ2GIqlHFz5WHcObHw3nDIDDBwMBvhjoKtltheiJDzHARPvFWS8CqKYhWCAeac1UdNQBa8dhu4BNw2tZ8BtY4xbpSxg3PcnX78M/vTL9ctz5AT0Dub4oFNE/66IiEiZUjg0zcq58sIPPIoqLNrcjblwL2b+bsyFe2FDu7t+QzsA9jvN2G1LMg6GMpHqucl0mtpau5877X1stfdwp72PtXZ/Tsbkn7+cKt+kCOgT/tH8htRRr1LID3aMiYcCMe++gIHaSuqrKth/5AQc63NhgvUqjsCtuDU47KYp5bJRtbWuEqacrTsHtvw5PP5FWv/2ffwAuPPx3fH7/ac6at3zD/GVz04MuNe0OgxNDTCnnmN7OuBAl9vWGDcVMHFFOj+cCZjchz+pWFxQ1TfkxlRbOQ0nzREFRCIiUoYUDuVBOQdEiZLDooILjjZ3Y27qwLRGMBZ3eUf3yOs3dcDmSZYHJgVQqQKd5OckVeXRWrufG9lBM30EgGb6uJEdOQmIkqeXiRQMBUQj+RU+sVg8gPErRkzCF7iQqKqCBasWsWrhHH7ym1d4/FcvQ3v3yACprsprbJxDmVSvlAsLF56zEIDP3vmY60WU2Fg8EnWBkL9c/aHjLsSrr4pPM6sMURMOMRyJuubVfmPqvkHo6qXtqX3sfHKPO0ZlaHqmd/nniMRcgFXoU8pERETKnHoOSUHKJAyZambTEUz/yI9Xk/+2Nf0WNh3BTrRyyA+g/PO0RriRHWDhYbMI8B775m7MhUfgQARaQtiNc9jatpOb583iWnbRZPrwuoKMUEWUa9nFwyya2PjSUP8hKSgXf0Wf9Pv8f7ISp20NR13/Gb9ixK8k8oOkcJCFK1uYXVfJPdv2MPzUPt56yTJoqI7vl2vG+08kpuXVDZx+/R3xqx/4NoO3f4rwmfOZ965vcOh4P+cunM0Z82Yyq7aSnfuP8MSrHdj7/tI9f3WVYIKE583gwMvttLQehTl1rrooFIRQkPmL5jB/Tt3IxtXpRBOqzSarttIFWkOR4muCrf5DIiJSZorsN3XpKLfm1JO1bv6qEcvET0tYdCCS2+1SSBVAVRFl44Kn2LCty92QIkDipg7sk/1srH7F3W7TfyjbSN+ExydSNBQQxQW9sMWfTuT/ExUIuKlGJuBCocSgKByi9rQmLu8f4rnWLhcoDQ57FSzecug5X7HMul40s2pdQORXx5RhWPTSlj9337Qdw3zg21Re+z3ef/5iVpwyk3ecvYC/+/0L6RuM0NHdz1vPmE/PvU9xeO9h5s6udVPLDFBfTVNDFYRDcLw/XvVVEXTP61DEhYbJIY3FVScFAvHKspxNHyT+81iMFBCJiEgZKeLf2MVPlRfZSZxyNi1Tz1oyzE5TbZfcqyjd1LN0wVLC7SkrmPqtm+KWdHsqndSMu81EKNyUiTDG7DPGPGuMedoYs927bbYx5gFjzG7vctaEDq43ce5NfTDgKkYqXNUI1rpgwJ+eFPV60iT2ngmHoLKC+YvmsGLeDFr3H+H43sNw8Bgc7XX7RmM57lVjoLrCHXvI+zevDIOhEebPxN73l/zFunP48VP7ePSlg9z8wTWcelYLZ57dwlsvPo2P/P6F7PrBnzJ3Vo173oYjLgw6eoJgwLB3/xH3nPqvvR/ODEVS93nyb/d7SuUyGIpEXSAVyvG0RBEREcm5Mv8rLP8UEE1ccn+iXPcrshvnYKtH/pGc/L7IVhvsxjkjb0zVqyhdb6J0AVTi7WkCpEz+fLfA7azMYMuJUUAkE7TWWnuetXa1d30D8JC1djnwkHd9Yso9IEqs/AgGXEAUDrk355EYDHkVQX7jaX8lK79xdWMDp56/mAXnLmRGUz3DA8N07j/CsRfaXK+bnn53jKidWFCU0OuagHH9cmI2vsy6wKxabvn87528+hd3Pe41GY+5y95B91VZ4aqu+oegdwB6BghUhZk3s4ahQa8nkbWuIbQBZta4fkSJvFXPqKrIXU+gmHU/IwPDJ/ta5b7qbBqpIlFERMqE/hKTojdlIdH6Buw3m0YuWf/JhpHXv9k0aqWytJU+m46MOkXKACo5cMq0gimFbipO9i7KtXz0gZKS9T7Ab7pyB/D+SR2tnAMiP3w5uUKZcc2kgwF3GfJWrxqKxKuI/JAIb79QwG0/u46KBbNpXNbMzKYGujp6aH/pEPu274XWo9Dd744RSxEUjRUcxRJWUauqiIdCuehxUypqwgzcfi0AP9q2h73PH4BB74MCY6B/GE4McOxAFz2vH3XPZ30V1ISpqq/CYuna3e62rQzF903mvw65ELNuXAPDrmIoHHSVYaXwsiogEhGRMqBwqACoemhq5Cwg2rYE27bcLVl/c/PI66kaUWcwVWzE8ZMDqKTAKVWAlIkBgtzKG7LeL1P+86vqIcmSBe43xuwwxlzn3dZsrT0I4F02Je9kjLnOGLPdGLO9s7Nz/LOUc0AUi41sSI1xgU/iVDOIB0SRmNvn5ObG2y7gQqOhCAxFsNZSFQ5RV1VBR2ePq1jxVzWLWVdN5FcUGVIHRAav75Hxto+5ECEcyq66ZDCSPvAoEZVntWD/+7MALL3pB5x1410ceO2IV70Vgdl1zJw3k/qmBmie4ZpQz6yB+ioqKys43j/E1vt3ub5O1eHRJ4hZF+LkSsybvuiviKZKMBERkaKi39wFQgHR5G1t2znqK2OZ9gjKZJ9MpoolSg6gkgMnP0BKMwy/UOA4YY5TQQxop4ZbuGDKqoaSKSCSLFxqrT0feBdwvTHmLZnsZK29zVq72lq7urGxMbMzlWtAZPECHxu/wQ9x/Dfu4ZALfqIx17NmKBqfauZXHBmvJ5F3nOFojJ37jzAcjdHUMstVquDfb+NT2hKnQKWqKoJ4NYkxI69nqiLoztE/nOM+SAVm/kzsXZ/hnAWzeaHtGOf97T389sWDxAaGvN5RxJe4j3i9pBqqYeEcFq9sYWlTA4/u2Efs+7+C938L3vhVuOpbsPVZr6os5BpR54LBhU3WTwdLjKqHRESkxGm1MikZmVQKpQyM0q0GBi6o2dztpoQlLCMPpN3HbpwDifeRpjdRNtY3wPXtKe8yuDDoGnPlxI+fZK3dz7XsopE+OqnhdlamDZq2tu1k3fxVWt5eMmKtbfMuO4wxW4CLgHZjzDxr7UFjzDygI2cnLMdVzIwB41WFBIx78x8gHhAZ46qHjIk3DPaDBguErLePcT1qwiGor6J5dh311WGqKgLx6WnGAAG3X4CEJdD9cMnGq5L81bDwzuNXCkW9EGm8RsgWV+kUDsaDDb+vTWLVUd+Q68dTzH1uEi1t5Jm7P8Nr2/dyxdd/xre++mNuCQVojsQwc+rgQ2vgHWe716OqIt6/qbaS5ZcuZ/lPniL2b4+65w5c76ibf+q+v+xs73Wb5HMVibmfiYZqN92wKgY1laWXEWn1MhERKWEKhwqIlrefeokBkh8UpesRxN90uK8uG//71guBqCL1PpuOuCoggORAKdUUtGwsCLnzp5BuuXo/5GmijxgGg6WTGn7DKVzCoRHhD8D1PE0DQ+7xeMdopo8b2QGWUQGRHwyJZMIYUwsErLU93vdXAF8F7gU+CdzsXf4kpycut4DIrybBr96JxoMXP0jxtwkGXEDkTy+KxWDYgonFew8F46FOTXMDfe3dVBzuoaI67MKImrCbtmRtPCTCC4ZOhkWMnLqWWK3ih0iZJAmJ4ZK1rrl2Tz/UVrov/zGWmoogp168jJc+egmRb99PKOI9l0dOEL3tlwQt8LYzXKPqk6+H99rd8WsCw0nTxwaG4as/hq9sgaYZ8Jm3w7pzJjY2/7XwVyWrrHABXUXIBXmlRgGRiIiUKIVDBUYB0fTxQ42tdnfqtySJoVAC02+x/WkO6vcVWt+AnWwYlMRunAM3tGNSvO8xC0Js3TYyqFlr93MjO6jCvSkIenMvmunjfewZEf7cxDYCQCjN/IwqolzLLja0dY24PTkYUvWQjKMZ2GJcWBAC7rbW/twYsw34oTHmWmA/oH8EJ8OvEPKXJfend/kCCU2BEvsLRWLxpc6t1z/GryDyp6NVhqipq3KhzNFeV43ir3bm95iJetPMggkVRD5r3XkCNqGqKaH/EN6YrTdef99oDPqGiLR1EaqqgMZ6Vzk0s/ZkTyT6hlxQUREsnaqhRAa46zfxYMgTjETp/vdfc3z5KSw8pcEFRF29LqgJBaDjeOrj+T8T7QmVRBMJiAzxxuThkDtnRRAOHIW59fHph6VEAZGIiJQghUMFSAHR9Flr96e9b0JvLVL1FdrQjrmzG6JAEOw1DXBzc/bHXt+AfbIf/r17RECUOGXNr4ZaN38V17LrZDCULPmxhTNo2tFIq5pQvAAAIABJREFUX0ZVQgqIJB1r7R5g1A+RtfYIcNmUntx/I1fqFUR+tQ5ev6CT/YO8aWZ+sBDy54B5/KlmfiWRv7y89ZpMk3CscCi+AtqJAdeY2uKthhaMhzuRWDygMia+v39eaxNWS0uYcuavoGXx+iFF3G0VQUILZ7vKlFBC5VF1hfvyQzC//07iNqUiTdBT1zvI9fc8ybpzFvD+85dQ11QPj74E//FY/EdiLAPD8N1fjB8ObX3Wbdd+3DXB/vTb4fKVLjj0pykOR+FoL3c++iLXvGtVaYZDIiIiJUjhkJS1a9mVsiv7uBMcZgWwA3b8vkIb2jF3dMePFQXu6HZRzEQCopubsRdVp5yyljxlzszvy2mj1k5qMt5WAZEUrFKfYuY3kTYAduTUropgvMF01LrAyA9ljJ/GJGwXicanDEWi8V42fthTHXbhTmePq9qpr3LTzPyAKBiIh0sB3Pn8fw39qiL/+NZCjJHBVizmKpRqKjMLevxQyl+FqxTDoeYZrmdQkp6aMMGA4cHnDvDs60d5Z+8gb/vNbsxQNPMPOtrTVBj5tj4Lm37qXhOAQ8eJfe1e9jzzOvfVhPnwmmU0nTqHo68fJRgMcM3H3uh+lpKVSr9qVQ+JiEiJKcG/nEqD3lhPj3S9esZiqw327xrHXYIewNzZPepvYOPdPiGpmmN75xzVbDvd6mgTMEDwZF8ikaJX8m/orKu48SuHTgYwxKeQ+VlQLBYPZyAeFBnjesaEAvHr1rqqkOGEJtZe36GX9x9h4NBx2p9v46UnXqX3xYOuqsif4haNxaet+YGRH1SRcE4/lPKrharD2S2JHgy4gGo4Ojocj9p4U+Zi9em3u8eXqLKCGZ96Kx+/dAWRaIzHdh/i3EdexAxl91h7airpeeZ111B60JumF4m55zEag39+KB4MeQLDUUKbt/Pndz5O85/9O6uuv4PP3f04T+3rTN3k2j9WqSjloFlERMqOKocKmKaXTb1OamhOFRClqgwCuglza/95bFjveu+M21co3d/mE3l/Ms6qauvmrxoxrWytnTWi59BEnHzMnJd2tbJ0VD0kBa1UK4hiXllGOBSvoAFXzUPCFDFwwY9fRZS8YpX/bSAAYW9lMj/YiVnoG4wft6GaFUsaGT4xwIyaMAe6evnJU/uoe/4AFyxupKmhioqaSq9HUA1YLwiKEQ9w/Kofv+IoSLxCKVtB73H1D7kgxW/GHTQQLPIGyf60r6SpXebLWwC4cEkjazqOMzvLw1rgkXMWcN0X/ouGqgo+8aYVfPJNp9MyfybMqHZBUXvqDzVOBdac1sRvX+1g562fZLCzm7sff4Vju1qZuXLByAquxN5UIiIiUlAUDklZu52VowIUvzIIcKuPtUZGLef+8Hx3OapaJ1mQ1EHQBN6fpF9VrRM2HeF+dtMxz1t5zHgri1lOrlYGqSv501X4W2ATF2UdCiVSQCQFrRQDopOzwwyE/HAloTLIb1DtVwv5AUzUS2oCCVVCif8yhPzG0TE3PSxmXdWR1wuIGdVUzKihYjhC/dFeljXPIGAM37jvaS5dfgoXLm3i1Dl1GGuhrireoyaxUskPDfxVzQImXl2UTYPpWMxVvzRUx0OyxClnoWBxT2tad86o3kB2adPJ1yF2zXcxvUMpd03377256gLe/bl30tZ+nLt/8Tx3Pb6bv75nGw9+4d1c9oZT3UpwjfVuCmGS7pow//1nV9CyegkEDZWN9fzh8lNcOOeHe371WkLLqZwZjqaevjZdNL1MRERKhD6+KXB6Yz21HjaLuIULaKeGGNBOzcnpYetueDN22xI4uJzGgy1sOOiqhdba/dxp72OrvYfOeQe4ed6stI2a7TUNo2Y2WO/2jG3uxpz1atpl7OmKYVojGOtWHtvAk1xvnzr5+G5nJYMEs/57vIOaSQVDIpIHfuBysiLIWxHMX5Ie4kFRxJvq5YdCQT9A8af92PiXxW3jr4AVDrrgJRp1U5Dw9g8FWb5qIRddsJjV5y7kQ2uW8e+PvcySv7ybf/j5M266mC9mE6Z+2YRV1hKaZSc2qc5UKAgzalyQdWIAegbgWB/0D4+/b+LDfu0I5pKvwuETmZ87X+bPhIWzoSpMoKs35SYW2LG0kd6K4MmH2AnwpffD53/PvaYLZ/PRT76Jn/2fj/Hcpg+y6X9+x9c3b+dY61H40MUu1EsQDQXZceESWs5bFK/QAhfW1FS6aWjRhCAykKNkyALDUXqeeZ0nH3re+xnMo1ILmUVEpCypcqgIaHrZ1HrYLOJhEkKQG7wvEpa79yqE1tr9bKzedrKCp5k+bmQHG77zGpAi8Lm52b33mehqZZu7MZ9rx4zxnib5T+0A8F728Lydy8Nm0ZirloGbNlZJdMQ2uewxpOohKWilWD3kl4f4PYXwAqOAiVcA+T2DrIWhmAuODF4w4zesJv5m3iYFRAHjgiW/MicUjPcvMsZN57JwyZuW8z9nt7D9xYPc8euXOHq0l9l1VW4coeDIRtQm6UHEiAcKsZg3HS2Dx2+AmrBb0r06HH+8g8Mu6KmuiFcrhQLxaWjDURgYiodW1vLtay4t/NW2Ypb/+/1f82fvfYN7vI0N0DF6CthwdQWnf+HdHO4dZOszr/PbPR109Q5ybc8Av/fjp+COX42YqnbWFefw4PJmPv+vj/KJ2x7mU289gzXrL6T5gWfhyAlobODHi2az9oMXe9MWkwQD3qpxMRc85UprF798ah/Hege5ctUiLlp75qjQSkRERLJnrM3i07hpsnr1art9+/Z8D6PgKCDKv/tbfoxJUcHTTg2NB1tyfj5z4d6U58tEOzVcY65kq70nbYngAEFu4QLATT9rpG/UFLpcUUCUG8aYHdba1fkeRz5Nye+IUgmIEits7Mn/eJLCFf8uv2G0z1/mPhaLB0IjQhnvur9cvT9tyz//UCQevERjrqojFqN1Tyc9A8Ocefo892Y+HBq5+hkJ5/KnxqUKgzKpPolZd15/+lpi/yX/sfoNtq11vYiCCZVRwUBxraoVtdDT78b96EvwjZ+NqKaxFUGevWIlXecu4i2nz8PMqSN2pIc7H3+FJ+98jG/0D1Gd+KNSVQEb3g3vOBv2H+H7v3iBp147zPvOX8wZ82bSWF9FuCbMH9z6IMGA4eiJQbb81XugZdbIcfkNzKvDuXsuewZ44IHnuPzS5dDUUFivUYFML9PvCWgws+0ac1m+hyECwGd2v5LvIXDL5z+W7yFw49fvyuv5f3Lk/LyeH6Dt4gkuSpRjD9p7Uv6e0EctIlmwrZGUf4c2mT6uSKoyyokDEy+VbzJ9bkn7C0Mpp6RFMdzCBfE+SkztFDJVEElBK6UKIr/Kx//Hyn/T74cuJNxvvGlnfnVNzAt8jFexE8C7nfjy8/7+EG82HEnoE1RVEV+dDLyqkSALTmui99BxXt93mKi1NNZXUdvU4AKciqALZhIbZge9kMofmz/+TEIbv+IoEnOhSWWFqxjKphlyIYUO4wkaqK921VHvPNe9Lv/8C+g4Do0NmD98C+eet4gX9h3m/l2tLG1qYPmSuXzi3efx0R/9llBfUo+igWHX9Brgn3/BH3Qc5+Mza3mgb5Cn15zGmtOamRsKMH9mDTNrKvn8e98A82aOHpffCyiXz+VQhAsWz4W5dYX3Gqn/kIiIFDGFQ0VE08vyL+3qZi0htm7bmbb30IS1pA52EqV9n+QtZb+p9fwRU+H4/9k78zi5qjrtf8+9tXX1nnR3tiYLJOQNhD0BFFACSNhkiaOC4jIyKhodHQdnyPuq47jBzDAOvCOjMjC+KqOOC7ggGlHCpqDsEGQnELJ2J+n0Vuu997x/nHvqnqquXpL0mj7fz6foqruce+6tIl336ef3/IgcQzZTyGIxOFgEoqBC/NFBzDr0WS8rdSkjKiHTAda67EyLRIJINNIiUZnYpLOCUM9jrjpW3gudOgJqk9QuaKG2J8uTz29n45Y9XFATh2Q9FPwoWNhxolwkoY9lKF1BeKyhhAHdXS3mQtxwKR3MaOeT50eh1X5YThdIiDksS8VxX9jB5t19dPXnWdzWwIy9VX6nAXJHN+IrvyjlRLld/Zzz8CZuL3g8GY9xelc/n7n3OVK9Ofp++QR1a8+CS04YONBohEVLVMB1f56+7d001yYP/vfTYrFYLJZxxv5mnWJY58XEcgvLkTXldySyRiDXzQSUa2g0BSK5biYyPoLtKl7ncJHrZrJ67jFsEPO5JruyLHR7JMLQ+m1Plh6jhRU3LZOeg+Gv/o4TdffSZWZuKLgUvShrqIQsd+RoEUlnC2mhRecS+aFwZJaw6TFdEW0HKvsnlVCD6zDqhhqOOXY+sxpruOupLWp5IhSFvECJRPqYnh+VJvlB9I+ddjlVexR9NU48LBFrSI1u5s1kJl9EnP4VFaRd8NX7kYwp11Q8Bs21HH70IZx6+CwSMYeXOnrwm2qrDhVAeYB4+PqSV3fRtnELwS33ku7N4QB1/Xm8r/4KfvzwwF9I+0qmoMLD/UCdQ64ImbwSGvtyPL1lD6Kt/gAPMgT+AZ7AwSAwWywWi2VaYsUhi2Uf0EKLbI8hBeqn7m5mlJWNmkC0pgF5/Sxks1O6d6tEAMRR24RzStzYAmuigOwNYj6Xi/NYLf6Cy8V5VYUhUwwyBaHRdkNZgcgy6ZnKAlHJ6WN0KdMuG925TIsoWoApyykySsZKIdaGSKTzhUod0ULBJhWPsopKxxJRq/t0UpV2gepwVpNg5aJWfvLIK7z2coeaVyquHD4QHUvnIWn3kzk/85z1Qwdix5yoXG46YQZR7+qFQHLvL59EnPol6OxR74WE5Jxmjj3qENpn1PLim5YSVLp7krHBvyB29nLUn14mbuZUATEvYOu//Rpe312eYbUv5Ioc8Z5v8JPbH4H+ghIHnVCodNVnMO46qhvdWFD0BwpiFovFYrFME6w4NAWx7qGJZYOYz9lbL0ZuW6Ja3YciTKXDZlQFoj8fhty+BHnjLGSVP4CLIlDrlM1pJMcfzh006mVyIVYgskx6prJAVPQjEacUJm2KRaFgYgo8nl8hwBhCjCB0EoVfGbQgpMcxHT2R/SiajxZ6Yg7UJpRIlC9CU5rTl82lL1dUJUO5glruhCJPwYuCqb1wjlrQ0s4nMwg7CCLXU2X49nThsDbkg59DPvg5aKmHvhxvPmUJAH98/DV17epTytHVUMPchS3Unb6MF88/Bq85jQT85lp+eMQ8cg3Vu7RtcQRBZ2/VdXMDycf/8x4e+92foS8ffS7WPw0XXY98wxfgkhvU62pkCuzszvIX/36X+kykE8r15TqQKbBlVx/L5jaNneinhdBcEQY5xxHx0D9YB5HFYrFYphwjFoeEEK4Q4nEhxB3h62OEEA8KIZ4WQvxCCFGljzcIIc4RQjwvhHhJCHH1aE18umMFooln9dxjBhVPtNgy6uLKmobQ61+FivDqwYSfsSgXs1gOSqaqQBR3o4yeUui0IZY4oYBiCkVadKkkCKKOZJKom5feh3B8M4Ba+xz1MU0XkReE4dAJEILZjWklDpnOJFBOo7TRhl6fh3Z26NIxPwi7qsnoPMvyiqYxCRfCbJ6Or72Xk44+BLqzqmwrkGr9zHrajz6Epe96I9s/ewkP/8Ml/M87T+JLnb30nn/cwBbxyRjtf3suufrqwlGmLsUX1qxkd1+e6791H7d++366v343XPML6OhRH8Ed3fD52+Ff7hw4wIxaHvvi29Tzvpx6f3tysLsPEjF6ckWyRb+sE9uIyHvQkx1+O6HOsSSaWiwWi8UyjdgX59AngGeN1zcDV0spjwJuBz5duYMQwgVuBM4FjgAuE0Icsf/TtZhYgWhiGU5gGTOBaN4gOfKDLa9gX+YzlgKSdQ9ZpgRTTSCqdNRAeeC0Li0riUdh2VgstCSWys1keRt7YYyvl+tAYH3MUmt4yoOktWijw7C9cJuCx9Y9/eSKPnJnj3Ka5ItK7MkWwsyc0LEkiXKPJOo4+aJyeBR8I3ybcufTdMd1oClN6wkLYXYTNKVVeVauoESiXEG9920NHHJ0OyceM58TD2ujJuHyt6/t4p6TF+M116rrO6tBtbc//1jSf3U6siLHKe8I/rUuyd9+/0Fe3NnN3OZaXu7ogR//qbqYc9sjVR1Ed23cAoC4/BvUnXktS9/3TdLv+HfE+f/KaV/6mfriWvT37To4IswxGuFnIpBKWDvQj5B1D1ksFotlCjEicUgI0Q6cjxKENEuB+8LndwFvq7LricBLUspXpJQF4AfARfs/XYtl8lJNdBkLh45cN3PIUOzRpjKLyAZUW6YdU0kgEoaYo7uJmTe4erlriDvamaOFGB0orTuamdtotPDiGhlH+uGF7h5pOH60o8cUlRIxin7A4XOaEMlYJGjlvXKHUsyN2t3rMfW2pkiUL4alZ0ZJXOXNvRa9BkOG5+obD53RNNVxwyym2qTK7Ekn1LXV189Xgsjiw2bx0TOP5J0nHUZ+5aH84B0ncuffnMufP3kO+fkzYcseOG4B4orTobUeBOQaarj50FYem9NEa32KbNGjsyfLrMY0DZnCoFPK/OuvkM/vUEKVL6Hg855TlvC+Uw8HoDmd5LMXHc81bz+J0//XHPb055m59tsqoHqkSEph1iMqNcx7avx0YnRKE61AZLFYLJYpwkhb2V8P/B1gtofYCFwI/Ax4O3BIlf3mAa8br7cAJ1U7gBDiQ8CHAObPt+21R4ptbz9+jCRoevXcY4bM7xkVYWVNg7q3uWa3KiWbF1PC0JqqlZ3AwLkf6FxGuyubxWIZRXTsD1QIOsbzUth02CHMEeXrnVBF0tk/cTcSSEpCEJGLSIaD6m30frp0TQtIWiCKu5CK87aVi3hpZw9z2hqUaFHwIrGiLqWWCUA4kbCky898GXVh8wK1bzIeuaP0cj1ffd6D3fB7RmmdnjOo4whR7saa6uj313FUCLguHyx4UJtgQUs93ZkCF5xyGPgBvXsz7OjO8syWLrqzBeKuQ9vCFhb90zuJJ2IkvYCzO3pIPreNe5/bzpauftYUPc55bseQ00j15nDefxMJ1+Ev37SUS09eTEMqzt+cczTffuAFXr/hcljUyrlrv809z203dkzAr5+Gb94NO7thViNceQasPmrgQaSEHXvD58NcF+1Kq0lEIqbFYrFYLNOEYX/zCSEuADqklI9WrPoAsFYI8ShKNKr2p6Fq36Sq/mqWUt4kpVwhpVzR2to63LQsBra8bHwwu5GZrwfbbkxZ04B8eNGAUOzBmMxCjhU3LVOCqeQegsg1Yz6qbaNFk5KgpIUUo+ws7kaCixZIdHg1RCHW4e6lbbWQk/eikGyI3Eg1CeqXzmFecy2vbu1SwkQqDg1pdXPen1cOjmxRuT5M907MjbKVSmVrbpR9pMOrS3OT0ZxNN5DOUvKlEhk8P+p0Zo5dFrp9EFFykjnKnVWvBLmrvv8gF9+wXglHNQnqWxuYP7OWpXOaWLW8nSWzGskUPJ7cvIede/rxgoAlcxr5wJuW8tV3vYG/bUzz1sdeo7YvN6T5ZlfCZeWiVgp+wDc3PMuqa37B9776K9o+fxs+4H/iVnpuvpcLMnk2AT6QT8XxLvhX+MfbVX6RRP388s/Ly9R0KWJHDxs371ECkjvIbCTqc5PJq/e7rnqm0n5j3UMWi8VimQKMxDl0CnChEOI8IAU0CCFulVJeDpwNIIQ4HFV2VskWyh1F7cC2A5uyxTLxVHPdVLZ/rxYEXSkwjRejmX+kz200x3yL83YrclomP1ogmuw3ema7eY0pbEjDWqSFgSBQL8w/GZXcNo6R3yNASLWMaFFpfGEESBPmGekuYp4fiS0669cRtC2bQ+/LHXS+vofW+TOV6JOMqYd2IOVDF1JJECLKMNIiUGAEZBe8qPMURsZSqQSu4poIoK2h3GVkXgd9jaaDm8R1eOT7H1XPtcDi+SQba0hmC9Cbo7k2iRCCfNFnb6ZAb75IvugzqzFN66wGWh/dNGwpXuAImlyHP27qJFuf4sFjF/CbjVv4XFc/tXqjrn6S376fK1BfQAESuWL1AYs+nZ+/nbbP3w7AzLoku/vyLJvbxO+ufuvg750XwGu7YE6TciQNJiBZLBbLfnLTJdVuk8eZJRM9Abi4tn9Cj/8P/7F8Qo8PMJs/TPQUhmTYbzlSynVSynYp5ULgUuBuKeXlQog2ACGEA3wG+EaV3R8GlgghFgkhEuH+Px+12VtK2Bvr8WOoLmWV21UyZiHVI2S0coP0/EfzPKyDyDJlmMwuIu3ckYYrplQOph/CCJE2lsHAG3rTCRTosqsqykpJcJHRan08XeIlRLnjSB9PCOoXtdLalKZz827VUavgqWNpsckVkQvJ80v7AZFoFHfVAYteNB9Q++nxvLC7WWWGktm1rdTdzXBVQeQ4OhgdRIMhUE6qVFxl8MysgwUtJBa20Da7EccROI6gqz/PjLoU+Xufo/DB/0IO0wZem9kS2SICSPfmOPn+57nKFIZCkkTC0HC0GM9396lcome37WXuX38Xcea1vPLAi+U7FH26n9nK9Xc+qZxqYykMTXZR2WKxWCzTngP5E9hlQogXgOdQbqBvAQgh5goh7gSQUnrAx4D1qE5nP5RSPnNgU7YMhhWIxpf9FUYmWiA6ECpziyyWactkFoiCoLxkzBRldEmX6Y4RFeKOFo70OFroibuh+0KrI4bopJ+USspEJExJGZWkmSKMFoxiripfmtdMa1OaFzfvZtNz21UHs1xRHTOViPKICp5aXgydSDEjjyjuhu4PR52rF0RilC41K/phkHUYPFxyHoWPolmSpq+pVMf2g3JRKVAhygwRujwhDFZKOBoEAV3P7wBHMO/QVpYc2sqi1gby9z3HnNsfIbE3M6IcZ7dCiKylXNzZX3pqk2y79GQ2f/XdPHftO5EfewvBzDoC4NDr7lSlZ34AO3t4YP3TbNyyh0++7zRoqR927APGCkQWi8VimcSMNJAaACnlPcA94fMbgBuqbLMNOM94fSdw54FM0jJybED15KDSmWOWmVWWnE1UqdmBMhah1La8zGIZDYQSPVLxqPRKiz9arKkMnxaADMvKdFmaDpiuDKuubDFf5j4SkRijhSRdmmTmEmlHEDJy6TgCWupZkk5S3NvPn7d1cWTMUY4VNwzOLrmD/KjMK+aWB1/rlveOkTNUDB1B6J8ycg/5gSqTixut2cNSKjxf5R3pkGIpQSaMPKWiOr+aeCQmSRl1fJsIvIB77nyS009fBg2jnJ0D4Do0L58XXXMJLbMaabn6ByNuMT9Wl0YA9f156n/wEPzs0dA15kfH29lN/+dv54Ofv53vA0tmNfDCd69UWUteEJUrjiUP/cPkFpctFovFMm3ZJ3HIYrGUY4ojZeLObT0Io5vYKtk8qAikRZZhRaKKMYfrUDYWDJWjNFpYgcgyZTj5HyenE8ARkE5GGTxlJWRioPOlJBjp0q9wY/3cLMFyhMoCKvpRpzLTAaK7kekd9FhaSAoMwUYfzxSJUnGIu8Rr4hzpOmSyRdLJWPlci54ay3XVz0L4OhmLhCJJlKMkpHpA6JySUbmcmVmk3Ue6m5p+nk5Cbw66M6obWm8uCnAWQoUXl4Sy0JkUm9hsoqde383pRW/sDmAKaQL43TPQuw/t5YfAo/zLaUFAXO6HoJStnk1UC3wF+D7w4s4exNn/XFp3/jHzmddcyzc/d7F6fw+mDnUWi8VisQzDNEhWnH7YG+sJ5rYexFUdiC2euifZ4vEpHmWV3AwMnkVkOosGbFNlTHFVB9zWM+anMxS2tMwy7ZmMDgCd+1PqMmaUl4HhcCEURAxRxGzZrkUlV0TjaFFFiwOljCGjHE0/1+KP7vxlHivmGIHR4T66jbwbtlef3Ui6LklXZ58q20qE28fDsGon/AqjXUTZAvSHLh99bIz56uM7xnXQwpAfRE4g3V1Nz6XgqYsVj0FjDTTXgpRsf6kjdJsYIlreg7395Z3QxhsBP3l4E6+/sHP8jvmNu0dlGAG4QCcqs7yrJs77h7iM+1s9t1CA/N5HypZd/+438t8fOZP/vOdZ3DOvYcW7v84la7/DX376B1z6yf/m6i/9XOUSjQaTUVS2WCwWy7THOocOUmx52fijy8XENbsR2fKvqyl8rmAjG5hftm0lleVmq+RmrmAjbQzMcBBZCdfsRo6ze8ic31hh3UOWKcVkcxAJEYZHYwRKV5SG6dIyXeIloJQkrUumBBg97ssxBRZdtlYt5NpcpucWVneV3Em6jEcLNKCEo5gLc+I0N6VDp5BxDuFUS0HUOkfICyCXUfOqSymRSYtXjlClc2ZZnXY9aXdUELqM8l7kJNLHzxWVeyihvjr1ZArc9YvHee+ph6tOV0GA3N6NqE2E7iF3gv4EJ7j1w2fQPqMy2nkM2dk9akMJoB9oA8gWuWyYbfeLWY2wqBX54OcGrAoe/Bz05Xn+0VfZ0Z1BCMGLO7v5q1vu5bpfPcljX3gbR5++jId/uxEpIVvwePNR7Zz82Z9w0fELWfepc0Y2B1teZrFYLJZJhnUOWSyjyOq5x6iyryq0khmwbTWBRS9bJTdzFQ8zq4owVGKQY1Xlth7Eyk2IuS8iVm46INfReDiGrLhpmVJMppu8kmATvvZlFEQNkVhU5qCpKA0z/9Hx5eAWjZIriXLnkRaGhLGs1EZev3bKy9J0oLTrRHk/Uqqsn5q42l5nF2nXkxcKSnEXapMqnyidVAJOvqgEHS04mQKVdjOZ7iI9B9NhpV1EUkJtAuY2KbGoqZalh7XythWLVAevuAvxGGJesxIeYu6EtkQ/ZGYtIhUfvwPOatznXYZy/MwPf16G6nYy6l9Wrzxj6PV1SZa+eSlvvvA43vTWY7nir96MvPf/cOuVZ3DR9b9GnPpFTvz87SRiDv15j4/cfC9Pvb6bI+c179s8JpOobLFYLJZpj3UOHcRY99DE0CHTzKoQggBEewy2Dty+0kXkeEWkAAAgAElEQVSkn69rfhjRNYxhft4I/xfWZWna0bTFg6s61Jfz/XAejVfHNesgskwpJpODyOwSJoxl2tajS8TMdu1SKkdPWTv68Il+rkulTOFDb2tScg5BWYlXab2I9nNLBwuFG1TZmS9DgSiIxKBUPCoDKwVcO5EApDORErGwM5mnysJcNyot09vq7XVpmZ6360SB0ro8zPOgQFQKJ4CGNLW1SWhIVziwJpiufvZmCjQdNmv8jnnlGXDtHUqMGyGiJj5oLtBmYFZDDV/pyZIcpSmWWLMCVh+17/slXC59zylc+p5Tos+LIzgWOO/tK/n6qE7SYrFYLJbxxzqHDnLsjfX4cwvLyeGWLZM1QgVID0JVkWUYYUgC9AcjcgBVK3UTWakCrveH0IX0G37MrfLOUp6SxTIcQghXCPG4EOKO8PUxQogHhRBPCyF+IYSYmDrJ0WIyOIhMh44u1dJlUzqnJwgiRxCUO39KJWYMzBEyRR1fRuKJKYro9ZVlbILoQcX2WjwSxmudPVQKmdYZRiJsX++ojXVeUcFTcyp46vxijmprH3PV+lwhEi/MsXVpGkTCU9FTrqSCp17HQsGo6KncmaKv9tEdzPQ5TlTOkEYCe/tpmlE7vs6l1UfB1RfA7Eb1/s1uVCLM7CEcRbkiXHjcgMV+zOGG+hQ7e7IlB9GosWYFfPq84bcbDtN5d6BMFkHZYrFYLNMe6xyaBlgH0fiyQcwHCVewkVYydJKm9bp65dD5WPm2lRlDI2lnX4rzAOgKRuYAGqz8bF/K0jQVLqRZZPgUj4IMz32Use6hg45PAM8C+gN7M3CVlPJeIcQHgE8Dn52oyY0KE+0gEoT5PEK1cy9lEBnumlKIsu7ohSH8GAKQMAaVUMoJgkh80GHM5gTM/U13kp6fHrfUrSxEVjibJOC4oXsnPI7rlHdY09lD2iUlws5hXhC1lXdEKIZJVW6my8eSqjuacjA5UQt7LawVQpEoFY/W685kuUJ4HRwVVh13jXI6RrfT1SDRTwPY3Ue+4JE8tG30jj1SVh9V3ZFzyQ2wY2AmUVcqzusnHsbRC1vo+MbdtBR8etMJtp21nA+9YTHvKXj0XPsLmgZxFw1LTUK9J71ZaGuED6+Cc4/ev7EsFovFYpkGWHHIYhkDNoj5pfBpgPVrqrd8N9vXVwpDPSRopDBgnwF/dB9JMPW8mColq7Z8HxlJ4PZoYwWigwMhRDtwPvBl4FPh4qXAfeHzu4D1THVxCCZWIKrmXHEd5RWu7CQmnEhgkaFQhAyFJVHuitFCjhY+ysSeCgEokOXCjxZ7zG30a8fcNxRlvCDKJYKoy5h29ghUqZgjoOAr4UeLSUEQOTscwyCt5xkPS9M8X7lXEm4UgG2Wv4kw78gUxTypQqdb6iIxyg/AK0BeKDeSzk4aqaAzEoKgvEQQyrOl/EB1dOvqJzlvxuRqwf6+05D/9itEwS8tkgmXPx6/kHM/8+PSsp99cjVnHjGPZfOaQ1HOh/e/ieCmu3H8ah/qiCIQb6hRQtCsRlXmtvqoSETU76v+bJmlg5MBG05tsVgslkmALSubJtgb64lnKFdQNeGo/sZmChX/iw729VhWE37M9etmImvKvwUPV+o2KCMM3LZYqnA98HeEvapCNgIXhs/fDhwy3pMaMybDzZ52smihQgsmenkpm0hEN8s6J0gaKo7p9EGWO3cq27aXStTC13obXYZmLqswHJWW6w5mQRCND2q5FnL0TX8qrnKG4q5y9QiUAKDDqr1Aba+FGy9Q42pHUcGDbEH91JlFrhOVksWcqLQtFVeh0/GYEqzyXpSHJKXqZra7Tz16c6F7yTjH/S450+4vGc4/vPadPbB9r3LmdPbCwhZorNnfg4wNb1rKPW9YQl9dEgnsScV5d8Hn3N+/yL9cejIAxy9s4cI3LaV28Sz1PjXUQGsDvPVYnCvPpFgTH/TSecD7AH76Sfjt1fCjjythSFKeLyWler8CqUSk/H46ksYKW15msVgslgnGikPTCCsQTQyVYdODiUQDupetaSB2YyuyPaa6L7fH6KF69xnRPowDaE0D8rq2srHkdW37FUY9mNuok/S+j7UP2NLIqY0Q4gKgQ0r5aMWqDwBrhRCPAvVQxS6n9v+QEOIRIcQjnZ2dYzzbUWQiBCJTB5bGLbUZuqwzfZywxsuvUnam99GlXmXrqtyqV4pFgSk6CSWmlEKtZeQC8oPy7RMuZblJA9wy4TzMLCJQopfOAaqJR8cqekoU8ENBSIs/ep0umSt4yklUCMvKtAhVNpfwfBKqOxn1KXW8RBiC3ZhWXdM8XwlF27pgx17oyUHR57PX3sGMM68t7w4HUUc5P7wmeU/Nw8xA6s9DRw/s6VM/X9+jtp9ZD7ObYFFrlJ00mZhRy6ovv526717Jq9ddxsxcke+Hq9b96I8AvOsNi/nrr9/Nlic3h5+B8D1KJ+Gtx7H3Xy7jhkUt5CqCzfuB9wI9x8xn02Ovwp5+9T4W/DB7Kvzs6g51+jPTmFai4mTDCkQWi8VimUBsWZnFMsaYgk81YWiV3FyWT3TtnOVsEPPVtmsaysrF6m/rIbd2FykMe37oAKrmPirLNPrYaSPKNBoOuW4mmJ3Pwjnckl1+wGMPhy0vm9KcAlwohDgPSAENQohbpZSXA2cDCCEOR5WdDUBKeRNwE8CKFSv2238xIYx3iZnWciSRQwj9nIrMHyKBSAsRmlKr+aBibGFk6xg367LK22KKSyVRKjye3tzzozn5MKDURwtEWhiqPCfd2azkkBJRiZjnQ19OOYNm1BkCmAPJUCDSTiLT1eT5UelRzI3mWQzFrGQsum6C8vK1mKuEB89XIkV/XglErsPxC1voyhT49rfv532XnhwJYjoXqjenxKV46IzK+dF1jLkwo1bNaU+fCnuejAJHNVwBsxtZNKOutOjHKxZx6pObafV95PqnuXd5O2f90x2sPLSN735ytTq/UOxqfbmDy7ftJSElHuovm5uB/w3c3VDDzic388snN/PTT6zmojcshoaUCiN3HfU+mJ+RmFP+ubVYLBaLxQJY59C0w95YTzymiLNKbuZTPMosMjhE4c6r5ObqGUUfO42vcgI7SauO08M4gLQbSY81VOv5EbelH8SFdPX2rpIzajREqMGwDqKpiZRynZSyXUq5ELgUuFtKebkQog1ACOEAnwG+MYHTPHioLGUqiTQMNP1oYcYV6mZcl3T5oWhiCkqBjMSisnweopwfUe2h52E4hbTTKOYa3ceqnEdl+DTGsfTxtIijW9N7fjRHvU7nGJlja6dSKq7cKlKqTmS6LK3oq/Ijz1dj1ITB1LmiWidlVL4mpRKhujPqkfeU46cuCfU1EHM5bkELb195KH+x8lA1B+04SsRUOLaU0J2NHC76HFPxSAja06cEpOQUEYZMEi5y/aeRnzqHS57czKyijwO4e/o5/cGXeOSco1kyu5FVn/sJf/jNRr730f8Hq/8Z/u3XtOQ9HNRfNQtCMOeKN3Peh1exsycLwGXAMTesR156I8V3fR1+/lhUPpYPy6F1MPlkxrqHLBaLxTJBWOfQNMR2L5sY1m8bGEp9BRvLXECgwp3/nocHzaYww67XP1xesgaDizzDCTb71DmtwtE0EsrGu62HzrW9JbfULSwfk05nlinBZUKIteHz24BvTeRkxoxxdw/pUisxeLA0RAKPaaRwBEhDqAmMEN9SuZl2ExnlZvplNecPRIHTQRAGQceiEjXt+DHnao6rO4BBuYOoNF8ioUl3FdOunPqUUa4VOkd0R7K4G4lPOl8oCAWigqeW6238QO2TjKmH1CVfPjiC/NY9ZPI+ze3NkaAWc1VQdCYPMZe+XJGrzjuG2ubaSPzR4pkQ0NbAH3+zkZPmNEYuJX38bEG5oBrTqpRtqppfGmrgew/iFMt/9zlFn7pfPkHmxMM4Zv5Mvvalnw36j0FKSvjhQ1z+w49z+ZlH8sA37ua4u/9Mbbg+vjeD97W7iBV9lT/Um1PvVyIGG56l48bf0lL0cWY1wkfOqN5lbSKxAdUWi8VimQCsOGSxjBGVYks1BgtxdpGqPfxtLWWuIFMAqiYGjbVrZ6QMWUp3Ww/iqg5mheqXdkshGbFAZMvLpjZSynuAe8LnNwA3TOR8xo3xFIh06ZVvlC1BuWhkRAkNCIY2y7hch1KGeEmkMYQcXZpldtOq1i2rVMomopItjdlBzZxPqXxMn5MxRmlf4xh6DD3nYqCygWLh+AUP+vLRmF4salOvc490yHXcjUrDCp4qU4qJyJUUC3OHHAeKHsnaFMmGcJyEqw4gQ0dRKsGujh5mNdbQWp9S4lLBi0Qe7WISgpPOOVodQweFF32CrXvYvLufhScsVPtOdXYObG0PwK5entm6h02dvdwBJIcaozcP9zwHbzmSU59+fcDqmBfQ9+376V4ym3kLZqr39K6N8B+/o00LUzu7Cb78c2Wjn2wCkcVisVgs44wtK5um2BvrsaeaQFIplAwV4pzCR1yzu+q6sRaA9nX8EZekAeKa3WV5RaDO9Qo27tMxrfvNMiUZDzeAFmw0gS4PC5fpzmEQOm6IBCFT3BFGNpDuPlYSk2RUVlVqGS+j8cwOZpUPIZSDA+OYpoAlRPl8ygSscBBd8mZSOQZEeUE6kDgRU4JMINVrz1edxXqyUWczLXrpYOR0QgkLBU9tlyuq7fJhu3VHqBKv1nqoTajXXhDlFqUTkC/S0pCiPhWPxpdE5WlBuL0+L13a1pOFjh6c5loWLm8/OIQhUO3mB1n+i4tPYGPRZ8FIxvnWfeq93dVbdXVtX57fPP06X771D2qbW/+g3kcDp+jDN+7et/mPB7a8zGKxWCzjjBWHpjFWIBp7KjN/oNzdcwvLyTFEd5ktHmLlJsTcF+mcsxVu6xmwyXjk/AyHLpkbkYtpqzdwGYO7qCyWg47xEIiE8VMLJb7ZVl1G5VO6RbpuNV9qHa8f4UBa2NHdn1wRCUFlB9W7SsqyjsxNzM5mpWMHhpBlHLpMwDKOL0QkpmghSo8N5Zk9erkwMoZ0p7SaMLvHCc+p6Jeflw641rk/gVSijs6yyRVV6VjeK++GJtRYsqO3VPLmB8pFVHo/ip4SK4p+NB8tLPXl1Dxa6lQpWX1qpO/+5OfKMwaEaftxF05eDNfeAR09I6ua6+xR17m1vupq0dbA8YtaeXFnNyd//nZkx8DfocBAJ9P6p+Hi6+GNX4BLblCvJwIrEFksFotlHLHi0DTHCkTjR6W7Zv22J9kg5vNVTkAO8S1YbPEQUpVfFdbu4to5zWMynwMdyxSCBi2pu61n0H91hnJRDYZ1D1mmLGMpEFXrxKSdM9qlo5/rIOiSS8co6zIFGjNIumwFUXevkjupIhNIyoECjs4X0oKMG5aCORXzqRSzysKpZST8lAVdE21nOpNiTlTOlk4oESdXhK5+Hn9is+oqFnOjNupa4CnN04mCo+OuOl5fTgVEZ/JhcHXoZkrF1fiZAsIV0FQLqTi1NeHyREyNoYUiXZqnS9YQao7pRHkXtYOF1UfB1RdAWwMSKDTWsPsvToQ/vKjekxFSaKjhid+/wI0NNfRXrkzG4F1voD4Z59sPvMDZy9vprxukUK21Ifrs/OopuOYXsLNHvd7RrQSriRKILBaLxWIZJw4Sf7JltBgu1NhyYKyeewzrv3a/Khfb6vGbeZtUa/hPQzXjTOUtXgqfde2PIVl0QHPY13UjDbPWDqJq4ds6a0j4A/fP4XILy4ee+CDY/CHLlGWsMojMrmSV7earbUPFNqU8IhFl8AgHHMNJY/4MpFqnj+cb27mmm0g7kkR0jMHmpfcv7V6RmVQaS0YCjh++1m3gpaEuSamyhmrC9uZumCvUWg+ZAsctmAl7M0pASsQiJ1NNIixNM9xNMScSv2KhSJQtKPdPrqj2iQeQKRDkCjjNtep4GU8dU7dSN91TQRAFYvuhwKRDwCEqP9MlabVDpvFMDVYfBWctR+SLJLqzbLt5A61DOIa0eUzTD3ywO8v3P/sTAGqPmMe7t3UR785ASz2882S2LGhhS2cvN1/xZo45ZCY7HMGiO57ANcKwZcJl2+nLmNfZq1xk//arqLuZJldUpWcTkUtkw6ktFovFMk5YccjCXcGPCHYcXrZsJGHKln1nldysBBKdubPFI79219Chm5UMUpY1UqoKN0TvczWBx1w3EqptXy1rCMBH8FVOOKBuZVYgskxZxkIgMsOmtYvGdPFUVoFpUUdjhk6XOomFr6UhyGCOY5RvadFGogQV00WkHT363tw1QoX8QL02ty07BuXClGOcjwBcVHi02QFMl6A5QF2KUnc0fU1irhJaAqlcQDlPCQPJeBgKXYjKucxr54TOHj3vdEIJCP15VSrWn2fP3gwz6pIlo5XMFhG1ifC4jjpOaR6i/HwdUd4RzQvU+RQD6C9EDqcR1V5NYtzwOt7/PIff+dSQp2N+XHcBnwCemtfMhveeyulHtav3NxG6srr6eemVDhbPaaL98Nm8yfOh6OMtmEmuKU3ND/+I6C8AkEUwr62BTGcv6Y1bkL356vMYLER7PDhIBCIhhAs8AmyVUl4ghDgG+AZQB7wKvFtKOaD2TwhxDqpxgQvcLKW8dvxmbbFYLNMHW1ZmAcCZ/cKg6yZD96uDhSvYWDWMeZ+YN1DTrcz7GQmDCYBmiVjle19ZIqbzkMTKTVXzkMoYRNQSSNvG3jK9Ge2bvlJ2UBC5U3Q2jxaCdHt6cxutfJREFT2gLP9ZGqPCSWR2QtNuHi32hDfnSKIyMkE0T2S0bemh5yIjYSaQVXKGdAlZmCfk+WGZm94n3L6UUSSjjmZ6P0dQLHh43Rm1b29WnUsQQHdGOYP8UKDRXco0Xnhdk3HVpt11kUWfZNwlr7ODdveRLfoE/XklIgXSyEOSStDQoddmGZnjqGPFY+G5CiVmZQqqPXumEB1/KvONu0dcWCxQrqFDzjuGs46cR03cxdPZTz1Z8ps62LSpk8XzZ0L7DGhKQzoJrQ3ElsyidnYTjheUPl7pgkfhP35L+onX8L/zwOAC1WAh2pZ94RPAs8brm4GrpZRHAbejfNRlhILSjcC5wBHAZUKII8ZhrhaLxTLtsOKQpcRQAoMViEaHAw1dljVClaFVsK/lgCN1hg3mMLp2TjOFtbtKeUhii4e4qgNu6ykLoi77TFURtWD/soaqYfOHLJZKwttfM5RZL3accpeOLgfz/Kh8SgdYVwZV64cOUNYCkynclNxIIhJUXCcSiXwtyjiR06k0dpV8IlMscgRl4pUpGnlhiReEr03xXUa5RnpfffzaJPFZjWQLHs+/3Am5It7uPuVEkqgxvaD8eumyOy2u6bnli4hkjNrmNMmmWoi7vL6rj1zBoydbpHtnj+qcpUUgP1CZRb4h2Jm5UDIUzhyhSuYSblR+lg/dSn05yOrOZ1NQKdpHV8584KGXO3jwpZ186/7n+eqvnuLOh17m1df3sH1vlkWbd8FffxdO+yK87f/CXRvZs3GLur7/de+AsrFEINn69d8hugYkF0VcecZ+nNgoMsXDqYUQ7cD5KEFIsxS4L3x+F/C2KrueCLwkpXxFSlkAfgBcNJZztVgslumKLSuzlLgr+JG9wR5jOkkzq5pA1CyQXbLqXyylCwTAvJgShtY0VB17LNrPD1aCdgUbBzieRFbCNbuR4fzM8jSAa7Ycz7qah8ucU7JGcEt2/7KGqmHLyyxTltEqL6ssERuw3qg3E9HTUrmVWS6mx9LCkZmTUyr9Crc14n1Kbh+M/QmPEXeNkGvDrVSam4zm4lTMxfwH0ixZM1vYAzSko3HzXiSOSUAG5dtql1QYUl2fTvDQI5u4/dFN1CbjfHzNCjVnJ6aEFy106UmEwhKBVK4fHaaciIV5RQ70ZGlIJwgCSU3CxfMlhWyRWKEHJx6LunblPHWsQIadz8KytyBQ45nd37RQpMU0PxTH8uH84m6UlWSKYZOVWY0q+HmEbAbue347AH96pZO3HreAH/zxJR5/bTfr2+pZ0NEbnfKObvjXO5lRk0T25YDql2NuIHkNWFjtgI01E5M3VMnULi+7Hvg7wGwttxG4EPgZ8HbgkCr7zQNeN15vAU4aozlaLBbLtMY6hyxlDHVjPdHt0g8GbmE5sqb8a6msEVzTtZKfcWjpPqu0DqDBQX5tFvLhRYMKQ4MxUgGo2n5D5QwN6oCqKB0z9716exfyujZkewwpUD+va+Pq7V3DztFimRaMxk2fme8jKh5mS/eyNvJGpy/zYQo7pfENIaZSSDKdPBLl0ih6kZul5EYyxg6CqAxMr9OOn6KvnDu6zbtnPPS2pS5fQXQcz4+WuU7Ubl6XpvlBxfmEz70AYi5rVizCEYL1T7+O7MmGolh4/fK69bwXzTlXVI/enFpXn4ra2Rd9Ct1ZGtMJmuuSJGMu6USMRENKCUNSKteQHjOTVy4lXfYHSujJF8OSvLDMTDuxtLPIC6LzK3ph2VlWzak/r66jdoANhV/5W2icuPIMdc1GgEzE+N/h83+59GTkzz7Jz699B499+e1856h2zjKFIY0voS8XVSpWG1fAHTCw61kqDn9zzghPxFINIcQFQIeU8tGKVR8A1gohHkWJRoVqu1dZVvWDLIT4kBDiESHEI0XyBzRni8VimY5YccgyAGf2C0OKClYg2n82iPlck13JTtJlAskGMZ8bxfFcy4l0Ey/7A7zoChBXdexXC/vB3qvK7KDfzPspv+HH3CrvZJXcXH07Y0zRXv1LfIccpkRsTQPy4UXIbUvKxK7R/ExZ95tlSnOgApH5j4es8qhsVS8IS60Mwaeyu5mZWRQE1V+brhbPV2KHH0QullJIdMXYOvtHrzedT3pcLfwUDRFJi1n6uV5e9Acuc4QaI1OIxKpcUc2xYIhX4dx7skV+8cRrHH3IDCRSiSz5YpQ5BEaZXHicfFFtF0glxOix80USzWm1X10KapN09ubAdVU+UX0qamefD8OwtRBW8FTeUX9YctafjwKvzfyjZFw5i8zSOn38Qig49WTVWH2hWJQtRscKpLoGu/voemoLdB1Y+fN+sfoo+PsLYFbDkC6nQIC48gx2HTkPgJ8+ukl9vtJJWNjKezp69vuLrSNhLWD4zpRj6OoLJodrSDM1y8tOAS4UQryKKgs7Qwhxq5TyOSnl2VLKE4DvAy9X2XcL5Y6idmBbtYNIKW+SUq6QUq6I71urD4vFYrFgy8os+4lteb//bBDz2cD8ckHkY9G6K+RGGimW7SOyknXtj3H23IuBkYkp5ntTmTFUKvn62v0U1u4qlYjNIsOneBQkZSHRle+3XDcTzK5rKAdU63X1+9zZbLBjHAi2vMwypTmQErNM6BJBqjIlXablOlEgs3YX6WWeHzloHCdysyDCrBvHCI82SspiTnmZk87+cYUSIHQotc710ccteFHmkR5biKiUSqsDWrgqhE6YSmGmNGcRrhflYdjmTwGyO4tIxiKxCiJhyHAQ7ezOkIg5FP2AP2/dS1tDDXWpOHFX4Ib7Obq0TEr29OTY3Z8n5gjqaxIU/YCGVJxEzKHgBXRnCxT9gMa+PNmiT3++SLqzF4DefJGd3Rn6ckWEEGQLHpmCR67o44TXWUpJ0Q+QEpJxl5gjSMRcMuG2xdDt4wCOI5ASYq7AEYJ0Isa8GbW0NdQwr7mWWF3oaoq55flNGeWyaF46OypzG2/OPVo9ANY/DdfeEZXpASRiPLrqf5Fsa6C5Vt34//7FnXg7e4ilk0og6ximMcIwDNCl+q37ZDSQUq4D1gEIIU4HrpJSXi6EaJNSdgghHOAzqM5llTwMLBFCLAK2ApcC7xqfmVssFsv0wopDlqpUa29fjcEyaSz7z1AlW/p6m1k+1RjW+RW6hVjrkapYn8LnCjaygfmskptVJ7KtHsyLsf5rPUCDcgABXLO7tO6aLcdz9Zou1q+ZHM4yKxBZpjT7KxAJgezooT/v4UuJlBLXEWFVlywJCb6UJYeFEIKiHxB3HRxHlJZ7gSQRc0qCSK7ok4q75IpKTA7CEqW46+AFMjQGqf316yCQJGIucddBIvF8iRcECCHU3IQgQAkbQgi8IMANRZEAJYzszRSoTcYQ4fK465CMOfiBDOOQJDFH0B+GDOt5ua4ovRZCIATkCr5qahZzkVKSKah9Yo6DELCrL8ez2/bS2ZPjt89sZUZdijcsnsWshhoKfoDvB7iuE56Xg5Tw0s5uurMFmmuTNKWTxF2HVNwtiUO7w5ybmXUpenNFin5Ab7aojteboz/vEXcdWupTOMk4cdeh6OfIF32yRR8HNd+6VJzZjTV0Zwu8uKOb57bvZU5TGkcIVcXnBxT9gHzRpz/vsTeT5/U9/aVsHoC3rVzEZScv5tyj55NuTkfiUMyB/jzPPredZUtnQ0t9JCBNBNqp8427VVh1WyO85xTmts/g+vVP0dZQU9r0737wEF99/2kwp0ltN5ot571AzWEyOYdgqmcPmVwmhFgbPr8N+BaAEGIuqmX9eVJKTwjxMWA9qpX9f0kpn5mY6VosFsvBjRWHLIPizH5hxAKRxgpF+4cpsg0aWh12+xqJI6dStCtz89zWg6hw/VTSJjKs/9r9arst4XZbPLiqQ4lCa0KBaE2DGnsrIOBqovyg4QSskc7dYpm27I9AlIghmtLUZQvg+XieElu8IAgrtCRBIKmrTYCUBH5AJu+TLXq0NqmyUOkFpe2VUKFEh5p4DCEgFVeuGd8Qh2bUxFVH+UDiS1kSopJxF0eobf1A4ghJAoe46xBIJbDodX4gSQqXdMIlU/DD5mYSL5DUxF2EEMQcQbbok8l7OI7ADUWmfDFQYpOrxpZS4jgCL1BiVNH3ibsONckYvh9Q8FSr+Zp4LFyvxLGZtSnmNqV509I5PLN1D4++2kkm77H6qHZqk3EyBQ8/FLz8QLJtbz+7+nKkEzFqQ2EnU/DY3ZfDdQQ79mYQQpAMRbVMwaPoBc5P1ncAACAASURBVOzN5MkXfbaHriEtTmnRam+mgOcrEa07W2Dzrl7+vK2LnmyRPaGbZWZdkt19wztbUnGXS4o+XwHmP7yJrY9s4n/XJfljWyOXnryYd71hMa0zVFe1lvoUfZ191DmOcljFXeUkMksPx4NfP6VEmY4eJfh8aBUUPeZ96af8064+9qbi7EbVIXVnC/z28ddoen4HK9asgK//bnTnMppi02gyRQUiKeU9wD3h8xuAG6pssw04z3h9J3Dn+MzQYrFYpi9WHLIMyb6WCI1GadB0FAfM872F5XyKR8u6gQ3Wwn64saqtu1XeyazqWY4lOmSatmt2DxCQRFbSsbaXViMYe/22J7l2TjNXsBHmZBDtsRHNtdpna7Tfd+seskx59lUgijnK9YEKJ46FIcXx8L4+YWb6CIEjoE5CHZRKvwQQ19sJQUKE++kwYzfM8Alk1AlM6H0FsfDfl0Qp/0i5iUppahXh1Y6UxM3lQpAO55iUUCulEifC+dSW8pSMMjezI5p+7vlRmZy5HlSpmi6fCsrXHxpIzjpuAd+79zk27+5jxaJWjjpkJo1zm9Q8/AB6c3zv/ufZtjfDecfMZ/6MOuqa06pDWX9eZRDFXfZ2ZwEo+gFN6QR+oNxTvpRKIAsk2aJHfSoeTkNSm4pH18d1wlyhGPiS7o4etnX1s2xhC8++uotlh8xQWUY9WV7f1VcqNXMd5Sba1tVP/eOvcuK9z5Hw1TkeIuHLvXnWBd28uquXU774U7qzBU5bOodjDpnJu9+wRJ2L50N/eH0TYfezVFzNRZfsjQXrnyb3jz+NXK07u+FLPy1lPDnAjFyR/wTam9J0BpKzTjqMF1/u4M/zZ7KMUTY8tTWO5mgWi8VisUxarDhkGZL9bW9fKSYNd9M/3UKuhxJENoj5IFW7+FYMsWUfO5UNxqBlayE6O0iu3Vn1C3YrGa6d08y69sdUSVmTw7p4gNDREKHDaP3X7mf13NMGfW9NIXEsBUErEFmmPPsqEOlSISlBDBLPqztb6f/JS0HVMmr5LlGB0zFX131FIopusS4pzw1Sg5Ufp4SRJVSWuh8uKM1HRMfSeUZ6jgLVeUqP6zrlY5vPdS6RORfzfIUoL53SY+aKUJfi/GMX8MAL23lxZzeN6QQr65J4RZ9bf/8iH//uA/SFZWznHn0IPbkCiV6HRDKmhJS4C56PEII9fTla6lPkiwEdPVm8IKCxJkEsFaeuLknXzjzb9mYoegErFrVA6MIiEYO6pAqcLnh0PqsyeJcdO58XnnqdZYvbVKh1X55C0WdOc5pYzKFQUA6lbNFnZm2S2f91LzG/XOivBT7Vn2fR+qe56LgFHDGvmdpknM/85GE+d9sjvOPEQ1k2t5l3nnQYyw6fra55IQzEjoUd02ri6qdrZEqNBl/91YBy51L3u4pz+OjeDIm3nwg1CZYsmMmuXX0jOoSH+kjvBmakEzgZ1SSr8hQ8AbEPnLbPpzBu6H8XpqCDyGKxWCyTjxGLQ0IIF3gE2CqlvEAIcSwqOC6F+j37USnln6rs9yrQC/iAJ6VcMRoTt4wfI80fGorBxKJ9FZEOBkbirtKh1QDrH973YOehxh6sbE0ChELU6o+dxq3cWXW7XhJ8ikcRW0JnU9fA1sfaYbSKzXTO2UobGQIEDrJ0DHRJGtPjfbdYxg0tvpjCyIAbdy3IyMgRVNmmXqDEpcp9S44hp7w9vV4GlHWaNsc3kUa4tSlWGZnUA8QrCEOrw42HEiRKbiCM6yGV0OEI5R6KueF2xkFTqvtXY8Ll7ITLgy918PArnfz6qdf53oMv8cKOvRy3oIVHX93F6qPa2dTZy5zGNIm6pOoCloyhxa6auEsi5vJKZy/pRIzXdvXSlyvS1Z9nZn2KZXObeeK1XcxpSnPq0jnQUqfmGXeV8CKlCtvu6MYLJHMWzKR3826WzGpUwlBBdUpLtNSpeWeLJLIFEgmXdMFXYtXe6n8QWACsu+BYvvv7F7n/hR2lcrUPr1pGIuayYGYdH7j5Hp7YvJu3HruAY+bP5NyjD+F4LRZ5fvT5ScSUiBVzosDwaoKRFhwlZY6z0jopoSc3xJs68BzEoW3qeM21tMRdghm1iD0DGtKX4aLEobOXt7P+cxezp7OX3H3PM/uOxxHhdfBq4nyhNsUXzzt2xPOxWCyWavgbn5/oKdD9ljdO9BRY/MMPT+jxD79r94QeH6DvkpMmegqK235cdfG+OIc+ATwLaPvCPwP/KKX8lRDivPD16YPsu0pKuWsfjmWZZIw0f2ikTEdRCCZHl7eqZWsAzaIkDA26XY2ArCxbNhhtYeczva2rbxaN7KLxKiG07iHLlGdf3ENalzFdO9UqSXULc6hyM286bozxtBCkBSJpjBVI9a1C3/RrXBEJMKZYIA0LkeOUlZUNFHaM/YRggOpgCgwaRxjCVLjOFCUKqLbwiVj5fINwzukk8QUtvKmlnjflF0AyxmfXrGBvb46tXf1cd+eTNNcmqU3GyPuBGjceo3dnNzWJGLGGGhLpBPOSMeY2pUvB2gCNNQkCKalNxujPq8whmmvVXOLhV7Ni2Iq+L09fpsic5jT9O3sQCERbfVjelo1ygfrz0JdX1y4RI5EKnT2tDVW7eInWer7yzpM49fA5vLBjL7v6cmzryvDNDc+Wtrnh8jfyf99zCulEjNX/8ks+85OHOap9BmtWLOIty9s55Yi5UF+jPkeFrJqT7kTnOpBKQH9OubFCAYfeHMwIz9X87PgBZPID3u6hEK31kE5En5nGNM7Ji5F3PjnkGL21SejPs/asI6GtgRnzmuHY+fDxs6Arwzu/+FPee+rhfPGERUrwmuxM0fwhi8VisUwuRiQOCSHagfOBLwOfChdLIqGoEdg26rOzTCr2N2B4JEyXUOt9PbexuN66bG0tT9BAIfp7eZeMysE+dlppO10+1iHT3JJdztUMMAhWx4WUX11EElkJ1+xGjlKp3EiwApFlyjNSgUgLK7L0nwpM4cYxRBRpbF5FfBEChFtengWRcyQIIveQdgKVjWk4hbRjRJcKxRxl5dD5PyVHieEs0suqndJgy4QwxiQStAQqH8gPlLCRiEXn5IjICZWMKydMb1a9bkrTlIzjhu3k73xyMzHXobk2SSrm0txaTyruki14xLv6CSSk0wmEK+joybKwtQ7XcYjNrFOijuNw7Iw6yBaUyAFqTkE4r6IPno9Esr0rQ10yRv3sRnBd1X5eouaY95TDyBHQEHYh8wNIJ+EjZ8I/VbSFT8Vh7VmweBbnHdrGed0ZyBTo2N7NR888khd3dvOrpzbziVv/UNrltKWz+ff3nEpTOsEjmzr5+u+e4b3fvJudPVnef+pS1r31OObNC3OZYmE2U09WiVz6c7KnX13bO5+Ebz8AnT0qI+s9p8CpS3nt1U7akjFqwpK9oZBxF3Hpyeq80ypgnaKP/OPLQ4tLiRhXhu6gC88/Vr0PGseB5jTf+fAq/vDiTvr789T6MnSrWSwWi8VycCOkrPaNqmIjIX4MXAPUA1eFZWXLUG0lBeor3RullK9V2XcT0IX6CvNNKeVNwx1vxYoV8pFHHtmnE7GMD6PpHhoJB7NQNFIORCAa6vqpYOoq5WXtMc7eenHV40eB1sPnFsls1Cq76jYC5LYlw85ztJnqApEQ4tHpXp477X9HDCcQ6fKckjg0VHmZsc9Qng1RsV2pHIxyl48XhE4LUT2wuEz4CbcveJEQY34nMUvOBsxnuJIyIsELQ1TSriftBpFSCSuJWDRffY5mOZwfls/lirA3w2s7u/nKzx8n7/msXNTKzPoUKxa1srh9BhQ8ckWfrkyBjp4sS2Y1km6rj9wzybgSJHQJlp6HRAkpMVddk1wRij753qwSlPT86lNKxOjLqedCKHEp5qocINeJBLGGGvV6/dNRW/hZjXDlGYO3Z9clYz1ZXnxuO//9hxfZ0Z3hpZ09/O7PWwFYtWwuX3vvqRzR3szzW/fy8xvv4l1bu5gTSHa4Dg8dtwD5hsXMbkxTm4yRiLnMaUyTTrq4D76M+50HEIVIANLvcLY+xZYFM1mycWvVt1eGDiG/uZb+C46l9/iF1NckiLsOz23r4qGXO/jIdx4Y/KPRWk/3OUfT9N3f84mzl3P95y6OcqvyngoST4bXcFcvP7/3OS686PgovHyyE7qH7O8JaBAz5EnizImehsUyadjxNxNfVta/YGAUxXhy+M1dw280xvQtaZroKQDw+9s+XfX3xLDOISHEBUCHlPJRIcTpxqqPAH8jpfyJEOIdwC3AWVWGOEVKuU0I0QbcJYR4Tkp5X5XjfAj4EMD8+fNHdFKW8We0y8uGo5owMt0Eo33tGGcyVNlWm8hU/Wu73OINedPVemM98qpsWSczGYeeYpx6inSiHEZXsHFoEWme+udnur2fFssBM5yDSDtuTJdNpZhi6kCS0LWD4ewZ4vgBkZBihkOXgqQBER6gMleo9DNc5jrKwRMEA0vAtAhVLavI3LYydNrsYqbdQuZ1cQndOTJalzGcO5rKLmi6M1tzmgW1ST59/jE8/uou+nJF6pJxfvnEZt6RjNOfL/Kt+55HCFjePoNj5s9UgkOpJXwoRmiHknZP6c5kUipxKLxMruMQ023lEzEo+mzdtoe5TWmEFpLirirh0icqZRgUHp7f6qMGF4MqcQW4MWitZ8mMOj6/chFs38tLW7vY1NHLE5t3ce9z2zly3Q8BuAz4LyFIhddqrh9wySObeP6VDj5zWBu5gs/JHd18eHc/8YKnPmYVh9Sv07055m3cyjMNKY7oyZX/gSEZQ/z12bDqCGJS0ugHNPblePylDr51//M8s3UPHzx9GUFzLW7XwMwhObMWcfMVXPz3/wPA9e891Qg0B/nqLr7yi8dZe9aRNM2fySNPvs6cpnTZNpMWW1JmsVgsllFgJGVlpwAXhrlCKaBBCHEr8FZUDhHAj4Cbq+0spdwW/uwQQtwOnAgMEIdCR9FNoP4qvI/nYRlHDkSssOw/+3vdBxOIOmT1YOpO0gOOq4+pfjao+7ZrdqtuZfNUwHR9WCLWCmyYq0rSKjOLNLJGlNrdj1fukMaWl1kOCkYkEFHu8BkOUxgazESkx/GNW3zd3t0xBJTAOKAWiLQQY5aOldY5pbIgpFRiSuW5VM6zEl+LKlIJMYO5i0q5Q+HrZAxihghkOqG00OIKEKHrJwggGWfxke0sXtTKzm17+fXTr7Nxyx5+9KeXeWZrF9mCR94LuOkv36REJ31dtFOo6EdinHnNCl7obHKhUKBzZw9Swoy6JLF0srTv3KY0ojYRiWDpZMW8pcobOtB2864ANw6LWlm8sJXFns9b9ma47OUOnnp9D89u7+LSH/6RVEU3NAEs3dPP/xydQhw3G/H/7kMU1O+C4aSWWmBWT44uYEa4LKhN0ve2lbiLZ1PX1Q/JGPnuDD/+0yayRY+/WHkoN7z7jYiGGnAc5E0bBjiTeg6fQ9+mXdzz3Ha1sKEmcpP5AWJOI//nHSdCzGXbSzuZUZvk0BMPjd6zyYQVgywWi8UyBgwrDkkp1wHrAELn0FVSysuFEM8CbwbuAc4AXqzcVwhRCzhSyt7w+dnAF0Zt9pYJoWp7+9t6EBViwWi1Xq9kuuQTVWOVDLt/icw+Xedq1+xBZnMRrwwwEzzI7Kr7lolTaxqGzAsqCT4SrmAjbWTABekr8an1uvqyeVuByGLZD0acQUS5W6fMzVOxjd7ODJuu1Bf0clN4cAxRRbvGtUBkikKVLh+NdgDFHOWkCYLIsTFEtVsJPywhS0Sdwko/zayisnI74+RcysvIRMX8vDCIW+fpCKGe54vMamvgfact5S3L2+nPe6RiLm2NKZJ1KSPoW0TlaXocLQp5uuMbSiRzUK8FPLO1i3QiRlM6oUq9gL6+PHWzGlRJmadCnEvXxw2vn57faCJQDqXWetpb6mk/YSHnZQrwg4cG3Vze8yyv3/MsC/fxUC2Uv+VOf57dj7/Kk01pHnt1F3/e1kVLXYoTFrZQX5Ngw7Nb+emjm2iqTZJwHc6Y18TKTbvKmt7VPfQSDxtCT2ZnD+k/vgz/eY/KPmprhPefBqcezq7eHEfMaw7f632c/FhhBSGLxWKxjDH70q2skg8CNwghYkCOsCRMCDEXuFlKeR4wC7hdqC+DMeB7UspfH9iULZOBsvb2t/UgruqIyoyMblRjJRBpDiahaJXczBVspJWMKs1iuQqFNtaXnDiS/b7O+pqJlZsRW8rXCeCi9s38eUtLNJc5aVpvrN8vAWeDmM8GwnMwawk+Fj4Yf2HIYjmoqCYQmTlAZcuNMi3zrlnv4wgVBhYEUccuHdtTNgZquzIMu9GgJWzhwqHcLEJEXazKjlsxV00goyDskvNpiEBsPU8tGkUt0kL3TjDwGHpeeu56PyGgKa3mkCkwtzZJsTvD1q4Mz2/vZk5TkebaJLE5YUhzrqgcTbp0TMpISNOXR89BSp55YQd/94OHeOtxC1jUWk+L6+BLSXe2QF3cVeOIUFzSOUra+TTWbhctFDXWqAyjHd1VN3NR7eb3Z/jK14ue3sK3d3bzlc7eMmMawLzmNJm8hy8lfiC5Iu8NGMOVcMz9UTvpv/rgLfy3KxDa9bSzG667E/yArVJy9MIWtXwk4uRYYMUgi8VisYwz+yQOSSnvQTmFkFI+AJxQZZttwHnh81cAe9d3kKLLjcQ1u8vyZ4AJ6UZVWW41lQSHVXIz62oeLl3HWWEbeCQlgegKNg4o0Tqg67x1kG4wW7wBc8mtzfPVtQtG/AV5sGtf6UIaatuxxrqHLAct2oVTyWBOoAH7i6hsrNK9o8c1M38gKosaTPgxM4AMjaZqyVi1ddIYQws/gW6b7gziDtITlBXPRUkPikQuGTl3kMqJVBLQdGmZU553FJYjqY5hMSh6PP7abuY0pjlkRh2+lMRqEpFLqC6pysYEUR5QyW0lI2Go6PPSs9tY96M/sac/T0+2wC+f2MyCljp6skUuOmEh1CbDTmte5HIqvW8M/hkYC648Az5/e9VVwhHQ1jCoeLQvCOCvOvt44aTD+PhZy3nj4bOVOFeTUEKVH4SuKeDsf646Rgtw8fEL+eljr3IDRMKQxg/gprtZ8ZV38PNHNnHhgpZIJB0vgciKQhaLxWKZIA7EOWSZ5ujyst8MrChUDCY+jBNTyVV0BRsHCGwpfNbyeMl50zpYuPP+Xud5MeU+qsSl6lyuYGPkAtoPqglDE40ViCwHBZXuoWo3smVuIlFFlDF/iuhmW4sWlYJNSdgwhBYzQ6gaZmlbUGUueo6mAGWOVRKGwn21KGSOV1V0EkOfu97BYaDopeekBQJ9TUrnQ1SCF3M58fDZ9PfnyRZ9fD9QeUMxJxLPYq4qDzNFNy1mhQ6g+/7wEl+7ayPL5zWzbE4Ty+Y2c8rhs4m7Di11KSWI5IvRPNIJ6M9HLdljbtg1rop7bCxYfRTc8QQ8smnguotPgKMPgWvvUM4pjXY5FSsy6RproDs76KEOQfL9f760XBDU6NMNgmq7Amrz2695O3RlkJfeWH2jnhyt82eSem0Xrzz0MoceN18JUGZZ5oFmOVViBSGLxWKxTAKsOGQ5IO4KfsROUVu9K9W8yfPxmuyuosE6hzVQZK18jBvF8XRSPUB6f6+zXPf/27v/ODvq+t7j78/uEgKBBcEEG5ZfvSAtoKRlQSxVIYjRELDGq8Urtn3UmlsKjwdC0ZJL60OvP2hrvBf746HNxd778HLrDzS0KNRUC6ktt1USIAgE5aeabC8JFA2SIE32c/+Ymd3ZszPnzDln5sycmdfz8ViyZ87MOd+ZswnnvPfz/XyPlOLTASX5QSbtTf5AscT29LxiWuvtKl1/AiLUQmtA1FrFE5/W1bFySMG/R6Mjc/9dSjou62fkOcvEazZQiYKR1h5IpnDKVUtYFE1rk83tHRQfR1IwFn/8+HNE1UPT0y1hQ/j9vnDK2v5w+loUDkWNpOMrsYVTxhYtGNUBu1/Qgmj1s2g1spGR2Wv77/vDx52eXfHNpH97Ypc+9pV79MZXHKNTjz5Ci8cX6qSjDtOiA8e0b9o1tvQlQZXS/nAsCw8IHuu5F+aOq2Pb55z96bukj98u/fWW2Z+9XzlDet/K2X0+fUcwdeuow4Jqo6RtK14hvbpNa8ol43NDwpmphPEA09JDpsMOkvZN68KP3aqvtjufffv1htOP1Sf/9n698pnndN7kCcEYx0YS/1/dNcIgAEAFVefTO4bWUTetl7/nXfNChmg1qiqqXBVLShWPSXqzPa6L/+x5SQlLyPdznVcnrzpm1z+TXFHURQjVWrU1TFVcKJaZjUraLGmHu68ys2WSPq1gNcx9kn7H3b9d5hiHVhQQzVTxKL1CqJN4VdC8+6LHankOC6tj4tU7M1PBYo8b3ReFK3OWpdfcAGnmmNhx8ZAraZn7OWOM7TcTLGl2mlpUBTLnccIqoqgJ9Yv7Z8c4OiItCCt0on5HUR+nEQUVPYsWasGCsdkV16LnjyqwFowFTaT3vijteTFYaWzfft378L/qjm07NDYyol847qU66WWHaXraNR2OcWz8oLlVMdH1WTAWfP3kBenwRbPTzKIx5l3lkuZ9K+eGQXErXhF8JW1v9bL0Hka67PwgoJn2ln5Y0bRAC0LF966QPvaVuZVJB4xKV71RevZ5XXb+qdqzbUqLWiuXIv/rH6VLztaVa86TnnlO39z8pPZPT+u8s08MAqrwKbuaakYgBACouAH/agm19M536vq9Z8onxjQt6SkdLF+3pPBm1HnYOLV15qvMMfjaI1N/GWmuILBZPS5ft0Q+MSY3BX/2e51Xj8vvPkE+dZL87hOC22uPDCqIYroNoVYsPX3mK7odKfNap5m3+h6KdKWkbbHbfyzpQ+6+TNIHwtvoVesH0F6rHOKrj7WaCXFioU181+lYQBNVdERfIxYuC29zQ6bp2If7mbGH/X2iPjojFq7mFZ+K5i1f0rzwyVrOxRU23Y4qT2x2nNED7wv7+eybnh3XAaOz08Giyp2owuqAsDn0gQcE08YWHRhU9YxY0Gfop/tm95+ejjWRHtUL//YTbbr3+/rFD3xZn/r7h7Ty9GP1ulOP1tIjD9GiAw/QoYculA5ZGD7eSOx6jMxe6wVjs1VOFr9OAwqG8vTby4NzbbV6cjZMin6Ooq9o5TcpeK3f8ArpuouDoMkU/HndxcHx/2GJVr31TPk7fyl9DLfeK931SPDYP3O4Xrv853XEIQv1N3du04+2/kB69vkg3ItWmtsTBn1xZ39o7hcAABVH5RBy8V/8W8EH7Oh9aLgaVRWDgDTxaqJuK116mSo159qsHpd+f5f0bEqvhKivUIcl5HORUlGUZ9hXtallEtPLBsHMJiRdKOmjkq4ON7uk6IfrMElTJQytXuIVRJ2mhbWT1mQ6KTOKdummeW98Sln0XPHHip4nCmTiU9Hi98criOLVR9F+c6aUxaYjjUahkc9+yD9gNKg2eTH8N3dsdG7IJAXhj2s2LDIFU8NGR2Z7DM0MYiwMiPYHVUVRNdHBYaPq/dN6/Ild+vht9+nIQw7UK485QuefcnSw3/S0XnLIgcHzHLIwXPY+rHSZdmk6HMfYSDBdau+LwfkcMDp7Lac9OM9hEgVASVPOshgbmX2ctGNGTYe89cygQijNF78lXRoGSIsO1Onn/bxO//Ee3X/P9yVJP3jmJzph8aGaenaPbt/6A/1474v64C2bs40RAIAKIhxCbqIG1XGtAcAwhEVJPXKk5OqX1iCp1+eQJP/IYumKp2RJH7yK7N+0YXdQmdQSBOUdQrVeqyoGRCjcDZLeL+nQ2Lb3StpoZusUVLO2+XU+MosHRP1ImpIUn0aWJMtzJk0JyxJitB6XFkS1Tj+bmSKXsHNUkRSFRiMjQV+fmWbb4TFR4BJVCUUPNe3SmM3dJgW9jBaMBFU9B7UEVdJMdcwpZ/2sbvvZxcF9By0IVjWLlruPjokCsiiYiiqjIqMjQXVM1NdotM21GQbtgp28vPSQIHh6KmUK29PPzb09YtJLFumVrz5RGh3R4funpeXX62RJJxc7UgAABoJwCAPVa5hSBUnj7qV3Uer5rx6Xf3uv9NndcwKiQvs3bdgtizel3r5PumZn8Av5nMOhdtdPqkYvIqqHimNmqyTtdPctZnZu7K7LJF3l7l82s7dL+oyk1yccv0bSGkk69tjeV81DzbQLPyzl+zn7xCuOwu9HLfmAtPAqrafPnOdPCtnCbQcvkI7N+G98tCJZktGWTgHDGgwN0mXLpQ/eknzfUYclb19+fXHjAQCgROZpPQVKNDk56Zs3U5o7rLrp3zKsQVGSXAIiKbWSpwh25hOyhObTPjEW9CDqUtYpeP0+RtGqHBCZ2RZ3nyx7HN0ys+slvUtB0+mFCqaSbZB0kaTD3d3NzCT92N3b/sDz/4guxVcxAzDXx2+XNrT8e7LwAOnaVUH10hD2CxrW/0/kadyO8FfZ+WUPA6iM/3dV+YXZzx+X0j5jQF5+47OlPr8k/eSkw8segiTprg3vS/z/BA2pkbtuPljHmxYPu6TG1mnb2kpoEl2YHQmrkrXbnoNOr3kVAkMaVOfP3de6+4S7Hy/pEkl3uPulCnoMvS7cbbmkR0oaIoAmet9K6YNvmdu8eoiDIQAAesW0MlTCME83a5V2HpXssXP0WN/L1ke6Pbd2r3nVppuhUO+R9EkzG5P0gsKpY8hR1H8IQLKoxxFhEACgwagcQiF6mZZTiyqiDbuDqVpLH5Gd+YS0YbckaeOf/aPszCf0d/rSnO1ly2PZeqn3AKfKrzfVQ8Vx903uvir8/p/c/Qx3P93dX+XuW8oeXy3xoRdIxnLzAABIonIIBUpavSyLKDAotJKoiL4+ac2dv71X9sXnBtL0uWsDWrY+krTiW5XRoBq1QgURECAIAgBgHsIhFKrXgEgqMEhICnEuf0r6/V3BcvI9BiN2/TOzROulrwAAH7xJREFUjxlt2+vSTbtl+zV/+/XP5L5cfE8KWLY+mkLXdb+lFlWuLAKGEgERmogwCACAjphWhqGQNuWsl/AgMcSRZM9Oy654Srr2qc4PkjR9LK2J8/7kzUU2fR6ktNellyCo9bGqUF3E9DIAGEJMFwMAoCtUDqFw/VQPtWoXEGUOEtqEMuaSPrtbftZB6RVEKdPHdPiI9Oz8JRpt3pZQD02fC9XDVLu0cC7La5F2bPz1rErlENPLUCtUD6GOCIEAAOgLlUMYiEF8sM7c0LpDKGMeVBel3p82fUw+r7lzmtSmzykNrQsXBV7b9wXnv32f7JqdmZ9/49TWOV+dZFnKvirBUIQKItQKH6RRB1QHAQCQG8Ih1E6nkChpha552k35SrvvRy5ft0Q+MSZP3kMuBfevWzK/KqfPgKYfaYFXu5Cs3TXuFOwkBUjxYKmboAlAj/hAjWFFIAQAQO4qNq8FdZbn9LIsUhtaRyt0/f5O6VlPnvbVrrro6LFgKlnS9rC5s535RPI+E2Pyu09IfNjUiqRBNK5OC7zahGTxCp9O/aBaQ55cpgeWgOllqB2mmGFYEAYBAFAoKocwUJX5YL16XP7QifJfH5e3pEOpU76i+xMqj1qPybLPPGlBzPZ9xVcPHZ5SSZW2PaaXMKdT5VC/j18kppehdvjQjapi2hgAAAND5RAGbtAVRFKbqpQ/PCpoPt1NI+ao8qjdMVn2aZVSkWSSdM3O4PEKqyAyKXEyXLYeShuntkobdmvX5c9psfZolw7WZ3Sa7rRjJc1fvayKq5IBjUYFEaqCIAgAgFIQDqFREqc6hVPBupLlmG4ed8Nu6flpuZLjmMKnl/1o/iprbbe3uvYp2Wd366jw5lHao6u1RXLNBERxVVqJrBdML0MtERChDIRBAABUAuEQSlFG9VCrLOFEa7VLu+qXtOM6ihpR701rYx1q1yS7X+36KHWyYbfss7tlLcNfqP1aO3GP7twxPxxq1XptW++rIgIiAOgRgRAAAJVDOITSVCEg6qQ1mMgaVCRNY0tr0JzUiDpRlqAmYYxZgipfe6TUElB17JEUsuufmRcMzdixb6ZKqPWaDHv1kERAhBqieghFIRACAKDSCIeAAqWFH3OCIn+kY2efLEFNu+fqGBBl6JGU9Pgbp7a2r2gKA60oCGo3jmFYrQxoBAIi5IEwCAAq4/njMraKKNDSf8jwy/ACfe+3XlLq80vSa85+qOwhSJLu2pC8ndXKUCqqLqRdOjhxu49KbpJPjMnXLZkJapKqmTpV4ET7tN1v9bj87hPkUyfJ7z6hYzA0sz2loslNM4FWdHxSVVPrKmXD1qy66tVvQE/4YI9esLoYAABDi8ohlG4YppcV6TM6TVdrixZq/8y2FzSq/7b/jKCZ8w5JV4RfoUFNx8rSV8k3JExJC/9j1z8jV7bKoXb9naqO6WUAGo9ACACAoUY4BJTsTjtWcundeiBxGfi8xYOXtPCnU1+gOcFNOCVt5+XPaYn2SKbZHkTb90nX7NT1lx8nWXJT73aNqOPjaDdeAAVgehnaIQwCAKBWmFaGSmh61cWddqwutZVaYf9Rl9rKwoIhqXPQEr8/SxizYunp0upxXWortVMHz2tObXtd79YD7Y/vcoxV1OTqN9QYAQDimDYGAEBtEQ6hMpoeEA1K1sbPrb2Akh4j2i+q8lmsPYn7L9Eenec/mHds0uPFHzfpuaqMgAi1RBDQbARCAAA0AuEQKoWAqHytFUOdApn4vmnNtU3S1dqi8/wHmQKn+HO37l/1gAioJYKBZiEQAgCgcTKHQ2Y2amb3mtlXw9vLzOxfzOw+M9tsZmelHPdGM/uumT1qZtfmNXAAvcs6XaxdlVHSfZ/RafKDLPHxFmq/3q0H2k5Xy7SqmtpXNZWN6iHUFkFBvREIAQDQaN1UDl0paVvs9h9L+pC7L5P0gfD2HGY2KunPJb1J0imS3mFmp/Q+XDQB1UODkbWfT6eAKP44d9qxun7vmfJ5ewaiaWdZgp1uQqKqhUUERKgtgoN6IRACAAChTOGQmU1IulDSjbHNLmk8/P4wSVMJh54l6VF3f9zdX5T0eUlv7n24aAoCourotjn0tf/6rDSRvBBiNO0sa7PrVu2OqVrzagIi1BZBwnAjEAIAAAmyLmV/g6T3Szo0tu29kjaa2ToFIdMvJRx3tKQfxm5vl/SqpCcwszWS1kjSsccWt1ITgO5lqeCJ/+lrj5Su2SnbO1tD5AeZFq87VLqi93FUqToIaDSWuR8uBEEAAKCDjpVDZrZK0k5339Jy12WSrnL3YyRdJekzSYcnbEucceLu69190t0nFy9e3GlYaACqh6oh83L2cavH5euWyCfG5Kbgz3VLpNXjXQc8VasI6gbVQwBKRYUQAADIKEvl0DmSLjazlZIWSho3s5skXaSgD5Ek3ay5U84i2yUdE7s9oeTpZ0Cir0/fzAfsITGvN9Hqcfnq8dT95uyr4Q6B2rlg5G0EnagnqoeqiTAIAAD0oGM45O5rJa2VJDM7V9I17n6pmW2T9DpJmyQtl/RIwuF3SzrJzE6QtEPSJZL+Uy4jBzCU0iqHNk5tbbuKWdpxdQ2VgKFAQFQNBEIAAKBPWXsOJXmPpE+a2ZikFxT2CzKzpZJudPeV7r7PzK6QtFHSqKS/dPcH+x00moXqIQx7hRHVQwByRyAEAABy1FU45O6bFFQKyd3/SdIZCftMSVoZu327pNv7GSRAQDS82lX9NAkBUcDMRiVtlrTD3VeZ2RcknRzefbikH7n7stIGiO5FIQUVRMUjEAIAAAXJtJQ9AHSj1+qepBCpTsESAaekoFfdtuiGu/+quy8LA6EvS9pQ2sjQH4KLYrD0PAAAGIB+ppUBA0X10HBoDYaicCdte6u0/dO2YXiY2YSkCyV9VNLVLfeZpLcr6GEHNBtBEAAAGDAqhzBUmJYzPPoNcjZOba1V1VCk4QHnDZLeL2k64b7XSHrK3ZMWN5CZrTGzzWa2edeuXUWOEf0g1OgdFUIAAKBEhEMYOgREwyNt2frW+5qmiQGRma2StNPdt6Ts8g5Jn0s73t3Xu/uku08uXry4kDEiJ4Qb3SEQAgAAFcC0MgCFqWPlD3p2jqSLzWylpIWSxs3sJne/NFz1crUSFjkAaotACAAAVAjhEIYS/YeqqcgwqG6VRk1bvczd10paK0lmdq6ka9z90vDu10t62N23lzQ85O3sD7F6WSvCIAAAUGFMK8PQatIH62GR1kS632CnbsFQhIBzxiVqM6UMQ4owhD5CAABgaBAOAcisU2VQu75CK5aezjQzSJLcfZO7r4rd/g13/3SZY0JBmhqKEAgBAIAhw7QyDDWml1VHlobTrdui22mhUV0rhuKaNr0MDdSEKWYEQQAAYMhROYShxwfrYrUGN0lBTnxbHtPImoaAE7VX1/CECiEAAFATVA4B6Eu/U8WaXDUURwURaq8uFUSEQQAAoIaoHEIt8KG6WBuntmYKgXoJdNKOoT8RUEPDGqzQWBoAANQc4RBqg4CoGGl9gqT508nyeo5+H29YMb0MqBgCISA3ZjZqZvea2VfD28vM7F/M7D4z22xmZ6Uc96SZfSfab7CjBoDmIBxCrRAQDVYU4BQR5GStVqobAiLU3jCELYRCQBGulLQtdvuPJX3I3ZdJ+kB4O8157r7M3SeLHCAANBk9hwC0lVYd1E9wk/XYJlYPAY1Qtf5DBEFAocxsQtKFkj4q6epws0saD78/TNJUCUMDAISoHELtUD2UnywrlXUb4HQTKjWxckiieggYGCqEgEG5QdL7JU3Htr1X0sfN7IeS1klam3KsS/o7M9tiZmuKHSYANBfhEGqJgKgYURDUT2jTTZjU5MohAiLUXlmhDM2lgYEys1WSdrr7lpa7LpN0lbsfI+kqSZ9JeYhz3P0XJb1J0uVm9tqU51kT9i7a/O/6aV7DB4DGYFoZgETtAqAVS0/XxqmtPYc3eYRMTcDy9qi9QU0vIwgCynSOpIvNbKWkhZLGzewmSRcp6EMkSTdLujHpYHefCv/caWa3SDpL0jcT9lsvab0kjdsRnvdJDJvR004uewja/8B3yx4CKuSQR35U6vMvOv7IUp9fkqbeu7vsIbRF5RBqiw/V/Wm3SlkkS7gTNZZOajDdLlxqctUQ0ChFBTdUCAGV4O5r3X3C3Y+XdImkO9z9UgU9hl4X7rZc0iOtx5rZIjM7NPpe0hskPTCQgQNAw1A5hFr7+vTNTM/pUV79hqIqo6T9qRzqjOohNEJeFUQEQcAweY+kT5rZmKQXJK2RJDNbKulGd18p6ShJt5iZFHxu+St3/1pJ4wWAWiMcAtCTrJU9acFQtK2baqKmIiBCI/QaEBEIAUPD3TdJ2hR+/0+SzkjYZ0rSyvD7xyXxxgAABoBpZag9PlT3b8XS02dCm/j3WY9Nk1Q5RDURgMwIhgAAAHJBOIRGICDKJq03UKSfBtTdTEsjIJqP6ZFohCxhD72EAAAAcse0MgCS5gYyRUzt6rahdbvpaE1FQIRGSJteRhgEAABQGCqH0BhUDyVrVymU9/O0yhL8DGp8ACokCoKoEgIAABgIwiE0CgFRZ1TqAKgEAiEAAICBYVoZGofl7WcVPZUs6Xmi2/EG171WFTVJFGyGy/k22pYtW542s+/n9HAvlfR0To9VN1ybZFyXZFyXdIO+NscN8LkAADVBOARgIFrDHqaKdUalWzJ3X5zXY5nZZnefzOvx6oRrk4zrkozrko5rAwAYBoRDaCSqhwbT8Jmm0r0hFAIAAAAwSJnDITMblbRZ0g53X2VmX5B0cnj34ZJ+5O7LEo57UtJzkvZL2sdvTlAVTQ6Iiq7ayRoKpU0payICIQAAAABl6aZy6EpJ2ySNS5K7/2p0h5l9QtKP2xx7nrszDx2okCKreagUyo5QqHTryx5AhXFtknFdknFd0nFtAACVlykcMrMJSRdK+qikq1vuM0lvl7Q899EBBWti9VAVKnWaHh4RCFWHu/OhLQXXJhnXJRnXJR3XBgAwDLIuZX+DpPdLmk647zWSnnL3R1KOdUl/Z2ZbzGxND2MECtXED+pVC2eSxlOFECtvX5++uZE/bwAAAACqrWM4ZGarJO109y0pu7xD0ufaPMQ57v6Lkt4k6XIze23K86wxs81mtnnXrl2dhgWgB1ULXOLjqXNARCgEAAAAoMqyVA6dI+nisLH05yUtN7ObJMnMxiStlvSFtIPdfSr8c6ekWySdlbLfenefdPfJxYtzW6UYyKRJH9yrUjWUtWn1sAZEUSDUpJ+tqjOzt5nZg2Y2bWaTse1nmdl94ddWM3tL7L4zzOw7Zvaomf1JOJW6VtpclwvCqt/vhH8uj933jnD7/Wb2NTN7aTmjL06P12WBma03s++Z2cNm9tZyRl+sXq5NbJ9bzeyBwY54MLq9LmZ2sJndFv6sPGhmf1je6AEATdcxHHL3te4+4e7HS7pE0h3ufml49+slPezu25OONbNFZnZo9L2kN0iq5RsCDL86f4jfOLVVG6e2ViYYSpM0vqqPuRWBUKU9oOAXGt9M2D4Zrrj5Rkl/Ef7yQ5I+JWmNpJPCrzcOaKyDlHZdnpZ0kbu/QtKvS/rf0swvhj6pYLGJV0q6X9IVgxvuwHR1XULXKai2frmkUyT9wyAGWoJero3MbLWknwxkhOXo5bqsc/efk/QLks4xszcNZKQAALToZrWyJJeoZUqZmS2VdKO7r5R0lKRbwl+0jkn6K3f/Wp/PCWDIZa0aGhaEQcPB3bdJUmvxj7vvid1cqKBXnszsZySNu/s/h7c/K+lXJP3tIMY7KG2uy72xmw9KWmhmByroP2iSFpnZMwpWMX10MKMdnG6vi7v/VNJvSvq5cL9pBaFA7fRybczsEAWLmqyR9MVBjXWQergueyTdGe7zopndI2liQMMFAGCOrsIhd98kaVPs9m8k7DMlaWX4/eOS6vHpD41Qx9XLqhjErFh6euJ0sWGbQkYoVB9m9ipJfynpOEnvcvd9Zna0pHhl7HZJR5cxvgp4q6R7wwBEZnaZpO9Iel7SI5IuL3FsZZq5LmZ2eLjtw2Z2rqTHJF3h7k+VNrpyzfmZkfRhSZ+QtCf9kEZovS6SpPDn5yIFVXkAAAxc1tXKgMao4wf+XoOhosKatCluaU2pqxYaMXWsuszsG2b2QMLXm9sd5+7fcvdTJZ0paa2ZLVRQHTNv1yLGXbRer0t47KmS/kjSfw5vHyDpMgXTYJYqmFa2tsDhFybP66LgF24Tku4KF+L4Z0nrCht8wXL+mVkm6UR3v6XgYRcu55+ZaPuYgkr8Pwl/sQoAwMD1O60MqKW6VBD1E6oUHcgkBURVC4HiCIOGg7u/vs/jt5nZ85JOU1ApFJ/iMSFpqp/HL0uv18XMJhQsJvFr7v5YuHlZ+JiPhft8UdK1eYxz0HK+Ls8oqIqJApCbJb2770GWJOdr82pJZ1iwuMmYpCVmtsndz81jrIOU83WJrJf0iLvf0O/4AADoFeEQUFO9TCfrtLR8HjqNK769CmERoVD9mdkJkn4YTiU7TtLJkp5096fN7DkzO1vStyT9mqQ/LXOsgxROc7lN0lp3vyt21w5Jp5jZYnffJekCSdvKGGMZ0q6Lu7uZfUXSuZLukHS+pIdKGWRJ2lybTylo7i4zO17SV4cxGOpVm79LMrOPSDpM0m+VMTYAACJMKwNS1CEUyBLwRNO24qFN2vSuQSsrKGIZ+noys7eY2XYFVQy3mdnG8K5flrTVzO5T8Jv933H3qJHwZZJuVNBw+THVrBm11Pa6XCHpREl/YGb3hV9Lwt6CH5L0TTO7X0El0cdKGXyBur0u4X2/J+mD4XV5l6TfHfjAB6DHa1N73V6XsJroOgUr290TbickAgCUwtyr1z5hcnLSN2/eXPYwAEka2ullnSp0eqkSSusV1O3+rWNLGmtrGDSIhtrDEAaZ2RZ3nyx7HAAAVNW4HeGvsvPLHkapRk87uewhaP8D3y17CAg9dsPZZQ9BL7/x2VKff8cFR5b6/JL0sv/+f8segiTpG/6lxM8TTCsDaq5dyJJH0JP2XGkVSK3PHX/cdtVB3QZT3RiGUAgAAAAAikI4BHRQh+bUWVcG6/QYWfsYtbu/jIqgJARCAAAAABAgHAIyGMaAaNChS7c9gaJKoCzH5XkuhEIAAAAAMBfhEICudGpW3U3VUKdgKK9QiEAIAAAAANIRDgEZDWP1UF7iIU673j9Zppt16jFElRAAAAAADBZL2QNdaGLY0BoGtYY3rcFR1scsEsvQAwAAAEB2VA4B6ChrmJOl6iepsXXWyqR2CIMAAAAAoDeEQ0CX6j69rJtpY2k9hNqtbJb0+FmWsk9DKAQAADDX/ge+W/YQKmH0tJPLHoJ2XHBk2UOQNF32AEp32JP7yx6CHrvh7LKHELjyS4mbCYeAHtQ5IOqmaidttbGkvkJZq4qyIBACAAAAgPzQcwjoUV0DitawJwps2lUTtQt14iFRr1PGIvQSAgAAAID8EQ4BmNEuvGk35Suteqi1aqjXYIhQCAAAAACKw7QyoA91nl6WJC08ytqLqBuEQQAAAAAwGFQOAX2qY4iRNrUsab+k4Cc+1aybnkMSVUIAAAAAMGhUDgGYJ2o03Wmp+Wi/dkvTZ0EYBAAAAADloXIIyEEdwo0o0OllCpjUWzBElRAAAAAAlI/KISAnw9p/qLUSKNoWvz/SWikU354m6T4CIQAAAACoDsIhoMFag5/WbVmWnm+9PwqQWhEIAQAAAEA1EQ4BORqm6qGk4Kdd8+ho/7TwJ358ZMXS0wmFAAAAAKDi6DkE5GxYwpCkYCitb1C7CqK0aWYjL/ve0FwLAAAAAGgyKoeABkmaRpZF0tSx1seTpJGXfU+S9PXpXkcIAAAAABg0wiGgAFWcXpalf1C3+0f3B6uO9TU8AAAAAEBJCIeAglQtIOqlsXQ7TBkDAAAAgHogHAIaKq35dLuG0xKhEAAAAADUDeEQUKCqVQ/FdTPFTCIUAgAAAIC6IhwCClblgChNvJcQAAAAAKDeMi9lb2ajZnavmX01vP0FM7sv/HrSzO5LOe6NZvZdM3vUzK7Na+DAMBm2kCVoMD1cYwYAAAAA9KabyqErJW2TNC5J7v6r0R1m9glJP249wMxGJf25pAskbZd0t5nd6u4P9TNoAPkjDAIAAACAZspUOWRmE5IulHRjwn0m6e2SPpdw6FmSHnX3x939RUmfl/Tm3ocLDK+qhi9UCQEAAABAs2WdVnaDpPdLmk647zWSnnL3RxLuO1rSD2O3t4fb5jGzNWa22cw279q1K+OwgOFSpRCGUAgAAAAAIGWYVmZmqyTtdPctZnZuwi7vUHLVkCRZwjZP2tHd10taL0mTk5OJ+wDoD2EQAAAAAKBVlp5D50i62MxWSlooadzMbnL3S81sTNJqSWekHLtd0jGx2xOSpvoZMDDsyli9jFAIAAAAAJCm47Qyd1/r7hPufrykSyTd4e6Xhne/XtLD7r495fC7JZ1kZieY2YLw+FtzGDcw1AYV1jB1DAAAAADQSTerlSW5RC1TysxsqaQb3X2lu+8zsyskbZQ0Kukv3f3BPp8TQBuEQQAAAACAbnQVDrn7JkmbYrd/I2GfKUkrY7dvl3R7rwME6irv6WWEQgAAAACAXvRbOQSgD/0GRARCAAAAAIB+EQ4BQ4hQCAAAAACQl44NqQEUq5ughwbTAAAAAIC8UTkEVEC76WWEQQAAAACAIhEOARVFKAQAAAAAGARz97LHMI+Z7ZL0/YKf5qWSni74OaqOa8A1aPr5S8N3DY5z98VlDwIAgKoa0GeJIg3be5Ne1P0cOb/hV+dzTPw8UclwaBDMbLO7T5Y9jjJxDbgGTT9/iWsAAACqpQnvTep+jpzf8GvCObaiITUAAAAAAECDEQ4BAAAAAAA0WJPDofVlD6ACuAZcg6afv8Q1AAAA1dKE9yZ1P0fOb/g14RznaGzPIQAAAAAAADS7cggAAAAAAKDxGhEOmdnbzOxBM5s2s8nY9gvMbIuZfSf8c3nsvgVmtt7MvmdmD5vZW8sZff96PP9NZvZdM7sv/FpSzujz0cs1iO1zq5k9MNgR56/Hn4OvmdnW8LhPm9loOaPvX7fnb2YHm9lt4d//B83sD8sbPQAAGHZFvic3s7Vm9mi474pBnE/CGAp5v21mx5vZ3tg1+HSR59FOkZ8phvU1zPJ5oSqvYZGfh6rw+vVrrOwBDMgDklZL+ouW7U9Lusjdp8zsNEkbJR0d3nedpJ3u/nIzG5F0xMBGm79ezl+S3unumwc0xqL1dA3MbLWknwxslMXq5Rq83d13m5lJ+pKkt0n6/KAGnLNezn+du99pZgsk/b2Zvcnd/3ZwQwYAADVSyHtyMztF0iWSTpW0VNI3zOzl7r4/19F3VuT77cfcfVmeg+1RIec45K9h1s8LVXgNCzm/Cr1+fWlEOOTu2yQpeD3nbL83dvNBSQvN7EB3/6mk35T0c+F+0wp+YIZSj+dfK71cAzM7RNLVktZI+uKgxlqUXq6Bu+8Ot49JWiBpaJuU9XD+eyTdGe7zopndI2liQMMFAAA1U+B78jdL+ny4/xNm9qiksyT9c/+jzq4J77cLPMehfQ2H6fNCgedXidevX42YVpbRWyXdG/4FPjzc9mEzu8fMbjazo8oc3ADMnH9s2/8My/7+wFr/BtVT6zX4sKRPSNpT3pAGbt7PgZltlLRT0nMK0vI6S/p7oPDfhIsk/X0powIAAE3Ry3vyoyX9MHZ7u+ZWHlVJr++3TzCze83sH8zsNYWOsH+9nOMwv4ZZPy8My2vYy/kN0+uXqjaVQ2b2DUkvS7jrOnf/mw7HnirpjyS9Idw0pqBC4C53v9rMrpa0TtK7chxyrnI+fykoX91hZodK+rKCc/9sXuMtQp7XwMyWSTrR3a8ys+NzHmphCvg5kLuvMLOFkv6PpOWSvp7TcHNXxPmb2Zikz0n6E3d/PK+xAgCA+inpPXlSYFRI9UZJ77f/VdKx7v6MmZ0h6a/N7NRYRUeuSjrHoXwNIxk+LwzsNSzp/Ab2+hWpNuGQu7++l+PMbELSLZJ+zd0fCzc/oyDZvSW8fbOkd/c9yALlfP5y9x3hn8+Z2V8pKIurdDiU8zV4taQzzOxJBX9PlpjZJnc/N4+xFiXvn4PY475gZrcqKJmsbDhU0Pmvl/SIu9/Q7/gAAEC9lfSefLukY2K3JyRN9TKOTsp4vx1WcPw0/H6LmT0m6eWSCumNWtJnimF9DeOPm/p5YZCvYUmfhwb2+hWp0dPKwqkit0la6+53Rdvd3SV9RdK54abzJT008AEWLO38zWzMzF4afn+ApFUKmnfVTpufgU+5+1J3P17SL0v6XtWDoV61+Tk4xMx+Jvx+TNJKSQ+XM8ripJ1/eN9HJB0m6b1ljA0AANRfDu/Jb5V0iZkdaGYnSDpJ0reLH3k2/b7fNrPFFq4QZWY/q+D8KlXNncNniqF8DbN+Xqj6a5jD56FKv35ZNSIcMrO3mNl2BcntbeGcQUm6QtKJkv7A5i8P+XuSPmhm9yso3/zdgQ88Jz2c/4GSNobnfp+kHZL+Rxljz0uPPwO10sM1WCTp1vDnYKuCebalLR3ar27PP/ztwXWSTpF0T7j9t8oZPQAAGHZ5vic3s4vN7L9Kkrs/qKDR8UOSvibpci9hlaQ832/Hz0/SayXdb2ZbFfR7+W13/7eCTqOtos5xiF/D1M8LVXwNizq/qrx+/bKgSAYAAAAAAABN1IjKIQAAAAAAACQjHAIAAAAAAGgwwiEAAAAAAIAGIxwCAAAAAABoMMIhAAAAAACABiMcAgAAAAAAaDDCIQAAAAAAgAYjHAIAAAAAAGiw/w9X04K/bMi/RAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 1440x432 with 3 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, axs = plt.subplots(1, 3, figsize=(20, 6))\n",
"axs[0].pcolormesh(mesh.glamt[0], mesh.gphit[0], mesh.tmask[0, 0])\n",
"rio_show(dataset, ax=axs[1], cmap='copper_r', vmax=5);\n",
"npoints = 100\n",
"lat = numpy.zeros(npoints)\n",
"lon = numpy.zeros_like(lat)\n",
"for ii in range(npoints):\n",
" lat[ii], lon[ii], x_index, y_index, data_out = ros.get_lat_lon_indices(\n",
" geotiffs_dir,\n",
" spill_date_hour.month,\n",
" watermask,\n",
" mesh,\n",
" random_generator=random_generator,\n",
" )\n",
" if ii < 10:\n",
" print(lat[ii], lon[ii], x_index, y_index, data_out)\n",
" for ax in axs[:-1]:\n",
" ax.scatter(lon[ii], lat[ii], color='r')\n",
" axs[2].hist2d(lon, lat)"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"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.3"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment