{ "cells": [ { "cell_type": "markdown", "id": "ca7ff999", "metadata": {}, "source": [ "# fO2" ] }, { "cell_type": "markdown", "id": "37d17ca2", "metadata": {}, "source": [ "In this notebook we benchmark MagmaPandas fO2 codes by reproducing calculations and/or figures from literature. When possible, we compare MagmaPandas results to results from codes or data provided with the original publication of the model." ] }, { "cell_type": "markdown", "id": "9df81e31", "metadata": {}, "source": [ "### QFM" ] }, { "cell_type": "markdown", "id": "f8de67cd", "metadata": {}, "source": [ "fO2 at QFM is calculated according to the method explained in van Gerve et al. (2004). This method is similar to the one implemented in the [online fO2 app by Michael Anenburg](https://fo2.rses.anu.edu.au/fo2app/). Here, the Anenburg implementation is compared to the one in MagmaPandas" ] }, { "cell_type": "code", "execution_count": 1, "id": "250b2737", "metadata": {}, "outputs": [], "source": [ "from MagmaPandas.fO2 import QFM\n", "from MagmaPandas.fO2 import IW\n", "import pandas as pd\n", "import numpy as np\n", "\n", "import matplotlib.pyplot as plt\n", "import geoplot as gp\n", "\n", "mm = 1 / 25.4\n" ] }, { "cell_type": "markdown", "id": "2a10ee7f", "metadata": {}, "source": [ "We import data calculated with the Anenburg app at 1 bar and 10 kbar, over a 600 - 1800 K temperature range" ] }, { "cell_type": "code", "execution_count": 2, "id": "a9622710", "metadata": {}, "outputs": [], "source": [ "fO2_anenburg_file = \"./data/fO2/QFM_anenburg.csv\"\n", "\n", "fO2_QFM = pd.read_csv(fO2_anenburg_file)\n", "T_K = fO2_QFM[\"T_K\"]\n", "P_bar = fO2_QFM[\"P_bar\"]" ] }, { "cell_type": "code", "execution_count": 3, "id": "7509cbe2", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/Users/thomas/Dropbox/research/python/packages/MagmaPandas/src/MagmaPandas/EOSs/tools.py:55: RuntimeWarning: invalid value encountered in sqrt\n", " Q2 = np.where(t > tc, 0, np.sqrt((tc - t) / tc0))\n" ] } ], "source": [ "fO2_QFM[\"fO2_magmapandas\"] = QFM.calculate_fO2(T_K=T_K, P_bar=P_bar, logshift=0)" ] }, { "cell_type": "code", "execution_count": 4, "id": "f866acfd", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtAAAAFaCAYAAAA3lGq/AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAqhZJREFUeJzs3Xd4VFX6wPHvnZn03nsjgYROIBB6lS66goKyioCu2Ls/XdzC2nt3bWsHC1YEld577yGEkN57n2Rm7vn9cUMUKRJImExyPs+TBzK3zHtmknfenHvuOYoQQiBJkiRJkiRJ0gXRWTsASZIkSZIkSbIlsoCWJEmSJEmSpGaQBbQkSZIkSZIkNYMsoCVJkiRJkiSpGWQBLUmSJEmSJEnNIAtoSZIkSZIkSWoGWUBLkiRJkiRJUjPIAlqSJEmSJEmSmkEW0JIkSZIkSZLUDLKAliRJkiRJkqRmkAW0JEnSZZKcnMz06dPx9fXF0dGRnj178uabbyKEAGD9+vUoinLG1/LlywEYOXIkiqLw0UcfnXbempoa3NzcUBSF9PT0y90sSZKkDkcW0JIkSZfBkSNHSExMBGD58uUkJyfz8MMP85///Ic777zztH2zs7PJy8tr+ho9enTTttDQUD7//PPT9v/+++/x9PRs9TZIkiRJGkWc6vqQJEmSWs3o0aOpra1l27ZtKIrS9Pgvv/zC5MmT2bRpE2azmVGjRmEymTAYDGecY+TIkfTu3ZsPP/yQo0ePEh4eDsDYsWNJSEjgueeeIy0tjcjIyMvVLEmSpA5J9kBLkiS1suLiYtatW8f9999/WvEMMGnSJGJjY1m8ePEFncvNzY2rrrqKhQsXApCTk8O2bdu47rrrWjxuSZIk6exkAS1JktTKUlNTAYiLizvr9tjYWE6cONH0vaenJ66urri6utK7d+8z9p81a1bTMI6FCxcyZcoU3N3dWyFySZIk6WzOvEYoSZIktagLGSlnZ2fX9P89e/ag1+sBsLe3P2PfsWPHUl5ezq5du/j888958cUXWy5YSZIk6U/JHmhJkqRWFhMTA0BSUtJZtx87dozo6Oim76Ojo4mJiSEmJqZpnPPv6fV6Zs6cyUMPPURxcTHjxo1rncAlSZKks5IFtCRJUivz9fVl1KhRvP7662f0Rv/yyy8cP36cmTNnNuucN998M5s2bWLmzJlNvdWSJEnS5SELaEmSpMvgzTff5NixY8yYMYO9e/eSkZHBZ599xqxZs7j//vtJSEho1vl69epFcXExzz77bCtFLEmSJJ2LLKAlSZIug+7du7Njxw5AG8McGRnJzTffzIMPPsirr756Uef08fHBwcGhJcOUJEmSLoCcB1qSJMkKzGYz06ZNIyUlhU2bNuHj42PtkCRJkqQLJAtoSZIkK6mvr+e1116ja9euXHXVVdYOR5IkSbpAsoCWJEmSJEmSpGaQY6AlSZIkSZIkqRlkAS1JkiRJkiRJzSALaEmSJEmSJElqBllAS5IkSZIkSVIzyAJakiRJkiRJkppBFtCSJEmSJEmS1AyygJYkSZIkSZKkZpAFtCRJkiRJkiQ1gyygpQ6jqqoKV1dXbr311lZ/rvXr15OQkNDqzyNJknSxIiMj8ff3x2QyNT22du1aFEXh4YcftmJkGkVR6NWrF3369KFXr1588803LX7+6urqFj2n1HHIAlrqML766iv69u3Ld999ZzNJ02w2WzsESZLasfDwcH766aem7z/66KM29cf/1q1b2b9/P59++ik333wzxcXF1g5JkgBZQEsdyIcffsijjz7KsGHDWLx4cdPjn3zyCePHj+eGG26gZ8+eJCQkcPLkyabtn3/+OYmJifTt25cRI0Zw+PDhCzrOZDIxZ84c+vXrR0JCAgcOHADO7J0+fPgwkZGRAKSnp+Pr68sTTzzBsGHDePPNNzl69CiJiYn06NGDmTNnMnDgQJYtW3ZG+/Lz8xk1ahT9+vWje/fu3HvvvQghAFiwYAHTp09n0qRJ9OjRg6uuuoqysrI/3SZJUvs2d+5cPvroIwAqKirYvn07EyZMAODQoUMMGzaMvn370q1bN5599tmm43JychgzZgzdu3fnyiuv5Morr+Stt94CYPbs2dx+++2MGTOGiIgI7rvvPtatW8fw4cOJjIzklVdeaTrPI488Qv/+/enTpw8jRowgJSXlrHHGx8fj6upKenr6OY85lT//9a9/0a9fP2JiYvjll1+azvH9998TFxfHoEGDePLJJ087/4033khCQgK9evXiyiuvpLCwEICUlBSGDBlC79696dmzJ//4xz8u9SWX2gshSR3A4cOHRXBwsDCbzeLHH38UgwcPbtr28ccfCw8PD5Geni6EEOLRRx8Vt912mxBCiM2bN4tJkyYJo9EohBBi48aNolevXn963Lp16wQg1q1bJ4QQ4uuvvxbdunVr2tavX7+m5z906JCIiIgQQgiRlpYmALFo0aKm7X379hWff/65EEKI3bt3C51OJ5YuXXpGG+vq6kRVVZUQQgiz2SwmT54svvnmGyGEEP/+979FYGCgyM/PF0IIcccdd4g77rjjT7dJktR+RUREiEOHDom4uDiRnZ0t3nnnHfHYY4+Jf//73+Khhx4SlZWVTbmvtrZW9OnTR+zatUsIIcTUqVPFk08+KYQQIiMjQ7i5uYk333xTCCHEzTffLIYMGSKMRqOoqakRfn5+Ys6cOcJisYjs7Gzh4uLSlKuKioqa4vnyyy/F5MmTm74HmvZbtWqVcHd3F+Xl5ec85lT+/PHHH4UQQvz666+iS5cuQgghCgoKhLe3tzh27JgQQojnn3/+tPP//pzPPvusuOuuu4QQQtx7773i6aefbtpWUlJyCa+41J4YrFq9S9Jl8uGHHzJr1iz0ej2TJ0/m9ttvJykpia5duwIwdOhQIiIiABg0aBBvvvkmAEuWLOHAgQMkJiY2nauoqIiGhobzHgcQExPDyJEjAZg+fTq33XYbubm5fxqro6MjN9xwAwCVlZUcPnyYmTNnAtCvXz969ep11uNUVeXRRx9l8+bNCCEoLCykT58+XHvttQBceeWVBAQEAHDbbbcxffr0pmPPt02SpPbtpptu4tNPP+XHH39k0aJFLFq0CIC6ujruvPNO9u/fj06nIysri/3795OQkMC6det44403AG0YyJgxY04751/+8hccHBwAiI2NZdKkSeh0OkJCQvDy8iI7O5u4uDhWrlzJm2++SVVVFaqqUllZedp5Bg8ejE6nw8vLiyVLluDh4cEXX3xxzmNcXFy4+uqrAS0np6amArB9+3b69u1LbGwsoOW5Rx99tOm4RYsW8fnnn1NfX09dXR2BgYEADB8+nEceeYSamhpGjBjBFVdc0WKvu2TbZAEttXsmk4mFCxdiZ2fHl19+CUBtbS0fffQRL774IqAVrafo9fqmscdCCObOncsTTzxx1nOf67hzURQFg8GAxWJpesxoNJ62j4uLC4qiND2/oihN35/PK6+8QklJCTt27MDR0ZEHH3zwjHP/MZaL2SZJUvsye/Zs+vbtS5cuXejcuXPT4/PnzycgIIB9+/ZhMBiYOnXqaTnlfHnij7nxbLkyMzOTe++9l507d9KpUycOHjzI6NGjTzvP1q1bcXV1bfr+z4754/OcyrWicTjb2WzevJm33nqLrVu34ufnx08//dSU86dNm8bgwYNZtWoVb731Fq+99tppw0KkjkuOgZbavSVLltCpUydycnJIT08nPT2dLVu28Nlnn5129/nZTJkyhc8++4ysrCxA6+XdvXv3BT3viRMn2LhxIwDffvstISEhBAUFERUVRVpaGiUlJYA2xvpcPDw86NatW1Phv2/fPg4dOnTWfcvKyggMDMTR0ZGCgoIz7lj/+eefm8b1ffjhh6f1pJxvmyRJ7VtwcDDPPvsszz///GmPl5WVERoaisFgIDk5mVWrVjVtGzlyJJ988gkAWVlZrF27ttnPW1FRgb29PYGBgQghmsZQt/QxoPVG79u3j+PHjwPwv//9r2lbWVkZ7u7ueHt709DQwHvvvde0LSUlBX9/f2bNmsULL7zA9u3bm9lKqb2SPdBSu/fhhx/y17/+9bTHevToQXBwMEuXLj3vscOHD+eZZ57h6quvxmKxYDKZmDx58gXdpd6nTx+++uorHnzwQYQQfPHFFwCEhITw8MMPk5CQQGRkJMOHDz/veT777DPmzJnDyy+/THx8PL1798bDw+OM/e69916uu+46+vTpQ0hIyBlF8JgxY7jllltIS0ujU6dOfPrppxe0TZKk9m/OnDlnPPaPf/yDm266iUWLFhEZGXlaT+/rr7/OrFmz+Prrr+nSpQtDhgw5a146n549e3LdddfRvXt3wsPDGTt2bKscA+Dv78/777/PlClT8PHxaRraBjBx4kQWLlxIXFwcoaGhDB48mBUrVgDwzTffsGjRIuzt7RFC8O677zarjVL7pYjzXdeQJMnqampqcHZ2RlEUjh49ysiRI0lOTsbLy+uCz7FgwQKqq6t56aWXmrVNkiTpbOrq6rCzs8NgMJCXl0f//v1Zs2ZN0xhjSWrvZA+0JLVxW7Zs4ZFHHmkaw/fBBx80q3iWJElqaSkpKcyaNQshBCaTiX//+9+yeJY6FNkDLUmSJEmSJEnNIG8ilCRJkiRJkqRmkAW0JEmSJEmSJDWDLKAlSZIkSZIkqRlkAS1JkiRJkiRJzdDhZ+FQVZXc3Fzc3Nzk6muSJHUIQgiqqqoIDg5Gpzt7P4rMjZIkdTQXkhtP6fAFdG5uLmFhYdYOQ5Ik6bLLysoiNDT0rNtkbpQkqaM6X248pcMX0G5uboD2Yrm7u1/wcUII6urqcHJyssneGVuPH2y/DTJ+67P1Nlxs/JWVlYSFhTXlv7ORudE24wfbb4OM3/psvQ2tmRtP6fAF9KkX1t3dvdkfEnZ2djb9w2XL8YPtt0HGb3223oZLjf98x8jcaJvxg+23QcZvfbbehtbMjafImwglSZIkSZIkqRk6fA/02ZhMJiwWy3n3EUJQX1+Poig2+9fZhcSv1+uxs7O7jJFJkiRJkiS1bbKA/gOTyURqaioXssK5qqp/epdmW3Yh8SuKQnR0tCyiJUmSJEmSGskC+g8sFgtCCIKDg3FwcDjnfkKIpgLUVnug/yz++vp6cnNzsVgssoCWJEmSJElqJAvoc3BwcMDR0fGc2ztCAS1JkiRJkiSdyXbHH7QTVVVVJCYm4urqyuHDh8/YbjabmT17NsOGDeO+++47Y/uCBQtYtmzZ5QhVkiRJkiRJQhbQVufk5MSyZcu49tprz7p96dKlhIaGsmnTJmpra9m6deslP6eqqpd8DkmSJEmSpI5KFtBWZjAY8PPzO+f2bdu2MW7cOAAmTJhw1gL622+/ZcKECYwfP57KykpUVWXcuHGMGDGCsWPHUllZCUC3bt2YNWsWjzzySOs0RpIkSZIkqS0QrdtZKMdAn4fJIig1ms++UQhUoaJTdHABY4i9HQ3Y6Zs/1ri8vLxpEQMPDw9KS0vP2Cc4OJhPPvmEDz74gA8++ICHHnqIJUuW4OTkxBtvvMHXX3/N3/72N7Kzs9myZQteXl4XNMuIJEmSJElSWyXqSqH0OKIsFcpOIBq/KEtF3/MWGP1Mqz23LKDPo9Ro5oujxS1yrpndfAlwaf5MFl5eXk09yOXl5Xh7e5+xT79+/QAYMGAA7777LjU1NcybN4/MzEzKy8uZNm0aADExMXh5eV1CKyRJkiRJki4vUVcGRYcRxUcQxUehqPHfmvzfdnIJAK9oFO8u0GkiatCQVo1JFtDn4e1oYGY337NvvIge6IsxcOBAVq5cyfDhw1mxYgVz5849Y599+/Yxbdo0du/eTUxMDMuXLyc4OJiFCxfyxhtvNPVa2/Kc1ZIkSZIktW9CCKjIQBTsRxTsh8LGfysytB0UPXh3RvHtjtLnbyh+3bWC2SsGxcHttPOIurpWjbXdFtApKSncfPPNFBcX4+npySeffEK3bt2adQ47vXLOXuOWnAZu0qRJ7N+/n+TkZObNm8fs2bOZN28e7733HlOmTOHHH39k2LBhxMfHM2jQoDOOz8rKYvz48QB88803VFVV8fTTTzNp0iSCgoIICwu7pPgkSZIkSZJaklDNUJKMKNiHyN8HBfsQBfvAWK7t4OSLEhiP0nU6SkAfFL9e4NMFRW9v1bhPabcF9Lx587jtttuYPXs23377Lbfccgvbtm2zdlhn9csvv5zx2HvvvQdoNxl++umn5zx2wYIFZzzm7u7O3r17z3h89+7dFx+kJEmSJEnSRRDmeig+gsjf2/RF4UEwN/YSe0ahBMSjG/AQBPRBCYwH1+A2vU5FuyygCwsL2bt3LytXrgRg2rRp3H333aSnpxMZGWnd4CRJkiRJktopYTZC4UFE3h5E/h6tWC46DKoJFB34xKEE9kXpOh0C+2q9y46e1g672dplAZ2VlUVwcDAGg9Y8RVEIDw8nMzPznAV0RUUFQggaGhpQVRVVVS94pgpbn9HiXPGfelwI0Sbb+Pv4bJGM3/psrQ1CtUBZqtaTU3wEUXQEXcgIRMK85p2nGe09lRtPcXBwwMHB4U/PbSuv6R/Zevxg+22Q8VvfhbahqVjO3wN5u7ViufgIqGbQGcC3O0pgX+h9C0pAPPj3QrF3OefzXe74LyWOdllAA2d0+//ZixIeHg5AREQE7777Lq6uroSGhv7p89j6oiTni//UHxJGo7FNJwKj0WjtEC6JjN/62lwbhAqVWSglR1FKjmj/Fh9FKTuOYqnXdnH0Rvh0hdDmx1/XjJtrTuXGU+bPn8/jjz/+p8e1ude0mWw9frD9Nsj4re+0NlhMKCWHUQr2oRTs1b5KjqCoZoTOgPDphgjoi+gxF+HfB+HbAwyOp5/QArTyzX3njP8CNCc3tssCOiwsjOzsbMxmMwaDASEEWVlZZ3wQ/F5mZibu7u40NDRQWFiIn5/fn85aceomQlv1Z/HrdDp0Oh2Ojo44Ojqecz9rEUJgNBpxdHRs0+OkzkXGb31toQ2iphCKDiGKDjf+e0TrwWmo1nZw8AC/HiihA6HPLSh+3cG3uzZlE1xU/CaT6YL3PZUbT7mQHmhrv6aXwtbjB9tvg4zf+lTVQn3eQRxKDkL+bkTebijYD5Z6bSYMv+4oQf2g720ogQng3xPlj8WyFV3se9Cc3NguC2h/f3/i4+NZuHAhs2fP5rvvviMyMvK84589PDxwd3fHaDRSXFz8p7Nr/L5H1hZ/QS4k/lOPK4rSptvY1uP7MzJ+67scbRCmOig+iig8iCg8qBXLhYegtlDbweAIvt1Q/HqgxF2rFc1+PcAt5Jyxnfo9bm78zdn3VG5sLlv/ubD1+MH22yDjv3xEVS4idycibyfk7kTk7cK+vhIB4N0FJag/SrcbUIL7azf52TlbO+QL0pq5sV0W0KDNYjF79myeeeYZ3N3dzzuThbU9/PDD7Nixg/DwcD7++GPs7X+bomX9+vXcdNNNREdHo9frWbNmzRnHv//++yxcuBBFUXBycuLFF1+kZ8+edO7cmZCQEKqqqliwYAFTpky5nM2SpA5JCAFV2YiCA1B44LeCufR449KyCnh1QvHridL3du1f/57aPKY6vbXDlySpnRMNNdp45dwdTV9U5WgbXYNQghNRBj5Kg08vHCKGoHOSC7CdTbstoGNjY9vstHW/t2/fPvLz89m0aRNPP/003377LTNnzjxtnxkzZvDSSy+d9fjVq1ezYcMG1q5di8FgoLi4mMJCrUfLw8OD9evXk5+fz+jRo2UBLUktTFgatF7lU5P+FxxAFB4AY5m2g6Mn+PVCiRqLMuBBFP9e2qVPe1drhi1JUgchhAqlKYicbYic7Yic7dqMGMICdi4QlIDS/a9a0Rw8AMU9tPE4bSESxdHJyi1ou9ptAW0rtm3bxrhx4wCYMGECH3/88RkF9HfffceOHTu49tprue+++07b9sUXX/DYY481zTji6+uLr+/pqycGBgZiNBo5dOgQ9957L/X19cTHx/PGG2+0YsskqX0RdWWnrYwlCvZD8VHtbnMU8I5B8e+NbsCDENAbxb83uIfZzCVcSZJsn6iv1HqUs7dpRXPujsY/6BVtiFjIQJR+d6GEJGozZMirXhdNFtDnIcy1UH7s7NsQoAqETgEu4APSMw7FcOaYofLycoKDgwGtx/jUstunJCQkkJycDMDVV1/N0KFD6devX9P2vLy8puOfe+45fvzxR6666irmz5/ftE9ycjJeXl7ExMSwdu1aFEVh6tSppKSkEBsb++exS1IHI6rztOmY8vci8htXx6pI1zYanLQbZoIHosTfjhLQR/te9ipLknQZCSGgLBWRsxWRrX1RdBgQ4OiFEjIQ3YAHIGQgStAAFEcPa4fcrsgC+nzKjyGW9j/n5lNl84VM8KZM2QW+fc943MvLi8rKSu3pysvx9vY+bbur628fyldddRUHDhw4rYAODg4mJycHHx8fHnvsMQYOHMiyZcsAbf7WkSNHotPpePvtt0lPT+fBBx+ktraWtLQ0cnNzZQEtdWhCCKjMQmQdRW0qmPdCTb62g6MXSkA8Stw07d+AeG0pWZ1MnZIkXV7C0qDlqKzNiOwtiOwtUFukbfTthhI6GGXAAyghg8AnVl79amXyU+B8POO0wvcsBAKhChSdgnKBPdBnM3DgQF5++WVmzZrFihUrGDJkyGnbKysrm+6A37RpE7fffvtp22fOnMnzzz/Pp59+isFgwGw2N207NQb6lHvuuYd77rmHSZMmMXXq1DY9t7MktbSmYjl/DyJvT+PUTHuwrytGBXAJQAnsh9LnVm2VrIB48IiQH0KSJFmFMFZovctZmxBZWyBvJ5iNYHDSxiz3uU0rmkMGocgb/S47WUCfh2JwPmuvMQBCIFQV5U+mu/sz8fHxBAYGMmzYMMLDw3nkkUcAmDdvHu+99x6LFy/m/fffx2AwMGTIEIYPH37a8WPGjOHEiROMGjWqaX7WBQsWnPW5pkyZwgMPPMD//ve/0wptSWqPRHUBIm8XIm+Xdsd53q7femtcAlGCElD63UWDd08cIgahcwu2bsCSJHVoojoPkbmpsWDeBIUHAaH9cR86FGXkMyihQ7Vp5PR21g63w1NEB++GrKysxMPDg4qKiqZ5oNPS0oiKijrv4iFCiKaFSGyxh+pC4r/Q18JahBDU1dXh5ORks++BjL+FYqmv1HqV83Y2zmW6CyqztI1OvtrcpYEJKMH9tV7mxmK5LbXhYlxs/H/Mexe7T0vG1FbYevxg+23oKPGLigxE5gZE5kZE5gYoO6Ft8IpBCRuGEjYUJWyY9v1lfh06ynvwR83Je7IHWpIkmyJUMxQe+t0cpjuhOAkQYO+KEpigTfgf1B8lKEEOw5AkyeqEEFCRjshYj8hcj8hYD5WZ2ka/niidxqOEPY0SPgzFNciqsUoXRhbQkiS1aaIqV5u/NLdxDtP8PWCq1ZaTDeiNEjYcJfFhlOAB4BMnp2WSJKlNEOXpWrF8WsGsQGDjjcnhI7ShGc4+1g5VugiygJYkqc0QlgYo2I84NYdpztbfhmK4h2k3zgx/AiU4EQL72sxyspIktX+iKheRsQ6Rvga79LWolRn8VjBfixIxUhua4ehp7VClFiALaEmSrEbUFmvzl+Zs1e4yz9+t3WWud4CgfihdZ2h3mIcMbBq3LEmS1BaIujKtdzl9DSJ9DZQ0rhvh1xO102Tsosdqvcxyhox2SRbQkiRdFtqk/yca5zDdrBXMpdoiQbgGoYQOQRn5DErIYO0uc4ODVeOVJEn6PWGu1/7gT1+NSF8NebtBqOAVjRIxBmXYf1AiRoKzHw11ddjb6A140oWRBbQkSa1CqBYoPKBNy5S9WZuWqaYAULSV+yJHowz9J0roEHmjnyRJbY4QAoqTEGkrta+M9WCuA2c/lMgx2jzMkWNQPCPPPE5q92QBbWVVVVVcccUVHDlyhO3bt9OjR4/TtpvNZm699VZSU1Pp27cvr7/++hnneP/991m4cCGKouDk5MSLL75Iz5496dy5MyEhIVRVVbFgwQKmTJlyuZoldUDC0gB5e7RpmbI2IrI2Q0OVNhwjeABK77na+L+QQXIMoCRJbZKoK0WkrUKcXIFIWwlVOaB3QAkbhm74f1CixoJ/LxRFZ+1QJSuTBbSVOTk5sWzZsqYFVP5o6dKlhIaG8sknn/C3v/2NrVu3Mnjw4Kbtq1evZsOGDaxduxaDwUBxcTGFhYXAbysR5ufnM3r0aFlASy1KmOtRcrai5m+FrI2I7K1a74y9K0roEHSD/67NYRqUgGJoe/OIS5IkCdWiLbR0cjnqyeWQu1MbluHbXbsHo9M47Q9/ecOy9AeygD6POrOF1PKas24T8NtKhBdwrmhPF5wMZ06vZTAY8PPzO+dx27Zt48orrwRgwoQJZxTQX3zxBY899hgGg/ZW+vr64uvre9o5AgMDMRqNHDp0iHvvvZf6+nri4+N54403LiBySdIISwPk7mqckmkdInsrduY6hIMHSvhwdCOeRAkbrt1xrpOpRZKktknUliDSViBO/Iw4uQLqSsDBAyVqLMrEW7U5md1DrR2m1MbJT7nzSC2v4aqlu1rkXD9N6U8P3wtfzeuU8vLyptVwPDw8KC0tPW17Xl4ewcHa7ATPPfccP/74I1dddRXz589v2ic5ORkvLy9iYmJYu3YtiqIwdepUUlJSiI2NvYRWSe1Z0xjm9LXaXeZZm8BUAw7u2tzLI56kIWAgDmGJ6PQylUiS1DYJIbRcduJn1NRfIGe71sscEI8Sfzu66IkQkij/8JeaRf60nEe0pws/Tel/1m0X0wN9Mby8vKisrAS0Ytrb2/u07cHBweTk5ODj48Njjz3GwIEDWbZsGQAVFRWMHDkSnU7H22+/TXp6Og8++CC1tbWkpaWRm5srC2ipiTZLRqo2/i99NSJjHRjLwOCkjf8b+i+UiFHa/Ms6PUIIRF2dXLhEkqQ2R5iN2h//KUsRJ5ZBVbY2vCxyLMrE91GiJ8qpMaVLIgvo83Ay6M/ZayyEQFVVdDpdq84eMHDgQFauXMnw4cNZsWIFc+fOPW37zJkzef755/n0008xGAyYzeambafGQJ9yzz33cM899zBp0iSmTp0q7xSWtEuZ6WsQ6asQaaugIgN0BghOREm4B13kGO3/cko5SZLaOFFTiDixDJHyk5bPTLXgGYUSOxWl8xSU8OEoentrhym1E7KAbgMmTZrE/v37SU5OZt68ecyePZt58+bx3nvvMWXKFH788UeGDRtGfHw8gwYNOu3YMWPGcOLECUaNGoWDgwMODg4sWLDgrM8zZcoUHnjgAf73v/+dVmhLHYewmCBnO+rJFYi0FZC3BxDg0xWl81XaGMDwESgOzR9uJEmSdLmJ0hOI4z+iHl8C2Vu0B0MHa1fMYqaAb1c5RabUKhTRwbshKysr8fDwoKKiAnd3d4xGI2lpaURFReHoeO6ZAy5XD3RruZD4L/S1sBYhBHV1dTjZ6GT1lyt+UZGBSF2OOLkckbEW6ivByUcrlqPGaf9exA0ztv76g+234WLj/2Peu9h9WjKmtsLW4wfbb8P54hdCQME+1GPfI47/AMVHweCo5bIuV6PEXIni4m+lyH+L0ZZff7D9NrRmbjxF9kBLUjsjLA2IrC2I1F8Qqb9C8RFQ9BAyCF3iIyidxjeNY5YkSWrrhFC1K2fHvkMkfw8V6eDopQ3LGPGUVjzbX9x9RpJ0sWQBLUntgKgu0ArmE8u0sX8NVeASqN0oM2wBStQVcvESSZJsh2pBZG5APfatVjRX54FLAEqXa1DipmlDzfR21o5S6sBkAS1JNki7jLkfcWIZ6oll2uT/KBAyEN2gR1GiJ0FAb7laliRJNkMIFbK2oCZ9jV3St6i1BeAWitJ1Brq4aRAySF45k9oMWUBLko0QlgZExnrE8SWIlJ8ap2Vy0yb973cXSqeJKC7nXpRHkiSprRFCQO4O1CNfIo59o/U0u4Wixl6HXY8bUEIHyo4AqU1qlwX07NmzWb16ddOKfGPHjuXFF1+0clSS1HzCWIFI/RWRsgSR+ot2A6B7OErsNdqsGXJapg6rwaJytKQKvWqih5OTtcORpGYRhYdQj36JOPoVlKeBazBK1+nouk5HBCfSYKzH3kZvYJOsq8GicqSkEjvVQvdWzI3tsoAGeOyxx7j77rutHcafqqqq4oorruDIkSNs376dHj16APD111/z2muv4eTkxKeffkpYWNhpx7366qt8++23eHl5sWjRIjw8PE7bvmXLFv71r39hsVhQVZX77ruPadOmMXv2bI4cOYJOp2PcuHE88cQTl62t0oURNUVawZz8PSJtNagmCIhHN+BBlC5Xg39v+aHSARXV1rO3qIJ9hRXsLargUHEV9RaVGdEBPBPgY+3wJOlPiYoMxJEvUA8v0m5udvRCibsWpftMlLBhvw3P6NiTg0nNVFhbz97CCvYVVbC3sIJDJVU0WFT+2jmQJwK8//wEF6ndFtC2wsnJiWXLlvHII480PWYymXjllVfYtGkTu3bt4sknn+T9999v2l5UVMTSpUvZvHkzX3zxBW+//fZpS3eXlJTwyCOPsHTpUnx8fDCZTOza9duS5B9//DFdu3YlISGB22677YziXLr8RFUO4tj3WtGctVF7MGwYujEvoXT5C4pHuHUDlC4riyo4Xl7N3sIK9hRqHwqZVXUABDk7EO/vwSO93eijP0m0e7l1g5Wk8xDGcsSxb1EPL4TMDWDnjNL5apRRz6F0GievoEnNcio37imsYE9BBXsLy8mqNgIQ5OJAXz8PJoa6Ea9PpZNHRavG0m4L6FdeeYX333+f8PBwnnrqKfr06XPe/SsqKhBC0NDQgKqq2ldDDZQcO+cxQlUROh0X9LeyTxyKnfMZD+v1+qahJkIIhBAcP36cbt26YWdnx+DBg3n44YdPWzVw586djBgxAoDx48cze/bs07YvW7aM66+/Hm9vb4QQGAwGBg0a1LSPEAKdTkdsbCw5OTnccsst1NfXY29vz7fffts09+Hv92+L04X/Pj5bpFZmozv8FZbUxgUAdHYQeQXKhHe14Rm/m8u0LbbR1l9/aDttqDVZOFBcye7CcvYUaD0p1SYLBkWhm48ro0N9iHerI14cJrhiMxRuhfQjAJg734EIHdys52tOe0/lxlNOLdj0Z+e29mt6sWw9frB+G4RqhtTlqIc/g5Sl2lW0iDEoV36qdQg4uJ0R62nH2/h7YOvxQ9tpQ63Jwv5irVjeU3hmbhwT5ktftzri1UMEVWyCom2QfhQAc5e7ESEDm/V8zWmvTRbQw4YNIykp6azb9u3bx9NPP01QUBA6nY4ffviBiRMnkpKSgqur6znPGR6u9fBFRETw7rvv4urqSoihED4dcN5Y1AsN+uadKIF9z7n51MImqqpSWlqKu7s7qqqd3Ww2N/0foKysDDc3N1RVxc3NjZKSktO25+bmEhMTg6qqLF++nBdeeAF3d3d+/PHHpuepqanh0KFDREZG8v333+Pk5MSbb77JV199xa233qq1rTEeo9Fo9V+i8zEajdYO4cJV56FL+QFdyvfocreh19ujRlyBOv4D1KhJ8Pup5urqrBZmc9jU638Ol7sNJUYTe4sr2V9Uxd7iKpLKajALgbudnt6+bsyNDSTePp9ept24lm1Bd3IbSn0RAgXVoyuqTyK10feQ59QPe49ogpoZf10zfrZO5cZT5s+fz+OPP/6nx9n6z4Wtxw+Xvw1K8RF0RxeiO/YVSm0hqm9P1MH/Ro2dDq5B2k4qF5zbbP09sPX44fK3odjYwL6iKvYVV7G3uJKkshosAtzt9fTxceOW2ADiHQro1bATl9KtjbmxuDE3dkP1GUht9H3kOvfD0aMTga2YG22ygN60adMF73vNNdfw2GOPkZycTL9+/c65X2ZmJu7u7jQ0NFBYWIifnx86gy/M2X3OY06t5HdBfOJQzrOvoijodDp0Oh0+Pj5UVlY2ndtgMJz2PN7e3qSmpqLT6aisrMTHx+e07SEhIeTm5qLT6Zg0aRKTJk2if//+TasO3nrrrTg5OXHffffh6urK7bffTmZmJuXl5UydOrXpXKficXR0bLMrERqNRhwdHdv0mGBRW4JI/k67WSZzA+gM0Gk8ypWf0BA2DkePgDYd/7nYyut/PpejDUIIcqqN7CosZ2d+ObsLKjhZWQtAiIsjCQEeXBftS4IhlZiaDegKt8CxnWCuBb0j+A5AdLmVCs9EMh3jyTa6kFvdQHW9CvXQW28mytetWfGbTKYL3vdUbjzlQnqgbfnnwtbjh8vbBlFXhjj6BeLgp5C/G5x8tTHNvWajD+hzcee08ffA1uOHy5cbs6uN7Coob/yqIK0xN4a6OpLg78GMGF/66VOJqV6PrmgLJO0ESx3oncBvACL2Nso9Esly6ENWvZYbaxpzY7zeQqRP6+VGmyyg/0x2djahodrSxNu3b6ekpISYmJjzHuPh4dG0lHdxcbFWPNo7QtDZi24hBIqqorTgUt6KoqAoCp07dyYpKalp7HKvXr1Oe47+/fvz0ksvoSgKK1euZMiQIadtnzx5MlOmTOGvf/0rvr6+mM3mpvODNga6e/fuqKrKDz/8QHBwMAsXLuSNN96gtLS0ab/f/9uWk0BbjE80VGvTzR39EnFyBQgVJXIMyuQPtRk0HD21Xv26ujYZf3PYevzQsm0QQpBWWcuOfK1g3llQRl5NPQBdPF0YFOzFPd19SVCOEFTxCxRshrR9ICzg4AMBQ1D7LKDELZFMQzeyayG3uoGGMoFOgQAXC7E+ToS42hPkYgfm+mbH35x9T+XG5rL1nwtbjx9arw1CCMjciLr/A8Sxb0E1o8RMRhnyvfZvC41rtvX3wNbjh5bPjakVtezML2NHQTm788vJq9VyY6yXK0OCvbivhy8J4jCBFT835sb9jbnRV8uN8U9Q4tafDEM3smsU8hpzo74xN3b1cSLYrTE3mlo3N7bLAnr27NkUFBSg1+txcnLim2++OWOWirZk0qRJ7N+/n+TkZObNm8fs2bO5//77GTFiBI6Ojnz22WcAPPfcc8yYMYOoqCimTJnCkCFDmmbh+D0fHx9eeOEFpk+fDmg9yffff/9Zn3vgwIE888wzTJo0iaCgIHlD4SUQqhmRvgZxeKG2cpapFkIHo7viVZS461BcA6wdotQKVCFIKa9hR14ZOwvK2VlQTnFdA3pFobuPG5MjA+jvJein7sez9HvI3wTHtPHLuIRDwFAsMXMocB1EBpFkV5soqDFhMYK9voFgV3v6B7oS7GZHoIs9Bt1vCV4IQZ3ZSg2XOhxRXYA49Cnqgf9BaQp4xaAb/h+UnjejuAZaOzypjVGF4HhZDTvyG3NjfhklRhN6RaGHjxtXdgqgv6egn7oPj5LvGnOjNn4Z1wgIGIq58y0UuiaSThQ5p+VGk5Ybg1wJdrMn0MXuzNx44Z3JF0URbXlw62VQWVmJh4cHFRUVTT3QaWlpREVFnXfYwqmxxLoW7IG+nC4k/gt9LaxFCEFdXR1OVpwrVAgBhQdQD32OOPIF1OSDdyy6njehdP8rimfkeY+1dvyXwtbjh4trg2gsmLfnl7E9r4yd+eWU1puw0yn08nVnQKAnAzwt9DXvxaV4PeRvhMrj2sHuXSBwGGa/YeS5JJJhCSS7uoHCGhMCcDboCHGz175c7fF1NqA7T1wX+x78Me9d7D4tGVNbYevxQ8u2QQihLeC09x3E8R9A0WtLaff5m7acdiv1cNvye2Dr8cPFteH3nQnb8rXcWPa73JgY6KXlRstunIvWQ94GqDqhHewRCwFabsx1GUiGxZ/sqgaKas+SG93s8XWyXm48pV32QEtSaxM1RYgjC1EPfAxFh8DZD6XbDeh63gSB/Ww2aUpn0oZk1LEtr5TteWVsb+xFOfWhcENsCIleKvHmXTgVfQ3pG6AyRTvYsxsEj8bc+9/kuiSSbvIhu6qeogozVICrXR2hbvb09HMmxNUeL0e9/NmR2gRhLEcc/BR137vabFQ+cehGv4TS8yYUp9abW1eyHaeGq23N0zoTtueVNXUm9PZz569xISR6Wog378Kx6CtI+13B7NkdQsZi9nuCHJdEMkxeZFc2NOVGN/s6Qtwc6OXvTKibPZ4ObS83ygJaki6QsJi0VQEPfow4sQxQtOnmRj6jLaett7N2iFILyak2sjWvlG2NHwr5tfXoFYXevu7M6BLCQB+FePMenIq+g4x1cDBZO9CzGwRfgbmP9qGQ2eBFVlU9RWVmKAM3eyNhbvbEB7gQ6maPu33b+1CQOjZRsB9191vaFTXVhNLlGpTx/0WJGCl/ViVyquvYklum5cb8Mgpq6zEoWsF8Q2wIid4KfS27cSz8FtLW/Xb1zbM7hIzD5P8UuS6JZDR4kl35W250t68n1N2e+EAtN3o4tP3ytO1HKElWJkpTUPf/D3HoU6gp0FYFHPOydqe5s6+1w5NaQHm9ie355WzNK2VLbinplXUoQDdvN66MCmCgvyP9xQFcCpdB1jo4eFA70L0zBI1C7fUvct0Gk9HgSVZVA4WlJkSp1osS6ubQVDDbwoeC1PEI1YxI/hF19xuQtQncQtENno/S5xaUU9PPSR1SmdHE9oIytuSWsTW3lIwqLTd293Hjqk4BDPR1IEEcwKVoaWNuPKQd6N4FgkZh6fMf8lwHkt7gSXbl73NjY2eCDRXMf2R7EV8m9fX1593eEcZA/9lr0J4JsxGR/ANi3/uIzPXakrM9bkTXey7KRU7NJLUd9RaVPQXlbM4tZVN2MUfLahBApLsTg4O8eSTenURDCp7FP0PeWji6E4QZXMIgaBSi+4MUug8h3exHVmUDeeUNqGXgbFdHeOOQjDA3ezza4GVHSTpF1BYj9v8Pde9/oTJLW/30mm9QYv+CopPlQUdUb7awu7CCzTklbMopIakxN0a5OzM0xJtHA9xINKTgUbwMctfCkZ3aLBku4Y258WEKPIaQYfIls7Ke/DITaim42NUR5mbfNCSjPeRG+RvyB3q99qbm5ub+6b7Nmge6DbqQ+BVFQa/XX6aIrE+UJKPufRdx+DOoK4XwEeiuWqjdNGNoezdSShdGCEFyWQ2bc7Ue5h35ZRgtKj6Odgz09+CmrmEMcS0huHwdIncVbNkI5hpw8IbAkZD4GuXeI0lTQ8msMpFT2YCpXOCgryHMzZ4RYe6Eudvj7Wiw+Q8Fqf0TJcmoO1/VrqoJgdJ9JrqEe1AC460dmnSZCSE4VlbNltxSNueWsjO/HKNFxdfRnsQAd27uFsZgl2ItN+asghMbtTnqHbwhaBREv0GZ90jS1RAyq0xkVzZgPpUb3R0YGe5EmJtDu7y/QxbQf2BnZ0d0dDQWi+W8+9n6ROkXGr9er8fOrn2P7RWqGXH8J8Te/yLS12gLAfS+BV3vW1F8ulg7POkilRlNbM4tYWOO9sFQUFuPg15H/wBPHojvxBAf6FK7BZH1Pvr9G6A2F6F3AP+hKL0fp9ZvNOmGrmRWmsisbKA2U0WvVBPiZk9isCth7vb4O9ud905wSWorhBCQtQl1x8uIlJ/AJQDd4MdR4uehuPhZOzzpMio1NmhX33JK2ZRTQmFdAw56HYmBnjzYtzE3Vm9CzVqOfu96qMvTcmPAUJQ+/6LWbxTpujgyq87MjQODXQl3d8DvT2YQag9kAX0WdnZ2f1o0CiEQQth0AW3L8bcEUZWL2P8B6r73oToXQoegu3oRSuw0FMO5V1qT2iaLKthfVMGGHK1oPlRciUBbvGRKVADDgtzprxzFPv8bSFsFu/cCoHr0gKjrsQRdQa5zf9JrDKRX1FOabQaq8Hc20M3XiQgPB4JdT5+HWZLaOqGaEce+Q93xEuTtBt/u6CZ/qN3DIa+qdQhmVWV/USUbskvYmFvC4eIqBNriJVdHBzIs0J0E5Qj2+YshdRXs2geA4tkTomdiCbyCHOf+pNfoyKhooDTrVG60o7uvE+EdNDfKAlrqUIQQkLMdddfriOTvQO+A0uMmdP3uQPHvZe3wpGYqqqtnU04p67O18XoVDWY87A0MC/Hmr3EhDPWsI7BkFSJ7OSSvA3M1OPpB8BXQ9R5KvUeQUuNObq1KTkEDFlGPq52JSA+Hpp4UJzvbHaYldWBmI+q+zxA7XoKyVJTIK1Bm/KrNGNRBO006kqLaejbklLAhu4RNuaVUNpjxdDAwNNiHm+JCGepRQ0Dpai03HlunDVlz9IeQsdD9Pko9R3C81o3cWpXcggYswoirnU7LjSGuhLvJ3CgLaKlDEJYGRNI3qLteh7xd2gpaY15G6TUbxaH5yxRL1qEKwaHiStZml7A+u5hDxVUA9PR1Y1bXUIYHu9PbvB9d7idwaDlUHEMoevAfjNLrMRoCx5Jh6EpGpYn0inqqy1UMipFQN3uGhroT4SHHMUu2TdRXoe59B7udryJqClDirkV3zWKUwL7WDk1qRaoQHCyuZF1WMeuySzhcUoUC9PTVxjGPCHSjl2UfupyPGnNjcmNuHILSe76WG/VxpDfmxpqyxtzo7sDQUHciPdrnOOZLIQtoqV0TtcWIve9qd5lX56FEjUWZvgwleiKK0rH/erYVlQ1mNueUsK6xaC4xmnC3NzA8xIfZXcMY6tWAb/FKRPYvcGSN1pPiHAIh4yH+CUq9hnOy1om0CiN5mSYElXg7Guji7USEuz3eBgtuLs7yg0GyaaKuFHXna4g9b0JDDWrXmdgN+Ts631hrhya1ksoGM5tySliXVcyGnBJKjCY87A0MD/VhbvcwhnrW41O8EpH1Cxxeo9385xwCoROg71OUeA4nrdaRkxVG8jO03OjjZCDO20nrTNBbcJW58ZxkAS21S6L0hHaX+cGPAVB6zkKXcC+KXzcrRyZdiMyqOtZkFrEmq5id+eWYhaCLpwvTYoIZFepNXyUZfc7XcORnKD2AUHRaL3Pv+ZiCxpOpjyOtop60inpqis3Y6aoJd7dndIQHkR72uDfOOXpquVdJslWitgR15yuI3W+CsGg3BfZ/gAY7H+ydnKwdntTCMiprWZNVzJqsYnY15sZYL1eu6xzMqBAv+ijJ6HO+hEO/QNnB367A9f4HDcETyFQ6k17ZQFp5PTVFJux0ZsLd7RkT6UGEuwPuDtqsWzI3/jlZQEvtisjZjrr9JUTy9+Dsi27QYyj97pQLnrRxqhDsL6pkdWYRa7OKOV5eg51OYWCgF/8Y0JlRQU6EVmxEZL0P638FY6E2jVLIBJSej1DpO4aTdS6cLDeSndGAKsrxctQT6+1EpIcDIW4d7wYXqX0TtcXajBp73tKmout3F7rEh1Bc/LV7PWTx0y6cujl6dVYxa7OKSSmvwV6nMCjIm38mdmFUoCMhFRsQWe825sYicPCB0AkovR6lwncMJ2udSaswkp3egCoq8G7MjVGeDoS42qOXufGiyAJasnlCCMTJ5ahbn9VW0fLujG7COyg9Z6HYyR6YtqrebGFrXhmrMotYnVVMcV0D3g52jArz5b74Tgz1NOKa/zMi8yfYtR6hNmhLZcfcDGGTyXNO4GSFiZPl9ZSWmNEplYS62TMs1J1Ong54Osr0JrU/oq4UdfsLiN1vAaAk3I1uwENyKrp2pN5sYcup3JhZRInRhI+jHSNDfbUpOD2NuOb9jMhcAjvXI1STtlR257kQOolc536crDCR1pgb9Y25cXiYO1EeMje2FPkqSjZLqBaUlB9Qd78MBfsgOBHdtB9Qulwlxze3UZUNZtZlFbMys4iN2SXUmC1EuDnxl06BjA33Jd4uE33WD4hDP0Hxbu3yY+AIlITnMYdMJEOEkVpuJC2vnjpzOU4GHVGeDgwOcSXCwwF7vXzfpfZJ1Fchdr2mTUenWlD634tuwIPy6lo7UVlvYm12CSszCtmYU0qt2UKkuxPTYoIYE+ZLvF0G+qzvEAd/gpK9CMUAQSNQ+r+IKXgSGWowJ8vrOZlnxGgux7kxNw4JdSPc3V7mxlYgC2jJ5giLCXFkEerW57ArTYbIMehmrkGJGCVvdmiDiusaWJNVxPL0IrbmlWJSBT193bi9VwRXhPnSueEAZL4Hm5dA1UmEnZs2NKPbvdQGjOdknTOpZUYyT9RjEWV4Oxro4etEtJcjgS528j2X2jVhNiL2voO65RloqETpeye6wX9HcfG3dmjSJSqqq2d1ZjErMgrZlleGSRX09nXnzt6RjA31Jrp+P2S9C5t+hOp0hJ27NjSjx4PUBowjtcaJk+WncmM5Pk4Gevk508lT5sbLQRbQks0QFhPi4CeoW5+GigzofBWmse/i2GmETBRtTH6NkRUZRSzPKGRXQTkAAwI8md+/M1eEehJctRWR8TGsXAJ1BeAUAOFXo4RfTYXXMFKrBCfKjOQmmVCoINjNniGhbkR7OsrLj1KHIFSzlu82/UebQaj3XHRD/4niHmbt0KRLkFdjZHl6IcszithdUI6iQGKAF4835sagys2IjP/Bip+0ez2cAiHiL1pu9BzKiUqVE+VG8o42oNBASOMUnHLY2uUnX22pzRMWE+LQZ6hbnoKKdJS469BdtxT8etAgb5RpM3KrjSzPKOTX9EL2FFZgp1MYHOTN04O7MibYDZ/S9Yj0N2HfT4iGCnCNgk5/hfCrKXFN4ESFiRNlRorzKtErEOHhwLhIDzp5Onb4CfuljkMIgUj9BXXto1B8BKXb9eiGP4Hi3dnaoUkXKafayK/phSxPL2RvkZYbhwR78+yQrowJdsW7ZB0i4zXYt0zLjW7REHMThP+FEtd+pJQ3cKLMSEluBXoFIj0cGBfVmBsNMjdaiyygpTZLqGbEoc+1wrn8pLYgwHU/ofj31LYLYeUIpYLaetaeLOKX9EL2FVVir1MYGuLDS8O6MTrIFY/i1Yj0l2DnMoSpCjzioOs9EP4XCh27caKsnpRCI+VZZdjrFTo1rgAoxzNLHZHI34u65mFExjqU8JHo5uxCCUqwdljSRcirrWfdySJ+TivkQHEl9nodw4O9eWV4N0YHueBWuAqR/gLs+BlhrtZuAux2L4RfQ6F9HCnl9aQUGKnILMVBrxDl6cigYFciPRywk7mxTZAFtNTmCKEiji5G3fhPKDuBEjsV3bU/yKW224iiunp+TS/k57QCdhdovSnDQ3waPxjccCtarX0wbFuKMNeAVy+UHg8hIqZSYNeZ46VGTuQYqWwoxdGgEO3pyMhwR8LcHeRUc1KHJCoyUTc8jji8EHy6ortuKUrMZDk0zcYU1dbzS3ohy9IK2FNYgb1Ox8hQH14d3p1RQS5abkx7DrYuQ5hrwbs3Ss//Q0RcQ4FdDMdLjaRkG6lqKMXJoCPa04HREY6EuTnIqebaIFlAS21G03R06+dDwX6UmCvRTf0GJaCPtUPr8CrrTfyaUcTSk/lszy9Dh3YJ8qkB0UyK8sO9ZB0i7WXY9pPW0+zVE6Xno4jIayk0RGkfDJlGKhtKcDboiPFypLO3I6Fu9uhkkSB1UMJUi7rtBcT258HBA92Ed1H63IKikx/NtqKi3sTyjEJ+Si1gR4GWG4eGePNMYgwTI31wL16LSHsRtv6k9TR790bpNR8ROY0CfWRjbqyjqjE3dvZ2JMZL5kZbIH9LpTZBZG/Fsu7vkLURQoeiv2kzStgQa4fVoRnNFtZmFfPTyQLWZxdjUgUDg7x4alAc48J88CzfgjnlDQz7f0I0lIFnt8ae5mkU23cmubSO4xla0exk0NHZy5Eu3o6EyA8GqYMTQiCOfYu65mGoyUcZ8AC6wY+jOLhZOzTpAtQ15cZ81meXYBGCgYFePD1Yy40epZsxp7yGYd9SREM5ePZA6fkIIvJaiuyitdyYbjytaO7i5UiwzI02RRbQklWJkmTUdY8iji8B/17opv+MEj1RXrq0Eosq2JpXypLUfFZmFlFtstDL153/S4hhUqQ/AbUHESdfgn3fQF0eOpdIiJuH0ukGyhzjSC6tIzm7jjJjMY4Ghc5eTnSRPc2S1EQUHsSy8l7I3IDSeQq6MWtQvGOsHZb0J07lxh9T81mZUUSN2UJvX3cea8yN/jX7ESdfhL3fQF0+OpcoiLsTpdP1lDp2IbnUSHJmHeX1xbJDoZ2QBbRkFaK2BHXzfxB73wG3EHRXL0Lpdr1cAMVKjpVW82NqHktOFlBQW0+kuxO3dA9nSqdAopRcSP0E8esiRNUJbVqlqOkQNYNS++5k1EJyjpGiuiLs9Qoxp8Y0y3F7ktREGCtQN/wDsfe/4N0F3fXL0XUab+2wpD+RVFrFDyfy+elkPoV1DUS5O3NbzwimdAogQuQgTn4Mvy5CVKU25sYZjbmxG+mNubG4rljLjV6OjIpwItxdFs3tgSygpctKmOsRe97SZtYQKroRT6H0vw/F4Gjt0Dqcorp6lqQW8ENqHkml1Xg6GJgSFcg1MYH0cm1ASf8WsWkRomgbGFwhchrK4Lep9x1BSoWJYyV1ZFdVo1egk6cjA0PciPSQNwJK0u8JIRBJi1FX3Q+manSjX0RJuAdFb2ft0KRzKKqtZ8nJfL4/kc+xsmq8Hey4slMA10QH0dPViJL2DWLjIkTRDrBza8yN71DvO5yUChNJxbXkVGu5MdrTkUEyN7ZLsoCWLgshBCL5B9S1j0BFBkr8PHTDFqC4+Fk7tA6lwaKyLruYb1PyWJ9dgl6B0WG+3N+nE8OD3LDPW444+HfI+hmBgJDxKCMWoYZOIa1GT1JJHWkHi7EICHezZ2SwE10D3HA06K3dNElqc0TpCdQVdyHSVqLETkM39jUU91BrhyWdRYNFZU1WMd+dyGNDdgl6ncKYMF8e6tuJYUFu2OX+ijjwKGT/igAIHY8y4ovG3KgjqbiWtNxiVAFh7vaMDHGiq7/Mje2ZzRbQH330Ea+++ipJSUm89tpr3H333U3bVFXlvvvu45dffkFRFB588EHuvPNOK0bbsYniJNSV9yLSV6NET0I3fRmKb1drh9WhJJVW8U1KHj+l5lNab6Knrxv/TOzMlKgAPKsPI048DVu+RNQXg09flAEvISKnU6h6kVRSx7EjVRjNAn9nA0NC3Yj1dsLFTkddXR0Ock5SSTqNMNcjtr+AuuVpcA1CN30ZupjJ1g5LOoujJVV8m5LLkpMFlNWb6O3rzr8HduHKSH88qg8hTjwJm79C1JeAbwLKgFcQUddRaPHkaEkdyYcrMVq03Dg01I0uMjd2GDZbQPfr14/Fixfz7LPPnrFt4cKFHD16lOPHj1NRUUHfvn0ZPXo0cXFxVoi04xL1laibn0Dseh08IuSHyGVW2WBm6cl8Fqfkcqi4Ch9HO6bGBDE1JohYZyOkLkIs/xhRdkhbSjtmFkrMLGpcu3G0pI6k1DpKjSW42Ono7utMNx8nfJ1/u+wsF7KRpDOJnO1Yfr4FSo+jJD6sLb9t52ztsKTfqaw38VNaAYuP53K4pApfJ3uu7RzEtJggOjvVwYmFiF8/QZQfBqcg6DwHJWYW1S5xJBXXkXSitik39vBzpquvE75OMjd2NDZbQPfu3RsAne7Mv/C+/vprbr/9dvR6Pd7e3kyfPp2vvvqKBQsWnPN8FRUVp/3QOzg44ODgcM79T+1rq78orRm/EAJxeCFi3aPQUIkybAHKgAdRDI4t+nzyPTj7OfcUVvD18Vx+SS+kQVUZFerLu6MjGRHsiV3+Gtj/OCLrJ+2AsKug39NYAseSVmnhSFEtGamF6HUQ4+nIiDA3wn53w8vvY7X11x9svw0XG39z9pe58QKPM9UiNv4TsfM1CEpAN2eP1VZN7ajvwZ+dc1dBOV+n5LE8vRCTKhgV6sM9vSMZHuyJXf4q2DcfkbUUUCD8Kkh4FkvgFZysNHOksI7Myt/nRneZG9uwy5EbbbaAPp/MzEwiIiKavo+MjGT37t3nPSY8PPy07+fPn8/jjz/+p89lNBovLsg2osXjLz2GYc296HK2YOk8DcvwZ8AtFEwCTHUt+1yN5HugTea/JL2Ib1ILSKsyEuriwG3dQrg60o9AkY/+5Gvoty2EuhxUj25Yej2JJXIGxaoXyeUNpBwqpd4iCHDSMyzIiWgPO+z1CqBS/yfx2frrD7bfhubGX1d34b+LMjf+OSVrI4bVd0F1LpahT6L2vQd0BmjG69waOtJ7cC7lv8uN6VVGwlwdmNctlKuj/Aiw5KJPewX9tkVQl4vq0QNL76exREynWPXkWHkDJ36fG4OdiHaXudGWtGZubLMF9LBhw0hKSjrrtn379hEWFnbe438/j/CF/EWRmZmJu7t70/cX0stiNBpxdHS0yTmLWzp+YTYitj6D2PY8eESiu2EV+sgxLRDpeZ6zg78HQgj2FlXyZXIOP6cVIhCMDffjicFxJPq7osv+Ffa+Dzkrwc4Vom6ALnMwe/Yluayew7m1FNZW49x4GbK7jxPeTheeEmz99Qfbb8PFxm8ymS54X5kbz7NvfRVi3aOIfe9C2DB0M37B4NPlMkV6nrg60HtwruP3FFbwZXIOv6QXIRCMj/Dj6cFd6e/vgi77Z9j9AeSu0mbR6HQDdJ6D2SOeY425sai2Ghc7HT39tOFrMjfalsuRG9tsAb1p06aLPjY8PJz09HT69+8PQEZGxhm9KH/k4eFx2ofEhVIUxSZ/uE5pifjV9LWoy2+H8nSUQY+hGzL/sk5L19Heg2qTmR9O5PNFcg7JZdWEuTpyf3wU13YOxsecizj+Gmz8GOrywHcAypD3EZHTyW+w53BRLcmZxVhUQaSHA1fFuBHl6XBJc5La+usPtt+G5sbfnH1lbjw7kbUZdeksqClEN+4tlH53tLl57Nv7e/BHVQ1mfkjN44tjORwvryHCzYkH+nZiWkwQPuYcxPGXG3NjAfgNRBn6ASLiOvLr7TlUXMvxxtwY5enA4Map52RutO02tGZubLMF9KW47rrreO+995g6dSoVFRV8/fXXLF++3NphtTuithh1zcOIQ59qy29P+xHFr5u1w2q3Ustr+PxYNt+fyKPOrDImzJe/949hSKAHutyViM33ILJ/1XpUov+K0uVWGjx6kVRSx8HjtZTUVeFur6d/kAvdfZ1xs5fTK0lScwlzPerGfyG2vwihg9HfsBrFq5O1w+rQUsqr+Twpmx9O5GO0qFwR7svjAzozONADXe4KxKY7ETkrGnPjjSixf6PevYeWG5NrKTVW4e6gZ0BjbnSVuVG6ADZbQC9cuJDHHnuMsrIylixZwnPPPcfSpUuJj4/npptuYteuXXTpol1Ke+SRR+jaVU6b1pLUY9+jLr8D1AZ0kz5A6T23zfW+tAdmVWVtVgmfJWWxNa8MH0c7bu4axg2xIQQZqiHlI8TW9xDVaeAdjzLkXYi6gcIGew4W1nAsvRCLKoj2cmR4mDsR7vY23ZsgSdYkCg9h+elGKE5CN+pZlMSHUXSy2LIGs6qyOrOYz49lsy2vDF8ne+Z2D+P62BAC9VVw/MPG3JgOPv1QhrwHUddT0GDHocJajqUVogpBtKcjI8PdCZe5UWommy2gb7zxRm688cazbtPr9bz99tuXOaKOQdQWo668B3H0K5QuV6Ob8C6Ka6C1w2p3KupNLD6ey6dJ2eTWGIn3c+eV4d2ZGOmPfdk+xN4nEemLtZ0jp6OMXITFuz/JZUYOptSSX1OJq52O/oEu9PCTPSqSdCmEUBE7XkHd8Dh4d0E/ZydKQB9rh9Uhldeb+Pp4Lp8lZZFXU08/fw9eG9GdCRH+2JXuQez5DyLtG1B0EDUDJe4rLD4JHGu8EldQW4mbvextli6dzRbQ0uXX1OsszOiuXoTS7Qb5F3sLO1lRwydHs/juRB4WVTClUyCzuobSw9sJ0r9DLH8bUbgVXCNR4v8DnedQpXhyoLCWwwcLMZpPjW32uuSxzZIkgajOR/3pRkT6WpTEB9GNeOqy3uMhaVLLa/gkKYvvf5cbZ3cLo5unA2R8h/j1LUTRdnDrhNL3Keg8m0rFgwOFtRzZX4DRIojycOCqzl5EXeLYZkkCWUBLF0DUlqCuvFv2OrcSIQRb8sr46Egm67NL8HG047YeEcyMC8GXCkh+E7HmXajNhaDRKKO/R4ROJqvGwv6sGk6WF2GvV+ju60xvf2c8HeWvtSS1BPXkCtSfZoGiQ3fDKnRRrTuzkHQ6IQSbc0v56EgWG3JK8HWyZ17PCG6IDcGXckh+HbHqXe2G6eAxKGN+RIRMJLP6t9zooFfo4etML5kbpRYmf5qk81JPrkRdNhss9bLXuYWZVJWl6UV8djyfpLJqunm78uLQblzZKQD7qmOIvfcgUhcCOu2mwK53Y/HswdHiOvYdLaPUaMbHycDoCA+6+jhiJ5eNlaSWYTGhrluA2P4iSqfx6KZ8huLib+2oOowGi8qStCI+S8kjuayG7t5uvDysG5OiArCvPIrYcxcidZE213b0TShd78Ls3o2kktqm3OjrZOCKSA/ivJ2w08vPLKnlyQJaOithNqKuewyx63WUqLHorvwExS3Y2mG1C1UNZr4+nsNHR7LIr61nRIgPjw/ozMBAT5T8tYg1t2h3jDsFofT5F3T5GzU6Tw4U1nAwowCjWRDj6cDoCHdC3eSNL5LUkkR5GoYfrkcU7kM3+gWUxIfkDdKXSWW9iS+P5/LJ0UwKahsYFerDvxNjGRDggZK3GrFmLiJnJTgHa0PYYm+lRvFgf0ENh9IKqLdoN0yPiXQnxFXmRql1yQJaOoMoPIjlx5lQdgLd2NdREu6WHyAtoLC2no+OZPJFcg71FpWrOgVwY3QAPf3dUNK/Rvz0KqLsIHj3Rhn2CUTNoNCosDenmuNlhegVhR5+TvTxd5GXIiWpFajHl6AuvRnFwRPdjZvQhQ60dkgdQkFtPR8ezuSr4zk0WFSujg7kr9H+9PBzQ0n7ErH9VUT5YW2moeGfQeR1FBhhb3YNKWWFGBSFHn7O9AlwxsNB5kbp8pA/aVITIVTEztdQ1/8dfOLQz9mN4t/D2mHZvMyqOt47lMF3J/Kw1yn8NS6U2d3C8DM0YDryDmx6G1GbDaETUQa8jAgcSXplA3tSqsiqasDdXs+wUDe6+zrjYJB/yEhSSxOqRZvbeesz0OUaTKPfwuAVZO2w2r30ylreP5TB9yfysNfruLExN/oa6rXcuPFtRG0OhE1GGfg6ImA4aY25MbuqAXcHPcNC3enu54SDHMImXWYtUkDv3LmTe+65B4PBwL///W/GjRsHwDXXXMMPP/zQEk8htTJRU4S6dBbi5HKUxIfQjXgaxXDu5XqlP5dUWsW7hzL4Oa0ALwc77u0TxY2xIbiJcsTRpyHpvxhMVdDpBpSeD6N6dCeptI49R0soqTMT4GLH5GhPYrwc5R3jktRKRG0J6pKZiPTV6EY9B4mPgNFo7bDataMlVbx7KJ1f0gvxdrTn/vhOzIwLxc1Sgjj6FBz7LwZTtXbvR8+Hsbh35VhJHXuOlFBqNBPoYseV0Z5Ey9woWVGLFNAPPfQQH330EQ4ODtxxxx1kZ2czd+5cysvLW+L0UisTmZuwLLkBLA3orl+BrtM4a4dk0w4WV/LW/jRWZxUT4urIvxK7ML1zMA7GHMS+hxDHP9LmKO1yC/WdbkfvFcOhojr2ZhRSY1KJ8nBgdLg7IXJ8syS1KpG3B8v306ChWst9UVcghLB2WO3WgaIK3tifxrrsEkJdHVkwMJZrY4JwMGYj9tyPSPkYFAPE3kp9p9vReXbiUFEd+9K13NjJ04ExkR6EuNrJ3ChZXYsU0Hq9nu7duwPw888/M3PmTCoqKuQPeFsnVNStzyI2/gvChqK/+gt5o+Al2FdYwZsH0lifXUKUuzMvDevGlE4BGKrTEDvuQJz4DOzcUHo9Cl3vok7nwa6cCo5kFWFSBV19nOgb6IKvk521myJJ7Z564GNtXnv/nuhv3IjiEW7tkNqtPQXlvHkgjY05pUR7OPPK8G5cGRWAvioVsX0e4sTnYO+B0ms+dL2TWsWd3TkVHM4swqwKujXmRh+ZG6U2pEUKaLPZTHV1Na6urtjb2/P1119z0003sXPnzpY4vdQKRG0Rhh//ishYgzLkcXTD/o2ik0PiL8bugnLe3J/GptxSYjyceW1EdyZHBqCrTEJs+j9E2lfg6IfS92mIm0eNcGJPfg2HiooQQtDDz5mEQFfcHOSKWJLU2oRqRl3zsDbDUJ9b0Y17Uy6M0kp25pfxxv40tuaV0cXThTdG9GBipD+6iiOITQ8j0haDoz9KwnMQexvVqiN7Cqo5VFQIAnr6O9Mv0BU3uVqg1Aa1SMX0xhtvUFNTg6urK6D1SC9atIjFixe3xOmlFiayt6F+fy2KpQHdjF/RRY+3dkg2aX9RBa/sPcnm3FJivVx5a2QPJkT6o5QfQWy4D5H+LbiEoSS+Dp3nUGmxY1dODUeLCzHoFOIDnIlz1+Ht5iKv1kjSZSCMFahLbkCcXIlu3FvoEu6ydkjt0p6Ccl7Zd5JteWXEebny9qiejI/wQyk7jFh/NyLje3AJR0l8AzrPocJsYFd2NUdLKrHTKfQLcCFW5kapjbvoAnrr1q0sXryYEydOABATE8OMGTMYNGgQAIqiMGPGjJaJUmox6r73UVfcDcEDME34FINfJ2uHZHOOlFTxyt5U1mWX0MXT5bcPh4okxPrGwtk1HGXwuxBzMxUmHbuyqjlaUo6DXsegEDd6+Ttjr1Ooq6uzdnMkqUMQZSexfDMFqnLQXf8ruqix1g6p3TlUXMkre0+yIaeEWC9X3hndk7Hhflqnwvq7EenfgWsUypAPIPpGKkw6dmZVc7S4DieDjsEhbvT2d8ZO5kbJBlxUAX3XXXfxzjvvEBsbS2xsLADLly/nzTff5I477uCtt95q0SClSyfM9agr70Hs/wCl310oY16GBou1w7IpyWXVvLbvJCsyiohy//1QjWOIDfdrlyNdwlAGv9NUOO/IrCap8cNhaKgbvfycm1YMlDcrSdLlITI3YfnuGnD0RH/zdhTfOGuH1K4cK9Vy48rMIjp5OPPGiB5MivJv7FS4t7FTIQJlyPsQM4vyBoUdGdUcK6nDyU7H8DA3evq5NK0YKHOjZAuaXUC/9dZbLF68mJ9//pmJEyeetu2XX37h5ptvplu3btx5550tFqR0aURVrnanef5edJM/Qtd7TmOCkn/hX4jsqjpe3XeSH1PzCXF15IWhXflLdCD6mnTE5tnakrIuoSiD3tYuR5p07MioJqnpw8Gdnn7OcjlZSbIC9eAnqL/chhI6BN3Ub1GcfawdUruRWVXHK3tTWXqygDA3J14a1o2rOgWgr05DbHwYcfLLMzoVtqdrhbOznY4R4VpuNOhkbpRsT7ML6A8++IBXXnnljOIZYNKkSbz88su89NJLsoBuI0TWFizfXws6PfqbNqEED7B2SDaj1NjA2wfSWXQsG3cHOxYMjGVGl2Ds6gsQO+5FHP8fOPigDHwDutxCtcXAjqxqjhTX4qiXHw6SZE1CCMTmJ1E3/Vu7WXD82yh6e2uH1S6UGBt4e38ai5Jz8HKw48lBsVzXJRiDMR+x/W7E8Q/ByR9l0FvQeS5VZj07sqo4UqQVziPD3ekhc6Nk45pdQB8/fpyRI0eec/vIkSOZN2/epcQktRD1wMeov86DkIHor/kGxTXA2iHZhFqThY+OZPL+4QwA7ukTxZxu4TiplYh9/0AcfQP0jijxT0C3u6kVDuzKqeZgYS12eoUhIW709neRPc6SZCVCtaCuuAux7z1tUajBf5c3o7WAGpOZD49k8sGhTHQ6hfviOzGnWxiO5nLEnr8jkt4CgzNKv6eh613UqPbsyq7mUFEt9nodw8Lc6OUnc6PUPjS7gHZ1daW4uJiwsLCzbi8pKcHFxeWSA5MunhAq6oZ/ILY+i9Lnb409L3L+zD9jUQXfpOTy6r6TVNSbuLFrKHf2isTLYIGkVxAHnwO1Abo/gNLjIRr07uzOr2ZvQSU6BQYEuxIf4CKXlJUkKxKmOm1lwZSlTUPWpEtjVlUWH9dyY5XJwqy4UO7oFYmnwQRHX0QcfB6EBXo8jNLjQep1buzKq2Z/YTl6RSGxMTfay9wotSPNLqBHjx7NCy+8wJdffnnW7c8//zxjxoy55MCkiyNMdahLb0Yc+xbd6BdREh+SPS8XYHNuKU/vTCG5rJqrOwXycL9OBLs4QOpCxN5/QW0exN6G0ucfWBz8OVBYw87cQsxCEB/gQkKgK44G+eEgSdYk6kqxfHOVdr/HdUvQxUy2dkg2b2NOCc/sTOF4eQ1TowN5oG80wc52kPq5lhuNhRB7O0rv+Vgc/NhfUMPOvEIsAvoGuNIv0EXmRqldanYB/eSTT5KYmMjgwYN54IEHmmbhSEpK4vXXXycpKYnt27e3eKDSnxPVBVi+vRoKD6Kb9h262GusHVKbd6K8hmd3pbAuu4SEAA9+uLI/vf3cEbmrEWsehdL9EDEVpd8zCPcYjpbUse14EdUNFnr4OTMw2BVXOcm/JFmdqMzC8tUEqClA/9d1KCGJ1g7JpqWUV/PMzhNsyClhQIAnS6b0p6evOyJnBWL1Y1B2ECKv1XKjWyeOFtexNbmQWrNKD19nEmVulNq5ZhXQK1asYOzYsWzevJl7772XGTNmNPVuCiEYMWIEmzZtaiqqpctHFB7G8s2VYGlAf9NGlKAEa4fUppXXm3ht30kWHcsh2MWBt0b1YGKEP5QfQV35COSsBL9BKJM2oQQMJr3CyMYjxZTUmYnxcmRqF2+8neTKjZLUFoiS41i+GAM6A/pZW1B85GfQxSozmnhlXypfJecS6urIO6N7Mi7cD8oOo654GHJXQ8BQlCu3gu8A0ivq2dSYG7t4OTI41A0vR5kbpfavWT/lkyZNIi8vj+7du3P99dfz5Zdfkp6ejhCCmJgYfHzk9EDWoGasR/32avCIRD99GYr72cenS9o456+P5/Dy3pOYVJVH+kVzc7cw7E1liO33QvJ72kT/o76BiGsoqTOzMbmUjMp6Qtzsub6rD0Gu8k5+SWorROFhLF9eAU4+6G9YheIWbO2QbJJFFXx5PIeX96SiAo8mRHNT1zDsTaWI7XdD8vvgFoMy+jsIv5riOjMbj5eSWdlAqJs9N3T1IVDmRqkDaVYB7e/vT05ODv7+/tx+++1cffXVDBggp0WzJjX5B9Qfb0AJG4Zu2vcoDm7WDqnN2ltYwYLtyRwuqWJaTBD/1y8aX0cdHHsbse8/ICwoCc9B17upVfVsy6jkcFEtHg56psR4Ee3pIMeTS1IbIvL3YflyLLiHob9+JYqLn7VDskm7C8pZsD2Zo6XVXNc5iEf6xeDjoEDSm4j9TwACpf+LEHcntaqerekVHCmuw9NRz1WdvejkIXOj1PE0q4CePXs248aNY8qUKQBs2bKFoUOH4ucnk5Y1qPv/h/rrPJS4a9FN+QzF4GDtkNqkorp6Xtidyncn8ujh48a3kxPo6++hjeXb8RBUHIMut6D0fRKLgx97C2rYlVuNosDwMHd6+zujl/OVSlKbInJ2YPlqPHh3QX/9ChQnL2uHZHMKa+t5fvcJfkjNp6evG99fmUAfPw9E9q+InQ9BZQp0+RtK3/9gsfdlT341u/Jq0CswItydXn4yN0odV7MK6GeffZaRI0fy888/I4Rgzpw5VFVV4e/vT8+ePenVq1fTV58+fVopZEkIgdj2HOr6+Sh970Q37g0UnbxZ448sqmBRcjYv7UnFoFN4alAsM7qEoKvJQF0zBzKXQOAIlBELwbs3J8vrWZ+i3SDY29+FgcFyZg1JaotE5kYsiyeDfx/0M35GcXC3dkg2xaIKFh7L5uW9qdjpdDw7JI7rOgejVKejrp4FWcsgaDTKyK/Aqyep5fVsSC6i2mShj78LiTI3SlLzZ+EYP34848ePZ/HixRw4cACTycShQ4c4ePAgBw8eZMWKFaSkpGA0Glsj3iYfffQRr776KklJSbz22mvcfffdTdsWLFjAf//7X4KDtbFw3bt3Z9GiRa0az+UihIq6+iHErtfQDVuAMvRf8tLZWRwtqeLxrcc4UFzJ9V2CeaRfDJ52Khx6FnHgGXDwRhn5JUReR5nRwvrjZWRU1hPh7sA1neUNgpLUVqlpq1G/uQolZBC6635CsZfrDjTH4eJKHt92jEPFVdwQG8Ij/aLx0FvgwFPaXPcOviijvoaIaZQazaxvHOcc5eHA1FhveYOgJDW66N+E/Pz8pv+HhoaetrS3yWS6tKguQL9+/Vi8eDHPPvvsWbfPmjWLl156qdXjuJyEakZdNgdxeBG68W+j6yeXS/+jGpOZ1/el8fHRLKI9nFk8qR8JAZ7acI1t90J1OnS/H6XPP2nQObMju4p9BTW42uu5KsaLTnKcsyS1WWr6OtRvpqBEjEI39TsUOydrh2Qzqk1mXt17kk+Tsujs6frbULbsXxHb74OaTOj+IErvx2lQnNieVcX+whrc7PVc3dmLTp6O1m6CJLUprfKnpJ1d669617t3bwB0upa5jFRRUYEQoul7BwcHHBzOPab41L6/P6Y1CYsJ8dONiOM/oPzlS5Su0y/puS93/K3hj21Ym1XMv7cnU2I08WB8FHO7h2NnzEVdeytk/ACBI2DM9+DRlWOlRjZmF9FgUUkM0ib7N+iU0853ueO3NbYeP9h+Gy42/ubs31Zyo8jeivrNFAgbjjL1OzA4tsr7Zus/E3BmG1ZlFrFg+3HK60083Deaud3DMNTl/DaULWgUjFkCHrEklRrZlFVEgyoYGOxK3wAtN17O18PW3wNbjx9svw2XIze222sxX375JatWrcLHx4d//vOfjBo16rz7h4eHn/b9/Pnzefzxx//0eVp7qAoAlgYMv9yMkvYr5smLEJFXQl1di5z6ssTfynLLq3huXxq/ZJYwJNCDD0d2I8zZDo68iTi4AAzOmAZ+hBp+LRUNKpuOlZBTYybKzcCgQBfc7HWY6o20/nWTs7P198DW4wfbb0Nz469rRv5oC7lRKdiH4btJCP94zJMWgUmAqWVy4LnY+s8EQE55Fc/sTWNFVgnDgjx5vG83Qp3tEIdeRRx6EgyumAZ9ihp2DeUNKpuOFZNbY6GTux2DAh1xtZO58VLYevxg+21ozdzYZgvoYcOGkZSUdNZt+/btIyzs3HMd33777Tz++OPY2dmxZcsWrrnmGnbt2kVERMQ5j8nMzMTd/bcbUS6kl8VoNOLo6Niql/yFuR512SxIX4Fu6ncYOl/ZMue9TPG3JlVV+TElh2f2piMQvDKsG1d1CoDyw7BuHhTvhC63Qb9n0Nt5sDevmt35NbjY6bk6xosoT+vOWmLr74Gtxw+234aLjb85w+ysnRtF4SHUH64C367oZ/yMXStP1WnrPxOg5cbvj2fz3L4MAF4d3o0pUQFQegDW3Q7FeyDuduj7FHqDO3vyqtmTrw1l+0tnLyI9ZG68FLYeP9h+Gy5HbmyzBfSmTZsu+tjAwMCm/w8ZMoT4+Hh279593gLaw8PjtA+JC6UoSqv9cAmzEfH9NEhfg+7aJeiiJ7T4c7Rm/K2psLaef247xqrMYiZG+LFgUCy+dipi3z/h0Evg3hll0gaUgKFkVNSz9ngxlQ0WEgJdGBDkhp2+7bTZVt+DU2w9frD9NjQ3/ubsa83cKIqPoX45Fjwi0F+/HMXx8s22Yas/E/k1Rv65LZk1WcVMjvTn3wNj8bEzI/bMhyOvgEccyqSNjSus1rM2o5iqBgv9g1wZEOTaNJStLbDV9+AUW48fbL8NrZkb22wBfSmys7MJDQ0FICUlhf3799OzZ08rR9U8wlSH+u1fEFkb0U1fii5qrLVDahOEEHx/Ip8ndx7HXqfj1cFduKpLKBRsQmz+G9RkovT5J/R8BKOwY8PJcpJK6gh1s+eqzl74OLX++HxJki6dKEvVlud29kN/w0oUR09rh9SmCSH4JiWPp3el4KjX8fqQWK7sHAL56xFbboPaHJT4BdDjYYyqgfWpZRwrNRLmZs9f5MxDktRsF/0bc/DgQXbs2EFhYSGgrVKYmJhIr169Wiy481m4cCGPPfYYZWVlLFmyhOeee46lS5cSHx/P448/zp49ezAYDOj1et5++226dOlyWeJqCcJsRP3mKkTOVnTTf0YXOdraIbUJxXUNPL41iVWZxfwlOpB/9O+MY0O5dgf5sbchYCjK2J/AvQvHS42syyxCFYKxkR5093Wy6b+iJakjEdUF2gqD9i7oZ65Gcfa1dkhtWlFtPX/feoy1WcVMiwni8f4x2NeXIbbdDcnvQsAwlLE/g3tnkkuNrM8sQxWCcVEedPORuVGSLkazC+iCggJmzJjBxo0biYyMbFqFsKioiPT0dEaMGMHXX3+Nv79/iwf7ezfeeCM33njjWbd9+umnrfrcrUlYTKg/zEBkb0Y341d0ESOtHVKbsCqjiPlbtTHx747uxbgIP9Tctaib/wbGApTE16DrXVSbBGtPlHGyvJ4YL0dGR7jjYicXmZEkWyHqq7AsngTmevQ3r0NxDfzzgzqw5emFPL71GDoFPhjTizHhfqg5q7Urcg3FKIlvQNc7qG4QrEkpI62ini5ejoyUuVGSLkmzC+h58+ahKAqpqalERUWdti0tLY1bb72VefPm8cMPP7RYkB2FUC2oy2YjUn9Fd91PsngGqhrMPLnzON+m5HFFmC/PDOmKj6EBddvdcOwd8BsC45eDewyHimrZlFWFQa8wJcaLGC85b6kk2RKtA+E6KE1Bf9MmFI9z37fS0VU2mPnP9mR+SM1nXLgfTw+Jw1tnRN16ByS/j+o/DP2w1eDWiYONudFer3BVjBfRMjdK0iVrdgG9atUqtm3bdkbxDBAVFcUrr7zCkCFDWiS4jkQIgbr8TsTRr9D95etWuWHQ1uzML+OhTUepqDfx/NCuXBsTBAWbEZtmg7EQEl+nIWI2Jr0jqxtXy+rh68SwMHe5zKwk2RghBOovf0Okr0U341eUgN7WDqnN2pZXysObjlLVYOalYd24JjoQCjYiNs0BYzEMfAtT+E3U6h1ZlVxKVlUDPf2cGRbqhoPMjZLUIppdQHt6epKWlnbOsc5paWkXdcd2RyaEQF37f4j976Ob/BG6rtdaOySrMqsqb+xP4+0D6fQP8OTLCX0JcdYj9jwOh16AgCEoE1YhXCJJzqtka0Ex9jqFa7p4W336JUmSLo668V+IQ5+iu3oRuqgx1g6nTTKpKq/tO8m7BzNIDPTkxWHdCHbSIXY/Codf0cY6T1iDcIngWF4l2/KLcTAoTO3iTYTMjZLUoppdQN99993cdNNN3H///YwaNapprHNhYSHr1q3j9ddfZ/78+S0eaHsmtjyN2PESurGvo+s9x9rhWFVOdR33bzjC/qJKHuzbidt7RqKrOIpYNgvKj6D0exp6PEytBVafKOdkRT1dfRwZGe4he50lyUape99FbHkK3egX0HWfae1w2qSsqjru23CYQ8VVPNwvmnk9I1DKDiNW3wQVx1ASnoPuD1BjhlUnykmvqKebjxMjw91lr7MktYJmF9B///vfCQwM5K233uKZZ55BVVVAW1K7d+/evPbaa8yZ07GLwOZQd72BuvGf6EY8ha7/vdYOx6p+SSvg71uP4W5v4KuJfenr7w5H30Ds+Tu4RaNcuQ3FJ57kkjrWZlSgUxTGhznTNcBd3kUuSTZKPf4T6oq7UBLuQUl82NrhtEnLThYwf2sSXg52LJ7Ujz5+bnDkVe2qnHsXlCk7ULx7c6wxN+p1ChPCnYnzl7lRklrLRU1jN2fOHObMmUNDQwMlJSUIIfD19cXe3r6l42vX1KRvUVfdj5L4MMrgjttrX2e28OSO43x1PJdJkf48MzgON0sJYuV0yF0D3e5D6fc0DYoD6xrnde7s5cjocHcw11s7fEmSLpIoOIC65AaUzleju+JVWez9Qa3Jwn92JPNNSh5XRgXw1OA43MyFiOXTIH8ddH8Ape9TNGDP2sZ5nWO9HRkZ7g4mmRslqTVd0szp9vb2BAUFnfZYTU0Ne/bsYfjw4ZcUWHsnsrei/nQjSrfr0Y1+vsN+cKSW13DXukNkVtXx7JA4pncOhtyViI2zQdGhjFuOEjKW3OoGfj1ZjNGkMiHKgzgfJwDqzNaNX5KkiyNqirB8ezV4d0F31ecoOjml2u+llFdz59pD5NYYf7uJOme5lhv19ijjV6EEjyanqoHlJ4sxWlQmdvIkzscJIQR1F74isSRJF6HFlx46ceIEo0aNwmKxtPSp2w1RchzLN1dByEB0V36MonTM8Wk/pxXw2OYkglwcWHJVf2LcHbTlZg+9ACHjUIZ9inD0Y0dOFdtzqwl0sWNaF288HbUfWyGElVsgSdLFEJYGLN9fC6Za9DduRLF3sXZIbcpPJ/OZv+UYoa6O/DRlAJ3c7BC7/k9bijt0IsqwjxEOvmzLqWJnbjVBrnZc28kbDwe5mqAkXS7N/m07Neb5XGThfH6iphDL1xO15Wmn/YBi6Hh3RjdYVJ7blcInSdlcGRXAs0PicK7PQfzyVyjehZLwPPR4kMoGleXHSsirNpEY7EpisCu6DtpTL0ntibryPsjZhv6va1E8wq0dTptRb1F5eudxFh7L4S/RgTw1KA4nY4aWG0v3ofR/CbrfR0W9yvKkEvJrTAwMcWVAkMyNknS5NbuAtrOzO+92IUSHHY7wZ4SpVut5NtVoy9M6eVk7pMsut9rIPesPcbikigWJXbipayhkLkFsvgXsPVEmbUDxH0hKaR2r0itw0Ou4Ls6HEDc5vl6S2gN1zzuIfe+im/Q/lLCh1g6nzcipruPudYdJKq3iyUGxzIwNgYzvtBUFHXxQJm1C8evP8cbc6GjQMb2rD8GuMjdKkjU0u4B2c3Pj73//OwMHDjzr9pSUFObNm3fJgbU3QrWgLrkRig6hv3EDiueZC9G0d5tzS7l//WEcDDq+ntSP3j4uiN1/h8MvQsQ1KEP+h8XOg00ZFewvrCXGy5GxkXJ6OklqL9T0dair7kVJuBddn1usHU6bsTGnhPs3HMbFzsDiyQn09HZC7HoEjrwKkdeiDHkfi8GdjRkVHCispYuXI1dEeeCgl7lRkqyl2QV0QkICqqoyYsSIs2739PSUY1PPQl37CCJlCbprl6AEJVg7nMtKCMFHR7N4dlcKQ4K8eXVEd7xEGWLFNCjY2HhZ8n4qGyz8fKyEoloTo8Ld6e3vLK9mSFI7IcpOov5wHUr4SHRXvGztcNoEIQQfHM7khT0nGB7iwyvDu+NhKUEsvxoKt6IkvgZd76ai3sLPScWU1JkZHeFOLz+ZGyXJ2ppdQN91113U1taec3t4eDgff/zxJQXV3qgHP0XsfBXduDfRdb7S2uFcVvVmC49vPcb3qfnc1iOCR/pFoyvejlg3A4QZZcJqlMDhnCgzsjKtHAeDjhldfQh0kZclJam9EGYjlu+ngaMnumu+RtHJm92MZguPbUnip5MF3NkrkgfiO6Er2qrlRkCZuAYlYCgppXWsTK/AyaBjRldfAlzOP4xSkqTLo9lZ7Jprrjnvdi8vL26++eaLDqi9Ebk7UX+dh9L7FpR+d1k7nMsqv8bI7WsPkVxWzavDu3NVpwBI+i9i10Pg2x9l1NeojkFszqxkb0ENMZ4OjI3ylEM2JKmdUVfdD8VJ6GfvQHHytnY4VpdbbeT2tQc5UV7DGyN6MDnKH46+qQ3b8B+MMvJLVMcANjYOZ+vcOJxNrigoSW2H7AZoRaI6H8t3UyEwHt34tzvUJbe9hRXcsfYgBkVh8aR+9PC0R2yaA6mfQ9d7UPq/QK2q55fjpeRWNzAizJ34AHlZUpLaG/XIl4h976Gb+B5KQG9rh2N1uwvKuWPtQRz0er6dnEBXDwNi4yw4+YW2MErCs9RadCxL1mbZkMPZJKltuugCeuPGjWd9XFEUHBwc6NSpE76+vhcdmK1rmudUqOinftehpqv7LiWPx7cm0dPXnf+O7omvKEUsnwClB1CGf44SPZP86gaWnihGFTAt1ptQt47z+khSRyFKklF/vQ2l+0yUPn+zdjhWt/h4Lv/cdow+fh68PaonPqIY8etUKDuMMuILlE4zyK1uYNmJYgRwXawPwXIGIklqky66gB45cmTTX8Snbhr8/feKojB+/Hi++OILPD09Lz1SG6OuvBfydmkzbrgFWzucy0IIwav7TvLWgXSmdw7miUGx2JXtQ6y5BhAoE9ej+PXnUFEt6zIq8HO2Y0qMF672cgUySWpvhKkO9fvrwDVE633uwD2oqhC8vDeVdw5mcEOXYP49MBa70j2ItVMBHcqkjeATz8HCGtZlVhLgYseV0TI3SlJbdtEDqpYtW0b//v1Zvnw5JSUllJSUsHz5chITE/nuu+/YsGED2dnZPPDAAy0Zr01Q976nXbKc8A5KyNmn+2tv6i0qD2w8wlsH0vm/ftE8OyQOu4xvEL+MAOdglCk7sPgksDq9gtXpFXTzdea6OB/5ASFJ7ZRYdS+UnUA/9RsUe1drh2M1RrOF+9Yf5t2DGczvH8NTg+OwS/8K8etIcAnXcqN3PKvTK1iTUUlPP2eui5W5UZLauovugX7wwQf5/PPP6d+/f9NjY8eOxd3dnVmzZpGcnMybb77JzJkzWyRQWyGytqCuvAel393oes+1djiXRZnRxLy1BzhYXMVbI3swMdIPsW8BHHgKOt2AMuQD6oQ9PyWXUFhjYmykBz38nK0dtiRJrUSX9CXiwIfoJn+I4t/T2uFYTYmxgXmrD3K0tIq3R/VkfIQvYu8/4OBzEH0TyuB3qRV2/NQ4fee4KA+6+8rcKEm24KIL6MzMzLM+rtPpyMrKAiAiIoKKioqLfQqbI2qKsPwwHUIGorviFWuHc1mkV9Yyd9V+KuvNfDGhL/HeDoj1N0D6dyj9noGe/0dJnZkfU0qwqILr4nwIkitnSVK7JYqT0K+5F6XHLJRec6wdjtWcrKhh7qoD1JjMfDmxL7287BDrpkPGjygJL0CPBymuM7MkRbsXZHqcD4EyN0qSzbjoAnr8+PHccsstvP766/Tt2xdFUdizZw/3338/EyZMAGDfvn3ExMS0WLBtmRAq6rKbwdKA/i9foejb/1ydewrKuW3NQbwc7fh+Sn/C7GoRK6ZAyV6U0d+iRPyFk+VGfkktx9NBz9VxPrg5yMuSktSeicOfgXs4Sgebeej3duaXMW/NQfydHVg4vj8hdtWI5ZOg7CDKmB9QwqeQWmbk15PleDnqubqztxyyIUk25qIL6I8//pj77ruP8ePHY7FYANDr9cycOZPXXnsNgJCQEN55550WCbStEzteQaT+im7Grx3ipsE1WcXcve4QvX3deXdMLzzqMxE/T4aGcpSJa8B3AHvyq9mYVUW0pwMTOnliL5edlaR2TxnxDKZed2Gwd7F2KFaxKqOIezYcpp+fB++M7ombMR2xbDKYqlAmrgOffuzKq2ZzdhUxXo5MiPLATuZGSbI5F11Ae3p68umnn/LWW29x8uRJhBBER0fj5ubWtM+AAQNaJMi2TuRsR13/d5SB/4cueoK1w2l135/I49HNSYwJ8+X1Ed2xL9uNWH012HuiTN6C6tqJtekVHC6uo3+QC0NC3DpsT5QkdTSKooCTj7XDsIrFx3OZvzWJ8RH+vDK8O/YlO7Tc6OiHcuVWLC6RrEmv4GhxHYlBrgwKcZW5UZJs1CUvpJKdnU1aWhoATk5OxMbGXnJQtkTUlaH+eD0EJaAb8ZS1w2l1HxzK4NndJ5jRJZgnB8Wiz1qK2PBX8IlHGfMDDQZvljYujjI+yoNu8oYYSZLaOSEE7x3K4IU9qcyMDeE/A2PRZX6vLZDi2x9lzPfU6z1ZeryUvOoGJnTypKuPk7XDliTpElz0daPCwkLGjx9P9+7dmTt3LnPnzqVbt25MnDiRoqKilozxrObPn0/Xrl3p3bs3AwYMYO3atU3bVFXlnnvuITo6mpiYGP773/+2ThBCoP5yK9RXtPtxz0IIntuVwrO7T3Bnr0ieGRyH/vj7iLXTIGwSyviVVOu8WHyshMJaE1NjfWTxLElSu6cKwTO7TvDCnlTu7RPFk4Ni0R17G7FuBoRfjTJuOdV4sPhYCcW1JqbF+sjiWZLagYsuoO+8804qKytJSkqitLSU0tJSjhw5Qnl5OXfddVdLxnhWw4YNY+/evRw4cIAPPviAadOmYTQaAVi4cCFHjx7l+PHj7Ny5kxdeeIFjx461eAy6A+/B8R/QTf4YxSOixc/fVphVlf/bnMT7hzP554DOPNS3Exx8FrHtLuh6N8rIryhpMPDV0WLqLYIZXX0IlatnSZLUzplUlUc2HeWjI5ksSOzCfX2iEPufQOy4T1uWe8RCihv0fJlUjMkimNHVlxCZGyWpXbjoAnrFihW88847pw3ZiIuL4+2332b58uUtEtz5TJw4EScn7a/4nj17YrFYKC4uBuDrr7/m9ttvR6/X4+3tzfTp0/nqq69a9PlF/l70m/6OknAPuti/tOi525IGi8rd6w6zJDWfV4Z3Z3a3UMSuhxF7/4kS/x+UxFfJrjKx+FgJTgYd13f1wcep/fbES5IkgbZ41N3rDrEsrYDXRnTnpq4hiB33w/4nUPo9g9L/BbIac6OLnY7ru/ng7XTJoyYlSWojLvq32c7OjpqamjMer62txWC4vEni448/Jjo6mtDQUECbozoi4rce4cjISHbv3n3ec1RUVDQtSQ7g4OCAg4PDOfcXJ1cifLrByOdPO85WnIr5fLHXmy3ctf4wW3JL+e/onowO8URsvhVOfAqJb0DXO0kurWNlWgUhrvZMjvHEQa+7bK/HhbShLZPxW5+tt+Fi42/O/s3OjR3gNTWaLdyx7hDb88p5Z3RPRgZ7IDbOhpNfwKD/QuxtHCupY2V6BWFu9kyO1mYhkrnxwsj4rc/W23A5cuNFV7pTp05l7ty5vPHGGyQmJgKwfft27r//fqZNm3axp20ybNgwkpKSzrpt3759hIWFAbBmzRr+85//sGrVqtP2+f2dzRfygoSHh5/2/fz583n88cfPfUD8fdDzNsxmAea6Pz1/W3Vq2Msf1Zkt3Ls5mb3Flbw5LI7BXnosa65Fl7sc08APUSOmcyi7nK35Rjp72DEi2BG1oR5rvBLnaoOtkPFbn623obnx19Vd+G9qs3PjRcbU1pwr/lqzhXs2HeNASTX/HR7HQA8dltVT0eWtwjToY9TwaRzIKmd7gZEunnYMD3bEInPjRZHxW5+tt6E1c+NFF9BvvPEG999/P1OmTDltHui5c+fy8ssvX+xpm2zatOlP99mwYQNz5sxh6dKlpw0lCQ8PJz09vWmZ8YyMjDM+BP4oMzMTd3f3pu8vpJfFaARHR0ebnIZIi9941virTWbuWp/E4ZJqPrqiD4m+Blg7FQq3wejvsA+dxI68GrbnG+kX4MLQUOtMxXS+NtgCGb/12XobLjZ+k8l0wfteXG5sn69pVYOZO9cdJam0hk/G9SHBSwdr/gLFu2DMj9iHjGdbbjU7C4z0D3RhsJWmqWvP74EtsPX4wfbbcDly40UX0M7Ozrz//vu88sorp80D7erqerGnbJaNGzdy0003sWTJEnr37n3atuuuu4733nuPqVOnUlFRwddff/2n47I9PDxO+5C4UIqi2OQP1yl/jL+ywczclQc4Xl7Np+P60NdLh1g5CcqPoIz7FQKGsTGrir0FNQwJcWNA8OV5v8+nvb0HtsbW4wfbb0Nz42/OvjI3airqTcxeuZ+TlbV8Nj6ePh4CsWoiVBxDGb8C/AezPrOS/YW1DAt1IyFI5sZLJeO3PltvQ2vmxmYV0P/6178ueN8nnniiOadutltuuYX6+nrmzJnT9Njnn39Oz549uemmm9i1axddunQB4JFHHqFr166tGk97UF5v4uaV+8iorOPz8X3p5a4iVk6AimSU8SsRvgmsaVwgZXSEO739O+ZKY5IkdSxlRhOzVuwjp6aOhePj6eFm0XJjZSrK+FUIn76sTKsgqaSOMRHu9JK5UZLavWYV0BcyrAKaV8FfrJSUlHNu0+v1vP32260eQ3tS2WDm5pX7yK4y8sWEvnR1Nf/uA2Ilqndffk0t50SZkQlRHnSVczxLktQBVNZrxXN+rZFFE/oS59yAWD4eajJRJqzG4tWbX0+UkVpez8ROnsTJOZ4lqUNoVgG9bt261opDsqJqk5m5K/eTUVnHogl96epy+geE2bM3y06UkVVZz5UxXsR4OVo7ZEmSpFZX1WDm5pX7yampayye6xErxkFNjpYbPXrwU0opOVUNTOnsRbSnzI2S1FHISSk7uFqThVtXa2OePx/fl24uRsSK8b/7gOjJkpRS8mpMXN3ZmwiPc988JEmS1F7UmMzMXXWAk5W1LBofT5xTHWL5OKgrQJm4BpN7N35MKaOwxsQ1XbwJc5e5UZI6EllAd2BGs4V7Nx7kcEkVn47rQy+3esSvV4CxEGXiWszuXVmSUkp+4weEXF1QkqSOoM5s4e4NSSSXVfPZ+Hi6u9Qhll8B9aUoE9dicottKp6ndvEmWOZGSepwZAHdQdVbVO7fcpy9RZV8PK4PfT2FNmyjvkgrnt1iWZJSRn6Nib/I4lmSpA6ivnEO/EMl1Xw8rg+93cyNubFMK55dO/Pj8TKKamXxLEkdmSygOyCTqnLv+sPsLKzgf1f0ZoCPoXHYRnZT8fxjShkFjT3PIfIDQpKkDqDBonLX+sPsLa7kwyv6kOClIFZMgto8lEnrTiuer4n1JthV5kZJ6qhkAd3BqELw6OYkNuSU8PqQWAb7OyJWXQnlSdqY58ZxfbJ4liSpI1GF4OFNR9mSW8qbw+IY6GuPWD0ZKk+gTFiDyS1OFs+SJDWRBXQHIoTgmV0pLEnN57UR3Rnu5wBrr4WS3SjjlmP2ipfFsyRJHY4Qgid3HOfntALeGtmDIb522uqrJftRxq/A5NlLFs+SJJ1GFtAdyHuHMvjoSBYLErswOcIby5rrIH8DytilWPwG81PK/7d359FR1Xcfx983ewLZyMKWjYQkkBATkFV2AZHNhRaqVVRc+rhh7VNtK4hgq7YVW7UtauuCVK2iFaVqVRBZw45ACAQIWQiBBLKQjSxMkt/zB4d5QNYgZDLk8zpnzmHmd+fm85sk3/ly85t7S9U8i0ir80paLvMz8nlmQDyjIwJp+GYiHE7Fuu6/NAT34z97So+veY5vR0c1zyICuDg6gDSPD/cc5PnNWUxLjmJKt06w+m5cDn4Fwz/EdLyWL7OOcKDyGDfGBqp5FpFW4/3dB/jTd9n8omc0t8Z1hJV34FL4DQz/N42hQ/gi6wgHq45xY1ygmmcRsVMD3Qp8k1fE9DUZ3BrfmZ+ndMFs+AXkLMA24C0IG8vinHKyy49fJEXnMhWR1uKr3MPMXLuLKd3CeOiqSMy6h2HfQmwD3obOo1mcU0ZueR0TugYS5qvaKCL/Tw30FW5D4RGmLU/nuogQfts/Hit9DmTMhf5/oyHsJpbvrySjpIbruwQQratoiUgrsbaglEdXpDM2KpRZ/eOw0p6D3f+Aa16jofMEvs2rYHdpLWOiA+ii2igi36MG+gq250gV9y1No1eIPy8O7YFL1j8xm6dDylNY8T9j4+E6th2uZkSkP/FB3o6OKyLSLHaVVvE/S9Po0yGQOYMTsTLfxGyZhdXrd1ixU9lwuI7tRTWMjPInrp1qo4icTg30Faqopo57vtlGpzZevDbiKjwKvsak3gdx92KlPMWmwqNsKa5jcJgvV4X6ODquiEizOFxdxz3fbCXC15tXr03C48AXmDUPQLcH4Kon2FBQxdbiOoaE+9IjRLVRRM5MDfQVqLa+gf9ZmsaxhkbeGJlM2/LvMMsmQ9hYrAFz2Vlcw+r8SnoFe3J1hzaOjisi0iyqbQ3ct3QbjQZeH5lMmyObMMtvhYgbsfq9THpxDWsOVHF1iCe92qs2isjZ6TR2V5gTFwPYVVrFB2OvplPDfsySCRCYjDXsX+yrbGBJbjk9gr3pHeru6LgiIs2i0Rj+d+UO9pYdZcGYq+lgy8V8MwGCr8Ya8g45FTaW5pZzVYgPV4forVFEzk1HoK8wf/4umy9zD/Pi0ESS2tRilowFr2CskYs4XOfO53uPEOXvybWRfliW5ei4IiLN4vlNe1mSV8TLQ3uQ6FN9vDZ6d8Qa8SmH6lz5IquM6ABPhkX4qjaKyHnpv9lXkI8yD/JKWi6/6d2V6zr7Yb4eBfVHscZ/Q4Xlz6eZJbTzcmNcTAAueoMQkVbig90H+Ed6Hk/2jWVEp7aYr0ZAwzGsMV9QbnxZtKeEEG83xkQH4qLSKCIXQA30FWJtQSkzUndxS1wn7k0Mx6y+C0o2Y41ZRp1XOJ9mFOPuYnFjXCDuri4YYxwdWUTkslt9sJSZa3dzW7fO3NU9DLNqCpSmYY1dTq1nZz7JKMbD1eLG2Ha4u1qqjSJyQbSE4wqQW1HNg99up1/HQJ4eEI+1/Y+Q9R7WoLdoCOrLosxSqusbuTmuHW3cXR0dV0SkWWSXV/PQsu0M6tSOWf3isNKehewPsIa8TX3g1SzKLKWuwXBzXDu83fV2KCIXThXDyR211XP/0jQCPN3527AeuOV9gvnuSUiZCV1+wlfZZRw6auOm2HYEeukPDiLSOlTZ6rl/6TZCvD14eVgPXPf9G7NlNlbP32Iif8SX2UcoqqnnpthAAlQbRaSJVDWcmDGGX6/OIL+qloXje+NXmYZZeSdETcJKeYo1B6rIPFLLhK6BdGzr4ei4IiLNotEYHlu5k8LqOj4Z3wff8u8wq6ZC9E8heTqp+ZVkl9UxITaQDqqNInIRdATaif0jPY//5h5mzuAEYj3KMUtvgsAeWIPnsbu0jg0FVQwK86VroC5DKyKtx6tp+1icV8SfBicS7VaKWXoztOuJNfB1Mkpq2FR4lMHhvsToEt0icpF0BNpJrTpQwpzNe3ngqkiuD/PF/Hco4II14hMO1bqyOKeE7kHe9NaFUkSkFVmeX8yfv8tiWkoXRnZug/nienDxwBrxMYW1LnyTe4TEYG9dKEVEfhA10E4or7KGR5anM7hTEL9IicasuQfKdmCNW81Rt1D+s7OYUB93Rkb563ymItJq5FZU8+iKHQwPD+aR5KjjZyOq2IM1LpUql2D+s6uY0DbuXBup2igiP4zTLuGYPn063bt3Jzk5mb59+/Ltt9/ax2bPnk1oaCgpKSmkpKRw2223OTDppVVta+D+pWn4e7rz0tBEXDP/AXv/iXXNqzQEpvCfvUewgAmxgbjphKYi0kqc+EB1Oy93/jw4AZddrx4/G9HAf1Dvn8R/9pbi6mIxoatqo4j8cE57BHrw4MHMnDkTb29vtm3bxrBhwygoKMDL6/iatjvuuIMXXnjBwSkvLWMMT6zJYF9lNR+P64NfxXeY9Y9CtwcgZgpLsssoqbExuVuwTlcnIq2GMYbfrM7gQFUtCyf0xrdsI2bD/0LCI9DlFr7OKqO0toGfdA9SbRSRS8Jpj0CPGTMGb29vAJKSkmhoaKC4uNjBqS6v93cf4LPsQ/xxYALx3tWYbydBUG+svn9mY8FRdpXWMrpLAO3buDs6qohIs3lnVz5f5B7m+cEJdPWowiybDKH9sfo8z/qDx89GdH0Xf0J9VBtF5NJw2iPQJ5s3bx4xMTGEhYXZH3v//fdZsmQJQUFBzJw5k+HDh59zH+Xl5adcgcrT0xNPT8+zbn9i2+a6atWu0ip+uz6T2+I7MzayHWbxaGi0wbAPyK1sJPVAJf06tiE20OuCMjV3/svB2eeg/I7n7HO42PxN2b6l18YdJZU8tyGTO7qHMTo8APP1KMDA0PfJLq9n7cEqBnRqS1fVRqeh/I7n7HNojtpomRb66gwePJiMjIwzjm3ZsoXw8HAAli5dytSpU1myZAnx8fEAFBYWEhQUhLu7O6mpqdx8881s3LiRyMjI0/ZVUVGBv7//aY9Pnz6dGTNmXMIZXbxqWwM/WZKGu4sL749Kok36U7ju/hvHhn1ORcA1fJxdRXtvV66P8NEHY0TkvCoqKujYsSPl5eX4+fmddZuWXhuP2hqYtDiNtu6uvDuiBz5p03HN/DvHhn9BuX9/Ps6qpGMbN0aHqzaKyPldSG08ocUegV61atV5t1mxYgVTp07ls88+szfPAB06dLD/e+DAgfTs2ZNNmzadsYE+IS8v75QX60KOstTW1uLl5XVZC7Mxhic3ZXCoxsaiCb3xL/4Kdr0MfebgFj6SpbtK8XR1YWzXdni5XfiKnObKfzk5+xyU3/GcfQ4Xm99ms13wti25Nj6xYSeldTbmXZeC3+HPYM9c6PsibmHXsnRXCd7uroyJUW10NsrveM4+h+aojS22gT6flStXMmXKFBYtWkRycvIpY/n5+fblHJmZmWzdupWkpKRz7s/f3/+8/9s4E8uyLusP18d7C/g0q5A/D0kk2qUQs/qe41caTPwFK/dVUFxjY3L3ILwv8oMxlzt/c3D2OSi/4zn7HJqavynbttTa+GHmQT7LOcTLQxOJIh+z5mcQfStWwjSW7yunpKaeW7oHqzY68RyU3/GcfQ6XszY6bQN9zz33UFdXx9SpU+2PvfPOOyQlJTFjxgw2b96Mm5sbrq6uzJ07l7i4OAemvTiZZVXMWrebSbEduTGqHeaLG8C7A9ag19lVUkNaUTUjIv3o0EaXohWR1mP3kSpmr9vDrXGdGB8ZiPl8HLQJw7rmNXaW1JBeVMOoKH9C9YFqEblMnLaBzszMPOvY/PnzmzHJ5VFT38C0ZemEtfVmVr94zObfwJE0rHGplNi8+GZfMd2DvEkK8XF0VBGRZlNta2Dasu1E+nozs18cZtNjULYTa/waim2eLM0tJjHYmx6qjSJyGTltA32l+936PeRV1vDphD54H1qM2fFnrD5zOBbYk893FhPg6cYIXU1LRFqZ2et3c+BoLYsm9MWz4EvMzr9g9X2ROv9kPttZTDtvN66NPP3DjyIil5LTngf6SvZNXhEf7DnIU/3iiPWsxKyaCp1HYxJ+zpKcco4ea2R810DcXdU8i0jr8VXuYf6dWcDs/vHEuJdhVt0N4eMw3R9mcU4ZtfWNjI/RlQZF5PLTEegWpqT2GE+kZjAiPJjJsR0wS8aA5Yo1+G22F9eReaSWcTEBBHrpWycirUdRdR1PrtnFdREh/Cg69Pi58F09sQa9RVpxDVllddzQNZAA1UYRaQaqNC2IMYYZqbswBp4b2A0r/QXMwW+xRn/FEdqxYn8RPUK8iWvn7eioIiLNxhjDE2t2YVnw7MBuWNv/iClcgXX9EkqMPyvyikkO9SEm0MvRUUWkldASjhZk4d5CFucV8ezAbgRXbsF8NxOSfkVDhxH8N/sIvh6uDAtv+umkRESc2YeZBXy7v5jfD+xOu/KNmK1PQ/J0GtoP48usMvw93Rii2igizUgNdAtxoKqGp9fvZmJMB67r5I1ZcTsE98Hq9TSp+ZWU1NQzNjoQd1d9y0Sk9cirrOGZ9XuYFNuRER08MSunQEh/rJSnWLW/giO19YyNCdC6ZxFpVlrC0QI0GsPjq3bi6+HGU/3jMRunQe0hrNFfk1fZyHeHjjIk3FfnNBWRVqWh8XhtDPRy58m+cZgND0BtCdb1S8mtqGfr4WqGRfgR4qPaKCLNS4czW4B5O/ezrrCMOYMS8C1aCrv/gdXneWq8o/gqu4xIPw96tW/j6JgiIs3qzR15bDpUxguDE2h7aDFkvoXV9wWqvSJYnFNOlL8nKaE637OIND810A6250gVczZnMTUhnAFBLpjV90GnEZi4/2FxTjkGGB0doPM9i0irklFayZ+/y+LeHhH0CTSY1J8dP51n7D0szikHYHQXnQtfRBxDDbQD1Tc28tiqnUS09eLxq2Mw6x8FWyXWwDdIK6ohp7yO67r408bd1dFRRUSaja2xkcdX7STKz4f/7RmNWfcINNRgDXydrYdryC2vY3S0Pz6qjSLiIFoD7UBv79zPjpJKPh7fG88Dn2Oy3sEa9CZlbp1Yub+I5FAfogN0WiYRaV3eTM9j15EqFo7vg0f+p5js97GGzKfUJZTV+cX0bO9DlL9qo4g4jo5AO0heZQ1//i6bOxPCSfa1Ydbcf/yKWjF38HVOGW09XBkcptMyiUjrkltRzctbc7g7IYKkNrWYNQ9CxI2YLj9lcU45fp6uDFJtFBEHUwPtAMYYnlyzi3ZeHvyyVzRm7cNg6rGu+TtbDldTUGXjui4BulS3iLQqxhhmrNlFqLcHP0/pgln7IADWNa+y+VA1hUeP10adsk5EHE1LOBzg06xCVh8s5a1Ryfjs/xiT+xHW0PcocwkhNb+Inu19CPP1cHRMEZFm9e+9BawtOML861Lw2b8As+8TrOELKKUdaw8Uc3X7NnRqq9ooIo6nI9DNrLjmGL/bsIcbotszNJjjH46J+hEmarJ96cbAzvrzpIi0LkXVdTy3IZOJMR0YFNR4/EPVXX6CifyRfenGNWG+jo4pIgKogW52z2zYA3D8ogAbf3V86Ua/v7BVSzdEpBX77fo9uFoW0/vGYjb8EgCr/8tsLjzKIS3dEJEWRks4mtHy/GL+k32IFwYnEFS2BrN3PtY1r1LmEsxqLd0QkVZqaV4RX+Qe5sUhiQSWrMBkvYc18HVKTQBrDxTTq4OWbohIy6Ij0M3kqK2eJ9fsZlCndtwUGYBZ8wCEDsTE3qOlGyLSalUeq+eptbsZ2jmICRF+mLUPQYehmK53/f/Sjc5auiEiLYsa6Gby4pZsSmuP8cw13SD9j1CVgzXwNbYertHSDRFptf78XRblx+p55pp4SHsOju63n3VDSzdEpKVSA90Mdh+pYv7OfB5J6UJ4Qy6k/QGSfkWldxypB6pICdXSDRFpfTJKK3lnVz6P9uxCJ1sWbJ+DddUTlHvGsPZAJb101g0RaaG0BvoyM8bw23V7iPD1ZmpCGGbxSGgTiXXVdJblVODlajFQnywXkVbGGMPsdXuI9vPhju6dMV8NB78YTNKvWJZdgY+7KwM6t3V0TBGRM1IDfZl9mXuYtYVHmDcqBY/s+ZhDq7BGLyarErLL6hgfE4CHq/4QICKty3+yD7HxUBn/HN0T98w3MIfXYI35lswKQ255HTfEBuKu2igiLZSq02VUbWvg2Y2ZjAwPZkhQA2bjryFmCrb2w1m2r5wu/p50DfRydEwRkWZ11FbPHzbtZXRkCAMD6jCbp0PsVI6FDGZ5XgUxAZ7EBKg2ikjLpSPQl9Fr23MpqbXxZL84zMb7wXLB6vsC6w5UUdvQyPBIPyxLH44Rkdblb9tyKauzMaNPLGbDPeDigdXnedYeqORYg2FYhL+jI4qInJPTHoGeMWMGSUlJpKSkkJKSwoIFC+xjjY2NTJs2jZiYGLp27corr7zS7PnyKmv4R3oeP+sRQXjNVsh6F6v3cxQ3+vPdoaP06+SLv6f+/yIirUt2eTVv7cjj/qRIOh/dCDkLsHr/gcP1vmw9VM2Azm3x83R1dEwRkXNy2g7u8ccf59lnnwXg4MGDdOvWjeuuu47AwEDeffdddu7cyZ49eygvL6dXr15ce+21dOvWrdnyPbN+D0Fe7tyfFIH56qcQ1AsTcxdLd5cR6OXG1e3bNFsWEZGWwBjD79bvob2PJz9LDMN8OQmC+2JiprB01xGCvN1ICVVtFJGWz2mPQAcEBNj/XVlZiWVZNDY2ArBgwQLuv/9+XF1dadeuHZMnT+aDDz5otmwr8kv4Zn8xM/rE4p37LhRvwur3Euklxyg4amNElB+uOq+piLQy3+4vZsWBEmb0jcUr520o3YrV70W2F9dy6KiNEVH+qo0i4hSc9gg0wF/+8hfmzp1Lfn4+b731FkFBQQDk5eURGRlp3y4qKopNmzadc1/l5eUYY+z3PT098fT0POv2J7Y9+TkAxxoa+e36PfTvEMDojp6YT2ZAl1uoDuzH6vRiEoK86dzW47TnNbez5Xcmzj4H5Xc8Z5/DxeZvyvaXqjbW1Tfwuw17GNQpkJGhbphPZkLM7RwN6M3q9GJ6BHvTsY27w78Xzv4zAc4/B+V3PGefQ3PUxhbbQA8ePJiMjIwzjm3ZsoXw8HAeeeQRHnnkEbZt28btt9/OyJEj7U30yR/Ou5AXJCIi4pT706dPZ8aMGed9Xm1t7Sn339p1gLzKal68JpaGLb/F1VZJXY/ZLM8twwJ6B7tRU1Nz3v02l+/nd0bOPgfldzxnn0NT8zelBl2q2viPnfkcrKrjb4O60fDdbFwbaqhLfIpluWW4AL2D3VUbLzFnn4PyO56zz+Fy1sYW20CvWrXqgrdNTk6mc+fOLF++nB/96EdERESQm5tLnz59ANi3b99pbwLfl5eXh5+fn/3+hRxlqa2txcvLy96sH6m18frOg/w0Powe3iWw5xVInkG5ZwR7yku5NtKPdr4+Fzyvy+lM+Z2Ns89B+R3P2edwsfltNtsFb3spamNxzTHezDjIlO5hJHgegsy/Q8/ZHPEIY295KaOi/Ahoq9p4qTj7HJTf8Zx9Ds1RG1tsA30+GRkZdO/eHYCsrCy2bNlCQkICAJMmTeLvf/87EydOpLy8nAULFvDVV1+dc3/+/v6nvElcKMuy7N+cV7bnYjA80rMLrP4x+HSCxF+yIrOSYG83kkJ8WtwP4sn5nZWzz0H5Hc/Z59DU/E3Z9lLUxrlpubi4WDx0VRSsuhHaRkDCo6zMrCTUx43EYNXGy8HZ56D8jufsc7ictdFpG+jf/OY37N27F3d3d9zc3Pjb3/5mb6inTJnCxo0biYuLA46fsePE2OWSX1nDuxn5PJzchXbF32L2f4E1fAF7KiwKjtr4UXw7XJz4h1BE5GLkVlTzr10H+EWvaAKLF2MOLMa69mN2lRsKj9qY1K2dU79Bi0jr5LQN9KJFi8465urqyty5c5sxDfzpu2wCPN2Z2r0j5otx0GEYDeETWZVeREyAJxF+Z/+Tp4jIlepPm7MI9vbgrvj2mM9GQ6cR1HeewOr0YmIDvQjzVW0UEefjtKexa0l2lFSyKLuQn/eMxmfv36EiE6vfi2w+dJSjtkYGhzf9z58iIs5uW1EFX+Qe5tGe0XjteQWqcrD6vsimQ9XU1DcyONzX0RFFRC6KGuhL4I+b9hLj78OPI9tgtj0DsXdztG0iGwqO0rN9GwK9nPZAv4jIRTHG8MdNmcQFtGFiuBcm7TmI/xlVPt3YVFhFr/ZtdDVWEXFaaqB/oNUHS1l9sJTHr+6K644/QX0NVs9ZpOZX4u5i0a9jW0dHFBFpdssPlLCusIzHe3fFZcccaDyGlfwkq/Mr8HBxoU8n1UYRcV5qoH+ARmN4ftNerg71Z2SwDXa+BAmPcMgEs7Okhms6t8XTTS+xiLQuDY2G5zdl0bd9AMMDa2HnXyDxFxQ2tmNXaS0Dw3zxdFVtFBHnpQr2A/x3XzE7Sqv4Te+ukPYsuHhAj8dYkVdBsLcbPUJaxnlNRUSa02f7ithTdpTf9OkK254BtzaQ+L8sz6sgxMeNhGBvR0cUEflB1EBfpLqGRv6avp/rIkLo5V0Mu9/Auuo3ZFZ7c7DKxtAIP522TkRandr6Bv66fT9jo0JJ9iiEzLewrnqCPVWeFB61MTRctVFEnJ8+wXGR3tuVT2F1HW9fF4357mfgHUpj/EOk7qqki79OWycirdP8jHxKam38slc05rup4NOJxvj7Sd1ZSUyAJ+GqjSJyBdAR6IuUHOLPL5Mj6VK/G3IWYKU8xY4yQ3ldAwPDdGomEWmdrg7155cpkUQe2wG5H2P1nM320kYqjqk2isiVQ0egL9LVof4k+HrAqmngH09DzF2sSy+lWzsvQnzcHR1PRMQhercPINHPE1beDwEJ1Efdxvr0UroHeRPkrdooIlcGNdA/gMuhlXBwMdbwD9laVEdNfSMDOusIi4i0bi6F30LBt1jXLmRrcR21DY0M6KzT1onIlUNLOC6SMQa3tFkQ3JtjYTexsaCKxGAfAnTRFBFpxYxpPF4bQ/pT12k8GwuqSArx0UVTROSKoop2sfZ9gkvpJhi9mO8OVWNrNPTXhQFEpLXL/TcuR7bC9d+y+VA1DcbQVxeUEpErjI5AXyxbBQ3hE6kJHsbmwqOkhLahrYero1OJiDiWrZKGiElUBw1ky6Gj9FRtFJErkI5AXyQr9i5sYT9hU2EVFtBHR1hERLDi7sEW/lM2FhzFxYLeqo0icgXSEegfoMrWyLbD1fTq0AZvd72UIiIAlcca2V5UTe8ObfFyU20UkSuPKtsPsLmoFndXi14d2jg6iohIi7G5qBZPVxdS2vs4OoqIyGWhBvoiHamtZ/cRG307tMXTVS+jiAhAaW09e8ps9OnYBg/VRhG5Qqm6XaTtRdX4uFlcFaojLCIiJ6QdrqaNu0VSiGqjiFy59CHCizQozJc4XxfcXCxHRxERaTGGhPvSzU+1UUSubDoCfZFcLAt/T52aSUTkZKqNItIaqIEWEREREWkCNdAiIiIiIk2gBlpEREREpAnUQIuIiIiINIHTNtAzZswgKSmJlJQUUlJSWLBggX1s9uzZhIaG2sduu+22S/716+rqePbZZ6mrq7vk+24Ozp4fnH8Oyu94zj6Hlpi/JWZqCmfPD84/B+V3PGefQ3Pkt4wx5rLt/TIqKysjICAAgIMHD9KtWzf27dtHYGAgs2fPpqqqihdeeOG8+6moqMDf35/y8nL8/Pwu+OuXl5cTEBBAWVkZ/v7+FzsNh3H2/OD8c1B+x3P2OVxs/gupe6qNzpkfnH8Oyu94zj6Hy1kbT3DaI9AnmmeAyspKLMuisbHRcYFEREREpFVw6gup/OUvf2Hu3Lnk5+fz1ltvERQUZB97//33WbJkCUFBQcycOZPhw4efc1/l5eWcfDDe09MTT0/Pcz7H19cXACc9iO/0+cH556D8jufsc7iY/E3ZVrXROTn7HJTf8Zx9Dpe7NrbYJRyDBw8mIyPjjGNbtmwhPDzcfn/btm3cfvvtLF++nKCgIAoLCwkKCsLd3Z3U1FRuvvlmNm7cSGRk5Gn7OnGY//t++ctf8thjj12y+YiItBSVlZXExcWd88+bqo0i0tpcSG20M1eI0aNHm3//+99nHLvuuuvOOrZ//34D6Kabbrq1utv+/fvPWlNVG3XTTbfWejtXbTzBaZdwZGRk0L17dwCysrLYsmULCQkJAOTn5xMWFgZAZmYmW7duJSkp6Yz76dSpE/v378fX1xfLsponvIiIAxljqKyspFOnTmfdRrVRRFqbC6mNJ7TYJRznc+ONN7J3717c3d1xc3Pj17/+NZMmTQLgzjvvZPPmzbi5ueHq6soTTzzBj3/8YwcnFhEREZErgdM20CIiIiIijuC0p7ETEREREXEENdBnUVdXx8MPP0xsbCyJiYncfvvtABw+fJjrr7+e2NhYevTowerVq+3Pqa6u5tZbb6Vr167ExcWxcOFCR8UH4Ouvv+bqq6+mZ8+e9OjRg/nz5wMtdw6PPPIIUVFRWJZFenq6/fGLzdvY2Mi0adOIiYmha9euvPLKKw7Jf/fddxMfH09KSgpDhgxh69atTpX/hPnz52NZFp9//nmLzH+uORhjmD17NnFxcfTo0YNhw4a1yDmcLf+mTZsYMGAAPXv2pHv37jz//PMOy6/aqNp4qfKrNjo2v7PUxXPNwaG18bwfM2ylHn30UTNt2jTT2NhojDHm4MGDxhhjpk6dambNmmWMMWbDhg0mIiLC2Gw2Y4wxTz/9tLnzzjuNMcZkZ2eb9u3bm9LS0mbPbowxjY2Npl27dmbbtm3GGGNycnKMp6enqaioaLFzWLFihdm/f7+JjIw027dvtz9+sXnnz59vrr32WlNfX29KSkpMZGSkycjIaPb8ixYtsuf97LPPTGxsrH3MGfIbc/yMDAMGDDD9+/c3n332WYvMf645vPTSS2bixImmrq7OGPP/v88tbQ5ny5+SkmIWLVpkjDGmpKTEhISEmB07djgkv2qjauOlyq/a6Nj8zlIXzzUHR9ZGNdBnUFVVZfz9/U1lZeVpY23atDGHDx+23+/Tp49ZtmyZMcaYhIQEs2HDBvvYpEmTzLx58y533DM68SaxYsUKY4wx27ZtM506dTJ1dXUtfg7f/wW52Lxjx441H374oX3s8ccft7/ZXE5nKrInFBUVGQ8PD9PQ0GCMcZ78Y8aMMevWrTNDhw495U2iJeY35vQ5dO7c2WRmZp5x25Y4hzO9ScyfP98YY0xeXp7p3LmzKSgoaPb8qo2qjT+EauM8Y0zLye9sddGYllUbtYTjDLKysggKCuKZZ56hd+/eDB48mKVLl1JSUkJjYyMhISH2baOiosjLywMgLy/vlIu1nDzW3CzL4sMPP2TixIlERkYyaNAg5s+fT2VlpdPMAfhBr3lLmwvAyy+/zNixY3FxOf6r5wz5X331VRITE+nXr99pY86Qv6KigqKiIj755BP69+9P//79WbBggX3cGeYwb948Zs6cSUREBHFxcfz+97+nQ4cO5814qfOrNraMOYBqY0vI78y18Uqoi+DY2ui054G+nGw2G9nZ2SQkJPCHP/yBbdu2MXLkSNLT0087H6r53klMTh7//lhzqq+v5/e//z2LFi1i4MCBbNy4kZtuuom0tDSnmcMJPyRvS5rLu+++y4cffsiqVatOebwl58/JyeH1118nNTX1rNu05Pxw/Pf52LFj1NTUsG7dOvLy8hgwYACJiYn06NHjvDlbwhzmzJnDnDlzmDx5MtnZ2QwbNoy+ffsSHx9/3oyXMr9qY8uYwwmqjWceaw7OXhuvhLoIjq2NOgJ9BpGRkbi4uHDbbbcBkJycTJcuXeyXFi8qKrJvu2/fPiIiIgCIiIggNzf3jGPNbevWrRw8eJCBAwcC0KdPHzp16kRaWhrgHHMACAoKAi4ub0uay4IFC3j66adZsmQJoaGh9sdbev61a9dy8OBBunfvTlRUFOvWreOee+7h9ddfP2/GlpAfjv8MtW3b1v5ht4iICAYOHMimTZvOm7MlzKG4uJhPPvmEyZMnAxAdHU2/fv1Ys2bNeTNe6vyqjblnHHME1cYzjzUXZ6+Nzl4XoQXUxiYt+GhFRo0aZb744gtjjDG5ubkmODjYHDx40Nx5552nfGgjPDzc/iGIWbNmnbJgPTQ01JSUlDgiviksLDS+vr5m165dxhhjMjMzTWBgoMnPz2/xc/j+GqeLzTtv3jwzYsQI+4cEIiIizM6dO5s9/4IFC0zXrl1Nbm7uads6Q/6TfX+dX0vMb8zpc7jvvvvM3LlzjTHGlJaWmsjISLN58+YWO4eT89fX15vAwECzfPlyY8zxtaJhYWH2tX3NnV+1UbXxUuVXbXRsfmeri9+fg6Nroxros8jKyjJDhw41PXr0MMnJyWbhwoXGmOPFd9SoUaZr164mISHB/o0z5vgHbCZPnmxiYmJMbGys+eijjxwV3xhjzL/+9S/To0cPc9VVV5mkpCTz/vvvG2Na7hwefPBB07lzZ+Pq6mrat29vYmJiflDe+vp68+CDD5ro6GgTHR1t/vrXvzokv5ubmwkLCzPJycn2W3FxsdPkP9n33yRaUv5zzaGoqMiMHz/eJCYmmsTERPPaa6+1yDmcLf+SJUtMr169zFVXXWW6d+9uXnrpJYflV21UbbxU+VUbHZvfWeriuebgyNqoKxGKiIiIiDSB1kCLiIiIiDSBGmgRERERkSZQAy0iIiIi0gRqoEVEREREmkANtIiIiIhIE6iBFhERERFpAjXQIiIiIiJNoAZaRERERKQJ1ECLiIiIiDSBGmgREWmScePG8e67717w9s8++ywPPPDAZUwkItK81ECLiNOzLOuct+XLlzs64iUTFhbG22+/7bCvn5qayo4dO7j11lsByM3NxbIs9u7da98mKyuLyMhIbr31Vurr63n44YdZsGABOTk5jootInJJqYEWEadXUFBgvz366KMMGDDglMeuueYaR0c8r2PHjjXb12psbKS+vv6invvKK69w22234erqesbx3bt3M2TIEIYPH857772Hm5sb/v7+XH/99bzxxhs/JLaISIuhBlpEnF6HDh3stzZt2uDh4WG/HxISwu9+9zvCwsLw9fVl2LBhpKWl2Z87e/ZsBg0axMsvv0zHjh0JCAjgueeeo66ujvvvvx8/Pz+6du3KkiVL7M95++23CQsL45133iE8PJy2bdty3333ndIEV1dX89BDDxESEkJAQADjx48nNzfXPn7XXXdx22238atf/Yrg4GAmTZoEwKOPPkp0dDQ+Pj4kJiayYMEC+3OGDRvGgQMHmDp1KpZlMWzYMPvjTz755CmvSVRUlL1hPXGU+KOPPqJPnz54eXmxfft2GhoamDlz5llfm++z2Wx8+umnjB079ozj6enpDB06lAkTJjBv3jxcXP7/LWbcuHF89NFHZ923iIgzUQMtIle0p59+mv/+97+8//77bNmyhYEDBzJq1CgqKirs26SlpbF161aWLVvGSy+9xIwZM7jhhhtITExk8+bNjB49mjvuuOOUBrmkpIQ333yTzz//nE8++YQvvviC5557zj5+//33k5mZyZdffsn69esJDQ1lwoQJNDQ02LdZtGgRdXV1pKam8sILLwAQFBTEBx98QHp6OtOmTWPKlCls374dgIULF9KxY0deeuklCgoKWLhwYZNei6eeeopnn32WnTt3Eh0dfUGvzcm2bdtGTU0NPXv2PG1sy5YtDBs2jFtuuYXXXnsNy7JOGe/duzeZmZkcPny4SZlFRFokIyJyBZkxY4YZOnSoMcaYmpoa4+3tbbZv337KNrGxseadd94xxhgza9YsExgYaGpra+3j8fHxZty4cfb7BQUFBjBpaWnGGGPmzZtnAJORkWHf5vXXXzdBQUHGGGNycnKMh4eHKS0ttY8fO3bM+Pj4mFWrVhljjLnzzjtNdHS0aWhoOOd8Ro8ebZ5++mn7/c6dO5t58+adss3QoUPNjBkzTnksMjLSvP766/Y8gHn77bft4xfy2nzfxx9/bPz9/U957MS+PTw8zC233HLWeZSXlxvAbNy48azbiIg4CzeHdu8iIpdRVlYWNTU19O/f/5THa2pqyM7Ott+PjY3F09PTfr99+/YkJiaech+gqKjI/pivry/dunWz3+/bty8lJSWUlJSwY8cObDYb4eHhZ/y6gwYNAiA5OfmUZQ4A8+fP569//Su5ubnU1tZSV1d32n4u1slHji/0tTlZXV3dKa/TycaNG8dnn31GamoqAwcOPG3c29vbvn8REWenBlpErlhVVVUALF++nICAgFPG2rVrZ/+3u7v7KWOWZZ3y2InlCI2Njac9drav6+3tzdatW08bCw0Ntf/bx8fnlLFVq1Zx3333MWfOHIYOHUrbtm2ZNm0aNpvtrF8LwMXFBWPMKY+d6Tknf70LfW1OFhQURHl5+RnHnn/+eTp06MDYsWP55ptv6NOnzynjpaWlAAQHB59zLiIizkANtIhcsbp3746HhwcFBQX07t37ku67oqKC3bt3Ex8fD8DGjRsJCgoiKCiI5ORkqqurqampISkp6YL3uX79ehISEvj5z38OHG/Ys7Ky7EfA4Xizf/I6aoCQkBAKCwvt94uKik65fyYX89okJydTV1dHTk4OXbp0OW187ty51NbWMnr0aJYtW0ZycrJ9bOfOnfj4+BAbG3tBX0tEpCXThwhF5Irl5+fHww8/zAMPPMDHH39MTk4Oa9euZfr06ezYseMH7dvLy4sHHniAbdu2sXTpUmbNmsVDDz0EQLdu3Zg4cSK33HILX3/9NTk5OaxcuZJp06ZRUlJy1n3GxMSwe/duPv/8c3bv3s20adNOa4QjIyNZuXIlhYWF9qPBQ4YMYeHChSxbtozt27dz7733nnWpxQ95bdq3b09SUhKpqalnHLcsizfeeIMxY8YwatQodu7caR9LTU1l6NChuLnpuI2IOD810CJyRZszZw4PPvggjz32GPHx8UyePJn9+/cTFBT0g/YbFBTEHXfcwdixY7nhhhsYPXo006dPt4+/9957XH/99dx9991069aNu+66C5vNdtqyjZPddNNN3HfffUyZMoVrrrkGX19fJkyYcMo2s2fPZv369YSHh3PjjTcCcO+99zJx4kRuvvlmxo0bx6233nrKUpGzuZjX5q677jrn6ehcXFz45z//yZAhQxg5ciSZmZkAfPTRR0ydOvW8mUREnIFlvr9wTkREzuntt9/mySefJD8/39FRml1lZSVxcXGkpqYSHR19Qc85sbY7PT1dR6BF5IqgI9AiInLBfH19eeutt5r0n4fy8nLefPNNNc8icsVQNRMRkSYZM2ZMk7YfP378ZUoiIuIYWsIhIiIiItIEWsIhIiIiItIEaqBFRERERJpADbSIiIiISBOogRYRERERaQI10CIiIiIiTaAGWkRERESkCdRAi4iIiIg0gRpoEREREZEmUAMtIiIiItIE/wcv9ksJEdvE9AAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "gp.layout(colors=gp.colors.hollywood)\n", "\n", "fig, axs = plt.subplots(1, 2, figsize=(180*mm, 85*mm), sharex=True, sharey=True)\n", "\n", "for pressure, df in fO2_QFM.groupby(\"P_bar\"):\n", "\n", " label = f\"{pressure} bar\" if pressure < 100 else f\"{pressure / 1e4} GPa\"\n", "\n", " axs[0].plot(df[\"T_K\"], np.log10(df[\"fO2\"]), \"-\", label=label)\n", " axs[1].plot(df[\"T_K\"], np.log10(df[\"fO2_magmapandas\"]), \"-\")\n", "\n", "axs[0].legend(frameon=True, fancybox=False)\n", "axs[0].set_title(\"Anenburg app\")\n", "axs[1].set_title(\"MagmaPandas\")\n", "\n", "fig.supxlabel(\"Temperature (K)\")\n", "fig.supylabel(\"log10 $f$O$_2$\")\n", "fig.suptitle(\"QFM\")\n", "\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "94d383d7", "metadata": {}, "source": [ "Note that there is a small difference in the pressure dependence visible at 10 GPa between the Anenburg app and MagmaPandas. This is most likely a result of the codes using different root-finding algorithms to solve the pressures of phase transitions. At the other calculated pressures, results are near identical." ] }, { "cell_type": "markdown", "id": "c6a22fcd", "metadata": {}, "source": [ "### IW" ] }, { "cell_type": "markdown", "id": "96ba16b0", "metadata": {}, "source": [ "fO2 at IW is calculated according to Hirschmann (2021). Here we compare Magmapandas results to those from the Matlabscript provided by the Hirschmann paper." ] }, { "cell_type": "code", "execution_count": 5, "id": "0519191b", "metadata": {}, "outputs": [], "source": [ "fO2_hirschmann_file = \"./data/fO2/fO2_IW_hirschmann.csv\"\n", "\n", "fO2_IW = pd.read_csv(fO2_hirschmann_file)\n", "T_K = fO2_IW[\"T_K\"]\n", "P_bar = fO2_IW[\"P_bar\"]" ] }, { "cell_type": "code", "execution_count": 6, "id": "4d4146cd", "metadata": {}, "outputs": [], "source": [ "fO2_IW[\"fO2_magmapandas\"] = IW.calculate_fO2(T_K=T_K, P_bar=P_bar, logshift=0)" ] }, { "cell_type": "code", "execution_count": 7, "id": "27859835", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAtAAAAFaCAYAAAA3lGq/AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjEsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvc2/+5QAAAAlwSFlzAAAPYQAAD2EBqD+naQAAmAFJREFUeJzs3Xd4FNX6wPHvbN/0DikklNAJHWlSFRAUCyoogoL92utPL1i4eu1eC4q9SxFs9CJSpFdpQugJISGd9J6d8/tjIYq0BEM2m7yf5/GR7MzOvmd3886bM2fO0ZRSCiGEEEIIIUSlGFwdgBBCCCGEEO5ECmghhBBCCCGqQApoIYQQQgghqkAKaCGEEEIIIapACmghhBBCCCGqQApoIYQQQgghqkAKaCGEEEIIIapACmghhBBCCCGqQApoIYQQQgghqkAKaCGEEEIIIapACmghhHCh/v3788wzz/Dwww/Tt2/fU7Z9/PHHaJrGTz/9dMrjrVu35r///W9NhimEEOIvpIAWQohaoG/fvmzatImSkpKKx1atWkVERASrVq2qeCw9PZ29e/eeVmwLIYSoOVJACyFELdC3b19KSkrYuHFjxWOrV6/m8ccfP6WAXr16NVarle7du7siTCGEEEgBLYQQtUJwcDCtW7euKJbj4+PJzc3lzjvvJDY2lpycHMDZK929e3esVqsrwxVCiHpNCmghhKgl+vbtW1FAr169mt69e+Pl5UWHDh1Yu3ZtxeP9+vVzZZhCCFHvSQEthBC1RL9+/Vi/fj3l5eWsWrWKSy+9FIBLL72UVatWkZuby44dO2T8sxBCuJgU0EIIUUv069eP/Px8fv/9d1avXk2fPn0A6NOnD6tWrWLdunVomkbPnj1dHKkQQtRvUkALIUQtERYWRrNmzfjhhx84cuQI3bp1A5w90Fu3bmXx4sV07doVT09PF0cqhBD1mxTQQghRi/Tt25cPP/yQbt26VdwoGBgYSNOmTfn8889l+IYQQtQCUkALIUQtcnIYx8nxzyf16dOH/Px8KaCFEKIW0JRSytVBCCGEEEII4S6kB1oIIYQQQogqkAJaCCGEEEKIKpACWgghhBBCiCqQAloIIYQQQogqkAJaCCGEEEKIKpACWgghhBBCiCqQAloIIYQQQogqkAJaCCGEEEKIKpACWgghhBBCiCqQAloIIYQQQogqkAJauIXGjRvzxx9/nPJY//79mT9/PgDPPfccM2fOrLbX0zSN/Pz8ajtebZadnc3rr79+ymN/fW/P5UyfixDCPTRu3JiQkBDKysoqHlu+fDmapvHEE0+4MDInTdNo3749HTt2pH379nz//ffVfvz6kudF9ZMCWtQJL7zwAqNGjTrtcV3X0XXdBRG5jzMV0EKI+iEyMpK5c+dW/PzFF1/QtWtXF0Z0qnXr1rF9+3a+/vprbrvtNjIyMlwdkhCAFNCijhg3bhzvv/8+AJMmTWLs2LGMGDGCjh07kpyczAMPPECrVq3o0KEDXbp0obi4GIAFCxbQrVs3OnToQMeOHdm4cWPFMadMmUL37t1p0qQJX375ZcXjjRs35rnnnqNXr15ERkYydepU3n33XS655BKaNWvGypUrASgvL2fIkCF07dqVtm3bcsstt1BYWAjAV199xZAhQ7j55puJiYmha9euHD58GICVK1fSsWNH7rvvPjp06EDbtm3ZsmXLWdt97733ctlllxEVFcXDDz/MihUr6Nu3L40bN+att96q2PfJJ5+kW7dudOzYkX79+nHgwAEA7r33XrKzs+nYseMZT5zTp0+ne/fudOrUiY4dO7Jw4cJTtk+bNo2+ffsSHR19yusJIWq/22+/nS+++AKAnJwcNmzYwBVXXAHArl276NOnD507d6ZNmza88sorFc9LSkrisssuo23btlx11VVcddVVFTm4OvLS33Xq1AkvLy/i4+PP+pz4+HiCgoJ47rnn6NKlC9HR0afkq59++olWrVrRs2dPXnzxxVOOP2bMGLp27Ur79u256qqrSEtLA+DAgQP07t2bDh06EBMTwzPPPPNP33JRVygh3EBUVJRq2bKl6tChQ8V/np6eat68eUoppW677Tb13nvvKaWUev7551V4eLhKTU1VSin1+++/q1atWimHw6GUUio7O1s5HA61b98+1aBBA7Vv3z6llFKlpaUqOztbKaUUoN555x2llFJ79uxRXl5eqqysrCKWJ554Qiml1KZNm5TdbldTpkxRSik1c+ZM1bNnT6WUUrquq4yMjIp/33vvveqNN95QSin15ZdfKl9fXxUfH6+UUuqpp55Sd999t1JKqRUrViiTyaQ2b96slFLqww8/VIMHDz7j+3Lbbbep3r17q+LiYlVQUKCCg4PV+PHjlcPhUImJicrT01Pl5eUppZRKT0+veN6MGTPUlVdeqZRSKi4uTgUGBp5y3H79+lW8txkZGUrX9Yp9Q0NDVWlpacV7MX78+IrjR0ZGqg0bNpzjkxRC1BZRUVFq165dqlWrVioxMVF9+OGH6umnn1bPP/+8evzxx1Vubq4qLi5WSilVWFioOnbsWJGXRowYoV588UWllFJHjhxR3t7eFTm4OvKSUs48fHK/pUuXKh8fH5WdnX3OXAao2bNnK6WUWrRokWrRooVSSqnU1FQVEBCg9u7dq5RS6rXXXjvl+H895iuvvKLuv/9+pZRSDz30kHrppZcqtmVmZv6Dd1zUJSbXlu9CVN4PP/xAu3btKn7u37//Wfe96qqrCAkJAaBp06aUlZVx++23M2DAAK688koMBgNLly5l2LBhtGjRAgCz2Yyvr2/FMW655RYAWrdujclkIiUlhYiICICK4SKdO3emqKiIkSNHAtClS5eKnmSlFG+//TYLFiygvLycnJwc+vbtW3H8Sy+9lKioKAB69uzJe++9V7GtZcuWFb3BPXv25M033zxrW6+99lqsVmvF84YNG4bBYCA8PBx/f38SExNp1aoVv/zyC++99x55eXnouk5ubu5Zj/lXcXFx3HLLLSQmJmIymcjIyODIkSNER0cDcMcddwAQFBTEddddx7Jly+jevXulji2EcL2xY8fy9ddfM3v2bKZNm8a0adMAKCoq4r777mP79u0YDAaOHj3K9u3b6dq1KytWrGDy5MmAcxjIZZdddsoxqysv9erVC4PBgL+/P3PmzMHX15fp06ef9Tmenp5cc801gDN3Hjp0CIANGzbQuXNnWrZsCcDdd9/NU089VfG8adOm8e2331JSUkJRURENGzYEoG/fvjz55JMUFBTQr18/Lr/88mp734V7kwJa1EleXl4V//b19WX37t389ttvrFixgn//+9+sWrXqvMew2WwV/zYajZSXl5+2zWg0nvbzyf2mT5/Ob7/9xqpVq/D29mby5MmnvG5ljn+mbeeL80zPTUhI4KGHHmLTpk00bdqUnTt3MnDgwPO+BwA33XQTb775Jtdeey0AAQEBFUNgzkTTtEodVwhRO4wbN47OnTvTokULmjdvXvH4hAkTaNCgAdu2bcNkMjFixIhTfvfP9bteXXlp3bp1p+Tz8z3n76/jcDgAZ4fG2axZs4b333+fdevWERwczNy5c3nhhRcAuP766+nVqxdLly7l/fff55133jltGJuon2QMtKjz0tPTKSgoYPDgwbz88ss0btyYPXv2MGTIEBYtWsT+/fsBKCsrIycnp9peNysri8DAQLy9vcnLy+Orr76qtmNXVU5ODhaLhYYNG6KUqhirCODj40NhYeFZi/SsrCwaN24MwNSpU8nKyjpl+8nx4cePH2f27Nmn9UQJIWq3sLAwXnnlFV577bVTHs/KyiIiIgKTycS+fftYunRpxbb+/ftX5LSjR4+yfPnyKr/uufJSdT4HnL3R27Ztq8j3n332WcW2rKwsfHx8CAgIoLS0lI8//rhi24EDBwgJCeHWW2/l9ddfZ8OGDVVspairpAda1HlHjx7lrrvuoqysDF3X6dWrF0OHDsVsNvP5559z8803U1ZWhtFo5OOPP+aSSy6plte99dZbmTNnDm3atCE8PJw+ffqQlJRULceuqpiYGG688Ubatm1LZGQkgwYNqtgWEBDALbfcQkxMDJ6enqfdsPjuu+9y3XXXER4eTs+ePYmMjDxle1RUFH369CE5OZmHHnqo2t4/IUTNGT9+/GmPPfPMM4wdO5Zp06bRuHHjU3p63333XW699VZmzpxJixYt6N279ylD4CrjXHmpOp8DEBISwieffMLw4cMJDAzkhhtuqNg2dOhQpk6dSqtWrYiIiKBXr14sWbIEgO+//55p06ZhsVhQSvHRRx9VqY2i7tLUua5rCCGEEEL8TVFREWazGZPJRHJyMt26dWPZsmUVY4yFqOukB1oIIYQQVXLgwAFuvfVWlFKUlZXx/PPPS/Es6hXpgRZCCCGEEKIK5CZCIYQQQgghqkAKaCGEEEIIIapACmghhBBCCCGqQApoIYQQQgghqqDez8Kh6zrHjh3D29tbVlATQtQLSiny8vIICwvDYDhzP4rkRiFEfVOZ3HhSvS+gjx07RqNGjVwdhhBC1LijR48SERFxxm2SG4UQ9dW5cuNJ9b6A9vb2Bpxvlo+PT6Wfp5SiqKgIu93ulr0z7h4/uH8bJH7Xc/c2XGj8ubm5NGrUqCL/nYnkRveMH9y/DRK/67l7Gy5mbjyp3hfQJ99YHx+fKp8kzGazW3+53Dl+cP82SPyu5+5t+Kfxn+s5khvdM35w/zZI/K7n7m24mLnxJLmJUAghhBBCiCqo9z3QZ1JWVobD4TjnPkopSkpK0DTNbf86q0z8RqMRs9lcg5EJIYQQQtRudbKAHjduHL/++itBQUEADBo0iDfeeKNSzy0rK+PQoUNUZoVzXdfPe5dmbVaZ+DVNo1mzZlJECyGEEEKcUCcLaICnn36aBx54oMrPczgcKKUICwvDarWedT+lVEUB6q490OeLv6SkhGPHjuFwOKSAFkIIIYQ4oc4W0P+U1WrFZrOddXt9KKCFEEIIIcTp6mwB/dZbb/HJJ58QGRnJf//7Xzp27HjO/XNyclBKUVpaiq7r6LpeqWEcQKX3O5O8vDwGDRrE7t27Wb9+Pe3atTtl+1dffUV+fv4F9aZX1tniP/m4UuoftfFi+Wt87kjidz13b8OFxl+V/U/mxpOsVut5r85dSEy1hbvHD+7fBonf9dy9DTWRG92ygO7Tpw+xsbFn3LZt2zZeeuklQkNDMRgM/PzzzwwdOpQDBw7g5eV11mNGRkYCEBUVxUcffYSXl9d5J9EG5zjif8JqtTJnzhyeeuqpisL9r04Wr1V5naqMzT7XcU/GU1xcXKt/iYqLi10dwj8i8bueu7ehqvEXFRVVet+TufGkCRMmMHHixGqPqbZx9/jB/dsg8bueu7fhYuZGtyygV69eXel9r7vuOp5++mn27dtHly5dzrpfQkICPj4+lJaWkpaWRnBw8HmL0Oq4idBisdCgQQM0TcNgMJx2PE3TWLduHYsXLyY3N5cZM2YQERHBmDFjSEpKory8nGnTphEZGUmXLl3o2bMnubm5fPPNN+d97fPFfzIem812zuEsrqKUori4GJvN5pbDUCR+13P3NiilKC4qwlbFuU7Lysoqve/J3HhSZXqg3f49deP4wf3bIPG7nru3QSlFcWE+Ng+vi5Yb3bKAPp/ExMSK3uMNGzaQmZlJdHT0OZ/j6+uLj48PxcXFZGRkYDAYKNfheHH5mZ+gFLrSMWgGqMSHE2AzYTaee78zTSl3srBeuHAhS5cu5fXXX+e9997j008/xcPDg7lz5/LJJ5/w0ksvkZWVxaOPPnretjrD/7NH+WxfrpOP1/ap+mp7fOcj8bueu7RBFaShUrdB6nZUyjZU6jaMza9HG/hSleKvyr4nc2NVuct7ejbuHj+4fxskftdzhzYopUPWoYqcSOp2Z25sd/tFzY11soAeN24cqampGI1G7HY733//Pb6+vlU+zvHicqbvyaiWmEa3CaKB54XNZHGy5/ySSy7hrbfewuFw8NRTT7F9+3ZKSkpo27YtAP7+/pUqnoUQtZtSCnKPolJ+R6X8Dqm/o1K2Qf4x5w4Wb2jQEa3JYPTwXq4NVgghaojSyyFjLyr197/kx+1QmufcwTscrUFHtA53okcMvKix1MkC+tdff62W4wTYTIxuE3TmjRfQA32htm3bBsCWLVuIjo5m+/btpKamsnr1aubOnctPP/0E4NZzUgtRXymlIOcIKmUrKnkLpGx1nhSKMp07eISghXZBa38bWoNOaA06gX9TNM3gvEeiCmP2hBDCXShHGWTsqciNKuV3SNsB5SfGNQc0R2vQGS36SmjQyVk4e4Y4n1sDubFOFtDVxWzUztprXJ3TwA0bNozt27ezb98+7rnnHsaNG3fK9tLSUq644gry8/P57rvv8Pf3Jzk5mUGDBtG6det/9NpCiJqjlIK8JFTyZmexnLwFlbL1z2LZOxytYWcMXR+Ehp3RGnYGr7BafwlVCCH+CaU7IHNvRW5UyVv+UixrENjK2ZHQ5iZnXmzQEc1a9aFl1UkK6Fpg4cKFZ902bty40wpqOPONlFu2bKnOsIQQ/5AqzEQlb4bkzc4Tw7HNUJDi3OgVihba1Vksh3ZFa9gFzauhawMWQoiLrOKq27FNJwrmzZCyFUrznTsEtEQL7XqiWO4CDTuhWc4+i5qrSAEthBDVQJWXOG9eObbxxH+bIOugc6PNDy20G1qH253/D+uG5h3u2oCFEKIGqOJsZ5Gc9JfcWJjm3OgT6cyJvZ+B0G7OjgRb1e9ZcwUpoIUQoooqelCSNqCObUAlrYeUbaCXgdHiHI/XbBha2CVoYZeAf7QMwxBC1HlKdzjHLSetd+bHpPWQude50errzImd7nb+P/QSNK8Grg34H5ACWgghzkOVF0PyVlTSOlTiOlTShj+HYvg3QwvrgdZ2DFp4Dwhpj2Y6+zzJQghRV6jibGexnLgOEtehkjc5h2JoRmcujBqA1vNpZ24MaI6m1Z3JDqSAFkKIv1H5qajEtRUFMylbwVEKZg9nr0mH8WjhPdHCulfc9S2EEHWZUso533LiGmdHQuI6yNjt3OgRjBbeE0OviWjhPZ33dVg8XRvwRSYFtBCiflMKlbnPWTAnrkEdXfPn2GXfKLTwXmhtR6NF9HL2qBgkbQoh6gG93DkjRuJa1NE1qMQ1UJAKaBDcFi2iN1qPJ9EiejuvxNWzYWpyJqgFnnjiCTZu3EhkZCRffvklFoulYtvKlSsZO3YszZo1w2g0smzZstOe/8knnzB16lQ0TcNut/PGG28QExND8+bNCQ8PJy8vj0mTJjF8+PCabJYQtZLSHZC2A5WwCj1hFeajq9GLMkAzQEgHtGZD0SIuRWvUW270E0LUG6q82Hmj39ETuTFpPXpZARitENYdrcMdaI0udV59s/m5OlyXkwLaxbZt20ZKSgqrV6/mpZde4ocffmD06NGn7DNq1CjefPPNMz7/119/5bfffmP58uWYTCYyMjJIS3Pe3err68vKlStJSUlh4MCBUkCLekk5ypyLkxxZiTq6CpW4FkpynSeF8B7oMXdgbjIALaKny+cVFUKImqJKC5zD1I78hjq6Co5tdA5Vs/pCoz44uj+NpclA5/zLcl/HaaSAdrH169czePBgAK644gq+/PLL0wroH3/8kY0bN3LDDTfw8MMPn7Jt+vTpPP3005hMzo8yKCiIoKBTV09s2LAhxcXF7Nq1i4ceeoiSkhI6derE5MmTL2LLhHAN5ShzLlCSsNJZNCeuhbICMHuiNboUQ4+n0CL7Qmg3MFooLSrCYrfXu8uPQoj65c+CeSUqYSUc2wR6uXP8cmQ/tIFvOnNjcDvQDJQWFaFJbjwrKaDPQZUXQvbeM29Dga5QBg2oxJfLrxWayeO0h7OzswkLCwOcPcbHjx8/ZXvXrl3Zt28fANdccw2XXnopXbp0qdienJxc8fxXX32V2bNnc/XVVzNhwoSKffbt24e/vz/R0dEsX74cTdMYMWIEBw4coGXLluePXYhaTOkO5/zLR5aj4pejjq52FswWL7RGfTBc+hxaZD/nyn7GU1cWVUq5KGohhLi4VHkJJG1AP5EbObbROdWmRwhaVH+0QWPQovo7V/n7W5EsufH8pIA+l+y9qHndzrr55NetMl8zbfhmCOp82uP+/v7k5uY6Xy47m4CAgFO2e3n9ufrO1VdfzY4dO04poMPCwkhKSiIwMJCnn36aHj16MH/+fABycnLo378/BoOBKVOmEB8fz2OPPUZhYSFxcXEcO3ZMCmjhdpRSkLkPFf8rKn6ZsyelOBtMdmfB3PsZtMYDnQWz3PAnhKgnlO6AlN9P5Mblzpv+yovBHogW2R/t8redufEMBbOoOjm7nItfK2fhewYKhdIVmkFDq2QP9Jn06NGD//3vf9x6660sWbKE3r17n7I9NzcXHx/nuMzVq1dz7733nrJ99OjRvPbaa3z99deYTCbKy8srtp0cA33Sgw8+yIMPPsiwYcMYMWKE/IUp3IbKT3EWy3FLUfG/Ql4SGMwQ3gNDt0fQogZCeHc0o+X8BxNCiDqgYlq5+F+dufHIcmdngsULLbIfhn4vOQvmkPZ1av7l2kIK6HPQTB5n7DUGnFNf6TqawfCP/pLr1KkTDRs2pE+fPkRGRvLkk08CcM899/Dxxx8za9YsPvnkE0wmE71796Zv376nPP+yyy7j4MGDDBgwAKvVitVqZdKkSWd8reHDh/Poo4/y2WefnVJoC1HbqLIi55Ryh5egDv8C6bucG0I6oLW5Ca3x5WiN+tT5eUaFEOKvVHG2s3f58BJU3C+QE+9ctCS8B1rXhzE0GQRhl5w2XE1UP03V827I3NxcfH19ycnJwcfHh+LiYuLi4mjSpAk2m+2sz1NKoes6hn9YQLtKZeKv7HvhKkopioqKsLvpTQ4S/6nHIiMWdXgxKu4XVMJvzkuPXmFoTQejNRmEFnVZtS/7Wl8/g7/nvQvdpzpjqi3cPX5w/zZI/H85lu5w3hR9eAl63BJI2gjKAQEt0JqczI39q30Gofr6GVQl70kPtBDCJVRJnnNYxqFFqMOLITcBTDa0Rn2dlx6bDHZO1u+GyVsIIS6UKkhHxS05kRuXQFEmWH2dV96u+MBZOPs1dnWY9Z4U0EKIGlFx89+hBaiDC5wr/ullENAcrcW1zgVMIvuhme2uDrVWyy8rZ1taDlvSctiSms2lDXy4t2MzV4clhLhASumQ8jv6gfmoQwsheQugoEEntE73YGg2zHmPh9wUfU5phSVsSctma6ozNw4K9+f+ztEX7fXk0xBCXDSqvMS5eMmB+ahDCyDrkLOXOWoghsvfchbN/lL8nUtaYQmbU7PZkprNltQcYrPy0BX4mxx0thwjwsMEyHsohDtRJXnOG/9OdChQkOrsZW46BK3LfWhNr0DzaujqMGstpRRxuYUncqOzYD6SVwRAI2sJXYxxtCryAKSAFkK4CVWY4exhPjDXefmxrAB8GqFFX4U26Eq0qAFo5tPnRBfOk0JCXhGbU7PZlJrN5pQ/TwqR1hK6mg4z2riaruUbaVKeRIG9PTn2210ctRCiMlROAurAPGduPLLCeQUusDVau7EYoq+CiF5y899ZOHTFvqx8NqVmsykli82p2WQWl2FA0cpeQH8tli6G5XTRtxFcXkSOvSv59rEXNSYpoIUQ/5jK3I86MBf9wFxIXAtKQXh3DL0moDUfDsHtZCzzGTh7UYrYmJLFxpQsNqVkk1JYgga0tBXQ17iXrobldNO3ElTmPCkk+HRjp3UES+ydsdt9iQkw08jVDRFCnEYpBanb0ffPQR2YC6nbwGBCi+yP4bI30aKvlCtwZ1Gu6+w5ns/G5Cw2pmazOTWbvNJyLJqivT2bkYaddOVXOrMLe7mN4z7dORJwKSssD5Nuj8HHZiMmwEzkRYxRCmghRJUppVDJWzDu/h7H4fmQGetcyKTJILShnzhPDHL58TQnC+YNycfZkJLNxpQs0otKMWrQ1p7HVfxBN20pXdV2vMoNHPfoQbxvd1bYHiDN3h4/u50IbwttvS0M9rLgZTFQVFTk6mYJIU5QejkqYTXGPd+jH17gvDna6ovWbBhaj/9zDluz+bo6zFrHoSt2H89jQ3IWG1Ky2JKaTX6ZA5tB0dmeyR3G3+nGL3RUsRgcQaT79ORw4HXMtr7CcWtzgj3MhHtb6OptIczLgqf54udGKaBdLC8vj8svv5zdu3ezYcMG2rVrd8r28vJy7rzzTg4dOkTnzp159913TzvGJ598wtSpU9E0DbvdzhtvvEFMTAzNmzcnPDycvLw8Jk2axPDhw2uqWaIOOnliUPt/Ru2fDblHMdgC0JpfjTbgFWfxLEMzTqGU4mh+MeuTsypODKmFJRg1iLHncp22i+7aErqoHdjKbGT69iLOfxALbM+QYWtNkKeVCG8L3bythHtZsJsNpx1fCOFaqrwYFffridw4F4oyMHiFOW+ObjkCLbKvDM34G105h2SsT85ifXIWm070MNsNis72DO42baV72RJi9H2gh5Lm3ZtDQbcxw9adXEsUDbwshHtb6HOiYLaZaj43SgHtYna7nfnz51csoPJ38+bNIyIigq+++oq77rqLdevW0atXr4rtv/76K7/99hvLly/HZDKRkZFBWloa8OdKhCkpKQwcOFAKaFFlylGGSliJiv0Btf9nKEwH7wi0ltehtbiW4qCu2D29ZXjGX6QXlrAuOYt1ycdZl5xFUn4xBqCtRz5Xa7vpYVhEF33bnwVzwJXMsb1AhrUVIZ4WGnlb6eljIfwMJwUhRO2gyoqci5ns/R51YB6U5kFAS7SOdzpzo19b7B6ekhtPUEpxJK+IdceceXFDchbHS8qwGBSd7VncYfqdHqWLaK/Hgh5Kuk9vDgTfyVRbDwqsjWjgaaaRt4XLvK2EepmxGF2fG6WAPoeicgeHsgvOuE3BnysRVuJYzfw8sZuMpz1uMpkIDg4+6/PWr1/PVVddBcAVV1xxWgE9ffp0nn76aUwm50cZFBREUFDQKcdo2LAhxcXF7Nq1i4ceeoiSkhI6derE5MmTKxG5qG+Uo9Q5P/PeH5w9zUXHwa8pWvvxGFrdAKFd0TTN+Re+DB8gv6ycjSnZrD12nHXHjrP/RM5obi9mkDGWHsbFXOLYhFepgeO+vYjzH8wc6yQybK0I9rAQ4WOll7ezN0UKZiFqL1VWiDq40JkbD8533iAdHIOhx5NoLa9HC27j3E9yIwCZxaWsPXb8xH9ZHCsoxqhBe3seN5u206N0IZ31nRjLA0j37cPBoNuYZutJnjmSBl4WIrwtXO5jJayWFMx/JwX0ORzKLuDqeZur5Vhzh3ejXVDVVwrKzs6uWA3H19eX48ePn7I9OTmZsLAwAF599VVmz57N1VdfzYQJEyr22bdvH/7+/kRHR7N8+XI0TWPEiBEcOHCAli1b/oNWibpC6eWoIytRsTNR+35yFs0BzdE63essmht0lJ6UExy6YldmLquTjrPm2HG2peVQrhRh1nJ6WQ5zr3kZPctWElRSSI5vd+L9erPQ8uiJMcxWGnlb6eHjLJztUjALUaup8hLnCql7vnP2NJcVQINOzhukW12PFijn0JNKyh1sTsthTdJx1hzLZM/xfMDZmTDYuJeexkVc4tiIR6mZ436XcjhwBDMtb5JljSbIw0ykj5UBPs7OBGstLJj/zm0L6C+++IK3336b2NhY3nnnHR544IGKbbqu8/DDD7Nw4UI0TeOxxx7jvvvuq/JrNPPzZO7wbmfcdiE90BfC39+f3NxcwFlMBwQEnLI9LCyMpKQkAgMDefrpp+nRowfz588HICcnh/79+2MwGJgyZQrx8fE89thjFBYWEhcXx7Fjx6SArseU0uHoWvQ936H2/gCFac6e5k73YGg9CkLaS9F8wrH8YlYfy2RVkrM3Jbe0HC+joodHKhOtq+lTvIDIkmSKrB1JDL6UtdYPSLR3xWb1JNLHSgcfC5E+Vrwtp1+FEkLULspRhor/1Vk0758NJbnOnuZeE9Baj0QLuHhzC7sTpRSHcgr5LSmT1UmZbEzJpsShE2TW6W1LYLzlV3qX/kpwST65vt1J8LuUhZbHSLXH4GOzEunjHK7WyNuCh9n9cqPbFtBdunRh1qxZvPLKK6dtmzp1Knv27GH//v3k5OTQuXNnBg4cSKtWrar0GnaT8ay9xkopdF3HYDBc1CKjR48e/PLLL/Tt25clS5Zw++2nzvk6evRoXnvtNb7++mtMJhPl5eUV206OgT7pwQcf5MEHH2TYsGGMGDFCbkCqh5RSkLYTffd01J4ZkHsUfCLRYm7F0GYUNOwiRTPOnpSNqdmsSspkVWImB3MKMQAdPPO5zbyNPmVzaO/4A1UWRopff/banmeJrRcOSwAR3hYa+1jp62sh0GaS91MIN6CUgqT16H9MQ8XOgqIM55jmbo9iaDMKLai1q0OsFXJLy1l77LgzNyZlklxQgsUA3TyO86h1A30K59CiLI4Sa2uSgvqx0TqZBHt3TGYPGvlYaetjZaiPBT+b25afFdy2BR06dADAYDi9m3/mzJnce++9GI1GAgICGDlyJN999x2TJk2q4SgrZ9iwYWzfvp19+/Zxzz33MG7cOO655x4+/vhjhg8fzuzZs+nTpw+dOnWiZ8+epzz3sssu4+DBgwwYMACr1YrVaj1rO4cPH86jjz7KZ599dkqhLeo+lR2P2j0dffd0yNgN9kC01qMwtB0NET3RtNp/uexiS8grYmViBr8lZrI+OYtih05Dq4M+1jgesiymd+lyvIt1sv0v5VDgSKZZ3yPb0pQQTwuNfS0M97US6mnBaJCCWQh3oTJinUXznumQHQfe4WjtxzlzowxdQynF3qx8fkvMZGViJr+fGLLWzF7GENNe+pjmckn5BiwlnmT692N/8EOstl1KkSWUUC9nZ0J3XwshHmYMdey9dNsC+lwSEhKIioqq+Llx48Zs2bLlnM/JyclBKUVpaSm6rqPreqV7aP9pT+6CBQtOO95HH32EUgqj0chXX311zte6++67ufvuu087xubNm0/Zf9CgQfzxxx8VP+u6fs74Tz6ulKqVvdV/jc8dXez4VXEOat8PqF3fwtFVYPZAa34N2oBXocngU6ZVupAY3P39L3XobE7N4tf4NNak5HA4txCTBl09c3jYtpF+BT/SvCSOEksbEoMHsNLyGUft3bBa7UT5WOjh47wE6eHCqeUu9DOoyv4nc+NJJ/9Qr+6Yagt3jx/cvw0XPTcWpjuHZ+z6BlK2gs3PeRPgsFsgsu8pHQr1MTcWljlYcyyTZUfSWJ2SQ2phCXYj9PRI5zn7KvoV/kxYcQYFvl1JaNiP+dbHSbZ1wMdmIcrHykBf5w2Afx/HXNdyY60toPv06UNsbOwZt23bto1Gjc699tZf/2qszBsSGelcryYqKoqPPvoILy8vIiIizvu8k0WouzpX/Cf/kCguLq7ViaC4uNjVIfwj1Rq/Xo6WsBxD7HQMB+eBowQVOQB9yKfoza4Gi5dzv9JyoHquQrjT+59ZXMaa5Cx+O5bF2pQcCsodhFgc9LMe5nHLInqVLsezyECmf18OBt7Hals/CixhNPAwEulloouXmUDbX4ZtlZdQVAsu5lT1M6jKAgMnc+NJEyZMYOLEidUeU23j7vGD+7ehWuMvL0GLW4Qxdjpa/BIAVJOhOK56HNX4CjCd+KOwuKTaXtKd3v+kgmJWHcvmt2NZbErLoVRXNLaVMtS0l/7GOXRzbMBY4kua30D2hUxinq0P5WZ/wjxNNPYy0dfLhK/15DhmhV5aQm2Yh+Ri5sZaW0CvXr36gp8bGRlJfHw83bo5bwA8cuTIaSeBv0tISMDHx4fS0lLS0tIIDg4+4/CQvzo5BtpdnS9+g8GAwWDAZrNhs9lqMLLKUUpRXFyMzWZzy8ts1Rm/ytyH2vkl6o9vIT8ZAluj9Xkere0taD7n/0Pwgl7TDd5/pRQHsgtYdjSDZUcz2J7uvCG3vWcRd9u3MiB/Jq1L91Fqbk5S8OWssH5Bgq0bVoudxr5W+vlaaeRTe6eXu9DPoKysrNL7nsyNJ1WmB7q2fy/Oxd3jB/dvQ7XmxpRtqJ1foHZPh+IsCO2GdtlbaG1uQvMI4mIsb+IO77+uFDszcvn1aAbLEjLYn12ASYNLvHJ40mMdA/Jn0bg4iULfTiSEXc4cy/+RbOuAr81ME18rQ32dizyZaumQtZrIjbW2gP4nbrzxRj7++GNGjBhBTk4OM2fOZPHixed8jq+vLz4+PhQXF5ORkXHemwP/2iNbW39BzqUy8Z98XNO0Wt3G2h7f+Vxo/KokD7X3e/QdX0DiWrD5o7UdjaH9eGjYucbek9r2/pfrOltSc1iakM6yoxkk5BXhYYQ+nmm86rGcfoU/EViYT15Abw5FjuFrS39yrE1o4Ok8MfTws9LAw1yr2nQ+Vf0MqrLvydx4sWOqbdw9fnD/NlxwbizMdN7zsfMLSN0Ong3ROt6FIea2irmaa0Jte/+Lyx2sTc5i2YncmF5Uip8Z+nsc5UHbQi4tXoRXIWQFDGB/0OMstfan2BJCuJeF5n5WBvva8LcZa1Wbzudi5ka3LaCnTp3K008/TVZWFnPmzOHVV19l3rx5dOrUibFjx7J582ZatGgBwJNPPknr1nIHrXB/Sik4thF92yfOO8XLCp1LaF/7HVqLa9BMte9KQU0oKnewOuk4vxxJY9nRDHJKywmxKC6zH+Z5y2x6lC7HVOxDRsAgdjV8nb3WS9HNvkT5WOjoYaBFsBdeFrdNh0LUe0rpqPgVqO2fOldNVTpa9HC0vi+iNbsCzVA/f79zS8pYnpjBL0fS+S0pk6JynSi7g6stsVxeNotOZZugNJSUgMGssX3CIVsPzGY7Tfys9LIbaB7khc0Np5irCW77jRozZgxjxow54zaj0ciUKVNqOCIhLh5VnI36Yyr6tk8gfRf4RmHo+RRazG1ovucenlRX5ZSUsfxoBkuOpLMqKZNih05zj3Juse/kcsd3tCvdgcPUmKSQoSy23s0RWxc8LBaa+dkY5melkY8Vo+Yc82aXE4QQbkkVpKF2fom+/VPIOgSBrTH0fxmt3Vg0zxBXh+cSaYUl/JKQzi9H0tmQnEW5UnTwKuZ+j80MyptG06LDlJjbcTR0CD9ZnyfZ2g5/u4mmfjZG+NkI8zKj4cyN1lo6fK02cNsCuq7Iy8vj8ssvZ/fu3WzYsIF27doBzqn43nnnHex2O19//fVpN02+/fbb/PDDD/j7+zNt2jR8fX1P2b527Vqee+45HA5HxcIy119/PePGjWP37t0YDAYGDx7MCy+8UGNtFVVTMS/pto+dvc16uXMWjcveRGtyeb2cei6zuJSlR9JZFJ/G+hMnhk7exTzktZlBud/SpDCeIt9OHIkYyneWN0m3tCTQw0y0n42e/qcPzajNN8cKIc7M2du8HLXtE+dCJ5oBrfVIDFd9BRG93WqIQXVJyi9m8ZE0Fsen8XtaDkYNunvn84zXai7Pm0qDgkwK/HtxKOp2vrYMJMcSSainmWb+Ngb72Qiwn1oOSm48PymgXcxutzN//nyefPLJisfKysp46623WL16NZs3b+bFF1/kk08+qdienp7OvHnzWLNmDdOnT2fKlCmnLN2dmZnJk08+ybx58wgMDKSsrIzNm/9ckvzLL7+kdevWdO3albvvvvu8M5qImqVK853j97Z+AGk7wK8phkufR2s/Ds2roavDq3FphSUsPpLGkvh0NqZmAXCJdwHPeK9mUO43hOQfP3FiuIsvLJeTZwknzMtMa38bV/nZ6sSE/UIIUEVZqF1fo//+IRzfD0FtMFz2Jlq7MWj2gPMfoI5JyCtiUXwai+LT2JmRi8UAfbyzed1rKQPyv8M3v4ScwAHsD5rAXOtASsyBRPpY6epnpamfDS9ZGfUfkTPLOaiyQsjce+ZtSqF0HWUwQGX+2g1shWb2OO1hk8lEcHDwKY8dOHCAtm3bYrFY6N27N0888cQp2zdv3kz//v3RNI0rrriC22677ZTtCxYs4OabbyYwMBAAs9lMr169TtnHYDDQsmVLkpKSuOOOOygpKcFisfDjjz9e0A1D4p9TGbHov3+I2vU1lOajNR+ONvD1etnbnF5UwuL4dBbGp7IpJRujBj198nnBawWD8r4lIC+fvMD+7Gv6b+ZYLqfEHEgjbwuX+Nto5m/DU4ZkCFFnqJRt6L9/gPpjmvNKXKvrMQz7DBpdWu96mxPzilgYn8b8uFT+yMzDZtTo553JeM/F9C+YhVe+IitoMLtDXmentT+62ZvGPlb6Bdhp4muttTMKuSMpoM8lcy+OL7qcc5fKzgJtvH0rNOxcqX2zs7NPKWIdDsdZt/v6+nL8+PFTticnJxMdHQ3A4sWLefXVV/Hx8WHu3LkV+xQVFbFr1y6aNm3KnDlzsNvtTJ48mZkzZ3LXXXdVslXiH9MdqANz0be8h4pfBh7BaF0ewNDpbjTfqPM/vw45XlzKkiPpLIhLZUNKFhrQyyefl72XMSjvG3zziskNHMi+4P/wvfUyyk1+RPla6eNvo5mfTU4MQtQljjL0PXNQW96DpPXgHYGh90S0DnfUuytxyQXFLIhLY2F8KtvTc7EaNAb4ZHKX5wL6F3yPR56B48FD2NbwXXZZ+6GZPGnqZ2Wwv43GvlbMRsmNF4MU0OcS2MpZ+J6BUqpiHuVK/QUc2KrSL+vv709ubm7Fz0aj8bTtBw8eBJzFdEDAqZeuwsLCSEpKAuCKK67giiuuoGvXrhXbx48fj91u5+GHH8bDw4O77rqLhIQEsrOzuf766ysdp7hwqjgHfccXmLe8h54TB+E9MFwzzbkalunsc+zWNfll5Sw9ks7cw6msPXYcHUVPn0Je8F7J4Nwv8c8rJDdwIHuDX+Z360AcJl8a+1oZEGCnqZ/1tJWuhBDuTRWko2//BPPWD1D5x9CiBqJd/5Pzilw9mknjeHEpi+LTmHc4lU2p2VgMGv18snnbezED8qbjmadxPPgKtjScwh/WvmgmD5r5WRka4JzDvrbOz1yX1J9v4wXQzB5n7zVWCk3X0SpbQFdBdHQ0e/bsobS0lM2bN9O+fftTtnft2pU33niD5557jiVLltC7d+9Ttg8bNozhw4czevRogoKCKC8/dam0L7/8krZt26LrOj///DNhYWFMnTqVyZMnn9abLaqXOn4AffNk1K6voLwEvcUIzNfOwBDe3dWh1ZgSh87KxAzmHU5l2dEMShw6XX1KeMZnNUNzPiMgN4f8wL7sa/Yftlgur+hpHihFsxB1lkrbib75XecwDU1Db3UT5u6PYmjQ/vxPriMKyspZmpDB3MMprEk6jkLRy6eA132WMij3K7xyHeQED2JHyDtstw5AM3lK0exCUkDXAsOGDWP79u3s27ePe+65h3HjxvHII4/Qr18/bDYb33zzDQCvvvoqo0aNokmTJgwfPpzevXtXzMLxV4GBgbz++uuMHDkScI53fuSRR8742j169ODll19m2LBhhIaGyg2FF4FSCo6uQd/0P9T+ueARhNbtUbTO91Jq9MNit7s6xItOV4otqdnMPpTCwvg0ckvLae3l4GGfLVyZ8wnhuUcpDOjF/qZP85NlMEWmIKJ8rfQLsBHtZ5OplISog5RSqLilqI1vouKWOodp9JkEHe6kVPOoF7mxXNdZeyyLnw8lszQhnaJynS4+pUz0XcOw7E8JyM0mP2ggsc1eY6vlchwmb5r42RgSYKOJrw2zUYpmV9FUPZ+rJDc3F19fX3JycipWIoyLi6NJkybnXL66ykM4apnKxF/Z98JVlFLOOXzt9lr5GSi9HLX3J/SNb0LyZuf8pN0fR2t3C5rJVuvjP5/KxH8ou4DZh1KYcziFxPxiwmxwjUcsV+d9SvOyPyjxjSEu4Do2WK4kxxxOuLeFlgE2mvvb8KiBGwHrw2dwJn/Pexe6T3XGVFu4e/xQ+9ugHKWo3TPQN/7POa99w84Yuj+B1uoGNKO51sd/PueLXynF7sw8fj6Uwry4VDKKSon2VFxj+Z3huR8T4YinyL8HB/yvZZP1CgqNzg6FlgHOm6Rr4ipcXf8MzqYqeU96oIWoZqqsELX9c/RN/4OcI84xfCMXOFfDqgezaWSXlDHvcCo/HUxmR0Yu3maNod7HuNY+ja5FS1GGCBJDb2CWfTLJ5hYE2U3EBNppFWDH2yqzZwhRV6mSXNTvH6Fvfgfyk9Gir0QbPBktsp9bFmlVlVZYwpxDKfx4MJn92QUEWTWu8ozjWv0L2hZsoNzQkvhGt/Ct7SqOmxoR6mnmkkA7LQJqpkNBVI0U0EJUE1WUhdo6BX3zu1CchdZmFIbrf0Zr2MnVoV105brO6qTj/HgwmV8T0nEoRT+/At7zmc/A3G+w5NrIaHANv3j+wD5zN7ytJloG2Lk80E6Qh9nV4QshLiJVkOYc37x1CpQVosXc6rwaF9Ta1aFddCXlDpYdzeCHg8msTjqOUYPBflk85fUDvfN/wEAAKSEjmOvxPPEm54qArU50KMgc9rWbfDpC/EMq7xj6prdR2z5yzlHa/nYMPZ5A82vi6tAuurjcIubtTuLnQymkF5XSwkvxuN8mrsmeQmB2GgVBl7Gt+ftsNQ9EM3nQ3N/GDUF2Irwt9aLHSYj6TGXHo298E7Xjc9CMaJ3vxXDJo2je4a4O7aKLzSpg7o4E5h5OJae0nE4+Dib5/caVWR/inV1ETshQ1oV+xQ5LH8wmC60C7dwcaKeBp1lyo5uQAvosSkpKzrm9PoyBPt97UN+p7Hj09a+idn4JJhta1wcxdH0YzauBq0O7qArLHCyKT2Pm/iS2pOXgazZwjW8CI9TntM1fS7l3C+Kj7mae7RpyjQ2I9LFwWaCd5v42mY9UiHpAZe5HX/cy6o+pYPPD0GsCWpf76/xqgbklZcw5nMqsA8fYnZlHsNXAKJ8DXJ/7Ic1yd1Li15GDzf7NesuVFJn8aeJrZVignSZ+NplBww1JAf03RqMRTdM4duzYefc9WYC6q8rEr2naafNQ13cq67Dz5LDra7D6YegzCa3zfWg2X1eHdlHtzMhl5r4k5sWlkl/moLdfKW/5LmRIzqdYcuxkNBjBQs9/c9DUAV+biXZBdtoEesi4ZiHqCZW+B33dS6g934FnQ+cy2x3vQrN4ujq0i0YpxabUbL7bd4zFR9Io13X6+xXyoPeP9M+bhoFAUhreyE8eb3HU1IpAu4muQXZaBdplXLObkwL6b8xmM82aNTtt9b+/U0pRXFyMzWZz2x7oysRvNBoxm2WMKvytV8UeiKH/K2id70WzeLk6tIsmv6yceYdTmbEviT8y82hoNzDOdx/X506mUfY+Cv17srv522wwD0GZPGjhb+PGIDvhMkRDiHpDpe1CX/tfVOz34B2OYfB7aB1uRzPVvtmbqkt2SRk/HUxmxr4kDuUUEuVh4CH/7VyX9S7BWYnkBw5gc/NP2WIegNlspWWAnUuD7TTwkCEadYUU0GdgNpvPWzQqpVBKuXUB7c7x1ySVdRh9zX+chbNnAwyX/Q+t013OhXbqqN2ZeUzfl8TcQykUORz09y/mE9+f6JfzDZoeSEroKGZ5fEKysRkNPMz0C/GgRUDNTK8khKgdVEYs+qrnUXu/B98oDEM/Qou5rc6upqqUYmtaDjP2JbEgPg2lFIP8c5nkPZ0eeT+iq3ASI0czzXYDmcZwIrwtDA523vshQzTqHimghTgLlXsUfc1/UTu/AHsQhsvfcRbOdbRXpaTcwcL4NL7dm8j29Fwa2oyM99/HyJy3CTt+gMKg/mxp8QmbTAMxm620CrDT28dAhL+X/BEmRD2isg6hr/4Pavc05+Inwz51Fs7Gunm1sqCsnDmHUpi6N4m9WflEehp4xH8H12e9RcDxZPJDhrAmbCrbTH2wmc20CbRzubdGqJ/kxrpMCmgh/kblp6Cve8U5q4bFB8OAV9E6/6vO9jgn5hUxbV8S3+8/xvGSMi71L+dD/7kMyPoEo/InOXQ0szy+JNnYmFBPM5eHeNDC347JAEVFRa4OXwhRQ1ROAvraF1E7vgTPEAyDJqN1vLPO9jgfyi5g6t5EfjyYTGG5g4EBxTzlO5PeOdNBb0hi5BgW224kwxhGpI+FYcEeNPOzYdAkN9YHUkALcYIqzkbf8Dpq0ztgsmK49Dm0rg+hWb1dHVq1U0qx5thxvolNZPnRDLzMBm7wP8bovMk0ydpCSUAPdrWYwjrTIDSTndaBdi4L8SD4L3M21/NFTIWoN1RhBvral1C/f+DsVBj4+olOhbq31LZDV6xIzOCrPUdZl5xFgNXIWP/D3JzzP8Iy91IYPIDNLT9nk2kAZpOFtkEeXBXigf9f5myW3Fg/SAEt6j1VXoza+gH6upegvBjtksec8zjb/FwdWrUrKnfw88Fkvo5N5EB2Aa28DbwYvJ6rM97EnlnG8dCRzPd5mUOGNgTZTfQN8aBVoF3GNgtRD6nSAtTmd9DXvwaAofezaJc8UidvnM4tLef7A8f4NjaRhLwiOvnC/wJ/5YrMyViyzKSHjuZHr09INDSl4YkrcS0D7DK2uR6TAlrUW0p3oHZPQ//tWchLQut4F4Y+z6F5hbo6tGqXlF/Mt7GJzNyfRF5ZOZcHljDJdyqX5MxE1xsT3+wJVlqup8DgR3N/Gzc28CTcS+4WF6I+Uo4y1I4v0FdPgqJMtC73Y+g9Ec0jyNWhVbvDOYV8E3uUHw8kU6rrDAvM420+o0POQsp9WnGw+Qv8Zr6aMqMXLQPsjG7gSQPPujnWW1SNFNCiXtIP/4K+7AlI34XW6gYM/V5CC2zh6rCq3Y70XD7ffYRF8el4mAyMDEhibM6bRGTspCi4HxtbfcMmYz9sZjMxwR60D/HA2yJzkwpRHymlUAfmoi//Pzh+AK3dLRj6vlDnVlU9OXfzZ38ksOxoBgFWI+MDDzM6+3VC0g9Q0GAoqyK+Z5uhB14WI11CPGkXLPM2i1NJAS3qFZW+B335E6hDi6BRH4y3bUAL7+7qsKqVrhTLj2bw2R8JbErNJtLTxDPBvzMi41U80gvJDhvJPJ+3OWxoSYiHicENPGkhlyKFqNdU6nb0Xx9DHVmB1mQQhutmojXo6OqwqlW5rrMoPo1P/0jgj8w8mnsbeSV4HVdnvI4lUyMj/BZ+9PqSRK0x4V4WrmzgQbS/DYNciRNnIAW0qBdUYYZzvtJtHzvnKx3xI1rL6+rUEIXicgc/HUzm891HicstpLO/xgdBC7ks4z0MjiCONvkXP9hGka0F0szPyo0NvWSYhhD1nMpPRv/tGefMGoEtMYxcgNZsaJ3KC/ll5czcf4wvdx/lWEExvQMUXwT8QJ/jn6IcjYhv9m9W2K6nUPOhhb+N0Q29ZJiGOC8poEWdphylqC3voa95EVDOKem6Plinpl3KLSlj6r4kvtydQFZJGUOCSnnd7ys6Zf2Aw6cl+1q+xUrTlWC00zbITqcGnvjZ5FdfiPpMlRWhNv4Pff2rYLI5Vw/sdHedmss5s7iUr/cc5ZvYRArLHAwPLuB240e0Pr6YMr+O7Gj1MWuMgzGZLLQP9qBjA0+8ZAibqCQ5i4o6Sz/8C/ovD0HWQbRO92DoMwnNM9jVYVWb1MISvtydwPR9SZQ6dK4PzuaO/LdonL6ekqA+bGw9jQ2GPnhaTPRo4Em7YA9sJplNQ4j6rGKc89JHnDdPd3sIQ+9n6tSsQ4l5RXy2O4FZ+49h0OCmwDTG5b5KWNpOihoMZnWjH/mdS/C1mejbwJO2QXbMMtOQqCK3LaC/+OIL3n77bWJjY3nnnXd44IEHKrZNmjSJDz74gLCwMADatm3LtGnTXBWqqGEqO945lm//zxDZD+OI79FCYlwdVrVJyCvio53x/HQwGYvRwJjAY9yW/RLBqfvID7uGXxsvYrcWQ6DdxOCGnrQKsGOU8c1C1Hsqcz/60odRhxejNR2C4abFaIEtXR1WtTmYXcAHO+OZdzgVH4uBe4OPcMvxF/BPSyQnfCTzfd/lkNacEA8zw0I9aS7jm8U/4LYFdJcuXZg1axavvPLKGbffeuutvPnmmzUclXAlVVaE2vAG+vpXwBaA4ZoZaG1G1ZmxfIdOnBzmHk7Fz2Lk4QaHGJ3xAl6paeRE3MQcn0+J15oQ7m3hmoaeNPG11pm2CyEunCrNdy6EsvF/4B2O4YbZaM2vrjP5IfZ4Hu/viGdxfBoN7CYmNNzDjWkvYE8v4HjEWL73votjhNHY18oNDT2J8LbUmbYL13HbArpDhw4AGAxy2aW+U0pH7fkefeW/nZckuz/unLO0jkz2v/d4PlN2xLHwxMlhYoPdjEyfhDW9lMyI21jkdQfJNKSpn5VRoV6EeVlcHbIQohZwznU/Ff23Z6Ao05kXe/xfnVlBcEd6LlN2xPHr0QwaeZp4seE2RqS9gDnDSGrUnfzscRsZBNLC38YtYV6EeNSd8d3C9dy2gD6fGTNmsHTpUgIDA3n22WcZMGDAOffPyck5ZflNq9WK1Xr2G81O7uuuS3a6e/zgjF1LXI2+5hlI2QLNr8YwalHFJcna3rbzfQaxx/OYvD2eXxLSifA082LDnVyXNglLpkZK1D0s9biNbPxpEWBjYEPPimW2a6rddeU79Nf/u5sLjb8q+0tudD9KKbT4pehrn4P0nWitR6L1f6ViPufa3rbzfQbb03N4d3scq5KO09TLzBsNNzE89UWMZZ4kNn2MX22jyde8aR1o58qGnhXLbEturDx3b0NN5MZaW0D36dOH2NjYM27btm0bjRo1Outz7733XiZOnIjZbGbt2rVcd911bN68maioqLM+JzIy8pSfJ0yYwMSJE88bZ3Fx8Xn3qc3cNv7MWExrnsUctwi9QRccNyxGRfRxbisqcm1sVfT3z+BAdiEf7D7K0sTjNPIw8XLIZq5NfwlzqYXEqPv41eM28g2+tPA1MyTIiq/VCJRTVFReK+J3R+7ehqrGX1SF3xHJje5FS9uBcc0zmBOWo4f1wjFqJSq0m3Ojm+fGP47nM+WPo6xOzqaZl5G3glcxLP11DGU+HGnyFMvsoyk2eNHK30KHQCveFgOoMoqKympF/O7I3dtwMXNjrS2gV69efcHPbdiwYcW/e/fuTadOndiyZcs5C+iEhAR8fHwqfq5ML0txcTE2m80tx1K5a/yqOAe18mnU9k/BJ4ryoV9jbj8as8H9ph76+2dwILuAydvjWBSfRrinmVcbbuOa1EmYyu0kNXuMpfZbyMebtkF2ujX0wsfq2ja763for9y9DRcaf1lZ5QsKyY3uQRVmopY/idr1NQS0oGz4TCxtrsfshsMc//4Z7M7MY/J251CNpt5m3mm4jqEpL2EoD+RIi+dYahlFiWanfbAHXRt64uniqejc9Tv0V+7ehprIjbW2gP4nEhMTiYiIAODAgQNs376dmJhzz8Lg6+t7ykmisjRNc8sv10nuFL9+cAH6onugJBfDZW9Cp3+hl+kYDEa3acOZJOQV8e72eOYeTiHUw8x/w3YxIuU5TGVmkqKf4BfrGAo1T9oGedAt1PWF89+503fobNy9DVWNvyr7Sm6s/fS9P6EvuQ8cJRiGTIEOd6BKyzEYDG7ThjM5lFPI29sOs/hIOo29zPwvdDNXpbyAodSHIy2e5xfrTZRpNtoHe9I11BPPWrbUtjt9h87G3dtwMXOj2xbQU6dO5emnnyYrK4s5c+bw6quvMm/ePDp16sTEiRPZunUrJpMJo9HIlClTaNGihatDFhdIFR13Tr30x1S0pldgGPYJmk8j51ilMve6JPlXqYUlvLP1MD8dTiPQZmJS2F5uTH0Wc7oitemDLLKPIx8vYoI96NbQC+9aVjgLIVxLFaShL3kAtfd7tBbXYrjiAzSv0BPjOF0zpKs6JOUX8daWg8w9kk6Yh5nXw7ZzdcrzGEvtJEU/zWLrLZQYPOgQ7EmXWlg4i/rBbQvoMWPGMGbMmDNu+/rrr2s4GnGx6Pt+Rl/8L2fPylVfosXc5tZ/DQNkl5Tx0c4jfBN7FItB4/HQeMam/RtrWhHpje9licftZClf2gba6R7mhY/VbX9NhRAXgVIKtWcm+i8PAqrOTNmZUVTKBzvimL4vCW+zgWdCDzAqdQKWUkVK00dZbL+NAjydnQqhXrJqoHApOTOLWkkVHUdfcj9qz3doza/GcMWHaN5hrg7rHykud/DlnqN8tOsIDl3njobp3JHxb7xTksiKuoufvf5FmvKnlb+Nq8O9K+4cF0KIk1RBOvqie1D7f0ZrdSOGIe+jeYa4Oqx/pKCsnE/+SODzPxIwavBA6DFuS3saj9QsMhr/i8Ued/6lU8G71g1jE/WTnKFFraMfXoI+/3YoL8Rw9VS0tqPdumfFoStmH0rhrd8PkV5Uys2hhdyf8zxBSTvICR/F7MDHSNDDiPa1MiTcmyCZq1QIcQb6gXnoC+4EpWMY8QOGVte7OqR/pFzXmbX/GO9sjyO3pIzbQnO4J+sZ/I4dICt8DAv9H+aYCqGVv41rwr3xk04FUYvIt1HUGqqsEH35/6G2TkFrMsg5ZMM73NVh/SOrkzJ5dctBYo/nM7QBPG55g8bHllISeiW/NHubWL054R5mbmrkQ6gsgCKEOANVkoe+7DHU9s/Qoq/CMOwzNK8Grg7rgimlWHY0g9e2HORwTiHXhDp4rOC/hCWtpTB8BAuiP+egHkVjbwtjGvlUzHEvRG0iBbSoFdSxzTjmjoHcBAyD30Prch+a5n7TL520LyufVzYfYFXScToHWJgVMovOqZ9RHtCNdTFz2aw6E2Q1MTTYQotgb1lRUwhxRuroWhzzboWCVAzDPkXrcIdbX5HbmZHLy5sOsCk1m15BZt4K/pK2yTMoC76U39r/wna9DQ1sZq4KsRAd5O3WbRV1mxTQwqWU7kCtfQl9zQvQsBPG27ehBbVydVgX7HhxKW9vO8yMfUlEelmYEr6RwUmTwLMRu9t+xjKG4G02MSTci5YBNkqKi+UEIYQ4jXKUoa+ehFr/KoT3wHjzL2j+zVwd1gVLKyzhja2H+PFgMs19rHwetpK+x15BeUezo923/Kb642s2cWWEN9F+VrdfwEPUfVJAC5dRuYk45twCiWvQej+DofczaEb3vFRXput8G5vI5O1xKKV4OjyRMSmPYy7WSWo5iUWW0ZQbbFwa6kXHBp6YDJrbLpEqhLi4VHY8jjk3Q/IWDH1fROv5FJobLhYFUFLu4PM9R/lwRzwWo8YLEYcYmfwkxhIbCa3fYJHpBjCY6RfuTftgD4ySG4WbkAJauIS+fy76gvFgsmO8ZQVaZF9Xh3TBVhzN4KXNB4jPLeSmsDIezp5IYNIushvfxULvB8jQ/Wgf7EGPMC88ZL5SIcQ56LHfoy+8C2z+GMeuRgvv4eqQLohSisVH0nl18wGSC0oYG17IA8efwvfYYY43uY8FXveSrXvRMcST7mFe2EwyjE24FymgRY1S5cXOGwW3vIfW4hoMwz5H8wh0dVgX5EhuIf/ZuJ+ViZn0DLYwOfgrWiVNp7TBQJY0W85evRlNfawMi/AhwC6/akKIs1NlhehLH0Ft/xSt9UgMQz9Gs/m5OqwLcjC7gOc37GN9chb9G5j5zPYhzRJ/oiTsKhYGf8UBRwTR3laubuQj03UKtyXfXFFjVOY+HD+Pgsy9GAa/f+JGQfcb/1tc7uDDXUf4eNcRgqwmpkRsZXDisyivSHa0+5aVqj+BVjPXR/oQ6WN1dbhCiFpOpf2BY/YoyI5z6xsFC8rKeW97HF/sPkqYp5nPwtfQ/9h/0X1asK3996zSexJkMXFDpA+NJDcKNycFtKgR+p5Z6AvvAK9wjOM2ojXo4OqQquzk1EsvbNxPWmEJd4bn8q+0h7GlZJHc4hnmW29FN1gZEO5N+xAPDG54AhRC1Cx917foi+4B/2iM47egBbdxdUhVppRiYXwaL206QFZJGQ9GZHBn6iNY0opJavkSCyw3o2tmBkZ5ExMsuVHUDVJAi4tKOUqdQzY2v4vW5iZn74rFy9VhVdnRvCImbdjHisRM+oRY+Mr6Po2PzqEw4gZ+DJhIkh5M+yAPeoV7YzfLWD4hxLmp8hLnkI1tH6G1H4dhyBQ0s4erw6qywzkFPLd+H+uSsxjU0MzE0jeIOLqEgqgxzPV7imMOf9oHetArwhu7jHMWdYgU0OKiUXlJOH4eCcc2Yxg0Ga3rA253WbJM1/n8jwQmb4/D32rig8htDEqYiPJuwuaYH1mnuhPmYeaWKF9CZLJ/IUQlqOx4HD/fCGm73HbIRolD58Od8Xy0M56GHhY+i1hP/6T/oPu0ZmOHeWxwdCLMbuaWSF9CPCU3irpHCmhxUehxy9Dn3AxGK8axq9zyTvLt6TlMWLuX/dn5jI8o4aGMh/FMSia5xQTmWsahGW0MbuRNm0C72538hBCuoR9ahD53DFh8MN62Dq1hZ1eHVGUbU7KYuG4vCblF3NUon/tTH8CWkkNSy/8y3zIaTbMwuInkRlG3SQEtqpVSCrXhDfSV/0aLGojhmulonsGuDqtKckvL+d/WQ0zdm0hbfxs/NZhJu6OfUxI6jDkh3xHvCKN9sAe9I7xl6iUhRKUopaPW/Bd99SS06GEYhn+DZg9wdVhVkl1SxiubD/D9gWQ6B9p4P/hLWiRMpzjiOmYHPceR8mBiAj24VHKjqAeqpYDetGkTDz74ICaTieeff57BgwcDcN111/Hzzz9Xx0sIN6DKCtEX3InaMwOt1wQMfV9wu8n/lx5J57kN+8grLWdiVBpjk+7HUGTjj7ZfsIxBhFjN3BzlS0Mvi6tDFUK4CVWajz7vNtS+n5x5sfdENM19CkylFPPjUnlh435KHTovRiUyKvFBNKs/u9pPZ7nelxCLiZua+xIquVHUE9VSQD/++ON88cUXWK1W/vWvf5GYmMjtt99OdnZ2dRxeuAGVexTHD9c6p6i7bhaG1je6OqQqySwu5T8b9jM/LpUBDSz8p+x1wo4sIbfx7cz2foJ8zZf+Ed50kNk1hBBVoLLjcHx/DeTEYbhhDoYWV7s6pCpJKyzh2fV7WZqQwdBQM88Wv0jIkd/IaXovs70epQBP+jXypmMDyY2ifqmWAtpoNNK2bVsAFixYwOjRo8nJyZGxT/WEOroWx48jwGTDeOtatAYdXR1SpSmlWBCXxqQN+9CV4n+R+7g64VGUTzPWxcxls+pMEx8r10f54m11r950IYRr6fEr0H++EWx+GG/bgBbc1tUhVZpSip8OpvDipv2YDRrvR+7kiqNPofu2Zm2HxWxxtKWJtzM3+khuFPVQtRTQ5eXl5Ofn4+XlhcViYebMmYwdO5ZNmzZVx+FFLaZv/wx98X0Q3hPjiB/carxzWmEJz63fxy8J6QwNNfN84TMEHd1IWvRj/Gz/FwajjWFRPrTwt8kfg0KISlNKobZOQV/6CFrUAAzXzXSr8c7H8ouZuG4vvyVlcm24iYl5T+KfuIPU5hOYY7sDTbMwrKkPLQIkN4r6q1oK6MmTJ1NQUICXl3N+X6PRyLRp05g1a1Z1HF7UQkp3oC973Dm/c6d7MQx+F83oHmPflFLMi0vl+fX7MBk03ovczdCjT+DwbcuvMb+wW29Ju0A7fRr5yI0wQogqUXo5+pIHnfM7d3sYw2Vvohnc4359pRQ/HEzmxY378TQZ+SRyGwOP/huHf0d+jfmV3Y5mtA2w01dyoxAXXkCvW7eOWbNmcfDgQQCio6MZNWoUPXv2BEDTNEaNGlU9UYpaRZUWoM8ZjTo4H8OQKRi63OfqkCotq7iM59bvZUF8GleFmXg+///wT9xOSvRTzLbdjdVk5fomvrIEtxCiylRJHvrPI1FxSzEM+xRDxztdHVKlpReVMHHtXn49msH1ESYm5j6KT+JuUlo8y8+W27EZLYxo5kuUr+RGIeACC+j777+fDz/8kJYtW9KyZUsAFi9ezHvvvce//vUv3n///WoNUtQeKj8Fx6yrnDcL3jgPQ/QwV4dUaSuOZvD02lhKHTrvRu7hyqOP4fDrwNJ2S9mjR9M+yIM+jbyxGKVnRQhRNSo30Zkbc+Iw3LQIQ5NBrg6p0pYcSeOZdXtRwIdROxmU8H84/DvzS8yvxDqaSm4U4gyqXEC///77zJo1iwULFjB06NBTti1cuJDbbruNNm3acN997tMrKSpHpf2BY9aVoJdjvHWN29wsWFBWzsubDjBj/zH6NrDySul/aHB0LSnR/+Zn213YpNdZCPEPqNTtOGZeCQYjxrFr0ULauTqkSsktLefFjfv58WAylze08lLxBAITNpLc/BlmW+/AZrRwfbTkRiHOpMoF9Keffspbb711WvEMMGzYMP73v//x5ptvSgFdx+hxv6L/dD34NsY4cgGaT4SrQ6qU7ek5PPrbbtKKSnixcRI3JdyL8mrCiphF7NRbS8+KEOIf0Q8uRJ89CgJaYhw5D80r1NUhVcqmlCweW7WH3NIyXmtylBFH7kN5N2V5zBJ26S0lNwpxHlUuoPfv30///v3Pur1///7cc889/yQmUcvou75BX3AHWuPLnXeTW31cHdJ5OXTFR7uO8M62w7Tzt/C59X0ax88mq+kD/OD5CJrRgxHRfjKeTwhxwfTtn6Evugct+koM18xAs3i6OqTzKtd13tsRz5QdcXQJtDHD413C4+aS1fTBE7nRLrlRiEqocgHt5eVFRkYGjRo1OuP2zMxMPD1rfxIRlaNv/B/6sifQOtyBYehHbnE3eXJBMY+v2s3GlGzubZTPQ8l3YLJ6s7n9z6zTu9HC38ZlUb5yF7kQ4oIopVDrXkH/bSJa5/swDJ7sFquuJuYV8eiq3WxLz+GhyHz+lTQSg+UvudHPxmWNJTcKURlV/i0ZOHAgr7/++lm3v/baa1x22WX/KKjKmDBhAq1bt6ZDhw5ccsklLF++vGKbrus8+OCDNGvWjOjoaD744IOLHk9do5TCsfwpZ/HcawKGYZ+6RfH8y5F0rpyzkfjcAr4JW8LjR0dQGjqEb5v9wlbtEoY29ePKZv5yghBCXBCldPRfH0P/bSKGPpMwDHnfLYrn+YdTuXLuJpILipgRupgHjoygJHQwU5svZat2CVc09WNYMz/JjUJUUpUrohdffJHu3bvTq1cvHn300YpZOGJjY3n33XeJjY1lw4YN1R7o3/Xp04dnn30Wu93Ojh076N+/P8nJydhsNqZOncqePXvYv38/OTk5dO7cmYEDB9KqVauLHlddoPRy9EX3oHZ8geHytzFc8oirQzqvgrJyXt58kBn7khjUwMTLBQ/jl3aY/W0/YpF2FY08LVzfxE9WExRCXDDlKENfcDvqj2luM4Vnflk5L2zczw8HkhkWauK/eQ/inRbH/jYfOnOjh5XrmvjJaoJCVFGVCuglS5YwaNAg1qxZw0MPPcSoUaMqViFSStGvXz9Wr15dUVRfTH+9iTEmJgaHw0FGRgYRERHMnDmTe++9F6PRSEBAACNHjuS7775j0qRJZz1eTk4OSqmKn61WK1br2ceAndz3r89xJ2eLX5UVoc+5GQ4tRBv+DVq7MbW2jSfj2pSSxVNrY0krLOWFRvHcfPQ+9MBOLGi2jDg9nEvDvOjS0BNN02pVW+rqd8iduHsbLjT+quwvufHE46UF6D+PhPhf0a6dgdZ6ZK1t48m41icf5+m1ezleXMorUfHccOQ+9ICOLGz2K4f1CHqHedFVcmO1c/f4wf3bUBO5sUoF9LBhw0hOTqZt27bcdNNNzJgxg/j4eJRSREdHExgYWKVAq8uXX35Js2bNiIhwzgyRkJBAVFRUxfbGjRuzZcuWcx4jMjLylJ8nTJjAxIkTz/vaxcXFFxBx7XFK/CU5mOaOREvdSvnwmagmV0BRkeuCO4/icgfv7TrKN/uT6eBv4VPf92l6dA4ZTR7iR69HsBqsXBNpJ8TDWKs/p9ocW2W4e/zg/m2oavxFVfi9ltwIFGdhmnM9WsYflF/zIyrqslqfG9/ZmcDUAyl0DbDylfV9oo7MIaPJg/zk+Qhmg5WrIz1oILnxonL3+MH923Axc2OVCuiQkBCSkpIICQnh3nvv5ZprruGSSy6pUnCV1adPH2JjY8+4bdu2bRU3MS5btoz//Oc/LF269JR9TvaMQ+X+okhISMDH58/ZJSrTy1JcXIzNZjvltdzF3+NXRcfRf74asg9iuHkppoherg7xnLan5/Dk6lgS84t4qnE24xPvwWDxYGv7n1ird6Olv42BUT5Ya/EUTHXtO+SO3L0NFxp/WVlZpfet97mxMAP956sgJwHD6GWYwi7OOa+6/DU3TmiSzbijd6NZPNnW/gdW691pGWBjYKQP1lo81rmufYfckbu3oSZyY5UK6HHjxjF48GCGDx8OwNq1a7n00ksJDg6uymEqZfXq1efd57fffmP8+PHMmzfvlGEjkZGRxMfH061bNwCOHDlyWi/K3/n6+p5ykqgsTdPc8st1kqZpUJSJPmMQ5B7FOHo5WsNOrg7rrMp1ncnb4/hgZzzt/O38HPQtLeO/oShqNLP9nyNL92ZwEx/aBNrd5nOpC98hd44f3L8NVY2/KvvW69xYkIY+43IoSMM4ZiVaSIyrwzqrv+fG2UHf0iLuG4qjbma2/3Nk6j4MbuxDmyDJjTXF3eMH92/DxcyNVSqgX3nlFfr378+CBQtQSjF+/Hjy8vIICQkhJiaG9u3bV/zXsWPHqhy6ylatWsXYsWOZM2cOHTp0OGXbjTfeyMcff8yIESPIyclh5syZLF68+KLG465UQSr69MuhMB3jLStr9QpamcWlPLzyDzamZPNw42LuOTYWIw72dfiWJY4BBFvMjG7qR4C99s8WIoSo3VTeMWfxXJLjLJ6DWrs6pLM6XlzKQydy40ONi7n32FiMlLO/wzcs0QcSZDFzi+RGIapVlX+bhgwZwpAhQ5g1axY7duygrKyMXbt2sXPnTnbu3MmSJUs4cODARR83c8cdd1BSUsL48eMrHvv222+JiYlh7NixbN68mRYtWgDw5JNP0rp17U1+LpOf7Lw0WZKDccxvaEG1d5aSnRm53Ld8F8Xl5XwV+gs9496gNPxq5gS9RJLDn+5hXlwS6oXR4L5/KQshaom8RPSfroTyYmduDGju6ojOaldGLv+qyI1L6Rn3OqXhw5kX/BJHywO4JMyL7pIbhah2F/znaEpKSsW/IyIiTpkVoypjSC7UgQMHzrrNaDQyZcqUix6DO1M5CZi/Hwx6aa0/QXx/4BjPrt9HKx8z75smEZa2hYT2HzJXvwovo5FRzf0I9ZZVs4QQ/5zKjnPmRk3DOHYVml8TV4d0Vj8cOMYzf8uNR9t/yFx1FZ4GIyNb+REmuVGIi+KiXM8xm80X47CimqicI+jT+oNSGMb8hubf1NUhnVGpQ+eFjfuZvi+JUeGK59JvwWLxZHnbhexyRNMhxIMugSZ8PC2uDlUIUQeorMPO3KiZMdyyHM0v6rzPcYVSh86Lm/YzbW8SI8MVz5/Mje0Wsau8Ge2DPegaJLlRiItJBkTVMyovCce0gaAZKLthEaZa2ruSXljCv5bv4o/MXF6KOsKohHsobXAZ00PepUDzY0QLPyJ9LFWackYIIc5G5STgmD4QjDbKrl+IyffcN567SnpRCfct38WujFxeikpgVMLdztzY4F3y8eW6Fn5E+VglNwpxkUkBXY+o/FQc0y4DvQzDmFVgCXF1SGe0JzOPu5btwKHrTG/wMx2PfMDx5k/wne0h/G0WRkf742M1ue0E70KI2kXlHcMx/TJAwzD6VzAHuTqkMzqZG8sdOtMazKHTkfcrcqOf1ZkbfSU3ClEjpICuJ1RhBo4Zl0NpLsYxq8A3qlYuBLA4Po3HV++mmbeZj/RnaZi+jX0xX7FYXU5LfxuDGvthNsrNMEKI6qEK0pzFc3kxxrGrwKdRrcyNS46k8diq3TT1MvOR+XlC07eyr/1XLNYlNwrhClJA1wOqOBvHjMFQkHrihsHoWtdDoZRiyo543tp2mGGhRl47fgs2s5nl7RbxR3lT+kR4VyzHLYQQ1UEVHccxYxAUZ1fcMFgbc+MHO+P53+8ncmPWGGxGIyvaLWRXeTPJjUK4yAUX0Dt37mTjxo2kpaUBzlUKu3fvTvv27astOPHPqZI8HDOHQu4RjLesqJVzmRaXO/i/NbHMj0vlkag87j86GkdgV2aFfUSWCuDaFv409pU7yYUQ1UcV5zg7FvKO1dqZiIrLHTy1NpZ5h1N5OCqf+4+ORg/ozPfhH3NccqMQLlXlAjo1NZVRo0axatUqGjduXLEKYXp6OvHx8fTr14+ZM2cSElI7x9fWJ6qsEMesqyBjj3OFwZDa98dNRlEpd/26g31Z+bwXtZuhRx6lIGoM03xewGa2M7p5AH42uVAihKg+qrQAx6xhkH3Y2bEQ3MbVIZ0mo6iUu5ftYO/xfCZH7WHYkUcojLyFab4vYpXcKITLVfm375577kHTNA4dOkSTJqfO4BAXF8edd97JPffcw88//1xtQYqqU3o5+uybIGULxpuXooV2cXVIpzmcU8D4X7ZT7HDwXcj3tDvyCamtXmCWaTzhXlauivbHZjK4OkwhRB2iHGXoP90AaTudHQsNOpz/STXscE4Bty/dQWFZOTNCfiDmyEektZrETNMdhHtbuaqZ5EYhXK3KBfTSpUtZv379acUzQJMmTXjrrbfo3bt3tQQnLoxSCn3RPahDizDcOA8toperQzrN5tRs7lm2g2CrkanmlwhPX0ds+2/5Re9Pu2A7AyN9ZeUsIUS1UkpHX3AHKn4ZhlEL0cK6uTqk05zMjUFWI9/YXyYifS1723/LEn0AMcEeDIj0kdwoRC1Q5QLaz8+PuLi4s451jouLw8fH5x8HJi6c/tszqB1fYBj+DYZmV7g6nNMsiEvl8dV76ORv5sOif+FdlsXqdvP4vbwlfRt507mB3BAjhKh++oqnUX9MxXDtDAxNLnd1OKdZGJfKYydy4wfF9+NTlsHqmHn8Xia5UYjapsoF9AMPPMDYsWN55JFHGDBgQMVY57S0NFasWMG7777LhAkTqj1QUTn65smodS9juOxNDDFjXR3OKZRSfPpHAq9uOcg1YUZezhiJ2e7PvMbzSNRDuTraj2b+NleHKYSog/SN/0NteAPDoMkY2oxydTinUErx2R8JvLLlIFeHGXnlRG6cHz2Po44wyY1C1EJVLqD//e9/07BhQ95//31efvlldF0HwGAw0KFDB9555x3Gjx9f7YGK89N3f4e+9BG07k9g6P64q8M5ha4UL2zczzexidwXVcKjR2/EEdiN78I+ocDgy43NA2jgKUvACyGqn/7HVPRlT6D1+jeGbg+6OpxT6Erx4sb9fB2byL+iSnj06Ej0wC7O3Kj5cWMryY1C1EYXdAvv+PHjGT9+PKWlpWRmZqKUIigoCIvFUt3xiUrS435Fn3crWrsxGAa+5upwTlHq0Hli9R4Wxqfy3ybJ3BQ3juKIEUwLeAOz2c5NLQLwtcrd5EKI6qcfWow+fzxah9sx9HvJ1eGcotSh8+SaPcw/nMqLTVK4Oe42iiOuY1rAm5Ibhajl/tFvpsViITQ09JTHCgoK2Lp1K3379v1HgYnKU6k70H+8Dq3xZRiu/BxNqz13ZxeWObhvxU42JGcxOWIHV8Q9SU70w0y1P0Gw3cI1zQOwm2tPvEKIukOl/I7+0w1oTa/AMPTjWjV+uLDMwf0rdrE++TiTI3cyNO4JcqMf5Fv7/0luFMINVPuftgcPHmTAgAE4HI7qPrQ4A5V3zDnXs39zDCO+RzPWnkt92SVl3LF0O/uzCvg8dBk9j75KSptXmWUYQxNfK8Oa+svSs0KIi0LlJuKYNRyC2mC4biaaofb05GaXlHHnUuf895+FLadXwiuktnmFmYaxkhuFcBNVzignxzyfjRTONUeV5juLZ8A4cj6axcvFEf0ppaCYcb9sJ6OohKlBs4g59hmHYz5inrqKmGAPBkb5YKhFvUFCiLpDlebj+H44GEwYb5yLZvZwdUgVTubG9KISvg3+gfZJnxDX7kPmMlxyoxBupMoFtNl87h5OpVStukxWVyndgT5nNGQdwDh2DZp3mKtDqhCXU8itv2xD13Vm+HxEs7Q57Gn/DUsd/egZ7kX3UC/5jgghLgqlO9Bn3wxZhzDeuhbNq6GrQ6oQl1PIbb9so1zX+c7nY5qlzmZPzDcs1fvRM8yL7mGSG4VwF1UuoL29vfn3v/9Njx49zrj9wIED3HPPPf84MHFu+rInUAcXOBdKqUUraR3IzueWxdvwNWt8aXmZsOz1bI2ZxRpHF/pH+tCpgaerQxRC1GH6ssedi0iNnI8WEuPqcCocyM5nzOJteJs1plteJSxrDb+3n8nq8q6SG4VwQ1UuoLt27Yqu6/Tr1++M2/38/FBK/ePAxNnpW6agNr+DYfD7GKKHuTqcCrHH8xi7eBvBNgPfqKcIKDjA2nY/sbWsNUOa+NImqPZcRhVC1D3O3PguhiEf1KpFpGKP53Hrkm0EWQ18w9MEFOxnXcxPbClrI7lRCDdV5QL6/vvvp7Cw8KzbIyMj+fLLL/9RUOLs9IML0Zc+hNbtEQxd73d1OBV2ZeRy65JtRHgY+ar8Ifz0DFa0mcue8iiuivYnWhYBEEJcRBW58ZJHMXT5l6vDqfBHRi63/rKNcA8jX5U9jJ+ezoo2c9hT3lhyoxBurMoF9HXXXXfO7f7+/tx2220XHJA4O5W+B332KLToqzBc9qarw6mwLS2HcUu308zLyBcl/8KbAha1nE2cI4xrmvsT5Wt1dYhCiDrMmRtvcubGgW+4OpwK29NzuO2X7TT1MvJl6X14qzwWt/yZw45wyY1CuLnaM6+POCdVnI3jx2vBtzGGa6ahGYyuDgmATSlZ3LF0B619zXxafBdeWhnzo38kSQ/l+hYBhHnL4jpCiIunIjf6RNaq3Lg5NZvbf9lOK18znxXfjZdWwvzmP5Kkh0luFKIOuOACetWqVWd8XNM0rFYrTZs2JSgo6IIDE39SSkefOwYK0zGO31Jrpqtbn3ycO3/dQQd/M58UjMdu1JjX9AeOqQZc3zJQlp8VQlxUtTU3bkjO4o5ft9PB38LHBbfjYVTMa/qj5EYh6pALLqD79+9fMd3OyZsG//qzpmkMGTKE6dOn4+fn988jrcf0Vc+jDi7EMGohmn8zV4cDOHue7/x1B10CzHyYdxs2s4U5TWeRqoK5oWUgIXKCEEJcZLUxN25OzeaOX7fTOcDMR3m3YTOZmNPsB8mNQtQxF7xO6Pz58+nWrRuLFy8mMzOTzMxMFi9eTPfu3fnxxx/57bffSExM5NFHH63OeCtMmDCB1q1b06FDBy655BKWL19esW3SpEmEhITQsWNHOnbsyC233HJRYqgJ+t6fUGv/i6H/y7XmrvJtaTncsXQHHfxMfJg7FpvFzuymP5BKMDe0khOEEOLiq425cXt6Dncs3U4HfzMf5d6KzWxldjPJjULURRfcA/3YY4/x7bff0q1bt4rHBg0ahI+PD7feeiv79u3jvffeY/To0dUS6N/16dOHZ599Frvdzo4dO+jfvz/JycnYbM47mm+99VbefLP23Gh3IVT6HvT5t6G1uhGt51OuDgdwzrYxbul2Wvua+CT/VqxWb35q+h2ZKpAbWgYQ7CEnCCHExVVbc+Ntv2ynpY+Jj/PHYbV68nPTmWRIbhSiTrrgHuiEhIQzH9Bg4OjRowBERUWRk5NzoS9xTkOHDsVutwMQExODw+EgIyPjoryWK/x502ATDFd9UStWpzo5l2kzLyOfFt3l7F1pIsWzEKLm1NbceNsv22jmZeKzoruwm8zMbvKdFM9C1GEX3AM9ZMgQ7rjjDt599106d+6Mpmls3bqVRx55hCuucF5O27ZtG9HR0dUW7Nl8+eWXNGvWjIiIiIrHZsyYwdKlSwkMDOTZZ59lwIAB5zxGTk7OKQvAWK1WrNazTzF0ct+LsWiMUjr6nDFQmIFh3CYwe1b761Q1/v1Z+YxdvI1Gnka+KL0fL0MZs5v9SAaBXN/Sn0C7qcYX0LmYn0FNkPhdz93bcKHxV2V/yY3ndiC7gLGLnfM8f1F6319yY5Dkxgsk8bueu7ehJnKjpi7w3cnOzubhhx9mxowZOBwOAIxGI6NHj+add97Bz8+PTZs24XA46NmzZ5WP36dPH2JjY8+4bdu2bTRq1AiAZcuWMX78eJYuXUrLli0BSElJITAwELPZzNq1a7nuuuvYvHkzUVFRpx0rNzcXX1/f0x6fMGECEydOrHLc1cGw5S1Ma56l7JqfUE2GuCSGv4rPK+K25bsJsmh8qz+CnyODBc1/IkkLZ3hjTwJttWPaKCFE5eTm5hIaGkpOTg4+Pj5n3Udy47kdOZEbAywa3+qP4u9IO5EbIyQ3CuGGKpMbT7rgAvqkvLw8Dh8+jFKKZs2a4e3t/U8OVyW//fYbY8eOZd68eXTo0OGs+w0ZMoS7776b66+//rRtJ08SCQkJp7xZlellKS4uxmazVeslRJW4Dn1qP7QeT2Do/0q1Hfe016lk/MkFxYxcuBW7EaZrTxFQGs8vLWdzWG/EiJYBNHThTTEX6zOoKRK/67l7Gy40/tzcXPz8/CpVQEtuPLOUghJuXLjlRG78NwGlh/ml5c8c1iMlN/5DEr/ruXsbLmZuPOkfL6SSmJhIXFwcAHa7vaIX+GJbtWoVY8eOZc6cOacVz4mJiRXDOQ4cOMD27duJiYk55/F8fX3P+2adiaZp1fblUoWZ6HNuhvAeGPr9t0a+tOeKP6u4jPFLd4BSfGX4DwFFB1nRZjYHHRFc1yKAUK/asRBAdX4GriDxu567t6Gq8VdlX8mNp8suKWPc0u0opfjS+CIBhftY2WY2Bx2NJDdWI4nf9dy9DRczN15wAZ2WlsbYsWNZunRpxTzPOTk5DB48mG+++Ybg4OALPXSl3HHHHZSUlDB+/PiKx7799ltiYmKYOHEiW7duxWQyYTQamTJlCi1atLio8fxTSin0BeOhrADjtTPQDK5dJLKwzMGdv24ns6iEmV6TCc3bxuq2P7OnvDHXNA8gQlbREkLUgNqYG+9Y6syNM7zfJyxnC2va/czu8iaSG4WoRy44E913333k5uYSGxtb0eu8d+9exo8fz/3338+sWbOqLcgzOXDgwFm3ff311xf1tS8Gtelt1IF5GEbOR/Np5NJYSh06963Yyf6sAqYFfUeT9MVsivmJbWXNGR7tT5Tv2S/fCiFEdaqtuXFq0EyapS9kc8yP/F7WQnKjEPXMBRfQS5YsYfXq1acM2WjVqhVTpkyhf//+1RFbvaGSNqKveAqt++MYoq90aSy6Ujyxeg8bkrP4PGwF7ZK+YFfMVNaXt2doUz+a+dtcGp8Qov74Mzc+USty45NrnLnxs7CVxCR9zh8x37KuvIPkRiHqoQsuoM1mMwUFBac9XlhYiMnk2kts7kQVZeGYPQoadrmoN8ZUKhal+M+G/SyMT+W9yN30PPIy8e0+YLneh4FRPrQKtLs0PiFE/VGRG0O7Yuj/smtjUYoXNu5n/uFUJjeOpVf8S8S3m8Iyva/kRiHqqQteSGXEiBHcfvvtLFmyhOzsbLKzs1m8eDF33nnnGWe7EKdTSqEvvBNKcjBe+x2a0bWT7X+46wjf7k3kxcbJDDnyKKmtX2AOV9M9zIsOIZ4ujU0IUX/8mRtza01u/CY2kRebpDI0/mHSWv+HOVwjuVGIeuyCu4onT57MI488wvDhw0+ZB/r222/nf//7X7UFWJepXV+j9v2EYcQPaH6NXRrL3MMpvLn1EA83LmRU3Dhyox9kpmE87YLs9AzzcmlsQoj65c/c+COa7+nz99ekeSdy44NRhdwUdyu5ze7nO8PttAuW3ChEfXbBBbSHhweffPIJb7311inzQHt5SUKpDJUdj/7LQ2gxt2Fo5doe+00pWfzf6j1cFwYPHLmRosib+Nb+fzTxtXFZY1+3nsJGCOFeKnJj+3EYWo1waSybU7N5cvUerg3TeCjhRoojR/Gtx9PO3BgluVGI+qxKBfRzzz1X6X1feOGFKgdTXyjdgWPebWDzxzDoXZfGEpdTyL3Ld9I5wMxL6ddT3qAvU/1fIdhuYVhTfwxyghBC1JCK3GgPqBW58Z5lO+gUYObl9OspD+nDVL9XCfaQ3CiEqGIBvXr16krtJ3+Vn5va9DYcXY3xlhVottOXyq0px4vLuH35boKsRqYU3YfJowHTG3yAzWzj2uYBmI3yOQohas4pudFa9cVbqsvJ3BhoNfJB8f2YPEKY3vBDrJIbhRAnVKmAXrFixcWKo95QabvQf5uI1v0xtKh+LoujuNzBg2v2UlBWzjeeb+BTmsG86AUUG7y5qUUANtMF318qhBBV9mdufNylubGk3MFDa/ZRUFbOD55v4FOQxvzoBRRrkhuFEH+S+eZqkCovwTF3DAS0wNDvvy6LwzmfaSz7sguZFvIzEWkrWBMzh6N6KCObB+BjNbosNiFE/fNnbmzp8tz4xJpY9mYXMC3kZxqlLWdtuzkk6GGMbCW5UQjxJymga5C++nnIiMU4fjOayXWT7r+77TCL4tOYHLGdDokfsLf9V2wtb8PwaD8aeLp2uighRP1zam503Wp+726P+0tunMK+9l+zxdFWcqMQ4jRyLaqGqITVqPWvY+j3IlqDDi6LY2F8Ku/tiOexqGyGJv4fqa1fYIl+GX0ivImWlbSEEDWstuTGRfFpvLc97s/c2OoFFktuFEKchRTQNUCVFeJYMB4ieqF1f8JlccQez+PJ1Xu4MlTj3iOjyY28jVnGcbQLstOloSwGIISoWX/mxt4uz41PrN79Z25sdCuzTOMlNwohzkqGcNQAffV/IDcR48gFaAbXjKE7XlzK3ct20tTLzKvHR1Me0pvv/P5DuKeVgTKfqRDCBf7MjQtdmhvv+UtudIT04jv/FyQ3CiHOSQroi0yl/I7a+D8MfV9AC2zpkhjKdJ0HVuyiqKycGaYXsZlMzAqdgtVo4cpmfhgNcoIQQtSsitzY779ogS1cEsPJ3FhYVs70k7mx4RSsJsmNQohzkwL6IlJ6OY4Fd0JwW7QeT7osjv9uPMCW1By+DV1IWOoGVrRbQBb+XBfpIVMyCSFq3J+5sR1a98ddFsdLm07NjSvbzSeLAMmNQojzkgL6IlIb34K0HRjHbUQzuuYO7u/2JfHt3kRebJxEt/i32R/zGbscLbi2uR++Ft0lMQkh6rfakBtn7k/im9hEXoxKpNuRt9nf7lN2OlpKbhRCVIr8iX2RqOMH0Fc/j3bJo2ihXV0Sw7a0HJ7fsI/RjRQ3H7mDrOhHWaSuoHeEN1G+rpsqSghRf6njB0/kxsdclhu3p+fw3Pp93NJIcXPCnWRFP8IihkpuFEJUmvRAXwRKKfSFd4NXKIY+/3FJDFnFZTywchft/Cw8k3oDpQ0u4zv7ozT3s9FV7ioXQriAUgp90d3gFYahr+ty4/0rnLlxYuoNlDUYyHf2xyQ3CiGqRAroi0Dt+ByVsBLDzUvRLDWfkHWleGzVborLHbxreh6z1YfvGkzG22JhcBPnXeVKqRqPSwhRv6kdX6COrHDmRrNHjb/+ydxY8rfc6CO5UQhRRVJAVzOVn4y+7Am09uMwNLncJTF8sCOeVUmZfB76K2HpW1nRbjE5yofR0f5YjDJqRwhR81R+MvryJ9Daj3d5bvzilNzoK7lRCFFlkjGqmb70ETBaMQx80yWvv/bYcd7ZfpgHG2XTN/l1DrZ6m53lTRnWzA8/m/y9JIRwDX3pI2CwYLjM9bmxT/LrHGz1luRGIcQFk6xRjfT45ajYWRiGf4PmEVjjr59SUMwjv/1BryAz9yeNIy9qHAu4kt7h3jT2laVohRCuUZEbr/4WzR5Q469/Mjf2PCU3XiW5UQhxwaSAribKUYb+y4POJWnbjanx1y/TdR767Q/MBvhf0ePg04yZPs/QxMdKt1C5MUYI4Rp/5sZL0dreUuOv/9fc+NZfcmNTX8mNQogLJwV0NVFbp0BGLMbbt7pk6df/bT3E9rRcpjWcS0DGARY1W4LB6MGQpn6yFK0QwmXU1imQubd25Mb0/Sxq+4szNzaR3CiEuHAyBroaqPwU57ymne9Fa9ipxl9/ZWIGn/yRwP9FptIl+T32tXqXQ3okw5r5Y5fVtIQQLlKRGzvdi9agY42//p+5MY0uye+x/y+5UVYaFEL8E26bQSZOnEhMTAwdO3akY8eOzJw5s2Kbrus8+OCDNGvWjOjoaD744IOLGou+4mkwmDH0++9FfZ0zSS8q4cnVe+gfYmJ8wu3kNr2HJWowvcO9CfOy1Hg8Qghx0p+58cUaf+2TubFfiInxCePJa3I3ixkiuVEIUS3cdgjHk08+yUsvvQTAsWPHaNWqFYMHD8bf35+pU6eyZ88e9u/fT05ODp07d2bgwIG0atWq2uNQietQu77GMPTjGr85RinFU2tiAcVrBY+i+7djptdTNPG20kUWBBBCuNCfufETF+fGx9D92vKd99OSG4UQ1cZte6D9/Pwq/p2Xl4emaei6DsDMmTO59957MRqNBAQEMHLkSL777rvqD0J3OG+OadgFrcMd1X/88/gmNpGViZm87reEwNLDLIn8EIPRJuOehRCudUpuvL3GX/5kbnzN/xcCSw6yJOojyY1CiGrltj3QAJMnT2bKlCkkJibyxRdfEBjonDouISGBqKioiv0aN27Mli1bznmsnJycU1agslqtWK3Ws+6vlMLwx1eQug3DretAM9ToClb7svJ5ZfNBbgsvol/Sa+xv8yEHy8O5sZUfNuP5V9M6ud2dV91y9zZI/K7n7m240Pirsv+F5cYva0Vu7J/0KvvbfCC50c1I/K7n7m2oidxYawvoPn36EBsbe8Zt27Zto1GjRjz00EM89NBD7NixgzFjxnD55ZdXFNF/7WWozBsSGRl5ys8TJkxg4sSJZ39CUSbmdZNwtBlLaUAHKCqqRKuqR3G5g4dW7CLK08iTqXeSF3odi7Sr6N7Ahr/RQVEVYikuLr6IkdYMd2+DxO967t6GqsZflRxxYbnxPy7LjQ+v/HtuHC650U1J/K7n7m24mLmx1hbQq1evrvS+HTp0IDw8nJUrV3L99dcTGRlJfHw83bp1A+DIkSOnnQT+LiEhAR8fn4qfz9fL4lj7Nkp3YLrsNSx2e6VjrQ6vbdhPQn4xs/0/wuqw8n3QSzTystIjwrfSlyeVUhQXF2Oz2dz2kqa7t0Hidz13b8OFxl9WVlbpfaucG9e86bLc+PrG/RzJK2a2/8eSG924DRK/67l7G2oiN9baAvp8YmNjad26NQCHDh1i27ZttGnTBoAbb7yRjz/+mBEjRpCTk8PMmTNZvHjxOY/n6+t7yknifAyXPktJ1BBsXg1r9Mu1/GgG3+5N5PlGCbQ4+hMb2/9MgebLiKZ+GAxVH9KuaZpb/nL8lbu3QeJ3PXdvQ1Xjr8q+Vc6NvSdSEjnYJbnxm9hEno88SouEH9kY8xMFBsmN7twGid/13L0NFzM3um0B/fTTT3Pw4EHMZjMmk4n333+/oqAeO3YsmzdvpkWLFoBzxo6T26qLZvNDRfSp1mOeT2ZxKU+t2cOAEBNjEv9FZvQjbNC7MqypL94WY43GIoQQZ6J5BKMiB9Toa/6ZG82MOXovmc0eZoPqxrAoyY1CiIvDbQvoOXPmnHWb0WhkypQpNRhNzXh+/T50pXil6P/Q/drxg/0hWvnbaBlYs5dJhRCiNpn099zo8bDkRiHEReW2BXR9szAulYXxabwbupGgtN38EvMrZoONAVG+rg5NCCFcZmFcKgvi03gnbBNBqZIbhRA1w23nga5PMopKeW79Pq4IUQxLfpaEli8SWx7FFU18ZTlaIUS9dTI3DglRXHnsGcmNQogaIxnGDTy/YR+g+E/+Y5SF9GeOYRRdG3oS4XP2O+GFEKKuO5kbX5DcKISoYTKEo5ZbEJfKovg03g3dQEB6PD83+JJAm5me4d6uDk0IIVzm77lxdsMvCbRKbhRC1Azpga7F0otKeG79PoY2gCuTnyMu+j8kEcYVTf0wGdx3WhkhhPgnzpQbE5XkRiFEzZECupZSSvHc+n1oKCblP0ZpcF/mG0fSPcyLIA+zq8MTQgiXOCU35j0uuVEI4RIyhKOWmh+XypIj6bwXvomA1MPMa/oZgTYL3Rp6uTo0IYRwmZO5cXL4ZgJSDzEvejmBVsmNQoiaJT3QtdDx4lL+s2E/wxpqDE16hoTo54hXEQxu4odRLk8KIeqp48WlTDqRG4clTXTmRl1yoxCi5kkBXQu9uuUg5UrnubwnKQvuzTzTaLqGetLAUy5PCiHqr1e3HMRxMjcG9XLmxoaSG4UQNU+GcNQyG1Oy+OFAMv8N3UlQ2l4WRK/A12ymR5jcWS6EqL/OlBt9zGZ6yKwbQggXkB7oWqTEofPMur109jcyMvkpkpo/w0FHIwY1ljvLhRD116m58f8qcuNgyY1CCBeRHuha5JNdRziSW8Rcn7dRfm2ZbRpDpwaehHlbXB2aEEK4zJ+58S2UX7sTudFDcqMQwmWkgK4l4nIKmbIznjsbZtAieRHL2y3E02qld7jcWS6EqL/+zI3ptEhezPJ2C/GwWuktQzeEEC4kQzhqAaUUz67fSwObkfvTHyAr6nZ2EcNlUb6YjfIRCSHqp1Nz44MVufFyyY1CCBeTHuhaYM7hVNYlZ/F5g8XYHEZmej9OywAbUb5WV4cmhBAuM+dwiuRGIUStJAW0i2WXlPHSpv0Ma6DTL/Ut9rT+mCKDL30b+bg6NCGEcBlnbjzwl9z4keRGIUStIdfAXOyNLQedd5jn/5vSkP4sNQyjV4Q3Xhajq0MTQgiXef2vuTG4H0sNV0puFELUGtID7UK7MnL5bv8xng07SEjKbuY3XU6I3UyHEA9XhyaEEC6zKyOXmfuP8Wz4QUKSnbkxWHKjEKIWkR5oF1FK8cLG/UR7mxid/CipTR/hkGrMwChfDJrMayqEqJ8qcqOPidHH/syNl0luFELUItID7SLz4lLZmpbDNwEzMHhGMMd+NzFBHoR6ybymQoj668/c+B0GzwjmSm4UQtRCUkC7QGGZg9c2H2RQUBm9Mr5ka9vpYPTg0giZ11QIUX+dzI2Dg8volf4FW9tORxk96C25UQhRy8gQDhf4+I8jZBaX8u+iFyhucDlr6EOfRt7YTPJxCCHqr5O58enCv+TGCG/skhuFELWM9EDXsKT8Ij7ZdYQ7QpKJTN3C3MbLCfO20CbQ7urQhBDCZc6aG4MkNwohah/5s76GvbL5IL5mA/dkPE5G5J3Eac0YEOmLJjfHCCHqMcmNQgh3IgV0DdqYksXC+DSe9FuDp6Gc+d4P0DbIToin2dWhCSGEy/yZG9dKbhRCuAW3LaAnTpxITEwMHTt2pGPHjsycObNi26RJkwgJCanYdsstt7gwUieH7pyaqaO/kWtSJnG48dMUGgLoFS43xwgh6q8/c6OJa1Kel9wohHALbjsG+sknn+Sll14C4NixY7Rq1YrBgwfj7+8PwK233sqbb77pyhBPMXN/ErHH8/kh4AuUTwsWW0ZySainrKolhKjXZh04diI3flmRG7uHekluFELUam5bQPv5+VX8Oy8vD03T0HXddQGdQ0FZOW9vO8y1DUrpmPodG9vOxG6x0aWhl6tDE0IIlykoK+et3w+dyI0z2Nj2O+wWG50bero6NCGEOCe3LaABJk+ezJQpU0hMTOSLL74gMDCwYtuMGTNYunQpgYGBPPvsswwYMOCcx8rJyUEpVfGz1WrFarWedf+T+/71OWfz+R8J5JaW82jecxQ1vIINWm+GRnhh1Cr3/IuhKvHXVu7eBonf9dy9DRcaf1X2v6i5cfdfcmODIWzQLpXcWA3cvQ0Sv+u5extqIjdqqpa+O3369CE2NvaM27Zt20ajRo0qft6xYwdjxoxh5cqVBAYGkpKSQmBgIGazmbVr13LdddexefNmoqKiTjtWbm4uvr6+pz0+YcIEJk6c+I/bcby4jCsW/M5Iv0QmZt7J3DbLKfJozjVNPOXuciGES+Tm5hIaGkpOTg4+Pj5n3UdyoxCiPqlMbjyp1vZAr169utL7dujQgfDwcFauXMn1119Pw4YNK7b17t2bTp06sWXLljMW0CclJCSc8mZVppeluLgYm812zmT/+c6jGDS4P3sCx6PuJF5ryqgoXzw8XLssbWXjr83cvQ0Sv+u5exsuNP6ysrJK73uxcuMXZ8iNN0lurBbu3gaJ3/XcvQ01kRtrbQF9PrGxsbRu3RqAQ4cOsW3bNtq0aQNAYmIiERERABw4cIDt27cTExNzzuP5+vqe96+NM9E07awfztG8IqbvS+LBoH34Hc9lqtf9tAqwEeZ99pNPTTtX/O7C3dsg8bueu7ehqvFXZd+LlRun7UviwaD9p+TGUMmN1crd2yDxu567t+Fi5ka3LaCffvppDh48iNlsxmQy8f7771cU1BMnTmTr/7d370FRnWkawJ8GBFHuLaByacJVBIMaURAi5GKpoJkMFVktY8RMrFIjiZsaM7OajDIzMbMxtWvG0WTKZIBoJkE3GFYziSFZbyHRqKPcY5CLYMDhpoAREeh3/3DooeWiDUj3aZ5fFVWe853T/XxHffvlq9PdZ8/CysoKlpaW2LlzJwIDA4c943/9vQxO1hZIqt+Iyz7r0GzhjARPw1+IiIjMyX+f66qN/4Eq1kYiUiDFNtBZWVl9jqWnpw9jkt4VN7bgf8uu4LduZ2DbZI3PbJMw3X0s7G340UxENHIVN7Ygq/QKUv5ZGz+3TcJ0dzvWRiJSFMV+kYqpe+NMKTRjLfFU7WZc9H4JWit7fmwdEY14XbVxce1mXPS6XRtn8GPriEhhFLsCbcpO1lzFsR8bsMP1C1jKBHxh/W+ImGCH0Vb8fYWIRq6u2vhH1+zbtdHmdm20YW0kIoVh1RpiIoL/PHMRDzpaYH7dmyjyehk21raY6sYVFiIaubrXxgV121DI2khECsYV6CH2+aU65NY3Y4/6Q3Q6BOP/rBbhkYl2GGWp3HexEhEN1uE7auMR1kYiUjA20ENIK4K3zpUhWi2IbEjF6eB02I+2Rui4McaORkRkNFoRbNerje/DgbWRiBSMt3AMoc8ravHDtZ/wQsc7uOUSjm8sHsXsiXawtOAKCxGNXPq1cSa+sXgEkR72rI1EpFhcgR4iWhH88Xw5op3bMf3qARyffADqMaMQpLY1djQiIqPpWRs/uV0bXUYbOxoR0YBxBXqIHL5Uhx+u/YTk9rdx0+0xnLMIR5SHPSwU/A0+RESD1b02tro9jnMWM1gbiUjxuAI9BLQi2HG+HLOd2/HQ1UP40vsQJowdBV8n0/laWiKi4XZnbcxmbSQiM8EV6CHwZWU9vr96Hckdu3HDNRaFFlMR5Wmv6O+PJyIarOzKum618REUsTYSkZngCvQgyT9XWCKcOxF+9RNkexyAl4M1vBy4wkJEIxdrIxGZM65AD9JXVfUobGxBcud7aFVHo8gqHBET+ZXdRDSyfVlVj6LG6/+sjQ+zNhKRWWEDPQgigj/mVmCWkxazmv8HJ11fhIe9NTztucJCRCNX1+pzV238lrWRiMwMG+hBOFZzDYUNLVgnaWhziUCeZQRmTeAKCxGNbMeqr6Kw8bquNuZbzmJtJCKzwgZ6gEQEbxdUIdxJENH0EU65rcd4O2t4O1gbOxoRkdGICHYVXmZtJCKzxgZ6gI7+2IDCqz8hGe+j3fkhnLOMwqyJdnx3ORGNaEcvN6Do6k9Ixh60O89gbSQis8QGeoCqWm4i0lmLyGsf4Iz7v8N17Cg84Mj7+4hoZLvYdAOznLSIvLYXZ9zXszYSkVliAz1AzwR7ItVqKzqcwnDaMgazJnCFhYhoVag3Uke9ztpIRGaNDfQAyT9OwKruOP4+/iW42I6Cv/NoY0ciIjI6+cfXsK47xtpIRGaNDfRAXTmBWw6hOGn5CO/vIyLq8uMXuOUQwtpIRGaN30Q4QKqwjfhi1NNw6hyFQBeusBARAYBq+m+RPfpZ1kYiMmtcgR6g+tZ2lN6wRviEsbDgCgsREYDbtfHijdGsjURk1thAD1BxfSvsRqkwycXW2FGIiExGQR1rIxGZP97CMUDRnvYIcrCEpQVXWIiIuszxssckBwvWRiIya1yBHiCVSgV7a14+IqLuLFQqONpYGjsGEdF9pfgO8OjRo7C0tMSf/vQn3T6tVovk5GT4+fnB398fu3btMmJCIiIiIjInir6Fo6WlBb/61a+wYMECvf179+5FUVERfvjhBzQ1NWH69Ol49NFHMWnSJCMlJSIiIiJzoegV6JdeegkbNmzAuHHj9PZnZGRg9erVsLS0hIuLCxITE/HRRx8ZKSURERERmRPFrkB/9tlnuHbtGp566ikcOnRIb6yyshIajUa37ePjgzNnzvT7eE1NTRAR3baNjQ1sbGz6PL7r2O7nKInS8wPKnwPzG5/S5zDQ/IYcz9qoPEqfA/Mbn9LnMBy10WQb6IcffhjFxcW9jp07dw6//vWvkZ2d3ef53b/96l4uiLe3t972xo0bsWnTprued/PmzbseY8qUnh9Q/hyY3/iUPgdD87e2tt7zsayNyqX0OTC/8Sl9DvezNppsA33ixIk+x77++mvU1NRg5syZAID6+nocPHgQdXV1SElJgbe3NyoqKhAeHg4AuHTpUo8XgTtVVlbCwcFBt323VZbm5mZ4eXmhqqpK7zylUHp+QPlzYH7jU/ocBpq/vb39no9lbVQepc+B+Y1P6XMYjtposg10f6Kjo1FbW6vbTkpKwowZM7Bu3ToAwOLFi/HnP/8ZCQkJaGpqQkZGBj7//PN+H9PR0dHgfyQtLS0A9Fe7lUTp+QHlz4H5jU/pcxhIfkOOZW1UJqXPgfmNT+lzuN+1UZEN9N0sX74cp0+fRmBgIABgw4YNCA4ONnIqIiIiIjIHZtFAp6Wl6W1bWlpi586d93Ru1/3Rzc3NBj1n1/HNzc2K/O1M6fkB5c+B+Y1P6XMYaP6u8/p7fwhrozLzA8qfA/Mbn9LncD9rYxeVKPUtlkPk8uXL8PLyMnYMIqJhV1VVBU9Pz17HWBuJaKTqrzZ2GfENtFarRXV1Nezt7RX5WxYRkaFEBC0tLZg4cSIsLHr/OgDWRiIaae6lNnYZ8Q00EREREZEhFP1NhEREREREw40NNBERERGRAdhA96GtrQ3r1q1DQEAAQkJC8PTTTwMAamtrMX/+fAQEBCA0NBRff/217pwbN25g6dKl8Pf3R2BgIDIzM40VHwBw+PBhPPTQQ5g2bRpCQ0ORnp4OwHTn8MILL8DHxwcqlQoFBQW6/QPNq9VqkZycDD8/P/j7+2PXrl1Gyf/ss88iKCgIU6dOxZw5c3D+/HlF5e+Snp4OlUqFQ4cOmWT+/uYgItiyZQsCAwMRGhqK2NhYk5xDX/nPnDmDyMhITJs2DcHBwXjjjTeMlp+1kbVxqPKzNho3v1LqYn9zMGptFOrV+vXrJTk5WbRarYiIVFdXi4jIypUrZfPmzSIi8t1334m3t7e0t7eLiEhKSoqsWLFCRETKysrE3d1dGhsbhz27iIhWqxUXFxfJzc0VEZHy8nKxsbGR5uZmk53DsWPHpKqqSjQajeTn5+v2DzRvenq6PProo9LR0SENDQ2i0WikuLh42PNnZWXp8h48eFACAgJ0Y0rILyJSVVUlkZGREhERIQcPHjTJ/P3NYfv27ZKQkCBtbW0i8q//z6Y2h77yT506VbKyskREpKGhQVxdXaWwsNAo+VkbWRuHKj9ro3HzK6Uu9jcHY9ZGNtC9uH79ujg6OkpLS0uPsbFjx0ptba1uOzw8XI4cOSIiIpMnT5bvvvtON7Z48WJJTU2933F71fUicezYMRERyc3NlYkTJ0pbW5vJz+HO/yADzRsXFyf79u3TjW3YsEH3YnM/9VZku9TV1Ym1tbV0dnaKiHLyL1iwQE6ePCkxMTF6LxKmmF+k5xw8PDykpKSk12NNcQ69vUikp6eLiEhlZaV4eHhITU3NsOdnbWRtHAzWxlQRMZ38SquLIqZVG3kLRy9KS0uhVqvx+9//HjNmzMDDDz+Mr776Cg0NDdBqtXB1ddUd6+Pjg8rKSgBAZWUlNBpNr2PDTaVSYd++fUhISIBGo0F0dDTS09PR0tKimDkAGNQ1N7W5AMBbb72FuLg43cfjKCH/22+/jZCQEMyaNavHmBLyNzc3o66uDgcOHEBERAQiIiKQkZGhG1fCHFJTU/Hqq6/C29sbgYGBeP311zF+/Pi7Zhzq/KyNpjEHgLXRFPIruTaaQ10EjFsbzeKbCIdae3s7ysrKMHnyZPzhD39Abm4uHn/8cRQUFPT4PFS541MAu4/fOTacOjo68PrrryMrKwtRUVE4ffo0nnzySeTl5SlmDl0Gk9eU5rJ3717s27cPJ06c0NtvyvnLy8uxe/du5OTk9HmMKecHbv9/vnXrFlpbW3Hy5ElUVlYiMjISISEhCA0NvWtOU5jDtm3bsG3bNiQmJqKsrAyxsbGYOXMmgoKC7ppxKPOzNprGHLqwNvY+NhyUXhvNoS4Cxq2NXIHuhUajgYWFBZYtWwYACAsLwwMPPIDi4mIAQF1dne7YS5cuwdvbGwDg7e2NioqKXseG2/nz51FdXY2oqCgAQHh4OCZOnIi8vDwAypgDAKjVagADy2tKc8nIyEBKSgqys7Ph5uam22/q+b/99ltUV1cjODgYPj4+OHnyJH7xi19g9+7dd81oCvmB2/+G7OzsdG928/b2RlRUFM6cOXPXnKYwh/r6ehw4cACJiYkAAF9fX8yaNQvffPPNXTMOdX7Wxopex4yBtbH3seGi9Nqo9LoImEBtNOiGjxFk7ty58umnn4qISEVFhYwbN06qq6tlxYoVem/a8PLy0r0JYvPmzXo3rLu5uUlDQ4Mx4suVK1fE3t5evv/+exERKSkpEWdnZ7l8+bLJz+HOe5wGmjc1NVUee+wx3ZsEvL29paioaNjzZ2RkiL+/v1RUVPQ4Vgn5u7vzPj9TzC/Scw6rVq2SnTt3iohIY2OjaDQaOXv2rMnOoXv+jo4OcXZ2lqNHj4rI7XtFPT09dff2DXd+1kbWxqHKz9po3PxKq4t3zsHYtZENdB9KS0slJiZGQkNDJSwsTDIzM0XkdvGdO3eu+Pv7y+TJk3V/cSK332CTmJgofn5+EhAQIPv37zdWfBER+etf/yqhoaHy4IMPypQpU+TDDz8UEdOdw9q1a8XDw0MsLS3F3d1d/Pz8BpW3o6ND1q5dK76+vuLr6ys7duwwSn4rKyvx9PSUsLAw3U99fb1i8nd354uEKeXvbw51dXWycOFCCQkJkZCQEHnnnXdMcg595c/Ozpbp06fLgw8+KMHBwbJ9+3aj5WdtZG0cqvysjcbNr5S62N8cjFkb+VXeREREREQG4D3QREREREQGYANNRERERGQANtBERERERAZgA01EREREZAA20EREREREBmADTURERERkADbQREREREQGYANNRERERGQANtBERERERAZgA01ERAaJj4/H3r177/n41157DWvWrLmPiYiIhhcbaCJSPJVK1e/P0aNHjR1xyHh6eiItLc1oz5+Tk4PCwkIsXboUAFBRUQGVSoWLFy/qjiktLYVGo8HSpUvR0dGBdevWISMjA+Xl5caKTUQ0pNhAE5Hi1dTU6H7Wr1+PyMhIvX2zZ882dsS7unXr1rA9l1arRUdHx4DO3bVrF5YtWwZLS8texy9cuIA5c+bgkUcewQcffAArKys4Ojpi/vz5ePfddwcTm4jIZLCBJiLFGz9+vO5n7NixsLa21m27urrid7/7HTw9PWFvb4/Y2Fjk5eXpzt2yZQuio6Px1ltvYcKECXBycsLWrVvR1taG1atXw8HBAf7+/sjOztadk5aWBk9PT+zZswdeXl6ws7PDqlWr9JrgGzdu4Pnnn4erqyucnJywcOFCVFRU6MaTkpKwbNkyvPzyyxg3bhwWL14MAFi/fj18fX0xZswYhISEICMjQ3dObGwsfvzxR6xcuRIqlQqxsbG6/a+88oreNfHx8dE1rF2rxPv370d4eDhGjx6N/Px8dHZ24tVXX+3z2typvb0dn3zyCeLi4nodLygoQExMDBYtWoTU1FRYWPzrJSY+Ph779+/v87GJiJSEDTQRmbWUlBT87W9/w4cffohz584hKioKc+fORXNzs+6YvLw8nD9/HkeOHMH27duxadMmPPHEEwgJCcHZs2cxb948PPPMM3oNckNDA9577z0cOnQIBw4cwKeffoqtW7fqxlevXo2SkhJ89tlnOHXqFNzc3LBo0SJ0dnbqjsnKykJbWxtycnLw5ptvAgDUajU++ugjFBQUIDk5GcuXL0d+fj4AIDMzExMmTMD27dtRU1ODzMxMg67Fb37zG7z22msoKiqCr6/vPV2b7nJzc9Ha2opp06b1GDt37hxiY2OxZMkSvPPOO1CpVHrjM2bMQElJCWpraw3KTERkkoSIyIxs2rRJYmJiRESktbVVbG1tJT8/X++YgIAA2bNnj4iIbN68WZydneXmzZu68aCgIImPj9dt19TUCADJy8sTEZHU1FQBIMXFxbpjdu/eLWq1WkREysvLxdraWhobG3Xjt27dkjFjxsiJEydERGTFihXi6+srnZ2d/c5n3rx5kpKSotv28PCQ1NRUvWNiYmJk06ZNevs0Go3s3r1blweApKWl6cbv5drc6eOPPxZHR0e9fV2PbW1tLUuWLOlzHk1NTQJATp8+3ecxRERKYWXU7p2I6D4qLS1Fa2srIiIi9Pa3trairKxMtx0QEAAbGxvdtru7O0JCQvS2AaCurk63z97eHpMmTdJtz5w5Ew0NDWhoaEBhYSHa29vh5eXV6/NGR0cDAMLCwvRucwCA9PR07NixAxUVFbh58yba2tp6PM5AdV85vtdr011bW5vedeouPj4eBw8eRE5ODqKionqM29ra6h6fiEjp2EATkdm6fv06AODo0aNwcnLSG3NxcdH9edSoUXpjKpVKb1/X7QharbbHvr6e19bWFufPn+8x5ubmpvvzmDFj9MZOnDiBVatWYdu2bYiJiYGdnR2Sk5PR3t7e53MBgIWFBUREb19v53R/vnu9Nt2p1Wo0NTX1OvbGG29g/PjxiIuLw5dffonw8HC98cbGRgDAuHHj+p0LEZESsIEmIrMVHBwMa2tr1NTUYMaMGUP62M3Nzbhw4QKCgoIAAKdPn4ZarYZarUZYWBhu3LiB1tZWTJky5Z4f89SpU5g8eTJefPFFALcb9tLSUt0KOHC72e9+HzUAuLq64sqVK7rturo6ve3eDOTahIWFoa2tDeXl5XjggQd6jO/cuRM3b97EvHnzcOTIEYSFhenGioqKMGbMGAQEBNzTcxERmTK+iZCIzJaDgwPWrVuHNWvW4OOPP0Z5eTm+/fZbbNy4EYWFhYN67NGjR2PNmjXIzc3FV199hc2bN+P5558HAEyaNAkJCQlYsmQJDh8+jPLychw/fhzJycloaGjo8zH9/Pxw4cIFHDp0CBcuXEBycnKPRlij0eD48eO4cuWKbjV4zpw5yMzMxJEjR5Cfn4/nnnuuz1stBnNt3N3dMWXKFOTk5PQ6rlKp8O6772LBggWYO3cuioqKdGM5OTmIiYmBlRXXbYhI+dhAE5FZ27ZtG9auXYtf/vKXCAoKQmJiIqqqqqBWqwf1uGq1Gs888wzi4uLwxBNPYN68edi4caNu/IMPPsD8+fPx7LPPYtKkSUhKSkJ7e3uP2za6e/LJJ7Fq1SosX74cs2fPhr29PRYtWqR3zJYtW3Dq1Cl4eXnhZz/7GQDgueeeQ0JCAn7+858jPj4eS5cu1btVpC8DuTZJSUn9fhydhYUF3n//fcyZMwePP/44SkpKAAD79+/HypUr75qJiEgJVHLnjXNERNSvtLQ0vPLKK7h8+bKxowy7lpYWBAYGIicnB76+vvd0Tte93QUFBVyBJiKzwBVoIiK6Z/b29vjLX/5i0C8PTU1NeO+999g8E5HZYDUjIiKDLFiwwKDjFy5ceJ+SEBEZB2/hICIiIiIyAG/hICIiIiIyABtoIiIiIiIDsIEmIiIiIjIAG2giIiIiIgOwgSYiIiIiMgAbaCIiIiIiA7CBJiIiIiIyABtoIiIiIiIDsIEmIiIiIjLA/wMyLS0/KzE5AQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "gp.layout(colors=gp.colors.hollywood)\n", "\n", "fig, axs = plt.subplots(1, 2, figsize=(180*mm, 85*mm), sharex=True, sharey=True)\n", "\n", "for pressure, df in fO2_IW.groupby(\"P_bar\"):\n", "\n", " label = f\"{pressure} bar\" if pressure < 100 else f\"{pressure / 1e4} GPa\"\n", "\n", " axs[0].plot(df[\"T_K\"], np.log10(df[\"fO2\"]), \"-\", label=label)\n", " axs[1].plot(df[\"T_K\"], np.log10(df[\"fO2_magmapandas\"]), \"-\", label=label)\n", "\n", "\n", "axs[0].legend(frameon=True, fancybox=False)\n", "axs[0].set_title(\"Hirschmann matlab\")\n", "axs[1].set_title(\"MagmaPandas\")\n", "\n", "# axs[1].set_ylim(-0.1, 0.1)\n", "\n", "fig.supxlabel(\"Temperature (K)\")\n", "fig.supylabel(\"log10 $f$O$_2$\")\n", "fig.suptitle(\"IW\")\n", "\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "py310", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.10.17" } }, "nbformat": 4, "nbformat_minor": 5 }