Skip to content

Instantly share code, notes, and snippets.

@arbennett
Last active March 11, 2022 18:05
Show Gist options
  • Save arbennett/1e8ef6fb9155caa6c5e302ed7e51be64 to your computer and use it in GitHub Desktop.
Save arbennett/1e8ef6fb9155caa6c5e302ed7e51be64 to your computer and use it in GitHub Desktop.
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "markdown",
"id": "42488a49-7587-424d-accd-dc25a492e6cc",
"metadata": {},
"source": [
"# Example of using `bmorph`s functionality for 2d pdfs\n",
"\n",
"First, just some imports and plotting setup."
]
},
{
"cell_type": "code",
"execution_count": 1,
"id": "a1be0fea-541b-460b-ac67-142662bc4aaf",
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"/pool0/data/andrbenn/spatially_consistent_bias_correction/bmorph/bmorph/core/workflows.py:6: TqdmExperimentalWarning: Using `tqdm.autonotebook.tqdm` in notebook mode. Use `tqdm.tqdm` instead to force console mode (e.g. in jupyter console)\n",
" from tqdm.autonotebook import tqdm\n"
]
}
],
"source": [
"import numpy as np\n",
"import xarray as xr\n",
"import seaborn as sns\n",
"import matplotlib as mpl\n",
"import matplotlib.pyplot as plt\n",
"import bmorph\n",
"from bmorph.util import mizuroute_utils as mizutil\n",
"sns.set_context('talk')\n",
"mpl.style.use('seaborn-bright')\n",
"mpl.rcParams['figure.figsize'] = (12, 8)\n",
"mpl.rcParams['font.size'] = 30"
]
},
{
"cell_type": "markdown",
"id": "03a2abb1-e9a7-4385-8223-776fb0b94c6e",
"metadata": {},
"source": [
"## Setting up of data\n",
"\n",
"This is the same as in the tutorial, so you can copy and paste the code following this cell into the binder instance to play around with this. All this is doing is loading in the data and setting up some metadata. For more details see the bmorph tutorial: https://notebooks.gesis.org/binder/v2/gh/UW-Hydro/bmorph/tutorial?filepath=tutorial%2Fbmorph_tutorial.ipynb"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "3495fad4-bce6-4267-8e17-e4ad25c01d94",
"metadata": {},
"outputs": [],
"source": [
"site_to_seg = {'KEE' : 4175, 'KAC' : 4171, 'EASW': 4170, \n",
" 'CLE' : 4164, 'YUMW': 4162, 'BUM' : 5231,\n",
" 'AMRW': 5228, 'CLFW': 5224, 'RIM' : 5240,\n",
" 'NACW': 5222, 'UMTW': 4139, 'AUGW': 594, \n",
" 'PARW': 588, 'YGVW': 584, 'KIOW': 581}\n",
"\n",
"seg_to_site = {seg: site for site, seg in site_to_seg.items()}\n",
"ref_sites = list(site_to_seg.keys())\n",
"ref_segs = list(site_to_seg.values()) \n",
"\n",
"yakima_topo = xr.open_dataset('./topologies/yakima_huc12_topology.nc').load()\n",
"\n",
"yakima_met = xr.open_dataset('./input/yakima_met.nc').load()\n",
"yakima_met['hru'] = (yakima_met['hru'] - 1.7e7).astype(np.int32)\n",
"\n",
"yakima_raw = xr.open_dataset('./input/yakima_raw_flows_full.nc').load()\n",
"yakima_raw['seg'] = yakima_raw.isel(time=0)['reachID'].astype(np.int32)\n",
"\n",
"yakima_ref = (xr.open_dataset('./input/nrni_reference_flows.nc')\n",
" .rename({'outlet':'site'})[['seg', 'seg_id', 'reference_flow']])\n",
"yakima_ref = yakima_ref.sel(site=ref_sites).load()\n",
"\n",
"yakima_met_seg = mizutil.to_bmorph(yakima_topo, yakima_raw, \n",
" yakima_ref, yakima_met)"
]
},
{
"cell_type": "markdown",
"id": "40a0e00a-2008-461c-8f4f-5fadc656e98f",
"metadata": {},
"source": [
"## Using `bmorph.hist2D` and `bmorph.marginalize_cdf`\n",
"\n",
"Below are examples of how you can use these functions. Note that the `temp` and `flow` variables are just the timeseries as pulled from the dataset directly. You can pull out any other variables you would like if you want to experiment.\n",
"\n",
"The first calculation using `bmorph.hist2D` will just compute the joint PDF for the x and y variables. It returns the bin centers for each of the x and y variables, as well as the PDF as a 2d array. I will plot this up in a subsequent cell.\n",
"\n",
"The second calculation you set up which values you want to actually select to compute the univariate CDF for bias correction. I've chosen 3 values to do this with. The calculation needs the bins which are used to find the index to slice out of the 2d PDF as well as the 2d PDF itself."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "87791291-ce08-4e19-914b-258eabbf2c1c",
"metadata": {},
"outputs": [],
"source": [
"site_to_test = 'KIOW'\n",
"n_bins = 50\n",
"\n",
"# Grab your variables out, feel free to experiment here\n",
"temp = yakima_met_seg['up_tmax'].sel(seg=site_to_seg[site_to_test])\n",
"flow = yakima_met_seg['IRFroutedRunoff'].sel(seg=site_to_seg[site_to_test])\n",
"\n",
"# Note - sum over `pdf` will not be 1, see notes section for details:\n",
"# https://numpy.org/doc/stable/reference/generated/numpy.histogram2d.html\n",
"flow_bins, temp_bins, pdf = bmorph.hist2D(flow, temp, n_bins, n_bins)\n",
"\n",
"# Note - this needs to be a sequence even with a single element!\n",
"marginal_temps = [10, 15, 20]\n",
"marginalized_cdf = bmorph.marginalize_cdf(temp_bins, pdf, marginal_temps)"
]
},
{
"cell_type": "markdown",
"id": "20d00326-a72e-407d-a806-c5db5762e2f2",
"metadata": {},
"source": [
"## Plotting the joint distribution"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "51832eec-db34-4884-b174-d2835d5af4e1",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 1.0, 'Joint PDF at KIOW')"
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAwgAAAIKCAYAAACHumioAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAABl1ElEQVR4nO3deZhkZXmw8fuZfZiB2RRxR1CjLIn7FhBQcI2yuHxxBdQsiluiH26ooFGjQUWDGhcUJBqXfAyEZIhiYBAjiEaUAAqyKCDIMjDDMMzez/fHOd0UNdXd76mu6u4a7h9XXUWf857nvH26uuc8590iM5EkSZIkgBlTXQFJkiRJ04cJgiRJkqQRJgiSJEmSRpggSJIkSRphgiBJkiRphAmCJEmSpBEmCJL6KiIyIk6e6npIkqQyJgiSxhQR+9c3+e+covMfEhHHNjxmZV3n4dfmiPh9RPxLROzZVva3bWXviojrImJFRLw1IhaPco4c4/Xu7r/jjuc6NiIOaVB+17oeJ3bYt2d9LdZGxLPqbUfU5V/aofzc+jr8OCJWR8SGiLgqIr4QEbu1lV1YX+vzO8SZVZ8zI+LZHfa/p973ktLvU5LUH7OmugKStnvzga0TOP4Q4HDg2IbHbQTe0FKHp9ZxXhgRT87MK1rK3gC8p/7/ecCDgP2BzwDvi4hXZOY5Hc7xC+CTHbZf3LCu4/kgcApw+kSCRMRTgLOABJ6dmReNU/4BdfnHA2dT/QzuAv4EOAI4vL42ZwBk5l0R8TPgKRGxQ2be3RLuycBCYAtwAPBfbafbv67XeRP4FiVJPWCCIKmvMnPDFJ16S2b+c8vXX46IXwHHA28FjmrZt6atLMCHImI/4N+AMyLi8Zl5VVuZ33c4blqqWwvOAO4EnpOZl41TPoDvUiUHf5WZX2rb/2mqm/l/qROu4XjnAk8DngH8oOWQ/YG1wPfq/2+NNasuf2lm3tbN9ydJ6h27GElqrO4u8q6IuLzucrIqIpZHxN4dym4zBmF4W0Q8PSLOi4h1EXFbRHwlIha2lFtJ9dS/vUvPEV1W/Xv1+yNLCmfmecA7qJ5897rb0Jsi4vt1d59NEXFTRPxzROzaUmbXiMj6y8Nbr0HDcx0MrABuBv50vOSg9mfAvsB32pMDgMy8BvgrqtaZ41p2nVu/H9B2yP7Af1O1HDwlInZo2TfcunAukqQpZwuCpG58A3g5VbeTLwC7UD2RvyAi9s3Mki42jwP+Hfga8E2qG8jXA0PAX9ZlPkL1IGNf4DUtx/64y3o/qn5v8pT6VOBE4AUd9s2OiPu1bRvKzNsL4r4TuBD4LHA7sBdVl6hnRcTembkKuJXq+z4VOB/Y5kZ9PBHxGuCrwOXAczPzD4WHDo9H+PIYZc6i6p71woiYm5kbqZKATbS0ErS0EHyEqtVhNvCnVJ8fWsqaIEjSNGCCIKmRiDiIKjn4DvDnmZn19m8DP6e64d23INQfA8/IzAvrr78YETsBR0bE32bmXZl5dkS8Cti3m648LTfvw2MQPl1//fXSGJm5MSKuBPaOiB0zc23L7udQ3cS3upkqYRrP3pm5rq2+/0bVLef1wCfq/f8cEacC13RxDV4IvAm4AHhhZq5ucOxe9fvPRyuQmRkRv6BqbXgUVRehuyPiIuCpEbGg/h6GWwhWZuavIuJmqhaG1gQhgR82qJ8kqU9MECQ1dWj9/pHh5AAgMy+JiH8HDo6I+2dm+41zuwtakoNh51A9qd8VuHSC9VzAtjfvfwCOyMzvdSg/ljvr952o+tEP+wlwTFvZjSUBh5ODiJgB7Ej1VP2XwBqqZKYXHgAE1VP+teOUbbdT/b5mnHLD+xe1bDsX2IeqleD7VAnAOuB/6v0/rLe1ti78srDlRZLUZyYIkpp6BFU3oF912HcpcHBdZrwE4ZoO21bV78u6rt09NgAvqv9/C9WT/Ssyc6iLWMM3y3e2bb8tM3/QXrhEPWj4A1TJwLy23Uu6idnBKVTJx6vqc74qM7cUHjv8vS6i6gI1mk6JxLnA+6laCYYThB9n5uZ6/3nACRGxgKolaSGwsrBekqQ+M0GQ1FT0KM5YU5/24hxbu715v1dFIuYCjwZuauteNJGYT6a6cb6KavDztcB6qm4236J3E0hsBV5bx301MKOelrQkSbgUeEL9Gus6Pp6q1eQ3LdsuqLftX7cQ/CnwsZb951H9+7NPHR8cfyBJ04YJgqSmrgaeCzwWuKRt3x71+7U9PF+jGXv64DXAXOA/ehjzlcBM4PmZOXKt6ifqvWo9ACAzhyLicKpWn9dWp4k/L0gSTqvLv4FREoSIeB7wEOC0eoDy8Dk3RMQFVAnAAVTdvVrXN7iMaqD4AVQJxhCOP5CkacNpTiU1dXr9/p56rnwAImIv4MXAjwrGHzRxVx1/aQ9jFqnXQfgkVf/9j41TvInh1pP2lpL30vnv8l1A199/3a3qSKrB2S8Bvh0Rs8c57N+oZiT6PxHxuvad9XSsX6TqyvXBDsefS/UQ6hiq1pGRRdnqsSvnAwdRtS78ouEAaklSH9mCIKmRemah7wB/DiypByYPT3O6gWoRsl66EHgz8PmI+A9gM/CT1ifvPbAoIl5d//9cqpWUD6DqO38L1WxNncZMdGs58DfAioj4EtW0oAdR9cfvNAXrhcCBEfEu4Dqqe+xvNTlh3ZJwJFWLzOHAdyLi5S3jAtrLZ0S8jGoq05Mi4uVUaymsq+t5JNW/Ia/IzE4Dys+lWh/hmcC5mbmpbf95wAktZSVJ04QJgqTxDD/lbh0z8Cqq6S+PoHrCvo7qhu/9mfm/PT7/v1B1Q/lz4GVUT9iPpLfdmB5CtdYAVE+7V1H1wX878PVeP93OzP+OiJdQDeT9cH3OHwD70bmrzZuAzwHvoxp0DNVYhabnHapbA4ZbFL5b3/iPVv6miHga8NdU1//DVAnUjVRrV/xDZl49yuE/qb+v+dy7e9Gw1m0rG34rkqQ+ipZZCiVpG/UqvKcDr8vMr01xdSRJUp85BkHSeJ5Wv/e6ZUCSJE1DtiBI6igiXgE8CXgL1Qq5TxjnEEmStB0wQZDUUUTcQTWg9WzgbzPz91NcJUmSNAlMEApFxBaqLlntK6lKkiR1aydgKDOndOKYiPgd1crpvbYmMx/eh7jqI2cxKjcDiFnM7scvjyRJug/awmaYHmNCFwGLFu3Uu6qsuXOoZ7E0uUwQyt05i9mL9o+Dp7oekrY30eAf5PQfXGl7sjLPYAubp0XvhEU7zeD2K3bvWbylf3S1ScKAMkGQpKnmTb+kaWII/x5pejRpSZIkSfcSEQsj4rMRcVNErI+In0XEiwuP3T0iTo+INRGxNiJWRMQeo5R9a0RcGREbI+LqiDg6Ytum3ZKYEXG/iPjXiLgqIu6KiDsj4uKIeHNEzJxIPSeTCYIkSZIA2JpDPXv1wHLgVcAxwAuBy4HlEfGCsQ6KiJ2B84FdgcOBVwBLgfMi4iFtZY8BPk21Ov1zgZOAjwAf7TLmPGAD8DHgEODlwA+BfwQ+1209J5tdjCRJkjSt1EnAgcBhmbm83nYusBvwSWDFGIe/E1gCPCkzb6yPvQC4Fngf8MZ627L66xMz8wP1sSsjYgFwdEScmJk3NIlZl391W33+MyIeABwZEW/OzC1NYk4FWxAk6b4uZpS/JG3HkqEevqqldLp2KLAGOGOkdtXc/KcAjxmnG86hwNnDN931sauAM4HDWso9j+qJ/yltx59M9RC9tTtTaczR3AYM1a9exewb/9pLkiQJqAYp9+q/2qKIWD3Wa5Sq7AVcnrlNX6VLWvZvIyLmA7sDl3bYfQmwc921ZzhGApe1FsrM3wDrh8/RMOZwPSIiZkXEkoh4OXAE8A/D3083MSeTCYIk3dflUPlLkibHMuD2Dttvb9nfyRIgCo9dBtydmRs7lL2jpVyTmMOOAjbX+78FnNDSjanbmJPGMQiSJEkCYGtOqFtQJ2syc3GXx45VmfEqWnpsk3M0Kftt4EKqROAA4J0RsSgz3zKBmJPGBEGSJEnTzSo6P0FfWr93evIO1ZP/LDx2FbAgIuZ2aEVY0lKuSUwAMvNW4Nb6y7MjYhVwfER8NTMv7ibmZLKLkSRJkkjo6SDlCT7+vgx4bIf1CPau3zv13Scz1wPX0HmMwt7ArZl5S8s5AtiztVBEPBKYP3yOhjFHc1H9/ugexuwbEwRJ6gdnBpI0gLaSPXtN0HJgMfCitu2vBa7IzMvHOfagiNhleENELK1jndZS7ixgI/CatuMPB7ZQzSbUNOZoDqjfr+phzL6xi5EkSZKmmxXAucBJ9XoF11LduO8DHDxcKCJWAvtlZrQcezzVTf+KiDiO6mb/mPp9ZAG0zFwVER8D3h8Ra+rzPR14F9Wg4uubxoyIdwKPBc4Bfg8sAp4D/BXw/zLzf5rGnAomCJLUD874I2kADU3duNh7ycyMiEOobpQ/StWacDnVwmlnjnEomXlzROxLdQN+KlWPmfOBZ2bmdW3FP0S13sJRwHuAG4EPAh/vMuYvqBZ4O55qfMFG4NfAO2hbSblhPSdVZO9Hq2+XImL1LGYv2j8OHr+wJElSgZV5BlvYPJGZfnoiIlbvtFMsuupXD+xZzEc+9ibuvDOn/HtTc7YgSJIkCejLNKcaQCYIkiRJAsDOkQJnMZIkSZLUwhYESZIkAfRielJtB0wQJElTq8laEM4OJUl9t913MYqIL0XETRFxZ0T8b0T82VTXSZIkabpJYGv27mVbxOC6L7QgfAp4S2ZujIgnA2dHxK6ZuXqK6yVJgmatArY2SFLfbfcJQmb+uuXLLcBc4MHA6impkCRJ0jRlWi2YZl2MIuIhEfGZiPhRRNwVERkR+49SdmFEfLbuPrQ+In4WES8epeznI2I98HPgv6hW4pMkSVKLrUTPXhpc0ypBAB4JvAK4i+pGfizLgVcBxwAvpLrpXx4RL2gvmJlvAhYCBwHfT5ePlqTBlEPlr0ESM8pekjQJplsXox9m5s4AEXEIMFqLwAuAA4HDMnN5ve1cYDfgk8CK9mMycyvwg4h4e0RckZnfa4u5epy6LWr2rUiSJA2WIR+himmWIGQWP/I5FFgDnNFybEbEKcCXImKPzBytG9FMYPeJ1VSSxuFgWjVR+hnoVyuCn0FJLaZVgtDAXsDlHRKKS1r3R8RC4BDgdGADVWJxAPCu9oCZuXisE9YtDLYiSJKk7ZZjBwSDmyAsA67ssP32lv1QTcH7OuBEIICrgFdk5iUdjpWk3unXE9nSJ8g+Ed4+9WtK2H7E9DM4kEwQBIObIMDY628kQGauA541OdWRJEmSBt+gJgiruKeVoNXS+v32DvskSZI0igSGsnctCI53HlyDmiBcBrwkIma0jUPYu36/dArqJEn9N9XdNqa6i8lUn3+Q9KM70n39mkr3EYM6qfJyYDHworbtrwWuGGMGI0mSJI3ChdIE07AFISJeWv/vk+v3/SLifsC6zDyr3rYCOBc4KSKWAdcChwP7AAdPZn0lqSOfdPfHVLdKDJIm12qqr6u/A9PG1oF9dqxemnYJAvDdtq+Prd9/B+wKI2seHAJ8tH4tplpJ+bDMPHOiFYiIlR02L5xoXEmSJGm6m3YJQmbZ6JjMvBN4c/2SJEnSRGRvByk7SnlwTbsEYTrIzP3bt7lQmqRGBqnLxFR3BZnq1YGnes0KIGaU3ZTlUPkdV8yaXVw2t2wuLlsetE9rNgzS75Y0oEwQJEmSBD0fXBzYjDCYTBAkSZIEwNZ0kLJMECRtp2LmzOKyuXVrH2syAAapy8ZUd3HaTq9VaXek+/zvinQfYYIgSZKkaiXlHk5zaueiwWWCIGm71OhJ5/b6BHmq9WN13n78rJo8ae9by1RZ3NLBzNBwQHPptzUNPv+lPwNbO6TumSBIkiQJwBWQBeByeZIkSZLuYQuCJA3SwNd+dNvpl8I6zJgzpzxkg24jOVR2rfrWbacP6xA06eLU6DPQy8Wxao26YzW4rnYd6i9nMRKYIHQUESs7bF442fWQJEmaTEN2MRJ2MZIkSZLUwhaEDjJz//ZtEbEaWDTplZE0mPrVxac07lR3cWoQd2jzlvKYDZR2HWrSvaWRPlzXZjMT9aGLT5OffxPToUucANjqs2NhC4IkSZKkFrYgSJIkCXCQsiomCJK0PS6U1q9uQ1Ot0fdVOItQls0g1FTpLEpN9GvGpX58rnOqJxvaHn+v+8yVlDVsgP5VkCRJktRvtiBIkiQJgK19WBNDg8cWBEmSJEkjbEGQpEHSj5WUp7j/dZN+9VA+dWfxNKcNzj7V/dqbjGtoNl7hvj1eQsOix9Oc2hoxqEwQJEmSBMCQsxgJuxhJkiRJamELQgcRsbLD5oWTXQ9JA6xfXVH60b2iQV0bdRvZWjbPZcwqnI4UyC3lU5KWrtA8Y/784pgUfk/QrNtM6bWaMbv8n+3SmE30q9tQk7gxc05RuX6t0L29cyVlgS0IkiRJklrYgtBBZu7fvi0iVgOLJr0ykiRJk8RpTgUmCJLUHwM000rfuo006DrUD6XdcWJWf7rt9KPbTKPzzyyf8QkK4zbojjZjTpPzlyvtZtav2bH60XVLmm5MECRJkkQCQz3sfd5oCmFNKyYIkiRJAmCr05wKBylLkiRJamELgiT1Q7+mOW0Sty8x+9D/usH332hK1NK+4lvKp8OMfk2JWjolbJMBpE1+rlHWGaRRv/4mY1s2bSov248xAI2uVR9WM59Ghlz9WNiCIEmSJKmFLQiSJEkCHIOgigmCJEmSAFdSVsUEQdLg6Fe//qnuK9xobvnez5ffRLO59Qs1iNlkzYLIwnUQFuxQHDM3lveVjx3KxyuU9vrOTWVrAACNxlYUjxdosg5Gg5/rUB/W1+jXmhXFcYvHKhSfWpo0JgiSJEkCgqGerqTsgOdBZYLQQUSs7LB54WTXQ5IkSZpsJgiSBsdUdwVqok9Tlw41mA6y1Ix5c8sLN+gKUtptpckUl1HYxQrKuyPlho3lMXds8KyoydSd6zeUFZzd4J/tJt3BSqdZXdBgmtcG3aEafQYL65oNetllk891edDex+yzpLdjEOw9NbhMEDrIzP3bt0XEamDRpFdGkiRpkgw5i5FwHQRJkiRNQxGxMCI+GxE3RcT6iPhZRLy48NjdI+L0iFgTEWsjYkVE7DFK2bdGxJURsTEiro6IoyO2bdotiRkRj46IT0XExXW5VRFxfqd6R8SxEZEdXn8ovUb9YguCJEmSANg6vQYWLweeABwNXAscASyPiBdl5orRDoqInYHzgVuAw4EtwDHAeRHx+My8oaXsMcBxwEeAc4Bn1P+/FHh3FzGfAzwfOBX4KTAbeA1wRkT8TWae0KHKBwF3tXzd+76kDZkgSNo+NZkSdao16Kvcl2lGm0wHuUODKUEL+9XPmN+gX3sDWTjNZ5PpSJlZ/rkqPT9ALFxQWLDBdJx3ry8uO2PpkrKY6+4ujtloDMTmBlOyFv5u92XqUk2aiHgBcCBwWGYur7edC+wGfBIYNUEA3gksAZ6UmTfWx15AlWS8D3hjvW1Z/fWJmfmB+tiVEbEAODoiTmy58S+KCXwL+Fxmtg7BWBERu1AlFCd0qO/PMnP1uBdlEg3Qv6CSJEnqp6Gc0bPXBB0KrAHOGN5Q33SfAjxmtO5CLceePXwjXx+7CjgTOKyl3POAeXXMVidTPURv7RZUFDMzb2tLDob9FFgWEf15ItJjJgiSJEnql0URsXqs1yjH7QVcnrlNE+slLfu3Ud+A7w5c2mH3JcDOdXeh4RgJXNZaKDN/A6wfPkfDmJ3qFMABwDWZ2alp71cRsbUea/HlsWJNFrsYSZIkCZhWYxCWAVd22H57y/5OllCt0HZ7h32tx95Sv9+dmZ3mOr6j5RxNYnbyNuBJwOvatl8NvBe4mGrcwZ9Sjbd4dkQ8MTPvGCVe35kgSJIkCejLNKdrMnNxl8eOtZTCeMsslB7b5ByN6xMRhwDHAydn5tfudUDmqW3Fz4mIC4HvA0cBfzfG+frKBEFSuSYDf/u0UFhfzj/FGg08Lh2g2SBmowGapQt6NdGHBb2gfFGzXFa+xE2sb7Co25LFxWWLBx+vvWv8MsMh55YvKpdr7iwr1+T6z59XXrbBAnBDhQvbNaprv35fNBGr6NxKsLR+7/Q0H6on/1l47CpgQUTM7dCKsKSlXJOYIyLihcC3gdOAN4xS33vJzLMj4ibg6SXl+8UEQZIkSQBsnT4LpV0GvCQiZrSNQ9i7fu80HoDMXB8R19B5jMLewK2ZOdwV6DKqrkN7Aj8fLhQRjwTmD5+jYczhGM+nSgzOAl6V2WRtb2YAU/qUa9p8CiRJkqTacmAx8KK27a8FrsjMy8c59qB6alEAImJpHeu0lnJnARup1iloNbzOwZldxCQinluX/wHw8szcPEZd7yUingM8ALiw9Jh+sAVBkiRJJDDUw0HK4w0SGMcK4FzgpHq9gmupbtz3AQ4eLhQRK4H9MrO14sdT3fSviIjjuGdRsy3AR0fql7kqIj4GvD8i1tTnezrwLuCEzLy+acyI2IcqOfg98AngCXHvLoQXD3dnioiLga8DVwCbqRZpeydwFfC5Zpert0wQJJXrV7/+qR4v0I+xFX1aqK20/3OTPt0z5jVY/Kyw/zeU90FvtFBZx+nFRwtcdqPTZFwBGxuUbbCoGqXXdUuDXgoNxiCUjgOZ0eBnNdRgvEQ2WSitD7+DOTTBW9ntyHTpYpSZWQ/w/Wj9WgxcTrVw2pljHEpm3hwR+1Ld1J9K1WPmfOCZmXldW/EPUa23cBTwHuBG4IPAx7uMeSBV96TdgJUdqvcI4Lf1//8aeBPwIKoVl68HvgJ8eKoXTjNBkCRJ0rSTmXcCb65fo5XZf5Ttv6GlpWGM45NqdeMTCsqOGzMzjwWOHS9WXfYVJeWmgglCB3VzVbuy6TAkSZIGUjCUvVwHYdqsqaCGpkc7kiRJkqRpwRaEDjo1V9VLgZdPli1JkjRgtvrsWJggSJoOSgcTTvXia000qGvMarCg1IzCgbezyv+8x4LyQcoxb25xWUrLNhl422DwdS4sG1Abd60vP3+D69qorjstKCoXt4y2NlSHmBvKF7WLhWXnZ0v5YOKY2+Cz0kTxQmX9WdAsh/rw92q6SHrbxcix3wPLNFGSJEnSCFsQJEmSBMCQz46FLQiSJEmSWtiCIEmSJBLY2sMxCA5BGFwmCJKIwlVUS1fxbWyKVyfuh9JrCkDhwGOA3FS2km+j8zcYeNpokG6T1XFLNfgMxMbC8ze4VtlgdeJosOpz8WrOTQafz55dXLZ4her5DVa9brLqdpPPVenfqwbnbzSgeBAHH0sNmSBIkiQJ6PEsRhpYJgiSJEkCYCgHp6VW/eOnQJIkSdIIWxAkSZIEwFbsYiQTBEn0cfBxr/VpJeXS1YmrKhTGbXL+JgOKS8s2iTm/fCXnofstKi4747Y1ZQXnlA/8zfnlZbfsVLaS7+zb1pWff4cGg5TvbLBCc5R9BnOnheUhb19dfv7SAc3rG3xPMxt0UmhQNteV1aHJ71U2+hM4IH8vpQkwQZAkSRLgIGVVTBAkSZIEOEhZFT8FkiRJkkbYgiBJkiQAhhykLGxBkCRJktTCFgSpX0pnsWkyM892qnS2kX7NttQobuHPtcnMSE3E/Pll5XYoKwfAho3l59+4uTxug9mRSg3NKf9na+b6LUXlNu1cPjNQE03+gZ1xx11F5XJh2cxMANxvcXHRuGNtWcFZ5d9V7LRjcVmazLg0r+waRIO6Dt1VPpNV6d/sgZkdrkUCW3s4SDl7FkmTzQShg4hY2WFzf/4FkSRJmhaix4OU7a40qOxiJEmSJGmELQgdZOb+7dsiYjVQvkKQJEnSgHEdBIEtCJIkSZJa2IIgSZIkwGlOVTFBkPple5ydqHRmJmj0/edQ4VwXDc4/Y3b5n7cms42UzrhUOttQY0OF17XB95T3X1JcNtbeXR53UdncDjmrwedqZvnNS6zbVFRuxpyynylAlH5Wgbi77PwAubDs8zJjffksUo2u65ayGZ+YPbs85rryz0ojMwu/r1nlP9ehTeU/qxlz5hSVG8RZjKRhJgiSJEkCHIOgigmCJEmSAHo8zakGlZ8CSZIkSSNsQZAkSRJgFyNVbEGQJEmSNMIWBKlPM/NMuX58X02+/z6cf0aDmYGGNmwsLttkxqPi7yv69BSucBalfMCy4pCxdn1x2dKZiQCG5hVe1ya/Vlk+i9DGXXYsKjd7bYMZbNaVf66GdppXXHbmbWsLY5b/DjSZ8Yj5hXXdXDjbEZB3l3+umsiNZT+DXL+hOGaTvwFRWDa2Fs741ODHNBmc5lRggiBJkiQg6W0Xo/J0XtONXYwkSZIkjbAFQZIkSYCDlFWxBUGSJEnSCFsQJEmSBNiCoIoJgiRJkgATBFVMEKRBmrq0iSn+vqJwOk6A3FpYrsEUizGj/B+53FpYAQDKysaM8ukoY3bhdIgAhVO9bt1hTnHIXDi3uOyMDeVzMpZOsznUoK5b55Z/rqJwStTNO5aff0aD889aV36tcoeyn0FsafB7vXZdednCz+DQA5YUh5zxhwZz2DSYFjjmlNW1yTSrTf4GbF13d3FZaVCZIEiSJAlwHQRVHKQsSZIkaYQtCJIkSYLs8RgEV0obWLYgSJIkSRphC4IkSZIAZzFSxQRBkiRJQPQ4QTDZGFQmCFI06GnXj6lDm5y/Xwq/ryZTlzZRGrfJVIQzF+xQXLZJ3JhV9mczCqcjBcglO5aXnV12rWJz+fe0eWl5XWfO7P0/+NlkStqZ5b8vc26+qyzmnPJ/CptMMxp3lU+zyabCKVELp24FYFaD39fCn8GMOxpMndpg6tImU5JS+Pva5Pc6h+wsL7UyQeggIlZ22LxwsushSZI0WZLedjEy7Rpc0+DRpSRJkqTpwhaEDjJz//ZtEbEaWDTplZEkSZok6SBlYQuCJEmSakNEz17qv4g4OyL+T0TM6WVcEwRJkiRpMD0R+CZwY0ScEBF79yJo112MIuLRwJ7AzlTjUG4FLs3M3/SiYpIkSZpcroMwcHYBDgNeD7wFeEtE/Az4CvAvmVk2nVubRglCRDwW+GvgZcADhjfX71mXuRn4DvDFzPxVN5WSJEmSNLbM3AR8C/hWRDwceB1wOPBF4FMR8R3gq5n5303iFiUIEbE78HHgUGA9cD5wAXA1sIoqSVgKPBJ4GvAGqgzmNOBdmXlNk0pJkiRp8jlIeXBl5u+AD0bEscBBwJuBI4AjIuIKqqThK5k57oImpS0IlwP/W5/ktPECR8QC4KXAW+tj5xWeR7rv6cfia9BoAbbiBdD6tKjcjB0LFwprsqDW+g3FZWPu3PKyCxcUlRvaeXFxzKG5DRbq2lp2XbcsKB+vtmVBgwW1Gtw7zNhUVtetc8s/V3NvK19Qa2he2TWYeVf5Z4Ut5YtvsXFTedkdyz5X2WChuljTYFGz0sUKF5b/cx6rt5SXnd/gNmGo8G/LuruLQ+bm8roW/22bDotgdsEuRtuFxwEvBval+qt9FTAEfBp4d0S8JDN/PFaA0k/vyzPzSZl5aknWkZnrMvOUzHwi8H8KzyFJkiSpoYhYHBFHRcTPgZ9R9eb5HnBgZj46M/cCDgTuBj43Xryix1aZeUa3FZ7IsZIkSZo8djEaLBHxLKoByodS9di5EjgaODkzV7WWzcxzIuLv6VWC0KEy/wg8KjOf17JtFvAnwO8z8w/dxJUkSZJU7AfARuA04EuZed445a8Cxh2w3G0HuWcB1w1/EREBrAQuAm6IiI91GVeSJElTZCijZ6+JioiFEfHZiLgpItZHxM8i4sWFx+4eEadHxJqIWBsRKyJij1HKvjUiroyIjRFxdUQcHbHtIJKSmBHx6Ij4VERcXJdbFRHnj1bvJvUcxTuAB2fmqwqSAzLz3Mw8YLxy3SYID6EatDzsBcAzgG8A/wwcHRF/1mVsSZIkaTnwKuAY4IVUE98sj4gXjHVQROxMNePmrlRTfr6CarbN8yLiIW1lj6EavPst4LnAScBHgI92GfM5wPOB71JN2PMa4AbgjIh4e7f1HMOOwING2xkRe0bEBwpjjeh2obSkGuQw7GDgt8DhmZkRsRR4G/DvXcaXJEnSJMuc6hpU6iTgQOCwzFxebzsX2A34JLBijMPfCSwBnpSZN9bHXgBcC7wPeGO9bVn99YmZOXwTvbKejfPoiDgxM29oEpMq0fhc5r2u5IqI2IUq0TmhaT3H8UGqbkOXjrJ/r7rMhwpijei2BeF3wB4w0r3oRcAZLRfje8P7JUmSNBiGiJ69JuhQYA0wMtlNfZ95CvCYcbrhHAqcPXzTXR+7CjiTatXhYc+jGth7StvxJ1M9RG/tFlQUMzNva0sOhv0UWBYR87uo51jGu9DzgAbz+Fa6TRC+AxwZEc8G3gPszL0zuSGqJhJJkiSpqb2AyzO3WXjikpb926hvwHen8xP1S4Cd6649wzESuKy1UGb+hmph4L26iNmpTgEcAFyTmesnGjMidoqIh0XEw+pNy4a/bns9jqqL1vWj1W003XYx+izwcuD7VJnLL6hGUQ97JHBrl7F1X9NkQa8mC2VtbbCgUXkFystuh4vpNLn+0GDxrdI27QZrys3YcWF54Z0KF2oDhnYsW9ApZ5X/XLfOK/9TvPZhZYt/RYNrNe+O8t+V2NL7/gfzr141fqHapgcvLi47a93monKNFh/bUH6tNu/2gOKys29aU3b+zQ2uf4O+Ijm/8HO1tsHiY4sb/A42qGvcsbYs5KYGC9U1Ufo3u1+LYPZZH6Y5XRQRq8c+Zy7usHkZ1ZSd7W5v2d/JEqp709s77Gs99pb6/e7M3Nih7B0t52gSs5O3AU8CXtdlPdv9DTDcJSqpui2dMMq5g2ra00a6ShAyc21EPJlq7MFi4F/bmlMOBn7ZTWxJkiSJ6ua3m31Njm1yjsb1iYhDgOOp1iX4Wulx4+xbORyeKlFYzj0tK63H3gVcON6qyZ1024JAZm6iGqF9LxFxP+BH3LtFQZIkSdNYQk+mJ22NB6wZpYVgPKvo3Eow3IW905N3qJ78Z+Gxq4AFETG3QyvCkpZyTWKOiIgXAt+mWqPgDROo573U05meV5/j4cA/ZeZPOpXtVnGCEBHXA6fXr5WZ2bGdNTNvA47oQd0kSZI0iabLLEZU4wJeEhEz2sYh7F2/d5y1JzPXR8Q1dB6jsDdwa2YOd9u5jOop/J7Az4cLRcQjgfnD52gYczjG86kSg7OAV7XfN3cTc5Tv98jxynSjSefnfwMOAc4Gbo2IUyPi0IjYoR8VkyRJ0n3Wcqpu7C9q2/5a4IrMvHycYw+qpxYFoJ6C/0VUN+3DzqJahfg1bccfTjXzz5ldxCQinluX/wHw8swcbUBUccyW/a2Dk+kwMLnja5Tzj6q4BSEzjwKOioinUE3LdAjVyOgNEXF2/U2eWU/PJEmSpAHTh0HK3VoBnAucVK9XcC3Vjfs+VGNdAYiIlcB+ee+KH091078iIo6jutk/pn4fWQAtM1dFxMeA90fEmvp8TwfeBZyQmdc3jRkR+1DdE/8e+ATwhGoSoxEXt3RnKorZ5rfAUETsUHf3/y3jj8eARjOHdDEGITMvAi4C3hMRj+GeZOEkqgr/iOrCnJGZv2saX5IkSfdt9cK7h1DdKH+UqjXhcqqF084c41Ay8+aI2JfqBvxUqh4z5wPPzMzr2op/iGq9haOopu6/kWphsY93GfNAqu5Ju3HPYOJWj6C6qW9az9b6JvesbTD8dU9F57UcuggU8WDuSRaeSZWp/BJ4b2b+Z09OMoUiYvUsZi/aPw4ev7CameppTvs1zegATXNael1jZqMHEOXnnz9//EIAUf7zjzmzyyvQj2lOZ5dfqy0LyqaYhKmf5nTWugZTog6V/fsy54Y7imP2Y5rT2FA+HWas6zQbYmebH1q+HFDpNKeNOohvKK9r6ZSkcdf68pg7NuiB3IdpToduK+/QkIWfVWjwb0vhvwEr8wy2sLnbgbw9ExGrZyyYu2iPb76rZzEvf+XHGVq3ccq/NzXXszuTzPx9Zp6YmQdSLZx2JFWG1HEhC0mSJE0jWc1i1KtX759ra7L05dFlZq7OzK9n5mGZeXw/ziFJkiTdl0XEIyPieW3bnhoRZ0bEf0fEX3YTt+t1EKQx9WPFYaDz5LqjVKGwO0xfVlzu0/ljVnm3mUarHhf+vGLe3AYxG3QHmlvYxWZL+bXKncu7dzRacbbw57rx/oXdphravLDsus7cUB5z3QPK/ymYu6b8d3vB78u6o2y9X3kXrznXjTb1+bZyx7KfQWzcMn6h4ZgLyn8HZt5d1sUJKP5s58KyLm4AQ0vKVzKeUdh1av0fla8O3WSF7JxX/rct15Wt5hwLF5THvPOu4rLbu2k0zanKfJxq3YT/hJH1yM4CFgLrgS9ExC2ZeXqToGP+pY+IWRHxvoj4WkQ8ISIeFBFnR8S1EfGFiCj/SyVJkiSpl57EvRcnfgWwE/AE4P7AT4C3NQ063qOgTwDPBh5GlZm8CfgS8F6qgcjHNj2hJEmSpqfM6NlLk+L+VDMvDXse8N+ZeWk9Deq3gD2aBh2vXfllwOOoEombgVMz8wqAiLiUajrTdzc9qSRJkqQJW0c1BSwRMZNqnYjPtuxfT9Wi0Mh4CcJOwwufRcRdw8kBQGb+b0Ts3PSEkiRJmo56/eTfVoRJcBnwmoj4OtWD/YXA2S37Hw7c2jToeAnCmoiYn5nrgY+07oiIRUD55NGSJEma1hyjPHD+ATgDuKX++mKqxdaGPQf4edOg4yUIZ1JlHr/OzI+37Tu0mxNKkiRJmrjM/I+IeBZwMNWK0CdmvQpyRCwDbgC+3jTumAlCZh41xu7vAt9pekJJkiRNTw4uHjyZ+UPghx22rwIO6yZm1+sgZOa6bo+VJEmSND11nSDUzRZ/DGwBrsjMW8Y5ZGBExMoOm8tXnJEkSRpEDkIYOBHxdODNwKOAZWw7Ojwzc/cmMbtKECLiYOCfgR1att1INSZh5JWZv+8mviRJkiafXYwGS0S8FvgasBm4EriuF3G7bUH4OLABeB9wJ/BYqhXb/hR4UV0mgZkTreBUyMz927dFxGpg0aRXZlDlUF/Czpgzp7js0OYtfalDr83cccfiskPrNxSXjfnzy8vOKvxTMKvBr3Q2eAw1Y7w1G2v3L/8V3LxsQXHZnFH+D+LQ7LKyQ7PKY27cqfD7B+avKvvdmrG5OCTrHlB+/tnry3+3Z91yZ1G53GFuccycX1421hdOtLdla4OY5b+D+YAl5WUXlH1fQ/PK/wauf1D534A5q8v+Bsy9bX1xzE0PXlxcNhr8uZi9vvC6NvhZzVhYXoHcuLGs4FBhzPJqSp28D7gCODAzbxyvcKluE4SHAu/PzM+274iIhwNPpFpgTZIkSQOiybMdTQsPB/5vL5MD6D5B+BXV6srbyMzfAb8DTuu2UpIkSZLGdQNQ3rRaqLxd+d4+A7yqXtJZkiRJ24HM6NlLk+Kf6MM9eVctCJl5akQ8BTgtIl5Xz7MqSZKkQeaN/aD5H+AlwEUR8TngWmCbgVX1WgnFup3F6MHAw6gGJN8YEecCF1LNXvQ/zl4kSZIk9d1/tfz/V9h2otqgi4mDuh2D8DXgQOAXwB3A44HnDFcqIm6jShRe0GV8SZIkTTIHKQ+cI/sRtNsE4RnAlzPzr4Y3RMRDqKY6fWL9esLEqydJkiSpk8w8pR9xu00QVlN1JxqRmTdQjaT+twnWSZIkSVPBFgTR/SxGpwH79LIikiRJmjpJb2cxMteYHBHx0Ij4akTcEBGbIuJZ9fb719uf3DRmtwnCPwJ/EhH/p8vjJUmSJE1ARDwC+BnVTEaX0TIYOTNvBZ4EvKFp3G67GF0B3A18s04SvgNcmJm/7TKeJEmSppqP/QfNR4AhYC9gPXBL2/4VVLOONtJtgvAPwOOoZi86pH5lRKymGpswPN3pd7qML0mSJGlsBwL/mJnXR8SyDvt/BzykadBuF0p71/D/12siPIEqWRh+PZsqBzVB2N5Et73SeiO3brP2xxiFh4qKzZgzpzjk0OYt5ecvlFvKY858eIPf8ejDYjezyqdRzgbnz/mzi8pteOCC4pibdiz/rN71oN5/rrc2WPh+1t3lZTcuLbuuTWLOa7DUZTa4VFuXLCwr2CDmzFV3FZfNOWWfq6HF5Z+rnFP+OzDr9nXFZYd2KPs7FA3moJy9tvzv5VDh97Vlx/IP9salZdcfYN5tm4rLbtllcVG5mXeU/xLEps3FZZlVdus0dHeDX8JpxBWQB85OwE1j7J9DF/f73bYgjKgXRfs9cObwtohYQpUoSJIkSeqP64E9x9j/NOCqpkH78jg4M+/IzHP6EVuSJEl9kj18aTKcBrwuIvZq2Ta8cPFLgJfRRY+eogQhIp7dNHDLsQd2e6wkSZImU/TwpUnwEap1yH4C/DNVcvDuiLiAKjH4JfDJpkFLWxD+MyLOiYg/i4hxOypGxOyIODQizqMaPS1JkiSphzLzTuDpwFeopjQN4CDgj4DPAwdk5oamcUvHIDwe+BTVKsm3RcTZwEXA1cDtdWWWAo+i6uv0LGAJ8H2q2Y4kSZI0nfW6a5DdjCZFnSS8DXhbRNyf6r781swGsxq0KUoQMvNS4DkR8XTgTcDBwCvY9kcfwJ1U/aG+kJk/7bZikiRJksrVi6NNWKNZjDLzAuCCupvRE4E9gPtTJQq3ApcCF2cWzi8pSZKk6cOn/tNaRDyzm+My84dNyne7DsJWqi5GF3VzvCRJkqYh10GY7lbSXRpXvogLPVgHQZIkSdKkOLLt6wDeAjwa+AZweb1tD6rhAFcCJzY9iQmCJEmSAOh+WKsmQ2ae0vp1RLyVqrv/H2XmjW37PgxcAOzY9Dx9WShNkiRJUt+9Bfhie3IAkJk3AF+syzRiC4IkSZIqtiAMmocCd4+xf11dphFbECRJklTJ6N1Lk+Fa4NURMa99R73ttcBvmwa1BUGSJEkaTJ8EvgT8NCI+B1xB1Q70WOCo+v2vmwY1QeiTmFk+m1QOlbfnxYyyjDy3bi2O2Q+l9eynGXO3SaYnHnN2+a/MjCWLi8rl0rJyAFuW7FB+/o1bisuuf2BZ3Bmbyz+r2aB9cvUjZxeV27CsPObWueVl5/zRmuKyT3zgDUXl5s/cXBzz0Qv+UFz2xIueVVRu6y1l1xRg7h3lv6+bdyj/wW59WNnnas7a8r9XG3aeX1x23q0bisrFlvKle2JD+c/19iferzxu4a/WvFXlv9ebFzb4d6jwxzpzY/m1avL3IraWl51R+DOIzeXXijnlvy9bb7+jPO4AKv0sanrIzK9ExELg74DPc08nsQDWA/83M7/cNG5XCUJEXAr8F3AusDIzV3cTR5IkSVL3MvOEiPga8BxgN6rk4Grg7G7v0bttQVgHvIlqVPTWiPglcE79Oj8z13UZV5IkSVPFFoSBlJlrgO/2Kl5Xg5Qz86nAUuBg4B+pVmd7B7ACuD0ifhQRx/WqkpIkSZImR9djEDJzLXBm/SIilgLPB94DPAN4OvDBHtRRkiRJk8HZh8QEBylHxAzgycCzgGdTJQbzgD9QdTeSJEnSoLCLkeh+kPJbqRKC/YCdgDuA84D/C5yTmb/qWQ0lSZIkTZpuWxBOALYC3wQ+A1ycmeackiRJg8y7OdH9Sso/ADYCrwH+A/jniHh9RDyiZzWTJEmSNKqIKF8cpoFuZzF6DrCEqovRF4GHAJ8DroqIayPipIh4Ze+qKUmSpL7LHr40GW6KiC9ExBN7GbTbFgQyc3Nmnp+Zx2bmflQJwyuAu4EjgFN7U0VJkiRNiozevTQZfgy8AbgoIn4REW+OiMUTDdp1ggBVs0ZEHBQRHwNWAt8AHkuVN1480cpJkiRJ6iwzXwA8HPgAsAD4LHBjRHwjIg7oNm5XCUJEfCAizqOaveg/gXcBC4F/Al4C3C8zn9RtpSRJkjT5Inv30uTIzBsz8yOZ+SiqWUZPAw4BfhARV0fEeyPiQU1idtuCcCzVuINTgVcDD8zMPTPzLZm5PDNXdxlXkiRJUhcy89zMfDXwIKqePY8APgz8NiKWR8RTSuJ0O83prpl5XZfHSpIkaTryyf9Ai4j7Uc0y+npgD2Ad8G2q2UdfDbwoIt6YmV8eK063sxiZHEiSJElTLCrPi4jvAjcAn6RKCN4EPCgz35CZRwEPoxoz/P7xYnY9SDkiFkTEcRFxSUTcVb8uiYhjI2JBt3ElSZIkjS8iPgT8jmpdsucCpwBPzswnZuY/Zeba4bKZuabe/+Dx4nbVxSgilgLnU81YdBv3zFj0aKpR1C+LiH0z8/Zu4kuSJGlyBb0dXBzYY2kSHAP8D9U4g29m5rpxyv8c+NB4Qbsdg/Ah4DHAm4EvZuZWgIiYCfwl8I9UA5nf2mX8gZdDDX4lcqi4aMyaV1Ruxg47lJ9+0+bisswom9c4ZjX4aGWDazVndnHRKL0GW7eWn39+2fUHGFpc1pB290MXFsfcuKi80W+owY9g7cP6MF91g5CbH7m+qNz8HTYVx9zvwdcWl73/nLXjF6p9eOdLi8qdvq68IfVfbn5acdl5O24sKrdhY/lnZcsO5R+WrXPL484o/NOycXGDus4r/2Bt2qnsZ7BxUXnMoZnFRZm9rvxv28bFZXVY++Dyv4Ebdi4//+y1Zedf+qvikI3uDGfNLb+wObOsrjmj/HM184415WUXlP3bsnXd3cUx1VlELAQ+CrwMWAxcBnwoM/+t4NjdqbraHEDVY+Z84J2ZeXmHsm+luqd9OFUXnS8Cx2fe+watNGZEfAR4AvBE4P7AcZl5bIfzHgt8sEP1b87MXcb7HmtPyMxfFJYlMy+juo5j6raL0YuBr2Tm54eTg/qkWzPzC8BXqaZXkiRJ0iBIertQ2sSbD5YDr6J6Sv5C4HJgeUS8YKyDImJnqpv3XYHDqRbyXQqcFxEPaSt7DPBp4FtUXXROAj5ClZh0FRN4G7ATcHrh93kQ8PSW15jfX5tPRcSzR9sZEQdExDkN4gHdtyA8gLEXQvs51cWTJEmSGqmTgAOBwzJzeb3tXGA3qqf4K8Y4/J3AEuBJmXljfewFwLXA+4A31tuW1V+fmJkfqI9dWY+lPToiTszMG5rErO2UmUP1isZ/UfDt/mwCSwTsD3xljP07A/s1DdptC8LNwOPH2P/4uowkSZIGRfbwNTGHAmuAM0aqlplUg2wfExF7jHPs2cM38vWxq4AzgcNayj0PmFfHbHUy1UP0F3cRk/auSVNsMdWMRo1024JwJvBXEfFz4MvDFyIiZgBvAF5H1X9LkiRJg6L3o4oXRcTqMU+ZubjD5r2AyzvcbF/Sur/9oIiYD+wOfLdDzEuAV0bEzpl5Sx0jaeuTn5m/iYj19f6mMbvxq7oL0y3AvwPvGytWRPwx8LiWTftGRKd7+qVUU51uc53G022C8AGq/lKfB46LiCvq7X9ENRjjKjoPupAkSZLGswy4ssP221v2d7KEapqMTjNpth57S/1+d2Z2esJ+R8s5msRs4mrgvVTd9jcBfwocDTw7Ip6YmXeMctyh3HOfncBf1a9O1tLFpEFdJQiZuSoingS8i2ow8pPrXddQ9YP6RGbe2U1sSZIkTY1eTnNaWzNKC0GJsWozXk1Lj21yjonUZ9sDMk9t23RORFwIfB84Cvi7UQ49mWrBswDOoRpQfXaH+txF1QqzoWndum1BoE4A3le/JEmSpF5ZRedWgqX1+2hrbd1BdXNccuwqYEFEzO3QirCkpVyTmBOSmWdHxE1UsxmNVuZ3VIujERFHAj/MzPK5vQt0vZLyoIiIuRHxtYi4ISLWRMS5EbHnVNdLkiRp2pk+g5QvAx5bj29ttXf93nFRmsxcT9WjZa8Ou/cGbm3p338Z1VP4e90XRsQjgfnD52gYsxdmAEUDnTPzlF4nB1DYghARr+0meGZ+vZvjemwW1Q/1acBNVHPTng48agrrJEmSNP1Mn6WPlwOvB15Ey0xGwGuBKzoteNZ27JsjYpfM/ANARCytY/1LS7mzqGb4eQ3VFP3DDge2UE3K0zTmhETEc6iWE7hwlP3D9+SnZmaW3qM3vScv7WJ0MtVHpslyqwlMeYJQLzn94eGvI+JE4PiIWFZPTyVJkqTpZQVwLnBSvV7BtVQ37vsABw8XioiVwH6Z2XqPejzVTf+KiDiO6mb/mPp9ZAG0ekztx4D3R8Sa+nxPpxpje0JmXt80Zl2n/agm7RlednuPiHjp8PeVmXfX5S6mule+AtgMPINqvYWrgM+Ncl1OprrH/hbVwObhr8e6R298T16aIBzQJOhE1KvR/V+q5akfBywADsjMlR3KdrME9zOAW0wOJEmS7q0Pg5S7Uj8dP4TqPu+jVPd5l1MtnHbmGIeSmTdHxL5UN/WnUnXZOR94ZmZe11b8Q1TrLRwFvAe4kWqGoI9PIOZx3HtxspfVL4BHAL+t///XVNOQPgiYDVxPNdnPh8dYOO2Auj6bWr/utdIE4XDgi5n5E4CIeCbwq8y8tQ91eiTV8tU/B/6Ley9S0W458ASqKaGuBY6gWoL7RZm5zQp79Yp2X6SaUkqSJEnTVD0hzpvr12hl9h9l+29oaWkY4/gETqhf45UtjdmxTh3KvaKkXNsx5431da+UJghHAD8AflJ/fS5VM8s3+1CnH2bmzgB15tgxQWi6BHdEzKMae/DvmfnVDvFWj1OvRU2+CUmSpIGTTXqTa3tVmiDcRjVgYljfPj0NlqfuuAR3RJwCfCki9hgewBIRM6n6av2eqm+XJEmS2k2TLkbqrO7F01hm/rBJ+dIE4cfAMRHxMKq5YAEOq6eBGqMu+eEx9k9UkyW4v0w1XdXL6qakbYy3iEfdwmArgiRJkqbKSpqlcVGXn9nkJKUJwtuBU6iWah4+0WH1azRJy+xBfVC0BHdEPBw4EtgA3BEx0vjx/Mw8v4/1kyRJGijTZZCyRnXkZJykKEHIzN8C+0XEHGAXqtHXb+fe89JOhXGXvK5Xm7NDnSRJkgZaZp4yGecpbUEARqZUuq7u5/+T+uZ7qnS7BLckSZI6sQVBNEwQhmXmpDRvjOMy4CURMaNtHMKYS3BLkiRJg6geD8zw2gvDX4+nw1oNY+oqQZgmJrIEtyRJkto4BmHa+y0wFBE71D17fktZu09fBilvIyJeSbXq3KPo3NUnM7Or+C3LUT+5ft8vIu4HrMvMs+ptRUtwS5IkqZAJwnT3Iaqf0pa2r3uq2xv4Y6iWkb6ZagrUO8Y+orHvtn19bP3+O2BXmNgS3JIkSdKgycxjx/q6V7ptQXgT1Tysz8vMzb2rTiWzbBm/kiW4uxERKztsXggQs2YXxcitW8vPVxizUdwt5RM3NanrzKWdGos6aBCTKK/r0IPuVx53c1kdtizZoTjklh3Kf2U2LyxrzbvpGQ2+/6UNft3uKq/rjCUbiso9cdfri2Nev3ZxcdlbVy8sKrfvQ64ujvnQeeXPLZ6yQ3nc/7h7XlG5T17znOKYt61dUFx2401lZWcu21gcc+0jZhSXbWLWurK4c1eV/w5sKb9UrH1EWbnZa8tjbt6xP49XZ99VVm7j0tK1RMuvP1D8/HFNg8/K7HXlp59zZ/lnYPNOc4rK7XBd+e/A0N13l5fdUB53INmCILpPEHYCvtOP5ECSJElSuYh4CnAosFu96Rrg9Mz8STfxuk0QLgYe2uWx015m7t++zZWUJUnS9s5ByoMlImYCXwKOYNt1v46OiK8Db8jMBl07oNt25WOAv46IJ3R5vCRJkqSJOYZqdeUzgGdQjcldDPwp8G9Us3se0zRot+sgnBcRrwcujIgLqKZYas9MMjNf3018SZIkSeN6HXB2Zh7Wtv0C4NCIOLsuc1yToN3OYvRU4OT6+H3rV7ukWqdAkiRJg8AuRoNmZ+ATY+w/HTi+adBuuxh9BthMtd7A0syc0eHVaEEGSZIkSY1cCewyxv4H1mUa6TZB+GPg+Mw8MzNXdxlDkiRJ00VWg5R79bI1YlJ8DDgqIv6kfUdEPJ5qaYKPNg3a7SxGtwCbujxWkiRJ05E39dNaRHygw+ZrgJ9FxPeBX1P9FPcADgJ+CTy66Xm6TRC+Crw6Ik7MzC3jlpYkSZI0UceOse/59avVE4DHAx9ucpJuE4QfAX9GNYvR54Fr2XYWIzLzh13Gn1JjraQsSZK03bIFYborXCN+YrpNEH7Q8v9fYduPU9TbHKgsSZIk9UBm/m4yztNtgnBkT2sxzbiSsiRJui9yJWVB9wulndLrikiSJElqJiJmAYcATwWWsO0spY0XL+62BUGSJEnbG1sQBkpELAXOBfbini7+Ue/Olm2TlyBExJMYO1tpNGJakiRJU8cuRgPn74DHAG8AVgJXA88FrgPeDzyq/rqRrhKEiJgPnAY8h7GzFRMESZIkqT9eCHw9M78WEcvqbVsz8wqqJQlWUi2m9sYmQbtdSfkDVMnBR4ADqBKCw6nmXj0f+CnVAg2SJEkaFNnDlybDLlT33QDDa5PNa9l/OvDipkG7TRBeCnw3Mz8AXFpv+31mfg84EJgDHNFlbEmSJEnjux1YUP//WmAz8NCW/ZuphgI00m2C8FDgvPr/hxdImwNQr6z8L8CfdxlbkiRJU8EWhEFzJXWvncwcAi4GjoiIuRGxA/Ba4JqmQbtNENZyz/iFtcAQ8KCW/WuomjwkSZI0ICJ799Kk+D7w0oiYW3/9KaoJhG4HbgGeBHy6adBuE4SrgUcDZOZW4DKqbkdERACHAdd3GVuSJEnS+D4K7JKZGwEy8ztU9+TfA84CXpmZJzUN2u00pz8AXhcRb68ThC8CJ0bE1VSNSo8A3ttl7ClXj/hut3Cy6yFJkjSpfPI/UDIzgY1t206jmm20a90mCH8PnEo9tWlmfj4i5gGvphqT8GXgExOpmCRJkqRy9VIEZOb6icTpKkHIzLuAK9q2fYqq39PAy8z927dFxGpg0aRXRpIkabLYgjBwImJn4FjgEOAB9bZbgOXAcZl5c9OYE1pJua7AXOB+wK2ZuWmi8SRJkjQ1HFw8WCLiEcCPgAdSPby/kKqHz2OAvwYOjoh9M7PRTEbdDlImIp4QEedQzWJ0HbBPvX3niPiviDiw29iSJEmSxvVJYBlwWGY+NjMPy8xDM/OxwEvqfcc3DdpVghARj6NaMXl34Out+zLzFmA+1crKkiRJGhSugzBong18LjNPb9+RmcuBL9RlGum2BeFDwI3AnsC7qQcrt/gv4CldxpYkSZI0vgR+M8b+K+kiXes2QdgX+HI9WLnTSa/j3gunSZIkaZpzobSBcx5wwBj79wdWNg3abYIwj2q15NHs1GVcSZIkSWXeDjw1Ij5Zz2YEjIwJHl5V+e1Ng3Y7i9HVwBPH2P8s4PIuY09ruXVr72Nu2dzzmDFzZl/K5tq7ymLutGN5zLvLp+qdcV35TF25y7KymOu3FMfcumR2cVmy7PHJwuvL8/Q75zb4Wc0dKi4758r5ReUuv+TRxTFnrysuyg5zysqdc/UTimNG+Y+VkxaP9fDl3ubcUfbzuv/F5RV44Ibyn9WWBWWfq9lryz8ra3Yr/wzOvrv8seDWOWVlN+3Y3kt1dPNWFRdl4Q1l5TY2eKQ1NLu8rlt2LP+5ZpTFXXJp+c9qqMGfqzlry35W81aX/xs4Z3X5v22zVpf/OxCrxno+eY9cd3dxTOYU/hECWD+h6eWnP5/8T2sR0Wk2oh2okoC319PyJ7Ck3reKquv/7k3O022C8E3g/RHxHeDielsCRMQ7gOcBb+sytiRJkqaCCcJ0dx2T8FPqNkE4HjgI+B7wa6qKfjoi7g/sApwNfL4nNZQkSZLUcTHffuhqDEK9INpBwDuB9cAG4NHAbcDRwJ9lZnnbqiRJkqZc9PClwdW4BSEi5gMvA67IzE8Dn+55rSRJkiQViYjdgYOB3epN1wBnZObV3cTrpovRRuDLVGMMftLNSae7iFjZYfPCya6HJEnSpOn1AmeOZ5gUEfFhqnXJ2men+EREfDQzP9A0ZuMuRnXXoetxKlNJkqTtiusgDJaIeB3wPqqH9ocCj6pfhwAXAO+LiCObxu12kPIpwGsi4jOZubHLGNNWpwEg9bRRiya9MpIkSVJnR1ElB/tnZuv82ldHxArgfODNwNeaBO02QfgxcBjwi4j4PNUSz9tMOJyZP+wyviRJkiabT/4HzWOB97QlBwBk5paI+BbwsaZBu00Qzm75/8+w7ccp6m3lK/VIkiRJamITY4+T3bEu00i3CcLrMMeUJEnavnh3N2h+CvxVRHwlM29u3REROwN/SReTCnWVIGTmyd0cJ0mSpOnLwcUD58PAfwG/ioiTgMvr7XsCR1K1ILyqadCuFkqLiK9GxFPH2P+UiPhqN7ElSZIkja8e73sYsBZ4B3BS/frbetthmXl+07jddjE6AvgBozdZPAI4nKorkiRJkgaBLQgDJzPPjIj/AJ5IdQ8ewNXAz+vlCRrrNkEYzwJgc59iS5IkSfdpEbEQ+CXwj5l5AtV4hJ/2InZxghARDwN2bdn0mIh4ZoeiS4E3AldNrGqSJEmaTI5BGByZeVdELAPu6nXsJi0IRwIf5J6FuN9Xv9oFMFSXlyRJktQfFwJPAr7Sy6BNBimfTnXT/3qqJODLVGMMWl9HAi8FHpGZp/ayopIkSeqz7OFrgiJiYUR8NiJuioj1EfGziHhx4bG7R8TpEbEmItZGxIqI2GOUsm+NiCsjYmNEXB0RR0fENvfIpTEj4iMRcVZE3BIRGRHH9qKeo3g38PKIODIiosFxYypuQcjMX1L1cyIiHg78v8y8tFcVkSRJ0tQJetvFqAd3q8uBJwBHA9dSTZKzPCJelJkrRj1vNf//+cAtVJPmbAGOAc6LiMdn5g0tZY8BjgM+ApwDPKP+/6VUN9+NYwJvo7pnPh34i17UcwyfAu6gakH4RERcDdzdViYz89kFsUZ0uw7Cca1fR8Qs4CnAg4HLM/OybuJKkiRJEfEC4ECqaTqX19vOBXYDPgmMmiAA7wSWAE/KzBvrYy+gSjLeRzVWlrr//vuAEzPzA/WxKyNiAXB0RJzYcpNeFLO2U2YORcRixkgQGsYczW5U7TXX1V8/oOCYcRV3MYqI/etmnl3atj8C+B+qDOhbwCWugSBJkjSApk8Xo0OBNcAZI1XLTOAUqolyxuqGcyhw9vBNd33sKuBMqjUDhj0PmFfHbHUy1UP01u5MpTFpMLVocczRZOaumfmI8V6F9RnRZAzCEcDBmfmHtu2nAHsDPwY+TbWC2+ERcXjTykwXEbGy/QUsnOp6SZIkDZhFEbF6rNcox+1F1Sul/Wb7kpb924iI+cDuQKdu8JcAO9dde4ZjJHCvni+Z+Rtg/fA5GsYs0ouYEXH/iHhqROze5NwlmiQIT6bKaEZExGOAfYAfZua+mflOqq5GvwFe27NaSpIkqf+mTwvCMuD2Dttvb9nfyRKq4Q8lxy4D7s7MjR3K3tFSrknMUl3HjIgZEfFPwE1UD+ivjIgfRcT9G9ZhVE3GIDwQuLJt2/5UH4GRqZUyc31EfBN4y4RrN0Uyc//2bXWGu2jSKyNJkjRJ+rAOwprMXNzlsWPVZryalh7b5BwTqU83x422783AXwI3AhcAj6IaXP1FCrsmjadJgjCXqrml1ZPr9/Patl+PN9OSJEnqzio6P0FfWr93evIO1ZP/LDx2FbAgIuZ2aEVY0lKuScxSE4n5WuBXwNMycy1ARHwZOCIiFmfm6oZ12UaTLkbXAXu2bdsHuCUzr2/bvgOwegL1kiRJ0mSbPl2MLgMe22E9gr3r945T7WfmeuAaOo9R2Bu4NTNvaTlH0HZ/GxGPBOYPn6NhzCITjPlHwMnDyUHtH4GZwKOb1GM0TRKE84HXRsTeABFxKFWTxlkdyu4N/H7i1ZMkSdJ90HJgMfCitu2vBa7IzMvHOfag1pk3I2JpHeu0lnJnARuB17QdP7wmQevY29KYTXQbcwFV96JWN7bsm7AmXYw+BrwK+EVEDDf7bKKai3ZERMykmhbq//WigpIkSZockb0fhNClFcC5wEn1egXXUt247wMcPFyonmlyv8xsXZfteKqb/hURcRz3LEC2BfjocKHMXBURHwPeHxFr6vM9HXgXcEJbD5mimHWd9gPuT9WjBmCPiHjp8PeVmcMLmRXH7GC08RE9WU25yUrK19bf8AeBRwIXAX/XYVG0A6j6dJ2BJEmSBsc0yQ8yMyPiEKob5Y9StSZcTrVw2pljHEpm3hwR+1LdgJ9K1WPmfOCZmXldW/EPUa23cBTwHqon8R8EPj6BmMcB+7V8/bL6BfAI4LddxGz3gra1yXag+um9LCIe11Y2M/PT48S7l0YrKWfmz9i2qae9zA+4p3+YJEmS1Fhm3kk1Y8+bxyiz/yjbf0NLS8MYxydwQv0ar2xpzI51mkjMDl5Zv9r9VafTUK1VVqxRgiBJkqTtVx+mOVXvHdDvE5ggSJIkSQMiM9uXF+g5EwRJkiT1anrSe8fTQDJBkCRJEmAXI1WarIMgSZIkaTtnC4IkSZIqtiAIWxAkSZIktbAFQZIkSYBjEFSxBUGSJEnSCFsQJEmSVLEFQZggSJIkqWYXI4EJQkcRsbLD5oWTXQ9JkiRpspkgSJIkqZI2IcgEoaPM3L99W0SsBhZNemUkSZKkSWSCIEmSJMAxCKqYIEiSJKligiBcB0GSJElSC1sQJEmSBEAMTXUNNB3YgiBJkiRphC0ITeVgpNZDGzf0J/DGwnJ33dWf80eDnHbVqrKQM2cWh9zhf2cXl82tW4vKLZxd/mu4y6ZNxWWbfF+xww5F5XJj6QcAKPz+qwr04VlFg9/VmD+/uOzQurvLYs6I4pilnxWAOYXlmvz8d/7J3OKyQ+vL/7ZEg892sX58rpp8VuaU/gQgG/y+Dm3eUly2vAKD8e8VQIOfqvrNMQjCBEGSJEk1ZzES2MVIkiRJUgtbECRJklRxJWVhC4IkSZKkFrYgSJIkCXAMgiq2IEiSJEkaYQuCJEmSKrYgCBMESZIkAWSPuxiZbAwsuxhJkiRJGmELgiRJkipOcypsQZAkSZLUwhYESZIkAU5zqooJgiRJkiomCMIuRpIkSZJa2IIgSZIkwC5GqtiCIEmSJGmELQgdRMTKDpsXTnY9JEmSJtWQTQgyQZAkSdIw8wNhgtBRZu7fvi0iVgOLJr0ykiRJ0iQyQZAkSRLgIGVVHKQsSZIkaYQtCJIkSaqkTQiyBUGSJElSC1sQJEmSBDgGQRUTBEmSJFVMEIRdjCRJkiS1sAVBkiRJBBA9HKQcPYukyWYLgiRJkqQRtiBIkiSpMjTVFdB0YIIgSZIkoLddjDS47GIkSZIkaYQtCJIkSaqmOO1lA4KNEQPLFgRJkiRJI2xBkCRJUsUxCMIEQZIkSbUwPxB2MZIkSZLUwhYESZIkVexiJGxBkCRJktTCBEGSJEkAxFDvXhOuS8TCiPhsRNwUEesj4mcR8eLCY3ePiNMjYk1ErI2IFRGxxyhl3xoRV0bExoi4OiKOjoht7pF7HTMijo2I7PD6Q8n32E92MZIkSdJ0tBx4AnA0cC1wBLA8Il6UmStGOygidgbOB24BDge2AMcA50XE4zPzhpayxwDHAR8BzgGeUf//UuDd/YzZ4iDgrpavN41zXfrOBEGSJEmVaTIGISJeABwIHJaZy+tt5wK7AZ8ERk0QgHcCS4AnZeaN9bEXUCUZ7wPeWG9bVn99YmZ+oD52ZUQsAI6OiBNbbvz7EXPYzzJzdaML1Gd2MZIkSVIle/iamEOBNcAZI1XLTOAU4DGjde1pOfbs4Rv5+thVwJnAYS3lngfMq2O2OpnqIXprd6Z+xJy2TBAkSZLUL4siYvVYr1GO2wu4PDPbRzNc0rJ/GxExH9gduLTD7kuAnevuQsMxEristVBm/gZYP3yOfsRs86uI2FqPtfhyS6wpYxcjSZIkARDTpIsRsAy4ssP221v2d7IEiJZyox17S/1+d2Zu7FD2jpZz9CMmwNXAe4GLqcYd/CnVeItnR8QTM/OOTt/gZDBBkCRJUr+syczFXR47VrYyXiZTemyTc/Q0Zmae2rbvnIi4EPg+cBTwd2PE6SsTBEmSJFWmTwvCKjq3Eiyt3zs9zYfqKX0WHrsKWBARczs88V/SUq4fMTvKzLMj4ibg6WOV6zfHIEiSJKky1MPXxFwGPLbDegR71++dxgOQmeuBa+jc139v4NbMvKXlHAHs2VooIh4JzB8+Rz9ijmMGvbiCE2CCIEmSpOlmObAYeFHb9tcCV2Tm5eMce1BE7DK8ISKW1rFOayl3FrAReE3b8cPrHJzZ55jbiIjnAA8ALhyrXL/ZxUiSJElA9niQ8oRirQDOBU6q1xa4luomex/g4OFCEbES2C8zo+XY46lu0FdExHHcs6jZFuCjI7XLXBURHwPeHxFr6vM9HXgXcEJmXt/PmBFxMfB14ApgM9WCau8ErgI+18U16xkTBEmSJE0rmZkRcQjVzfdHqVoTLqdaOG3Mp/CZeXNE7Et1U38qVY+Z84FnZuZ1bcU/RLXewlHAe4AbgQ8CH+93TODXwJuABwGzgeuBrwAfnuqF00wQJEmSVJk+g5TJzDuBN9ev0crsP8r239DS0jDG8QmcUL/GK9vTmJn5ivFiTRUTBEmSJFWmUYKgqWOC0EHdn63dwsmuhyRJkjTZTBAkSZJUjSnu5eSaNkYMLBOEDjr1Z4uI1cCiSa+MJEmSNIlMECRJkgTQ42lONahcKE2SJEnSCFsQJEmSVLEFQZggSJIkaZgJgrCLkSRJkqQWtiBIkiSpYguCsAVBkiRJUgtbECRJklTp5UJpGlgmCJIkSQJcB0EVuxhJkiRJGmELgiRJkiq2IAhbECRJkiS1sAVBkiRJlSFbEGSCIEmSpGF2MRJ2MZIkSZLUwhYESZIkVWxBELYgSJIkSWphC4IkSZIqtiAIWxAkSZIktbAFQZIkSRWnORUmCJIkSQJIIId6G08DyS5GkiRJkkbYgiBJkiQgezxI2SaEQWULgiRJkqQRtiBIkiSp4iBlYYIgSZKkYa6DIOxiJEmSJKmFLQiSJEmq2IIgbEGQJEmS1MIWBEmSJFVsQRAmCJIkSRo21MOVlDWw7GIkSZIkaYQtCJIkSarYxUjYgiBJkiSphS0IkiRJqtiCIGxBkCRJktTCFgRJkiRVhmxBkAmCJEmSaplOcyq7GEmSJElqYQuCJEmSKnYxErYgSJIkSWphC4IkSZIg6e00pzZGDCwTBEmSJFWGHKQsuxhJkiRJamELgiRJkoDs8UrK9jEaVLYgSJIkSRphC4IkSZIASMcgCBMESZIkDetpFyMNKrsYSZIkSRphC4IkSZIqrqQsbEGQJEmS1MIWBEmSJFXSQcqyBUGSJElSC1sQJEmSBEA6BkHYgiBJkqRhOdS71wRFxMKI+GxE3BQR6yPiZxHx4sJjd4+I0yNiTUSsjYgVEbHHKGXfGhFXRsTGiLg6Io6OiG3ukac65mQyQZAkSdJ0tBx4FXAM8ELgcmB5RLxgrIMiYmfgfGBX4HDgFcBS4LyIeEhb2WOATwPfAp4LnAR8BPjodIo52exiJEmSJGD6dDGqk4ADgcMyc3m97VxgN+CTwIoxDn8nsAR4UmbeWB97AXAt8D7gjfW2ZfXXJ2bmB+pjV0bEAuDoiDgxM2+Y6phTwRYESZIkTTeHAmuAM4Y3ZGYCpwCPGacbzqHA2cM33fWxq4AzgcNayj0PmFfHbHUy1UP01u5MUxlz0tmCUG6nLWxmZZ4xfkn1Tz8ebGzpU9lSm/sQE5rVdVOUlcvp8WSp5zYWfv9Afz6EfdDk59+v779fn+2ptHGqK6DtzZbqF2Wnqa4HVHVZObS8p/GARRGxeqxymbm4w+a9gMsztxnMcEnr/vaDImI+sDvw3Q4xLwFeGRE7Z+YtdYwELmurz28iYn29fzrEnHQmCOWGgBlb2HznVFdkQCyq39dMaS22f725zgNyzzuF/DxPDq/z5PA6T47S67wT1T3GVFsDIzf1vdRt8rMMuLLD9ttb9neyBIiWcqMde0v9fndmdkr972g5x1THnHQmCIUy02vVwPDTglGeCqhHvM6Tw+s8ObzOk8PrPDkG7Tpn5sOnug4djPX4arxHW6XHNjnHVMacVI5BkCRJ0nSzis6tBEvr905P3qF6Sp+Fx64CFkTE3A5ll7SUm+qYk84EQZIkSdPNZcBjO6wdsHf9fmmngzJzPXANdV//Dsfe2tKv/zKqbj57thaKiEcC84fPMQ1iTjoTBEmSJE03y4HFwIvatr8WuCIztxmg3HbsQRGxy/CGiFhaxzqtpdxZVFMPvKbt+MOppno4c5rEnHT2q5ckSdJ0swI4FzipXlvgWqqb7H2Ag4cLRcRKYL/MbJ2O7XiqG/QVEXEc1Y35MfX7yGJlmbkqIj4GvD8i1tTnezrwLuCEzLx+OsScCpHb67SFmlKDNjhrUHmdJ4fXeXJ4nSeH13lyeJ0nLiJ2orpRfilVa8LlwIcy8/SWMivZNkEgIh5FdQN+AFWPmfOBd2bmZW3lAngbcBTwMOBG4EvAx9unWJ3KmJPNBEF94R/GyeF1nhxe58nhdZ4cXufJ4XXWIHMMgiRJkqQRtiBIkiRJGmELgiRJkqQRJgiSJEmSRpggSJIkSRphgiBJkiRphAmCeioi9oyIL0TERRGxISIyInYdo/wrI+KXddkbIuLvI2LeJFZ5IEXEwoj4bETcFBHrI+JnEfHiqa7XoIqIh0TEZyLiRxFxV/253X+UsgdFxIX1db8lIr4YEYsntcIDKiKeHREnR8QVEXF3/Tt/WkTs3aGs17lLEfGMiPheRPy+/tt6a0ScExHP71DW69xDEXFs/ffjFx32ea01MEwQ1GtPoloi/A/Af49VMCJeDXyjLvd8qsVQjgJO7m8VtwvLgVdRrbj4QqrFY5ZHxAumtFaD65HAK4C7gP8arVCdNKwArqf6nL8TeDHwHxHh39Px/TXVokGfpvqd/9v6659GxNOGC3mdJ2wJcAXwDuB5wF8CG6lWa/3z4UJe596KiD2pVsu9ucO+/fFaa4A4zal6KiJmDK8SGBFvp7oReERm/rat3EzgBuCizGxdMv0vqFYbfFpm/mSy6j1I6iTgP4DDMnN5vS2oVl9clpmPncr6DaK2z+0hVAnYAZm5sq3cRcBs4Ikt5Q8Cvg/8eWZ+ezLrPWgiYufMvKVt22LgWuCczHxJvc3r3GMRMYvqOv8mM59Vb/M690h9k/9j4KfA3sDizHxcy36vtQaKWat6qn0J8TE8DdgFOKVt+zeAzcBLelmv7cyhwBrgjOENWWX6pwCPiYg9pqpig6rkcxsRDwaeDJzaWj4zzwZ+j5/ZcbUnB/W21cBvgIeA17lfMnML1d+NzeB17oO/ofoMv699h9dag8gEQVNlr/r90taNmXk3cHXLfm1rL+DyDje1l7TsV+91/MzW/heve1ci4v5U1274unqdeyQiZkTErIh4UEQcBzyaqlUXvM49ExG7AR8C3pyZd3Yo4rXWwDFB0FRZVr/f3mHf7S37ta1ljH7dhver9/zM9ljdNe5LVP8WHV9v9jr3zneoWgx+D7wdeHlm/me9z+vcA/Vn+MvA9zLz9FGKea01cEwQNKqI2L+ejaHkdb8uTzPaIBgHx4xtrOvjtesvP7O98w/AIcBfZ+av2vZ5nSfuaOApVINhVwDfiYhXtJXxOk/MX1BNzvGWgrJeaw2MWVNdAU1rvwaOLCy7tmHsVfX7spb/H7aUajCdOltF5ydOS+v3Tk+pNHGtn9l2S/G6NxIRH6GaZedtmXlyyy6vc49k5jXANfWXZ0bEmcDnIuLbeJ0nrH4w9gngY8C6lilLZwEz66834LXWADJB0Kgy8w/0b8rRy+r3vYArhzdGxA7A7sCZfTrv9uAy4CWtM+/UhueS79TPVRPX+pn9ftu+valmMFGBiPgQ8F7g6Mz8bNtur3P/XAT8GXB/vM698BBgEVWC8LEO++8APg6cWH/ttdbAsIuRpsqFVGslvKZt+yuopoI7bdJrNDiWA4up5tJu9Vrgisy8fNJrdB+QmTcAPwNe1TpveUQ8G3gwfmaLRMQHgfcD78/Mf2jf73Xuj7qv/P7AamCV17knrgIO6PD6JdVkGwcAX/JaaxDZgqCeqlsAhhfr+pP6/fkRcStwa2aeB9WUexHxbuDkiDgR+FfgsVRPW/41My+c5KoPkhXAucBJEbGMqjvW4cA+wMFjHajRRcRL6/99cv2+X92FYF1mnlVvexfVE8B/iYgvAQ+i+sz+BPjuZNZ3EEXEO4BjgX8HftC6OBqwMTMvrv/f6zwBEfEN4HfA/wC3AQ+k+hvxLOAt9ZSn4HWekMy8C1jZvj0iVtf7W/d5rTVQXChNPRURuzL6+IHzMnP/tvKvpvrD+Wiqf8i+AXwwM9f3sZoDLyJ2olp5+qVUrQmXAx8aYxYNjSMiRvtj+LvM3LWl3POA46gS4LXA6VRdZe7odx0HXUSsBPYbZbfXuUci4s1UK60/mqoLzBqqJ9gnZuaZbWW9zj1Wf87vtVBavd1rrYFhgiBJkiRphGMQJEmSJI0wQZAkSZI0wgRBkiRJ0ggTBEmSJEkjTBAkSZIkjTBBkCRJkjTCBEGSJEnSCBMESZIkSSNMECRJ24iIvSJiS0Qc1OfzHBIRmyLiUf08jySpnAmCpDFFRDZ47TrV9Z2OIuJxEXHsgF2fTwH/nZlnt26MiAMi4qKI+GFE/PFoB0fEDhHx9og4LyL+UCcBt0bEWRHx6oiYAZCZpwP/C3y8n9+MJKncrKmugKRp7zVtX+8L/CXwJeD8tn23TkqNBs/jgA8CK4HfTmVFSkTE04GDgEM67P4E8DxgEfBl4Nkdjn8K8K/AQ4CzqJKNVcCuwMuAU4GdgM/Xh3wGOCUi9szMy3r4rUiSumCCIGlMmfnPrV9HxCyqBOGC9n33BRExE5ibmXdPdV2G9aFOb6K6oV/RYd8sYD0wE5jdoS6PB34A3A3sk5k/btt/LPC3wC9aNp8GfAH4a+AtE669JGlC7GIkqaciYm5EvDciLouIDRGxOiLOrG8c28seUXdNenZEfCAifhcR6yPiJxHxtLrMfhHxo4hYFxE3RcT7R4lxYN2N53cRsTEiLomIP59IHdtivz8irgY2AC+PiB0j4u/qut5Wn/OqiPj7iNihJcaxwNfqL89t6Y518vD+0bpnRcRvI2JlaZ2aXv9Rrs0sqpaDszNzc4cifw9cBPwn8N726wp8E5gL/Fl7cgCQmVsz8x9a92XmXVStUS8rqaMkqb9sQZDUMxExm+rG8RlU3UhOpOqK8hfAf0fEMzPzZx0O/XuqJ9KfAeYA7wC+FxGHAydRdWf6BtVN8Ici4toOrRcfBxZQPYlO4EjgXyJiXmaePME6Hk/1tPzLwJ3AFcCDgTcA/4/qpngLsB9wNPB44Ln1sacBD6Rqdfko8Kt6+9WjXcdC29RpAte/1ROBhVRJwDYy89vAt0c59kjgMcCnC87T7gLguRHxmMz8dcNjJUk9ZIIgqZfeDOwPPC8zvze8MSI+D1xKdVO7f4fjZgJPy8xNdfnLgTOo+rE/PTN/Wm8/CfgdcBTQniDcD/jjzFxTl/0n4BLgUxHx7cxcP4E6zgce39qFJyLmAA9te8r+uYj4MHBMRDwlMy/KzEsi4gKqBOHszFzZ4fvvRqc6/U0X31u7Per3bhKYv6JKzr7QxbHD59sTMEGQpClkFyNJvfRqqpu7/4mI+w2/qFoFzgb2iYj5HY77wnByUBse/HzhcHIAUJe5COg0JeYXhpODuuwa4J+AJdz7pribOn6hvX9/Zm4aTg4iYlZELKnj/KAu8tQOdeylbepE99e/1f3r99ubVKY+z58A/5uZv2lybG1V/b5zF8dKknrIFgRJvfRYqifbY81mdD/g+rZt17R+kZl3RATAtR2OvwNY1mH7rzpsu7x+322CdbyyU6GIeBPVwNo92faBy5Ix4vdCpzp1e/1bZf0eDevz8PqYKxoeN2z4fDlmKUlS35kgSOqloJrT/m/HKNPp5nXrKGVH295JpxvLTje53dRxm9mBIuJvgU8C3wc+C9wIbKIam3Ay5S20Y90Qj/U3utOMRd1e/077l45Trt3M+n2bmY0KDZ/PqXIlaYqZIEjqpd9QdVE5JzOHJvncewD/1rbtsfV7awtFr+r4Gqo1DZ7fGicinteh7FhJwHBXnqW0rJEQEfOoBjdf1aBOvfjeLq3fm65sfBUwBOwdEZGZTVsCHtl2fknSFHEMgqRe+jqwC6M8wY6IB/Tx3G+MiEUt51pE1f1nNXBeH+q4lerGf6SVop4i9N0dyt5Vv3d6Kj/cVejAtu1/Q/O/0b343i6mmhXpaU1OnJm3A/8B7A68dZTz7xYRfz1KiKcBN2dmt12UJEk9YguCpF76DNUKvP8QEc8CzqG62XwY1Yq7G4AD+nTu24CfRMRXqW7aj6zP+4a2wby9quO/Ah8DzoqI06hWBn4l0GntgJ9SPV1/X0QsAdYB12bmT6gGNf+aavrWZVTjLvahumG+rfzb7833lplb6+/n4IiYm5kbG5z/jVTjMU6IiOcD51J1GXog8EyqJOjv2g+KiIVUK3R/tcG5JEl9YoIgqWcyc3NEvJBqJd7XAMfVu26kmn3olD6e/l1UN5lvBh5A1d3mVZn5zT7V8R+oEpHXU92Y/4FqfYCvcc/g6OFzXhcRr6vr+AWqfvqnAD+pb8gPphrH8BaqcQzfp1pT4b9Lv/kef29fAI4A/oxqnYfS8/8+Ip5A1YJxCPCBetdNVIPIjwK+2+HQlwA7AF8sPZckqX+ieTdRSZo+IuIIqpvyA3q4xsB9XkT8J7AgM/edhHP9D/C7zDys3+eSJI3PMQiSpE7eATw9Ip7Tz5NExCHA3lStK5KkacAuRpKkbWTmZUzCvxGZeTrVQm6SpGnCFgRJkiRJIxyDIEmSJGmELQiSJEmSRpggSJIkSRphgiBJkiRphAmCJEmSpBEmCJIkSZJGmCBIkiRJGmGCIEmSJGnE/wciDRk9E16gXgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 864x576 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.pcolormesh(temp_bins, flow_bins, pdf, shading='auto')\n",
"plt.colorbar(label='Probability density')\n",
"plt.semilogy()\n",
"plt.xlabel('Temperature (${}^{\\circ}C$)')\n",
"plt.ylabel('Streamflow ($m^3/s$)')\n",
"plt.title(f'Joint PDF at {site_to_test}')"
]
},
{
"cell_type": "markdown",
"id": "5226204f-2c44-4b02-b282-9cd20f5427e0",
"metadata": {},
"source": [
"## Plotting the conditional CDFs"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "f441f0ef-0884-4542-bafe-682d12f0629e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 1.0, 'Conditional CDFs at KIOW')"
]
},
"execution_count": 5,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAugAAAIPCAYAAAAhGyuqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAACnrElEQVR4nOzdd3hUVf7H8feZkkIKCQm9hd4RBAUVlWZXFMuqu6Luuqvrb3dde+9d1rLqVte+trWBDRtSFBWwgID03lt6QsrM3PP74wYEDJCESW4y+byeZ54kt81nEki+Ofnec4y1FhERERERqR98XgcQEREREZGfqEAXEREREalHVKCLiIiIiNQjKtBFREREROoRFegiIiIiIvWICnQRERERkXpEBbqIyAEYY6wx5vkDbTvANVYbY6ZFOdpBM8YMr3gtF3udRUREXCrQRcRzxpgmxpgrjTFfGGNyjDEhY8wWY8wkY8zFxpiA1xmrwhhzpzHmDK9z1DZjTLox5nZjzDfGmDxjTLkxZr0x5i1jzJnGGLPbsasrfgHY+Sgyxqyt+NpeYYxJ28dz2P08bqyzF/tTnjOMMXdW85xpxpiiSrbHGWPeqHgtf9v5+ar4+P19XOuYinM2Vny+t1Z8Ds+o5NjXKq7VpZJ9/67Y999K9rWq2PdedV6niERfg/ihJyKxyxjTFfgA6A5MBh4AtgMtgNHAc0Bv4HqvMu5DIhDZa9sdwAvAxEqO7wE0+JXhjDGHA+/gfn3eBV4GCoC2wMnAW8AfgH/sdtp64KaK9xOANsBw4HHgFmPM+dbaKZU83VzgkUq2zznY11EDZwAXAXcezEWMMU2At4ETgHuttbdV4Zz7gJuBNcAzwCqgFfBLYEJFsf1ra+3Of49TgXOBEcCKvS43HAhX7Nvbzm1Tq/GSRKQWqEAXEc8YYxKB94HOwFnW2rf3OuQhY8xhwGF1Hu4ArLWl1Ty+rLay1BVjTCvgPdwi+1hr7Yy9DrnHGHMCkL7X9nxr7Ut7bbvbGHMsbpH/jjFmoLV2+V7HbKjkvAar4q8FHwBHAFdbax+rwjmX4Bbnk4HTrbU7dts3HrdgvxBYDdxesWtngT0ceHq341vh/iL8HPBrY0zXvT7nwyveTqvWCxORqFOLi4h46be4I8uPVFKcA2Ct/cZau/to7M52gy8r2iWKKt4/fe9zd/Z9G2N6GmM+MMYUGmPyjTFvVhQrex/fxxjzkTGmuKLV5iVjTIvKcu3eg26MyTLG7Bwdv2j3doy9s1Rynai/FmNMG2PMI8aYucaYXGNMqTFmoTHmBmOMv7LXU0XX4Y6c31BJcQ6AtfZja+1rVbmYtXY6cA2QDNSobcUYc6Qx5kNjzOaK17mhovVjaBXOPdwY87wxZqkxZkfF5/RLY8zYvY6bhjt6vnfrzcXVyNkSt/AdAvymisV5HHAvUAT8cvfiHMBaGwYuA9YC1xpjmldsXwps4Oej5Ds/vp/KR9GHA3m4f7kQEQ+pQBcRL51d8fapqp5gjPk/YALQDLd4uafi/YnGmEsrOaUtbmG0FrfAfAU4E3hxr+t2Ar4Ajgb+hjsa2Rz4qAqxtgHjKt7/ouL9cbttq9PXAvSv2D4FuBW3+F0HPMierSfVdRZQjtvGEy3/Bcpw22P2FjTGZO71aLZzpzGmB/Ap7qjw48D/AX/HbSU6pArPPRboCbwO/Bm4D/fz/7Yx5pe7HXcf7tcV9vzafl6VF2iM6Vhxfk/gHGvt81U5DzgKt5XlHWvttsoOqPhLzku4LVe7fw6nAW2MMd122zYcWFIxav49P42YY4xpjft5nG6tdaqYT0Rqi7VWDz300MOTB5ANFFTj+HTc0cTlQOpu21Nxe20LgbTdtq/GLdZ+sdd1dhZxPXfb9krFthG7bTO4BbQFnt/rGlXatleWaXX0WhIBU0mG/+L2zbfebdvwivMvPsDnPqXiuHnV/BqvBhYc4Jh5FddO2etzWdlj827HXFGx7fAa/vtLqmRbE2AJsHCv7c+7PzKrdf1puL/QrKv4eo46wPEWeH+3j/9Use3qA5x3VsVxD++27ZKKbZfutm0J8K+K9x/CbSHaue/8iuOvrMnnUg899IjuQyPoIuKlVNwbDKvqOCAJeMJau+u8ivefxG2VGL3XORutta/vtW3nDYldAYwxPuA04Ftr7a4b5Ky1FhhfjXzVUSuvpeIaJRXZd84Y0swYkwl8jPuX08E1yJta8bY6X6+q2nnN1L22z8L9PO3+OGe3/fkVb083xiRU90mttcU73zfuTEIZuAX6FKCXMWbvPDXhx20LKsQt1Ktj5/Pn7/eon/Y33W3b7n3oe4yQV2yfjjvC3n3349ANoiL1gm4SFREvFeCOzFZVp4q3P1ayb0HF2857bV9ZybHZFW8zKt62wC2IF1dy7MJq5KuO2notGHdayhtxbx7sivuXgN3tfRNnVewsoqvz9aqqfRX/2621k/dz3mvABbg3UV5ljJmJ+0vIa9baNQd60or7C+4FTsf9N7C3tEoyVVcJ7uj0m8A0Y8wIa+2SKp6787mb7veoSgp5a+1KY8xafiq8d77d2ZYzA/evKSOApRX7c3D/miEiHtMIuoh4aQGQaozZuxDdl70LzarYeyrEyq63821dToNYW68F4FHcfvbvgV/j9iYfB9xQsb/a3/uttYW40/z1rJh9JyqMMfG4I7ubKp6jOpnKrLXH4d54+QDu5+duYPHeN3pW8rwG+AT35s8XcaclPBH38/RKxWFR+RlprX0Ptw0lA7dI71nFU3f+onboAY7buX/+XtunAq0revWHA8uttRsqMhXg3gw6fK/+8wY/FahILFCBLiJeeqvi7W+rePzOOZ37VLKvd8XbykaZD2Qrbj94r/1cN9pq67VAxQ2M1trzrLUvWGs/rBiJPtjR4LeBOA5w82s1jQPicacfrBFr7Wxr7T0VxXpXoBh3ZHx/+uPeSPqgtfY6a+3r1p2BZjJuW8rPnqam+Soyvo9bpDcDphpjKvu3trevgC24LTyZlR1Q0dpzAVAKfLjX7p3tKiNwC/Tpe+2fXrF9+F7Hi4jHVKCLiJeexr1x7drKphYEMMYMqpjtBNwZO4qBPxljUnY7JgX3hrqiimOqxboLvLwPDDbG7Jp6rmKUtToLJBXhFmBVUSuvpUKEvUbojTFJwFU1vN5O43FnrBlvjDmisgOMMccbY86rysWMOw/6I7j92Q9UN8w+itb1FRkP9HXY+deIvT9PfXFnd9lbUcX+qn59f6aiSB+L22I01Riz31/+rDt3/u247Vcv7f2Xi4opM/8BdAT+Yq3dutcldhbc57Nn//lO03Fnifl9xcfTqvN6RKT2qAddRDxjrd1hjDkVd/R0ojHmE9yiNBt3isMRuCsujq84Ps8Ycz3uzCWzds5DDlyMO3J6mbX2QDfU7cutwEnA+8aYJ3ELvdMqclTVTGC0MeYG3KkQrd3HnOC1/FreBC4zxvwPd4GblsBv+KlfvUastZsrvl7vADOMMRNxe5oLcFcHPREYBly+16lNjTEXVLwfX3HszlHdrcB51tqa/LXgVmPM8bi/XK3CLbZPw53O8EA39y7C7f+/3rirey7BLWIvw20t2butZCbwR+AfxpgPgBAwy1q7qjqBrbWTKtpvJuAW6SOttZXdh7Dz+KeMMV1wf1FcaIx5EXdmnFa4hXc/3GkW76rk3LXGmJXAMRWb9i7QvwCciv3b+amlRkQ8pgJdRDxlrV1ujBmIWxidBdyCO2KYA3yL2yP8ym7H/8MYswl3HvA7Kjb/AIy11k48iBwrjDFH447o/gl3bu4PcVswtlTxMjvn4b6Fn26m3OeiPbX1WoCrcUelf4F7A+Q63Lnmv8Et2GvMWjvbGNMH93N0Gm7PdyJuoT0Td7XLd/c6rR3uFI/g3jSZjVsMXgm8aK3Nq2GciUBr3NfZsuLay4Df4a6wub/XETHGnAI8jPtvLKki00W4rS97F+ivAgOB83BnkvHh9vdXq0CveO4PKynS91kcW2tvMMZ8iPs5vxS3lz0f9//HHdbaCft5uqm4Nxuvttau3eu6ucaY+bivV/3nIvWI0f9HEREREZH6Qz3oIiIiIiL1iAp0EREREZF6RAW6iIiIiEg9optEd2OMCeP+0lIbS1mLiIiIiOyUCjjW2p/V47pJdDfGGAcwTZseaFVlEREREZGay8/PB3c63p91tGgEfU8FTZs2bZqXl+d1DhERERGJYWlpaeTn51fataEedBERERGRekQFuoiIiIhIPaICXURERESkHlGBLiIiIiJSj6hAFxERERGpR1Sgi4iIiIjUI5pmsZqstWzfvp3S0lIcx/E6jtQin89HQkICmZmZGGO8jiMiIiKNhAr0arDWsmHDBgoLC4mPj8fv93sdSWpRKBSiqKiIsrIy2rZtqyJdRERE6oQK9GrYvn07hYWFtGzZkmbNmnkdR+pATk4OW7ZsYfv27TRv3tzrOCIiItIIqAe9GkpLS4mPj1dx3og0a9aM+Ph4SktLvY4iIiIijYQK9GpwHEdtLY2Q3+/X/QYiIiJSZ1Sgi4iIiIjUIyrQRURERETqEU8LdGNMO2PM48aYGcaYImOMNcYMr8b5XYwxE40x+caYQmPMJGNM79pLLPXRa6+9xgknnMDf//53r6OIiIiIHDSvR9C7AucDRcBn1TnRGNMC+ALIAi6quE4zYLoxpl10Y8Y+Y0yVHqtXr67R9devX8+f//xnhg0bRnJyMsYYpk2bts/ji4qKuOKKK2jdujWJiYkMHjyYd99992fHbd68mRdffJE33niD2bNnM2fOnEqv9+WXXzJu3DiysrJISEggNTWVYcOG8be//Y3y8vIavSYRERGR2uD1NIufW2tbABhjzgDGVOPca4F0YLC1dmPFNb4GVgG3AJdHN2ps+/rrr/f4+IYbbmDp0qVMmDBhj+2tW7eu0fWXL1/Oq6++yqGHHsqoUaMqLbZ3N3bsWL7//nvGjx9Pp06deP755xk7dizvvfceJ5988q7jCgsLad26NampqXTu3Jns7Ow9rmOt5dprr+Wxxx5jxIgR3HnnnWRlZZGTk8M777zDlVdeSXJyMhdffHGNXpeIiIhItHlaoFtrD2ZqjLHApzuL84rrZRtj3gPORAV6tQwdOnSPj9PT04mPj//Z9po65phj2Lp1KwATJ07cb4E+adIkJk+ezNtvv83YsWMBGDFiBCtXruSaa67Zo0Dv1q0bTZo0YciQIXTv3p3bbrttj2vddNNNPProo4wfP57rrrtuj31nnnkml19+ORkZGVF5jSIiIiLR4PUIeo0YYxKBLsAbleyeB/zSGNPCWrt1r/PyDnDpptFJKHvz+areTTVhwgSaNm3K6aefvmubMYaLLrqISy+9lIULF9K790+3Gjz55JOVXmfmzJmMHz+ecePG/aw43ylav4CIiEjsc8rDhAuKiRSWEC4qIVK4o+JtCZHS6LRLRsodQsVhyovDFW8jhIpDECrDlJdiQqUVb8vARvY82YIFcCzWgq14627cmyUQHyKYWEZck3IC8aF9HRg1xljimpQRn1xKfFIpwcTyWn/OqsrJ/CPtLzjD6xi7NMgCHbe1xQA5lezbuS0D2FrJfjlI1loikciBDwQCger/E1uwYAG9e/f+WVHfv3//PfYfyH333YcxhrvvvrvaGUREpHGLlJZTsnITO5ZtYMeyDRQv30Dp2m3gwboYwSodZUlMLyKYWLbHVgP448IEE8oJJpYTSCyveL8Mn79+FMf1ga1iXVNXGmqBvtP+/mX9bJ+1Nm1/F6sYYa/2KHok7FCwvqS6p9Wa1HaJ+AO1d//v9OnTGTFiRJWOXbVqFVlZWdW6fnZ2Nt27d//Z9p0ruO7dZ16ZwsJCPv74Y0aPHl3t5xcRkdhkHYfSLfkUrdxOKKeIyI5SnIpHpKSs4v0SQltyKN+cjTv8vJ/rBYJYfxyOCRJx/ETCDk7I7nrrhJ2DGiA2PoM/zocv4CNCkJANErZxhB33rYOPQKCcdlnL6dR1IalpuTV6HsfGESaZ2p47xAJhm0rINqt4pGE9n6/ElTF0kNcR9tBQC/Rc3K9zZc3DzSreVja6XisK1pfw167v19XTHdCVy08lPSup1q4/aNAgvvnmmyod26ZNmxo9hzGmRvt2mj9/PqFQiAEDBtTo+UVEpO45EYeywjBlBaGKR5iy/BClBSFK88t3+9h93wlXPprtj5SQUL4ds6MQX1kxAWcHQVNKQqAcn696FXNxSRz5BYnkFzYhrzCRHSVxhEJ+QuEA1h7459Hu4pIDNMmMJykzniaZ8TRpXvF+83iSdns/uWUCSS0SiEv6eZlmrQOlBdiCTUSWfIKz5BMI7fjpgISmYPYsek1cE0hugUnKxCS3wCQ3dz9OboFJaQlxSVX62Sp1p0EW6NbaEmPMSqBvJbv7Adv27j+X6ElOTq5y4VuTFpeMjIxKR8lzctzfuXaOpO9PXl4eUPNZZ0REpHrCZRGKNpdStLmU0rxyyorClBeG3IK7sOL9ndsKwpQVhigrDFG+qyB3+65rwhiHZmnFtMwopEVmAanJpT/t3Ed/SDjsIxT2u4+Qf4/3S8uCuwryUKjyn2MJTYMkZsSR2CyOJs3iScyIc4vuiveTMuPdfRlu0d0kM55ggn+/r8OWFWELNmJ35MCabMI7crA7sqHYfWuLs2FHNjh7fZ4C8fi6jcLf93R8zX/+F2hpeBpkgV5hAvBHY0wra+1mAGNMM+A04NW6DJLaLpErl59al0+5X6ntEmv1+rXd4tKnTx/eeustHMfZow99/vz5APTtW9nvZXtq2bIlAGvXrq3Wc4uIxDprLeVFYYq3lbFjWxklOeVY58CjytZCWWGIos2lFG4q2VWMF2523y/JqZ01JQLxPuJTg8Q3DZLQNEh8ivt+k6QwSeGtJJRsJr50Gz67Z9HqmCCRxFRISsHfNIVAs6bEtUgnsW06TTpkEJfWpEZ54lICJKbH4Q/WvDXDWguFm3G2L8duX47NXoGzfTkUbqnWdUzTdvj6jsHf40RMQkqN80j943mBbow5u+LdwyreHmuMyQSKrbUfVhwzDTjW7vm3pIeBccAkY8xdQBi4teLt/XWRfSd/wFerLSX1TW23uIwdO5ZnnnmG9957b4+ZXF588UV69OhRpRtEBwwYQIcOHXjllVe4/fbbSUtL+9kxU6ZMYeTIkdXOJyJSH1lrKdpcyvbFBeStKaZoq1uAF28vo3hrKTu2l1G8tYzibWWES2v3hjif3xCfGiQuJUB8SoC45KD7NmXn2wDxyW6hHZ8aJCE14BbeqT894lICJDQNEoj/adS5bGM2uV8uIHfGD+xYsv5nz9ukaxtSB/eg6WHdSerZHuPf/4h1XbOREM6yz4j88AY2e+W+D/QFIakZpkkGpkkGNEl321OSMjBJmdDE3UdimlpTYpTnBTo/nyrxzoq3a3BXCa2UtXaLMeZo3EL9v7h3NnwBHGOt1bBpLUpJSWHw4MHVPu/NN98E2FXcT58+ne3bt5OUlMRJJ52067iTTz6ZESNGcMkll5CdnU2nTp144YUXmDFjBu+8806Vnsvv9/PUU08xZswYBg0axDXXXEPPnj0JhULMmTOHl156iT59+qhAF5EGJ1wWIWd5EduXFLBtcSHZSwvYvriQ7UsLKSsIVft67k2IVSvy4pICJLdOIKVVIsmtEkiueJvSOoHklgmktHY/jksORK1wLFm7lbwZC8j9cgElKzbtmT05gdRDu7lF+eDuBJvVz1FkW1pI5Md3icyf4Lao7OSPw2R0xmR2wZfZFZPRFZPeHuJTVXg3csYe4A7lxsQYk9e0adOmO/uX97ZmzRoAOnbsWIepvHHGGWcwd+5cVq9eHbVr7uubTceOHX/2PAUFBdx88828+eab5OXl0bt3b26//XbOOOOMaj3n3LlzefDBB3f9MpCRkUHHjh0ZPXo048aNo2fPnge8RmP6uouId5yIQ0luiJKcckpy3NaTktxydmSXk7+2mOylhWxfUkjuquL9tqQYAyltE90bDZsn7LoRMalFxY2IO7c1j9/njYheK9ucS/an35L7xQJK1+55S1kgPYX0o/qQNqwvKf2y6t0o+e5s/gYi894isuhDCFf0xRsfvi7D8fc/C9OiB8ZXf/NL7UpLSyM/Pz+/slkGVaDvRgW6VEZfdxE5GNZaCjeVsnlOLpt+yKNocwklOW7hXZJbUYjnlFOaV73R77ikABk9UsjskUJmj9RdbzO6JRNMrH9F94HYiEP+t0vY9v4sCr5duscUh8HMpqQP60P6sH4k9eqA8dePqfkqY63Frv+eyLy3cNbMZNc8i8Em+Hufgr//mZiUVp5mlPphfwV6w/sfLCIiUk9Za8lbXcymuXls+j6XTXPcR9GW0gOfvJdAvI/EDHcmkOSWCWT2rCjEu6eQ2TOV1LaJMdEGEcotZPvH37J90mzKt+bt2h7MSKXZ8ENIP7ovTbq3w1RjRWov2FAJzpJPicx/G5u75qcdyS3x9z8Tf6+TMfHJ3gWUBkUFuoiISA04EYfsZUVsnpvLxt2K8X2NhPuDPpr3SSW9UzKJzeJITK+Yoi8jjsRm8bs+3vmIaxK7P6KttRQtWM2292eS9+WP2PBPN62mDOxK81OGkDa0FyZQv9s/rBPBbl5AZMV0nKWfQlnRrn2mTX/8/c7C1+kotbFItcXu/34REZEociIO62dms/TDTaz5fBubf8jb57zdwUQ/Lfun0frQdFoPTKf1gDRa9Gm6x4wkjVHJ2q3kfj6P3OnzKF23bdd2f3IiGccPovnJh5PQrrmHCQ/MOhHspvlEVkzDWfkF7NhtXUR/0J2PvP9Z+DK7ehdSGjwV6CIiIvuwI6eM5R9vZumkjSz/eHOlc33HpwZpPSDNLcQPTaf1gHQyeqTgD9Tvloy6UrphO7nT55H7+XxKVm/eY19Sz/Y0P2UI6cf0xxe/jxWF6gknb73bV75iOpTk7rHPNO+Gr+sI/D1PwiSmeRNQYooKdBERkQrWWrbMz2fppI0sm7SJdTOzfzZjSos+qXQ9vjXthmTQemA6aZ2S8Pkafi94NJVtyiH383nkfD6fkhUb99gX1yKN9GP602z4ITTpWv21Muqas305ke9exln5OVhn13bTvAe+rsfi73wspmn9fx3SsKhAFxGRRq18R5hVn21h6YebWPbhJvLX7dhjfyDBT6cRLeh+cmu6ndSmUS1MVx2R0nJyp/7Ato9m/2wRoWBmKulH9yP9mP7uAkIN4OZWZ9N8tzBfO+unjSmt8Pcdg7/LcExqa+/CScxTgS4iIo1O7qoilk7axNJJG1k9bSvhMmeP/Wkdm9Dt5DZ0P7E1WSNaxPQNmwerZPUWtk2aRfbk73F2lO3aHkhPIf3ovjQ7pj9JvTvU/1lYyndg89djc9YQWfQ+duO8XftMsyz8A3+Jr9tI3fApdULfcUREpNHYsiCPD6+cw6ppey5+4/Mb2h+ZSfeTW9P95DY0762VHPfHKQ+T9+UCtk2aTdH8Vbu2m2CA9GP6kXn8IJL7dqq385U7WxbjbJ6PzV3nFuV566B4+8+OMy164h/0K3xZR2JM/XwtEptUoIuISMwrLQgx7e4FzHpyGU7E7SlvkhlPtxNb0f3kNnQ5rhWJ6XEep6z/yjbnsH3SbLZ//C3h/OJd2+NbNyPzlCFkHjeIQNP62wLkbFpA+Jvnseu/2/dB8cn4WvbBP+AcTNtD9YuaeEIFuoiIxCxrLQv+t5aPr5tL4SZ3saCMbsmc+PBAup7YCl89HeGtb8o257Du3x+QP3PRTyt8+nykHdGL5qcMIWVAl3rdwuJsXkj4m+ew6779aWNqG3wZnTBp7fd4kNBURbl4TgW6iIjEpK0/5vPBFd+zerrbzhJM9HPMzb058uoejX4+8qqykQhbJnzJxv9Oxpa5CzAFM1PJPOlwMk8YTFxmU48T7p+zZTGRb57f40ZP06In/sMuxtfhcBXiUm+pQJcG77XXXuO5555jzJgx/OEPf/A6joh4rKwwxLR7fmTmE0txwu5ob68z2nLiIwNJ61h/2y/qmx3LNrDm8bfZsdydJjGYkUq7S08hfVgfjL9+/4Lj5Kwh8vW/cdZ8vWubad4N/2G/xtdxqApzqffq79+jpE4ZY6r0WL16dY2uv379ev785z8zbNgwkpOTMcYwbdq0So+dNm3aPp9/8eLFexy7efNmXnzxRd544w1mz57NnDlzKr3ml19+ybhx48jKyiIhIYHU1FSGDRvG3/72N8rLf77wiIg0PNZaFry+lif7fMhXjy7BCVuadU3mgveO4bw3h6k4ryKntJz1/5nEoj//Y1dx3vzUIfR56iqaHdu/XhfnNlRKeObThF7/7a7i3GR2JXDSPQTP/jf+rCNUnEuDoBF0AeDrr7/e4+MbbriBpUuXMmHChD22t25ds3lfly9fzquvvsqhhx7KqFGjePfddw94zkMPPcQxxxyzx7asrKw9Pi4sLKR169akpqbSuXNnsrOz99hvreXaa6/lscceY8SIEdx5551kZWWRk5PDO++8w5VXXklycjIXX3xxjV6XiNQP2xYV8MEV37FqqtvOEkjwc8xNvTjymp4EE+pvQVnfFHy/jDVPTKB8s7tSZkKHFnT881iS+2R5G6wKnLWzCX3+VyjY5G5IbUPgiMvwdR6mGVikwVGBLgAMHTp0j4/T09OJj4//2faaOuaYY9i61f3BOXHixCoV6N27dz/g83fr1o0mTZowZMgQunfvzm233bbH/ptuuolHH32U8ePHc9111+2x78wzz+Tyyy8nIyOjmq9GROqLsqIQ0+9ZyNePL9nVztJzTFtOfHSgFhSqhnB+Mev+8wE5k92/Qpqgn1bnDqfVL4bji6vfpYIt2kb4y7/jrJjubvAF8A88H/+gX2EC8d6GE6mh+v2/TmKGrxbv7n/yyScr3T5z5kzGjx/PuHHjflac7xStX0BEpG5Za/nxzXV8fO1cCjaUANCsSzInPTaQ7idr2fWqstaSM3Uu6//9wa5pE5P7ZNHhz2NJ7NDC43T7Z50IkQUTicx6FkLu6q+mzQACx16FL72Dx+lEDo4K9CiwkQjl2wq8jrFLXPPUWu0RtNYSiUSqdGwgUPN/Ypdddhlnn302SUlJHH300dx1110MGjSoyuffd999GGO4++67a5xBROqfbYsLmPTn71n52RbAbWc5+sZeHHWt2lmqo2DuCja+8AnFi9YC4GsST7tLTiLzpMPq9ZSJALZ4O6FJt2C3LXU3JKYROPJyfN2PU4+5xAQV6FFQvq2ABReP9zrGLn2fv574Vum1dv3p06czYsSIKh27atWqn/WNH0jTpk258sorGT58OM2aNWPRokU8+OCDHHXUUUyfPp0hQ4Yc8BqFhYV8/PHHjB49utrPLyL1U3lxmOn3/sjXf11KJOQA0OO0Npz06EDSOyV7nK7hKFqwmg0vfkrRvJW7tqUd2Yf2fxhDXEaqh8mqxu7IJfTONdi8il8sep9GYOjvMAkpHicTiR4V6FJtgwYN4ptvvqnSsW3aVP9PzQMHDmTgwIG7Pj766KMZM2YMffv25ZZbbmHy5MkHvMb8+fMJhUIMGDCg2s8vIvXP4vc2MOmK78lf57YypHdK4qTHDqXHqWpnqariJevY+OKnFHy3bNe25L5ZtLnwOFL6d/YwWdXZ0gJC713nFue+IIGT7sHf8cCDNiINjQr0KIhrnkrf56/3OsYucc1rdwQkOTm5yoXvwbS47K5Vq1Ycf/zxVbq5FCAvLw+o+awzIlI/5K/bwaQrv2fxOxsACMT7GHZ9L4Zd35Ngon6EVcWOlZvY+OKn7iqgFZr0aEfbi44nZWDXBtMSYsuLCb1/AzZ7Bfj8BE68U8W5xCx9d4sC4/fXaktJfVPbLS774jhOlX+QtGzZEoC1a9dG5blFpG5Fwg6z/raMqXcuoLwoDECX0S055W+DyOiqVoaqKFmzhY3/nUzejAW7tiV2aU2bC4+j6eE9G0xhDmBDJYQ+uBm7dTEYH4HRt+DPOtLrWCK1RgW6VFttt7hUZvPmzXz66adVnnVlwIABdOjQgVdeeYXbb7+dtLS0nx0zZcoURo4cGZV8IhI962dn8/4fvmXTnDwAklrEc9KjA+l7bocGVVR6pWxzLhtf+IScaT+AdaeeTOjQgjYXHkfakb3r/Q2ge7PhckIf3obdNA+AwIjr8Xet2iCRSEOlAl2qLSUlhcGDB1f7vDfffBNgV3E/ffp0tm/fTlJSEieddNKu4371q1/RuXNnDj30UNLT01m8eDEPPfQQJSUlPPDAA1V6Lr/fz1NPPcWYMWMYNGgQ11xzDT179iQUCjFnzhxeeukl+vTpowJdpB4pzS/ns9vm880/l++sKxl8aRdG39efxPQ4b8M1EKXrtrHk2n/vmjIxvm0GbS4YTfox/TH+hlWYA9hIiPDHd2DXfwdA4Nir8Pc8weNUIrVPBbrUmXPOOWePj++8804AOnbsyOrVq3dt79evH6+99hpPPvkkxcXFZGRkMHz4cG699Vb69u1b5ec74YQTmDVrFg8++CD33HMP27dvJyMjg44dO3L66aczbty4aLwsETlIO+c0/+jqORRuKgWgZb+mnPaPwbQ/ItPjdA1H2ZZclt70DOH8YvwpibT73clkjBpYq9Pu1ibrRAhPvg9nzUwA/Ef9AX+fMR6nEqkbxu4cphCMMXlNmzZtuvMGw72tWbMGcAtKaTz0dRepPTkri/jgT9+x/OPNAASb+Bl+e1+O+HN3/MGGN+LrlVB2AUuu/Tdlm3LwNYmn+4O/Jal7O69j1Zi1DuHPHsRZ+ikA/iGXEBh0gcepRKIrLS2N/Pz8fGtt2t77NIIuIiJ1Llwe4atHlzD93oWES92Fz7qf0oZTnjiUtI5JHqdrWMKFO1h6y7OUbcrBxAXoetdFDbw4t4SnP/ZTcT7oAhXn0uioQBcRkTq1ZsY23v/Dt2z90V2BObVtIif99VB6ndFWN4FWU2RHGctueY7S1VswAT9dbr+AlH6dvI5VY9ZaIl/+A2fh+wD4+5+N//DfeJxKpO6pQBcRkTqxI6eMT2/8ge+fXQWA8RmG/LErI+/qR3xK0ON0DY9TFmL5HS+wY+l68Bk63XAuTQf38DrWQYnMfpbIPHdCAV/vU/Ef9X/6pU0aJRXoIiJSq6y1/PDSaj6+7gd2bC8DoM2gdE77x2DaDGrmcbqGyQmFWXHvyxTNd3/Z6XjlmaQf3c/jVAcn/N3LRL57CQBf9+MIHHuVinNptFSgi4hIrdm+pID3//gdq6ZuBSA+JcDIe/px+OVd8TXAaf/qAxtxWP2X1yn4ZgkA7X9/KpnHV3/q2/ok/MObRGY9DYCv89EERt6AMfr3IY2XCnQREYm6SNjh8/sX8sWDi4iUOwD0PqsdJz06kNS2TTxO13BZa1nzxARyP58PQJsLj6PFGUd5nOrgRBZ+QOTLvwPg6zCEwHG3YXwNc2pIkWhRgS4iIlFVWhDijfO+Yvkn7tSJaVlJnPLEoXQ/OTorCzdW1lrWP/UB2R9/C0DLs4+h1fkNe0XNyNJPCU97BADTdiCBE+/C+HU/gogKdBERiZq8NcW8fPoXbF2QD8DQK7oz6t5+xDXRj5uDtemlz9g64UsAMk8+nLaXnNige7QjK78g/NmDgMW07EPw5PswgXivY4nUC/qOKSIiUbHhm2xeOWMGRVtK8QUMY/41mIEXd/Y6VkzY8tYXbHr5MwCajTiEDn84vWEX52tmEf7kbrAOpnk3gqc+iAkmeh1LpN5QgS4iIgdt4YT1vH3hTEIlERLSgpz3xlF0GtHS61gxYfuH37D+P5MAaDq0F1nXnINpwDfYRlZ/Rfjju8AJY9I7Ejz1L5j4ZK9jidQrDfd/uEiF1157jRNOOIG///3vXkcRaXSstcx4eDH/O+dLQiUR0jsn8dsZo1WcR0n+t0tZ8+QEAFIGdKHzzedjAg33BsrIkk8Jf3gbRMoxTdsRHPMIJrGp17FE6h0V6AKAMaZKj9WrV9fo+p999hkXX3wxPXr0oEmTJrRr144zzzyT+fPnV3p8UVERV1xxBa1btyYxMZHBgwfz7rvv/uy4zZs38+KLL/LGG28we/Zs5syZU+n1vvzyS8aNG0dWVhYJCQmkpqYybNgw/va3v1FeXl6j1yTS2EVCDu9d/i2f3vgDAO2PyOB3X46mec9Uj5PFhtJ121j1wKvgWJp0a0uXO8bhi2u4N1CG571F+LP73baWzK4Exz6OScrwOpZIvaQWFwHg66+/3uPjG264gaVLlzJhwoQ9trdu3bpG1//Xv/5FdnY2V111Fb169WLLli2MHz+eww47jGnTpjF06NA9jh87dizff/8948ePp1OnTjz//POMHTuW9957j5NPPnnXcYWFhbRu3ZrU1FQ6d+5Mdnb2Htex1nLttdfy2GOPMWLECO68806ysrLIycnhnXfe4corryQ5OZmLL764Rq9LpLEqySvn9XO/YuVnWwDod14HTn/6cIIJDXd0tz4JF5aw/M4XiBSXEmyWQpc7xuFPbJg3UFpriXzzApFvXwDAtO7v3hCqthaRfTLWWq8z1BvGmLymTZs2zcvLq3T/mjVrAOjYsWMdpvLGGWecwdy5c2s8Yr63rVu30qJFiz225eXl0alTJ0aOHMlbb721a/ukSZM45ZRTePvttxk7dizgfoM/+uijyc7OZtGiRXtc509/+hOzZ8+me/fuvPDCC/h8P/1h6MYbb+Shhx5i/PjxXHfddT/LNXPmTDIyMujWrds+szemr7tIVeSuKuLl079g28ICAI69tQ8j7ujToG9arE9sJMKyW5+ncM5yTDBAj4cvJalHe69j1Yi1DpEZfyMy3x3s8XUcSuD4OzDBBI+TiXgvLS2N/Pz8fGtt2t77NIIudWLv4hzcf5jdunVj/fr1e2yfMGECTZs25fTTT9+1zRjDRRddxKWXXsrChQvp3bv3rn1PPvlkpc85c+ZMxo8fz7hx4yotzoGfjdyLyP6tm7mdV8+cQfHWMvxBH6f/5zAOuSDL61gxZd2/P6BwznIAsq4+q+EW55Ew4anjcZZ+CoCv2ygCI2/E+FV6iByI/pdEgXUiULTN6xg/SW5eq6uwWWuJRCJVOjYQ2Pc/sW3btrFgwQLOP//8PbYvWLCA3r177zESDtC/f/899h/IfffdhzGGu+++u0pZRWT/fnxzHW9fPItwaYTEZnGc9+ZRZB3z81++pea2TZrNtnfdlsNW5w6n2YgB3gaqIRsuI/zxXThr3Nfi6zeWwLA/YoxufROpChXo0VC0jfKXzj/wcXUk7oJXIbVVrV1/+vTpjBhRtdXrVq1aRVZW1s+2W2u59NJLcRyHa6+9do992dnZdO/e/WfnNGvWbNf+AyksLOTjjz9m9OjRlT6/iFSdtZYvHlrEZ7e6N3VndEvmV+8eQ0a3FI+TxZbCeStZ+/d3AGh6RG/aXHScx4lqxpYVEfrwFuzGeQD4B1+I/7CL1QIlUg0q0KXaBg0axDfffFOlY9u0qXxp7+uuu46JEyfy3HPP0atXr5/t39838qp8k58/fz6hUIgBAwZUKaeIVC5cHuH9//uOOc+vAqDj0c05782jaJLRMG9YrK/KNuWw4t6XIeKQmNWKTtf9AuNreKPNtrSQ0HvXYLctA8B/1B8IHHK2x6lEGh4V6NGQ3Nwdta4vkpvX7uWTk6tc+FbW4nLLLbfwyCOP8Pjjj1c6e0pGRkalo+Q5OTnATyPp+7PzRt+azjojIlCSW87/zvmSVdO2AnDIBR0Z8+/DCMRrppZoihSXsvzOF4kU7CDQNIkud16Iv0nD+wXIWofwZw+4xbnxERhxPf6eJ3gdS6RBUoEeBcbnr9WWkvrmYFpcbr/9du6//37Gjx/PFVdcUek5ffr04a233sJxnD360HfOmd63b98DPm/Llu4iKWvXrq1SThHZU86KIl4e8znblxQCMOLOvhx7S2+1KUSZjTisGv8/StdswQT8dL7tV8S3Svc6Vo1E5vxvV895YMR1Ks5FDoIKdKm2mra43HXXXdxzzz3cc889+5xVBdw50J955hnee++9PWZyefHFF+nRo0eVbhAdMGAAHTp04JVXXuH2228nLS3tZ8dMmTKFkSNHVul1iDQma7/cxqtnfcmO7WX443yMffZw+p2naUZrw4bnPyZ/1mIAOvzxdFL6dvI4Uc04G+YSmfU0AL7ep+DveaLHiUQaNhXoUm0pKSkMHjy4Wuc88sgj3HnnnZx66qmMHj2amTNn7toXHx/PwIEDd3188sknM2LECC655BKys7Pp1KkTL7zwAjNmzOCdd96p0vP5/X6eeuopxowZw6BBg7jmmmvo2bMnoVCIOXPm8NJLL9GnTx8V6CJ7mf/aGiZeMptwmUOTzHjOf+soOhxVu21zjVX25O/Z8sbnALQYexSZJx7mcaKascXZhD65u2KF0G4EhlX+11ERqToV6FIn3nvvPQDef/993n///T32dezYcY8FkYwxTJw4kZtvvpmbb76ZvLw8evfuzdtvv81pp51W5ec84YQTmDVrFg8++CD33HMP27dvJyMjg44dO3L66aczbty4qLw2kVhgrWX6fQuZeucCADJ7pPCrd4+hWRet9lgbihatZc3jbwOQOqgb7X57kseJasY6EUKf3gMluRCXRPCEOzCBOK9jiTR4Wkl0N1pJVCqjr7vEunBZhHd//y0//Hc1AJ2Gt+DcN44iMV2FVm0o35bHoiv+Tji3iPh2zen518sJJCd6HatGwl//h8icVwAInHQP/k7DPE4k0nBoJVEREanUjuwyXjv7S9Z84S62NvDiTpz6j0EE4jRTS22IlJaz/M7/Es4twp+cSNe7LmywxXlk9Ve7inP/gHNVnItEkQp0EZFGKntZIS+P+ZzsZUUAjLq3H0ff0EsztdQS6zisfvgNSlZsBJ+Pzjf/koS2mV7HqhFbsInwZw8AYFr3xz/ktx4nEoktKtBFRBqhoq2lPDdyCoWbSgnE+xj7/BD6ntPB61gxyzoOa5+YSN4Mt8e//e9PIfXQrh6nqhkbLif08R1QVgSJ6QSPvw3jVzkhEk36HyUi0shYa3nnd99QuKmU+NQg4yYdQ/uhDXMktyGwEYc1j71F9uTvAWhx+pE0P+0Ij1PVXPjLv+9ajCh4/G2YJP3bEYk2FegiIo3MN/9ewdIPNgJw6t8HqTivRTYSYfXDb5Az9QcAWpxxJO0uO7XBthFFlnyK8+O7APgP/zW+tgMPcIaI1IQKdBGRRmTbogI+vnYuAP1/2ZH+52t2otpiwxFWjf8fuZ+7qyC3PPsY2l5yYoMtzp2cVYSnPwqAr8MQ/If+0uNEIrFLBXo1+Hw+QqGQ1zGkjkUiEYLBoNcxRA5auCzCm+O+JlwaIa1jE0558lCvI8UsJxRm1YOvkffljwC0Om84bS46vsEW57Z8B+GP7oRwKaS0JDD6ZozxeR1LJGbpf1c1JCQkUFZWRk5OjtdRpI7k5ORQVlZGQkKC11FEDtpnt81n89w8jM9w5gtDSWiqec5rg1MeZuV9r+wqzlv/alTDLs6tJTzlIWzeWvAFCB5/JyYh1etYIjFNI+jVkJmZSVlZGVu2bCEvLw+/X/MEx7JIJEJZWRkpKSlkZqpHVxq2lZ9t4atHlwBwzE296DisuceJYpNTHmLFPS9T8I37uW5z4XG0/uVIj1MdnMjc13BWfg5A4Og/4WvZ0+NEIrFPBXo1GGNo27Yt27dvp7S0FMdxvI4ktSgYDO4qzhvqyJcIuIsRTfjNLADaHd6MY2/t43Gi2OSUhVhx138p+H4ZAG0vOZFW5xzrcaqD46z/jsjMpwHw9TwRX+/TPE4k0jioQK8mYwzNm2vkSUQaBmst713+LQUbSohLDnDWf4/AH1R3Y7RFSstZceeLFM5dAUC7S0+h5ZkNe2VNW7iF0Cf3gHUwzbsTOOZKDVaI1BEV6CIiMWzO86tY+PZ6AE5+/FCadUn2OFHsiewoY/kdL1A0fxUA7f/vNFqMOdLjVAfHhssJfXQHlOZDQirBE+/CBOK9jiXSaKhAFxGJUdnLCvnwyjkA9Dm7PQMuzPI2UAyKFJey7LbnKV64BoAOfzqD5qcM8TjVwQt/8QR22xJ3MaLjbsOktPI6kkijogJdRCQGRUIOb42bSXlxmNR2iZz6j0FqT4iycFEJy299juLF68AYOv55LJknHuZ1rIMWWfg+zqIPAPAf/ht87Qd7nEik8VGBLiISg6bd/SMbvs3BGDjzuSE0aab2hGgKF+5g2c3PsmPZBjCGrKvPIuO4QV7HOmjOlkWEP38CAF+no7UYkYhHVKCLiMSY1Z9v5YsHFwJw1LU96TSipceJYks4v5ilNz9DyYpN4DN0uu4XNBsxwOtYB82W5BH6+E5wQpi09gRG3aC/uoh4xNNb+Y0xycaYJ4wxm4wxJcaYb40xY6p47lnGmK+MMbkVj6+NMb+o7cwiIvVZSV45b188C2uh9aHpjLirr9eRYopTWs6yW5+rKM59dL7x/Ngozp0IoU/uhqKtEEggcOLdmLgkr2OJNFpez7U1AfgVcCtwCrAQmGCMOXl/JxljLgLeBDYCv6x4bAD+Z4z5Ta0mFhGpp6y1vP+H78hfu4Ngop+zXhxKIE4LqkWLdRxWPfyG29biM3S++XzSj+nndayoiMz8D3aDe0NxYNSN+JpleRtIpJHzrMWloggfDZxprZ1QsW0q0Bl4BJi0n9N/DawBfmGtdSrO/RhYCVwIPFuL0UVE6qV5r6xhwf/WAnDCIwNo3lPLsUfTxhcnkzdjAQDtLzuV9GGx8deJyIrpROb+DwD/gHPxd2nYiyuJxAIvR9DHAvnAOzs3WGst8ALQ0xjTez/nhoCincV5xbkOUASU1U5cEZH6K3dVER/88TsAepzWhsG/6+JxotiS/dkcNr82FYDmpw6h+ZgjPE4UHU7OGsJTHgLAtB2If+jvPE4kIuBtgd4XWLh7kV1h3m779+VvQC9jzC3GmExjTHNjzC1AD+CxWsgqIlJvRcIOb104k7LCMMmtEjj9qcN0c18UFf24mjV/fQuAlIFdaf/702Li82sjYcKT74NQCSS3IHj87RifWqJE6gMvZ3HJAJZWsj1nt/2Vsta+U3Ez6UvAvRWbi4FzrLUf7es8Y0zeATI1PcB+EZF654sHFrHu62wAxj5zOEnNEzxOFDvKNuey4u6XsKEICe2b0/mWX2ICsVHERua8ht2+DDAEj7sVk5jmdSQRqeD1TaK2JvuMMccBrwD/A44HTsK94fRVY8wpUU0oIlKPrft6O9Pv/RGAoVd0p+sJrT1OFDsixaUsv+MFwvnF+FMS6XLXhQSSE72OFRVOzioi374IgL//Wfhax8bNriKxwssR9GwqHyVvVvE2p5J9GPfvii8AU6y1v99t10fGmHbAk8AHlZ1rrU3bX6CKEXaNootIg1BaEOKtC2fiRCwt+zVl9P39vY4UM2zEYeWDr1G6Zgsm4KfLbReQ0CbT61hRYZ0I4SnjwQlBahv8Qy7xOpKI7MXLEfQfcfvI986w89f4Bfs4ryXQGvi2kn3fAp2MMfr7rojEvA+v/J7cVcUE4n2c9d+hBBNio/WiPlj/9CQKvlkCQIc/nUFK/84eJ4qeyA+vY7cuBiA44jpMUD8yReobLwv0CUAacNpe2y8EllhrF+7jvFygFDi8kn1DgWxrbWm0QoqI1EcL3ljL3BdXA3DcQ4fQsm+ap3liybYPZrF1wpcAtDzraDJPGOxxouhxctcSmf0cAL6+Z+BrO8DbQCJSKS9bXCYBU4FnjDEZwCrgImAYcPrOg4wx04BjrbUGwFpbZoz5F3ClMeZp3AWL/LiF/TDcRY9ERGJW/rodvHe5+0fErie0YsgfunmcKHYUzFnO2r+/C0DTob1o+5sTPU4UPbtaWyIhSGlF4IhLvY4kIvvgWYFurbXGmDOA+yseabgriZ5prX3vAKdfCywGLgXOBhzcGWHGAS/XUmQREc85EYe3L5pJaV6IpObxnPHM4TEx5V99ULpuGyvvexkch8TOrel0w7kYv9dzKURPZP7b2C3uDcXBEddigrFxw6tILPJyBB1rbQHwx4rHvo4ZXsm2CPDvioeISKPx5SNLWP35NgBO/8/hpLRSkRUN4cIdLL/zBSJFpQTSk+l654X4E+O9jhU1Nn8DkVnPAODrfSq+doM8TiQi+xM7QwMiIjFu43c5TLl9PgCH/b4LPU5t43Gi2OCEwqy892XKNmRjggG63jGOuBZpXseKGmsdQlPHQ7gMklsQOPL3Bz5JRDylAl1EpAEoLw7z5riZOGFLZs8Ujh8/wOtIMcFay7q/v0vhDysByLrmbJJ6dvA4VXQ5C97BbnQX6Q4OvwYTl+RxIhE5EBXoIiINwEfXziF7aSH+oI+z/3sEcU087VCMGVvfnsH2j74BoPUFo2g2/BCPE0WXLdhE+OunAPD1PBFfh8omQBOR+kYFuohIPbfonfV89x93hHfUvf1oPTDd40SxIW/WItY//SEA6cf2p/WvRnmcKLqstYSm/gXCpZCUSeCoP3gdSUSqSAW6iEg9VrCxhHcvdUd4O41owRFX9fA4UWzYsXITqx58DawlqWd7sq4+O+Zmw3EWvofdMAeAwLFXY+KTPU4kIlWlAl1EpJ5yHMvE38xiR3Y5ielxnPn8EHy+2CoivRDKLWTFnS/ilJQTbN6ULrePwxcf9DpWVNnCLYS/cic683U/Dn/WER4nEpHqUIEuIlJPzXxiKSsmbwFgzL8Hk9q2iceJGj7rOKy8/1XKt+bhS4ij610XEWyW4nWsqLLWEpr2MIR2QGI6gWH7nMlYROopFegiIvXQ5h9ymXyzO/PGob/pRO8z23ucKDZs/+hbiuavAqDTDefSpHNrjxNFn7P4Q+w6d6XZwLFXYRJSPU4kItWlAl1EpJ4JlbhTKkbKHZp1TebERwd6HSkmhHIL2fCse1Nos1EDSTuit8eJos8WbSP85T8A8HUdgb/z0R4nEpGaUIEuIlLPfHLjD2xbWIAvYDj7v0OJT46t/mivrH/qAyJFpfhTEmn3u5O9jhN1u2ZtKS+GxDQCR1/hdSQRqSEV6CIi9cjSSRuZ/fflAIy4oy9tD8vwOFFsKPh+GTlTfwCg3SUnEUyLvRlNnEWTsOvcGX8Cx1yJSUzzNpCI1JgKdBGReqJoSykTfzsbgI7DmjPs+p4eJ4oNTlmItU9OBCC5TxYZxw/yNlAtsIWbf2pt6TYSf5djPU4kIgdDBbqISD1greWd382meGsZCU2DnPnCEHx+fYuOhk2vTaVsUw4m4KfDFWdgfLH1eXVbW3abtUWtLSINXmx9lxIRaaCWTdrE0kmbADj174NI65jkcaLYULJ2K1ve+ByAlmcfTWLHlh4nij5n4XvY9d8BEBh+NSahqceJRORgqUAXEfGY41gm3+ZOqdhldEv6ndfR40SxwVrL2icnYsMR4lo3o/X5I72OFHW2YBPhL/8JVCxI1GmYx4lEJBpUoIuIeGzB/9ayZV4+AKPu7e9xmtiR/el3u+Y87/jH02NvtVDruLO2hEuhSQaBYX/yOpKIRIkKdBERD4XLI0y5Yz4Avc9sR9vBzTxOFBvC+cWsf9qd8zz92P6kDurucaLocxa8g90wB4DA8GswCbG1IqpIY6YCXUTEQ98/u4rclcUYn2Hk3f28jhMz1v9nEpGCHfiTEmh/2alex4k6m7+B8NdPAeDreSL+rCM8TiQi0aQCXUTEI+U7wnx+348ADLgoi+Y9tSR7NBT+sJLsyd8D0PY3JxJsFlsjy9Y6hKaMd1tbkpoTOOoPXkcSkShTgS4i4pFZf1tG4aZS/HE+ht/Wx+s4McEpD7OmYs7zpJ7tyTzpMG8D1YLIvLexm9ybioMjrsXEx96iSyKNnQp0EREPlOSWM2P8IgAOv7wraR00rWI0bH5jOmXrt4HPR4crxsbcnOdO3jois54GwNfrFHwdDvc4kYjUhtj6ziUi0kB8+fBiSvNCxCUHOPrGXl7HiQml67ex+bVpALQ8cxhNOrf2NlCUWSdCeMpDEC6D5BYEjrrc60giUktUoIuI1LHCzSXMfHIpAEde1YOk5gkeJ2r4rLWs/ds72FCYuBZptL5glNeRoi4y7y3sZveeheCI6zFx+quLSKxSgS4iUsc+v38hoR0RmmTEccRVPbyOExNypsylcO4KADr84XT8CXEeJ4ouJ3ftT60tfcbgaz/I40QiUptUoIuI1KGclUV8+5RbSB59U28SUmNr8RwvhAt3sP6pDwBIG9aXpkN6epwoutzWlgchEoKUVgSOuMzrSCJSy1Sgi4jUoal3LcAJW5q2b8Jhv+/qdZyYsOGZjwjnF+NrEk/7y0/zOk7URX54HbvFvaE4OPJ6TFwTjxOJSG1TgS4iUke2zM9j/itrADj2tj4EE/weJ2r4ihasZvtH3wDQ9qLjicuIrbnkneyVRGY9B4Cv31h8bQd6nEhE6oIKdBGROvLZ7fOxFjJ7pDDgwiyv4zR4TuinOc+bdGtL81OHehsoymwkRHjy/eCEILUNgaG/8zqSiNQRFegiInVg7VfbWfLeRgBG3t0Pf0Dffg/WlrdmULpmC/gMHa8Yi/HH1uc08s3z2OwVYHwER92ECSZ6HUlE6khsfTcTEamHrLV8dqu78mObQen0PrOdx4kavrJNOWx65TMAWpx+JE26tfU4UXQ5mxYQmfMaAP6B5+Fr3dfjRCJSl1Sgi4jUshWfbGb159sAGH1vf4wxHidq2Nw5zydiy8MEM5vSZtxxXkeKKhsqIfTZA2AdTGZX/Idd7HUkEaljKtBFRGqR41gmV4yedxregs6jW3qcqOHL/Xw+Bd8tA6DD/52Gv0m8x4miK/zlP6FgI/iCBEbdjPFrKk6RxkYFuohILVr41jo2zckDYNS9/TR6fpDCRSWs+9f7ADQ9ojdpR/bxOFF0RdbMxFn4HgD+IZfgy+jkcSIR8YIKdBGRWhIJOUy5fT4APce0pf3QTI8TNXwbn/+YcG4hvoQ4OvxfbM15bkvzCU/9CwCmTX/8h5ztcSIR8YoKdBGRWjLnhVVkLyvCGBh1Tz+v4zR4xYvXsu2D2QC0GTeauOZp3gaKImst4emPwY4cCCYSHHkTxqd58kUaKxXoIiK1IFQSZvo9PwLQ/4IsWvRp6nGihs0JhVnz+ASwlsQurWlxxpFeR4oqZ9lknBXTAQgM+yMmtZXHiUTESyrQRURqwex/LqdgQwn+oI8Rd2iKvIO1+fXplKzavNuc57EzumwLtxL+/HEAfFlH4ut5kseJRMRrKtBFRKKsNL+cGQ8tAmDwZV1Iz0ryOFHDVrJqM5tfnQpAy7HDSOrR3uNE0WOtQ2jqQ1BeDIlpBIZfqxuJRUQFuohItH316BJ2ZJcTlxTgmJt6ex2nQbORCKsffRMbjhDfNoM2F8bWnOfO/InY9d8DEDj2GkyTdI8TiUh9oAJdRCSKiraW8vVflwIw9IruJLdM8DhRw7bl7RnsWLYBjCHrqrPxxcfOnOBO7lrCX/8bAF/PE/F3HuZxIhGpL1Sgi4hE0RcPLqS8OExisziOuraH13EatNJ129j44mQAmp82lOS+Wd4GiiIbCROefD9EyiGlJYFhf/Q6kojUIyrQRUSiJG9NMd/8awUAw67vRULTOI8TNVw24rD6sbewoTBxLdNp++sTvI4UVZHvX8ZuWwIYgiNvxMTpPgUR+YkKdBGRKJl69wIi5Q4pbRIZ8oeuXsdp0La+9zXFC9cA0PHKM/EnxnucKHqcLYuJfPsiAP5DzsbXdoC3gUSk3lGBLiISBVsX5vPDf92C8thbexNMDHicqOEq25jNhuc+BiDzpMNIHRg7v+zYcBnhz+4H62CaZeEf8luvI4lIPaQCXUQkCqbcsQDrWJp1TebQX3f2Ok6DZR2HNY+/jS0LEcxMpd1vT/Y6UlRFvn4Km7cOfH4Co27GBNQGJSI/pwJdROQgrZ+dzaIJ6wEYeVc//EF9a62p7R9+Q+EPKwHo+Kex+JNiZxYcZ/33ROa/DYB/8MX4mnfzOJGI1Ff6KSIicpA+u3U+AK0GpNHnnNhZRKeulW/NY/0zHwLQbNRAmg7p6XGi6LGhEkJT/wKAadkL/6Hne5xIROozFegiIgdhxeTNrJyyBYDR9/TH59MqkDVhrWXN42/j7CgjkJ5M+9+f6nWkqAp//RQUbgZ/kMDIGzE+v9eRRKQeU4EuIlJD1tpdo+cdhzWn64mtPE7UcGV/+j0F3y0DoMMfzyCQ0sTjRNHjbJiDs2AiAP7Df4MvvYO3gUSk3lOBLiJSQ4smbmDDtzkAjL6vH8Zo9LwmyrMLWP/U+wCkH9OP9KP6eJwoen7W2nLIOR4nEpGGQAW6iEgNRMIOU253R8+7n9KGDkc19zhRw2StZe2TE4kUleJPbUL7/xvjdaSoCs/8DxRscltbRtyg1hYRqRIV6CIiNTDvpdVsW1SAMTDqnn5ex2mwcqfPI3/mIgA6XH4awbRkjxNFj7NhLs78CQD4D/s1vmYdPU4kIg2FCnQRkWoKl0WYevePAPQ9rwOt+qd5G6iBCuUVse6f7wHQdGgv0ocf4nGi6NmjtaVFT/wDfuFxIhFpSFSgi4hU0zf/XkH+2h34AoaRd2r0vKbW/fM9wvnF+JMS6PCnM2Kqhz8y82ko2Ai+IIGR16u1RUSqRQW6iEg1lBWG+OKBhQAM+m1nmnWJnZaMupT75Y/kTp8HQLvLTiEuI9XjRNHjbPzhpwWJDrsIX7NOHicSkYZGBbqISDV8/delFG8rI5jo59hbYme2kboULtzB2r9NBCB1UDcyjhvkbaAosqFSQlPGA2Cad8c/8DyPE4lIQ6QCXUSkioq3l/HVo4sBGPLHbqS0TvQ4UcO07t/vE84twpcYR8c/nxlbrS2zdm9t0YJEIlIzKtBFRKpoxkOLKCsMk5AWZNj1vbyO0yDlz15MzuQ5ALT77cnEtUjzNlAUORvnEZm3s7XlQnwZam0RkZpRgS4iUgX563Yw+x/uSpdHXduTxPQ4jxM1PJHiUtY8MRGA5P6dyTzpMG8DRZENlRKeOh6wFa0t53sdSUQaMBXoIiJVMO3eHwmXOSS3TGDon7p7HadBWv/0h4S252Pig2RdeSbGFzs/giKzn8HmbwBfgMBILUgkIgcndr47iojUku1LC5n7/CoAjr21N3FJAY8TNTwFc5az/cPZALS9+Hji22R4nCh6nE3zifzwFgD+wRfiy+jscSIRaehUoIuIHMCUO+bjRCzpnZI49BIVX9UVKSljzV/d3uykXh1oMeZIjxNFjw2XEZ5S0dqS2U2tLSISFSrQRUT2Y+P3Ofz4xjoARtzZl0CcWheqa8NzH1O+JRcTDJB19VkYf+z86InMehabvx58fre1xa+/rojIwYud75IiIrXgs1vnA9Cib1P6ndfB4zQNT+GCVWx792sA2lwwioT2LTxOFD3OpgVEfngDAP+gcfgyu3icSERihQp0EZF9WDVtK8s/2QzA6Hv64Yuhkd+64JSFWPOY25vdpFtbWp59tMeJoseGy36atSWzK/5Df+V1JBGJIfppIyJSCWstk29xl6JvPzSD7qe28ThRw7Pxv5Mp25CNCfjJuvpsjD922oMis5/D5q1Ta4uI1AoV6CIilVjy/kbWz8oGYPR9/WNqtcu6UPD9Mra8/QUArc4bTmKnVh4nih5n8497tbZ09TiRiMQaTwt0Y0yyMeYJY8wmY0yJMeZbY8yYKp5rjDGXGmO+M8bsMMbkGWNmGmNiZ3oAEfGEE3H47Da397zrCa3IOjZ2+qbrQtGitay4+yVwLImdW9Pq3OFeR4oaGyl3Z22xDiajC/5Df+l1JBGJQV7/TW4CcChwPbAKuBiYYIw5zVo76QDnPg2cBYwHvgKSgEEVb0VEamz+q2vZuiAfgNH39vc4TcNSsmozy297Dqe0nLiW6XS9+yJ8Qa9/1ERP5LuXsXlrwfgqWluCXkcSkRjk2XdNY8zJwGjgTGvthIptU4HOwCPAPgt0Y8xZuMX8MGvt17vt+qDWAotIoxAujzD1rgUA9DmnPa0HpnucqOEo3bidpTc/S6SolEB6Mt0e+A1xmU29jhU1TvYqIt+/AoB/wLn4mnfzOJGIxCovW1zGAvnAOzs3WGst8ALQ0xjTez/n/gn4fK/iXETkoH339EpyVxXj8xtG3tXP6zgNRvn2fJbd9Czh3EL8yQl0u+83JLTJ9DpW1FgnQnjaX8AJY5q2xX/YRV5HEpEY5mWB3hdYaK119to+b7f9P2OMCQJDgfnGmPuNMVuMMWFjzI/GmP1+x6zoU9/nA4idoR4Rqbby4jCf37cQgIG/7kRm9xSPEzUM4YJilt38LOVbcvHFB+l698U06dza61hRFVkwEbtlEQCB4ddgAvEeJxKRWOZlY2AGsLSS7Tm77d/XefHARcB64I9AHnAJ8LwxJs5a+5/oRhWRxmDmE0sp2lJKIN7Hsbf28TpOgxDZUcayW5+ndO1WTNBPlzvGkdy7o9exosoWbiYy82kAfL1Owdd2oMeJRCTWeX3njq3Bvp2j/gnAydbaNQDGmMm4/eu3A5UW6NbatP2F0Si6SOO1I6eMLx9eDMDh/9eNpu2aeJyo/nPKQ6y460V2LF0PPkOnG84j9dDY6su21hKe/lcIl0KTZgSOuMzrSCLSCHjZ4pJN5aPkzSre5lSyDyAXt3hfvLM4h1396x8B7YwxmhNNRKrly78spjQ/RHxqkKNv7OV1nHrPhiOsvP9VCn9YCUDHK88kfVilnYkNmrPsM5y1swAIHH0FJkFtTyJS+7ws0H8Eehlj9s6w866sBZWdZK0tAZbv45o7VxLZu69dRGSfCjaWMPPJZQAceXUPmmSov3h/rOOw+tE3yZ/p9mS3u/QUMo8f7HGq6LMl+YRn/A0AX6ej8Xc51uNEItJYeFmgTwDSgNP22n4hsMRau3A/576NW9xn7dxg3GX+TgJWWmu3RzeqiMSy6ff9SLg0QlLzeI74c3ev49Rr1lrW/et9cqbMBaD1L0fS8sxh3oaqJeEv/w6l+RCXROCYK7yOIyKNiJc96JOAqcAzxpgM3IWKLgKGAafvPMgYMw041lq7+zrbfwF+BXxkjLmLn24SHQScVxfhRSQ25Kwo4vtn3DaNY27uTXyKFp7Zn00vTWbbu+4Mt83HHEHrcaM9TlQ7nLWzcZZ+CkDgiMswSbEzZaSI1H+eFejWWmuMOQO4v+KRBizEXbjovQOcm22MORq3UP8HkAjMB8ZaayfWYmwRiTFT7piPE7akdWzC4Eu7eB2nXtsyYQabXp4CQLNRA2n/+1Nx/3gZW2yohND0RwEwbQ7B1/sUjxOJSGPj6Swu1toC3GkS/7ifY4bvY/tq4JxaCSYijcLG73KY/9paAEbc0ZdAvN/jRPXX9k++Zf2/3cWamw7tRdZVZ2F8XnZJ1p7IrGehcAv4g+6c5z+7VUpEpHbpu46INErWWj69yV0XrWW/pvT/VWzN3R1NuTMWsOavbwOQckhnOt98PiYQm7/MOFsWEpn3FgD+wRfhS2vvcSIRaYxUoItIo7Tik82snLIFgOPuPwSfX98OK1Pw/XJWPfQaOJYm3dvR5Y4L8cXFZp++jYQIT30YsJiMLvgHnOt1JBFppPQTSUQaHcexfHqzO3reaXgLup7YyuNE9VPRorWsuPu/2FCEhA4t6HbvxfibxO4UlJE5r2FzVoHxERhxHcbv9Vp+ItJYqUAXkUZn/itr2PxDHgDHPdA/Jm90PFglqzaz/LbncErLiWuZTrf7f0MgNcnrWLXGyV1L5Nv/AuDvfza+Fj08TiQijZkKdBFpVEKlET67fT4Afc5pT9vDKlvQuHEr3bidpTc/S6SolEB6Ct0fuIS4zKZex6o11jqEpz0MTghSW+M//GKvI4lII6cCXUQaldn/WEb+2h34AoZR9/Q78AmNTPn2fJbd9Czh3EL8yYl0v/83xLeJ7V9inB/fw25yf2kLHnsNJpjocSIRaexUoItIo1GSW84XD7jL0w++rAsZXVM8TlS/hAuKWXbzs5RvycWXEEfXey4msVNs9+fbom2Ev34KAF/PE/G1H+RxIhERFegi0oh88dAiSnLLiUsOcOwtfbyOU69EdpSx7NbnKV27FRP00+X2cST36uB1rFplrSX8+V8htAMS0wkc+XuvI4mIACrQRaSRyF+3g1lPLgXgqGt7ktwiweNE9YdTHmLFXS+yY+l68Bk63XAeqYd29TpWrXNWTMdZ/RUAgWF/xCTEbp+9iDQsKtBFpFGYetcCwmUOya0SOPIqzdCxkw1HWHn/qxT+sBKAjleeRfqwvh6nqn22tIDwF08A4Ot4BL6uIzxOJCLyExXoIhLztizIY+6LqwEYfnsf4pI0vzWAdRxWP/om+TPdvvx2l51C5vGNowc7/NU/oSQXgk0IHHuVptoUkXpFBbqIxLzJN8/DOpbMHikc+pvOXsepN9Y/9QE5U+YC0PpXI2k5dpi3geqIs+47nMUfARA44neY5OYeJxIR2ZOGkUQkpq2evpWlkzYBMPq+/vgDGpcAyPl8Hlsnuv3XzcccQesLRnucqPZZa3EWTXJHzwHTqg++PmM8TiUi8nMq0EUkZllr+eTGHwBof0QGPU9v63Gi+qFsSy5rH58AQOqgbrT//akx3+Jh8zcSmvYwdsMcd0NiGoER12OMfmETkfpHBbqIxKyFb61nwzc5ABz/4CExX4RWhY04rB7/OpHiUgJpyWRdew7GF7tFqnUiROZPIDLrGQiXAuDrNorAsD9hEjVri4jUTyrQRSQmRUIOk2+dB0DPMW3pcJT6jAE2vzaVoh9XA5B19VkE02N3sSYnZzXhqX/BblnobkjKJHDs1fizjvA2mIjIARywQDfG/AN41lr77W7b4qy15bWaTETkIHz7nxXkLC/C+Ayj7+vvdZx6oWjhGja+PAWAFqcfSdPDe3qcqHbYSJjInFeJfPtfcEIA+HqfSuCIyzDxyR6nExE5sKqMoP8emAF8C2CMyQC2GmOOs9ZOqc1wIiI1UVYYYvo9PwJw6G860bxXqseJvBcpLmXVQ/8DxyExqxVtLznR60i1wtm2lPCU8djsFe6G1DYER1yLr+1Ab4OJiFRDTVtc1MgpIvXWl48soXhbGcFEP8Nvj/1Fd6pi7d/eoXxLLiYuQKcbz8MXF/Q6UlTZcBmRb14gMvd/YB0wPvz9z8J/+G8wQa0aKyINi3rQRSSmFG4u4evHlgAw9M/dSW2T6HEi72V/NoecqXMBaH/pKSRmtfQ2UJQ5m+a7veZ56wAw6R0JjLweX8veHicTEakZFegiElOm3fMj5cVhmmTEMey62Oyxro6yTTms/fs7ADQd2ovMU4Z4nCh6bPkOwrOexpk/EbDg8+M/9Jf4B12A8cd5HU9EpMaqWqDbKm4TEfHM9qWFfP/0SgCOuaUPCU0bd5FmwxFWPfQazo4ygs1SyLrqrJiZatJZ+w2haY9A0RYATPPuBEZchy+zq8fJREQOXlUL9AeNMTdVvO/HLc6fNsYUV3KstdYeEpV0IiLV8Nmt83AilvROSRx2WRev43hu48ufUbzYbfvIuvYXBJomeZzo4NnSQsJf/QNn8UfuBn8Q/+G/xn/ILzA+v7fhRESipCoF+lrcgjxlr22+vbaJiHhm3cztLHx7PQCj7ulHIL5xF2uF81ay+bVpALQ8+xhSD234I8uRNTMJT/0L7HAXnzKt+7mj5mntPU4mIhJdByzQrbVZdZBDRKTGrLV8eqO7KFGbQen0+UUHjxN5K1xYwqrxr4O1NOnahjYXHed1pIMWWfyRW5xbB4KJBIb+Dl/f0zEmdldBFZHGSzeJikiDt/T9jayZsQ2A4+4/BJ8vNvqsa8Jay5rH3ya0PR9ffJBON5yHL9iwv9WH571FZMbfADDNexA88U5MSiuPU4mI1J5qf9c2xiQCnYFUoABYaa0tiXYwEZGqiIQdPr3FHT3venwrOo+KrSkEqyv7k+/Im7EAgPaXn0ZC++YeJ6o5ay2R7/5LZPZzAJg2hxA8+T5MXMPvpRcR2Z8q/23QGHOEMeYTIA+Yh7u66Dwg1xjzkTHm8NqJKCKyb3NfXM22hQUYA8c90N/rOJ4qXbeNdf94F4C0YX3JOGGwx4lqzlpL5Kt/7irOfR2HEjz1IRXnItIoVGkE3RhzLvAiEATW4BbmBbij6P2B44ERxphfWmvfqqWsIiJ7KN8RZuqd7mhx/191pNUh6R4n8o4TCrtTKpaFCDZvSsc/n9lgp1S0ToTw9EdxFk0CwNd1JIFRN2H8DbtVR0Skqg743c4Y0xJ4CtgC/Npa+1klx4wGngOeMcZ8Ya3dGvWkIiJ7mfnEUgo3luCP8zHyrn5ex/HUxhc+ZcfyjWAMna77BYGUhrmCqo2ECE++D2fFdAB8vU8jcMyfNYWiiDQqVWlx+Q2QCJxUWXEOYK2dDJwMJAEXRy2diMg+FG8vY8b4xQAM+UM30jo23taHgu+Xs+XNzwFode5wUvp39jhRzdhQKaEPb91VnPsHnkfg2KtUnItIo1OVAn0k8JG19sf9HWStnQ98CDT8+bxEpN77/P6FlBWESEgLcvRNvbyO45lwfjGrH34dgKSe7WlzwSiPE9WMLSsi9P712LWzAfAP+S2BIy5rsG06IiIHoyoFem/gyype78uK40VEak3uqiK++edyAI6+oRdNmsV7nMgb1lpWP/YWoZxCfE3i6XTDeZhAwxtttiV5hN65GrtpPgCBY/5MYNCvPE4lIuKdqtxxk4bbf14VW4DGe5eWiNSJKXcsIBJySG2XyJA/dvM6jme2fzCL/JmLAOjwf2OIb93M40TVZ4u2EXr3WmzeWjA+AiNvxN9Df4gVkcatKgV6IlBexeuFgMY5lCUidWLTnFzmvbIGgJF39iWY2Dhn9ihZs4V1T30AQLMRh9Bs1ECPE1Wfzd9A+bvXQOEW8AUJnHA7/k7DvI4lIuK5qv5ks7WaQkSkij696QcAWvRtyiHjsrwN4xGnPMSqB1/DloeJa5lOhz+e0eB6tZ3slYTevRZKciGQQPDke/G1G+R1LBGReqGqBfozxph/R/F6IiLVtuLTzayY7HbcHXdff3z+Kq+1FlM2PPsRJas2g89HpxvOxZ+U4HWkanE2LyT0wY1QVgjxyQRPeQhfK92+JCKyU1UK6s/RCLqIeMxxLJ/e7I6eZx3TnG4nt/Y4kTfyv1nC1olfAdD6VyNJ7t3R40TV46z/ntCkWyBcConpBE/7C77MLl7HEhGpVw5YoFtrh9dBDhGR/Vrw2ho2zckD4LgHD2lwLR3REMotZPUjbwKQ3CeL1ucN9zZQNUVWzSD8yd0QCUFKS+LGPIJp2tbrWCIi9Y5aUkSk3guXRfjs9gUA9Dm7Pe0Oz/A4Ud2zjsPqR94knFeEPymBrOt/gfE3nCkVI0s+JTzlQbAOJq09wdMexqS08DqWiEi9dMAGTmOM3xjzoDHm9wc47nJjzP2mMQ5riUitmv3P5eStLsYXMIy6t5/XcTyx9d2vKfh2KQAd/jyW+JYNZ0bbyIKJhD+73y3OM7sRHPuEinMRkf2oyh1WFwDXAd8c4LjZwA3A+QcbSkRkp5K8cj6/fyEAgy/tQkbXFI8T1b0dKzex4ZkPAcg4bhDNjunvcaKqC3/3MuHPHwfAtO5H8PRHMYlp3oYSEannqlKg/wKYbK39bn8HVez/GBXoIhJFM8YvpiSnnLjkAMfe2sfrOHXOKS13p1QMRYhvk0H7y0/zOlKVWOsQ/vrfRGY9DYDpcDjBU8dj4pM9TiYiUv9VpUAfBEyu4vWmAoNrHkdE5Cf563cw8wm3reOoa3qS3KJhTScYDeufnkTp2q3gr5hSsUn9XwvOyVtHaOJVROa8BoCvy7EET7oXE2x8Xz8RkZqoyk2izYCtVbzetorjRUQO2tS7FhAujZDcMoEjrurudZw6l/f1Qra9PwuAthcdT1KP9h4n2j/rRIj88AaR2c9BxF2A2tf3dALD/oTxNZwbWkVEvFaVAr0QyKzi9TKAoprHERFxbf0xn7kvrAZg+O19iE8OehuojpVvz2f1o+6UiikDutDy7KM9TrR/TvZKwlPGY7ctcTcktyQ4/Gp8HQ73NpiISANUlQL9R+B44JEqHHtcxfEiIgdl8i3zsI4lo3sKh/6ms9dx6pR1HFY//AaRwhL8KYlkXXsOxlc/V021kRCR714m8v3L4IQB8PU9g8DQ32HimnicTkSkYapKgf428Igx5nRr7Tv7OsgYMwa3QL86WuFEpHFa88U2lry/EYDR9/XHH6yfxWlt2fLWDArnrgCg45VnEZfZ1ONElXO2LCY89S/YnJUAmKbtCIy4Fl+bQzxOJiLSsFWlQP83cDnwujHmYeA/1trVO3caY7KA3wLXAksrjhcRqRFrLZ/c+AMA7YZk0OuMxrXSZPHS9Wx4/mMAMk8+nPSj6t/MNTZcRmT2c0R+eAOsA8aH/5Bz8B/+a0yg/t/EKiJS3x2wQLfWlhhjTgHeB24CbjTGFAIFQAqQChhgCXCqtba0FvOKSIxbNGE962dlA3D8g4fQmNY+i5SUseqh1yDikNC+Oe0vPcXrSD/jbJznjprnrwfANOtEYMT1+Fr29DiZiEjsqMoIOtba5caYAcDvgLOBPkAr3CL9C+At4GlrbUkt5RSRRiAScph8yzwAepzaho5HN/c4Ud1a98/3KduQjQn66XTjefgS4ryOtIst30F45n9wFkx0N/gC+AddgP/QX2L8jesGXhGR2lalAh2gYmT8yYqHiEjUfffMSrKXFWF8htH3N5zVMqMhZ+pcsj/5FoC2vz6RJl3aeJzoJ87abwhNewSKtgBgmvcgMPJ6fBmN6+ZdEZG6UuUCXUSkNpUVhZh+jzsJ1MCLs2jRu37eGFkbStdvY80TEwBIPaw7Lc440uNELltaSPirf+As/sjd4I/Df/iv8R9yjuY1FxGpRSrQRaRe+OrRJRRtKSWY6GfEHX29jlNnnPIQK+9/FaeknGBmKp2u/UW9mFIxsnIG4c8fgx05AJjW/QiMuA5fWv1eLElEJBaoQBcRzxVtKeWrR9wFboZe0Z3Uto1n/ux1//6AkpWbwOej843nE2ia5GkeuyOX8BdP4KyY5m4IJBA44lJ8fU/HGO9/cRARaQxUoIuI56bd+yPlxWGaZMQx7PrGMxtIzvR5bP9gFgBtLjyO5L5ZnmWx1uIs+4zwjCehtAAA034wweHXYFJaeZZLRKQxUoEuIp7KXlbId/9xF+U55ubeJDStPzOX1KbSjdtZ8/jbAKQO6karXxzjWRZbtI3w9Mdw1nztbohPJnDk/+HreWKjmuZSRKS+UIEuIp6afOs8nLAlLSuJw37f1es4dcIpD7HyvldxdpQRzEgl6zpv+s6ttTiLPiD81b+gvBgAX6dhBI65EpOUUed5RETEpQJdRDyzflY2C99yF7wZdU8/AvGNY2aQ9U9NomTFRvAZOt14HsG05DrPYAs2E5r2F+z6790NiWkEjr4CX5fhGjUXEfGYCnQR8YS1lk9v+gGA1oem0/fcDh4nqhu5X8xn2/szAWhzwWhS+nWq0+e31sH58T131DzsLvzs6zaKwLA/YRIbz9SWIiL1mQp0EfHEskmbWP35NgCOu78/Pl/sj9qWbcxm9WNvAZByaFdanTu8Tp/fFmwiNPUv2A1z3A1NMggMvxp/Vv2Yd11ERFwq0EWkzjkRh09vmQdAl+Na0WV07M8S4pSHWXn/Kzg7ygikp9Dpul9g/HXTd26tg7PgHcJfP/XTqHmP4wkc9UdMQkqdZBARkapTgS4ide6H/65m64J8jIHjHujvdZw6sf7pSexY7vadd77xXILpdVMY2/wN7qj5RrediKRMAsdejT/riDp5fhERqT4V6CJSp0IlYabcuQCAfr/sSOsB6R4nqn25Mxaw7V13CsPWvxxJyiFdav05rXVw5k8kPPM/P42a9zyRwFF/wMTX/U2pIiJSdSrQRaROzXxyGQXrS/DH+Rh5Vz+v49S6sk05rNnZdz6gC63PH1nrz+mOmo/HbnTbiEjKJDD8Wvwdh9T6c4uIyMFTgS4idWZHdhkzHloEwOGXdyU9y9tl7WubEwqz8oFXiRSXEkhPptP159Zq37m1DpF5bxOZ9TSEywDw9TyJwFH/p1FzEZEGpO5XxtiNMSbZGPOEMWaTMabEGPOtMWZMNa9hjDFTjDHWGPPXWooqIlHw+QMLKc0PkdA0yDE39/Y6Tq3b8OxH7Fi6Hoyh0/XnEmxWe33nTt56QhOvJPLl393iPKk5wVMfIjjyehXnIiINjNcj6BOAQ4HrgVXAxcAEY8xp1tpJVbzG74CetRNPRKIld3Uxs/+xHIBh1/eiSUa8x4lqV95XP7J1wpcAtP7lCFIH1s4qqdaJEJn/NpGZT0OkHABfr1MIHPl7FeYiIg2UZwW6MeZkYDRwprV2QsW2qUBn4BHggAW6MaYtMB64BHiz9tKKyMGacsd8IuUOqW0TGXpFN6/j1KqyzbmsftT9lpTcvzOtfzmqVp7HyV1LeOp47OYf3Q3JLQgOvxZfh8Nq5flERKRueDmCPhbIB97ZucFaa40xLwBPGWN6W2sXHuAa/wQ+t9a+paWpReqvTXNzmf/KGgBG3NmXYKLXf7yrPU4ozKoHXyVSVEqgaRKdboh+37l1IkTmvUVk1jM/jZr3PtUdNY+L7b5+EZHGwMufkn2BhdZaZ6/t83bfv6+TjTHnAyOAKjeyGmPyDnCI1rkWqQWTb56HtdCiTyoDLszyOk6t2vj8JxQvXlfRd/4L4jJSo3p9J3ct4SkPYbdUfHtMaUlw+HX42g+K6vOIiIh3vCzQM4CllWzP2W1/pYwxmcDjwC3W2nW1kE1EomTlZ1tY/slmAEbfdwi+Olo90wt5Mxex5a0vAGh17nBSB3WP2rWtEyHywxtEZj8LkRAAvj5jCBxxGSauSdSeR0REvOf135ltDfc9gXtT6d+q9WTWpu1vf8UIu0bRRaLEcSyf3OSuYNnx6OZ0P6W1x4lqT/nWPFY//AYAyX2zaDMuen3nTs4awlMfwm5xp6gkpSXBEdfha6dRcxGRWORlgZ5N5aPkzSre5lSyD2PMccC5wEggda/e83hjTBpQZK0NRy+qiNTEj6+vZdP3uQAc/+AhxOq9IjYccec7Lypx+85vPA/j9x/8dZ0IkbmvE/nmuZ9GzfueTmDopRo1FxGJYV4W6D8CZxljfHv1oe9cWnDBPs7rgzt/+7RK9v2+4nES8FGUcopIDYTLInx223wAep/VjnZD9tm11uBteP4TihetBSDrunOIyzz4P8Q5OasITxmP3brY3ZDa2h01bzvwoK8tIiL1m5cF+gTc6RFPY7eZXIALgSX7mcHlTWBuJdunAm/htr3Mq2S/iNShb/69gtxVxfgChtH39vc6Tq3Jn7WYLW9+DkCrXxxL08E9Dup61okQmfMakW9eAKdi1LzfWAJDf4cJJh50XhERqf+8LNAn4RbVzxhjMnB7yi8ChgGn7zzIGDMNONZaawCsteuB9XtfrOJP5+uttdNqO7iI7F9pfjmf3+fOzT3od13I6FZ7K2h6qXxbHqsecfvOk3p3pM1Fxx3U9ZzsVe4MLduWuBtS21SMmg84yKQiItKQeFagV8x5fgZwf8UjDXdaxTOtte95lUtEDt6MvyxmR3Y5cckBht/ax+s4tcLtO3+NSMEO/KlN6HzT+TXuO7fWcec1//o/u0bN/f3OxD/0txo1FxFphDydxcVaWwD8seKxr2OGV/FasXn3mUgDU7BhBzMfd2dQPfLqHiS3TPA4Ue3Y+OKnFC90F1/qdO05xDWvWd+5Ld5O6LOHsOu/dTektiE48np8bQ6JVlQREWlgvJ5mUURizNS7fyRUEiGpRTxHXnVw/dj1Vf63S9j8+nQAWp5zDE0P71mj60RWziA87S9QWgCAr9fJBIb9UaPmIiKNnAp0EYmarQvzmfPcKgCG39aH+JSgx4mir3x7Pqv/8lPfeduLjq/2NWyohPCXf8dZ+IG7IT6FwPBr8Hc5NppRRUSkgVKBLiJR4UQcPvjjd1jHktEtmUG/7eJ1pKizkQirHnyNcH4x/uREOt94HiZQvb5zZ+sSwpPvw+a5iyCbdocSHHkjJrl5bUQWEZEGSAW6iETF139dyurPtwFw0qOH4g/6PE4UfRv/+xlFC1YDkHXtOcS1SKvyuT9Nn/gcOBHwBfAP/R3+Q87GmNj7XImISM2pQBeRg7Z5Xt6uRYkGX9aFbie19jhR9BV8t5TN/5sGQMuzjiZtaK8qn2sLNxP67AHsRneJBpPekcBxt+LL7FobUUVEpIFTgS4iByVcFuHti2YSKXdo1jWZE8YP8DpS1JVnF7Bq/OtgLUk929P21ydU+dzIsimEpz8K5cUA+PqeQeCIyzDB2JzdRkREDp4KdBE5KFPuWMCW+fkYn+HM54cQlxRb31b27DtPoNNN51ep79yWFxP+/HGcpZ+6GxLTCYy4Dn/WEbWcWEREGrrY+kkqInVq9edb+eqRxQAcc1Mv2g/N9DhR9G16eQpF892ZabKuOYf4lukHPMfZNJ/Q5PuhcDMAvg5DCIy8HtOkWa1mFRGR2KACXURqpLQgxIRfz8JaaDMonWNjcMXQgu+Xs+nVqQC0GHsUaUf03u/xNhIm8t2LRL57GawD/jgCR12Or8/pGKO11EREpGpUoItIjXx41ffkrdlBIMHPmS8MjblZW0I5haz6y//AWpp0b0fb35y43+Nt/gZCk+/DblkEgMnsSmD0LfiaZdVBWhERiSUq0EWk2hZNXM/cF1YDcNyD/WneM9XbQFFmIw6rHnqNcG4R/qQEOt98Pr5g5d8urbU4iz8i/MUTEC4FwD/gXPxDfoPxx9VlbBERiREq0EWkWgo3l/Du778FoMvolhz+f908ThR9m16dQuEPKwHoePXZxLeqvHfclhYQnv4ozorp7oakTIKjbsTXblBdRRURkRikAl1Eqsxay7uXfsOO7WUkpsdxxjOH4/PFVm914Q8r2PTyFACajzmC9KMq76131n9P6LMHoHg7AL7OxxAYfjUmoWmdZRURkdikAl1Equy7p1eydNImAE79+yBS2zbxOFF0hXILWflgRd95t7a0++3JPzvGRsqJzHqWyNzXAQuBBAJH/wlfz5N0I6iIiESFCnQRqZLs5YV8fO1cAPqd14G+v+jgbaAosxGHVeNfJ5xbiK9JPJ1v/iW+uD2/RTo5awhPvhe7fTkApkVPgsfdimna1ovIIiISo1Sgi8gBRcIOEy6eRXlxmNS2iZzyZOz1WG/+3zQK57iFd9ZVZxHf+qe+c2stzo/vEP7ynxApB+PDP+hX+AddiPHr26iIiESXfrKIyAHNGL+YdTOzATjjmcNJTI+t2UkK561k40uTAWh+2lDSj+63a5/dkUt46nicNTPdDSmtCI6+GV/rfpVdSkRE5KCpQBeR/dr4XQ7T7l4AwNArutFldCuPE0VXKK+IVQ++Bo4lsUsb2v3up77zyJpZhKc8BCW5APi6H0fg6Csw8clexRURkUZABbqI7FOoJMzbF8/CCVua90pl9H39vY4UVU4ozMr7XiGUU9F3fsv5+OKC2HAZ4a//jTN/gntgXBKBY6/C322Ut4FFRKRRUIEuIvs0+ZZ5bFtUgC9gOPP5IQQTY+dbhrWWNY+/TdH8VYDbd57QJhNn+3LCn96LzV0DgGndn+DomzApsfWXAxERqb9i56etiETVismbmfnEMgCG396XNoMqX6ynodr86lRyJs8BoM1Fx5M2rA/hua8Tmfk0OCHw+fEfdjH+gedjfH6P04qISGOiAl1EfqYkt5yJl8wGoP3QDIZd39PjRNGVM3UuG1/8FICM0YfS8tQ+hN67Hrv+OwBM03YERt+Cr2VsvW4REWkYVKCLyM98cMV3FGwoIS4pwNjnh+AP+LyOFDVFP65m9aNvApDcvzPtxrYh9PrvoKwAAF/vUwgc9QdMMNHLmCIi0oipQBeRPcz/31rmv7oWgBMeHkBG1xSPE0VP2cZsVtz9EjYUIb5dJp3GxhH5+DbAQkIqgeHX4u98tNcxRUSkkVOBLiK7FGzYwft/+BaA7ie3ZtBvO3ucKHrChSUsu/0FwvnFBNMT6DYmDzvnfQBMy94ET7gTk9zc45QiIiIq0EWkguNYJl4ym9K8EE0y4xnz1GEYY7yOFRVOKMyKe16ibP02gilhup20GtavBsDX80QCx16F8cfW4ksiItJwqUAXEQC++edyVkzeAsCYfw0mpVVs9GBba1n7+ASK5q0kMb2QTsPX4isqBOPDf+Tl+PufFTO/iIiISGxQgS4ibFtcwCc3/ADAgIuy6HVGO48TRc/m16aRPfl70tpvpe3ANZhIBOKTCR53O74Oh3kdT0RE5GdUoIs0cpGQw9sXzSRcGiGtYxNOeuxQryNFTc70eWx88WNa9VlDZpfNAJi0DgROvg9fWuz8EiIiIrFFBbpIIzf93h/Z+F0uxsDY54aQkBr0OlJUFC1cw7rHXyFryGKSW+QD4Os4lMDoWzDxyR6nExER2TcV6CKN2LqZ2/niwUUAHHl1T7KOaeFxougo25jNur/8i05DfyA+uRQA/8Dz8Q+5RKuCiohIvacCXaSRKi8O8/bFs3Ailpb9mjLy7r5eR4qKcGEJGx4eT8cB3+EPRrC+IMGR1+PvPtrraCIiIlWiAl2kkfr4+rnkLC/CH+fjzBeGEohv+CPLkfIQ25+4mTZd5mIM2Lg04k57AF/Lnl5HExERqTIV6CKN0NJJG/n23ysAGHVPP1r1T/M2UBQ4oVIK/3UF6enLAIgkdiTxF49gkjI8TiYiIlI9KtBFGpni7WW8c+k3AHQ8ujlHXNnd40QHzxZto+TlP5MY2ARAabA/qeP+gglo8SEREWl4fF4HEJG6Y63lvcu/pWhzKfEpAcY+NwSfv2F/G3A2/0jpK5cQiGzCWsgvO5KUSx5TcS4iIg2WRtBFGpEfXlrNognrATj58UNJz0ryONHBiSz+iNDUR/DZMJFyP9tyjqb9rbfg8zXsXzpERKRxU4Eu0kjkrSlm0hXfA9BrbDsOGZflbaCDYJ0Ika//ReSHNzFAWWECm9YMocuD1+KL07c1ERFp2PSTTKQRcBzLhF/PoqwwTHLLBE7752CMMV7HqhFbWkjok7ux678FoHBLGhuX9KP7I38kkNqw/yIgIiICKtBFGoWv/7qE1Z9vA+D0pw4jKTPe40Q14+SsIfzhLdj8DQBsW9aGrcuz6PbAr0lom+lxOhERkehQgS4S47bMz+OzW+cDMPjSLnQ/pY3HiWomsvprwp/eC6EdWOtn/fedyN+QSdb1Z5PSt5PX8URERKJGBbpIDAuXRXjroplEyh2adU3m+PGHeB2p2qy1ROa8RmTmfwCLY1JY9XlHSvKSaX3BKDJGDvQ6ooiISFSpQBeJYVPvXMCWefkYn+HM54cQnxz0OlK12HAZ4al/wVn2GeAuPrRsYgbhsjiajRxA61+N8jihiIhI9KlAF4lRa77YxpcPLwbg6Bt70X5ow+rRtkXbCH14G3bbEgCcVsNY8oKDU+aQ3DeLjlee1WBvdBUREdkfTRYsEoNKC0K8/etZWAttBqUz/LY+XkeqFmfzQsrf/L1bnBsf9LuYpf+LwylziG+bQZfbL9B0iiIiErNUoIvEoI+unkPe6mICCX7OfH4o/mDD+a8eWfwRoYlXwo4ciEvCN/Iulr+YTTh/B/6URLrefbGmUxQRkZimISiRGLPonfXMeX4VAMc90J/mvVI9TlQ1uy8+BGDS2hM4/h6WP/wppeu2YQJ+utw+TtMpiohIzFOBLhJDiraU8u5l7gI+nUe15PA/dPM4UdXsvfiQr8MQ/KNvYe0/P6Fw7goAOl51Fin9NJ2iiIjEPhXoIjHCWss7l37Dju1lJKQFOeOZw/H56v9NlHsvPuQfeB7+Ib9ly5szyP7YLdhb/2oUGaM0naKIiDQOKtBFYsR3z6xk6QcbATj174Np2q6Jx4kOLLJmprv4UHkx+IMEhl+Hv8dx5H4+nw3PfgRAsxEDaH2BplMUEZHGQwW6SAzIWVHEx9fMBaDvuR3od24HbwMdwN6LD5GUSfDEe/C17EnRorWsevh1AJL7ZNHxKk2nKCIijYsKdJEGLhJ2ePuimZQXh0ltm8gpTx7qdaT92nvxIdOiJ8GT7sEkZVK2OZcVd/0XWx4mvo2mUxQRkcZJP/lEGrgv/7KYdTOzATjj6cNp0ize40T7tvfiQ74exxM49hpMII5wUQnLb3+ecF5RxXSKFxFoqukURUSk8VGBLtKAbfw+h6l3LQBgyJ+60eW4Vh4n2jdn80JCH93mzm9ufPiPuAz/IedgjMGGI6y892VK1279aTrFds29jiwiIuIJFegiDVSoJMzbF8/CCVsye6Zw3P39vY60T5HFHxGe9ig4IYhLInj87fg6HA64/ehr//bOT9MpXnmmplMUEZFGTQW6SAM1+Zb5bFtYgC9gOOuFoQQT699/Z3fxoX8T+eENoGLxoZPvw5fWftcxW978nO0ffQNA61+NJGN0/e6hFxERqW317ye6iBzQys+2MPOJpQAMv60PbQY18zjRz9nSQkKf3o1d585lbjocTvC42zDxybuOyf1iPhue2Tmd4iG0vmC0J1lFRETqExXoIg1MSV45Ey6ZBUC7IRkMu6GXx4l+zl186FZs/noA/APOxT/0dxiff9cxxYvXsuov7nSKSb07ajpFERGRCirQRRqYSVd8T8H6EoJN/Jz5whD8AZ/Xkfawr8WHdle2OZfld1ZMp9i6GV3vGIcvLuhRYhERkfpFBbpIA7Lg9bXMe2UNACc8PICMrikeJ/qJdSJEvv0vkW9f5KfFh+7G13LPEf5IcelP0ykmJ9L17os1naKIiMhuVKCLNBAFG3bw/h++A6DbSa0Z/LsuHif6ic3fQGjyfdgti4A9Fx/a47hwhBX37Tad4m0XkNBe0ymKiIjsTgW6SANgrWXib7+hJLecJhlxnP6fw+pFv7a1FmfRJMIz/gbhUgD8h5yDf8hvMYG4nx279u/vUPj9cqBiOsVDOtd5ZhERkfpOBbpIAzD7n8tZ8elmAE7712GktEr0OBHYknzC0x7GWTXD3ZCUSXDUjfjaDar0+C1vfsH2D93pFFudP0LTKYqIiOyDCnSRem7b4gI+veEHAAZcmEXvse08TgTO2tmEpjzkrgoK+LocS+DYqzEJqZUenztjARuedadTTD+2P20uPK7S40REREQFuki9Fgk5vH3xTEIlEdI6NuGkxwZ6mseGywh//W+c+RPcDcEmBI6+Al+P4/fZclO8ZJ07naK1JPXuSNY1Z9eL9hwREZH6SgW6SD02/b6FbPw2F2Ng7LNDSGgad+CTaomzbRnhyfdhc91ZZEyrvgRH34xJbb3Pc8q25LL8zhexZSHiWjejy+0XaDpFERGRA1CBLlJPrZ+VzRcPLATgyKt7knVsC09yWCdC5IfXicx6Fpww+Pz4D7sY/8Dz91h4aG/udIovEM4twp+cQLe7LyKYlrzP40VERMTl6QonxphkY8wTxphNxpgSY8y3xpgxVTjvt8aYd40xayrOW1ZxHc3XJjGhvDjM2xfPxIlYWvZrysi7+3qSwxZuIfTuNUS+fgqcMCatPcGxfyMw6IL9Fuc2HGHlfa9QumYL+H0V0yl68wuGiIhIQ+P1CPoE4FDgemAVcDEwwRhzmrV20n7OuwuYCtwEbAB6A3cAY4wxA6y1ebUZWqS2fXLDXLKXFeGP83HmC0MJxO+7GK4tkWWfEZ7+mLsiKODrM4bAkb/HBPc/g4y1lrX/eJeC75cB0PHPZ5JySP2Zs11ERKS+86xAN8acDIwGzrTWTqjYNhXoDDwC7K9AH2it3brbx9ONMQuBacA44MlaCS1SB5Z9uIlv/rUCgJF396NV/7Q6fX5bVkT487/iLPvM3ZCYRmDE9fizjqjS+Vve+oLtk2YD0Oq8EWQeX/m0iyIiIlI5L0fQxwL5wDs7N1hrrTHmBeApY0xva+3Cyk7cqzjf6ZuKt97PQSdSQzuyy5j4O7e47Xh0c468qnudPr+zYS6hzx6Eoi0A+DoOJTDiOkyTZlU6P2fqXDY8s/t0iqNrLauIiEis8rJA7wsstNY6e22ft/v+alxvZMXbBfs6wBiTd4BrNK3G84lElbWW9y7/lqLNpcSnBBj73BB8/rq5TcRGQkRmP0dkzmuAhUA8gSP/D1+f06o8JeLWiV+y7l/vA5DUq4M7naLP09tcREREGiQvC/QMYGkl23N2218lxphmwBPAMuD1g48mUvfmvbyGhW+vB+Ckvx5KelZSnTyvk7OG8OR7sduXA2Cadycw+hZ86R2qdL61lo3Pf8Lm/00DIKl3R7redaGmUxQREakhr28StTXct4sxpgkwEWgGHGOtLdvnBa1NO8C18tAounggb20xH1zxPQC9zmjLgAuzav05rbU4CyYS/upfECkH48N/6C/xD74Q469acW0jEdY8PpHsT74FoOmQXnS+6Tx8Cd7N1y4iItLQeVmgZ1P5KPnOZtecSvbtwRiTCLwLDAROsNbOO8ApIvWO41gm/Ho2ZQUhklsmcNo/B9f6Spt2Rw6hKQ9h17r97qS0JDjqZnxt+lf5Gk5pOSsfeJX8WYsByDh+MB3/fAbGX/czzoiIiMQSLwv0H4GzjDG+vfrQ+1W83WcvOYAxJgH3BtMjgJOttV/VTkyR2jXz8aWsnu7e93z6U4eR1DyhVp8vsmoG4akPQ2k+AL4exxMY9idMfNUXEQoX7mD57S9QvGgtAK3OG06bi46v9V8sREREGgMvC/QJwCXAaew2kwtwIbBkXzO4ABhj4nHbWo4GTrPWTq/FnCK1ZsuCPCbf4v7hZ9DvOtP9lDa19lw2VEJ4xt9xFn3gbohPIXDsVfi7jqjWdcq35bHslucoXbsVjKH970+lxelH1kJiERGRxsnLAn0S7mJDzxhjMnAXKroIGAacvvMgY8w04Fhr7e5Dc28CJwB3A0XGmKG77dtmrV1Ry9lFDlq4LMLbF80iUu7QrEsyJ/xlQK091/+3d+fxVVXn/sc/zzknEwkkhFGmMIWAI4gi4IQo0tahamtvbR3Qtr/WDrZ6vZ3U1rZqB1upQydbLWr1tr2tqK2t4IBWZRBnAUkIhDAmQEgCIdMZ1u+PfTIQEkjISc5J8n2/XnkdsvZea68d4OTJyrOfFSldR+iFu3CV2wGwkdNIOvc7WEbHNt+tKS5lwy1/JLinEgv4GfvNT5F9VvvTYkREROTI4hagR2ueXwLcFf3IwiureJlz7h9H6H5h9PV70Y/mHsHbkVQkoS37wVpK3qvAfMZli04jJSP2VU9cJEz4rT8RfvNRcBHwJeGf+Xn8J30Ss46VQKxaV0zh9x4hXFWDr18KE267kgHTJsZ8ziIiIn1dXKu4OOf2AV+NfrR1zpxW2pToKj1a8Wu7ef3uDwE481tTGD1rcMyv4Sq3E3zhLlyply1m2WO98omDOx5UV6z8kE0//l9cXZBAVga5P1pAv9yRsZ6yiIiIEP8yiyJ9Tu2+IE8uWIVzcMzJAzn7tmNjOr5zjsj6fxN67QEI1gDgP/GT+Gd+AQt0vPzhnqVvUvzLxRCJkHxMNpPuvI6UEe3epkBEREQ6SAG6SDd77r/foWLzAQKpfi5bdBqB5NiVJXS1lYRe/gWRTa96DemDSZr7LXyjT+n4WM5R8pdX2LFoCQBpE0aQe8cCkgb2j9l8RURE5FAK0EW60fpntvPOH4sAmPfjExl6bOz2xYpsWU3wpZ9CdRkAvvFnEZhzE5ba8Wu4SIRtv3uWXU971Uv7T53AhNuuxJ/etSUgRURERAG6SLepKq3lmS+uBmD83GHM+EpuTMZ1oTrCKx4k/MGTXkNSGoEzb8CXN/+o6pJH6kNs/sX/Uf6KV/5x4FknMPbmT+FL1tuFiIhId9B3XJFu4JzjmS+u5sDuOlKzkrjk4Rn4fJ1/1jmyp5DQC3fi9m4GwIYdR9J538Uyj66eeri6jo0/+hP73ykEYMjFsxj9pQsxX8cqvoiIiMjRU4Au0g3efngT+f/cAcCFD0wnc1S/To3nXITwu/9HeNVDEAmC+fCfeg3+kz+L+Y4upz1Yvp/C2xZRXejNc8SC+Qz/r7O1O6iIiEg3U4Au0sX2bqziuZveBeD4/xrDCZ/O6dR4bv8ugi/9BLf9HQAsc6RXPnHYlKMes25HGRtueZi6nXvBZ+TccCmDP3Jqp+YpIiIiR0cBukgXioQjPLlgJfUHQvQfkcYF95/cqfHCG14i9J+FUFcFgO/YCwic/hUsKe2ox6wu3MGG2/5IqLwKSw4w/jtXkDUrtqUfRUREpP0UoIt0odfuXs/WFV5VlUsfmkG/7JSjGsfVVRF69T4iBc97DamZBM65Gf+4Mzo1v33vbmTjDx8jUl2HPyONiT+4mozjxnZqTBEREekcBegiXWTnO+Usu30NAKd9NZcJ84Yf1TiRHe8RfOHHUFUKgG/MaQTmfhPrl92p+e39z/tsvvuvuGCYpMGZ5N5xLWljh3VqTBEREek8BegiXSBYG+bv16wkEnIMntyf8+46scNjuHCQ8OpFhN/+X8CBP5nA6dfjO+7jnX5wc9czy9n6m3+Cc6SOGUrundeSPCSrU2OKiIhIbChAF+kCL97yPrvX7cMXMC5bNJPkfh37rxYp30LohTtwuzcAYENyCZx7C77sTj5g6hw7Hnmekj8vAyB9yhgm/vAaAv07V1VGREREYkcBukiMbXqxlBX3FgBw9q3HMfKU9qeiOOeIrH2a0PLfQqgOMPzTrsA/YwHmT+rUvFw4TPF9T1G25E0AMmfkMf67n8GXmtypcUVERCS2FKCLxFBNRT2LP7cKgFEzsjnz2+0vfeiq9xJ66WdEtnj9yRhG0nnfwTfipE7PK1IXZNOP/5fKlR8CMGjedHK+fikWOLqa6SIiItJ1FKCLxEjpmgqW3Pwu+7bVkNTPz2WPzMQfaN8OnOGi1wm9/HOoqQDAN2kegTNvwFIyOj2v0P4aNt7+KFVrNwMw/L/mMGLB+dqASEREJEEpQBfpBOccm14sZfk9+RQuLWlsn3/3VAbl9j9y/2ANodd/TWTdP72G5HQCZ9+EP3duTOZXv7uSDbf+kdpirwLMqC9dyLBLTo/J2CIiItI1FKCLHIVQfZi1f93K8oX5lLxX0dg+aFJ/zvzWFKZePfaIY0RK1xN64U5c5TYAbMRUks79NtY/NqUOa7bsYsMtDxPcXYkF/Iy9+XKy53Q+XUZERES6lgJ0kQ6oqajnrd9vZOX9G9i/o6axfexZQ5h9Yx65F4zA5zt86oiLhAm//QTh1YvARcAXwH/a5/FPvRyz9qXEHEnVumIKv/8I4f01+NKSmXDbVQw4eWJMxhYREZGupQBdpB3KNx9g5X0FvP3wJuqrQgD4/MaxnxzN7Bvz2l2pxVXuIPjiXbiStQDYwBwC827FNzh2wXPlqvVsvOsJXF2QQGY6E++4lvTckTEbX0RERLqWAnSRw9i+uozlC/NZ+7dtuIgDIDkjwPTPjWfmDZPIyklv1zjOOSL5Swi9ej8EqwHwn3AZ/ln/DwukxGy+Zc+/xeaFT0IkQvLwgeTeeR2pIwfHbHwRERHpegrQRVqIRBwFz+5g+S/yKX5td2P7gJFpnPa1SUz//HjSstpfO9zVVhJ6ZSGRja94Df2ySZr7LXxjZsRszs45Sv/2H7Y/9BwAaROOIfdH15KUfeQHVUVERCSxKEAXiQrWhHj30c2suLeAsoL9je3Dp2Yx+8Y8jrt8NIHkjtUNj2x9k+BLP4UDewDwjTuDwJybsbTMmM3bRSJs+/2/2LX4dQD6nzSeCd+7Cn96asyuISIiIt1HAbr0eVW7aln9m0Le+E0h1XvqGttzP3IMs2/MY9zcoR2uGe5C9YRX/p7w+3/zGgKpBM78Gr7JH41p/fFIMMTmX/yN8pffA2DgmScw9n8+hS9Z/7VFRER6Kn0Xlz5r9/p9rPhlPu89tplQXQQAf7KPEz+bw+xv5DH0uKNb5Y7s2eiVT9xbBIANO5ak876LZcb2Qc1wdR0b7/gT+98uBGDIRTMZ/aWLMH9sKsGIiIhIfChAlz7FOUfxf3bz+j35FDy7o7E9bWAyp35pIjO+MpH+w9OObux9JYQ//Bfhd/4MkSCYD/8pV+OffiXm61hqzJEEK6oovG0R1Ru2AzDi6nkMv+Ic7Q4qIiLSCyhAlz4hHIqw7u/bWH7Pena8Vd7Ynj0hg5lfn8S0a8aRnN7x/w4uWENk4yuE1y/B7Xi36cCAESSddwu+4cfGYPYHqyvZy4bvPkzdjjLwGTlfu5TBHz015tcRERGR+FCALr1a7b4gbz+8iZX3FVC5pbqxffSsQcy+aTKTLx6Br4MpIc5FcNvfIbx+CZFNr0Kotulgv2z8eR/BP/2zWHK/WN1Go+qNO9hw6yJC5fux5ADjv3MFWbNi/0OAiIiIxI8CdOmVKrdWs/L+At76wybq9gUBMIMpl45i9o15jJ7V8drgkYqtRPKXEs5fClW7mg74k/CNOwNf3nx8o0+JeTpLg/3vbaLwB48Sqa7Dn5HKxNuvIeP4sV1yLREREYkfBejSq+x8p5zlC/NZ89ctRELexkJJ/fxMu3Ycs27II3tCRofGc3VVRApfIrx+Ka507UHHbNhx+CfPxzfxHCylY+N2VPmrH1D0s7/ggmGSBg8g945rSRs7vEuvKSIiIvGhAF16POcchc+V8Po96yla1rSynTE8ldO+ksspX5xAv+z279bpImEiW1cTWb+EyObXIRxsOpgxFH/e+fjyzseXNTqWt9Gm3f9cyZZfPQPOkTJqCJPuuo7koVndcm0RERHpfgrQpccK1oZ5/4liVvwyn93r9jW2Dz1uALNvzOOEK3IIpLQ/3SRStonI+iWEN7wA1XubDgRS8U04C3/efGzkVMy6p4yhc46df3qBnY+/BED65NFM/ME1BDLTu+X6IiIiEh8K0KXHqS6rY/VvC3nj14VUlTY9oDn+3GHMvjGPifOHt7vcoKupILzhRSL5S3C7Nxx0zEZM9VJYJpyNJR1d6cWj5cJhtjzwDHv+/QYAA07NY/wtn8Gfmtyt8xAREZHupwBdeoyywv2suLeAdxcVEawJA+ALGCd8egyzvpHHMVMHtmscFw4SKV5JJH8JkeKVEAk3HRwwAv/k+fgnnY8NiE+Od6QuSNFP/0zF8nUADDrvZHK+cRkW6JqHT0VERCSxKECXhLdl+R6W37Oe9U9vx3nPfZKamcT0L0zgtK/mkjnqyOUMnXO43QWE85cQ2fAi1DalxJCcjm/iHC+FZfjxcd3sJ1RVw8bbH6VqzWYAhn3qbEZeO18bEImIiPQhCtAlIUXCEdY/vZ3Xf5HPtlVlje1ZOf2YecMkTr5uPCn9k444jjtQRrjgeS+FZe/mpgPmw0afgj9vPr5xp2OB9j9E2lXqy/ax4ZaHqd1cCsCoL17AsEvPiPOsREREpLspQJeEUlcV5J1FRay8r4DyTQca20eeks3sm/KYctko/IHDP6TpQnVEil4nnL8Et/VNcJHGYzYwB9/k+fgnzcPSO14LvavUbt3Nhlsepn5XBfh9jLv5crLPmRrvaYmIiEgcKECXhLB/Zw2rHtjAmw9upKa8vrE976IRzL4xj5wzhxw2zcM5hytZ46WwFL4M9U3BPSkD8E0610thGTIp4dJFDqzfwobvPUJ4XzW+1GQm3PZZBkyfFO9piYiISJwoQJe4Kl1TwfKF+XzwxBbCQW+lO5DqZ+rVY5n19UkMzhtw2P5ufwnh/OeJ5C/FVW5rOuDz48uZ6e3umTMT8x85HSYeKlfns+mOx4nUBQlkpjPxh9eQntc99dVFREQkMSlAl27nnGPTi6UsvyefwqUlje3pQ1I49fqJzLh+IulDUtvuH6whsvEVwvlLcdvfOeiYDc7Fl3c+/knnYWlZXXULMVH2wttsXvh3CEdIHjaQ3DuvJXXUkHhPS0REROJMAbp0m1B9mDV/2cLyhfmUvl/Z2D44rz+zvpHHSVfmkJTW+j9J5yK47e96KSwb/wOhpvrnpA3EP2kevsnz8Q0a39W3ERMlf/sP2//wbwDSxh9D7o8WkDTo8L8tEBERkb5BAbp0uZqKet58cCOrHtjA/h01je1jzxrC7BvzyL1gBD5f63nhkYptRPKXEM5/HqpKmw74k/CNPd0Lykefivl6Ro1wF4mw/aHnKP37qwBknDieid+/Cn96278xEBERkb5FAbp0mfLNB1h5Xz5vP1xEfVUIAJ/fOPaTo5l9Yx4jT8lutZ+rqyJSuMyrwlKy9qBjNuxYrzTixHOw1P5dfg+x5EJhNt/zN/a+9C4AWWccz7hvfgpfcmLmx4uIiEh8KECXmNv2RhnLF+az7u/bcBFvZ6HkjADTPzeemTdMIisn/ZA+LhImsvVNb3fPotcgHGw6mDHUS2HJOx/fwDHddRsxFa6pY9Mdj7PvrQ0ADLnwNEZffzHmP3zJSBEREel7FKBLTEQijoJ/7mD5PfkUv7a7sX3AyDRO+9okpn9+PGlZyYf2KyvyUlgKnofqvU0HAqn4xp+Jf/J8bOQ0zHpuIBusqKLwe49QXeBVmRlx1XkM/8zchCv3KCIiIolBAbp0SrAmxLuPbmbFvQWUFexvbB8+NYvZN+Zx3OWjCSQfnB/uaiqJbHiRcP5zuN0bDjpmI6binzwf3/izsOR+3XIPXamupJwNtzxM3fY94DPGfPUShnxsRrynJSIiIglMAboclapdtbzx60JW/7aQ6j11je25HzmG2TfmMW7u0INWiF04SKR4lZfCUrwSIqGmwQaMwJ93vreR0IDh3XkbXap6004Kb/0jwb37saQA47/zabJmHxfvaYmIiEiCU4AuHbJ7/T5W/DKf9x7bTKjO21jIn+zjxM/mMPsbeQw9LrPxXOccbneBl8Ky4SWobSqtSHI6vglzvBSW4cf3unSP/e9vYuMPHiN8oBZ/eioTbr+a/ieMi/e0REREpAdQgC5H5Jyj+D+7ef2efAqe3dHYnpadzKlfnMiMr0yk//C0pvMPlBEueIFI/nO4vZubBjIfNmq6l8Iy7gwskNKNd9F9yl9bQ9FP/4ILhkjK7k/undeRNq73/GZAREREupYCdGlTOBhh3d+3snxhPjveKm9sz56QwcyvT2LaNeNITvf+CblQPZGi14jkLyWydTW4SOP5NjAH3+T5+HPPwzJ6306ZLhymdnsZNZtLOPDhFnY9vRwijpSRg8m98zpShg+M9xRFRESkB1GALoeo3Rfk7Yc2sfL+Aiq3VDe2j541iNk3TWbyxSPw+X0454jsXOPt7lm4DOoPNA2SMgBf7lwvr3xoXq9IYXHOEdxTSU1RCTWbS6nZ7L3Wbt2FC4YPOrdf3ihyf7iAQOahJSVFREREDkcBujSq3FrNyvsLeOsPm6jb59UhN58x5ZKRzL4xj9GzBgPg9pcSyl9KJH8prnJb0wA+P74xp+HLm49v7EzMf2hZxZ4itL/GC8CLSqgpLqWmqITa4lLCB2rb7ONLSyYtZxgZx4/jmCvPxZ/ac+9fRERE4kcBurDznXKWL8xnzV+3EAl5Gwsl9fMz7dpxzLohj+wJGbhgDeH1S7zdPbe/C7jG/jZ4Ir68+fhzz8X69ax0jkhdkNotuxpXwxteg2X72u7k95E6eghpY4eTNnZY9HU4yUMzMV/PrdcuIiIiiUEBeh8ViTgKn9vJ8oX5FC3b1dieMTyV076SyylfnEDawCTcjvcIvriEyMZXINRs9ThtIP5J53mr5YMnxOEOOsaFI9TtLPOC8Gar4nU7yyDi2uyXPGygF4CPG0ZazjDSxg0nZeRgfEn6ryMiIiJdQ1FGHxOsDfP+E8Ws+GU+u9c1rRIPPW4As2/M44QrcvDXlhBe/zj1BUthf2lTZ18SvnGzvaB8zAzM52/lCvHlnCNYtq/ZangJtZtLqdmyC1cfarNfIDP9oNXw1LFeQO7v1zsrzYiIiEjiUoDeR1SX1bH6t4Ws+tUGDuxq2lho/LnDOP2mPMafnYHb9ArhZxcSLllzUF8bdiz+vPn4Jp6Dpfbv7qm3KVRVQ21x6SGr4uGqmjb7+FKSSM0Z1mxV3HtNGpg49yUiIiJ9mwL0Xq6scD8r7i3g3UVFBGu8SiO+gHHCp8cw6+sTGTZoE+H1fyD4yGsQrm/qmD4Ef948b7V84Jg4zd4TqQ9Ru3XXQaviNZtLCe6ubLuTz0fqqMHNVsW91+ThA5UnLiIiIglNAXov5Jxj6/I9LF+Yz/qnt+OiKdapmUlM/8IEZl6dTL+KVwi/93OC1WVNHQMp+Maf5ZVGHDm121NYXCRCXcleaopKqW14aLO4lNpteyASabNf8tAsb1V8XFMgnjpqCL5k/fMWERGRnkcRTC8SCUf48KntLL8nn22rmgLvrJx+nH7DCE6aXYRt/g3u5XyaV+22ESd5KSwTzsaS+3X5PJ1zhMqrDq6cUlRC7ZZdROqCbfbzZ6QdFIQ3vPrTU7t8ziIiIiLdRQF6L1BXFeSdRUWsvK+A8k1NmwWNOjWT875cy8hBb+C2rITVoabiiANG4M87H3/e+diAY7psbuHqukNKGNYWlxKqPNBmH0sOeBVTmq2Kp44dTlJ2/16x4ZGIiIjI4ShA78H276xh1QMbePPBjdSUe/njZo4ZVxgzLygm/cBy2F+J2x/tkNQP38Q5XgrLMSfENNiNBEPUbdtzyKp4/a6Ktjv5jJQRDXniTaviKccMwvzKExcREZG+SQF6D1S6poLlC/P54IkthINebnbm4FrO+8IeJuZ+gL+qGBozXAwbPd1LYRl3BpbUuXQQF4lQv6vi0O3ut+2GcNt54kmDB3gBeMOqeM4wUscMxZeS1Kn5iIiIiPQ2CtB7COccm14oZfnCfAqXlgDgD4SYevZ2Zl1UTFbKesxFoMo737LG4Js8H/+keVjGkKO6ZrCiqjEIr20IxotLidTUt9nHn57aVEd8bFMwHujf9bntIiIiIr2BAvQEF6oPs+YvW1i+MJ/S9ysBx8ixpZx6fhF5x2/AT7V3ogNS+uPLneulsAyd3O4UlnBNnbfdfVFJsxSVUkIVVW32saQAqWNabnc/jKTBmcoTFxEREekEBegJqqainjcf3MiqBzawf0cN/bOqmD2vgJPnFDKg/96mE82HL2cmvrzz8Y2dhfmT2xzThcLUbt9z8MY+m0uo37m3zT6YkXJM9sGB+LhhpIwYhPkTbydRERERkZ5OAXqCKd98gJX35fP2Q0W4YA15JxVx4mX55OTuwKyxBgs2aIKXwpJ7LtYv+6AxnHPU76poSktplifuguGWl2wUGNjf212zeTA+Zii+1LaDfhERERGJLQXoCWLbG2UsX5jPuie3MnrcduZ9vIAp0zaRnNKsLnjaQPyTzvV29xw8EYDQvgPUvLeJmuISaopKG1fFI9V1bV7L1y/Fe1izRXpKIDO9q29TRERERI4grgG6mWUAdwGXA1nAWuCHzrln2tF3AvAL4BzAB7wK3OycW9dlE+4C4VCER+e/QsXaDZxwagHX31pA1qD9TSf4kvCNm42Nm0tdZAz7t5RR82Q+NZtfoXZzKcG9+9sc2wJ+UkcNaVwVT23Y7n5olvLERURERBJUvFfQFwMnA98EioAFwGIzu8g596+2OpnZULyAfBdwDRACbgVeMbNpzrltXT3xWPH5Yd7cxxnyiY0HtYfTcqiJTKZyxxCqV5dTt+M5cK6NUSB5+MDoSnjz7e4HYwHliYuIiIj0JHEL0M3sY8B5wGXOucXRtmXAeLyV8TYDdOBmYCBwinNuR7TvCrwg/xbg+i6cekyZ+UjPTgYHoWAq5VsGU148iPqqNKAy+tEkkJnetN19jvfAZmrOMPxpKXGZv4iIiIjEVjxX0C/Fiz6fbmhwzjkzewR40MyOPUy6yqXA8w3BebRvmZn9A7iMHhSgA1RsGU3pNseB3ZmAl3riS00mNWdotGpK0wY/SVkZ8Z2siIiIiHSpeAboxwPrnHMtt598v/nxlp3MLA2YAPxfK2O+D3zGzIY653a10rfiCHPKPNKku0LqiTNxA8aS2Sw9JXlYFubTdvciIiIifU08A/RBQEEr7XubHW/NQLxl5taKdzfve0iAnqiOueKceE9BRERERBJEvB8Sbfupx8MfO6q+zrmsww0YXWGPyyq6iIiIiAh45QnjpYzWV8kbdt1pa3vLcrwA/Gj6ioiIiIgktHgG6GuBKWbWcg4nRF/XtNbJOVcDbMLLUW/pBGB3a/nnIiIiIiI9QTwD9MV4mxNd1KL9aiD/CBsOLQbmmdnwhgYzy46O9WSM5ykiIiIi0m3iGaD/C1gGPGRm15nZOWa2CDgD+J+Gk8zsZTNrmVP+c7wSjf8ys4+b2QXAs3gbFt3VLbMXEREREekCcQvQnXMOuAT4M15Q/W/gRLyNi/5xhL6lwJnAVuAx4C9ABXCWc25L181aRERERKRrmTvM9vF9jZlVZGZmZlZUVMR7KiIiIiLSi2VlZVFZWVnZWpVB7YQjIiIiIpJAFKCLiIiIiCQQBegiIiIiIglEAbqIiIiISAJRgC4iIiIikkAUoIuIiIiIJBAF6CIiIiIiCUQBuoiIiIhIAlGALiIiIiKSQLSTaDNmFgEsMzMz3lMRERERkV6ssrISwDnnDlkwV4DejJmF8H6rsK+VwxnR16rum1Hc9aR7TpS5duc8uvJasR47VuN1dpyGn74rOzkP6RqJ8v+4O/W0e06E+Xb3HPRe23F6r22fAUDEORdoeUABejuZ2csAzrk58Z1J9+lJ95woc+3OeXTltWI9dqzG6+w4ZlYR7Z/VmXlI10iU/8fdqafdcyLMt7vnoPfao+pfEe2f1Zl59GXKQRcRERERSSAK0EVEREREEogCdBERERGRBKIAXUREREQkgeghURHpM/TgkohI19N7bedpBV1EREREJIEoQBcRERERSSBKcRERERERSSBaQRcRERERSSAK0EVEREREEogCdBERERGRBKIAXUREREQkgShAFxERERFJIArQRUSizCzFzP5oZtvMrNLMlpnZcfGel4hIb2NmD5rZTjPbZ2YfmNmF8Z5TIlGALiLSJABsAmYC2cA/gKfiOSERkV7qHmCsc24AcB3wJzPLiu+UEofqoIuItMHMkoFaYIhzrize8xER6Y3MbBqwHDjFObc23vNJBFpBF5Eey8xGmdm9ZvaamVWZmTOzOW2cm2Fm90V/pVpjZm+a2cVHuMRsYJeCcxHpy7rqvdbMfm1mNcDbwIvAui67iR5GAbqI9GQTgSuAKrw398NZDHwWuBW4AO8bwWIz+1hrJ0d/1fo74LuxmqyISA/VJe+1zrkvAxnAPGCpU1pHI6W4iEiPZWY+51wk+udL8L4xnOOce7nFeR8DngUuc84tjrYZ8CowyDk3pcX5qcBzwFvOuf/u6vsQEUlkXfVe26LvP4H7nXNLuuQmehitoItIj9XwDaMdLgUqgaeb9XXAI8BkMzu2od3M/MCfge3AzbGbrYhIz9QV77Wt8AMTjnqSvYwCdBHpC44H1rXyTeb9Zscb/B5IAxbo160iIh3SrvfaaJ76ldHXgJldDpwDvNaNc01oCtBFpC8YBOxtpX1vs+OYWQ5wLXAWUB59GKrKzM7snmmKiPRo7XqvBRxeacVtQBnwbeAK59z7rfTtkwLxnoCISDc53Gq4A3DOFQPWPdMREemV2vNeewCY2z3T6Zm0gi4ifUEZTSs3zWVHX1tb8RERkY7Re22MKEAXkb5gLTDFzFq+550QfV3TzfMREemN9F4bIwrQRaQvWAxkARe1aL8ayHfOaXMMEZHO03ttjCgHXUR6NDP7ZPSPp0ZfzzazwcAB59y/o23/ApYBD5nZIKAIuAY4A/h4d85XRKQn0ntt99JGRSLSo5lZW29ixc65sc3OGwDcBXwSb4VnHfBD59xTXTxFEZEeT++13UsBuoiIiIhIAlEOuoiIiIhIAlGALiIiIiKSQBSgi4iIiIgkEAXoIiIiIiIJRAG6iIiIiEgCUYAuIiIiIpJAFKCLiIiIiCQQBegiIiIiIglEAbqIiEgLZvZrM9tqZvvMbLuZ/dLMkuM9LxHpGxSgi4iIHOoBYLJzbgAwFTgJ+FZcZyQifYYCdBERAcDMxpnZU2a228ycmS0ysznRPy+I9/xaMrPjzSxkZvNiPbZzbp1z7kDzJmBirMY3s0vMrN7McmM1poj0HgrQRUQAMxtvZg+a2XozqzazcjNbZ2aPmNk5zc6bama3m9nYOE63qywCzgZ+ClwF/C6uszmye4DXnXPPd8XgZvZtM9sP7AJOBO4/zLm/iqbCWHvGds49BXyA97UWETlIIN4TEBGJNzM7BXgFCAKPAmuBNGAScBGwH1gWPX0q8H3gZWBz986065hZCnAm8IBz7ufN2ufEa06HY2azgHnAJV11DefcT4CfmNkUvB9YStqYiwEfB552zrkOXOJe4BEzO845t7bTExaRXkMBuoiIF3D3A6Y5595tfsDMvgoMP9qBzcwPpDjnqjs1w643DDBgb7wn0k5fBsqAf3X1hZxzH5rZu3i/YTivlVNOBUYCT3Vw6CeB3wBfAr529DMUkd5GKS4iIpALlLUMzgGccxHn3A4AM7sd+GP00LJobrYzs0XR4wuin59nZreZ2UagFvhU9HiKmX3XzNaaWa2ZVZjZP8xsWvNrmll/M7vDzFaZ2R4zqzOzQjP7iZn1a3FuwzXPNbPvmVmxmdVE+86MnnO2mb1mZgfMbKeZ3dZijEVAcfTT7ze7rzltfcHMbHA0rWNrNJd6a/TzQc3OyYmOc3uLvkuj7d9o0b7KzNa1dc1m5wXwVs6fd84FWzn+0+j4k8zsvmjqyQEze97MRkfPucrM3oqmM+Wb2SVHuKwf799Jay4FKmn6LQtmNsDMbjGz982sMloNZp2ZPdBwjnOuCngVuPxI9ywifYtW0EVEYCOQZ2aXOeeePMx5TwLHAP8PuAv4sFn/5n4OJAG/B/YB+WaWBDwHzAYew6sSkgl8AXjdzM5yzr0Z7T8S+Dzwd+AJIISXG/5NYBowv5W5/QQviLwXSAb+G1hiZtcADwEPAo/j/bDwQzMrcs79Kdr3d8C7wEJgcfQ+id7flJYXMrNMYDneQ5MPA29H53U9MNfMZjjn9jvnis2sCDgXuD3aNxk4HYhE238ZbR8ATAd+28q9tTQdyADeaOP4NKAmeh/vAXcAJ+N9TX9lZlvw0nkej87j28DjZjbWObfbzDLwgubFeIH38cBtwJI2rncp8GzDDwvRdKFXgRy8H+jW4f2G5gS8tKnmVgDzzWyyc259O+5dRPoC55w+9KEPffTpD2AWUI9XqaMAL+i8HpjSyrkLoufNOcyxfKBfi2M3Ro/Nb9E+ANgCvNysLRlIamX8H0XHmNHKNd8Gkpu1XxxtDwGnthh7J7Cixdhjo+ff3qJ9TrR9QbO2O6NtX25x7lei7T9q1vaH6Nc2Pfr5WdFzHsP74SUQbb8o2n5ZO/6+ro2ee3Ebx3dHj3+mRfs70fZ/Nv/6Al+Ptp8b/TwdeB4v3acK2ATc3fLvNHrulGjfTzZruzzadn477uXK6LmfiPf/A33oQx+J86EUFxHp85xzK/BWZR/BW9W+Fvg1sM7MXjWz8R0c8jfu0JzzK4H1wFvR9JDBZjYYL2B+HjjDzNKi86l3TauxATMbGD33hehYp7Vxzfpmn78afV3pnFvd7F7r8VaeO1Pe71K8IPjBFu2/A/ZEjzd4Ce+3CWdEP5+LVxXlXqA/Xv42wDl4q9kvt+P6Q6Kvh+TLm9koYDDwjHPuiRaHy/F+YLnOHZwasy/6GgJwzh1wzs1zzmU75zKcc+Odc//Tyt8peKk2dXi/HWkwMPo6w8yO9H22LPo69AjniUgfogBdRARwzn3gnFvgnBuGt5p8DV6QewbwtHVsF8mCVtqmAJPxAtuWH9fhpacMbjjZzL5sZu/jBX97o+e9HD08kENtanE/5dE/FrVybjkwqJX29hoH5DvnQi2uGcL77UHzH2heir7Obfa6DG/Fv7xF+3vOufY8pNpQKaW1koYnR1//0sqx44FlzrldLdob0njy23Htli7Fy4Wvatb2N7zUmh8BO8zs92Z2YRvBesM9dKT6i4j0cspBFxFpwTlXDDxqZo/hBemnAzOA19o5RGsrrYZX9/qmw/TbDWBmNwG/AJYC9wE78NJERuJVEmkt0Au3MWZb7d3COVdiZh/i5ab3w1v9/5pzLmJmrwDnmtlv8eqM39POYXdHX7NbOdbwwO3K5o3Rh0OHtGyPOhnY4ZxrtYxiW6Kr9afgPUfQyDm318ym41V8+RheCcbPAyvMbE6L33Q03MNuRESiFKCLiLTBOefMbBVegD6yofkoh9uAFyC+5JyLHOHcq/BqrH+0+blm9pGjvHasbcJ7qDbQfBU9Wl1lEi1W8/FW0a/HyzNPBl6Mtr+I90DtR/F+gHmJ9lkTfW0tTWcaUOGcazmHhpX1t9vos7yd127uErx/D8+0POCcC+M9VLok+gPXIrw0p5OA1c1ObdiddA0iIlFKcRGRPs/M5kWDy5btacD50U8byv81pDK0tnp7OI/i1VNvdQXdzIY1+zSMF/hZs+MBvGojieApvB82Pt+i/QvR9sUt2l/C+37zfWCLc25js/YU4Dt4+d+v0j7v4OWNz2zl2DRaD8KnR1/fat5o3o6w2W30OZJLgdecc42r32Y2xOzg3USjwXrD3+n2FmPMBEqdc0eTXiMivZRW0EVEvPKCg8zsGbw0lGpgNPAZvBXhR51zH0TPXY33MOMtZjYQOAAUOedWHeEa9+LtfHm3mc3FC073AWPwyg3W4j0oCV4O84+Bf5vZk3iVXj6Dt9NpIvgZXqWSX5nZyXgB8zTgc3h53D9rcf4yvK/ZFLyVZACcc+vMrAQ4Fq+qzP72XNw5F45+XT5uZinOuToAM8vG+3r+tZVuJwN7nHNbW2mHDgbo0b/7s/BKXzb3c7wHfp8GCvF+MJkPXAjc7aI19aNjZOCVe3y4I9cWkd5PAbqIiLeq/XG8B0I/AWTh1b9+H/gpBweVW8zsOuBbeLtAJuFVfzlsgO6cC5rZBXg7YF4F/CB6aAdeVZVHmp1+N97q+efwAvsSvIceG2pqx5VzrtLMTse7h4vxqt6U4tUw/37LQNs5V27eTpwnc2gay0t4P3y0N72lwW/wSkxeiFcvHpryz9taQW8rvaWtPodzEd730KdatL+I97Dvp/B+m7AX7+/sEufc0y3O/QReffTfdfDaItLLmXN6cFxERHoeM3sOr776mXG49mJgnHNuaifGeAsods5dFrOJiUivoBx0ERHpqf4bmGVm5x/xzNhbAXz3aDub2SV4O4t+K1YTEpHeQyvoIiIiIiIJRCvoIiIiIiIJRAG6iIiIiEgCUYAuIiIiIpJAFKCLiIiIiCQQBegiIiIiIglEAbqIiIiISAJRgC4iIiIikkAUoIuIiIiIJJD/D2O2Jnh1+2a5AAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 864x576 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"cmap = plt.cm.plasma(np.linspace(0.25, 0.75, len(marginal_temps)))\n",
"labels = [f'T={t}'+r'${}^{\\circ}C$'for t in marginal_temps]\n",
"plt.gca().set_prop_cycle(plt.cycler('color', cmap))\n",
"plt.plot(flow_bins, marginalized_cdf, label=labels)\n",
"plt.legend()\n",
"plt.semilogx()\n",
"plt.xlabel('Streamflow ($m^3/s$)')\n",
"plt.ylabel('CDF')\n",
"plt.title(f'Conditional CDFs at {site_to_test}')"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "a4f53318-f8e8-4541-b26e-83ec51cbb026",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "bmorph",
"language": "python",
"name": "bmorph"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.9.2"
},
"widgets": {
"application/vnd.jupyter.widget-state+json": {
"state": {},
"version_major": 2,
"version_minor": 0
}
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment