From 096f2055c9c743a0c40d72e0dbefa0bc4c8f8ba9 Mon Sep 17 00:00:00 2001 From: rstiskalek Date: Mon, 8 Apr 2024 10:39:41 +0100 Subject: [PATCH] Add MDPL2 HMF --- notebooks/diagnostic/hmf.ipynb | 49 ++++++++++++++++++---------------- notebooks/diagnostic/hmf.py | 18 +++++++++++++ 2 files changed, 44 insertions(+), 23 deletions(-) diff --git a/notebooks/diagnostic/hmf.ipynb b/notebooks/diagnostic/hmf.ipynb index 0fed805..5a7b0bf 100644 --- a/notebooks/diagnostic/hmf.ipynb +++ b/notebooks/diagnostic/hmf.ipynb @@ -2,7 +2,7 @@ "cells": [ { "cell_type": "code", - "execution_count": 8, + "execution_count": 6, "metadata": {}, "outputs": [ { @@ -28,44 +28,45 @@ }, { "cell_type": "code", - "execution_count": 151, + "execution_count": 22, "metadata": {}, "outputs": [], "source": [ - "bin_edges = np.arange(12, 15.2, 0.25)\n", + "bin_edges = np.arange(11.8, 15.2, 0.25)\n", "xbin = (bin_edges[1:] + bin_edges[:-1]) / 2\n", "bin_edges = 10**bin_edges\n", - "xbin = 10**xbin" - ] - }, - { - "cell_type": "code", - "execution_count": 152, - "metadata": {}, - "outputs": [], - "source": [ + "xbin = 10**xbin\n", + "\n", "paths = csiborgtools.read.Paths(**csiborgtools.paths_glamdring)" ] }, { "cell_type": "code", - "execution_count": 153, + "execution_count": 23, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ - "100%|██████████| 20/20 [00:08<00:00, 2.36it/s]\n", - "100%|██████████| 20/20 [00:08<00:00, 2.41it/s]\n", - "100%|██████████| 20/20 [00:07<00:00, 2.64it/s]\n" + " 0%| | 0/20 [00:00" ] @@ -125,12 +126,12 @@ }, { "cell_type": "code", - "execution_count": 164, + "execution_count": 24, "metadata": {}, "outputs": [ { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHVCAYAAAB8NLYkAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAABsIklEQVR4nO3dd3iUVd7G8e+UZNJ7SAihl0AIBAwBsQEKoiKCbdVdFdC1sui7WFZ9d7Gtuqsu6/qaVVcX0bUBKghWFEF6J/TeW0JCek9m5v0jkjUCIZmZzEwm9+e6uC7zPGfO84sauDnnOecY7Ha7HRERERFp8YyeLkBEREREXEPBTkRERMRHKNiJiIiI+AgFOxEREREfoWAnIiIi4iMU7ERERER8hIKdiIiIiI9QsBMRERHxEQp2IiIiIj5CwU5ERETERyjYiYiIiPgIs6cL8KTDhw9z2223ceLECcxmM3/605+48cYbG/15m83GsWPHCA0NxWAwNGOlIiIi0lrZ7XaKi4tJSEjAaGx4TM5gt9vtbqrL6xw/fpzs7Gz69etHVlYWaWlp7Nq1i+Dg4EZ9/siRI7Rv376ZqxQRERGpHZBKTExssE2rHrFr27Ytbdu2BSA+Pp6YmBjy8vIaHexCQ0OB2n/RYWFhzVaniIiItF5FRUW0b9++Lnc0xKuD3eLFi3nppZdYt24dx48fZ/bs2YwdO7Zem4yMDF566SWysrJITU3l//7v/xg4cGCTn7Vu3TqsVmuTRuBOTb+GhYUp2ImIiEizasxrX14d7EpLS0lNTeWOO+7guuuuO+3+jBkzmDx5Mm+88QaDBg3ilVdeYeTIkezcuZM2bdoA0K9fP2pqak777Pz580lISAAgLy+P22+/nbfeeqvBeiorK6msrKz7uqioyJlvT0RERMSlWsw7dgaD4bQRu0GDBpGens5rr70G1C5maN++PZMmTeKxxx5rVL+VlZWMGDGCu+66i9tuu63Btk899RRPP/30adcLCws1YiciIiLNoqioiPDw8EbljRa73UlVVRXr1q1j+PDhddeMRiPDhw9nxYoVjerDbrczfvx4Lr300nOGOoDHH3+cwsLCul+HDx92uH4RERERV/PqqdiG5ObmYrVaiYuLq3c9Li6OHTt2NKqPZcuWMWPGDPr27cucOXMA+M9//kOfPn3O2N5isWCxWJyqW0RExJXsdjs1NTVYrVZPlyIOMplMmM1ml2yd1mKDnStcdNFF2Gw2T5chIiLikKqqKo4fP05ZWZmnSxEnBQUF0bZtW/z9/Z3qp8UGu5iYGEwmE9nZ2fWuZ2dnEx8f76GqRERE3MNms7F//35MJhMJCQn4+/trs/wWyG63U1VVRU5ODvv376d79+7n3IS4IS022Pn7+5OWlsaCBQvqFlTYbDYWLFjA7373u2Z9dkZGBhkZGRr2FhERj6mqqqpbNBgUFOTpcsQJgYGB+Pn5cfDgQaqqqggICHC4L68OdiUlJezZs6fu6/3795OZmUlUVBQdOnRg8uTJjBs3jgEDBjBw4EBeeeUVSktLmTBhQrPWNXHiRCZOnFi3SkVERMRTnBndEe/hqv+OXh3s1q5dy7Bhw+q+njx5MgDjxo1j+vTp3HTTTeTk5DBlyhSysrLo168f33zzzWkLKkRERERaA68OdkOHDuVc2+z97ne/a/apVxEREZGWwKuDnYiIiDTdzf9q3H6urvLx3YPd+jxvNH78eAoKCuq2T/MUTcw7ICMjg+TkZNLT0z1dioiISIuUlZXFpEmT6NKlCxaLhfbt2zN69GgWLFhQ16ZTp04YDAYMBkPd6t8777yT/Pz8ujaLFi1izJgxtG3bluDgYPr168cHH3zQ4LMPHDhQ1+fRo0fr3Tt+/HjdnnIHDhxo9Pfzj3/8g+nTpze6fXNRsHPAxIkT2bZtG2vWrPF0KSIiIi3OgQMHSEtL44cffuCll15i8+bNfPPNNwwbNoyJEyfWa/vMM89w/PhxDh06xAcffMDixYt54IEH6u4vX76cvn378umnn7Jp0yYmTJjA7bffzhdffHHOOtq1a8d7771X79q7775Lu3btmvw9hYeHExER0eTPuZqCnYiIiLjV/fffj8FgYPXq1Vx//fX06NGD3r17M3nyZFauXFmvbWhoKPHx8bRr145hw4Yxbtw41q9fX3f/iSee4Nlnn+WCCy6ga9euPPjgg1xxxRV89tln56xj3LhxvPPOO/WuvfPOO4wbN67eNavVyp133knnzp0JDAwkKSmJf/zjH/XajB8/vt559kOHDuWBBx7g0UcfJSoqivj4eJ566qlG/htynIKdK1hrPF2BiIhIi5CXl8c333zDxIkTCQ4OPu1+Q6NeR48eZd68eQwaNKjBZxQWFhIVFXXOWq655hry8/NZunQpAEuXLiU/P5/Ro0fXa2ez2UhMTGTWrFls27aNKVOm8MQTTzBz5swG+3/33XcJDg5m1apVvPjiizzzzDN8991356zLGQp2rjDnPtg0EyqLPV2JiIiIV9uzZw92u52ePXs2qv0f/vAHQkJCCAwMJDExEYPBwNSpU8/afubMmaxZs6ZRe9r6+flx6623Mm3aNACmTZvGrbfeip+f32ntnn76aQYMGEDnzp35zW9+w4QJE84Z7Pr27cuTTz5J9+7duf322xkwYEC9dwibg4KdA365eCI3J5uatdNh1nhY9SYUZ3m0PhEREW91rm3MfumRRx4hMzOTTZs21YWiUaNGnfH0p4ULFzJhwgTeeustevfu3aj+77jjDmbNmkVWVhazZs3ijjvuOGO7jIwM0tLSiI2NJSQkhH/9618cOnSowb779u1b7+u2bdty4sSJRtXlKAU7B/xy8UTNiD+zydyHoyeLKN80Bz67Cxa+ADm7PFuoiIiIl+nevTsGg4EdO3Y0qn1MTAzdunWje/fuXHrppbzyyissX76chQsX1mv3448/Mnr0aP7+979z++23N7qePn360LNnT2655RZ69epFSkrKaW0+/vhjHn74Ye68807mz59PZmYmEyZMoKqqqsG+fznyZzAYsNlsja7NEQp2LhDftS/n3f5Xwm6dzs6oS9lfYKVwxyLsX/wevnoUDq6AJv4NRURExBdFRUUxcuRIMjIyKC0tPe1+QUFBg583mUwAlJeX111btGgRo0aN4q9//St33313k2u64447WLRo0VlH65YtW8YFF1zA/fffT//+/enWrRt79+5t8nPcQcHOhUKj29Hv+kdof89Mjnb/NTtLg8jZt4GaBc/C7Htgx5dQU+npMkVERDwqIyMDq9XKwIED+fTTT9m9ezfbt2/n1VdfZfDg+psdFxcXk5WVxfHjx1m9ejWPPPIIsbGxXHDBBUDt9OuoUaN44IEHuP7668nKyiIrK4u8vLxG13PXXXeRk5PDb3/72zPe7969O2vXruXbb79l165d/OlPf/LaLc908kQzMAeEkDx8HFx6Gwc3fMfxtbOIObSXqJOvEbDhA+h5FfS8GgIjPF2qiIj4IG8/CaJLly6sX7+e5557joceeojjx48TGxtLWloar7/+er22U6ZMYcqUKQDExsaSnp7O/PnziY6OBmpXnpaVlfHCCy/wwgsv1H1uyJAhLFq0qFH1mM1mYmJiznr/nnvuYcOGDdx0000YDAZuueUW7r//fr7++usmfufNz2Bv6luMUqeoqIjw8HAKCwsJCwtrsO3JAxs5suQjInLWEBlkJjQoCEPXYdB7LER0cE/BIiLiMyoqKti/fz+dO3cmICDA0+WIkxr679mUvKEROzeJ7pRKdKdUyvOOsW/JR/jv+57IonlE7PoWc/t06H0dtO177o5EREREzkLBzgEZGRl17wc0VWBUAr3HPIS96l72Lf+MnM1zCN+yhKj9qwiI7wG9r4VOF4PR1AyVi4iIiC/TVKwTmjI0elY2G1mbF3By1QxCivcSFexPSHRbDL2ugR5XgH+Qa4sWERGfoKlY36KpWF9hNBKfOoL41BEUH9rM4aUfEHBwNeHZrxOR+RGmpCsgeQwEn/2lThERERFQsPMqoR36kPzrv1BTcIz9Sz7i5O7vCDn5IdGbZmPpNgRSroPorp4uU0RERLyUgp0XMkck0H30Q1B1L0dWzeZA5mwC1n5F1PbvCel0HoaU66FdGhgMni5VREREvIiCnTfzDybx4lvhwlvI3/YDx1bOwLx5BRF71xIe3xVzn+ug6zAw+Z27LxEREfF5CnYtgdFEZMoIIlNGUHFkM4eXfkjBvpUEHfkrkdHvENBnDPQcBZZQT1cqIiIiHqRg18IEJPah+80vYC86zqFlH3F4+3z8F7xJ5OoPCU0ZiaH3tRCW4OkyRURExAMU7BzgzD52rmIIa0vHKyfDpXdzYt3nZG+YTc7SWUSu/5ywpIsx97ke4pI9Vp+IiHjQ9Kvd+7zxX7j3eU4wGAzMnj2bsWPHerqUZmH0dAEt0cSJE9m2bZt3HABsCaHNBb+h+30ziLv6j+QHdmD/+u/J+ngSZbN/DweWgbYqFBERLzJ+/HgMBgMGgwE/Pz86d+7Mo48+SkVFhadLa/E0YucrjCZCeg0npNdwqo9t4uiyjyjdsxK//ZmExXYgYsAN0G0E+GkTSxER8bwrrriCd955h+rqatatW8e4ceMwGAz89a9/9XRpLZpG7HyQX0JfOt34Al3v+RBz72vIO5nN3i/+Rs47t1C95h0oy/N0iSIi0spZLBbi4+Np3749Y8eOZfjw4Xz33XcAnDx5kltuuYV27doRFBREnz59+Oijj+p9fujQoTzwwAM8+uijREVFER8fz1NPPVWvze7du7nkkksICAggOTm5rv+f27x5M5deeimBgYFER0dz9913U1JSUnd//PjxjB07lueff564uDgiIiJ45plnqKmp4ZFHHiEqKorExETeeecd1/9LcoCCnS8La0vCyN/T5b5ZRF9yNyU1Rvb/MJ3j/76F0h9ehvyDnq5QRESELVu2sHz5cvz9/YHa47XS0tL48ssv2bJlC3fffTe33XYbq1evrve5d999l+DgYFatWsWLL77IM888UxfebDYb1113Hf7+/qxatYo33niDP/zhD/U+X1paysiRI4mMjGTNmjXMmjWL77//nt/97nf12v3www8cO3aMxYsXM3XqVJ588kmuvvpqIiMjWbVqFffeey/33HMPR44cacZ/S42js2Kd4JKzYt3JWkPFnkVkr5hB1YldBJhNBHZKJ+b8myGhv6erExGRJmjwrFgvXzwxfvx43n//fQICAqipqaGyshKj0cjMmTO5/vrrz/iZq6++mp49e/Lyyy8DtSN2VquVJUuW1LUZOHAgl156KX/5y1+YP38+o0aN4uDBgyQk1O4W8c0333DllVfWLZ546623+MMf/sDhw4cJDg4G4KuvvmL06NEcO3aMuLg4xo8fz6JFi9i3bx9GY+14WM+ePWnTpg2LFy8GwGq1Eh4ezttvv83NN9/ctH93P9FZsdJ0JjMBScPpmDQc27FNHF/xEcX7VlK4dxXBbToTm349ph4jwGzxdKUiIuLjhg0bxuuvv05paSl///vfMZvNdaHOarXy/PPPM3PmTI4ePUpVVRWVlZUEBQXV66Nv3771vm7bti0nTpwAYPv27bRv374u1AEMHjy4Xvvt27eTmppaF+oALrzwQmw2Gzt37iQuLg6A3r1714U6gLi4OFJSUuq+NplMREdH1z3bkxTsWiljQl/aXd8XCo+Su2YmpVu/Ze+8lwgNf4vIvlcR0GcshMZ5ukwREfFRwcHBdOvWDYBp06aRmprKv//9b+68805eeukl/vGPf/DKK6/Qp08fgoOD+Z//+R+qqqrq9eHnV//kJYPBgM1mc3mtZ3qOu57dVHrHrrULb0fM8N/T8d5ZJFx6L9X4cWjJhxyb9huKv3oKjm/0dIUiIuLjjEYjTzzxBH/84x8pLy9n2bJljBkzhltvvZXU1FS6dOnCrl27mtRnr169OHz4MMePH6+7tnLlytPabNy4kdLS0rpry5Ytw2g0kpSU5Nw35SEKdg7IyMggOTmZ9PR0T5fiOpZQQtJvocNdH9Pp2inYIjuTvWURBz/6H/I/vBt2fgM1lZ6uUkREfNSNN96IyWQiIyOD7t27891337F8+XK2b9/OPffcQ3Z2dpP6Gz58OD169GDcuHFs3LiRJUuW8L//+7/12vzmN78hICCAcePGsWXLFhYuXMikSZO47bbb6qZhWxpNxTpg4sSJTJw4se5lRp9iMuPf4zISe1yG/fgmsld8TOH+leTN/QvhEW8S2W80pl6jISTW05WKiMjZtKCTIE4xm8387ne/48UXX2TDhg3s27ePkSNHEhQUxN13383YsWMpLCxsdH9Go5HZs2dz5513MnDgQDp16sSrr77KFVdcUdcmKCiIb7/9lgcffJD09HSCgoK4/vrrmTp1anN8i26hVbFOaHGrYh1VeJT8tZ9QvPUbKsrLiAi2EJ50MZY+10J8yrk/LyIiLtfgqlhpcbQqVtwnvB2Rlz1I5EV3ULb5C06u+4yD678nZOsiwhN7Epw6FroMBZPfuXoSERGRZqRgJ41nCSVowC0E9b+Rmn2LObHyY7L2b8V8cBsRUW8R3ncU9LwagqM9XamIiEirpGAnTWcyY+5+KQndL4WszZxYOYOCfcvJ/eEdIld/THjSJZh6j4G4ZE9XKiIi0qoo2Ilz4vvQZmwfKDpGwbpPKN7yNSfXfEP41gVEtE/GP2UMdB4CJv2vJiIi0tz0p624RlgCEcMeIOKCCZRv/ZK8tZ9ycNdGgg9sJTL6LQJTroakqyAoytOVioj4FK2B9A2u+u+oYCeuZQkl8LybaZd6A9Z9izmxagbHju/EfOLfRK35iNCkoZA8BmJ7eLpSEZEW7dTJB2VlZQQGBnq4GnFWWVkZcPopF02l7U6c0Gq2O3FW1mZyV82gbO9yqq02ooL9CWufjCl5DHS6WNO0IiIOOn78OAUFBbRp04agoCAMBoOnS5ImstvtlJWVceLECSIiImjbtu1pbZqSNxTsnKBg10RFxyha/ylFm7+mvLyU8EA/IqPj8UseBUlXQmCEpysUEWlR7HY7WVlZFBQUeLoUcVJERATx8fFnDOcKdm6iYOegymIqt33FybWfUpafTZC/maiwIAK6D4VeYyCmm6crFBFpUaxWK9XV1Z4uQxzk5+eHyWQ6630Fu2aWkZFBRkYGVquVXbt2Kdg5ylqDbf8SclbNoOz4DowGA9Eh/oQk9oHka6DjhWA8+//oIiIirYGCnZtoxM6FsraQt3oGpXuXU1VjJSrYn/DoeIw9R0HSFRDgY2fyioiINJKCnZso2DWDomOUZH5G4aavKCv76T280BD8ug2FXqMhuqunKxQREXErBTs3UbBrRpXFVG3/mrw1n1CSf4JAfxPRwf4EJPatDXgdLgCj0dNVioiINDsFOzdRsHMDaw32A0vJXfUxpcd2YjBQ+x5eVFtI+mma1hLq6SpFRESajYKdmyjYuVnWFvLXzqJ0z1Iqq396Dy8kBGPXobWLLSI7ebpCERERl1OwcxMFOw8pOk5Z5qcUbv6K0tIywgLNRAX7Y048D1Kuh4T+nq5QRETEZRTs3ETBzsMqS6je/hV56z6jJC+LQD8T0SH+WNp0g97XQuch2i5FRERaPAU7N1Gw8xI2K/b9Szi58gOKj+/BZIToYAvBkXG1U7Q9rgT/IE9XKSIi4hAFOzdRsPMydjscWUP+qg8oPrSJGpud6GB/wkLDMSSNhOSxEBzj6SpFRESaRMHOTRTsvFjWForXfEjR3lWUV9UQGexPREggps4XQ+/rtB+eiIi0GAp2bqJg1wKc3Ev5uo8p3LGQovJqwgP9iAr2xy+xf23AS0zzdIUiIiINUrBzEwW7FqToGNUbZ1Gw6WvyS8oJtphrNzxu0+WnhRZDwWT2dJUiIiKnUbBzEwW7FqgsD+vmTynInEd+YRF+JgPRIRZCIttAr2sg6UrwD/Z0lSIiInUU7NxEwa4FqyzGvn0eBes+pSD/JDY7xIT8tNCix+W1Cy1CYj1dpYiIiIKduyjY+YDqCtj1DUXrZlGQe4yKahvRdQstLqrd8FgLLURExIMU7NxEwc6HWGtg30JK131EQdZBSiqqCQ/yJzrYH792qbUBr10aGAyerlRERFoZBTs3UbDzQXY7HFxG5fqPyT+yg4KyKkIsZmJCLATEdq5daNFlmBZaiIiI2yjYNbOMjAwyMjKwWq3s2rVLwc5XHV1PdeYMCvavJ6+0igCzkegQf0IiYv+70MIS6ukqRUTExynYuYlG7FqJEzuwbpxB0e7lnCypwGAw/HSiRQjGHiMheQyExnu6ShER8VEKdm6iYNfK5B/AvvkTirb/QF5xGdVWO1HB/kQGW2oXWvS+DmJ7eLpKERHxMQp2bqJg10oVZ8OWTynZ+jV5RSWUVlmJDPQjKsSCf0Kf2oDXfqAWWoiIiEso2LmJgl0rV54P2z6nfPNc8goKKSyvJizATHSIhcCYjrV74XW7DEx+nq5URERaMAU7N1GwEwCqSmHHl1Rtmk1eXg75pVUE+puIDvYnNCIGel4NPUdBgP4fERGRplOwcxMFO6mnphJ2f0fNpk/IzzlKXkkVJqOB6BB/wkNDMHYbXrtdSlhbT1cqIiItiIKdmyjYyRnZrLBvEbbNn1B4fC8nSyqx2vlpoYU/5o4X1G543KanpysVEZEWQMHOTRTspEF2OxxejX3TTIqPbCG3pJLKahsRQX5Eh/jjH98bUq6DDoO10EJERM5Kwc5NFOyk0Y5vgs0zKT2wlpMllZRU1hAa4Ed0sD9BMR3+u9DCbPF0pSIi4mUU7NxEwU6aLHcPbJ5Fxd4lnCypoLC8hkA/EzEh/oRERGPoORp6Xa0TLUREpI6CnZso2InDCo/A5k+o2r2AvOJy8kur8DMZiA6x1C60SLqydhQvJNbTlYqIiIcp2LmJgp04rTQXts7GuuNr8oqKySupBCAy2J/IkCD8ug2BlBsgsqOHCxUREU9RsHMTBTtxmYoi2D4X27Z5FBTkcbKkimqbnYhAP6JDLFg6DawNePEpnq5URETcTMHOTRTsxOWqy2HXN9i3zKYoL4vckioqqq2EWGpPtAhJTIE+N0D7QVpJKyLSSijYuYmCnTQbaw3sWwibP6H4xAFySyopq7IS4GciOsSfsLguGPtcD12Ggcns6WpFRKQZKdi5iYKdNDu7HQ6tgM2zKD22ndySSkoqavAzGYkK9iciJh5zyrXQ4wrwD/J0tSIi0gwU7NxEwU7c6vjG2q1SDq4jt6SSwvJqjAYDEUF+REVGYuk9GpLHQGCEpysVEREXUrBzEwU78YiTe2HzLKr2Lia3pJKCsmrsdjthAX5EhYcQ3Ovy2hMtwhI8XamIiLiAgp2bKNiJRxUeha2fUb3re/KKSskrq8JmsxPkbyI6JIDQpCEY+twIMd08XamIiDhBwc5NFOzEK5TlwbY5WLd/RV5hIXklldTY7Piba9/DC+8yAHPfG6HdeZ6uVEREHKBg5yYKduJVKktg51fYts6hIC+HkyVVVFltmIwGIoP8iUjsiaX/r6DjRWA0erpaERFpJAU7N1GwE69UUwV7vsO+5TOKcg7X7YVnMBgICzATGdeB4PNuhG4jwOzv6WpFROQcFOzcRMFOvJrNBgcWw+ZZFB/fQ25JFWVVNQAEW8xERMUScd510HMUWEI9XKyIiJyNgl0jFRQUMHz4cGpqaqipqeHBBx/krrvuavTnFeykxTi8pnYvvCOb6vbCA7CYTUSGhxHR72pMKddBcIyHCxURkV9SsGskq9VKZWUlQUFBlJaWkpKSwtq1a4mOjm7U5xXspMXJ3la7F96BVeQWV1BYXg2A2WggIiSIiN7Da9/Di+jg4UJFROSUpuSNVn0WkclkIiiodrf+yspK7HY7rTjnSmsQlwxxTxKQf4DEzZ/QZs8icovLKSirJreolJMr5xKx6RvCul9ASNrNte1FRKTF8OqlcYsXL2b06NEkJCRgMBiYM2fOaW0yMjLo1KkTAQEBDBo0iNWrVzfpGQUFBaSmppKYmMgjjzxCTIymoqQViOwElzyM/6+mkTDoRronRBMTYsFggPzSSg5mLuTQBxPJm/UAHFpVe7SZiIh4Pa8esSstLSU1NZU77riD66677rT7M2bMYPLkybzxxhsMGjSIV155hZEjR7Jz507atGkDQL9+/aipqTnts/PnzychIYGIiAg2btxIdnY21113HTfccANxcXFnrKeyspLKysq6r4uKilz0nYp4SEgbOP9e/FJvJm77PKK3f0F+fh55JZUUV9RQvHsDBQc2ERLXmejzf42p6zAwefVvGyIirVqLecfOYDAwe/Zsxo4dW3dt0KBBpKen89prrwFgs9lo3749kyZN4rHHHmvyM+6//34uvfRSbrjhhjPef+qpp3j66adPu6537MRnVJfDzq9r98LLPV63Fx6An8lISFQ80QN/hSX5KvAL9HCxIiKtQ1PesfPqqdiGVFVVsW7dOoYPH153zWg0Mnz4cFasWNGoPrKzsykuLgZqw9nixYtJSko6a/vHH3+cwsLCul+HDx927psQ8TZ+gZByHcYbphE14hG6dU8iMTKQAD8T1VYb+TnH2PfVPzj8r5soXjENygs8XbGIiPxMi51Tyc3NxWq1njZtGhcXx44dOxrVx8GDB7n77rvrFk1MmjSJPn36nLW9xWLBYrE4VbdIi2AyQ4/LMXQfQfjB5YRvnkXx0e11e+EVFeZTvOgdglfNILj3SGIG3QJhbT1dtYhIq9dig50rDBw4kMzMTE+XIeK9DAbodCF0upDQY5mEbp5F2cF15Py0F15JWRkla2ZTtOkLLF0vJuGi2zDEdPN01SIirVaLDXYxMTGYTCays7PrXc/OziY+Pt5DVYn4sIR+kNCPoNw9dNw8i4q9S+r2wiuvrKZ82w+U71pEQKd02l0yDkP82Ue/RUSkebTYd+z8/f1JS0tjwYIFdddsNhsLFixg8ODBzfrsjIwMkpOTSU9Pb9bniHilmG4w7HECbnyLxIFj6R4fQWSwPwaDgcoaG4V7VrH3vUlkffwA9sNrPF2tiEir4tWrYktKStizZw8A/fv3Z+rUqQwbNoyoqCg6dOjAjBkzGDduHG+++SYDBw7klVdeYebMmezYseOsW5a4kk6eEAHK8mDrbKq3f0VeQSF5pVXYfvptJcDPRFi7HsQOvg06XVQ7tSsiIk3iM0eKLVq0iGHDhp12fdy4cUyfPh2A1157jZdeeomsrCz69evHq6++yqBBg9xSn4KdyM9UFsOOr6jZOoe8kznklVZhtdX+9hLoZyKybSciB/0GumgvPBGRpvCZYOftFOxEzqCmEnZ9S83mT8nPOcbJksq6gBfkbyaqTQLhA26CHiPBrFXmIiLnomDnJgp2Ig2wWWHfIqybZpF/bC8nS6qosdVudhxsMRMTE0tI/xug59XgH+ThYkVEvJeCXTPLyMggIyMDq9XKrl27FOxEGmK3w6GVWDfOoODwVk6WVFH902kWIRYzMVGRBPcdA8ljIDDCs7WKiHghBTs30YidSBMd34ht40wK960ht7iy7riy0AA/YiNCCex9FaRcD8ExHi5URMR7KNi5iYKdiINyd2PbNJPCnYvJLa6gqqY24IUF+hEbFkxAz+HQ50YIS/BwoSIinqdg5yYKdiJOKjiMfcunFG6ZT25RGZU1VgDCA/2IDQvE0vUS6PsriOrs4UJFRDxHwc5NFOxEXKQkB/vWzyjc+AUnC4upqLZiMEB4oD+xoRb8Ow6EvjdBm16erlRExO0U7JqZFk+INJOKItg+l4L1s8nLO0l5tRWDwUBEkB+xIRb8EvrWjuC1O8/TlYqIuI2CnZtoxE6kmVSXw86vKVg7k/zcbMqqajAYDEQF+RETYsHcpkdtwOswWKdZiIjPU7BzEwU7kWZmrYY9CyhY8zEF2Qcprfwp4AX7ExPijzmyI/S5AboMBaPJ09WKiDSLZg92c+fObXJRI0aMIDAwsMmf82YKdiJuYrPBwaUUrv6QgiM7KamswWgwEB3iT3SwBVNYPKRcB90vB5Ofp6sVEXGpZg92RqOxSe0NBgO7d++mS5cuTX2UV1OwE/GAI+soXP0BhQcyKa6oxmQ0EB3sT1SIBVNQVO1Gxz1HgZ9v/UVSRFqvpuQNh0/izsrKok2bNo1qGxoa6uhjRETqS0wjPDGN8BPbKVz1AUV7VnCiuJKTpVXEhFQRtfYdjJs/gV6jIfkasOj3HxFpPRwKduPGjWvStOqtt96qES0Rca02vQgf/WfC8/ZTtOYjirb/QHZRBSdLqogJrSIy80OMW2dD0hXQ+zoIivJ0xSIizU6LJxyg7U5EvFBxFsXrZlCw8SuKy8oxGw3EhFiICPLDaPKDbsNrF1qExnu6UhGRJtGqWDfRO3YiXqgsj+INn1Kwfg7FJcWYjQZiQy1EBPphMBqh8yXQ51cQ2dHTlYqINIrbg111dTVZWVmUlZURGxtLVFTrmPJQsBPxYlWllG78nJNrZlJSmIefqTbghQf6YcAAHc6vDXixPTxdqYhIg9wS7IqLi3n//ff5+OOPWb16NVVVVdjtdgwGA4mJiVx++eXcfffdpKenO/RNtAQKdiItQE0lpVu/JmflR5TnHcfPZCQ21EJYoLk24LVNrd3suG2qpysVETmjZg92U6dO5bnnnqNr166MHj2agQMHkpCQQGBgIHl5eWzZsoUlS5YwZ84cBg0axP/93//RvXt3h78hb6VgJ9KC2KyU7/yBrGX/oSL3AP4mI21CLYQG/LTvXWxP6HMjtB+o0yxExKs0e7C75ZZb+OMf/0jv3r0bbFdRUcH06dPx9/fnjjvuaOpjvJ6CnUgLZLdTsW85x5a8S1XWjtqAF2YhxPJTwIvsCH1vhk4XKeCJiFfQ4gk3UbATadkqD2/g8I/vYT2yDn+zkdjQAEIsP+0CFd6+doq28xBo4qbsIiKu1JS84fTvVgsXLjzrvTfffNPZ7kVEmo2lfX+63fp3Ok94i5r2F3C0oJwDJ0spq6qBwsOw5G8w517Y/T3YrJ4uV0TknJwesbNYLDzwwAM8//zz+PnVTmXk5uYyYcIEli5dSn5+vksK9Sbax07EN1XnHWL/wnex7vmBYD8DbUItBPiZam+GxNW+g9dtOJgcPrRHRKTJ3DoVu3z5cm6//XZCQkL48MMP2b9/P3feeSdJSUm89957dOzou3tFaSpWxDdVnDzM7vlvYj6whNAAE23CLPibfgp4wTGQcgP0GAkmP88WKiKtgtvfsSspKeHee+/lk08+wWaz8eyzz/Loo49i8PEXjxXsRHxb4dHd7Jn/BoFZa4gM8iM21IL51Pt2gZGQcj0kXQlmi2cLFRGf5tZ37AB27drF2rVrSUxMxGw2s3PnTsrKylzRtYiIx4S3607ahL8RceNrHLQksTu7hBNFFVhtdijPhzVvwyd3wOZPoLrc0+WKiDgf7P7yl78wePBgRowYwZYtW1i9ejUbNmygb9++rFixwhU1ioh4VEK3vpx/7z8xXfVXdhk6sPtEMbklldjsdqgohHXTawPexhlQpb/UiojnOD0V27ZtW6ZNm8aVV15Zd626uponnniCV199lcrKSqeL9FaaihVpfex2O+tWLKB85XTiqg//9xzaU6+e+AdDr2sg+RqwhHq2WBHxCW59xy43N5eYmJgz3vvxxx8ZMmSIM917NQU7kdarssbKyoVfYMj8gLb2bNqEWgg7dQ4tgF8g9Lwaeo+FgHCP1ioiLZtb37GLiYlhyZIl3HrrrQwePJijR48C8J///AfTqVVkIiI+xmI2MWTEGFLvf5/dSfeyuSiYfTmlFFdW1zaoLofNs2qnaNf8u/adPBGRZuZ0sPv0008ZOXIkgYGBbNiwoW7qtbCwkOeff97pAkVEvFl4sD9XjbmFTr99jw2Jt7Il35/9uT9tcgxQUwlbZ8Mnd8KqN6Esz7MFi4hPczrY/fnPf+aNN97grbfeqtugGODCCy9k/fr1znbvlTIyMkhOTiY9Pd3TpYiIl0iMDuWGW35L6G/eYVn0dWzJM3Aor5SK6p9OrLBWwfZ5tSN4K/4JJTmeLVhEfJLT79gFBQWxbds2OnXqRGhoKBs3bqRLly7s27eP5ORkKioqXFWr19E7diJyJna7nWU7j7J70Yf0zV9AXEA1bUIt+Jt/9nqK0QRdL6s9jzY03nPFiojXc+s7dvHx8ezZs+e060uXLqVLly7Odi8i0uIYDAYu6pnIb377MMdHvMa35qFsy63heGE51VZbbSObFXbPh8/uhiVTofCoZ4sWEZ/gdLC76667ePDBB1m1ahUGg4Fjx47xwQcf8PDDD3Pfffe5okYRkRbJ32zk6gHdufGux9ly/t/42n4B23MqyT61yTGA3QZ7f4A598KPL0HBIc8WLSItmtNTsXa7neeff54XXnih7rQJi8XCww8/zLPPPuuSIr2VpmJFpCmOF5Yze+kmgnd+SnrVGmJDzEQH+2P85fGLHS+E1JsgSrMeIuKBs2IBqqqq2LNnDyUlJSQnJxMSEuKKbr2agp2IOGJHVhGzF6+nw+G5pNWsp02IP5FBfqefr91+EKTeAjHdPFOoiHgFjwS71kjBTkQcZbfbWbHvJN8sW0dKzjz6WbcQF+pff5PjU9ql1Qa8Nj09U6yIeFSzB7vJkyc3uu3UqVOb2n2LoWAnIs6qsdr4ZmsWy1avYWDBV/S176RNqIXQAL/TG7dNrQ148SnuL1REPKYpecPsyAM2bNhQ7+v169dTU1NDUlISALt27cJkMpGWluZI9yIirYbZZOTqvgkMTRrF7PWpvJm5mkuKv6Z3yQHiwiwE+f/st+njG2t/xfWuDXgJ/TxWt4h4J4eC3cKFC+v+eerUqYSGhvLuu+8SGRkJQH5+PhMmTODiiy92TZUiIj4uxGLmtsGdONE7ng9X92PJztVcWvANPc1HiQsNIMDvZ3vgZW+F+X+E2J6QejMkDvBc4SLiVZx+x65du3bMnz+f3r1717u+ZcsWLr/8co4dO+ZUgd5MU7Ei0lx2Zxfz/sqDWA+vYVjZtyT5nTh9k+NTYrpD35uhwyD3Fyoiza7Zp2J/+bCcnNOPxsnJyaG4uNjZ7kVEWqXucaE8PSaF1fsT+GhVP5bmrGZY3nd0s+QTG2LBz/SzbUhzd8MPz0JU59qA1/EC+OUKWxFpFZwOdtdeey0TJkzgb3/7GwMHDgRg1apVPPLII1x33XVOF+iNMjIyyMjIwGq1eroUEfFxAztHcV6HCL7fHs9/1vanc9Fqhp38ni4BJUSHWDAbfxbg8vbDohcgoiOcd7tG8ERaIaenYsvKynj44YeZNm0a1dXVAJjNZu68805eeuklgoODXVKoN9JUrIi4U1lVDbM3HGX+5qOklK1kWMVCOgZVnHmTY4A2yTBgArTp5f5iRcRlPLKPXWlpKXv37gWga9euPh3oTlGwExFPyCmu5OPVh1i1J4u0iuUMq/qRDkHVRASdJeB1GAxp4yA80f3FiojTtEGxmyjYiYgn7csp4f2Vh9h97ASDypcytHopiSE2ws+0ybHBCN1HQL/fQFCUZwoWEYco2LmJgp2IeIN1B/P4YNUh8vJOMrj8R4ZYV5IYaiDEcoZNjk3+0PtaSLke/IPcX6yINJmCnZso2ImIt7Da7CzYns0n645QU5bP0LL5XMJ62ob5YznTFikB4dD3Jki6CkxOr6MTkWakYOcmCnYi4m3Kq6zM3XiULzYdJ7zyOJeXf0m6eR+xoQH1V9CeEhJXu4K28yXaIkXESynYuYmCnYh4qxNFFfxn5UHWHMijS9Uurqz4iuTAPCLPtsAiuiukTdAxZSJeqFmDXX5+Pna7naioKHJycliyZAlJSUmnnTzRGijYiYi323SkgOnLD3Asv4z+lau5onoBXYKrCA04w/t3UBvs0ibUBj0R8QpNyRvGBu/+wttvv01aWhoDBgzg9ddf59prr2XBggXcfPPNvP32204VLSIirtc3MYIXr+/LbYM7syP0Al4OeYgPKgazN6+KiuozbLJ+LBPmPQiLX4aSE26vV0Sc06QRu759+7Jq1SrKy8vp0KED+/fvJzY2lsLCQoYMGUJmZmYzlup9NGInIi1JYVk1H605xI87cwix5nNZ2bdcbNpMXKgFs+kMf883mqHnqNpFFgH6PU7EU5ptxM5sNhMYGEhUVBTdunUjNjYWgPDwcAx66VZExKuFB/lx75CuPDs2hTbxicwJvYmX/e/lh4I25JRUYvvl3/NtNbDtc/j0t7BpFtRUeqZwEWm0JgU7k8lERUUFAD/++GPd9ZKSEtdWJSIizaZbmxD+PDaFe4Z0pSy0M9NC7+E1+69Ymx9IYXn16R+oLoP178Jnd8Gu+WCzub9oEWmUJk3FnhoC/OXo3IkTJzh48CDp6ekuL9CbaSpWRFq68iorn6w7zLdbs7FZqxlQsYKrrAvpEGIjyP8s+9uFt4e08dBhkFtrFWmttN2JmyjYiYivOJJfxnsrDrLpSAEBtjIuKf+BEYY1xIea8TvT+3cAbZJhwARo08u9xYq0Mm4Ndi+88AJxcXHccccd9a5PmzaNnJwc/vCHPzjTvVdTsBMRX7PmQB7/WXGQE8UVRFpPMrL8Ky7w20V0iAXT2d6l7jAY0sZBeKJ7ixVpJdwa7Dp16sSHH37IBRdcUO/6qlWruPnmm9m/f78z3Xs1BTsR8UVVNTa+2HSMOZnHqKqx0qF6H1dXfk3fgBOEB/lh4AwBz2CE7pdDv19DUJT7ixbxYW4NdgEBAWzfvp3OnTvXu75v3z6Sk5PrFlv4IgU7EfFluSWVvL/yICv3nQS7nZSqDVxT8z1dg8sJPtv7d2YLJI+FlOvBP8it9Yr4qmbb7uRM2rdvz7Jly067vmzZMhISEpzt3itlZGSQnJzc6haLiEjrEhNi4X+G9+BPVyfTITqYLZbz+GvQ75lWdjH7Cm1U1Zxhg+OaStg0o3YF7ba5YK1xf+EirZjTI3YvvvgiL774Ii+99BKXXnopAAsWLODRRx/loYce4vHHH3dJod5II3Yi0lrYbHbmb8vmk3WHKamsIdhWzGXl33GpOZPYYH9MxrO8fxcSB+fdDp0vAe13KuIQt07F2u12HnvsMV599VWqqqqA2unZP/zhD0yZMsWZrr2egp2ItDZFFdXMXHOYBdtPYMdObE0WV1d+Rbr/QSKD/M6+WX1019ozaBP6ubVeEV/gke1OSkpK2L59O4GBgXTv3h2LxeKKbr2agp2ItFb7c0uZvmw/O7OLAehStZOxNd+SHJhHiMXv7B9M6F+7RUpUFzdVKtLyeWwfu1NdtZbjxRTsRKS1W7o7lw9WHSS/rAqD3Ub/yjWMsS+kU3AVFrPp7B/sMrR2ijakjdtqFWmp3Lp4AuDf//43KSkpBAQEEBAQQEpKCm+//bYruhYRES92UfcY/n5TP65JbYfJZGZ9wCCes/wP75UO5GixlRrbWcYO9i2Cz+6G1W9BRZFbaxbxZU6P2E2ZMoWpU6cyadIkBg8eDMCKFSt47bXX+P3vf88zzzzjkkK9kUbsRET+K6uwgndXHGDDoXwAwqwFXFE1n0vMW4kK8sN4ttkcvyDocyMkX1O7XYqI1OPWqdjY2FheffVVbrnllnrXP/roIyZNmkRubq4z3Xs1BTsRkdOtP5TPf1Yc5HhhOQBtaw5zbfU39LccJTSggffvgqKg363QbTgYXTKhJOIT3DoVW11dzYABA067npaWRk2N9i8SEWltzusQyUs39OWWgR0I8DNx3Nyefwbexd+qb2BzURAV1WfY/w6gLA+WvwqfT4RDq9xbtIiPcHrEbtKkSfj5+TF16tR61x9++GHKy8vJyMhwqkBvphE7EZGG5ZdW8eHqQyzdnYsdO0a7lfSKFVxnWky7QCtmUwPjC3EpkP5biOnmvoJFvJBbp2InTZrEe++9R/v27Tn//POB2nNiDx06xO23346f33+H3X8Z/lo6BTsRkcbZlV3MO8v2sz+3FIAAWxmXVf3ASNM6YoJMZ3//DqDLEDhvPITEuqdYES/j1mA3bNiwRrUzGAz88MMPzjzK6yjYiYg0nt1u54cdJ/h4zWGKK6oBiLLmMLrmOy4w7yQ80P/sHzb5Qa9roO9NOoNWWh2P7WPX2ijYiYg0XUllDbPWHua7bdnYfvojqEP1Pm60fUtvywmC/M1n/3BAOKTeDElXgbGBffJEfIiCnZso2ImIOO7QyTLeXXGArccKay/Y7fSt3MCNpoV0CijDr6H378La1Z5g0eF89xQr4kFuCXZ33HFHo9pNmzbNke5bBAU7ERHnrdh7kvdXHeRkSSUAZnsVF1UtZaxxGW2CDA2/fxeXAul3Qkx3N1Ur4n5uCXZGo5GOHTvSv39/Gupi9uzZjnTfIijYiYi4RmWNlc8zjzFv4zGqrTag9v27X1XO4byAo4Q1tP8dQOchkDZeCyzEJ7kl2E2cOJGPPvqIjh07MmHCBG699VaioqIcKrilUrATEXGtE8UVvL/iIKsP5NVesNsZULGCsdbvSAw1EODXwHt1WmAhPspt79hVVlby2WefMW3aNJYvX86oUaO48847ufzyyzE0NHTuIxTsRESax6YjBbyz7EDd6RXh1nxGl35Kmt8B2oQGYDI28GeMJQz63aIFFuIzPLJ44uDBg0yfPp333nuPmpoatm7dSkhIiCu69loKdiIizafaamP2hqPMzTxGja12erZfxRqurPiSDiF2IoL8MNBAwAtLgLQJ0HGwmyoWaR5uPVKsriOjEYPBgN1ux2o9y3ExIiIijeRnMvKrAe356w196Rlf+4dZZkA6/xc2mUVlndmXU0pZVQNHVxYdg4XPwdd/gNzdbqpaxLOcCnaVlZV89NFHjBgxgh49erB582Zee+01Dh065POjdSIi4h7tIgJ5cnQy9wzpSojFjxJjGDPCxvNRwK/YmgdH8svqFlycUfZW+OL38ONLUJLjvsJFPMDhqdj777+fjz/+mPbt23PHHXfwm9/8hpiYGFfX59U0FSsi4l5FFdX8Z8VBluyuDWhBthKuKP2cPlWbiAm1EB3s3/D2KFpgIS2Q27Y76dChA/37929wocRnn33mSPctgoKdiIhnbDlayL+X7q9bXJFUuYWrSucQZSwhPiyA0HNtj6IFFtKCuCXYjR8/vlErX9955x1Hum8RFOxERDznl4srLLZyRpbOI7VyLSEWM/HhAVjM5whtWmAhLYCOFGuisrIyevXqxY033sjLL7/c6M8p2ImIeN7RgnLeXrKP7ceLAOhStZPRJZ8Sbi8kKtif2BBLw9ujAMT1hvTf6gQL8UrNuio2Pz+fvLzajSNzcnL47LPP2Lp1q2OVeonnnnuO88/XeYMiIi1R7eKK3nWLK/b5J/F6xGTWWM4nt6SKPSdKKCirwk4D4xj1FliccF/xIi7WpGD39ttvk5aWxoABA3j99de59tprWbBgATfffDNvv/12c9XYrHbv3s2OHTu48sorPV2KiIg4YVhSG6belMrF3WOpMgbwdci1vBd+DycMkRwtKGd/binl1Q1sjwKw/0eYfQ+snQZVpe4pXMSFmhTsXn31VbZu3cratWt55JFHmD17NhkZGSxdupTXXnvN5cUtXryY0aNHk5CQgMFgYM6cOae1ycjIoFOnTgQEBDBo0CBWr17dpGc8/PDDvPDCCy6qWEREPCkswI+Jw7rxx1HJtA0P5JBfF96M+D0rAy+mrMrGvpxSjhaU1214fEbWatjyGXx6F2ybCzbtzSotR5OCndlsJjAwkKioKLp160ZsbO1hy+Hh4c1yhFhpaSmpqalkZGSc8f6MGTOYPHkyTz75JOvXryc1NZWRI0dy4sR/h9H79etHSkrKab+OHTvG559/To8ePejRo0ej6qmsrKSoqKjeLxER8T4p7cJ58Ya+XHdeIpgC+C54NNPCJ5JrakNBWRW7T5RwsqQSW0OvmVcWwep/wZz74OAK9xUv4oQmLZ5IT09nyZIlBAQEUFhYSHh4OAAlJSVcfPHFbNiwofkKNRiYPXs2Y8eOrbs2aNAg0tPT60YLbTYb7du3Z9KkSTz22GPn7PPxxx/n/fffx2QyUVJSQnV1NQ899BBTpkw5Y/unnnqKp59++rTrWjwhIuK9jhWU89ZPiyuMditDyr7jgvJFGLFhMZuID7cQYjnH9ihQu8BiwJ0Q27jBABFXabZVsac6/OXo3IkTJzh48CDp6emOVdwIvwx2VVVVBAUF8cknn9QLe+PGjaOgoIDPP/+8Sf1Pnz6dLVu2NLgqtrKyksrKyrqvi4qKaN++vYKdiEgLsGjnCd5feYiSymriao5xTcks4muOAhAa4Ed8mAX/c22PAtB5CKSNg5A2zVyxSK2mBDtzUzo+NUL3S23atKFNG/f+D56bm4vVaiUuLq7e9bi4OHbs2NEsz7RYLFgslmbpW0REmtfQpDac1zGS91ceZPEueDv8d1xQ/iNDyr+nuKKaksoaooP9iQm1YGro9aL9P8LBZZB86gSLYPd9EyLn0KRg90tZWVnEx8e7qhaPGj9+vKdLEBGRZhYW4Mf9Q7txSfdY/r10P8sMl7LTvzejS2aRWHOI3JJKCsuriQuzEB7of/aObDW1Cyx2fw+pN0PPUTrBQrxCk/ex+7nLL7/cVXU0WUxMDCaTiezs7HrXs7OzfSZsiohI8zi1uOL68xIptLTlnfCJzA++mmqDH9VWG0fya7dHqag+x4pYLbAQL+NUsPPkoRX+/v6kpaWxYMGCums2m40FCxYweHDzHg2TkZFBcnJys75TKCIizcvPZOTGAe356/V9SU6IYFXgJbwZMZmDfl0BKKuqYW9OCcfOtT0KQNExWPgcfP0HyNnlhupFzsypI8X69u3Lpk2bXFlPPSUlJezZsweA/v37M3XqVIYNG0ZUVBQdOnRgxowZjBs3jjfffJOBAwfyyiuvMHPmTHbs2HHau3fNQUeKiYj4jrrFFRVVpFWu5LLSr7DYaxfMmYwGYkMtRAX5N257Ly2wEBdy21mxzR3sFi1axLBhw067Pm7cOKZPnw7Aa6+9xksvvURWVhb9+vXj1VdfZdCgQc1W088p2ImI+JaiiuqfFlfkEGYt4OrST+latbPufoCfifiwAIItjXhF3eQHqbdAyvV6/06c4jPBztsp2ImI+Katxwr595L9HCssp2/FWi4vnUegvbzufligH3FhAfibGvFGU1RnuOBBiOnWjBWLL2tK3nDqHTuTSX8DERER39M7IZy/3tCXG9Lasz14IK9HPswO/5S6+0Xl1ew5UcKJ4oqGT68AyNsPX/6+9vzZmsqG24o4yakRu9YqIyODjIwMrFYru3bt0oidiIgPO15YzttL9rP1WCHJlRu5ovRzgm0ldff9TEbiwwMIC2jE6RWhbeGCSdC2bzNWLL7GbVOxrZ2mYkVEWo8fd+Xw/sqD1JQVcEXp56RUZta7H2wxEx8WQIBfI2azeoyEAXdoc2NpFLcGuzVr1vDYY4+Rk5NDt27d6NevX92vDh06ONO111OwExFpXYorqnl/5SF+3HWCHlXbGFXyGSG2orr7BgNEBvnTJjQAk/Ecq2cDI2HwROhwfjNXLS2dW4Ndz5496dChA9dccw379+8nMzOTzMxM8vPziYyM5OTJk85079UU7EREWqdTiytO5uczomwe/SvW1LtvNhpoExZAZFADp1ec0ukiGHRPbdATOQO3Brvg4GA2bdpE165d610/ePAgmZmZjBkzxpnuvZqCnYhI61VjtTEn8xifZx4lsXwnV5d+QoQ1v16b0AAzCRFBmM81eucfAum/he7Dm7FiaanctioW4MILL+TIkSOnXe/YsaPPhjqdPCEiImaTkRvSEnnxhr4EdRrAGxEPsSbgAuz8N8QVV9SwL6eEsqqahjurKoFlr8D8P0JxVvMWLj7NoRG76667jr59+5Kamordbuef//wns2bNIjKydQ0ja8ROREROObW4IqJ4F9eUzCLKmlt3z2CANqEBRIf4Y+Aco3dmC/S/DZLH1H5QWr1mn4p95JFHyMzMZOPGjeTm1v6PGx0dzZgxYzj//PPp378/ffr0wd+/Ee8WtGAKdiIi8nOnFlcs3XmUoWXzOb98CUb+e85sSICZdhGBmI2NmDCL6QEXPgiRHZuxYmkJ3PqO3dGjR+sWTJz6tW/fPsxmM0lJST59MoWCnYiInMm2Y0W8vXQfNbn7GFXyGYk1B+vu+ZmMJEYGEuTfiGPJjCbo8yvoexOYGtFefJLH97ErKSmpG9GbOHGiq7v3Ggp2IiJyNtVWG28t2cfinSfoX7may0q/qjuWzGCA2NAAYhozNQsQ3r529K5Nz2auWryRx4Nda6FgJyIi5/LlpuN8sOogFmsJw0u/IrVyLQZq/+gNsZhpF9nIqVmAXtfAebeDX0AzVizeptlXxR46dKhJ7Y8ePerIY0RERFq8UX3b8ocremIMCGde6I28G34vJ0zxAJRU1rAvp5TSc62aPWX7XPj8fji6vhkrlpbMoWCXnp7OPffcw5o1a87aprCwkLfeeouUlBQ+/fRThwv0RtruREREmiK1fQTPjk2hbXggh/0686+IB/k++CqqDf5UW20cyC0lp6QSO42YRCs5Ad9NgSV/g8ri5i9eWhSHpmJPnjzJc889x7Rp0wgICCAtLY2EhAQCAgLIz89n27ZtbN26lfPOO48//elPXHXVVc1Ru8dpKlZERJqitLKGVxfsZuORAgDCrLXnziZVbQUcmJoNCIdB90Lni5upYvEGbnvHrry8nC+//JKlS5dy8OBBysvLiYmJoX///owcOZKUlBRHu24RFOxERKSpbDY7H6w6yJebj9dd6161jStKPyfCmo/ZZCQxIpBgSxNWwbYfCOdPhODoZqhYPM0twW7KlCmMGTOGtLQ0h4r0BQp2IiLiqEU7T/Dvpfupttbuc2e2V3Fx2QIGly/GhJXYUAuxoZbGrZoF8AuCAROgxxXa2NjHuOVIsSNHjnDllVeSmJjIfffdxzfffENVVZWj3YmIiLQqQ5Pa8Kerk4kIrN3Mv8bgz8LgK/lXxP9w0K8rOcWVHDxZRo3Ndo6eflJdBisy4JvHoVCLFlsrp6ZibTYby5YtY968eXz++eccP36cESNGMGbMGK6++mqioqJcWavX0YidiIg462RJJS/P38X+3JJ61/tUrGdE2ReEU0a7yEBCmjI1a/KDfr+B3tfWbnIsLZrH9rHbvn17Xchbt24dAwcO5JprruGWW26hXbt2rnqM11CwExERV6issfLGon2s2Jdb77rFVs6lZd8woGIlsaH+TZuaBYjqUruxcXRXF1cs7uQVGxTn5OQwd+5c5s6dy8UXX8zDDz/cHI/xKAU7ERFxpdkbjjBzzZHTtj1JqD7MqNLP6GrKpl1EIH6mJrxJZTBCynWQ+msw+/YZ7r7KK4KdL8vIyCAjIwOr1cquXbsU7ERExGXWHczjtR/2UF5trXfdYLcxoGI5l1d+R+dwY9OmZgHCEuCCByDet3es8EVeEewOHz7Mk08+ybRp05qje6+gETsREWkOh/PKeOnbnZworjjtXoitiMtLv2Co/3ZiQywYmroCtscVMOAO8A9yUbXS3Lwi2G3cuJHzzjsPq9V67sYtlIKdiIg0l+KKal75fjdbjxWe8X6Xql3cbPuKXsHFTZuaBQiKhvPvhw6DXFCpNLem5I0mjuP+19y5cxu8v2/fPke7FhERafVCA/x44qpeTF9+gO+2ZZ12f59/D/5q78JlFUu40X8FoX5NGKcpOwk/PAudLoZB90BghOsKF49yeMTOaDRiMBho6OMGg0EjdiIiIk76bls27y4/cNY97aKsJ5kU/D09rbubPjVrCYX030K3y1xQqTQHt2xQ3LZtWz777DNsNtsZf61fv97RrkVERORnRiTH8cRVvQgN8Dvj/TxTNE9X3MS7/jdT5R/etM4ri2Hp3+G7KVBywgXViic5HOzS0tJYt27dWe+fazRPREREGi85IYznrk2hQ9TZFz18U9qNyVX3caTtiKYfK3Z0Pcy5H7bNBf353WI5PBW7ZMkSSktLueKKK854v7S0lLVr1zJkyBCnCvRmmooVERF3q6i2krFwD2sO5J21jQEDv+5ew9WlszHk7mz6Q2J7woUPQEQHJyoVV/GKVbGtgYKdiIh4gt1uZ9baI3y24UiD7XrGhTK5037Ctn4AVSUNtj2N0QR9b4Y+N4LJ4bWW4gJueceuNcvIyCA5OZn09HRPlyIiIq2QwWDgV+ntefCyHvibz34W7I7sYiZntmVT+l+g66VNe4jNCpkfwFcPQcEhJysWd3FoxG7Tpk2kpKRgNDYuF27dupWkpCTMZt9K/BqxExERT9ufW8rL83dysqTyrG0MGLg6tS03dyjCtPKfUHi4aQ8x+cF5t0Py2Ka/uydOa/YRu/79+3Py5MlGtx88eDCHDinti4iIuFrnmGCev7YPPeJCz9rGjp15G4/x9GoDuZe9DGnjwdSEc2Ot1bDm3/DtE1CS43zR0mwcGrEzGo3cfffdBAU17jiSf/7zn2zbto0uXbo0uUBvphE7ERHxFjVWG28v3c+inQ1vWRJiMXPf0G6kRVfDqjfg8OqmPcgvEAbeA92HO1GtNEWzL54YOnRokzdA/PDDD2nbtm1TH+XVFOxERMTbfLX5OO+vPIitoQMEMHBV37bckt4e89E1tQGvtIkjcR3Oh8G/06kVbqBVsW6iYCciIt5o05ECXl2wm5LKmgbbdW8TygOXdSc2EMj8ELbNqV000VgB4XDBAzpztpkp2LmJgp2IiHir44XlvPTtTo4VlDfYrnZqtitpHaNqV7+uyIDsrU17WLfLaqdn/Rv3ipY0jYKdmyjYiYiINyurquHVBXvIPJx/zrZX9WnLrwd2wGwywp7vYe07UFHY+IcFx8LFkyG+jxMVy5ko2LmJgp2IiHg7u93OB6sO8cWmY+ds2zU2hAeHd6dNaEDtGbLrpsOub5v2wOQxP626PfO5ttJ02qBYREREgNrNjG89vyP3D+2Gn6nhP/b35pTw+Keba48rs4TCBZPgot/XnkLRWNs+h3kPwsm9TlYujnA42O3btw8N9omIiLQMl/SI5cnRyUQGNbx/XWlVDX+bv5Ppy/ZTY7XVvj834hnwD278wwoOwZeTYePHYLM5Wbk0hcPBrnv37uTk/Hdp9E033UR2drZLihIRERHX69YmlOeu7UOX2JBztv1maxZT5m6lsKwa2qbCVS9BSJvGP8xmhQ3vw1cPQ+FRJ6qWpnA42P1ytO6rr76itLTU6YJaAp0VKyIiLVVUsD9Pje7NhV1jztl2X04JT8/bSl5pFUR0gFF/g5juTXtg7i6YOwm2zwPN9DU7hxdPGI1GsrKyaNOmNr2HhoayceNGnztdoiFaPCEiIi3Z55lH+Xj1Yew0HAXahAbwx6t71S6qqKmExS/DoRVNf2BCP7jw9xAc7VjBrZRbFk8YDIbTTp9o6mkUIiIi4jlj+rXj4ZFJBPo1vDjiRHEFT8/bxvHCcjBbYNgTtatfm+pYJnx+P+xd6FjBck5OjdhdeeWVWCwWAObNm8ell15KcHD9lys/++wz56v0UhqxExERX3Akv4yXvt1JdlFFg+0ig/z531G9SIz8aSPi7V/A6jcdm2LtdBEMnli7+lYa5JZ97CZMmNCodu+8844j3bcICnYiIuIrSipr+Pt3u9h6rOFNicMC/PjfUb3oGP3TQM7h1fDjX2unaJsqMBIufBASBzhQceuhDYrdRMFORER8idVm570VB/h2a1aD7UIsZh67shfd2vy0uvbkXvj+KSg/9wkXZ9RjJKTfBX4Bjn3ex2mDYhEREWkyk9HAhAs7c9fFXTAbzx4RSipreO7LbezIKqq9EN0VRk2FiI6OPXjXtzD3d5C9zbHPSx2HR+zuuOOORrWbNm2aI923CBqxExERX7X9eBF//24XRRXVZ21jMZt4ZGQSKe3Cay9UlcGi52sXSTgq5XrofxuYzI734WPcMhVrNBrp2LEj/fv3b/AEitmzZzvSfYugYCciIr4sp7iSJ+duqd3H7iz8TEYmj+hB/w6RtRdsVliRAbvnO/7gyE5w8UMQ1dnxPnyIW4LdxIkT+eijj+jYsSMTJkzg1ltvJSoqyqGCWyoFOxER8XX7ckp4at42qmqsZ21jNhp5cHh30jv9LAdsmgnr33P8wUZT7chdyvXQyrdTc8s7dhkZGRw/fpxHH32UefPm0b59e371q1/x7bff6gxZERERH9ElNoT7hnTFwNnDVY3Nxivf72L5ntz/Xuz7K7jkETA6OKVqs8K66fD1H6C44cUc8l8uWxV78OBBpk+fznvvvUdNTQ1bt24lJOTcZ9G1ZBqxExGR1uLTdUeYte5wg20MGLh7SBeGJf3sTNnsrfDDn6Gy2PGHmwMg/beQdIXjfbRgHlkVazQaMRgM2O12rNazD9eKiIhIy3N9WiIXnON8WTt2/vXjPub/fLuUuN5w1csQ2tbxh9dUwIrXardUKctzvJ9WwKlgV1lZyUcffcSIESPo0aMHmzdv5rXXXuPQoUM+P1onIiLS2tw7pCtdYxv+892OnWnL9vPlpuP/vRjeDka9DG16OVfAkbXw+UQ4sNS5fnyYw1Ox999/Px9//DHt27fnjjvu4De/+Q0xMQ0neV+jqVgREWlt8kur+N85mxtcKXvKTentubZ/4n8vWKthyVQ4sMT5QroMgUH3gcX3B5Lctt1Jhw4d6N+/P4YGVqvorFgRERHf0piVsqdc278dN6V3+O8Fux3WvwubP3G+kKBouOh/IKG/8315sabkDYd3/7v99tsbDHQiIiLim7rEhnD/0K784/vd2Gl4fGj2hqNU1di4bXCn2gsGA6SNh5B4WPlPsNscL6TsJMz/E/S8GgZMALPF8b58hM6KdYJG7EREpDVrzErZU0Ykx3PHhZ3qDwodXQeL/gLV5c4XE9audlPj2B7O9+Vl3LIqtqioqFG/fFFGRgbJycmkp6d7uhQRERGPacxK2VO+25bFm4v31d/rtl0aXPli7ZSqs4qOwlcPwYb3wVrjfH8tlFPv2DU0FWu32zEYDD699YlG7EREpLWrqrHx9Lyt7M0paVT7C7vGcP+wbpiMP8sQZXm1W5nk7XNNUdFda0fvIjqcu20L4JbFEz/++GPdP9vtdq666irefvtt2rVrV6/dkCFDHOm+RVCwExERadpKWYCBnaJ44LLumE0/mzisLq+dlj26zjVFmfzgvNsheWyLP5LMLcHul0JDQ9m4cSNdunRxRXctgoKdiIhIrf25pTw5d2ujVsoC9GsfyeQRPfA3/yzc2Wyw6g3Y+ZXrCovvC8MeB0uo6/p0M4+cPCEiIiKtV+eYYO4f2vCZsj+XeTifF7/ZQeXPg6DRCIPvhwF3uK6wrE21R5rVNG40saVTsBMRERGXOL9LNDekJZ674U+2HCvkha92UF71i1G+lOtg6ONg8ndNYdlbYfFLtXvo+TiXBjvtayciItK6NWWlLMCOrCL+/OU2Sip/sZK104VwxQsQEO6awg6tgJWvu6YvL+bwO3bXXXddva/nzZvHpZdeSnBwcL3rOnlCRESkdWnqSlmATtHBPDGqF2EBfvVvFGfVrpgtPOKa4vrfBqk3uaYvN3HLO3bh4eH1ft16660kJCScdl1ERERaF3+zkYcvTyI6uPEnQRw4WcrTc7dRUPaLd+FC4+GqlyAuxTXFbfgP7P7ONX15IZ084QSN2ImIiJxdU1fKArQND+R/R/UiJuQXodBaA8tegX2LnC/MYITLpkDiAOf7cgOtihURERGPa+pKWYDjheU8PW8bJ4oq6t8wmeGShyH1ZucLs9tg4fOQs8v5vryMgp2IiIg0m6aulAXIKa7g6XnbOF54hjNk+98KF/0PGE3OFWat+undvaPO9eNlFOxERESkWTV1pSzAydJKnpq7lcN5Zaff7DYcRjwD/sGn32uKyiL4bkrtkWY+QsFOREREmt29Q7rSNTakSZ8pLK/m6Xnb2J9bevrNtqm1iypC2jhXWEk2fP8kVJ0hQLZACnYiIiLS7BxZKQtQUlnNn7/Yxp4TxaffjOgAo/4GMd2dKy5vPyx8rnaBRgunYCciIiJuERnsz8Mjk/A3N+39uNKqGp77cjvbjxedfjMwEq74C3Q437nijm+EpX9v8adTKNiJiIiI23SOCWZiE1fKApRXW3nh6x1sPlJ4+k2zBYb9LySPca64/T/Cmred68PDFOxERETErQY5sFIWoKrGyovf7mD9ofzTbxoMMPAuGHRP7T87atvnsKXlnpqlYCciIiJud31aIhc2caUsQLXVxtT5u1i17+SZG/QaDcP+WDuK56i101yzEbIHKNiJiIiIR9zjwEpZgBqbjX8s2M3S3blnbtBhEFz5Yu37d45a+nc4tsHxz3uIgp2IiIh4hKMrZQFsdjsZC/ewcMeJMzeI7gqjpkJER8eKs1lrT6c4udexz3uIgp2IiIh4jKMrZQHs2PnX4n18uzXrzA1CYmv3unN0O5Tq8to97orP0r8XUrATERERj3J0pSzUhrt3lu1n3sZjZ27gHwTDn67d884R5QW1p1NUnGE1rhdSsBMRERGPc3Sl7CkfrDrIp+uOnPlmQBhc/mcIiXOs86JjtefKVlc4XJ+7tPpg16lTJ/r27Uu/fv0YNmyYp8sRERFptRxdKXvKrHWH+Xj1oTPfDIqCkc85vqAidzcseqH23Tsv1uqDHcDy5cvJzMxk4cKFni5FRESkVXN0pewpczKP8u7yA2e+GRpfO3Ln72D/R9fB8lcdrs0dFOxERETEazizUvaUr7cc5+0l+7Cf6XiwyI4w4hkwBzjW+Z4FsO5dh2trbl4d7BYvXszo0aNJSEjAYDAwZ86c09pkZGTQqVMnAgICGDRoEKtXr27SMwwGA0OGDCE9PZ0PPvjARZWLiIiIo5xZKXvK99uzef3HvdhsZwh3sT3gsj+B0exY55tnwfYvHK6tOXl1sCstLSU1NZWMjIwz3p8xYwaTJ0/mySefZP369aSmpjJy5EhOnPjvnjb9+vUjJSXltF/HjtWunlm6dCnr1q1j7ty5PP/882zatOms9VRWVlJUVFTvl4iIiLieMytlT1m8K4f/+2EP1jOFu7apMPQxMDgYhVa9AQeWOVxbczHYzzhO6X0MBgOzZ89m7NixddcGDRpEeno6r732GgA2m4327dszadIkHnvssSY/45FHHqF3796MHz/+jPefeuopnn766dOuFxYWEhYW1uTniYiISMM+W3+EmWsPO9XHBV1jmHRpNwxnOkN2z4LaUyYcYfKDEc9CfIpT9Z1LUVER4eHhjcobXj1i15CqqirWrVvH8OHD664ZjUaGDx/OihUrGtVHaWkpxcXFAJSUlPDDDz/Qu3fvs7Z//PHHKSwsrPt1+LBz/6OJiIhIw647z7mVsgDL9+by+o97z/zOXbfLYNA9jnVsrYYfnoX8A07V50oOTi57Xm5uLlarlbi4+nvSxMXFsWPHjkb1kZ2dzbXXXguA1WrlrrvuIj09/aztLRYLFosThwqLiIhIk90zpCtZRRXszSlxuI/Fu3LwMxq565Iup9/sNRoqSyDTgXftq0rhuyfhqpdrT7rwsBYb7FyhS5cubNy40dNliIiISAP8zUYeGZnE/87ewsnSSof7WbAjGz+TgfEXdj79Zr9boKoEtn3e9I7LTtaeTnHVi2AJdbg+V2ixU7ExMTGYTCays7PrXc/OziY+Pt5DVYmIiEhziAhyfqUswDdbs/hg1cEz3xx4F3QbfuZ751J4GBY8CzVVjhfnAi022Pn7+5OWlsaCBQvqrtlsNhYsWMDgwYOb9dkZGRkkJyc3OG0rIiIirtU5JpjfDevm1EpZgHkbjzFzzVnek7/gAeh4oWMdn9gGi18CD65L9epgV1JSQmZmJpmZmQDs37+fzMxMDh2qPS5k8uTJvPXWW7z77rts376d++67j9LSUiZMmNCsdU2cOJFt27axZs2aZn2OiIiI1DewcxQ3DnD8TNlTPttwhNkbznC2rNEIlzwCCf0d6/jQClj5T+eKc4JXv2O3du3aeue3Tp48GYBx48Yxffp0brrpJnJycpgyZQpZWVn069ePb7755rQFFSIiIuI7rjsvkaP55Szbm+tUPzPWHMbPZOTqvgn1b5jMcOkfYf4f4cT2pne882sIiobUm52qzxEtZh87b9SUfWVERETEdapqbDw9b6tTK2VPGX9BZ65IOcP7+ZUl8M1jjm9ncsED0ONyp2qDVrKPnYiIiLRep1bKOnOm7CnvLj/Agu3Zp9+whMDlf4awhNPvNcaK/4PDTTvq1FkKdg7Q4gkRERHPc9VKWTt23l6yn8W7ck6/GRgBlz8HwQ5skmy3w6K/QM5Op+prCk3FOkFTsSIiIp63en8ef/9uF3acizRGg4HfXdqNC8500kXhEfj6D1BR2PSOLaFw1UsQ7tiiD03FioiISKvhqpWyNrudjIV7WHMg7/Sb4Ykw4hnwC2p6x5XFtRsYl52hXxdTsBMREZEWzxVnygJYbXb+8f1u1h/KP/1mdFe4bAqY/JvecckJ+P5JqCpzusaGKNiJiIiIT7hnSFe6xoY43U+Nzcbfv9vFpiMFp9+MT4FhT4DRgff68vbDD38Ga43TNZ6Ngp0DtHhCRETE+7hypWy11cbL83ex7VjR6TcTB8DFDznWcdYmWDq12U6n0OIJJ2jxhIiIiPfZn1vKk3O3UlVjdbqvAD8TT1zVix5xoaff3PkNrHjNsY6Tx9SeTdsIWjwhIiIirZarzpQFqKi28pevd5x5I+SkKyBtvGMdb/sctnzqVG1nomAnIiIiPmdg5yh+le78SlmAsqoaXvhqOwdPlp5+s88N0OdGxzpe+w7sXehccb+gYCciIiI+6dr+rlkpC1BSWcOfv9zO4bwzrGpNGwdJVznW8bJX4Oh6p2r7OQU7ERER8Vn3DOlKtzbOr5QFKK6o5rkvt3O8sPz0m+ffB52HNL1TmxUWvQAn9zpfIAp2IiIi4sP8zUYevtw1K2UBCsqrePaLbZwoqqh/w2CAiydDogM7ZlSX125gXJzldH0Kdg7QdiciIiIth6vOlD0lr7SKZ77YRm5JZf0bRhMMfRziUpreaUVhbbgrL3CqNm134gRtdyIiItJyrNh7klcX7Hb6TNlT4sICeGp0byKDf3ESRVUZfPsEnNzT9E5jusPIF8AvoO6StjsRERER+YXBXaO59rx2Lusvu6iCZ7/YRmFZdf0b/kG158qGt296p7m7a9+5szm2B5+CnYiIiLQaN6YlMqhztMv6O1ZYzp+/3EZxxS/CXUAYXP5nCGnT9E6ProNl/3CoHgU7ERERaTUMBgP3D+tK55hgl/V5OL+M57/aTmnlL86ADY6Gy5+DwIimd7r3B1j3bpM/pmAnIiIirYrFbOKhy5OICPQ/d+NG2p9byvNfbae86hdTqGFtYcSz4O/AliubZ8H2eU36iIKdiIiItDoxIRYeurwHfibXRaG9OSX85evtVFT/ItxFdYbhT4HZgS1XVr0JB1c0urmCnYiIiLRK3eNCueviLi7tc2d2MS99u5PKml+EuzY94dI/gtHc9E6Xv9ropgp2DtA+diIiIr7hkh6xjE5NcGmfW48VMnX+Lqqttvo3EvrDkEdrNzNuClvNudv8RPvYOUH72ImIiLR8drudl+fvZN3BfJf2m9Yxkt8P74H5l9O9e76Hpa80up+i8mrC75+vfexEREREzsVgMPC7Yd1pHxnk0n7XHczn1R/2YLX9Ygyt23AYeLdLn3WKgp2IiIi0eoH+Jh4ZmURogJ9L+129/yT/XLiH0yZIk6+Bfr926bNAwU5EREQEgDZhAUwe0QOz0bXxaNneXN74cd/p4a7fryF5jEufpWAnIiIi8pNebcO486LOLu/3x10n+PfS/affSP9t7dSsiyjYiYiIiPzMsJ5tuDKlrcv7/X57Nu8uP1D/osEAFzwAHQa75BkKdiIiIiK/cNv5HUlNjHB5v19vOc4Hqw7Wv2g0wpA/QEI/p/tXsBMRERH5BaPRwIPDu5MQHujyvudtPMbMtYfrXzSZYdgfIbanU30r2DlAGxSLiIj4viB/M49ckUSIxYHTIs7hs/VHmLPhaP2LfgEw/EmI7Ohwv9qg2AnaoFhERMT3bT5SyF++2X76fnQucNv5nRjV9xfv85Xnw1ePQvFxQBsUi4iIiLhMn8Rwbju/U7P0/Z+VB/h2a1b9i4GRMPI5CIpucn8KdiIiIiLncEVKPMN7xTVL39OXHeCHHdn1L4a0gcv/DJamzQgq2ImIiIg0woQLO9M7Idzl/dqx89bi/SzelVP/RkR7GPEMmBu/gEPBTkRERKQRTEYDvx/Rg7iwAJf3bcfOGz/uZcXek/VvxHSDoY81uh8FOxEREZFGCrGYeWRkEoF+Jpf3bbPbeW3hbtYcyKt/Iy650X0o2ImIiIg0QWJkEA9c1h0DBpf3bbXZ+cf3u9lwKN+hzyvYiYiIiDRR/w6R/HpQh2bpu8ZmY+p3u9h8pLDJn1WwExEREXHA6NQEhvRo0yx9V1ttvDR/J9uPFzXpcwp2IiIiIg666+LOJMWFNkvfVTVW/vrNDvacKG70ZxTsRERERBxkNhmZfHkSMSGWZum/otrK1O92Nbq9gp0DdFasiIiInBIe6MejI3tiMbt+pSxAeZW10W0V7BwwceJEtm3bxpo1azxdioiIiHiBDtFB/O7Sbs2yUrYpFOxEREREXCC9UxS/Sk/0aA0KdiIiIiIucm3/RC7sGuOx5yvYiYiIiLjQPUO60jU2xCPPVrATERERcSF/s5GHL08iKtjf7c9WsBMRERFxschgfx66PAn/ZlopezYKdiIiIiLNoGtsCPcO6eLWZyrYiYiIiDSTC7rGcF1/962UVbATERERaUY3DkhkYKcotzxLwU5ERESkGRkMBu4f1o1O0cHN/iwFOxEREZFmFuBn4uGRSYQH+jXrcxTsRERERNwgJsTCQ5cn4WdqvvilYCciIiLiJj3iQvntxc23UlbBTkRERMSNhvSIZXRqQrP0rWAnIiIi4ma/HtiB8zpEurxfBTsRERERNzMYDEy6tDuJkUEu7VfBTkRERMQDAv1NPDoyidAA162UVbBzQEZGBsnJyaSnp3u6FBEREWnB2oQF8PvhPTAbXRPJFOwcMHHiRLZt28aaNWs8XYqIiIi0cMkJYUy4sJNL+lKwExEREfGwy3rFcUXveKf7UbATERER8QK3D+5E38QIp/pQsBMRERHxAkajgQcv607b8EDH+3BhPSIiIiLihGCLmUevSCLY3+zQ5xXsRERERLxI2/BAHhzeHaPB0OTPKtiJiIiIeJm+iRHcPrhjkz+nYCciIiLiha5IactlPeOa9BkFOxEREREvNeHCTvSID210ewU7ERERES9lNhn53bBujW6vYCciIiLixZpylqyCnYiIiIiPULATERER8REKdiIiIiI+QsFORERExEco2ImIiIj4CAU7ERERER+hYCciIiLiIxTsRERERHyEgp2IiIiIj1CwExEREfERCnYiIiIiPkLBTkRERMRHKNiJiIiI+AgFOxEREREfoWAnIiIi4iPMni6gJbPb7QAUFRV5uBIRERHxVadyxqnc0RAFOycUFxcD0L59ew9XIiIiIr6uuLiY8PDwBtsY7I2Jf3JGNpuNY8eOERoaisFg8HQ5pKens2bNGk+X4ZY6XPkMZ/ty5PNN+Uxj256rXVFREe3bt+fw4cOEhYU1utaWylt+HqD5a9HPQ9Pb6efBc/Tz0PTP2O12iouLSUhIwGhs+C06jdg5wWg0kpiY6Oky6phMJq/4DcoddbjyGc725cjnm/KZxrZtbLuwsDCv+P+kuXnLzwM0fy36eXC8nX4e3E8/D4595lwjdado8YQPmThxoqdLANxThyuf4Wxfjny+KZ9pbFtv+e/vLbzp30dz16KfB9fU4cu86d+Hfh6a95maihVpJYqKiggPD6ewsNBr/uYu4in6eRBfpRE7kVbCYrHw5JNPYrFYPF2KiMfp50F8lUbsRERERHyERuxEREREfISCnYiIiIiPULATERER8REKdiIiIiI+QsFORERExEco2Im0Utdeey2RkZHccMMNddcKCgoYMGAA/fr1IyUlhbfeesuDFYq415l+Jk4pKyujY8eOPPzwwx6oTKTxFOxEWqkHH3yQ9957r9610NBQFi9eTGZmJqtWreL555/n5MmTHqpQxL3O9DNxynPPPcf555/v5opEmk7BTqSVGjp0KKGhofWumUwmgoKCAKisrMRut6OtLqW1ONPPBMDu3bvZsWMHV155pQeqEmkaBTuRFmjx4sWMHj2ahIQEDAYDc+bMOa1NRkYGnTp1IiAggEGDBrF69epG9V1QUEBqaiqJiYk88sgjxMTEuLh6Eddrzp+Jhx9+mBdeeMHFFYs0DwU7kRaotLSU1NRUMjIyznh/xowZTJ48mSeffJL169eTmprKyJEjOXHixDn7joiIYOPGjezfv58PP/yQ7OxsV5cv4nLN9TPx+eef06NHD3r06NEcZYu4nNnTBYhI01155ZUNTgtNnTqVu+66iwkTJgDwxhtv8OWXXzJt2jQee+yxRj0jLi6O1NRUlixZcsaXyUW8SXP9TKxcuZKPP/6YWbNmUVJSQnV1NWFhYUyZMsXl34OIK2jETsTHVFVVsW7dOoYPH153zWg0Mnz4cFasWNHgZ7OzsykuLgagsLCQxYsXk5SU1Kz1ijQ3Z34mXnjhBQ4fPsyBAwd4+eWXueuuuxTqxKtpxE7Ex+Tm5mK1WomLi6t3PS4ujh07dtR9PXz4cDZu3EhpaSmJiYnMmjULk8nE3XffXbdoYtKkSfTp08fd34KISznzMzF48GB3lyviFAU7kVbq+++/P+P1zMxM9xYi4iXO9jNxyvjx491TiIgTNBUr4mNiYmIwmUynLXrIzs4mPj7eQ1WJeI5+JqQ1UbAT8TH+/v6kpaWxYMGCums2m40FCxZoWklaJf1MSGuiqViRFqikpIQ9e/bUfb1//34yMzOJioqiQ4cOTJ48mXHjxjFgwAAGDhzIK6+8Qmlpad2KQBFfo58JkVoGu7aVF2lxFi1axLBhw067Pm7cOKZPnw7Aa6+9xksvvURWVhb9+vXj1VdfZdCgQW6uVMQ99DMhUkvBTkRERMRH6B07ERERER+hYCciIiLiIxTsRERERHyEgp2IiIiIj1CwExEREfERCnYiIiIiPkLBTkRERMRHKNiJiIiI+AgFOxEREREfoWAnIiIi4iMU7ERERER8hIKdiEgzGDp0KAaDAYPBQGZmpqfL4eGHH2bs2LFnvT9+/Pi6eufMmeO2ukTEtRTsRMQnDRkyBIPBwPPPP1/vut1uZ9CgQRgMBp555plmreGuu+7i+PHjpKSkuL22CRMm8Mc//rHu68zMTPr27XvW9v/4xz84fvy4U88UEc9TsBMRn2O329mwYQMdO3Zk8+bN9e69++67HDt2DIDzzjuvWesICgoiPj4es9ns1tqsVitffPEF11xzTd21jRs3kpqaetbPhIeHEx8f7/AzRcQ7KNiJiM/ZvXs3xcXFjBs3rl54Ki4u5vHHH2f8+PEApKWl+WRty5cvx8/Pj/T0dACOHDlCbm4uACNGjCAoKIikpCRWrVrl+DciIl5JwU5EfM66desICgrilltuYefOnVRVVQHw7LPPMmDAAGJjY4mPj6dt27Y+WdvcuXMZPXo0BoMBoO4dv4yMDJ544gk2btxIhw4deOyxx5z+fkTEu5jP3UREpGVZv349ffv2JSkpiYCAAHbs2EFgYCCvv/4669ev57nnnqub6jSbzfXegVuxYgWBgYFn7bup7Z2pzVGff/45f//73+u+zszMJCoqipkzZxITEwPANddcw5tvvunUc0TE+yjYiYjPWb9+Peeddx4Gg4G+ffuyefNmPvroI+677z66d+/OunXruPbaawGIiIho0qrVprZ3prZfmjNnDm+//TZVVVXcdNNN3Hnnnae12b59O8eOHeOyyy6ru5aZmcmYMWPqQh3A/v376datm8Pfh4h4J03FiojPORWeAPr168crr7zC2rVr+dOf/kRFRQU7duw456jYX//6V1JSUujTpw8ffPCBx2v74IMPmDlzJq+//jrvv/8+27Zt47nnnjut3dy5cxkxYgQBAQF11zIzMzn//PPrtcvMzKRfv34u+75ExDso2ImIT9m3bx8FBQV14ah///6sXbuWF154gdDQUDZu3EhNTU3d4oSCggL69etHv379+O1vfwvAmjVrmDlzJmvXruXHH39kypQpdatVz9S+uWr7uX/961+8++67tG/fnjZt2vC3v/2NRYsWUVxcXK/d559/zpgxY+q+Li4uZt++ffTv379eOwU7Ed+kqVgR8Snr1q3D39+/7j24cePGMXbsWKKjo4HaEbPY2Fjat28PnHlqddmyZVx//fUEBAQQEBDAZZddxpo1axgzZoxTU7FNre2UkydP0qFDB/z8/Hj77bdZunQp06dPZ9CgQezatasuCJ44cYK1a9cyd+7cus9u3LgRk8lEnz596q4dPHiQ/Px8BTsRH6RgJyI+Zf369aSkpODn5weAn59fvXfL1q9ff9rolbfXFhkZWbd58I033sioUaMA2LlzZ70QOG/ePAYOHFivz8zMzLqFGqds2LCBiIgIOnXq5NLvT0Q8z2C32+2eLkJExFNiYmLq9ng7Ze3atdx7770sW7aMsrIy0tPTWbJkCW3btj1j+zMZOnRo3Tt0rjBlyhQiIiKYPHkyADNnzuSrr75i+vTpdW2uueYaLrroIh599FGHn2MwGJg9e3aDx4+JiPfSO3YiIr8wYMAAbrzxRtLS0rjkkkt4+umnHdpX7p///CchISGnnTDhiClTpnDixAlSU1Pp168f3333HRkZGfXaXHTRRdxyyy0O9X/vvfcSEhLidJ0i4lkasRMRaQZHjx6lvLwcgA4dOuDv7+/hihp24sQJioqKAGjbti3BwcEerkhEHKFgJyIiIuIjNBUrIiIi4iMU7ERERER8hIKdiIiIiI9QsBMRERHxEQp2IiIiIj5CwU5ERETERyjYiYiIiPgIBTsRERERH6FgJyIiIuIj/h8aUjz8GWei1gAAAABJRU5ErkJggg==", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAnYAAAHWCAYAAAD6oMSKAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAACBN0lEQVR4nOzdd3iUZbrH8e/MpPfeE0JvgYReFOkCKmLbtdMsu4plRVF0XdzVo7hi16yuBbGLKCLqiihSpNeA9BY6IQnpPZmZ88dAICaE9En5fa5rrsO8z/M+7/26c/D2qQar1WpFRERERJo8o70DEBEREZG6ocROREREpJlQYiciIiLSTCixExEREWkmlNiJiIiINBNK7ERERESaCSV2IiIiIs2EEjsRERGRZsLB3gE0ZRaLhRMnTuDp6YnBYLB3OCIiItIMWa1WsrOzCQsLw2isvE9OiV0tnDhxgsjISHuHISIiIi3A0aNHiYiIqLSOErta8PT0BGz/oL28vOwcjYiIiDRHWVlZREZGluYdlVFiVwtnh1+9vLyU2ImIiEi9qsq0Ly2eEBEREWkmlNiJiIiINBNK7ERERESaCc2xExERqWdms5ni4mJ7hyGNlKOjIyaTqU7aUmInIiJST6xWK0lJSWRkZNg7FGnkfHx8CAkJqfW+uErsRERE6snZpC4oKAg3NzdtZi/lWK1W8vLySE5OBiA0NLRW7SmxExERqQdms7k0qfP397d3ONKIubq6ApCcnExQUFCthmW1eEJERKQenJ1T5+bmZudIpCk4+zup7VxMJXYiIiL1SMOvUhV19TtRYiciIiLSTCixExERkWZt4sSJXHPNNfYOo0Fo8YSIiEgDuumdNQ36vC/uHlDte5KSknj22Wf54YcfOH78OEFBQcTFxfG3v/2N4cOHAxAdHc3hw4cBMBqNBAcHM2bMGF588UV8fX0BWLZsGa+88grr168nKyuL9u3bM23aNG699dYLPvvQoUO0bt0ao9HIkSNHCA8PLy07efIkkZGRmM1mEhMTiY6OrtL7vPbaa1it1mr/c2iKWnSP3dGjRxkyZAhdunShe/fuzJs3z94hiYiI2NWhQ4fo1asXv/76K7NmzeL3339n0aJFDB06lClTppSp+/TTT3Py5EmOHDnCp59+yooVK3jggQdKy1evXk337t35+uuv2bZtG5MmTWL8+PF8//33F40jPDycjz76qMy1Dz/8sEyiV1Xe3t74+PhU+76mqEUndg4ODrz66qvs3LmTxYsX87e//Y3c3Fx7hyUiImI39957LwaDgfXr13P99dfToUMHunbtytSpU1m7dm2Zup6enoSEhBAeHs7QoUOZMGECmzdvLi1/4okneOaZZxg4cCBt27blwQcfZPTo0cyfP/+icUyYMIEPPvigzLUPPviACRMmlLlmNpu54447aN26Na6urnTs2JHXXnutTJ0/DsUOGTKEBx54gEcffRQ/Pz9CQkL45z//WcV/Qo1bi07sQkNDiYuLAyAkJISAgADS0tLsG5SIiIidpKWlsWjRIqZMmYK7u3u58sp6vY4fP853331Hv379Kn1GZmYmfn5+F43l6quvJj09nZUrVwKwcuVK0tPTGTt2bJl6FouFiIgI5s2bx86dO5kxYwZPPPEEX375ZaXtf/jhh7i7u7Nu3TpeeOEFnn76aX7++eeLxtXYNerEbsWKFYwdO5awsDAMBgMLFiwoVyc+Pp7o6GhcXFzo168f69evr9GzNm3ahNlsJjIysvo3H1oNFkuNnisiItJY7N+/H6vVSqdOnapU/7HHHsPDwwNXV1ciIiIwGAy8/PLLF6z/5ZdfsmHDBiZNmnTRth0dHbntttuYPXs2ALNnz+a2227D0dGxXL1//etf9O7dm9atW3PrrbcyadKkiyZ23bt356mnnqJ9+/aMHz+e3r17s2TJkiq8dePWqBO73NxcYmNjiY+Pr7B87ty5TJ06laeeeorNmzcTGxvLqFGjSo/lAIiLiyMmJqbc58SJE6V10tLSGD9+PO+8806N4ixe8RLMvxN2LIDi/Bq1ISIiYm/VXWAwbdo0EhIS2LZtW2lSdOWVV2I2m8vVXbp0KZMmTeLdd9+la9euVWp/8uTJzJs3j6SkJObNm8fkyZMrrBcfH0+vXr0IDAzEw8ODd955hyNHjlTadvfu3ct8Dw0NLZM/NFWNelXsmDFjGDNmzAXLX375Ze66667SzP/tt9/mhx9+YPbs2UyfPh2AhISESp9RWFjINddcw/Tp0xk4cGCN4sz7KAnzZQV4pv8X162fQ4dR0HkcuOsIGRERaTrat2+PwWBg9+7dVaofEBBAu3btSu999dVXGTBgAEuXLmXEiBGl9ZYvX87YsWN55ZVXGD9+fJXj6datG506deLmm2+mc+fOxMTElPv3+hdffMEjjzzCSy+9xIABA/D09GTWrFmsW7eu0rb/2PNnMBiwNIPRt0bdY1eZoqIiNm3aVOaHYzQaGTFiBGvWVG0pudVqZeLEiQwbNozbb7/9ovULCwvJysoq8wHw3rCdoJdWYn11P6eWHSd73edYv54MK16E0wdq9oIiIiINzM/Pj1GjRhEfH1/hYsKMjIxK7z97xml+/rnRq2XLlnHllVfy73//m7vvvrvaMU2ePJlly5ZdsLdu1apVDBw4kHvvvZcePXrQrl07Dhxouf/ubbKJXWpqKmazmeDg4DLXg4ODSUpKqlIbq1atYu7cuSxYsIC4uDji4uL4/fffL1h/5syZeHt7l35K5+Pdfjs4OeG2cz/B767C5bnfSfkxhfSN/8Oy8AFY9AQc2wgtZA8dERFpuuLj4zGbzfTt25evv/6affv2sWvXLl5//XUGDCi7J152djZJSUmcPHmS9evXM23aNAIDA0tHwJYuXcqVV17JAw88wPXXX09SUhJJSUnVWqh41113kZKSwp133llhefv27dm4cSM//fQTe/fu5R//+AcbNmyo+T+AJq7JJnZ14dJLL8VisZCQkFD66dat2wXrP/7442RmZpZ+jh49ait48004cgRmzICAABxPnCLoizW4LnVil1MMKfs3UrJ4Biy4F/YuBnPtDvgVERGpL23atGHz5s0MHTqUhx9+mJiYGEaOHMmSJUt46623ytSdMWMGoaGhhIWFcdVVV+Hu7s7ixYvx97dNRfrwww/Jy8tj5syZhIaGln6uu+66Ksfj4OBAQEAADg4Vzx77y1/+wnXXXceNN95Iv379OH36NPfee2/N/wE0cQZrE9mK2WAw8M0335TuQ1NUVISbmxtfffVVmb1pJkyYQEZGBt9++229x5SVlYW3tzeZmZl4eXnZLubnw0cfwSuvwLvvwqBBZKcc4dg3bxFwcAVe7bxw9vCDTldBpyvBxave4xQRkYZXUFBAYmIirVu3xsXFxd7hSCNX2e+lwnzjAppsj52TkxO9evUqszTZYrGwZMmScl3FDcrVFf7yF9i5Ey69FADPwCg678wn8N+rKXnjICnLD5KzZg7Mmwhr/gOZx+0Xr4iIiDQbjXpVbE5ODvv37y/9npiYSEJCAn5+fkRFRTF16lQmTJhA79696du3L6+++iq5ublV2h+nNuLj40vnIFyQ8Q85c0EBODjgvm0v7tugKCqMlEEROKZ9jfee/2GI7Addr4WQmHqNXURERJqvRj0Uu2zZMoYOHVru+oQJE5gzZw4Ab775JrNmzSIpKYm4uDhef/31i+56XVeq0zUKwNGj8MYb8M47kJkJgNnXm9Mju2C6zB0fVydMQR2g63XQ6pLyyaGIiDQZGoqV6qirodhGndg1dtVO7M7KzobZs+HVV+HQIXjkETIemMyJVZ8QkLoRH1cTTt4h0Plq2554jq719QoiIlJPlNhJdbT4OXZNmqcnPPgg7NsH8+bBgw/iE9mZLjc9i0/UXRTFH+T0sh3krXwLvpwAG2dDbqq9oxYREZFGrlHPsWv2HBzghhvKXHJ690OctuzBYwsUto4gdVAwTmmf47nzWwzRg2zz8Pzb2ilgERERacyU2NVAlRZP1NRrr0Hr1jB7Ns6Jx3BOPIb5B1/SBrXG2P97vA8sxRgaCzHXQXgvMBjqPgYRERFpkjTHrhZqPMeuKtLTbYssXn8dTpwAoKRrRw4+NhrfzB34uDri6NcKul4DbYeBybHy9kREpEFpjp1Uh+bYNXe+vvDYY5CYCB9/DHFxONz/EB1ufxX3G97ioEtPTq/YRMGyV2z74SV8DgVZ9o5aRESkSfvnP/9JXFycvcOoMSV2jZ2TE9x2G2zeDGfOyXMLaUfHvFb4v7oBXtpF2k+7yVk9+8yGx/Ha8FhERGolKSmJ+++/nzZt2uDs7ExkZCRjx44tcyhAdHQ0BoMBg8GAyWQiLCyMO+64g/T09NI6y5YtY9y4cYSGhuLu7k5cXByffvpppc8+dOhQabsGgwE/Pz8GDx7Mb7/9Vm/v25xojl1TYTCAyXTue0oKuLjgsvcILnuPUPKjP2mDIjD2T8drz48YteGxiEjjNOeqhn3exO+rVf3QoUNccskl+Pj4MGvWLLp160ZxcTE//fQTU6ZMYffu3aV1n376ae666y7MZjN79+7l7rvv5oEHHuDjjz8GYPXq1XTv3p3HHnuM4OBgvv/+e8aPH4+3tzdXXVX5P4dffvmFrl27kpqayrPPPstVV13F3r17CQ4Orv4/gxZEPXY1EB8fT5cuXejTp4/9gnjiCThyBP71LwgMxOHUafy+2orXU+vJ+vwwp3csp+R/j8H3D0Hib2Cx2C9WERFpMu69914MBgPr16/n+uuvp0OHDnTt2pWpU6eydu3aMnU9PT0JCQkhPDycoUOHMmHCBDZv3lxa/sQTT/DMM88wcOBA2rZty4MPPsjo0aOZP3/+RePw9/cnJCSEmJgYnnjiCbKysli3bl1p+ccff0zv3r1LY7jllltITk4uLV+2bBkGg4ElS5bQu3dv3NzcGDhwIHv27CnznOeff57g4GA8PT254447KCgoKFNusVh4+umniYiIwNnZmbi4OBYtWlRafraH8csvv2TQoEG4urrSp08f9u7dy4YNG+jduzceHh6MGTOGlJSUqv2PUAtK7GpgypQp7Ny5kw0bNtg3kMBAmDHDluC9+y507owxNw8fhyicxs3ioEsXkhO3U7TkOZh/J+xYAMX59o1ZREQarbS0NBYtWsSUKVNwd3cvV+7j43PBe48fP85333130dOfMjMz8fPzq3JM+fn5fPTRR4DtnPiziouLeeaZZ9i6dSsLFizg0KFDTJw4sdz9f//733nppZfYuHEjDg4OTJ48ubTsyy+/5J///CfPPfccGzduJDQ0lP/85z9l7n/ttdd46aWXePHFF9m2bRujRo3i6quvZt++fWXqPfXUUzz55JNs3rwZBwcHbrnlFh599FFee+01fvvtN/bv38+MGTOq/N41pVWxtVCvq2JrwmKBn36yJXy9ewNQ+PtazLfcQmF/b5x7BeHm6Q0dR9tOtXAPsHPAIiLN1wVXOTbiodj169fTr18/5s+fz7XXXltp3ejoaE6ePImjoyNms5mCggL69evHokWLLpgAfvnll9x+++1s3ryZrl27Vljn0KFDtG7dGldXV4xGI3l5eVitVnr16sWaNWtwdKx4F4iNGzfSp08fsrOz8fDwKD2W9JdffmH48OEA/O9//+PKK68kPz8fFxcXBg4cSI8ePYiPjy9tp3///hQUFJCQkABAeHg4U6ZM4Yknniit07dvX/r06UN8fHxpvO+99x533HEHAF988QU333wzS5YsYdiwYYCtZ3DOnDllhrLPp1WxUp7RCGPGlCZ1AM6fLsBteyK+7yXg+MwW0r/bTebKT7F+dQesmAWnD9gxYBERaUyq29czbdo0EhIS2LZtW+nCiiuvvLLCfV6XLl3KpEmTePfddy+Y1J1v7ty5bNmyha+//pp27doxZ86cMkndpk2bGDt2LFFRUXh6ejJ48GAAjhw5Uqad7t27l/45NDQUoHTIdteuXeV6GAcMGFD656ysLE6cOMEll1xSps4ll1zCrl27Lvics/MAu3XrVuba+UPF9UWLJ5q7+++3Lbx4+20cT6TgOzcFy48eZF7SCsugdLz3L8MUpg2PRUQE2rdvj8FguGCv0h8FBATQrl270ntfffVVBgwYwNKlSxkxYkRpveXLlzN27FheeeUVxo8fX6W2IyMjad++Pe3bt6ekpIRrr72W7du34+zsTG5uLqNGjWLUqFF8+umnBAYGcuTIEUaNGkVRUVGZds5PBg1n/h1nqYd55xU954/X6uO5f6Qeu+YuPBxmzoSjR+GNN6BNG4xZOfj8uAOfF7aSWBJG8t71FC36Byy4F/b+BOZie0ctIiJ24Ofnx6hRo4iPjyc3N7dceUZGRqX3m87s3pCff24+97Jly7jyyiv597//zd13312juG644QYcHBxK57/t3r2b06dP8/zzzzNo0CA6depUo96wzp07l1mQAZRZIOLl5UVYWBirVq0qU2fVqlV06dKlBm9S/5TY1UCjWBVbXR4ecN99sHcvfP019O+PccIk2t39Lp7Xv8ohj1hSVm4id+nL5214nGnvqEVEpIGdPTKzb9++fP311+zbt49du3bx+uuvlxmmBMjOziYpKYmTJ0+yfv16pk2bRmBgIAMHDgRsw69XXnklDzzwANdffz1JSUkkJSWRlpZWrZgMBgMPPPAAzz//PHl5eURFReHk5MQbb7zBwYMHWbhwIc8880y13/XBBx9k9uzZfPDBB+zdu5ennnqKHTt2lKkzbdo0/v3vfzN37lz27NnD9OnTSUhI4MEHH6z28xqCErsaaDSrYmvCZILrroM1a+CFFwBwDY+hQ9TVBL68AdOsHaQu2k7GyvexzJsEq9/UhsciIi1ImzZt2Lx5M0OHDuXhhx8mJiaGkSNHsmTJEt56660ydWfMmEFoaChhYWFcddVVuLu7s3jxYvz9/QH48MMPycvLY+bMmYSGhpZ+rrvuumrHNWHCBIqLi3nzzTcJDAxkzpw5zJs3jy5duvD888/z4osvVrvNG2+8kX/84x88+uij9OrVi8OHD3PPPfeUqfPAAw8wdepUHn74Ybp168aiRYtYuHAh7du3r/bzGoJWxdZCo1sVWxvvvQcPPABnus9LgvzIuDQMBgTg4+eGQ6uzGx53u0hDIiICOitWqqeuVsVq8YTY3HknXHMNvPUWvPEGDskpBMxPw7LYjYyBrSganotP4hpcgjvaErzoS8FoumizIiIi0nA0FCvnBATAP/4Bhw/DO+9Ax44Yc/Lw21eA29h/cdihNcf2byPn5+fg6ztgxzdQlGfvqEVEROQM9dhJea6ucNddcMcd8MMPYDbj1XkwXp0HU3BoK0W33kJSTDLOycfx3vIpxo6jofM48Ai0d+QiIiItmhI7uTCjEcaOLXPJZXkCLqt34rUaCqNDSRnojyE5Dd/t3+LY5jLbMG1AOzsFLCIi0rIpsauB+Pj40uXgLc7gwbZFFu+/j/OhkwQfOknJYl/SB4ZS3D8T372/4hoRa0vwIvtqw2MREZEGpFWxtdCsVsVWV3o6/Pe/8NprkJQEgMXdlYMzxuBoSsPP3QmPoFYYul4LbYeDg9NFGhQRaV60KlaqQ2fFin35+sL06XDoEHzwAXTtinHAJbSb9hX+42ZywrU9R7b8zulfXsb85UTY8inkZ9g5aBERkeZNQ7FSO87OMHEiTJgAaWlgMODRrj8dPaKx3h1FfodwTvVLxphyCr+tX+LUYQR0uQZ8Iu0duYiISLOjxE7qhsEAZ3YaB2DFCgxmC27bD+G2HYpahZDWP4DilK/x3/4/3Nqc2fA4NNZ+MYuIiDQzSuykfvz5z9C/P7z+OrzzDk6Hkwg5nIT5Z2/SBoaT3C8b331r8ArvaJuH1/oybXgsIiJSS5pjJ/UnKgpefBGOHrWdSxsejiktk8BfDuJ16X0kGYM5sHsrqT8+R8mXk+D3r6Ao195Ri4i0eBMnTsRgMPDXv/61XNmUKVMwGAxMnDixTF2DwYCjoyPBwcGMHDmS2bNnY7FYytwbHR1dWtfd3Z2ePXsyb9680vJ//vOfxMXFVRhTcXExjz32GN26dcPd3Z2wsDDGjx/PiRMn6uy9mwMldlL/vL1h2jQ4eBA++gj+9S/8htxEx798RMifZmFYks2Jlds5+etbFHw+Hta/CznJ9o5aRKRFi4yM5IsvviD/zBniYFu5+dlnnxEVFVWm7ujRozl58iSHDh3ixx9/ZOjQoTz44INcddVVlJSUlKn79NNPc/LkSbZs2UKfPn248cYbWb169UXjycvLY/PmzfzjH/9g8+bNzJ8/nz179nD11VfXzQs3ExqKlYbj5AS3317mkkcKeCxYhz+Q36UVp/t6Y075FL+E+Xh0GGKbhxfYwS7hioi0ZD179uTAgQPMnz+fW2+9FYD58+cTFRVF69aty9R1dnYmJCQEgPDwcHr27En//v0ZPnw4c+bM4c477yyt6+npSUhICCEhIcTHx/PJJ5/w3XffMXDgwErj8fb25ueffy5z7c0336Rv374cOXKkXLLZUqnHrgbi4+Pp0qULffr0sXcoTZ+fn21FraMjrjsPEz5nG6Gv7SL356McXP09GfPuw/K/R+HIWtCWiyIiDWry5Ml88MEHpd9nz57NpEmTqnTvsGHDiI2NZf78+Res4+DggKOjI0VFRTWKLzMzE4PBgI+PT43ub46U2NXAlClT2LlzJxs2bLB3KE1fu3YwZw4kJsJjj4G3N47Hkgn+agetZm0j45iF/b+vI+XbJyn+6i7Y8yOU1OwvABGRRiE398KfgoKq1z1viLTSurVw2223sXLlSg4fPszhw4dZtWoVt912W5Xv79SpE4cOHaqwrKioiJkzZ5KZmcmwYcOqHVtBQQGPPfYYN998c8s7JKASGoqVxiE8HJ5/Hv7+d3j/fXjlFUxWK9FPfU/6vhUkr/mc9F278DjxIn5+H+LS7WrodCW4eNs7chGR6vHwuHDZFVfADz+c+x4UBHl5FdcdPBiWLTv3PToaUlPL16vFaEdgYCBXXnklc+bMwWq1cuWVVxIQEFDl+61WK4Y/HC352GOP8eSTT1JQUICHhwfPP/88V155ZbXiKi4u5s9//jNWq5W33nqrWvc2d0rspHHx9IS//Q3uu8+22MLZBd+Yy/HtOBRLm2jygl1J7pMGye/it+kLPLqMss3D8wqzd+QiIs3S5MmTue+++wDbVKTq2LVrV7n5eNOmTWPixIl4eHgQHBxcLvG7mLNJ3eHDh/n111/VW/cHSuykcXJwgA7nLZpYuRLjsRN4HAOPTVDQIYL0vj4kp36Jf8J3eHW4FEO36yGos/1iFhGpipycC5eZ/rCfZ3IlOwQY/zCb6gJDnrU1evRoioqKMBgMjBo1qsr3/frrr/z+++889NBDZa4HBATQrl27GsVyNqnbt28fS5cuxf/8jfEFUGInTcXQobBzJ7zyCnz0ES57jxG69xglIf6kDQjmQOoifHcvxzuqOw7dr4eoAbbTMEREGht3d/vXrQaTycSuXbtK/1yRwsJCkpKSMJvNnDp1ikWLFjFz5kyuuuoqxo8fX63n5efnk5CQUOaap6cnUVFR3HDDDWzevJnvv/8es9lMUlISAH5+fjg5OVX/5ZohJXbSdHTuDO+8A888A//5D8TH45B0mqBvTuN09XOk5ieQunMjPolb8QluhXPsddBuBDg42ztyEZEm7WLDnYsWLSI0NBQHBwd8fX2JjY3l9ddfZ8KECRj/2LN4EXv37qVHjx5lrg0fPpz33nuPhQsXApTbxHjp0qUMGTKkWs9prgxWq/aQqKmsrCy8vb3JzMzUGL895OXZNjxeu9a2stZiIWffcnJn/Yt8j2wcWnnj5xeAW/erodNV4Opj74hFpAUpKCggMTGR1q1b4+LiYu9wpJGr7PdSnXxDPXbSdLm5wV//avsAGI14+MXg8ek6KCggv3s0ab1TSUlJxm/jXDy7Xg5drwPvcPvGLSIiUk+U2EnzkpcHY8fC11/juu0Q4dugODqE9H5+JKfOwy/he7zbD8TY/QYI7mrvaEVEROqUEjtpXlq1gi+/tG14/Oqr8P77OB5KIuhQEv6/eJN8SxdSsn/Fd89KvKNicOp+PbS6RAstRESkWdDJEzWgI8WagNat4bXX4OhRmDkTQkMxZecT+ujnhI56iFzHAA5t38TJb5+i4ItJsOs7KCm0d9QiIiK1osUTtaDFE01IURFs2ACXXGL7brViGdyf/IJU0vp4QVtffH398Ii9GjqPBVdf+8YrIk2eFk9IdWjxhEh1ODmdS+oADh7E+Nt63AH3DVDYKYL0vt6kpKTgt2EuXjGjMHS9Fnwi7RayiDQP6j+Rqqir34kSO2mZ2raF33+Hl1+GTz7BefcxQnYfoyQsgMwBARxM+RKfMwstHLpfDyHd7B2xiDQxjo6OAOTl5eHq6mrnaKSxyztzJvDZ301NaSi2FjQU20ycPAlvvglvvQXp6QDkP3kHKaGpFGSn4+3qiE9UV5xjr4dWl5Y/xkdE5AJOnjxJRkYGQUFBuLm5VftcVGn+rFYreXl5JCcn4+PjQ2hoaLk61ck3lNjVghK7ZiYnx7bR8RdfwJIlYIKC3YvJmvMa+ZYMHKK88QmKxD3uOmh/OThqzoyIVM5qtZKUlERGRoa9Q5FGzsfHh5CQkAqTfyV2DUSJXQtgNkOHDnDwIAWxUaT186G4vS++fv54nV1o4eZn7yhFpJEzm80UFxfbOwxppBwdHS94Di9o8YRI3cnIgLg4SEzEZesRwrYeobh1MJkD/DiQfGahRdeRmLpdDz5R9o5WRBopk8lU6b+4ReqKeuxqQT12Lcj+/fDKK/DBB5CfD4Al0IfT41qT0d4THzdnvNv1wyn2BgiNtXOwIiLSnGgotoEosWuBTp+Gt9+GN96AU6dgyU/keSaTvmk+eRmpeLg64hPZGde4GyB6EBj1X+giIlI7SuwaiBK7FqygAH78Ea65xnYcmbkY8503UrRrC6f7+mFt7Y1PYDiePW+A9qO00EJERGpMiV0DUWInpXJyIDTU9n+Boi7hpPfzJa+TP77+fnjHjcPQ5WqdaCEiItWmxRMiDc3DA1assG14/MUXOO08TvDO45gj/MkYGMCBU6n4rZ+Ld7fRmLpdB94R9o5YRESaIfXY1YJ67KRCx47B66/Df/8LWVkAFNx0Gan9HMgrLMHH3QmfDpfiFPsnCO5i52BFRKSx01BsA1FiJ5XKyoJ334X//AeWLwcPCzkbPiN36UKyLSW4h3riHdUdt55/hqj+trl6IiIif6DEroEosZMqsVjKHkM2agTWJUvJ7xHMqYEhOET74B3cCq9eN0LbYeDgZL9YRUSk0dEcO5HG5PykrrAQSqwYzBbcNp6k9caTFHbyJ/mSFE4fT8TH/318el6HofNV4Oxpv5hFRKRJUo9dDcTHxxMfH4/ZbGbv3r3qsZPq27QJZs2CefNsPXpAcaQXycOjyO/sh4+XJ76xV2GKuRY8g+0crIiI2JOGYhuIhmKl1hITbSdavP8+5OVheXwi6dFppGdkYLZY8fVwwbvzMJzj/gT+be0drYiI2IESuwaixE7qzOnT8N578OCDYDRj3bmQ3PhXKD54iuR+oXgFu+PduhduPW+EiF72jlZERBqQErsGosRO6o3FAp07w969WB2M5PQKJmlAKC6R3niHtcer95+h9WAwaZqsiEhzp8SugSixk3pjscC338ILL8DatQBYDZDfxZ+kS8IxtPfDOyAE395/wtBhNDi52TlgERGpL0rsGogSO6l3ViusWmVbaLFwYenl7CGRHBnZCmcHIz7ePvj2uNq20MLNz47BiohIfVBi10CU2EmD2rULXnoJPv4Y3n+OfKc9pJ88QFZaPkYHI96+Hvh0uxzn2BvAJ8re0YqISB1RYtdAlNiJXSQnQ2Cg7c+HV1Hyt/sx/LqT0/1DSesXileAO97t+uPe60YI6WbfWEVEpNa0QbFIcxYUdO7PUQNxOGiG7CKCfj5M4PJjpPcO5vjATFz3rMYrsgvefW6CVpfoyDIRkRZAiZ1IU2Y02jY7/vJLmDULw9at+K0+ge/ak2R3CyDl0nTSDu/AKygSvz5/xtB+JDg42ztqERGpJ8aLVxGRRs3REW69FbZsgcWLYcQIDBYrXltTCN9xGrPFQtKxRA58N4uUD26mZNMnUJBp76hFRKQeaI5dLWiOnTRaW7bAiy/C1Duwpq8ka/9qsnalYDiZQ25sIN4+Xvh0H4NL7A3gFWrvaEVEpBJaPNFAlNhJk5GWCOOugpU7KfFxJnVAGBl9Q/DydcO742DbQovADvaOUkREKqDFEyJSlm80XHk77H0Zh+QUQn5MJGjpEU73DeXYgGzcdi7Fs1UsPn1vgcg+9o5WRERqSD12taAeO2lyCgps++C98G/YfwAAq8lAymURpIxohZuTA15h7fDvd4vtyDKjpuGKiNhbdfIN/a0t0pK4uMBdd8GevbBgAQzoi8FsxTfMB183J/KLzSQl7uLA/H+R+tF4rLu+h5Iie0ctIiJVpB67WlCPnTQLq1dDhyg4tJjiHd+Tu3g3jquPkzIkkpJO/vj4B+HX5waMna8CJ3d7Rysi0uJo8UQDUWInzU5hNsR2gz2HASgIcyflskgKugfi4+uDX89xOpNWRKSBKbFrIErspFk6dgxefAHeeQfyCwEoDHAldVA4Ob1C8PX1xLf7aBy7/0lbpYiINAAldg1EiZ00a6mp8Nqr8MZrkJkDQFa3AI7e1AmT0YCfhwvenYfh3OPP4NfGvrGKiDRjSuwaiBI7aRGys+G//4VZz2P+a3/SAgtJyynEkluM0QDeAe74tOuPS8+bICTG3tGKiDQ7SuyqKCMjgxEjRlBSUkJJSQkPPvggd911V5XvV2InLUphoe34smMbMG/9gsL4BbgsPUpavxDSLonAI8Qdn1bdcet5E0T2BYPB3hGLiDQLSuyqyGw2U1hYiJubG7m5ucTExLBx40b8/f2rdL8SO2mxrFYY2AfWbgLA4mAgo2cwpy+LwDXME9/wdrj3vBHaDAGjyb6xiog0cTp5oopMJhNubm4AFBYWYrVaacF5rkjVGQywegP88AM8PQPjhi34rU/Cd0MSmd0DOXlZBs7H9uMbPAePnn+G9peDg7O9oxYRafYa9QbFK1asYOzYsYSFhWEwGFiwYEG5OvHx8URHR+Pi4kK/fv1Yv359tZ6RkZFBbGwsERERTJs2jYCAgDqKXqSZMxjgqqtg3SZYvhyGXYbBCj5bUwhcn0RWQTGHDx/i8P9eJvuT22HrF1CYY++oRUSatUad2OXm5hIbG0t8fHyF5XPnzmXq1Kk89dRTbN68mdjYWEaNGkVycnJpnbi4OGJiYsp9Tpw4AYCPjw9bt24lMTGRzz77jFOnTjXIu4k0GwYDXHYZLFkOmzbBdVfjOfUuwvy9cHYwUXw0i9PLd5K4+G0yP74NNrwPeWn2jlpEpFlqMnPsDAYD33zzDddcc03ptX79+tGnTx/efPNNACwWC5GRkdx///1Mnz692s+49957GTZsGDfccEOF5YWFhRQWFpZ+z8rKIjIyUnPsRCqSl4Z1xzcU//UfOG0+SX6YO6mDIymKDcLfxwPvmFEYYq4H73B7Ryoi0qi1iLNii4qK2LRpEyNGjCi9ZjQaGTFiBGvWrKlSG6dOnSI7OxuAzMxMVqxYQceOHS9Yf+bMmXh7e5d+IiMja/cSIs2Zmx+G3pNxuuwGrC7OuJ7IJfLz3US8uIHcxfs4sHwe6Z9NxrJ0Jpw+YO9oRUSahSab2KWmpmI2mwkODi5zPTg4mKSkpCq1cfjwYQYNGkRsbCyDBg3i/vvvp1u3bhes//jjj5OZmVn6OXr0aK3eQaTZMxjgldcxHD0GTzwOXu44p+YTPn8fUS9sIPfXRPavX0TaF3/FvOjvcHKrvSMWEWnSWvSq2L59+5KQkFDl+s7Ozjg7a2WfSLUFBMCzz8H0x+Gt/8CLL+CUkoZroYVMs4WTmQWkJKzEf/96fKJicIj9E0QN0F54IiLV1GQTu4CAAEwmU7nFDqdOnSIkJMROUYlIpTw94dHH4IEH4cMP8R8YhWn7N2Se2IdhcxIczSZxUDbeh3/HN6yt7TzaNkPA1GT/qhIRaVBN9m9LJycnevXqxZIlS0oXVFgsFpYsWcJ9991Xr8+Oj48nPj4es9lcr88RabZcXOAvfwHAJ2Y0PkfWYX7jKkyHT+O35gQZPYM5fFkWHscP4Bf4EU6x10OH0doLT0TkIhr1qticnBz2798PQI8ePXj55ZcZOnQofn5+REVFMXfuXCZMmMB///tf+vbty6uvvsqXX37J7t27y829qw86eUKkjlit8P338K8nYdM22yUjZHYL5PTgSFxb++Dv749zt2ug81hw9rRvvCIiDajZHCm2bNkyhg4dWu76hAkTmDNnDgBvvvkms2bNIikpibi4OF5//XX69evXIPEpsROpY1YrrFgB/3wSlq0svZwyJJKUy1vh5eKIv483rl2vgK7XgXvVjv8TEWnKmk1i19gpsROpR5s2wTMzsC78kdQ7Yklu7QGAocSCp7sT/t7uuHUaAd3+pL3wRKRZq/ezYhcuXFjte0aOHImrq2tNHiciLVGvXrDgBwz79xMY5I7zprlkbvsR5+9247krjZTBEVjTcwjYtRj3DkMh9ibw0d6SItKy1ajHzmis3vZ3BoOBffv20aZNm+o+qlE6f/HE3r171WMn0lCyU7C064gxOR2AwgBXUgdHUNI3lCBfd1w7DDmT4EXZN04RkTpU70OxRqORpKQkgoKCqlTf09OTrVu3NpvE7iwNxYrYQUoKvDwL65vxGHLyACjydSb1sgjMAyMI9HPDpd1ltgTPN9q+sYqI1IF6P1JswoQJ1RpWve2225T4iEjdCAyEmS9gOH4Snvs/LH7eOKUXEvbtAdwW7OVAcjbHExZTNP9eWPocpCXaO2IRkQajxRO1oB47kUYgLw/efgvLC89z4i89yHIpwGq14pBTjJePCwEB7jhGD4DYm8G/rb2jFRGpNq2KbSBK7EQaEbMZsJDz+w+krPkUv7eX474/ndOXRsDwVvgHuOHQqj/E3gIB7ewdrYhIlTV4YldcXExSUhJ5eXkEBgbi5+dX2yYbNS2eEGnkcnOwdO2E8fBxAEpcHUgfGAYjo/EL9MAU1ceW4AV2sHOgIiIX1yCJXXZ2Np988glffPEF69evp6ioCKvVisFgICIigssvv5y7776bPn361OglmgL12Ik0YiUl8NmnmP/5JKbEYwCYnU2kDwjDMKo1vsEeGCN62RK8oE52DlZE5MLqffHEyy+/THR0NB988AEjRoxgwYIFJCQksHfvXtasWcNTTz1FSUkJl19+OaNHj2bfvn01ehERkRpzcIDxEzDtOwSff465XTSmQjMBy45SsuwI+5JzSNu7BssPD8PiJyF5l70jFhGptRr12N188808+eSTdO3atdJ6BQUFzJkzBycnJyZPnlzjIBsr9diJNCEWCyxYQMmsZzh8Y3sK85OwWq24n8rDJ9gd71BPDGFxtkUWITH2jlZEpJQWTzQQJXYiTZTVSvauX0n67QMinv4W51O5ZPcJwXBlOzwjPDGExELczRDSzd6RiojU/1Ds+ZYuXXrBsv/+97+1bV5EpO4ZDHh2GU77G17DMbQ1RrMV77Un8XxqJdn/2UzO+t9g0ePw42Nwcqu9oxURqbJaJ3ajR49m2rRpFBcXl15LTU1l7NixTJ8+vbbNN0rx8fF06dKlWS8MEWkR/P1x2LgNli6lZEAvDBYrXhuScP/nSnLe2EDehtXw09/hf4/C8c32jlZE5KJqPRS7evVqxo8fj4eHB5999hmJiYnccccddOzYkY8++ohWrVrVVayNjoZiRZqZNWsonj4VxxVrAThySyesPUMI8nLG1dEBAjvZhmjDe9k5UBFpSRp8jl1OTg5//etf+eqrr7BYLDzzzDM8+uijGAyG2jbdqCmxE2mmNm6k8I1ZHLjEC07vwWq1ErQvA69wT5zb+UFAB4i7BSJ62ztSEWkBqpNvONTFA/fu3cvGjRuJiIjgxIkT7Nmzh7y8PNzd3euieRGRhtW7N84fzqULkHlgA8d/ehvfeR/jkFtMfhd/HMadxjF1LwS0t62ijexr74hFRIA6mGP3/PPPM2DAAEaOHMn27dtZv349W7ZsoXv37qxZs6YuYhQRsRvvtn3ocsNMLCNGYDUYcN15GseZaymYuZqS39bAkqdh4QNwZC1okwERsbNaD8WGhoYye/ZsxowZU3qtuLiYJ554gtdff53CwsJaB9lYaShWpIXZt4/CJx7B6ZvvMZgtABS19cF4e1ccWvmAX2tbD17UAGjmU1FEpOE06By71NRUAgICKixbvnw5gwcPrk3zjZoSO5EW6tAhCv8+Dce58zGaLRx8oCfubX0J8HDGZDSAbytbgtfqEiV4IlJrDbqPXUBAAL/99hu33XYbAwYM4Phx26HbH3/8MSaTqbbNi4g0PtHROH86D+Ohw+T++19k9RjM6dwi9iVnkzN/N+Zf1sKvM+HbKZC4QkO0ItJgap3Yff3114waNQpXV1e2bNlSOvSamZnJc889V+sAGyPtYyciAERE4P7oDOImvkzIbf/F7NAJ9x8OYnpzMyUzVmBZvBaW/hsW3AsHlynBE5F6V+uh2B49evDQQw8xfvx4PD092bp1K23atGHLli2MGTOGpKSkuoq10dFQrIiUkZZG/ozHcXp/DqaCIgDMIe4Yx7bD0C8MfCMh9kZoPQSMtf7vahFpIRp0KHbPnj1cdtll5a57e3uTkZFR2+ZFRJoOPz9c3/wvpuMnyX/ofkrcnDEl5WJ4dyvmx5dh/X0H/PYyLLhHPXgiUi9qndiFhISwf//+ctdXrlxJmzZtatu8iEjT4+eH68uv43DiFHmPPUyxhytkFXEIK9kFxZB1HFa8aJuDd2ilEjwRqTO1TuzuuusuHnzwQdatW4fBYODEiRN8+umnPPLII9xzzz11EaOISNPk7Y3b8y/ieDKZ3LlfkNJ2OIfTC0lMyab4v1vgh9W2RRYL77ftgyciUku1Pnli+vTpWCwWhg8fTl5eHpdddhnOzs488sgj3H///XURo4hI0+bhgde4G+jLDZxOOkLyS9NxXLsC1p7A8sN+jNclQVoiBLSDHrfrqDIRqbE6OSsWoKioiP3795OTk0OXLl3w8PCoi2YbNS2eEJEayc0l/Z9P4vGft3DMs+0kYI3ywnBDR4gJhKBOEHcrhPe0c6Ai0hg06AbFLZkSOxGplcxMTj10LwGffYmpsAQAa3tfDPf2BB8XCO4KPW6DkG52DlRE7KneE7upU6dWue7LL79c3eYbvfj4eOLj4zGbzezdu1eJnYjUSuHx46Tedweh3/9Msa8zOU8OwNfTBePZUytCukPP2yGos30DFRG7qPfEbujQoWW+b968mZKSEjp27AjA3r17MZlM9OrVi19//bW6zTcZ6rETkbqUtn07+79+CxennbgaLQS5OuI1fy+Gka0h3BPCetjm4AV2sHeoItKAqpNv1GjxxNKlS0v//PLLL+Pp6cmHH36Ir68vAOnp6UyaNIlBgwbVpHkRkRbJLyaGvjHx7D2wn0OL4/H4YSHey49iXXEUw4BwGJcHJ7ZARB/bEK1/W3uHLCKNTK3n2IWHh7N48WK6du1a5vr27du5/PLLOXHiRK0CbMzUYyci9cVqtbLt808Ie/7vBP5+1HbNaMAwOBLGtgdfF4gaAD1uBd9o+wYrIvWqQU+eyMrKIiUlpdz1lJQUsrOza9u8iEiLZDAYiL3ldtw3HmTtyzNJbx+IwWKFpUewPrYUvtgJh1bDt/fBsn9DxlF7hywijUCtE7trr72WSZMmMX/+fI4dO8axY8f4+uuvueOOO7juuuvqIkYRkRbLzcmB/g9Np2DdAVb+/QFyWvlgKLZQcCyLkrOVDv0G395rO80iq/mOkojIxdV6KDYvL49HHnmE2bNnU1xcDICDgwN33HEHs2bNwt3dvU4CbYw0FCsiDW3b/iNkv/oYgaZ9WELd8Pdwxr/IgnHdCRgWDS6O0HYYxN4MnsH2DldE6oBd9rHLzc3lwIEDALRt27ZZJ3RnKbETEXswW6ws35BA3up3aJ2/nfDvDuCz9iRWb2cMY9vB4ChwcoR2I6D7TeARaO+QRaQWtEFxA1FiJyL2lFVQzC+//kLHj/5F+8UJOKXbTrEgwBXGtYcB4eDoBB1GQ/c/g5uffQMWkRpRYtdAlNiJSGNwJDWXlT98Rv/5rxH5614cc2zTYgh1hxs6Qc8QMDlCxyug25/A1ceu8YpI9SixayBK7ESkMVm/7wT7fnyfIT9+RPDyQzjkl1ByeTQON5+3HZWDM3S6EmJuABf9vSXSFCixayBK7ESksSkqsbB4004Kfn2HQUsWktE/BM9Ad4I8nXE4lgP5xdDRHxxdofNY6HotOHvaO2wRqYQSu3qms2JFpLFLzy3i++Wr8dnxMW2K92I0GGg7ZztOe9OhawBc3xFa+4CjG3S9BrpcA05udo5aRCpSr4ldeno6VqsVPz8/UlJS+O233+jYsWO5kydaAvXYiUhjt+9UNr/88j+6HJpLt4Vr8N14yrbRMUDPYLiuo+0cWicPiLkOOl8Nji72DVpEyqi3xO69997jueeeA2DatGl8+umnxMbGsmLFCh588EHuvPPO2kXexCixE5GmwGq18tveZH5fNo9LD8yj1S+78U5IxmAFDEC/MLi2AwS5g4u3LcHrdJVtPp6I2F29JXbdu3dn3bp15OfnExUVRWJiIoGBgWRmZjJ48GASEhJqG3uTosRORJqSgmIz323cT+b6z7kk8UfCfjmA947TtsKJ3Wz7353l6mNbQdvxCtuKWhGxm+rkGw7VadjBwQFXV1dcXV1p164dgYG2TS+9vb0xGAw1j1hEROqdi6OJPw3oSHLX6Xy9YhS+kV/Q++AK/NedJLuLL0HFZlwcTXAwAwIKIf9d2P41dL8R2o8CU7X+lSEidlCt/y81mUwUFBTg4uLC8uXLS6/n5OTUeWAiIlI/grxcuOeqAWzv0YVvlq0gJnoeUcWJZKfk4OPkQNhbmzHkFMHI1jC6Dax9y5bgxd1qO65M/yEv0mhVayj2bBfgH3vnkpOTOXz4MH369KnzABszDcWKSFNnsVj5ZdcpNq/8kUsyFhKceozIz3fjevzMf7C7O8KYNjAiGpwdwKcV9LwdovrbNW6RlkTbnTQQJXYi0lzkFJbw9fpE0rcs4NKcxQTtOEbwL0dwTs6zVfByhmvaw2WRYDJCYEfoOQFCu9s3cJEWoEETu5kzZxIcHMzkyZPLXJ89ezYpKSk89thjtWm+UVNiJyLNzbH0PL5YuRPf/fPpk7sKv61JBP16BKe0AluFxwdAh/POnA2Lg16TwL+tXeIVaQkaNLGLjo7ms88+Y+DAgWWur1u3jptuuonExMTaNN+oKbETkeZq0+E0Fv62mW5J8+mStxXf9Ul4JOViujsON6cz07OzCm09eQDRl0KP28E73H5BizRT9bYqtiJJSUmEhoaWux4YGMjJkydr27yIiNhBr1Z+xEYM43/bu/DZ2t8YfNm3hJUchdRcvFwdCbaC099XQPcguL4TsBIOr4Z2I2yLLNz97f0KIi2SsbYNREZGsmrVqnLXV61aRVhYWG2bFxERO3EwGbk6Noxpt1/Hgb5Ps8DzZrKMPmTlF5Oy+jjWQjOsPwlPLIPPd0J2AexbDPPvhA3vQUGWvV9BpMWpdY/dXXfdxd/+9jeKi4sZNmwYAEuWLOHRRx/l4YcfrnWAIiJiX95ujvx1SDv2dwnmo9/6EXT0Ry7ptYyCEDdCfjqE+/4MWJwIvx2Fq9rByGjYsQD2/gRdr4Ou1+qYMpEGUus5dlarlenTp/P6669TVFQEgIuLC4899hgzZsyokyAbK82xE5GWxmq1smxPCt+u2U7v9O+JK9iIx740Qn86hPPJXFulIDd4djA4nBkUcvG2bXLc8QptcixSA3bZ7iQnJ4ddu3bh6upK+/btcXZu/mcMKrETkZYqt7CEeRuPkrBtCyNzFhJduB/vhGRCfjkCfUJwuDWm/E0eQRB3C7Qdrk2ORarBbvvYnW2qpRwvpsRORFq6I6fz+GB1IubDaxmZ+z3+BacwWq14+boR5OmMw+Es+GYv/KkTRJ75e9InyraCttUA+wYv0kRUJ9+o9eIJgPfff5+YmBhcXFxwcXEhJiaG9957ry6abpTi4+Pp0qVLiztpQ0Tkj6L83XhqbFdGjx7H3PDH+dHnWnKcPUnPLWJfcg5FX+6C31Pgqd/g/a2Qlg8ZR2Dps/D9VDi5zd6vINKs1LrHbsaMGbz88svcf//9DBhg+6+vNWvW8Oabb/LQQw/x9NNP10mgjZF67EREzikoNrNgy3F+STjAwNzF9M5fg+vpHEJ+Pozn76m2Sk5GuLw1XNEWXB1t18LioOdECGhnr9BFGrUGHYoNDAzk9ddf5+abby5z/fPPP+f+++8nNTW1Ns03akrsRETKS8osYM7qQxw5tJcRuT/QsWgnrkeyCPvpEC6HzmyB4uEEN3aGSyPO3djqEts5tN4RFTcs0kI16AbFxcXF9O7du9z1Xr16UVJSUtvmRUSkiQnxdmH6mE5sOhzMR2siWX96OyPbfE/+nZ547kojdPEhHFPyodhc9sbDq+DI+ZscB9jnBUSasFr32N1///04Ojry8ssvl7n+yCOPkJ+fT3x8fK0CbMzUYyciUrlis4Xvt53gm83H6ZK7jqF5P+FRnIHPtlQcLwnH39sVk8EAO1Jtw7Ttz5xDa3K0bY/S/UZw0d+v0rI16FDs/fffz0cffURkZCT9+/cHbOfEHjlyhPHjx+Po6Fha94/JX1OnxE5EpGpScwr5ZO1hNh04ycD8ZQzIX46jtRgHk5EgF0d8/m81htR86BUCN3SEEA/bjY6uZzY5vsb2Z5EWqEETu6FDh1apnsFg4Ndff63NoxodJXYiItWz/XgmH64+REbqCYbnLSKmcAum/GLCFx/Cc0MSBitgMsDgKBjXHrzO7Inq4g3d/wwdr9Qmx9Li2G0fu5ZGiZ2ISPWZLVYW70hi3qZj+OQmMjLvO6KKD+F8Kpfwn4/guuu0raKLA1zRBi5vA84m2zX3QNv8u7bDwFgnO3aJNHpK7BqIEjsRkZrLzC/m8/VHWL4nhU6FWxmR9z98zGm4H8wk/KdDOB7LtlV8tB90/sNCCu9I2wraVgMbPnCRBtYgid3kyZOrVG/27Nk1ab5JUGInIlJ7+5Oz+WDVIQ4lZ9C3YCWD8pbgbC7AZ3sq/sdycLojFuPZE41S8iDA9dyRZAEdoNcECI213wuI1LMGSeyMRiOtWrWiR48eVNbEN998U5PmmwQldiIidcNqtbJ0TzJfrD9KSV46g/N+plfBOgxYcHIwEuLlgmehGR5bBm194E+dIdr7XANhcdBzAgS0t9MbiNSfBknspkyZwueff06rVq2YNGkSt912G35+fjUKuKlSYiciUrdyC0uYt/Eoi3eewq84icvzvqdt0R4AAvemE/jpTgwlZ/61NSAcrusAAW7nGmg1EHqO1ybH0qw02By7wsJC5s+fz+zZs1m9ejVXXnkld9xxB5dffjmGs93kzZgSOxGR+nHkdB4frE5k18ks2hbtYWTu9wSaT+GYXkDE0qO4bTplq+hghJHRcFU7cDuzvZbBAG2HQ4/bwd3fbu8gUlfssnji8OHDzJkzh48++oiSkhJ27NiBh4dHXTTdaCmxExGpX6v3p/LJusOk5xTQs2Adg/N/xt2Sg9uJHMIXH8ZpX7qtorczvDAUnEznbnZwgdiboMs12iJFmrQGPVLsLKPRiMFgwGq1YjabL36DiIjIRQxsF0DPVr58s+U4P2xz4HfnHgzKX0K/8FXsm+COz/5MQn5KxNgtEMP5SZ3VCiUFsGkO7PsZ+v0Fwnva7T1EGkqdDcWuXLmSq666ikmTJjF69GiMLWB/IfXYiYg0nJOZ+Xy4+jAJR9PxMZ9meN6PdCncBhYrbgYDwYHuuDk5wKFM+Hg7/LkTdDxvKDaqP/S5CzyD7fcSIjXQIEOx9957L1988QWRkZFMnjyZW2+9lYCAlnVgsxI7EZGGt+lwOh+tOcSprAIiig8xKvc7wkqOAuDl6kj4h9sxbk2xVR4QDjd1PneChckJut0AMTeAg5Od3kCkehpsu5OoqCh69OhR6UKJ+fPn16T5JkGJnYiIfRSbLXy/7QTfbDlBUXEJ3Qq3MCxvEV6WDBxyiolccQzX1cdtR5S5O8KfOsGgSDCe+feVRzD0vcvWiyfSyDVIYjdx4sQqrXz94IMPatJ8k6DETkTEvlJzCvlk7WHWHjyNg7WIAfkruCR/GY7WIjxP5BD+7QFMZ0+waO8L47tBhOe5BsJ7Qd+7wTvcPi8gUgU6UqyBKLETEWkcth/PZM7qQxxLz8PDksXQ3J+ILdyIwWwhfNMpvBclYig0w61dYETrsjcbTdD1Wuh+Ezi62OcFRCpRr4ldeno6VqsVPz8/UlJS+O233+jYsSNdu3atVdBNkRI7EZHGw2yx8tOOJL7adIy8ohIiig9xdc48/M0puGQVEbEpGeebu5wbjs0ttg3TnuXmD33ugNaX2ecFRC6gOvlGtZauvvfee/Tq1YvevXvz1ltvce2117JkyRJuuukm3nvvvVoFbU95eXm0atWKRx55xN6hiIhIDZmMBq7oFsorN8YxpGMQxx1b81+fv7HadTD5Xi7sHxrB0Yx8is0WKDLD0yvhP5shvcDWQN5pWP4CLHoC0g/b92VEaqhaPXbdu3dn3bp15OfnExUVRWJiIoGBgWRmZjJ48GASEhLqMdT68/e//539+/cTGRnJiy++WOX71GMnItJ47TqZxTsrDnIyM5+w4qOMy5lLgDkZk9FAxNEcPN5OAIsVXB3g+o4wtNW53jyDETqPhbhbwcmt0ueI1Ld667FzcHDA1dUVPz8/2rVrR2BgIADe3t5N9gixffv2sXv3bsaMGWPvUEREpA51DvXi39d3Z2xsGElOUbzj8zdWug6lxGLgcLg7Jx7siaW1N+SXwCc74P9WweFM281WC+z8Fr65G/YvsW14LNIEVCuxM5lMFBTYuqyXL19eej0nJ6duozpjxYoVjB07lrCwMAwGAwsWLChXJz4+nujoaFxcXOjXrx/r16+v1jMeeeQRZs6cWUcRi4hIY+LkYOTWfq14elwMYX5eLHUfw/s+95HsEEJ6gCt77uhG7p86YnV1gMRM+NdK+HwnmC22BvIzYOUr8OOjcPqAXd9FpCqqldj98ssvODvbNnn09vYuvZ6Xl8c777xTt5EBubm5xMbGEh8fX2H53LlzmTp1Kk899RSbN28mNjaWUaNGkZycXFonLi6OmJiYcp8TJ07w7bff0qFDBzp06FDnsYuISOPRLsiDmdd14/qeEaQ4RfGu94OscBtBicHIobhAjj7cB3PvELACp3LPDcmelbwLvnsQ1vwHCrPt8g4iVdFktjsxGAx88803XHPNNaXX+vXrR58+fXjzzTcBsFgsREZGcv/99zN9+vSLtvn444/zySefYDKZyMnJobi4mIcffpgZM2ZUKSbNsRMRaXqOnM7jreUHSEzNIbjkBONy5hJcchKDwUD40Wy8WvtiCDwzry67yLbQwt/1XAPOXtBzPHQYBU10GpI0LQ22j11SUhIhISE1vb1a/pjYFRUV4ebmxldffVUm2ZswYQIZGRl8++231Wp/zpw5bN++vdLFE4WFhRQWFpZ+z8rKIjIyUomdiEgTY7FY+W7bCb7adAxzSTGX5C9lUN4STJhxdTQR5uOKi6MJ3t8KG07CtR1gRDSYzhvo8m8H/e+BwI52ew9pGept8cQfXX755bW5vVZSU1Mxm80EB5c9zDk4OJikpKR6eebMmTPx9vYu/URGRtbLc0REpH4ZjQbGxYXzwg3daR/iw29uI3jP5wGSHMLJLzZzMDWX5LQ8rMm5UGiGL3bB06vgYMa5Rk7vhx8ehpWv2ubiiTQCtUrsmsgobpVMnDjxoludPP7442RmZpZ+jh492kDRiYhIfQj1duWfV3dl4sBoMl0ieM/7Pn51G02J1UhKQTEHJ3ejaHxX20bGR7JsK2c/3g55xeca2f8LfPMX2LkQLBb7vYwItUzs7LnFSUBAACaTiVOnTpW5furUqXobHnZ2dsbLy6vMR0REmjaDwcDomFBm/ak7MRF+rHIbxrs+D3LCIZICs4V9HX1Jfrw/1gHhtsUVvx6GJ5aX7b0ryoX178B390PSdnu9ikjtEjt7cnJyolevXixZsqT0msViYcmSJQwYMMCOkYmISFMU5OnC36/swl8ua0uuWwSzvaewxG0MZoMDKQYr+69pS/5DvSHE3bavXbB7+UbSD8Oi6bB8FuSlNfxLSIvnYO8AKpOTk8P+/ftLvycmJpKQkICfnx9RUVFMnTqVCRMm0Lt3b/r27curr75Kbm4ukyZNqte44uPjiY+Px2w21+tzRESk4Q3tFERspA/vrzzI6sND2ePUlatzviSi5AgHA1zwe6QPQbklmM6eM2u1wtoT0CcUHM70lyQuh6PrIPYm6HINmBr1v26lGanVqtgePXqwZcuWuoynjGXLljF06NBy1ydMmMCcOXMAePPNN5k1axZJSUnExcXx+uuv069fv3qL6Xza7kREpHlbtT+VOasPkZ1fRP+C3xiS9xOO1mIcTUbCfFzwcHaENcfhnQQI94Dx3aCDX9lGvMKh318gvKdd3kGavgbb7qSlU2InItL8ZeYXM2fVIdYcTMXPnMLYnHlEFR8CwMfNiZC96Zg+3gE5RbYbLouEP3UCD6eyDUUNgL53gUdQw76ANHlK7BqIEjsRkZZjw6E0Zq9MJD23kL4FqxiW9yOO1mIcjEbCHIx4LjwAv53ZLcHTCW7qDAPCy25ibHKCbn+CmOvBwaniB4n8QYMmdhs2bGD69OmkpKTQrl074uLiSj9RUVG1abrROn+O3d69e5XYiYi0ELmFJXy89jDL9iTjaz7N2Jx5tCo+CICXqyOhp/Jw+HgHnDhzhvrgKJjYrXxDHsHQ926IapipQ9K0NWhi16lTJ6Kiorj66qtLFzckJCSQnp6Or68vp0+frk3zjZp67EREWqZtxzL474qDnM4uoHfBGkbk/Q9HaxEmo4EQdyd8lh+D7/bBQ32hk/+FGwrvZZt/5xXWcMFLk9OgiZ27uzvbtm2jbdu2Za4fPnyYhIQExo0bV5vmGzUldiIiLVd+kZnP1h/hl52n8DKf5uqceUQXHwDAw9mBMJMJRx+XczesPGo7c7ZzQNmGjA7Q9VrbCloH5wZ8A2kqGuxIMYBLLrmEY8eOlbveqlWrZp3UiYhIy+bqZOKOS1szY2wX3PzC+djrbn7wuI5CgzM5hSXszy8iLbcIK1ZIzbOdWPHCOng3AbLOnTuOpQR+n2c7veLoeru9jzQPNeqxu+666+jevTuxsbFYrVb+85//MG/ePHx9fesjxkZLPXYiIgJQWGJm3sZj/O/3k3iUpDM2Zx5tivcB4O7sQJijCadv98PSw7bTK9wdbStnB0WC8Q+nOLUbbpt/51TBBsjSItX7UOy0adNISEhg69atpKamAuDv78+4cePo378/PXr0oFu3bjg5Ne8VP0rsRETkfPuTc3h7+QGOpecRV7Cey3O/x9lagMFgINjTGb/kPAwfbredOwvQ3hcmdINwz7INufnDJQ/Y5uBJi9egc+yOHz9eumDi7OfgwYM4ODjQsWNHtm3bVpvmGyWtihURkQspMVuYv+U4CxNO4FqcxlW582lXtBuwDd+GeTrjsuwofLMXCs3gbIIXh5Xf9w6gwyjocyc4ujbwW0hjYvd97HJyckp79KZMmVLXzTca6rETEZELOXI6j7eWHyAxNYfYgo1cnvsdLtZ8DAYDgR5OBBRaMHy209Zbd33HCzfkEQSXPAihsQ0XvDQqdk/sWgoldiIiUhmLxcp3207w1aZjOBdncGXOfDoU7QTAxdFEmI8rribjuXl2RzLhl8O2zY3dHMs21ukq6D1JK2dboHpfFXvkyJFq1T9+/HhNHiMiItKkGY0GxsWF8+/ruxMWGs5cr4ks8LiJfIMrBcVmElNzOJVTiMVqBasVZm+znV7x5Ar4PblsY7u/h2/vg1M77fMy0iTUKLHr06cPf/nLX9iwYcMF62RmZvLuu+8SExPD119/XeMARUREmrowH1f+dXVXJg6MZq9HH97yfYTdTjFYrZCaU8jBlBzyis1wa1cIdof0Anh5gy3Ryys+11D2SfjxUdjwHpQU2e+FpNGq0VDs6dOnefbZZ5k9ezYuLi706tWLsLAwXFxcSE9PZ+fOnezYsYOePXvyj3/8gyuuuKI+Yrc7DcWKiEh1JWcV8M6Kg2w/kUnXwgRG536LmyUXgAAPZ4KcHDB8sxd+TrRtjeLrApO7Q0xg2Ya8I+DSqRDYoeFfQhpUg82xy8/P54cffmDlypUcPnyY/Px8AgIC6NGjB6NGjSImJqamTTdqWhUrIiK19evuU3yy9ggUZHBFzjd0LvodADcnByJ8XXE8kAHvb4XkPNsNj/SFrn9I7gwGiLkB4m4Fk0PDvoA0mAZJ7GbMmMG4cePo1avl7rGjHjsREamNtNwi3vvtIJuPpNO5cBtjchfgbsnBwWgg3NcVDwzw1R44mgWP9i+/mfFZvq1svXf+bSsulyatQY4UO3bsGGPGjCEiIoJ77rmHRYsWUVSk8X4REZGq8nN34tHRnbhvaHuOevXiLZ+H2e4cR4nFyuHTeZwqLMF6Sxd4pN+5pK7QDF/vgfyScw2lH4YfpkLCZ2Ax2+dlpFGo1VCsxWJh1apVfPfdd3z77becPHmSkSNHMm7cOK666ir8/PzqMtZGRz12IiJSVzLzi/nPsv1sPZpBh6KdXJkzHw9L1rmhWdOZvpjPdsDPh8Df1Tb3rktA2Yb82sCgqeAb3dCvIPXEbvvY7dq1qzTJ27RpE3379uXqq6/m5ptvJjw8vK4e02gosRMRkbpktVqZu+EoCxKO42wpYGTed/Qo2HBuaNbZEXafhtlbISXfdtPQKPhTZ3A9b46d0WSbdxdzAxhrPDgnjUSj2KA4JSWFhQsXsnDhQgYNGsQjjzxSH4+xKyV2IiJSH9YdPM1byw9QUGymTdFersz9Gh9zOoEezgR6OmMoNMO83fDrYdsNAWd67zr/ofcuoANc+hD4RDb8S0idaRSJXUugxE5EROrL0bQ8Xlq8h6SsAhythYzI/ZFeBWvwcDIRfnZodmcqfLANUs/03k3oBkOiyjZkcoSe46HLNbZVtNLkNMjiiYs5evQokydPrq/m7So+Pp4uXbrQp08fe4ciIiLNVKSfG89e243YSB+KDc786HENH3n/haNmHw6m5JBTWGKbX/f0ZbbhWHdHiA0q35C5GDa8Dz8+BlknG/5FpEHVW4/d1q1b6dmzJ2Zz812dox47ERGpb+fPuwNwsBYxJO9n+uf/RpCno21oFgNkFYLXeefIrjkOPYLB5by5dw7O0GsSdLpSvXdNSHXyjRrvZrhw4cJKyw8ePFjTpkVEROQMg8HATX2jiA5w561lBygsceIX9yvZ6dSdq3O+JLroNBE+rjicn9RtPQXvJECgm23uXSd/2/WSQlj3NhxeDZf+DTwq6OGTJq3GPXZGoxGDwUBltxsMBvXYiYiI1JGjaXm8uHgPp7IKADBazVyW/wuDC5cT6eOMh/OZ/ppdqfDeVkiz1WNENNzQEZzP689xdIXed0DH0Q37ElJtDTLHLjQ0lPnz52OxWCr8bN68uaZNi4iISAXOn3cHYDGYWOY2irc972Ndpg/J2QVYsdpWx/7fZXDZmdWwvxyCGb/B3rRzjRXnw5o34ecZkHu6wd9F6keNE7tevXqxadOmC5ZfrDdPREREqs/D2YHpoztxTdy5/WFPOYTxnvd9zDUPJTGtiBKLBVwdYVJ3mNoXfF1sZ84+vwYW7ivb4PHN8O29sP+XBn4TqQ81TuymTZvGwIEDL1jerl07li5dWtPmRURE5ALOzrv724gOODuYALAaTKxyG8bLLlNYlRlAbuGZI8e6Bdp67wZFgBWI8CzfYFEurHwVljwNeWnly6XJ0D52taA5diIiYm9/nHcHgNXKgMKV3GRaRpC7wbZqFuBIJkR5n6uXmAHhnuBkOnfN2RP6/QXaDGmA6KUqtEFxA1FiJyIijUFOYQlv/LqPrUczylz3Nacyie/o6XwMhz8eLZZRAE+uAE8n28rZ9n84373VJdD/HnD1qdfY5eLqffHEtm3bsFgsVa6/Y8cOSkpKavKoRkkbFIuISGNydt7duLiy57KnmwJ42TiRd4tGk2P+ww5np/PB0QhJuTBzDXyxE4rO28ni8Cr4doptaxRpMmrUY2cymUhKSiIwMLBK9b28vEhISKBNmzbVDrAxU4+diIg0NmsOnObt5QcoLCm73ZiPJYOHPX6iXcm+c0OzucXw+U5Ydcz2PcQd7oiFdr5lG209GPr/1TZMKw2u3odijUYjd999N25ublWq/5///IedO3cqsRMREWkAR07n8dLPf5h3d8a1nju5rmQRjub8cxcTTsGHv0NGIRiA0W3gT53Knk7h6gsDH4BIjVY1tHpP7IYMGYKhmkeRfPbZZ4SGhlb3UY2aEjsREWmscgpLeGPJPrYeyyhXFumSz2PeSwhI33reDUXw2U7bUWTDWsHtMRU33G449L0bnNzrJ3ApR4snGogSOxERacysVitfbDjKt2fOmT2fyWjg3lbHGJj6NYbC7HMFW5Oho9+5M2YzC8HNARzPWznrHmDrvQvvWc9vIKDErsEosRMRkabgQvPuAPqFOnCv6884H6tgkYTFCrPWQlaRbe5dG5+y5R1GQ587wdGlfgIXoIGOFBMREZGmYUBbf54ZF0OwV/kEbN3JEqYmj+ZI3EPltzZJzYcTuXAiB55dDV/thuLzksO9i+D7v0H6ofoMX6pBPXa1oB47ERFpSiqbd2cyGrglzp8r8r/DcODX824qgk93wNoTtu/hHrbeu9Y+593sZNvUuMOoeo2/pdJQbANRYiciIk2N1Wrl8/VHWbi1/Lw7gJ5RvtzXMRO3jW9Bbuq5gk1J8NHvtmFZowGuaAPjOoDDeYN/bYbAgPs0NFvHGmQo9uDBgygnFBERaVoMBgO39IviweHnzpk93+Yj6Uxb48j+AS/Y5tCd1SsE/m8w9A21zb3bkgx/zAMOLrMNzaYl1us7yIXVOLFr3749KSkppd9vvPFGTp06VSdBiYiISP06O+8uyLN879rpnEL+uSiR7z2ug1HPgWeIrcDTCe7pCVN6wp2xZVfKnpV5DH54GPb+VM9vIBWpcWL3x966//3vf+Tm5tY6oKZAR4qJiEhzEOXvxnPXdSM2wqdcWYnFwidrDzNrmxM5o1+DLuPOFfYOhWjvc9+XHIIfD5zrwTMXweo3YMUsKM5HGk6N59gZjUaSkpIICgoCwNPTk61btza70yUqozl2IiLSHFxs3l2AhzN/G9GedhyDVa/ZeuXOOpIF//wNrMCIaLi5i20O3lle4TBkOvi1rtd3aM4aZI6dwWAod/pEdU+jEBEREfu72Ly71JxC/rlwJ/875QNXvwHdbgDDmRQi0hP+3Nn2518OwVuboei8LVGyjtuGZvcsqvf3kFr22I0ZMwZnZ2cAvvvuO4YNG4a7e9kjRubPn1/7KBsp9diJiEhzc+R0Hi8u3kNydvlzZgH6RPvx18Ftcc85DCtfObeH3boT8N5WKLFAe194oDd4OJW9uc3gM6tmXev3JZqZBtnuZNKkSVWq98EHH9Sk+SZBiZ2IiDRHle13BxDo6cKDw9vTzt8FVr8OZ/e9230aXt8I+SUQ4g4P94UAt7I3a2i22rSPXQNRYiciIs3VxebdORiN3NovijExIbaFEvsW2wqOZ8NL6yG9AG7rCsOjy99scoK+d0PH0eXLpBwldg1EiZ2IiDR3lZ0zC3Btj3Bu7B0Ja9+CPf+zXUwvsA3NjmoNlc2/bz0YBmpo9mKqk2841PQhkydPrlK92bNn1/QRIiIiYmcD2voT7uN6wXl332w5jo+bE6MG3AtGE+z6DnxdYPR5u2TkFcPOVNs2KedLXA6n98GQxzU0W0dqtXiiVatW9OjRo9ITKL755psaB9fYqcdORERaipzCEl5fso9tFcy7M2DgbyPa06+NP2x4D3YsOFdYYoFX1sPO03BtBxjbrnwvnskR+v5FQ7MX0CA9dvfccw+ff/45iYmJTJo0idtuuw0/P7+aNiciIiKNmIezA9NHd+LzDUf4buuJMmVWrLy5dD+eLo506XMnGB3h93m2QqMBon1sid03eyEtH26PAdN5O66Zi2HNm5D0OwyYAk5/WHAhVVbjfezi4+M5efIkjz76KN999x2RkZH8+c9/5qefftIZsiIiIs2Q0Wjg1n6tKtzvrths4cXFezhyOg96TYDYm87cZIA/dYLbu4IBWH4U3tgEhSXlH5C4/MxZswfr/V2aqzpbPHH48GHmzJnDRx99RElJCTt27MDDw6Mumm60NBQrIiIt1YX2u/Nzd+LpcTEEeDjD1i9gyyfnCjcnwdtboNhiO5LsoT7g5Vy+cZPjmVWzY+r5LZqGBjl5olxDRiMGgwGr1YrZXPHKGREREWkezp4z28qv7MEEablFPPe/XeQUlth67XpNPFfYMwQe7Q8ejnAoE/6z+dz5suczF8OaeFj+AhTl1e+LNDO1SuwKCwv5/PPPGTlyJB06dOD333/nzTff5MiRI82+t05ERKSl83B24KGRHXBzKjtl/0RGPi8s2m3bIqXbDdDnjnOF7Xzh7wMhygtu7Vr5diiJKzQ0W001Hoq99957+eKLL4iMjGTy5MnceuutBAQE1HV8jZqGYkVERGB9Yhov/7yn3PVerXx5eGRHjEYD7Poe1r19rtBqLZvUZRaCdwXDsmAbmu1zF3S6oo4jbxoaZINio9FIVFQUPXr0wFBJtq2zYkVERJq/j9ce5odtJ8pdH9oxiL8Mbmv7smeRbfXrH+1Ns51WcWNnGNbqwg9pfZntrNkWtmq2QbY7GT9+fKUJXXMWHx9PfHy85hKKiIiccUvfKA4k57A7KavM9aV7kvF1c+LPfSJt+9QZTbDqtbI3b06CIjN8vB1O58MNHSseok1cAaf3w+DHwL9tPb5N06UjxWpBPXYiIiLnpOcWMX3+NjLzi8uVTb6kNZd3DbF9OfArrHzl3MIJqxW+3Wf7AAwIh8ndweECSwFa2NBsgwzFZmVlXbwSNOuER4mdiIhIWduPZ/Lc/3Zh+UN6YcDAQyM70Lf1mcMMDi6H314Cq+VcpRVH4cPfwWKFLgFwX09wdbzww6IHwcD7m/3QbIPNsatsKNZqtWIwGJr1cKUSOxERkfIWbDnOFxuOlLvuaDLyxBWd6Rx65t+Zh1bCillgOS9X+D0Z4jdDoRkiPc9sj+J04Yd5hTX7odkGmWO3dOnS0j9brVauuOIK3nvvPcLDw2vapIiIiDQD4+LC2Hsqm81H0stcLzZbePGnPfzz6q5E+rlB9KVgdIBlM88ld92CYPoAeGUDhHiAWyU9dgBZJ+CHh6HvXdDpynp6o6ajzubYeXp6snXrVtq0aVMXzTUJ6rETERGpWE5hCU/M/73cyRTwh9MpAI5ugKXPguW8Y8ZO54OXEziayt1/Qc10aNYuJ0+IiIiInOXh7MDfRrTH0VQ+1UjLLWLm/3bbTqcAiOwDw2eA6bwhV3/Xc0mdxQrvb4V15bdTKePQb/Ddg3D6QB29RdOjxE5ERETqRZtADyYMiK6w7HhG3rnTKQDCe8KIp8omd2etOQ4rj9nOmf3pIqdQZJ+0Dc3u/qF2wTdRdZrYtdR97URERKRiI7oEM6h9YIVle09l8/qSfVgsZ2aFhcbCyKfBwaVsxQHhMPzMxsVf7ILPd9p68S7EUgJr34Jlz7e4s2ZrPMfuuuuuK/P9u+++Y9iwYbi7lz0MWCdPiIiItGyFJWae/GY7R9MrTrKGdwrmrsvOm6OfvAt+ngHF+eeuWa2w6CB8udv2vXcI3B138Tl4nqEwZHqTXjXbIHPsvL29y3xuu+02wsLCyl0XERGRls3ZwcTUyzvgeoEkbMnuU3y58ei5C0Gd4fL/A8fzFkEYDDCmLfwlDkwG2JgEL66HnKLKH352aHbX97V/kSZAJ0/UgnrsREREqm7NgdO8tmTvBcvvvLQNI7oEn7uQuh8WPwlFOWUr7kqFNzZBsQUe7Qft/aoWQPSlMPCBJrdqVqtiRUREpNEZ0NafMTGhFyyfvSqRDYfSzl0IaAejngXnPyQznQPgiYFwT4+qJ3Vg2xD5uwcg83g1I286lNiJiIhIg7m1XxQdgj0rLLNYrby+ZB+7k847ttS/LYx+Dlz+ML0rwhN6hpz7fiwbdqRePIDsJPjpccg6WYPoGz8ldiIiItJgHExGHhzeHi+Xik+UOHs6xdG08xZa+EbD6Jng6lNxoxkF8PJ6eGU9rD528SDy0mzJXXZSteNv7JTYiYiISIPy93Dm/mHtMVDxNmk5hSU8/+NuTucUnrvoEwWj/w1uFQy9ujtCBz8wW+HdrfD9ftsq2srkpsKixyEnuRZv0vgosRMREZEG1y3Cmxt6RVyw/HRuIc+dfzoFgHe4LblzDyhb2dFk2/pkzJktU77eAx9vB7Ol8iByU84kdyk1e4lGSImdiIiI2MV1PcOJjfS5YPnxjDxe/GkPRSXnJWheobbkziOobGWjAf7cGW7tCgZg6RF4cxOcnxhWJOcU/PSEbXi2GVBiJyIiInZhMBi4b2g7/D2cL1hnd1JW2dMpADyDbcmdZ0j5G0ZEw5Re4GiEhGRb793FZJ+ERdObRXKnxE5ERETsxtPFkYdGdMDBeOGUZOPhNGavSix70SMQxrwAXmHlb+gVAtP6QSc/uKZD1QLJOgE//R3yM6oefCOkxE5ERETsql2QB7cPaFVpnV92neKrTX9Y8ermB6OfB+/I8je094NH+4PbeatvT+eXr3e+zKO2YdkmnNwpsRMRERG7G9U1hEvaBlRa56tNR1my61TZi25+tq1QfCtIDA3nrbr96SA8sQwSTpWvd76MI7D471CQVXm9RqrFJ3bR0dF0796duLg4hg4dau9wREREWqy7LmtDuE/lx329vzKRjYf+MBfO1QdGzQS/1hXfZLHCrtNQZIHXN8LSw5UHkn7YdpRZYXbVg28kWnxiB7B69WoSEhJYunSpvUMRERFpsVwcTTw0sj0ujqYL1rFYrbz+6372JP0h6XLxglHPgX+78jcZDXBfLxgUAVbgo+3w60WSu7SDZ5K7nMrrNTJK7ERERKTRiPB1465BbSqtU1RiZtZPuzmWnle2wNnTdrZsYMfyNzkYYVL3c3vdfbwdfjtaeTCnD8DP/4CivMrrNSKNOrFbsWIFY8eOJSwsDIPBwIIFC8rViY+PJzo6GhcXF/r168f69eur9QyDwcDgwYPp06cPn376aR1FLiIiIjV1SbsALu9SwVYm58kpLGHmj7tJyy0qW+DkDiOfgaDO5W8yGOBPnWxbogB8sA3WHq88mNR98POMJpPcNerELjc3l9jYWOLj4yssnzt3LlOnTuWpp55i8+bNxMbGMmrUKJKTzx0PEhcXR0xMTLnPiRMnAFi5ciWbNm1i4cKFPPfcc2zbtq1B3k1EREQubPyAVrQL8qi0zumcQmb+bxe5f9yE2MkNRj4NwTHlbzIY4JYuMDjKNiybXnDxYFJ2wy//hOIq1LUzg9V6scPUGgeDwcA333zDNddcU3qtX79+9OnThzfffBMAi8VCZGQk999/P9OnT6/2M6ZNm0bXrl2ZOHFiheWFhYUUFp47ty4rK4vIyEgyMzPx8vKq9vNERETkwlJzCpn+9e/kFBZXWq9TiBdPXNEZJ4c/9FeVFMKSp+Hk1vI3WaywMxViAqseUHAMjPwXOFx4Q+X6kJWVhbe3d5XyjUbdY1eZoqIiNm3axIgRI0qvGY1GRowYwZo1a6rURm5uLtnZtsmXOTk5/Prrr3Tt2vWC9WfOnIm3t3fpJzKygn1zREREpE4EeDhz37B2GDBUWm93UhZv/rqPcn1VDs4w/CkI61H+JqOhbFKXXwz70ysP6NR2+OVftoSxkWqyiV1qaipms5ng4OAy14ODg0lKSqpSG6dOneLSSy8lNjaW/v37M378ePr06XPB+o8//jiZmZmln6NHLzLpUkRERGolLtKHa3uGX7Te+kNpvL8ysXyBgxMMnwERvS98c14xzFoHs9bC3oscK5a0DX59BsyV9yLai4O9A7CnNm3asHVrBd2zF+Ds7Iyzc8N2v4qIiLR0f+oVwf7kHLYdy6i03i+7TuHr5sT1vSLKFpgcYeiTsPx5OLK2/I2ORvBwsu1z98p6eKQftPW98INOJMCv/wfD/gGmxpVKNdkeu4CAAEwmE6dOld1B+tSpU4SEVL6SRkRERJoOg8HAfcPa4e9+8c6VeZuO8uvuCk6XMDnAkMeh1SXlyxxNtn3uOvtDgRleXg+HMyt/0PFNsPRZMJdUXq+BNdnEzsnJiV69erFkyZLSaxaLhSVLljBgwIB6fXZ8fDxdunSpdNhWRERE6o6XiyMPjmiPg/Hiqct7vyWy6XAFQ6pGEwx+rOJhWScTPNgb2vtCXgm8uA6OXeTkiWMbbL2AFnMV36L+NerELicnh4SEBBISEgBITEwkISGBI0eOADB16lTeffddPvzwQ3bt2sU999xDbm4ukyZNqte4pkyZws6dO9mwYUO9PkdERETO6RDsya39oi5az2K18tqS/ew9VUFiZjTCZY9WfLasswM81Adae0NOsW3O3cmLnDxxZC0sf6HRJHeNOrHbuHEjPXr0oEcP22qWqVOn0qNHD2bMmAHAjTfeyIsvvsiMGTOIi4sjISGBRYsWlVtQISIiIs3DmG6h9G/jf9F6RSVmXli0m+MZ+eULndxsq2VdvMuXuTrC1L4Q6WXb885ShV3hDq+CFS+CxVKFN6hfTWYfu8aoOvvKiIiISN3ILzLz929+50RmBUnbH/h7OPPMuBj83J3KFybvhp8er3iFa1ahbUg2xL3qgbUZDIMesSWEdahF7GMnIiIiLZOrk4mHRnbAycF00boXPJ0CIKgTXPK3im/0ci6b1O1Ng4yLnDxxcDmsfAXs2GemxK4GtHhCRETEviL93LhrUOsq1T2ansesn/ZQbK5gqLTNYIi9ufIGdqbaFlPMWmfryavMgV9h9et2S+6U2NWAFk+IiIjY36D2gYzoXLV59buTsnjj1/3lT6cA6HErtB584ZsD3Wz73J3IsSV4OUWVP2zfz7Am3i7JnRI7ERERabImDIymdYBHlequTzzN7FWHKi685EEI7FhxWaAbPNoPvJ3haDa8tN52WkVl9i6CdW9XKa66pMROREREmixHk5GHRrbHw7lqJ0D8vDOJ77aeKF/g4ATDngT3wPJlACEeMK2frefuUCa8sgEKLrI58e4fYN07VYqrriixExERkSYtyNOFe4e2w0DVVqN+tu5IxRsYu/razpV1dK34xnBPmNYX3Bxgfzq8thGKLrJ/3a6FsOG9KsVVF5TY1YAWT4iIiDQuPaN8GRcXVqW6Vqy8vmQ/h1Jzyxf6tbadTnGhLUuivOHhfuDiAJ5OYKxCMrljAWz8oEqx1Zb2sasF7WMnIiLSeFgsVp793y52nLjIOa9n+Ls783/XxOBb0R53OxZU3tN2IhuC3cFUjT6y7jdCz9urXv8M7WMnIiIiLY7RaOCB4e3xdasgUavA6dxCXly8h8KSCoZTu14DHcdc+OYwz3NJncUKK45CRdupnG/bXEj4rEqx1ZQSOxEREWk2vF0deXBEe0xVGSIFDqTk8NayAxVvg9LvrxAae/FGPtoOH2yD2dsufgRZwmewdW6VYqsJJXYiIiLSrHQK8eLmvlFVrr/24Gm+3Hi0fIHRBEMeB++IyhvoFmiba7f6uC3Ju9gsty0fw+9fVTm+6lBiJyIiIs3OVd3D6BvtV+X632w5zm/7UsoXOHvA8KfA2fPCN/cKgbvjwAAsPwKf7bx4crdpDuz4psrxVZUSuxrQqlgREZHG769D2hLi5VLl+u+sOMiepOzyBV6hMPTvth68C+kXBpPPDNv+cgjm7b54crfhfdi5sMrxVYVWxdaCVsWKiIg0bkdO5/Hkt9spqmiBRAW8XBz5v2tiCKooIdz3C6x6tfIGlh62DccCXNMBxrW/+EP73wOdrrxgsVbFioiIiABR/m5MviS6yvWzCor596I95Fe08XD7ERBzfeUNDG0FN3cBkwEiKhm+Pd/at2DPoirHWBkldiIiItKsDekYxNCOQVWufzwjj1d/2YulohWuvSZC1IDKG7i8NTw3xDb3rqrWvAn7fq56/QtQYiciIiLN3qRLWtM6wL3K9bcey+DDNYfKFxgMcNkj4N+28gaC3M79OTUP1hy/+ENXvQb7l1Q5xooosRMREZFmz8nByEMjOuDu5FDle37akcRPO5LKFzg421bKulVh1W1WIcxcA+8mwKpjF6+/8hU4uKzKMf6REjsRERFpEYK8XLhnSFsMVG3zYoAPVx9i69GM8gVufrbkzsG58gY8naBHMFiB97fCuhMXf+hvL0Hib1WO8XxK7GpA252IiIg0Tb2j/bgqNrTK9S1WK68t2cex9Lzyhf5tYdAjlTdgMMAtXeGySFty904CbK6gF/B8ViusmAWHVlU5ztLHabuTmtN2JyIiIk2PxWLlmR92sutkVpXvCfJ04f+ujcHLxbF84e9f2TYcrvShVnhvq22unYMRHugF3S6yoOPMyRdZPl203YmIiIhIRYxGA38b3gFfN6cq35OcXcBLP+2h2GwpX9jtBmg34iIPNcAd3aFPKJRY4I1NsDet8nssZlg2E45trHKcSuxERESkxfF2c+SB4e0xGas+327PqWzeWXGw4sIB90FwTOUNmIy2o8figiDE3fa5GIvZNueuipTYiYiISIvUOdSLm/pEVeue3/al8M2WCla3mhxg6BPgeZH5ew5GuLcnPDYAvC6y8OIsS0mV41NiJyIiIi3W2Ngw+kZXYduS83y54RhrD54uX+DiBSOeAqeL9MQ5msD9vLl6K4/BkarP96uMEjsRERFp0e4Z0o5Qb9cq17di5T/LDnAgJad8oXcEDHkcDFVMsTactG2DMmsdHM+ucgwXosROREREWjRXJxMPjeiAk4OpyvcUlZh58ac9nM4pLF8YFgf976laQ10DINobcopsyV1SbpVjqIgSOxEREWnxovzduPPS1tW6Jz2viFk/7aGg2Fy+sOMY6DLu4o24OcLDfSHCEzILYdZa2xFkNaTErga0QbGIiEjzc1mHQIZ3Cq7WPYdO5/Lmr/upcFvgPndCRBVyBQ8nmNYPQt0hrQBeWAtp+dWK4ywldjUwZcoUdu7cyYYNG+wdioiIiNShiZdE0zrAo1r3bDycxqfrjpQvMBhg8KPg2+rijXg5w7T+EOQGKfm2Ydn84mrFAUrsREREREo5mow8NLI9Hs4O1brv+20nWLo7uYIGXWH4P8HF++KN+LrAo/3B3xX6hoJL9WIAJXYiIiIiZQR5ujBlaDsMVH3zYoD3Vyay40Rm+QKPQBg+A0wVHEf2R/6u8PQguLajrcevmpTYiYiIiPxBjyhfru0RXq17SiwWXvl5LyczK5gfF9gRLn2oag25nZcAFprh6z1VjkGJnYiIiEgF/tQ7gm7hVRhCPU9OYQmzFu0hp7CC0yJaXwZxt1QviHe2wJLDVa6uxE5ERESkAgaDgfuGtcffvYpHf51xIjOfV37ei9lSwUrZuFug9eCqNza2PVwaUeXqSuxERERELsDb1ZG/jWiPg7F6KdOOE5nMXplYceGlf4PATlVrKNobbulS5ecqsRMRERGpRPtgT27rH1Xt+5bsPsUP206WLzA5wrAnwT2wDqIrS4mdiIiIyEWMjgllQJuAat/3ydrDbDqcXr7A1QdGPGXbDqUOKbETERERqYK/DG5DuI9bte6xYuWNX/dx5HQFx4T5RsPgx2q0rcmFKLETERERqQIXRxMPjWyPi6OpWvcVFJv590+7ycgrKl8Y0Rv63FVHESqxqxGdFSsiItIyRfi6cfegNtW+73ROIbN+2kNRiaV8YZeroeMVdRCdErsa0VmxIiIiLdfAdgGM6hpS7fsOpOTwn2X7sVor2Aal318gLK7WsSmxExEREamm2/u3on2QZ7XvW3vwNPM2HStfYDTBkMfBO7JWcSmxExEREakmB5ORB0e0x9OlCue//sH8zcdYtT+1fIGTu22lrLNXjeNSYiciIiJSAwEeztw/rB0Gqr+q9e3lB9h7Krt8gWcIDPu7rQevBpTYiYiIiNRQ9wgfbuhV9SO/zio2W3jxpz0kZxeULwzuCgMfqFE8SuxEREREauG6nuHERfpW+76sgmJeWLSH/CJz+cJ2w6Hbn6rdphI7ERERkVowGAzcN6wdAR7O1b73WHoery3Zh8VSwUrZnuOh1cBqtafETkRERKSWPJwdeGhkBxxN1U+tEo6m89GaQ+ULDAYY9DD4VX3fPCV2IiIiInWgbaAH4wdE1+jeRTuSWLwjqXyBg7Pt2LEqUmInIiIiUkdGdglmUPvAGt374ZpDbDuWUb7Aza/KbSixExEREalDdw5qTaSvW7XvM1usvPrLPo6l59X42UrsREREROqQs4OJqZd3wNWx+nvR5RWV8MKiPWQVFNfo2UrsREREROpYqLcrfx3ctkb3JmcX8PLivZSYLdW+V4mdiIiISD3o18afK7uF1uje3UlZvPPbwWrfp8ROREREpJ7c0q8VnUJqdvbrir0pLNhyvFr3KLGrgfj4eLp06UKfPn3sHYqIiIg0YiajgQeHt8fb1bFG98/dcJSNh9KqXN9gtVor2OpYqiIrKwtvb28yMzPx8qpZNi4iIiLN344TmTz7wy4sNUm7ivOZe9/wKuUb6rETERERqWddw7y5sU9kje4tLqn6IgoldiIiIiIN4OrYMHq18q3XZyixExEREWkABoOBe4e0I8jTpd6eocROREREpIG4OzswdWQHHE31k4IpsRMRERFpQNEB7ky+tHW9tK3ETkRERKSBDe0YxJCOQXXerhI7ERERETuYfElrov3d67RNJXYiIiIiduDkYGTqyA64OznUWZtK7ERERETsJMjLhXuHtsWAoU7aU2InIiIiYke9WvkxNjasTtpSYiciIiJiZzf1iaRLqHet21FiJyIiImJnRqOBB4e3x9fNqXbt1FE8IiIiIlIL3m6OPDiiPSZjzefbKbETERERaSQ6hXhxS99WNb5fiZ2IiIhII3Jl91D6tvav0b1K7EREREQamXsGtyXU27Xa9ymxExEREWlkXJ1MTB3ZAScHU7XuU2InIiIi0ghF+rlx56Wtq3WPEjsRERGRRuqyDoEM6RhY5fotPrFLTExk6NChdOnShW7dupGbm2vvkERERERK3dKv6qtk6+7U2SZq4sSJ/N///R+DBg0iLS0NZ2dne4ckIiIiUsrRVPV+uBad2O3YsQNHR0cGDRoEgJ+fn50jEhEREam5Rj0Uu2LFCsaOHUtYWBgGg4EFCxaUqxMfH090dDQuLi7069eP9evXV7n9ffv24eHhwdixY+nZsyfPPfdcHUYvIiIi0rAadY9dbm4usbGxTJ48meuuu65c+dy5c5k6dSpvv/02/fr149VXX2XUqFHs2bOHoKAgAOLi4igpKSl37+LFiykpKeG3334jISGBoKAgRo8eTZ8+fRg5cmS9v5uIiIhIXWvUid2YMWMYM2bMBctffvll7rrrLiZNmgTA22+/zQ8//MDs2bOZPn06AAkJCRe8Pzw8nN69exMZGQnAFVdcQUJCwgUTu8LCQgoLC0u/Z2VlVfeVREREROpNox6KrUxRURGbNm1ixIgRpdeMRiMjRoxgzZo1VWqjT58+JCcnk56ejsViYcWKFXTu3PmC9WfOnIm3t3fp52xCKCIiItIYNNnELjU1FbPZTHBwcJnrwcHBJCUlVakNBwcHnnvuOS677DK6d+9O+/btueqqqy5Y//HHHyczM7P0c/To0Vq9g4iIiEhdatRDsQ3hYsO953N2dtZ2KCIiItJoNdkeu4CAAEwmE6dOnSpz/dSpU4SEhNgpKhERERH7abKJnZOTE7169WLJkiWl1ywWC0uWLGHAgAH1+uz4+Hi6dOlCnz596vU5IiIiItXRqIdic3Jy2L9/f+n3xMREEhIS8PPzIyoqiqlTpzJhwgR69+5N3759efXVV8nNzS1dJVtfpkyZwpQpU8jKysLb27tenyUiIiJSVY06sdu4cSNDhw4t/T516lQAJkyYwJw5c7jxxhtJSUlhxowZJCUlERcXx6JFi8otqBARERFpCQxWq9Vq7yCaqrM9dpmZmXh5edk7HBEREWmGqpNvNNk5diIiIiJSlhK7GtDiCREREWmMNBRbCxqKFRERkfpWnXyjUS+eaOzO5sQ6M1ZERETqy9k8oyp9cUrsaiE7OxtAZ8aKiIhIvcvOzr7oNmsaiq0Fi8XCiRMn8PT0xGAwlCnr06cPGzZsqPC+rKwsIiMjOXr0aLMbwq3svZvy8+uq3Zq2U937qlq/Lurp99z0nl8X7damDXv+ni9WR7/npvfslvB7tlqtZGdnExYWhtFY+fII9djVgtFoJCIiosIyk8l00b8UvLy8mt1fHFV576b4/Lpqt6btVPe+qtavy3r6PTed59dFu7Vpw56/56q2pd9z03l2S/k9V/VABK2KrSdTpkyxdwh2Ye/3rq/n11W7NW2nuvdVtX5d12tu7P3ejfn3XJs27Pl7tvf/pvZkz3evz2e35N9zRTQUawdaTSvNiX7P0pzo9yxNnXrs7MDZ2ZmnnnoKZ2dne4ciUmv6PUtzot+zNHXqsRMRERFpJtRjJyIiItJMKLETERERaSaU2ImIiIg0E0rsRERERJoJJXYiIiIizYQSu0bg2muvxdfXlxtuuKH02tGjRxkyZAhdunShe/fuzJs3z44RilRdRb/njIwMevfuTVxcHDExMbz77rt2jFCk6ir6PZ+Vl5dHq1ateOSRR+wQmUjFtN1JI7Bs2TKys7P58MMP+eqrrwA4efIkp06dIi4ujqSkJHr16sXevXtxd3e3c7Qilavo92w2myksLMTNzY3c3FxiYmLYuHEj/v7+do5WpHIV/Z7P+vvf/87+/fuJjIzkxRdftFOEImWpx64RGDJkCJ6enmWuhYaGEhcXB0BISAgBAQGkpaXZITqR6qno92wymXBzcwOgsLAQq9WK/ptSmoKKfs8A+/btY/fu3YwZM8YOUYlcmBK7WlqxYgVjx44lLCwMg8HAggULytWJj48nOjoaFxcX+vXrx/r166v1jE2bNmE2m4mMjKyjqEUqVp+/54yMDGJjY4mIiGDatGkEBATUcfQiZdXn7/mRRx5h5syZdRyxSO0psaul3NxcYmNjiY+Pr7B87ty5TJ06laeeeorNmzcTGxvLqFGjSE5OrlL7aWlpjB8/nnfeeacuwxapUH3+nn18fNi6dSuJiYl89tlnnDp1qq7DFymjvn7P3377LR06dKBDhw71EbZI7VilzgDWb775psy1vn37WqdMmVL63Ww2W8PCwqwzZ84sU2/p0qXW66+/vsy1goIC66BBg6wfffRRvcUsciF1/Xs+3z333GOdN29encYrUpm6/D1Pnz7dGhERYW3VqpXV39/f6uXlZf3Xv/5Vr/GLVJV67OpRUVERmzZtYsSIEaXXjEYjI0aMYM2aNZXea7VamThxIsOGDeP222+v71BFLqo2v+dTp06RnZ0NQGZmJitWrKBjx471Gq9IZWrze545cyZHjx7l0KFDvPjii9x1113MmDGjvkMWqRIHewfQnKWmpmI2mwkODi5zPTg4mN27d5d+HzFiBFu3biU3N5eIiAjmzZuH2Wxm7ty5dO/evXReyMcff0y3bt0a8hVEStXm92wymbj77rtLF03cf//9+i2LXdXm9zxgwICGDlekypTYNQK//PJLhdctFksDRyJSexf6PSckJDRsICJ14EK/57MmTpzYMIGIVJGGYutRQEAAJpOp3CTxU6dOERISYqeoRGpGv2dpTvR7luZKiV09cnJyolevXixZsqT0msViYcmSJerKlyZHv2dpTvR7luZKQ7G1lJOTw/79+0u/JyYmkpCQgJ+fH1FRUUydOpUJEybQu3dv+vbty6uvvkpubi6TJk2yY9QiFdPvWZoT/Z6lRbLzqtwmb+nSpVag3GfChAmldd544w1rVFSU1cnJydq3b1/r2rVr7RewSCX0e5bmRL9naYl0VqyIiIhIM6E5diIiIiLNhBI7ERERkWZCiZ2IiIhIM6HETkRERKSZUGInIiIi0kwosRMRERFpJpTYiYiIiDQTSuxEREREmgkldiIiIiLNhBI7ERERkWZCiZ2ISD0YMmQIBoMBg8FAQkKCvcPhkUce4Zprrrlg+cSJE0vjXbBgQYPFJSJ1S4mdiDRLgwcPxmAw8Nxzz5W5brVa6devHwaDgaeffrpeY7jrrrs4efIkMTExDR7bpEmTePLJJ0u/JyQk0L179wvWf+211zh58mStniki9qfETkSaHavVypYtW2jVqhW///57mbIPP/yQEydOANCzZ896jcPNzY2QkBAcHBwaNDaz2cz333/P1VdfXXpt69atxMbGXvAeb29vQkJCavxMEWkclNiJSLOzb98+srOzmTBhQpnkKTs7m8cff5yJEycC0KtXr2YZ2+rVq3F0dKRPnz4AHDt2jNTUVABGjhyJm5sbHTt2ZN26dTV/ERFplJTYiUizs2nTJtzc3Lj55pvZs2cPRUVFADzzzDP07t2bwMBAQkJCCA0NbZaxLVy4kLFjx2IwGABK5/jFx8fzxBNPsHXrVqKiopg+fXqt30dEGheHi1cREWlaNm/eTPfu3enYsSMuLi7s3r0bV1dX3nrrLTZv3syzzz5bOtTp4OBQZg7cmjVrcHV1vWDb1a1fm9hq6ttvv+WVV14p/Z6QkICfnx9ffvklAQEBAFx99dX897//rdVzRKTxUWInIs3O5s2b6dmzJwaDge7du/P777/z+eefc88999C+fXs2bdrEtddeC4CPj0+1Vq1Wt35tYvujBQsW8N5771FUVMSNN97IHXfcUa7Orl27OHHiBMOHDy+9lpCQwLhx40qTOoDExETatWtX4/cQkcZJQ7Ei0uycTZ4A4uLiePXVV9m4cSP/+Mc/KCgoYPfu3RftFfv3v/9NTEwM3bp149NPP7V7bJ9++ilffvklb731Fp988gk7d+7k2WefLVdv4cKFjBw5EhcXl9JrCQkJ9O/fv0y9hIQE4uLi6uy9RKRxUGInIs3KwYMHycjIKE2OevTowcaNG5k5cyaenp5s3bqVkpKS0sUJGRkZxMXFERcXx5133gnAhg0b+PLLL9m4cSPLly9nxowZpatVK6pfX7Gd75133uHDDz8kMjKSoKAgXnrpJZYtW0Z2dnaZet9++y3jxo0r/Z6dnc3Bgwfp0aNHmXpK7ESaJw3FikizsmnTJpycnErnwU2YMIFrrrkGf39/wNZjFhgYSGRkJFDx0OqqVau4/vrrcXFxwcXFheHDh7NhwwbGjRtXq6HY6sZ21unTp4mKisLR0ZH33nuPlStXMmfOHPr168fevXtLE8Hk5GQ2btzIwoULS+/dunUrJpOJbt26lV47fPgw6enpSuxEmiEldiLSrGzevJmYmBgc/7+9+2VVGArjOP5bmMWFCZYV8R0o+CcNsQu2Bd/BssE48EVoNlsMgsUoNlGrL8CyaNC4mxx4FcOul8nh+8kPZ8/aj7Oz59i2JMm27YezZfv9/mn36tt7K5VK6fDgIAjU6/UkSafT6SEELpdLtdvthzWPx2P6o8bd4XCQ67qqVqsffT8A+bOSJEnybgIA8lIul9MZb3e73U5hGGq73ep6varVammz2cjzvJf1r3S73fQM3SdEUSTXdTUcDiVJ8/lcq9VKs9ksren3+/J9X6PRKPNzLMvSYrF4e/0YgO/FGTsA+KXZbCoIAjUaDXU6HY3H40xz5abTqRzHebphIosoihTHsWq1mur1utbrtSaTyUON7/saDAaZ1g/DUI7j/LlPAPlixw4A/sH5fNbtdpMkVSoVFQqFnDt6L45jXS4XSZLneSoWizl3BCALgh0AAIAh+BQLAABgCIIdAACAIQh2AAAAhiDYAQAAGIJgBwAAYAiCHQAAgCEIdgAAAIYg2AEAABiCYAcAAGAIgh0AAIAhfgC+yAoG2mDbLAAAAABJRU5ErkJggg==", "text/plain": [ "
" ] @@ -146,7 +147,9 @@ "plt.fill_between(xbin, ylow, yhigh, alpha=0.75, label=\"CB2 Main\")\n", "\n", "ylow, yhigh = np.percentile(hmf_cb2_random, [16, 84], axis=0)\n", - "plt.fill_between(xbin, ylow, yhigh, alpha=0.75, label=\"Random\")\n", + "plt.fill_between(xbin, ylow, yhigh, alpha=0.75, label=\"CB2 Random\")\n", + "\n", + "plt.plot(xbin, hmf_mpdpl2, label=\"MDPL2\", color=\"red\", ls=\"--\")\n", "\n", "plt.yscale(\"log\")\n", "plt.xscale(\"log\")\n", @@ -156,7 +159,7 @@ "plt.xlim(xbin.min(), xbin.max())\n", "\n", "plt.tight_layout()\n", - "plt.savefig(\"../../plots/hmf_comparison.png\", dpi=450)\n", + "plt.savefig(\"../../plots/HMF_comparison.png\", dpi=450)\n", "plt.show()" ] }, diff --git a/notebooks/diagnostic/hmf.py b/notebooks/diagnostic/hmf.py index e485eb6..02503a2 100644 --- a/notebooks/diagnostic/hmf.py +++ b/notebooks/diagnostic/hmf.py @@ -17,6 +17,7 @@ import csiborgtools import numpy as np from tqdm import tqdm +from h5py import File def calculate_hmf(simname, bin_edges, halofinder="FOF", max_distance=135): @@ -46,3 +47,20 @@ def calculate_hmf(simname, bin_edges, halofinder="FOF", max_distance=135): hmf[i] = cat.halo_mass_function(bin_edges, volume, "totmass")[1] return hmf + + +def MDPL2_HMF(bin_edges): + """MDPL2 FoF halo mass function.""" + fname = "/mnt/extraspace/rstiskalek/catalogs/MDPL2_FOF_125.hdf5" + with File(fname, "r") as f: + mass = f["mass"][:] + + dx = np.diff(np.log10(bin_edges)) + if not np.all(np.isclose(dx, dx[0])): + raise ValueError("Bin edges must be logarithmically spaced.") + dx = dx[0] + + hmf = csiborgtools.number_counts(mass, bin_edges) + hmf /= 1000**3 * dx + + return hmf