"
],
"text/plain": [
" SiO2_Liq TiO2_Liq Al2O3_Liq FeOt_Liq MnO_Liq MgO_Liq CaO_Liq \\\n",
"0 49.1 3.22 14.4 14.8 0.14 3.20 6.72 \n",
"1 49.2 3.89 15.3 13.7 0.12 3.88 6.76 \n",
"2 49.6 3.79 15.8 13.0 0.14 4.26 6.59 \n",
"3 47.1 4.21 12.0 17.8 0.18 3.40 7.28 \n",
"4 48.1 3.88 13.2 16.4 0.16 4.02 6.51 \n",
"\n",
" Na2O_Liq K2O_Liq Cr2O3_Liq P2O5_Liq H2O_Liq Fe3Fet_Liq NiO_Liq \\\n",
"0 3.34 1.70 0.0 1.13 0.0 0.0 0.0 \n",
"1 3.44 1.22 0.0 0.83 0.0 0.0 0.0 \n",
"2 3.65 1.04 0.0 0.63 0.0 0.0 0.0 \n",
"3 2.93 2.02 0.0 2.32 0.0 0.0 0.0 \n",
"4 3.36 1.36 0.0 1.59 0.0 0.0 0.0 \n",
"\n",
" CoO_Liq CO2_Liq Sample_ID_Liq \n",
"0 0.0 0.0 Glass1 \n",
"1 0.0 0.0 Glass2 \n",
"2 0.0 0.0 Glass3 \n",
"3 0.0 0.0 Glass4 \n",
"4 0.0 0.0 Glass5 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"## Inspect them to check they read in right\n",
"display(Plags.head())\n",
"display(Liqs.head())"
]
},
{
"cell_type": "markdown",
"id": "c5cf7408-64a7-4bac-9117-0ef52dc8f19f",
"metadata": {},
"source": [
"## Please check your feldspars are actually plagioclases! If you apply plag-liq thermometers to kspars, you get stupid results, that mess up the iterative calculations"
]
},
{
"cell_type": "markdown",
"id": "1fc10d10-9af6-4ca4-afbc-b88b05691c57",
"metadata": {},
"source": [
"## Lets plot to inspect and check they are actually plag, else you will have convergence issues\n",
"-This function relies heavily on the ternary plot package from Marc Harper et al. 2015 - https://github.com/marcharper/python-ternary, if you use these figures, you must cite that (Marc Harper et al. (2015). python-ternary: Ternary Plots in Python. Zenodo. 10.5281/zenodo.594435) as well as Thermobar.\n",
"\n",
"You may have problems with this package if you have the separate “ternary” package installed (yes, there are python packages called ternary and python-ternary- Yay!). I (penny) got the error “module ternary has no attribute figure, so had to uninstall the ternary I had through pip (pip uninstall ternary), and re-install python-ternary through conda in the command line “conda install python-ternary”. If you have everything in pip, or conda, keep in 1 environment, don’t follow my bad example here!"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "859ffb88-f57c-4b71-8036-47bb3f26738e",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
""
]
},
"execution_count": 4,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAhEAAAIDCAYAAAC6gEyCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAACx2UlEQVR4nOzdd1RUV9cH4N8MHUTsIlYUDUXsBVtM1Kixt1iiUWM0MfbeFXsDFQtijyX23rvYUURQqiAiCNJ7nX6/P3iZTyIqzNyZO2U/a71rJcy95+yXwLBn33PO5jEMw4AQQgghpIz4XAdACCGEEO1ESQQhhBBCFEJJBCGEEEIUQkkEIYQQQhRCSQQhhBBCFEJJBCGEEEIUQkkEIYQQQhRCSQQhhBBCFEJJBCGEEEIUQkkEIYQQQhRCSQQhes7b2xvDhg1D7dq1YWJigkqVKqFjx47YsmULBAIB1+ERQjQYj3pnEKKfJBIJJk+ejD179sDCwgI///wz7OzskJWVhVu3buHdu3do1KgRrl69Cjs7O67DJYRoIEoiCNFTc+fOhbu7O1q3bo3z58+jZs2a8tekUilWrlyJlStXws7ODi9fvkT58uU5jJYQoonocQYheujt27fYvHkzKlWqhMuXLxdLIADAwMAAK1aswK+//orIyEi4u7vLX6tXrx7q1auHzMxMTJs2DbVr14ahoSEOHjyo5v8XhBCuURJBiB46ePAgZDIZ/vzzT1SvXv2L1y1duhQAcODAgWJfFwqF6NKlC65fv46+fftiypQpXx2HEKKbDLkOgBCifk+fPgUAdO3a9avX2dvbw8bGBh8/fkRsbCxq164NAEhMTESTJk3w5MkTmJmZqTxeQohmokoEIXooMTERAORJwdcUXZOQkFDs625ubpRAEKLnKIkghHxV0dprHo8n/5qpqSmcnZ25CokQoiEoiSBED1lbWwMAYmNjv3ltXFxcsXsAoFq1asWSCkKIfqIkghA91L59ewDA3bt3v3rdmzdvEB8fj5o1axZ79EEJBCEEoCSCEL00ZswY8Pl87N27FykpKV+8bs2aNQCAcePGqSs0QogWoSSCED3UqFEjTJ8+HWlpaejbt+9niyZlMhlWrVqFf//9Fw0aNMCcOXM4ipQQosloiychemrjxo3IysrCgQMH0LBhQ/Tu3RsNGjRAdnY2bt26hbdv36Jhw4a4du0anVZJCCkRHXtNiJ67c+cO9uzZgydPniAlJQUWFhZwcHDAkCFD8Pfff3+2jbNevXoAgOjoaPUHSwjRKJREEEIIIUQhtCaCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQQgghCqEkghBCCCEKoSSCEEIIIQqhJIIQUqKMjAz8/PPPePLkCdehEEI0FI9hGIbrIAghmsfFxQXPnz+HpaUlsrKywOPxuA6JEKJhqBJBCPlMZmYmQkNDAQA5OTl49uwZxxERQjQRJRGEkM9s3LgROTk5sLKyAgAsWLAAVLQkhPwXJRGEkGISExPh4eEBAFi7di2MjIzw8OFD3L59m9vACCEah5IIQkgxq1atQkFBAVq2bInff/8df/75JwBg4cKFVI0ghBRDCysJIXLv379Hw4YNIZVKcePGDXTu3BnJyclwdHREXl4eTp06hV9++YXrMAkhGoIqEYQQuWXLlkEqleLHH39E586dAQDVqlXDtGnTAABLly6FRCLhMkRCiAahSgQhBAAQFBSEJk2aAAAeP36Mli1byl/Lzs6Gvb09MjIysG/fPvzxxx9chUkI0SBUiSCEACjcgQEAAwcOLJZAAED58uUxb948AIXVCIFAoPb4CCGahyoRhBA8ffoUHTp0AI/Hg7+/P+zt7T+7pqCgAA4ODkhKSsKmTZswa9YsDiIlhGgSqkQQoucYhsH8+fMBAKNHjy4xgQAAMzMzLFu2DEDhDo6cnBy1xUgI0UyURBCi527duoXHjx/D2NgYixcv/uq1o0ePRv369ZGZmYlNmzapKUJCiKaiJIIQPSaTyeRrIf766y/Url37q9cbGhpixYoVAAB3d3ekpKSoPEZCiOaiJIIQPXbmzBm8evUKFhYWmDt3bqnuGTRoEJo2bYq8vDysWbNGxRESQjQZJRGE6CmJRCJf4zBjxgxUrVq1VPfx+XysXLkSAODl5YXY2FiVxUgI0WyURBCipw4ePIjw8HBUrlwZ06dPL9O9P/30Ezp16gSRSITly5erJkBCiMajLZ6E6CGBQIAGDRogPj4eGzZskJ9IWRbPnj3Djz/+CB6Ph5CQEDg4OKggUkKIJqNKBCF6yNPTE/Hx8bC2tpY32CorFxcX/Pzzz2AYBkuXLmU5QkKINqBKBCF6Jjs7G3Xq1EFWVha8vLwwduxYhccKCgpCmzZtAAAvXrxAq1atWIqSEKINqBJBiJ7ZtGkTsrKy0KBBA4waNUqpsZydnTFs2DAAha3CCSH6hSoRhOiRlJQU1KtXD/n5+Th69CgGDRqk9Jjh4eFo1qwZAODo0aP49ddflR6TEKIdqBJBiB5ZvXo18vPz0bx5cwwcOJCVMe/duyf/Z3d3d9DnEkL0B1UiCNETHz58gJ2dHcRiMS5fvoxu3bopPWZubi6cnJyQnJwMQ0NDSCQSnD9/HgMGDFA+YEKIxqNKBCF6wtXVFWKxGJ07d0bXrl1ZGdPT0xPJycmwtbXF1KlTARS2FJdKpayMTwjRbFSJIEQPhIWFwcnJCQzD4P79+2jbtq3SY6anp8PR0RFZWVn4559/0KNHD9jb2yM7OxuHDh3C6NGjWYicEKLJqBJBiB5YsmQJGIZB7969WUkggP/f5eHs7IyhQ4eiYsWKmDNnDgBg2bJlEIlErMxDCNFcVIkgRMe9ePFCfpaDn58fnJyclB7z48ePaNy4MQQCAc6ePYtevXoBAPLz8+Hk5ITExERs27ZN/oiDEKKbqBJBiI5btGgRAGDEiBGsJBAAsH79eggEAvmplUXMzc3lrcXXrFmDvLw8VuYjhGgmqkQQosPu3r2Lbt26wdDQEMHBwahbt67SY7579w7NmjWDRCLB7du30bFjx2Kvi0QiNG3aFNHR0Vi9ejUWL16s9JyEEM1ElQhCdBTDMJg/fz4A4I8//mAlgQCAlStXQiKRoEePHp8lEABgbGwsbzG+ceNGpKenszIvIUTzUCWCEB11/vx5DBo0CCYmJggPD0f16tWVHvP169dwcXEBAPj4+MhPqvwvmUyG1q1bIzQ0FHPnzsXGjRuVnpsQonmoEkGIDpJKpfIqxPTp01lJIABg+fLlAIAhQ4Z8MYEAAD6fj5UrVwIAtm7divj4eFbmJ4RoFkoiCNFB//77L96+fYvy5ctj5syZrIz55MkT3LhxAwYGBnB1df3m9b169UKrVq0gEonkCQUhRLdQEkGIjhEKhVi6dCkAYN68eahQoYLSYzIMI1/nMHbsWNjZ2X3zHh6Ph7Vr1wIA9u3bh8jISKXjIIRoFkoiCNExu3fvRmxsLKpXr46///6blTFv3ryJp0+fwtTUtEwtvzt16oRu3bpBKpXSLg1CdBAlEYTokNzcXKxevRoAsHjxYpibmys9pkwmkz++mDhxImrWrFmm+4seZZw6dQqvXr1SOh5CiOagJIIQHeLh4YGUlBTUr18fY8eOZWXM06dPIzAwEOXLl5cfa10WzZs3x5AhQwBAfhAVIUQ30BZPQnREWloabG1tkZOTg4MHD2LYsGFKjykWi9GsWTNERUXB1dVV4SQgMjISTZs2hUwmw8OHD9GpUyelYyOEcI8qEYToiPXr1yMnJwdOTk745ZdfWBnz0KFDiIqKQrVq1TBlyhSFx7Gzs8OYMWMAFFYj6LMLIbqBKhGE6ICPHz+ifv36EIlEOHfuXLF+For6tJmWu7s7Jk+erNR4cXFxaNy4MYRCIa5cuYLevXsrHSMhhFtUiSBEB6xcuRIikQht2rRBz549WRnTy8sLiYmJqFOnDsaPH6/0eLVq1cLEiRMBFDYFk8lkSo9JCOEWVSII0XJv376Fvb09ZDIZ7ty5gw4dOig9ZmZmJhwdHZGRkYG9e/di1KhRLERauG7DwcEBOTk5OHbsGEaMGMHKuIQQblAlghAtt3jxYshkMvz000+sJBBA4S6PjIwMODg4sPqHvnLlypgxYwaAwrjFYjFrYxNC1I8qEYRosYCAALRo0QIA8OzZMzRt2lTpMRMTE+Hk5IT8/HycOHEC/fv3V3rMT+Xm5qJRo0bIyMiAl5eX/BEHIUT7UCWCEC1WtOVy6NChrCQQALBhwwbk5+ejVatW6NevHytjfqpcuXLy0ytdXV1RUFDA+hyEEPWgSgQhWurhw4fo3LkzDAwM8Pr1azRo0EDpMaOjo9GkSROIxWJcv34dP/zwg/KBlkAoFMLR0RHx8fHYsGED5s2bp5J5CCGqRZUIQrQQwzDyKsTYsWNZSSAAYPXq1RCLxejSpYvKEggAMDExkR+lvWbNGmRmZqpsLkKI6lAlghAtdOXKFfTt2xcmJiYICQkpcz+LkoSEhKB169ZgGAaPHj1Cq1atWIj0y6RSKVq1aoU3b95g0aJFWLNmjUrnI4SwjyoRhGgZmUwmX1MwefJkVhIIAFi+fDkYhsGAAQNUnkAAgIGBgbwa4eHhgaSkJJXPSQhhFyURhGiZEydOyBtizZ49m5Uxnz9/jitXroDP58v/sKtD//790apVK+Tn52PVqlVqm5cQwg5KIgjRImKxGEuXLgUAzJw5E5UqVVJ6TIZh5InDqFGjYG9vr/SYpcXj8bBixQoAwO7du/H+/Xu1zU0IUR4lEYRokX379iEqKgqVK1dWqiHWp+7du4cHDx7A2NhY/phEnbp06YIff/wREolErVUQQojyKIkgREvk5+dj+fLlAAp7T5QrV07pMRmGwbJlywAAEyZMQJ06dZQeUxFF1YgjR44gODiYkxgIIWVHSQQhWmL79u1ITk5GzZo1WWmIBQDnz5+Hv78/ypUrx+lZDa1bt0bfvn0BAEuWLOEsDkJI2VASQYgWyMzMxNq1awEU7qIwNjZWekyJRCKvAEybNg3VqlVTekxlrFixAjweDxcvXsSzZ884jYUQUjqURBCiBTZu3Ijs7GzY29uz1hDr6NGjiIiIQOXKlTF9+nRWxlSGg4MDfv31VwDA/PnzQUfYEKL56LApQjRcYmIi6tevj4KCApw8eZKVfhYCgQCNGzfGx48fsW7dOnlnTa59+PABTk5OkEgkuHnzJrp37851SISQr6BKBCEabtWqVSgoKECbNm3k6waUtWfPHnz8+BE1a9bEX3/9xcqYbKhTp448nnnz5lE1ghANR5UIQjTY+/fv0ahRI0gkEtYaYuXk5MDR0RGpqanYuXMnfv/9d+UDZVFycjLs7e1RUFCA06dPY8iQIVyHRAj5AqpEEKLBli1bBolEgq5du7LWEGvbtm1ITU1Fw4YN8dtvv7EyJpuqVasmX6OxaNEiSCQSjiMihHwJVSII0VBBQUFo0qQJALDWECslJQWOjo7Izc3FkSNHNPZTflZWFhwdHZGeno59+/bhjz/+4DokQkgJqBJBiIYqOi+hqL8EG9zc3JCbm4vmzZtj0KBBrIypClZWVpg7dy6Awi2tAoGA44gIISWhSgQhGujp06fo0KEDeDwe/P39WelnERsbC2dnZwiFQly6dAk//fQTC5GqTkFBAZycnJCQkIDNmzdj5syZXIdECPkPqkQQomEYhsHChQsBsNsQa+3atRAKhejUqRO6devGypiqZGZmJu/lsWbNGuTk5HAcESHkv6gSQYiGuXnzJnr27AkjIyMEBwez0s8iPDwcLVq0gEwmg7e3N1xcXFiIVPUkEgmaNWuGd+/ewdXVVd47hBCiGagSQYgGYRhG3sNi4sSJrDXEWrFiBWQyGfr06aM1CQQAGBoayhMHd3d3pKamchsQIaQYqkQQokFOnz6NoUOHwszMDOHh4ahatarSY/r7+8vXV7x48QJOTk4sRKo+MpkMbdq0QUhICGbOnInNmzdzHRIh5H+oEkGIhpBIJPK1EDNnzmQlgQAAV1dXAMDw4cO1LoEAAD6fjzVr1gAAduzYgdjYWI4jIoQUoSSCEA1x8OBBvHv3DpUqVWKtIdaDBw9w584dGBoaYunSpayMyYXu3bujQ4cOEIvF8qSIEMI9SiII0QACgUD+7H/+/PkoX7680mMyDCNPHP744w/Y2toqPSZXeDweVq9eDaAw2Xrz5g3HERFCAEoiCNEIO3fuxMePH1GrVi38+eefrIx55coVvHjxAubm5liwYAErY3LJxcUFvXr1AsMw8q2fhBBuURJBCMeys7Oxdu1aAIW9IkxNTZUeUyqVyisbkydPhrW1tdJjaoKi/0/nzp2Dn58ft8EQQiiJIIRrmzZtQlpaGqsNsU6cOIHQ0FBUqFBBp056dHZ2xogRIwBAvgiVEMId2uJJCIdSUlJga2uLvLw8HD16lJV+FiKRCE2aNEFMTAxWrVqFOXPmsBCp5nj//j2cnZ0hlUpx7949/Pjjj1yHRIjeokoEIRxau3Yt8vLy4OzsjIEDB7Iy5v79+xETEwNra2tMmjSJlTE1ia2tLcaNGwegsBpBn4MI4Q5VIgjhyIcPH9CgQQNIJBLWGmLl5eXByckJSUlJ2Lp1K2uLNDVNQkICHB0dIRAIcOHCBfTv35/rkAjRS1SJIIQjy5cvh0QiQYcOHVhriOXp6YmkpCTY2tpi7NixrIypiWrUqIEpU6YAABYvXgypVMpxRIToJ6pEEMKBsLAwODk5gWEY3L9/H23btlV6zPT0dDg6OiIrKwv//PMPhg8fzkKkmisjIwP29vbIzs7GoUOHMHr0aK5DIkTvUCWCEA4sXrwYDMOgT58+rCQQQOEuj6ysLDg7O2Po0KGsjKnJKlasKF80unTpUohEIo4jIkT/UCWCEDV78eIF2rRpI//nxo0bKz1mfHw8nJycIBAIcPbsWfTq1UvpMbVBfn4+GjVqhLS0NGzfvl3+iIMQoh5UiSBEzYpOj/z1119ZSSAAYP369RAIBHBxccHPP//MypjawNzcXH60t6urK/Ly8jiOiBD9QkkEIWp079493Lt3DzweD2PGjMGmTZsgFAqxfv16AIXJQHx8PA4ePAg/Pz/cvn0bFy5cQFhYGDw9PZGdnV3s2vT0dKxatQoHDhwAAPTu3RuvXr3Cvn37kJKSUuza/Px8eHh4ICIiAmfOnIG3tzd8fHxw5MgRxMbGFrtWKpXCzc0NMTExOHbsGJ48eYKHDx/i1KlTiIyMxObNmyEQCIrdk5iYiAMHDsDPzw+3bt3CxYsXERISgp07dyIrK6vYtZmZmfDy8kJwcDAuXbqEGzduwN/fH/v370dSUlKxawsKCuDh4YHIyEicPn0a9+/fx9OnT3H06FF8+PABqampqFatGtLT0+Hh4aG2/5aEEHqcQYjaMAwDFxcX+Pr64s8//0S9evVYOU1yzJgxOHXqFHr06IELFy4oH6iWOXPmDOLj4+WNy96/f49KlSpxHRYheoGSCELU5Pz58xg0aBDMzMwQFhaGSpUqQSgUoly5cgqPGRgYKF+Y6ePjg2bNmrEUrXYoKCgAn8+HgYEB2rVrh+DgYMydOxcbN27kOjRC9AI9ziBEDaRSKZYsWQIAmDJlCqpXr46srCxcunRJqXGLGlINGTJE7xIIAAgJCcHLly9haGgIV1dXAMC2bdsQHx/PcWSE6AdKIghRg3///VfeEGvWrFkAgCpVqqBDhw7w9/dXaMwnT57g+vXrMDAwkP8B1SfJyckQiURo3749gML1IC4uLhAKhVi5ciXH0RGiHyiJIETFhEIhli1bBgCYM2cOKlSoIH/NysoK1apVK/OYDMPIxxw7dizs7OxYiVWbmJubo2bNmvJ/5/F48uRh3759ePfuHVehEaI3KIkgRMV2796NDx8+oEqVKvj777+LvVahQgUEBwcjICCgTGPevHkTT58+hampqV62xM7Ly8P+/ftRt27dYl/v1KkTunbtCqlUKt/6SQhRHVpYSYgK5ebmom7dukhPT/9iQyyGYRAfHw9ra2sYGBh8c0yZTIZ27dohMDAQM2bMwLp161QRukb78OEDateuDR6P99lrAQEB8kccAQEBerlWhBB1oUoEISrk4eGB9PR01K1bF7///nuJ1/B4PAQFBSEjI6NUY54+fRqBgYEoX768/NhnfcIwDB4/flxiAgEAzZs3x6BBgwAUHi9OCFEdSiIIUZH09HT5VsMVK1bAyMjoi9f27NkTV65c+Wb/B7FYLH/uP3PmTFSuXJm9gLXE/v37MWLEiK9es2LFCvD5fFy7dg2PHz9WU2SE6B9KIghRkXXr1iEnJwfOzs745Zdfvnl97969v5poAMChQ4cQFRWFatWq6WWfCIZh0L9//y9WIYrY2dlhzJgxAIC5c+eCntoSohqURBCiAh8/fsT27dsBFJ7lwOd/+1etatWqXz0kKT8/H2vXrgUAzJs3T6lDqrTV9u3bYWFhUaprFy1aBCMjIzx79gzXr19XcWSE6CdKIghRgRUrVkAoFKJdu3Zlaog1Z86cL54b4eXlhYSEBNSpUwfjx49nK1StERERgTFjxsDc3LxU19eqVQuTJ08GUFiNkMlkqgyPEL1ESQQhLIuMjJQ3xFq1atU3S++f4vP5SE9P/+zrmZmZ2LRpEwBg6dKlMDExYSdYLZKVlVXmRGDOnDmwtLREaGgoTp48qaLICNFflEQQwrKlS5dCKpWie/fu6NChQ5nu5fF4aNmyJfbv31/s6x4eHsjIyICDg8M3FxXqoocPH4JhGFSsWLFM91WuXBkzZswAUPjfRSwWqyA6QvQXnRNBCIsCAgLQokULAMo1xBIKhZDJZDAzM0NSUhKcnJyQl5eHEydOoH///ixGrPmkUinEYjEMDAy+ufC0JDk5OXByckJKSgp27dqFv/76SwVREqKfqBJBCIuKziUYPHiwUoccRUVF4eHDhwCADRs2IC8vD61atUK/fv3YCFOrfPjwATdv3lQogQAAS0tLzJ8/HwCwcuVKFBQUsBkeIXqNKhGEsOTRo0f4/vvvYWBggNevX6NBgwZKjRcSEoL09HT07t0bYrEY169fxw8//MBOsFpCIBDg0aNH+Omnn5QaRygUwtnZGbGxsdiwYQPmzZvHUoSE6DeqRBDCAoZhMHfuXADAmDFjlE4ggMK+Gnv27IFYLEaXLl30LoEAChea1qlTR+lxTExM5K3Y169fj6ysLKXHJIRQEkEIK65evYrnz5/D2NgYixYtYmXMrKwsnD17FkDhllF9I5PJsGXLFnz33XesjDdy5EjY2dkhIyMDbm5urIxJiL6jJIIQJclkMvkz98mTJxdrT62M5cuXg2EY9OjRA61atWJlTG2SkJDA6mMHAwMDrFq1CgCwadMmJCUlsTY2IfqKkghClHTixAmEhobC0tKStYZYz58/x+XLl8Hn89GpUycIhUJWxtUmd+/eZX3M/v37o1mzZhAIBFi9ejXr4xOibyiJIEQJYrFY/qx91qxZqFSpktJjMgwDV1dXAMCoUaMwe/ZsHD16VK/6Pxw7dgxDhw4t00FdpcHj8bBmzRoAwK5duxAdHc3q+IToG0oiCFHCvn378P79e1YbYt27dw8PHjyAsbGxfMtor169WBlbGzAMgy5dusDU1FQl43fp0gU//vgjJBIJli5dqpI5CNEXlEQQoqD8/Hx5W+758+ez0hCLYRgsW7YMADBhwgT5zoRq1arpzWLAI0eOQCKRqHSOooWq//77L4KDg1U6FyG6jM6JIERBGzZswIIFC1C7dm0EBQWx0s/i3LlzGDlyJMqVK4eQkBBUq1ZN/ppAIEBUVBQcHR2VnkdTxcXFwcjICNWrV1f5XMOHD8fFixfRt29fXLp0SeXzEaKLqBJBiAIyMzOxfv16AMCyZctYSSAkEon8E/K0adOKJRBFdH1HQW5uLvLz89Uyl6urK3g8Hi5fvoxnz56pZU5CdA0lEYQoYOPGjcjMzIS9vT1rDbGOHj2KiIgIVK5cGdOnT//sdVNTU9SuXRsXL15kZT5NEx4ejnfv3sHW1lYt8zk4OGDkyJEAgIULF+rVwlVC2EKPMwgpo8TERNja2kIgELDWEEsgEMDZ2RlxcXFYt26dvPPkf0mlUkgkEvB4PBgbGys9ryYRCAQAoLIFlSWJiYmBs7MzxGIxbt68ie7du6ttbkJ0AVUiCCmj1atXQyAQoHnz5qw1xNqzZw/i4uJgY2Pz1S6TBgYGePbsGYKCgliZV1NkZWXh2LFjak0gAKBu3bqYMGECAGDRokVUjSCkjKgSQUgZvH//Hg0bNoRUKmWtIVZOTg4cHR2RmpqKnTt34vfff//mPT4+PrC3t0fFihWVnp9rDMPg5s2b6NmzJyfzJycnw8HBAfn5+Th9+jSGDBnCSRyEaCOqRBBSBkuWLIFUKsWPP/7IWkOsbdu2ITU1FQ0bNsRvv/1WqnuqVKkCQ0NDVubnGsMwqFGjBmfzV6tWDdOmTQNQ2Mpd1dtLCdEllEQQUkpBQUE4duwYAMjPh1BWSkoKtm7dCqBwl0dpE4OGDRti//79KCgoYCUOLrm5uaFp06acxjBjxgyUL18eEREROHz4MKexEKJNKIkgpJQWLlwIABgwYABrDbHc3NyQk5ODZs2aYdCgQWW6d/r06Vrf0joxMZHVJluKsrKykv/3XbJkiXyRJyHk6yiJIKQUfHx8cPXqVfB4PHlfC2XFxsZiz549AAorG3x+2X4deTwerl+/rtWLAe/cuQOxWMx1GACAv/76C9WrV0dCQgK8vLy4DocQrUBJBCHfwDAMFixYAAD47bffYG9vz8q4a9euhVAoRKdOndCtWzeFxhg7diyOHDnCSjzqdu3aNXTv3l1jtqqamZnJe2msXr0aOTk5HEdEiOaj3RmEfEPRzgEjIyMEBwfL+1koIzw8HC1atIBMJoO3tzdcXFwUHuvjx4+wsbFhveOlKjEMg6SkJFStWhUGBgZchyMnkUjQrFkzvHv3DsuXL2et6kSIrqJKBCFfwTAMFi1aBAD4888/WUkggMIGUDKZDL1791YqgQAKn+fv2LGDlbjU5c6dO4iLi9OoBAIADA0NsXz5cgCAu7s7UlNTuQ2IEA1HlQhCvuL06dMYOnQoLCwsEBYWhqpVqyo9pr+/Pzp06AAejwdfX180btxY6TGzs7ORlpamtiOjlZGZmYmsrCzUrVuX61BKJJPJ0L59e7x+/RozZ87E5s2buQ6JEI1FlQhCvkAikWDJkiUACndCsJFAAJCXyIcPH85KAgEUHhmdnJzMyliqlpubi/T0dK7D+CI+ny9vhObp6YnY2FiOIyJEc1ESQcgXHDx4EBEREahQoUKJDbEU8fDhQ9y5cweGhobyRXxsqFatGhiGwdOnT1kbUxXS09Nx9+5dNG/enOtQvqp79+7o0KEDRCKRPKEghHyOHmcQUgKBQABbW1skJiZi/fr1rCQRDMPghx9+gK+vL/766y94eHgoH+gnBAIBeDwejIyMyrxdVF1EIhEkEgnMzc25DuWbfHx80KVLF/B4PISGhrK2K4cQXaKZ7zSEcGznzp1ITExE9erVv9oQqyyuXLkCX19fmJubY/78+ayM+SlTU1NcunRJY8vvUqkUO3fu1IoEAgDatWuHHj16gGEYVqtGhOgSqkQQ8h/Z2dmoV68eMjIySt0Q61ukUinatGmD0NBQzJ07l7Vjs0vi7e2NTp06aVxvjdu3b+Onn37iOowyCQoKQps2bQAAfn5+aNmyJccREaJZqBJByH9s2rQJGRkZZWqI9S0nT55EaGgoKlSogJkzZ7Iy5pdUrlwZUqlUpXMookKFClyHUGbOzs4YOnQoAMgPHCOE/D9KIgj5REpKCjZt2gSgcBcFG5/mRSKRvPIwa9YslbfvbtKkCbZu3apRx2G7u7ujWbNmXIehkOXLl4PP5+POnTvw9vbmOhxCNAo9ziDkEzNmzMDWrVvRrFkzPHnyhJUFil5eXpg1axasra0REhKiljUBDMMgNTWVtW2pykhNTUXFihU17mCpspg+fTr27NmDli1b4sWLF1p1OighqkSVCEL+58OHD/LGS4o0xCpJXl4eNmzYAKCwC6g6FxVeunRJbXN9zaNHjzT6XIjSWLBgAUxMTPDy5UuN+b4SogkoiSDkf5YvXw6RSKRUQ6z/8vT0RFJSEmxtbTF27FhWxiwNHo+H4cOH49y5c2qbsyQ+Pj5wdnbWiIqIMmrUqIGpU6cCKEwGNXHNCSFcoMcZhAAICwuDk5MTGIbB/fv30bZtW6XHTE9Ph6OjI7KysnDgwAGMGDGChUhLj2EYxMbGonbt2pyU3xmGQUZGBoyMjGBpaan2+dmWkZEBBwcHZGVl4fDhw6wtuiVEm1ElghAAS5cuBcMw6NmzJysJBFC4yyMrKwvOzs4YNmwYK2OWBY/Hg0wmw4kTJ9Q+N1C4PdLHx0cnEggAqFixImbPng0AWLZsGUQiEccREcI9qkQQvffixQv5WQAvXrxgpZ9FfHw8nJycIBAIcPbsWfTq1UvpMRWVkpICmUyG6tWrq21OsViM9+/fo1GjRmqbUx3y8/Ph4OCA5ORkbN++HVOmTOE6JEI4RZUIovcWLlwIgN2GWOvXr4dAIICLiwt+/vlnVsZUVFpaGjIyMtQ6Z15eHhISEtQ6pzqYm5vLf15WrVqFvLw8jiMihFuURBC9du/ePdy9excGBgby7prKevfuHf755x8AhX9ouN4OaG9vj8jISERFRallPrFYjGPHjqFz585qmU/dxo0bhzp16iA5ORlbt27lOhxCOEWPM4jeYhgGrVq1gr+/P6sNscaMGYNTp06hR48euHDhAitjKisvLw/GxsYwMjJS+VxisRhCoRDlypVT+VxcOX78OMaNG4dy5cohJiYGlSpV4jokQjhBlQiity5evAh/f3+YmpqydqRxYGAgTp06BaBwy6imsLCwwP79+5GTk6PyubZu3arTCQQADBs2DPb29sjNzZWfA0KIPqIkguglqVQqf7Y9depUWFtbszJuUeIwZMgQjTvmeeLEifDz81PpHA8fPsSMGTNUOocm4PP5WLVqFQDAw8NDJ9d/EFIalEQQvfTvv//izZs3qFChAmbNmsXKmE+ePMH169dZXV/BNjMzM8hkMpWNb2xsrNXHW5dF79690bp1a4hEIqxYsYLrcAjhBCURRO+IRCL5H/nZs2ez0l2SYRgsW7YMQOGaCDs7O6XHVIW2bdvCzc1NJWPv3bsXDRs25HwhqbrweDysWbMGALBv3z68e/eO44gIUT9aWEn0zvbt2zFt2jRWG2LduHEDAwcOhKmpKYKCglCrVi0WIlUNmUyGzMxMVhcDZmVlwdzcXC0LNzVNv379cPv2bQwbNoyzg70I4QpVIoheyc3NlX96ZKshlkwmk1c2Jk6cqNEJBAAIhUJcvnyZ1TFfv36NyMhIVsfUFkWPMk6ePInXr19zHA0h6kVJBNErHh4e8oZYv//+OytjnjlzBoGBgShfvjzmzJnDypiqZGZmhp49e+LevXusjBcREQFzc3M4ODiwMp62ad68OQYPHgwAWLRoEcfREKJe9DiD6I309HTUq1cPOTk5+OeffzB8+HClxxSLxWjWrBmioqKwbNky+Y4PTScSiZCcnMxK1SQ3NxcFBQVa36lTGW/fvkWzZs0gk8nw6NEjdOzYkeuQCFELqkQQvbF+/Xrk5OTAwcEBQ4cOZWXMQ4cOISoqClWrVpW3itYGxsbGiIuLU7oakZSUhFOnTul1AgEADRs2xOjRowEUPiajz2ZEX1AlguiFjx8/on79+hCJRKw1xCooKICTkxMSEhLg7u6OyZMnsxCpesXGxqJixYoKHw4VHh6O7777juWotFNcXBycnJwgEolw9epVTpuuEaIuVIkgemHFihUQiURo06YNaw2xvLy8kJCQgDp16mD8+PGsjKluHz9+RGZmpkL3ikQiREREsBuQFqtVqxYmTpwIoLAaocrzOAjRFJREEJ0XGRmJAwcOAADWrFnDyjkGmZmZcHd3BwAsWbIEJiYmSo/JBRcXF3h7eyMrK6vM9+7duxd9+/ZVQVTaa968eShXrhwCAwNx8uRJrsMhROUoiSA6b/HixZBKpejevTtrC948PDyQkZEBBwcH/Prrr6yMyZX+/fvDwsKiTPdIJBL89ttvKopIe1WuXBkzZ84EULhTQywWcxwRIapFSQTRaQEBAfKGWGwdTZyUlIQdO3YAAFxdXbX+mOfy5cvDw8OjTOX3rVu3ljnx0BdTp05FxYoVER0dLa+AEaKrKIkgOq1oy+Uvv/zCWkOsDRs2IC8vD61atUK/fv1YGZNrs2fPxuPHj0t1rZ+fH/766y+tT55UxdLSUn5exLJly1BQUMBxRISoDiURRGc9evQIN2/eBJ/Pl/e1UFZMTAz27dsHAFi5cqXO9Ing8Xjg80v3diCTyWBsbKziiLTbhAkTULNmTSQnJ8urVoToIkoiiE5iGEZehRg7dixrDbFWrVoFsViMLl264Mcff2RlTE3RunVreHh4fPWaixcvonz58pREfIOJiYk8cV23bp1CC1cJ0QZ0TgTRSVevXkWfPn1gYmKCkJAQ1KxZU+kxQ0ND0apVKzAMg0ePHqFVq1YsRKpZJBIJ8vLyYGVl9dlrBQUF4PF4MDY2LnXVQp9JpVK0atUKb968weLFi7F69WquQyKEdfROQHSOTCaTP5P++++/WUkgAGD58uVgGAb9+/fXyQQCAFJSUnD37t0SX4uMjISvry8lEKVkYGAgr0Zs2bIFSUlJHEdECPvo3YDonBMnTiAwMBDlypVjrSGWr68vLl++DD6fj+XLl7MypiaqUaMGmjRpgoCAgGJfT01NRXJyMr7//nuOItNOAwYMQIsWLZCfny/vHkuILqEkgugUsViMxYsXAwBmzZqFypUrKz0mwzDyT5SjRo2Cvb290mNqssqVK3/WmMvS0hINGzbkKCLtxePxsGrVKgCFJ5xGR0dzGxAhLKMkguiU/fv3Izo6GpUqVWKtIda9e/fw4MEDGBsbyxMUXVaxYkU8efIEYWFhAAChUIgdO3agTp06HEemnbp06YLvv/8eEokErq6uXIdDCKtoYSXRGfn5+ahXrx5SUlKwadMmTJo0SekxGYZBx44d4e/vj8mTJ8uPutYH79+/R926dfH+/Xs0aNCA63C0mq+vLzp37gwej4egoCA4OTlxHRIhrKBKBNEZ27dvR0pKCmrWrIk//viDlTEvXLgAf39/WFhYYN68eayMqS0iIiKQnZ2Nly9fch2K1mvTpg369OkDhmH0oppF9AdVIohOyMzMRN26dZGdnY09e/aw0tdBIpGgZcuWiIiIwKJFi7B06VIWItUuo0aNwj///AMjIyOuQ9F6oaGhaNmyJQDAx8cHLi4uHEdEiPKoEkF0gpubG7Kzs1G5cmW4uLhg06ZNEAqFWL9+PQBg/fr1iI+Px8GDB+Hn54fbt2/jwoULCAsLg6enJ7Kzs4tdm56ejgkTJiAiIgKWlpZwdHREQEAA9u3bh5SUlGLX5ufnw8PDAxEREThz5gy8vb3h4+ODI0eOIDY2tti1UqkUbm5uiImJwbFjx/DkyRM8fPgQp06dQmRkJDZv3gyBQFDsnsTERBw4cAB+fn64desWLl68iJCQEOzcuRNZWVnFrs3MzISXlxeCg4Nx6dIl3LhxA/7+/ti/fz+SkpKKXVtQUAAPDw9ERkbi9OnTuH//Pp4+fYqjR4/iw4cPWLduHVatWgU3NzfIZDJs2LABHz58wNGjR/H06VPcv38fp0+fRmRkJDw8PFBQUFBs/KSkJOzfvx/+/v64ceMGLl26hODgYHh5eSEzM7PYtVlZWdi5cydCQkJw8eJF3Lp1C35+fjhw4AASExOLXSsQCLB582ZERkbi1KlTePjwIZ48eYJjx44hJiYGbm5ukEqlxe6JjY3FkSNH4OPjA29vb5w5cwYRERHw8PBAfn5+sWtTUlKwb98+BAQE4Pr167hy5QqCgoKwa9cupKenF7s2Ozsbnp6eCAsLw4ULF3D79m34+fnh4MGDiI+PL3ZtgwYN4OjoCAB6V9UiuosqEUTrJSUlwdbWFgUFBdi0aRPq1q2L3r17KzWmQCCAs7Mz4uLisG7dOsyYMYOdYLWIl5cXRo4cid27d2Pu3Llch6P1oqKicOfOHcyePRsSiQS3bt3CTz/9xHVYhCiFkgii9aZMmQJPT0+0bt0ad+/ehUwmk5+sqKjt27dj3rx5sLGxQXBwMMzMzFiMWPO9efMGVapUQZUqVSAWixEQEIA2bdpwHZZWy8/Ph4GBARYtWoSdO3eiWbNm8Pf315n+K0Q/0eMMotXev3+P3bt3Ayhs9W1kZITHjx8jPDxc4TFzcnKwceNGAMCSJUv0LoEACteDFP1xMzQ0hEAg4Dgi7Zabm4ujR4/CxMQE8+fPh7m5OV69eoWzZ89yHRohSqEkgmg1V1dXSCQS/PDDD/KGWF27dkVmZiays7MVGnPbtm1ITU1Fw4YNWVmgqW2ePXuGjIwM+UFdPB4P3333HQ4dOsRxZNqJYRg8ePAAEyZMAABUq1YN06ZNA1CYpEokEi7DI0Qp9DiDaK2goCA0adIEAPDw4UO0bt1a/lpkZCRsbGxgbm5epjFTU1Ph6OiInJwcHDlyBEOGDGE1Zk0nk8lQUFAAExMTGBoayr/OMAzEYjHEYjEsLCw4jFD7MAyD4OBgODs7y7+WlZUFBwcHZGRkYP/+/Rg3bhyHERKiOKpEEK21ZMkSAEC/fv2KJRAAYGdnBy8vL4hEojKN6ebmhpycHDRr1gyDBg1iLVZtkZCQgJs3bxZLIIDCakRYWBidGaGAjRs3onHjxsW+ZmVlJV+sumLFCnpcRLQWVSKIVvLx8UH79u3B4/Hw8uVLODg4fHYNwzBISEiAjY1NqcaMjY2Fs7MzhEIhLl68iO7du7MdtkYTiUS4efMm+vbt+8Vr/Pz8UKVKFdSrV099gWmx+Ph41KhRo8TFkwUFBXByckJCQgI2b96MmTNnchAhIcqhSgTROgzDYP78+QCAkSNHlphAAIWfnm/fvl3qcdeuXQuhUIhOnTrp5dY7Pp//xe9lkRo1aqBSpUpqikj7PXz4EEKhsMTXzMzM5KdXrl69Gjk5OeoMjRBWUBJBtM7t27fx6NEjGBoafvMUyVGjRpVqQWBERAQOHz4MAFi5cqXebbtjGAbu7u6ws7P76nU1a9bEiRMnkJqaqqbItNf169fx/fffw9TU9IvXjBkzBvXq1UN6ejo2b96sxugIYQc9ziBahWEYNGvWDIGBgZgyZQrc3Ny+eU9iYiKqV6/+1cRg5MiROHfuHHr37o0zZ86wGbJWiI6ORt26dUudPEVHR9Mjja9gGAZpaWmwsrL65pHhp0+fxujRo2Fubo6YmBhUqVJFTVESojyqRBCtcubMGQQGBsLc3LzURwebmppi165dX3zd398f586dA4/Hw/Lly1mKVLv4+PiU6XpfX1/amvgVDx48wNu3b0vVc2Tw4MFo3Lgx8vPzsW7dOjVERwh7KIkgWkMikch3ZEyfPh1Vq1Yt1X0VKlTA8OHDER0dXeLrrq6uAIDhw4d/topeHxw9ehR9+/Yt0yOcoUOHYv/+/SqMSntlZ2ejdu3aaNeuXamu5/P5WL16NYDCk1JjY2NVGR4hrKIkgmiNQ4cOISIiApUrVy5zLwuBQICMjIzPvv7w4UPcuXOnVOsrdJFMJkPPnj1Rrly5Mt87ZMgQyGQyFUSl3fLz85GVlVWme7p3746OHTtCLBbrbTWMaCdKIohWEAgE8jfXuXPnonz58mW6v0aNGsjNzcWLFy/kX2MYRp44/PHHH7C1tWUtXm1x5swZZGZmKnSvlZUVLQb8j+zsbFy5cgUtWrQo0308Hg8rV64EABw8eFCpY9sJUSdaWEm0wubNmzF79mzY2NggJCTkqyvev0QoFIJhGJiYmIDH4+HKlSv45ZdfYGZmhpCQENSoUUMFkWuuuLg4MAyD2rVrKzxGXl4eoqKiip3GqM+KDo1S5OcTKFwfce3aNQwZMgSnT59mMzRCVIIqEUTj5eTkYO3atQCAxYsXK/wGbWJiggsXLiAhIQFSqVS+FmLKlCl6l0AAhWtMxGKxUmPw+Xw63+B/ZDIZdu3apfDPJwB5te3MmTN0OijRCpREEI23adMmpKWloX79+hg9erRSYw0fPhwhISE4fvw4QkNDUaFCBb08KTAmJgYvXrxA/fr1lRrHzMwMVlZWuHbtGkuRaa9bt25h+vTpSo3h7OyM4cOHAwAWLVrERliEqBQ9ziAaLSUlBfXq1UN+fj7+/fdfDB48WOkxX758iVGjRiE6OhorV66U9zDQJ3l5eTA2Ni7VFsRvkUql8tbhxsbGLESnnfz9/cu8FqIkUVFRaNKkCaRSKby9vfHDDz8oHxwhKkKVCKLR1q5di/z8fDRu3BgDBw5kZcwXL14gOjoa1atXx+TJk1kZU5vk5eXh+PHjrCQQAGBgYIBHjx4hMjKSlfG00ebNm+Hk5MTKWPXr15d39Vy4cCHocx7RZFSJIBorNjYWDRo0gFgsZq0hVl5eHpycnJCUlIRVq1Zhzpw5LESqPRiGwaVLl9C/f3/Wx/b29kbLli3LvHNG26WkpKBSpUowMDBgbcyEhAQ4OjpCIBDgwoULKvnvRQgbqBJBNJarqyvEYjGrDbE8PT2RlJQEW1tbWFlZsTKmtmnYsKFKxq1Tpw5r1Q1t8vz5c6SkpLA6Zo0aNeRVsgULFkAqlbI6PiFsoSSCaKSwsDAcPHgQALBq1SpWGmJ92uRo6dKlGD58OM6ePav0uNpkw4YN3+zUqagGDRrA09NT6R0f2uT58+ews7ODtbU162PPnj0b5cqVw5s3b3Ds2DHWxyeEDfQ4g2ikQYMG4fz586w2xFqyZAk2bdqExo0b4/nz5+DxeIiPj4eNjY1edO2Mi4uDjY0N+HzVfXZgGAZxcXFKnT2hLRiGQU5ODng8HiwtLVUyx8aNG+Hq6opatWrh3bt3er1wlWgmqkQQjePn54fz588DAGtHAMfHx8PT0xMAsGLFCvD5fPB4POTn5+tN104fHx8IhUKVzsHj8fDgwQOVzqEp3rx5g3v37qksgQCAyZMno0qVKoiLi8OePXtUNg8hiqIkgmichQsXAgBGjBjBWkOs9evXQyAQwMXFBT///LP86w0bNkTHjh2RnJzMyjya6tq1a2jXrh3MzMxUPteIESNw4MABlc/DJalUCplMhgEDBqh0HgsLCyxevBgAsHLlSuTl5al0PkLKipIIolHu3bvHekOsqKgo/PPPPwBKXl+RmZmJ7OxsVubSRAzDwMXFBdWrV1fLfAYGBujTp49Ob03My8tDamqqWuYaN24c6tWrh5SUFGzdulUtcxJSWpREEI3BMIy8CjFu3DjWGmKtXLkSEolE3inxvxwcHBAYGIiYmBhW5tM0jx49QlhYmFp3ThgZGels+V0qleLQoUPo3LmzWuYzNjaWJ9QbN24ssRstIVyhhZVEY1y4cAEDBw6EmZkZQkNDWVnxHhgYiLZt2wIoXBPQrFmzEq/Lz8+HoaGhzi1cy8rKQmJiIr777ju1z52WloasrCylj9bWNCKRCBKJBObm5mqbUyqVwsXFBcHBwZg/fz7Wr1+vtrkJ+RqqRBCNIJVK5c9+J0+ezNqWuaKFmUOGDPliAgEA5ubm2LdvH/Lz81mZV1MIhUIUFBRwMrdIJEJubi4nc6vStm3b1JpAAIWPiIoaxm3duhUJCQlqnZ+QL6EkgmiEo0ePIjQ0FOXLl8fs2bNZGfPp06e4fv16sTfgr5k0aRKePn3KytyaICcnBxcuXPhq8qRKNWrUQHp6uk51o3zw4AGmTZvGydy9e/dG69atIRAIsHLlSk5iIOS/KIkgnBOJRPIqxNy5c1GhQgWlx2QYBsuWLQMAjBkzBnZ2dqW6z8rKSmcWBBoaGmLMmDGcxtC+fXs4OjrqzPfU3Nycs1M5eTweVq9eDQDYu3cv3r17x0kchHyKkgjCud27dyMuLg5VqlTBpEmTWBnz1q1bePLkCUxNTeWLNUujVatWcHNzYyUGLslkMuzZswcmJiacxmFsbIwLFy4gKSmJ0zjY8M8//6BOnTqcHkz2/fffo0uXLpBKpfIkmRAu0cJKwqm8vDzUq1cPqamp8PDwwF9//aX0mDKZDO3atUNgYCBmzJiBdevWlfn+9PR0VKlSRelYuHL16lX06tVLY07ivH79Orp3785qkyp1ysrKgqmpKedJGQAEBASgffv2AIBXr16hadOmHEdE9BlVIginPDw8kJqaCltbW3n7Y2WdOXMGgYGBKF++vEJdOgsKCnDr1i1WYuFK7dq1NSaBAAqbc2nz55XQ0FBERERwHQYAoHnz5hg0aBCAwuZchHCJkgjCmfT0dGzcuBEAsGzZMlaeNYvFYvmisxkzZqBy5cplHsPCwgKdO3eGt7e30vFwYcuWLWjUqBHXYRTj5OQEd3d3rUwkitYeODs7cxzJ/1u+fDn4fD5u3LiBx48fcx0O0WP0OINwZu7cuXB3d4ezszOePXvGSmOo/fv3Y8qUKahatSpCQkIU7msgEAiQmZmpku6MqpSUlITKlSvD0NCQ61A+wzAM4uPjUbNmTa5DKZOCggLk5OSgWrVqXIdSzKRJk/DPP//AxcUFT58+1ajKE9EfVIkgnIiPj8f27dsBAK6urqwkEAUFBVizZg0AYP78+Uo1RjI1NUVkZKTWfcoLCAhASkoK12GUiGEY3Lt3j+swyiQtLQ1HjhzRuAQCABYtWgRjY2M8e/YM169f5zocoqcoiSCcWLlyJYRCIVxcXNCrVy9WxvTy8kJCQgLq1KmD8ePHKz1ex44dUatWLa05gOr58+eoXbs2atSowXUoJeLz+ejXr59WdU1NSkrCn3/+yXUYJapVqxb+/vtvAIVN62QyGccREX1ESQRRu8jISOzbtw9AyQ2xFJGVlYVNmzYBAJYsWcLaKvqEhAStaM7FMAycnJxQp04drkP5qvLly6Njx45asTZCIpEgKiqK6zC+as6cObC0tERgYCBOnTrFdThED1ESQdRu6dKlkEql6Nq1a4kNsRTh4eGB9PR02Nvb49dff2VlTABo164drl69qvEtmN++fYubN28q9QhHHXg8HrKysnDu3DmuQ/mmnTt3ok+fPlyH8VVVqlTBjBkzABQuThaLxdwGRPQOLawkavXq1Ss0b94cQOGx1EX/rIykpCQ4OTkhLy8PJ06cQP/+/ZUe81O5ubkwMTHh7KTCb5FKpQgJCUGTJk24DqXUPn78CENDQ7W1Jy8rsVgMoVCIcuXKcR3KN+Xk5MDR0RGpqanYtWsXK2etEFJaVIkgarVo0SIAwODBg1lJIABgw4YNyMvLQ6tWrdCvXz9WxvxUuXLlsHXrVo0twRcUFCArK4vrMMokJydHo9ea7NixQyMOlioNS0tLzJ8/HwCwYsUKzhquEf1ElQiiNo8ePcL3338PPp+P169fl7qfxdfExMTA2dkZYrEY165dw48//shCpJ9jGAb3799X2fiKkslk2L59O6ZPn851KGV26tQptGvXDrVr1+Y6lGJevHiB7777DuXLl+c6lFITCoVwdHREfHw8Nm7ciLlz53IdEtETVIkgasEwjPyNbezYsawkEACwevVqiMVi/Pjjjyr9A8/j8WBqaqqy8RUlEom0tnzdr18/VK1aleswPmNgYKA1VYgiJiYm8k61a9as0brKFNFelEQQtbh27RqeP38OY2Nj+SMNZYWGhuLo0aMAoJbWyE2bNsW2bdtUPk9ZeHl5ad0fvCKmpqb4559/IBAIuA5F7urVqzAyMtLK7+nIkSNhZ2eHrKwsnWgiR7QDPc4gKieTydC0aVMEBwcr1BDrS4YOHYrLly+jf//+OHHiBCtjfotEIkFubi4r7cqV5e3tjfbt22vlH7xP3bx5Ez169OA6DBQUFIBhGJiamrJy+BkXzp8/j19//RWmpqaIjo7W2IWrRHdo528K0SonTpxAcHAwLC0tFWqIVRJfX19cvnwZfD4fy5cvZ2XM0khISNCYUyytrKxgbGzMdRhKq1q1qkYsWo2JiYGvr6/WJhAAMGDAALRo0QICgQCrV6/mOhyiB7T3t4VoBbFYjKVLlwIAZs6cqVBDrP9iGAbLli0DUFjCtbe3V3rM0qpduzYaNGiAV69eqW3Okhw6dAjW1tY60S+hWbNmnJffMzIyEB0djR9++IHTOJTF4/Hkj/Z2796NmJgYjiMiuo4eZxCV2rVrF/7++29UrVoVoaGhrOy7v3v3Lvr06QNjY2MEBgaibt26LERaeunp6eDxeKhYsaJa5y2SmZkJU1NTjVzoqSiZTIaUlBTOyu9isRiJiYkat1NEUb169YK3tzd+++03HD58mOtwiA6jSgRRmYKCAqxatQpAYUMsNhKIT6sQEyZMUHsCAQCVKlXCnTt38PbtW7XPDQDh4eGcza0q+fn5uH//Pidzi8VibN26VWcSCADyR3z//vsvQkJCuA2G6DRKIojKbN++HfHx8ahVqxYrDbEA4MKFC/D394eFhQXmzZvHypiK+OWXXwBA7U2PIiMjIRaL4ezsrNZ5Va1cuXJo3749J10+379/z9paHU3Rpk0b9OvXDwzDYMmSJVyHQ3QYJRFEJTIzM7F27VoAhWf6s7GDQCKRYMWKFQCAadOmcd6eOSoqSu2nLtauXRvfffedWudUl6pVq6JJkyZqX2QZFBSk1vnUxdXVFTweDxcuXMDz58+5DofoKEoiiEq4ubkhKysLdnZ2rDXEOnbsGMLDw1G5cmWNOKGxR48eOHz4sNqqEenp6fjnn3808oAmNpiamiIkJAQ+Pj5qm3PXrl0YMGCA2uZTJ0dHR4wYMQIAWDubhZD/ooWVhHVJSUmoV68eBAIBjh8/zsqbtEAggLOzM+Li4rBu3Tp550KuZWdnw9zcHIaGhiqfKzAwUKuabCkqMjISNjY2MDc3V+k8EokE+fn5WnW8dVl9eiz8rVu38NNPP3EdEtExVIkgrFu1ahUEAgGaN2/OWkfNvXv3Ii4uDjY2Nhp1zLOlpSW2bNmi8nmkUini4uJUPo8mSE1NRW5ursrn+eeffyCRSFQ+D5fq1q0rX4+0cOFCjTiPg+gWqkQQVr1//x6NGjWCRCJhrSHWp62OPT09MW7cOBYiZY9IJIK/vz9cXFxUNsfWrVs14hGOuuzZswejRo1SWTUiPDwcFSpU0IsTHZOTk2Fvb4+CggKcOXMGgwcP5jokokOoEkFYtXTpUkgkEnTp0oW1hljbt29HamoqGjZsiNGjR7MyJpsMDQ0hlUpVNr5IJMKECRNUNr4mGjVqlMofEWnzyZRlUa1aNXkCunDhQp2vvhD10o/fIqIWwcHB8oZYbB1FnZqaCg8PDwCFuzzUsfagrPh8PurVq4d///1XJePv2rULRkZGKhlbU5mbm2PHjh0qKb+/fPkSsbGxOrtAtSQzZsyAlZUV3r59S4dPEVbR4wzCmv79++PSpUsYMGAAjh8/zsqY8+fPx7Zt29CsWTM8efJEYz89MgwDkUgEsVjMyqFaRXx9ffHdd9/BysqKtTG1BcMwuHfvHrp27cramDKZDLm5ubCwsICBgQFr42qDLVu2YNGiRbCxsUFUVJTWN24jmkEz35GJ1vHx8cGlS5fA4/HkJ0oqKzY2Frt37wYArFixQmMTCKCwZ0FISAjrZw7o2vHWZcHj8VjfOZGSkoJbt27pXQIBABMnTkSNGjUQHx8PLy8vrsMhOkJz35WJ1mAYRr4PfeTIkXBwcGBl3HXr1kEoFKJTp05asTWtRYsW4PF4+PDhAyvjXb9+HQzD6PUnRkdHR2zbto2VsSQSCZ4+fYohQ4awMp62MTMzw+LFiwEAa9asQU5ODscREV1AjzOI0m7duoUePXrAyMgIwcHBqFOnjtJjRkREoHnz5pDJZPD29lbpzgc2ffz4EVZWVko/0igoKIBUKoWFhYVOdOpUhkQiQXZ2NipVqqTUOFKpFNHR0WjQoAFLkWkfiUSC5s2bIzIyEsuXL4erqyvXIREtR5UIohSGYbBw4UIAhQ2x2EgggMLHFzKZDL1799aaBAIAatasiSNHjiAzM1OpceLi4uDn56f3CQRQmJj5+voqPY6bm5teJxBA4U6ioseN7u7uSE1N5Tgiou2oEkGUcvr0aQwdOhTm5uYICwtjpZ9FQEAA2rdvDx6PB19fXzRu3JiFSNUrKioK9evXV+jezMxMPH78GH369GE5Ku0VHBwMsViM5s2bK3R/VFQUbG1tKSlD4eJSFxcXBAUFYdasWdi0aRPXIREtRpUIojCJRCJfCzF9+nTWGmIVlViHDRumlQkEULiNUNGzI8qVK6fwH0tdVbNmTaWqCC9fvqTTGv+Hz+dj1apVAIAdO3bozUmoRDUoiSAKO3ToECIjI2FlZcVaL4uHDx/i9u3bxcqu2uiXX37Bnj17ynyfRCLBli1bULNmTRVEpb0qVqyIa9euISoqqsz3Hjt2DD/99JNG7+5Rt+7du8PFxQUikYi1M12IfqLHGUQhAoEAtra2SExMZK0hFsMw+OGHH+Dr64u//vpLfsiUtsrMzISlpWWZthOGhobC0dFRhVFptzdv3qBRo0alTgikUilycnJgZWVFjzL+w8fHB126dAGfz0doaKjOtpgnqkWpOVHIzp07kZiYCGtra9YaYl29ehW+vr4wMzPD/PnzWRmTS+bm5mXenhgREaGiaHRDbGwsCgoKSn395cuXkZiYSAlECdq1a4fu3btDJpNhyZIlXIdDtBRVIkiZ5eTkoF69ekhPT2etIZZUKkXbtm0REhKCOXPmyJ/Zarvs7Gx5O+Zv2blzJ/7880+NPNpbk2zfvh1Tpkz5ZmLw8eNHiEQi2Nraqiky7RMYGIi2bdsCAPz8/NCyZUuOIyLahioRpMw2bdqE9PR0VhtinTx5EiEhIahQoQJmzZrFypiawMDAoFSfnMViMUaPHk0JRCn8/vvvpWoiJZPJIJPJ1BCR9mrSpAmGDx8OAFiwYAHH0RBtREkEKZOUlBT5ljC2GmKJRCJ55WHWrFmoWLGi0mNqCgsLCxgbG+P27dtfve7ff/+FQCBQU1TazcLCAlu3bv3qNXFxcXj06JHenwtRGkuXLoWBgQHu3LmD+/fvcx0O0TL0OIOUycyZM+Hh4cFqQ6xdu3Zh5syZsLa2RkhICMzNzVmIVHNIpVKIRCIYGBjA2Nj4s9ffvHkDKysr1KhRg4PotJNQKMTLly/Rvn37El/PycmBmZkZVXZKadq0adi7dy/atGmDZ8+e0RoSUmpUiSClFhsbi507dwJgryFWXl4e1q9fD6CwnKprCQRQ+Ejj4cOHiI6OLvF1Q0ND+mNXRkZGRl/8+RMIBDhx4gR9T8tg4cKFMDU1ha+vLy5fvsx1OESLUBJBSm358uUQiUTo0KEDaw2xPD09kZSUBFtbW/z++++sjKmJevTogejoaOTm5hb7ekBAAN69e4eqVatyFJl24vP5sLGxwbFjxz577fr165gwYQIHUWmvGjVqYMqUKQAKEwpFD0oj+oceZ5BSefPmDRwdHcEwDGsNsdLT0+Ho6IisrCwcOHAAI0aMYCFSzRUVFYWaNWvKu3LKZDJkZ2ejfPnydBCSAhiGgUgkglgsljc8YxgG4eHhsLe35zg67ZORkQEHBwdkZWXh8OHD+O2337gOiWgBeucipbJ06VIwDIOePXuy1hBr8+bNyMrKQuPGjTFs2DBWxtRk9evXx7Zt2+Sf8tLS0nDnzh1KIBTE4/EQHByM0NBQ+dc2btyIRo0acRiV9qpYsaJ8Z5SrqytEIhHHERFtQJUI8k1+fn5o3bo1AODFixes9LNISEiAk5MTCgoKcObMGfTu3VvpMbVFTEwMatasifPnz+OXX37hOhyt9/jxY3n32Jo1a5bphFBSXF5eHuzt7ZGamoodO3Zg8uTJXIdENBx9BCLfVHR65PDhw1lriLV+/XoUFBTAxcUFvXr1YmVMbfH06VPw+Xy0adOG61B0QoMGDVC1alX4+/vTNlklWVhYYPHixQCAlStXIi8vj+OIiKajJIJ8lbe3N+7duwcej8daQ6yoqCgcOHAAALBq1Sq92042ZMgQDBs2DHXr1uU6FJ1Qo0YNzJ07F3Z2drCwsOA6HK03btw41KxZE8nJyd88j4MQSiLIFzEMg7lz5wIoPGefYRhs2rQJQqFQvi1z/fr1iI+Px8GDB+Hn54fbt2/jwoULCAsLg6enJ7Kzs4tdm56ejt9++w0SiQTNmjVDTk4OAgICsG/fPqSkpBS7Nj8/Hx4eHoiIiMCZM2fg7e0NHx8fHDlyBLGxscWulUqlcHNzQ0xMDI4dO4YnT57g4cOHOHXqFCIjI7F582YIBIJi9yQmJuLAgQPw8/PDrVu3cPHiRYSEhGDnzp3Iysoqdm1mZia8vLwQHByMS5cu4caNG/D398f+/fuRlJRU7NqCggJ4eHggMjISp0+fxv379/H06VMcPXoUHz58wKJFi7Bjxw6sW7cOMpkMGzZswIcPH3D06FE8ffoU9+/fx+nTpxEZGQkPDw8UFBQUGz8pKQn79++Hv78/bty4gUuXLiE4OBheXl7IzMwsdm1WVhZ27tyJkJAQXLx4Ebdu3YKfnx8OHDiAxMTEYtcKBAJs3rwZkZGROHXqFB4+fIgnT57g2LFjiImJgZubG6RSabF7YmNjceTIEfj4+MDb2xtnzpxBREQEPDw8kJ+fX+zalJQU7Nu3DwEBAbh+/TquXLmCoKAg7Nq1C+np6cWuzc7OhqenJ8LCwnDhwgXcvn0bfn5+OHjwIOLj4z/7fteoUQNxcXE4efIkHj58iMePH+P48ePyuCUSSbF74uLicPjwYTx//hx3797F2bNnER4ejm3btiE3N7fYtampqdi7dy9ev36Na9eu4erVq3j9+jX27NmDtLS0Ytfm5ORg+/btCAsLw/nz53Hnzh34+vri0KFDn8UtFovh7u6O9+/f4/jx43j8+LE87vfv38Pd3R1isfiz37VDhw7B19cXd+7cwfnz5xEWFobt27cjJyen2LVpaWnYs2cPXr9+jatXr+LatWt4/fo19u7di9TU1GLX5ubmYtu2bQgPD8fly5cxaNAgAMC6deuQkZGh1PsI0W20JoJ80YULFzBw4ECYmZlh7dq16NevH2xsbJQaMygoCG3btgXDMPDx8UGzZs3YCVaLnDt3Du3bt8fVq1fxxx9/cB2O1vPx8YFQKER6ejoGDBhAC1WVJJPJsHXrVhw+fBhv3rzBvHnzsGHDBq7DIhqKkghSIqlUiiZNmiA0NBRz5szB4sWLwePx5NsTFTV48GBcu3YNQ4YMwZEjR1iKVnscPXoUvXr1QsWKFZGcnIy8vDxqEKWE3NxcxMTEwMnJCQzDYOfOnbQYUElCoRAMw+DOnTv45ZdfYGJigvfv39OJqqRElLKTEh09ehShoaHyhlimpqbYt28fhEKhwmM+ffoU165dg4GBAVxdXVmMVjtIJBL07dsXFSpUAFDYdIsWAipHKBTKtyLyeDz5ozKiOE9PT5iYmKB3795o27YthEIhVqxYwXVYRENREkE+IxKJ5H/kP22INXnyZIUb9DAMI1+YOWbMGNjZ2bESqza5fv064uLi5AtJa9asidjYWLx+/ZrjyLRTfn4+Tpw4gebNm8u/ZmJiAk9PTw6j0m7379/HpEmTwOPxwOPxsHLlSgDAvn378O7dO46jI5qIHmeQz2zfvh3Tpk1D9erVERoaWqyfxcuXL9GiRYsy76i4efMmBgwYABMTEwQHB6NWrVpsh63RPn78CIFA8FlXyaITF83NzfVul4qy8vLyYGRk9FlTs8zMTMTExKBp06YcRaa9/P390bx582I/i3379sWdO3cwYsSIEo8ZJ/qNKhGkmLy8PKxevRpA4Rn6/22I1bx5c7i5uZVpTJlMJq9s/P3333qXQBQpKUkwNjbGxYsXkZqaykFE2othGOzfv7/ErqhGRkbU+0EBR44cQbVq1T77OS16lHH8+HEEBgZyERrRYJREkGI8PDyQnJyMOnXqlNgQi8/nY86cOUhOTi71mGfOnMHr169Rvnx5zJkzh81wtUJ8fDzu3buH+vXrl/j6r7/+iufPn0Mmk6k5Mu115coVTJ06tcTXLCws5AsDSelkZmZiyJAhJSb4LVq0kG/5LDqIipAilEQQuU/36i9fvrzET3lA4Yr40q6NEIvF8ueqM2bMQOXKlVmJVZtYWFjg119//eo1DRo0AD1ZLD1bW9uvPv5p1qwZ2rdvT/0fSunt27cIDw//4uuurq7g8/m4cuUKHj9+rMbIiKajJILIFR064+Dg8NWGWOXLl0fbtm1LlUgcPnxY3uq6qNWwPhEKhThz5sw3+zk4ODiU+TGRvtq2bds3t8UaGBjg0aNHiI2NVVNU2uv9+/coKCj46pktjRo1knf1XLRoESW8RI4WVhIAhcdbd+/eHRKJBAsWLEDLli0hk8nAMAxkMpn8f0X/LhQKIRAIYG5uXuyaT/9ZIBBg7dq1yM7Ohru7u17u3z9z5gyGDBlSqmsZhsHHjx/1ds1IaSQlJaFixYpfrJL91/Xr19GpUyd5q3DyOaFQiOzsbFStWvWr18XFxcHJyQkikQjXrl3Dzz//rKYIiSajJIIAABwdHREWFqaSsU1NTZGYmKj0QVXahmEYREZGomHDhqW6XiqV4uTJk9989KHP7t27h0aNGpU60YqOjkbNmjVhZGSk4si0U1ZWFo4ePYpJkyaV6vr58+dj27ZtcHZ2xqtXr+h0UEJJBAEiIyPRqFEjMAyDRo0aoXz58uDz+eDxeODz+fL/lfTvBQUFMDc3h5GRkXxvedHrUVFReP36NZycnODn58f1/02127hxI+bMmVOmN9r09HR4e3tj8ODBKoxMO7148QJGRkZlPip948aNmD17NrUIL0FQUBCcnZ1LfX1qaiq+++475Ofn4/jx4xg+fLgKoyPagJIIgmHDhuHUqVPo3r07Ll68WKZ7fXx80LBhQ1SpUuWz19zc3LBs2TKMHj0au3fvZitcrRATE4OaNWvC0NCwTPcxDIPU1FRUqVKFzo34BMMwEAgEkEgksLS0LPP97969++yMDn0nk8lw48YN9OrVq0z3rVmzBqtXr0a9evUQERFBVR49R7UoPffq1SucOnUKQOGOjLJq164dzp49W+Lxzbm5uQCgl8+jg4KCFDoinMfjISkpCVeuXFFBVNorOjoaFy9eVCiBAKCXlbBv2b59e5kTCACYNm0aKleujOjoaPzzzz8qiIxoE0oi9FzRvu8hQ4YUOz64LEaNGlViyV5fk4gbN27AyckJFhYWCt3fuHFjNGnSpExncegymUyGjIwMpUrnAwcOxL59+1iMSruJRCKMHz9eoXstLS2xYMECAIUfPAoKCtgMjWgZSiL02KNHj3Dt2jXw+XylGmJZWFhg+/btn2370sckQiqVomPHjqhZs6ZS4+Tn51Nzrv8RCATIy8tTagxjY2MMHjyYTrL8n927d5f5UdunJkyYgNq1ayMhIQE7duxgMTKibSiJ0FMMw2DhwoUA2GmINXPmTHh7exf7mj4mES9fvsSLFy9KvQXxSxwcHPDo0SMkJiayFJl2YhgGu3btQqdOnZQeSyQS6WX7+f/y8/PDqFGjlNotZWJiIq9irlu3DllZWWyFR7QMJRF66tq1a3jy5EmxNwNl8Pn8z5KFoiRC0efY2iY3NxcmJib48ccfWRlvyJAhsLKyYmUsbSUQCFg7X6R69ero0aMHoqKiWBlPWxkbG8PU1FTpcUaOHIlGjRohIyMD7u7uLERGtBElEXpIJpPJqxATJ05UuvRexMHBoVhpsyiJUHRtgLYRi8WsjmdiYoIDBw7o9dHNu3btUrqq8ymZTAaJRMLaeNrm1q1bEIvFMDMzU3osQ0ND+WLszZs30xoePUVJhB46efIkgoKCYGFhgblz57I2rqWlJSZOnIiMjAwA+lWJEAgE+Pfff1lvPz158mS9bSR17949TJw4kdWtrjVr1sS7d+8QHBzM2pjaQiAQoF27dgovoC7JgAED0Lx5c+Tn58u7/xL9QkmEnhGLxVi0aBEAYPbs2aw3xIqLi8OLFy8A6FclQiwW4++//1bJ2GxVirRNlSpVWCm7/9dPP/2EunXr6l3/h48fP8LX15fVUyZ5PB5WrVoFAPDy8kJMTAxrYxPtQEmEntm/fz+io6NRuXLlL7ZSVka9evVgbW2NwMBAvalEMAyDQ4cOKbXa/WsaN26sd825jh49igoVKqjkwC1DQ0NcvnwZ6enprI+tqbKzsxEaGoquXbuyPnbXrl3RsWNHSCQSpXZ5Ee1EJ1bqkYKCAtja2iIpKUmlDbHS09NhaGgIW1tb5OfnIywsDPXq1VPJXJrgwoUL6Nevn0r7CMhkMiQnJ8Pa2lplc2iKzMxMGBsbw9zcXKXzqOO/m6aQSCRITExUWXM3X19fdO7cGTweD0FBQXByclLJPETz6P5vD5Hbvn07kpKSULt2bYUPmimNSpUq4fLly8jPzweg+1s87e3tVf6HKCcnB0+fPlXpHJri/fv3iIiIUPk8jo6OKp9DE0ilUmzZskWl3WHbtGmDvn37gmEY+eNSoh8oidATmZmZWLduHQBgyZIlKu+o2bdvX/k/63ISsWPHDrWsWbCyskLz5s1x//59lc/FpejoaGRmZpa5yZYiGjVqhI0bN6p8Hq5FRESwuoD6S5YvXw4ej4dLly7h+fPnKp+PaAZKIvSEm5sbMjMz4eDggJEjR6p8vpycHACFz591tQV4QkICxo8fr7Y1H9bW1mjatKlOLwi0sbEpU1dJZc2fP1/nFwOGh4erZR5HR0f5e8v8+fPVMifhHiUReiApKQlbt24FACxbtkwtLZGLjinW5Q5/b9++RUpKitrmMzMzw8uXL3W2mVR2djb27t1bYkdYVZFKpXj27Jna5lO3Xbt2oXfv3mqbb8mSJTAyMsKDBw/0dmuyvqEkQg+sXr0aeXl5aNmyJfr376+WOYsqEZUqVdLJfgUvXryAhYUFateurdZ5u3XrBktLS53sqxEVFaWyxb5fYmhoiK5du+LMmTNqnVcdJBIJRo0apdZEvm7dupgwYQIAYMGCBTpdNSOFKInQcdHR0di1axcAYOXKlSrZMleSokqEpaUlPDw81DKnujAMA2dnZzRq1IiT+TMzM5VuSKVpZDIZZ31CKleujC5dukAmk3Eyv6ocPXpUvrhZnebPnw8LCwu8fPkS586dU/v8RL0oidBxy5Ytg0QiQefOndGlSxe1zVtUiShXrhymTJkCHx8ftc2tah8+fMD58+c5O//CxcUFJ06c0KlqxNatW9GjRw9O5ubxeIiNjcWNGzc4mV8VwsPD0a1bN1SrVk3tc1erVg1TpkwBACxdulQnK5Hk/1ESocOCg4PlXQuLTpVTl087eBoZGenMXvyi8xpGjBjBaRzjxo1TW1VJ1YRCIevHW5dV06ZNYW9vrzP9H4yMjFR2+FlpzJw5ExUrVkRYWBgOHz7MWRxE9XTjnZ2UqKg7Z9++fdG6dWu1zv1pEmFgYABra2ucOHFCrTGoglAo1IiGWGZmZvD09NSJZ8779u3TiIRIJBKx3kSNC4GBgXjz5g2qV6/OWQxWVlbybaWurq4QCoWcxUJUi5IIHfXs2TNcunQJPB4PK1asUPv8nyYRAFCnTh0MGDBA/nVtxDAMPD090aFDB65DAQDMmDEDt2/f5joMpfj6+mL48OEq6ZFRVvb29rhz545WVyMYhkHt2rXRs2dPrkPBxIkTUb16dcTGxsLLy4vrcIiKUBKhgxiGkX8KGDVqFBwcHNQew3/7ZvB4PAQGBqptz7oq5OfnY9q0aVyHIcfn81lvoKZu5cqVU/nx1mUxYsQIrW4Yl5aWhtu3b2vE40MzMzMsWbIEQOHj1KJ1UkS3cP+TRlh3+/ZtPH78GIaGhvJfYnUrqYNnmzZtkJ+fj48fP3ISk7L279+vcedeNGzYEJ6enlyHoZC7d+8iLy8PZmZmXIciZ2xsjEOHDmnlYw2pVApvb28MHTqU61DkxowZg3r16iE9PR1btmzhOhyiApRE6BiGYbBgwQIAwF9//YU6depwEseXOng2aNAAlSpV4iIkpdy+fRvjx4/XiGf3nypfvjz+/PNPpKWlcR1KmQgEArRq1QqtWrXiOpTPTJo0SSt3avB4PI37fhoZGWHlypUAgI0bNyI1NZXjiAjbKInQMWfPnkVAQADMzc0xb948zuIoqRIBFB5rvG/fPq1bG2FjY6NRn5g/FRcXh9evX3MdRpkkJCTAz89P45KyIra2tlyHUGYbN27UyLgHDx4MZ2dn5OXlYe3atVyHQ1hGSYQOkUgk8h0Z06dP52SPeJEvVSIAYOrUqUhISFB3SAo7fvw4ypUrp9F/8CpWrKg1iUR2djZevXqFrl27ch3KFzk4OMDNzY3rMEotMjKS0w8NX8Pn8+XVCE9PT8TFxXEcEWETJRE65NChQ4iIiEClSpUwY8YMTmP5UiWiSGBgoFZsT8zIyEDfvn1Rt25drkP5qnr16sHOzo7rMErFwsICbdu25TqMrzIwMMCsWbMQHx/PdSiloum/Tz169ECHDh0gEomwfPlyrsMhLKIkQkcIBAL5Vs65c+eifPnynMZTtBL7S6c6Dho0CDt37lRnSAqJjY1FZGQk12F8U8WKFXHhwgV8+PCB61C+SiqVYtOmTbCxseE6lG/KycnBixcvuA7jm44fP47OnTurpbGeong8nrwa8c8//yAiIoLjiAhbKInQEV5eXoiNjYW1tTX++usvrsOR93YoOifiv3g8HkaNGqXRR+LGxMQgKSkJzZo14zqUUhk5ciSys7M1+hNpaGioxpbd/6tChQpo3Lgx7t27x3UoXySVStGrVy+tWKzcvn17/Pzzz5DJZJztGiPsoyRCB+Tk5MiPtV6yZIlGLAD8tHfGlxgbG2t0NcLGxgYtWrTgOowySUxM1OieGtHR0VyHUCY2NjZo1aqVxiZmN2/exIcPHzR2vc5/FT3KOH36NF6+fMltMIQVlETogE2bNiEjIwP16tXDmDFjuA4HwLcrEUDhYTSjRo1CUFCQusIqtdzcXOzevVvrDnPq1q0b9uzZo5F/9Ly8vDhrsqUoMzMzPHv2DP7+/lyH8pmEhAQ0bNgQzs7OXIdSak2aNMEvv/wC4P+P5Sfajcdo4rsNKbXU1FTUrVsX+fn5OHLkCIYMGcJ1SGAYBpaWlpBKpXj37t1Xn3/n5OTg3bt3GvfIwM/PT+P23JdWXl4ejIyMYGxszHUockV9KbT1NMjg4GDY2dlpxPHcReLj45Gfn681C2qLREVFoUmTJvLDsX744QeuQyJKoEqElluzZg3y8/PRuHFjDBo0iOtwABQ2qSpa6/C1SgRQuPBSJBLB29tbHaGVCsMwyMjI4DoMhZmbm2P79u1ch1HMyZMntfrY47y8PBQUFHAdhlxSUhJu3bqldQkEANSvXx+///47AGDBggUaWTUjpUeVCC0WGxsLOzs7iEQiXLhwQWNKxSkpKfKTMnNzc7+5alwqlUIoFMLIyEgjjpXevHkzZs6cqTXPmUuSn5+PgIAAjWgW9ubNG1hYWKB27dpch6KUHTt2YMKECTAxMeE6FGRlZck75GqjhIQEODg4QCgU4uLFi+jXrx/XIREFUSVCi7m6ukIkEqFTp07o3r071+HIFa2HMDc3L9WbnIGBAe7fv68Rh9AUFBRg8uTJWp1AAICJiYlG/LEDAFNTU416tKKoCRMmQCaTcR0GRCIRTp06pbUJBADUqFEDU6ZMAQDMnz9fI76vRDGURGipN2/e4ODBgwCAlStXatQfvW+dEVGSXr16ITQ0FPn5+aoKq1QOHTqkE29oBgYGqFSpEk6ePMlpHMHBwQgKCkL16tU5jYMNJiYm2L17N9dh4NKlS5gwYQLXYSht9uzZsLS0xJs3b3Ds2DGuwyEKoiRCSy1ZsgQMw6BXr15wcXHhOpxivnVa5Zc4Ojpy+jjj+fPnGDJkiEZskWWDra0t+vXrx1mfEoZhUKNGDfTu3ZuT+VVh2rRpnDfnatq0Kafzs6VixYqYM2cOgML3M5FIxHFERBGURGghPz8/nD17FgA08gjZr/XN+BpbW1ts2bKFs0pAxYoVtXb3QEl4PB6CgoI4O3EzIyMD9+7dA5+vO28zfD4fNWrU4Gx+Nzc31K9fn7P52TZ58mRUrVoVMTEx2Lt3L9fhEAXozm+3Hlm0aBEAYPjw4Rq5R1zRSgRQeGQ3FwcSeXt7Iz09XWeqEEXatGmDjIwMfPz4Ua3zSqVS3Lp1S34mgC6pV68eJ4ekxcTEYMaMGVq9FuK/LCwssHDhQgDAqlWr5OupiPagJELLeHt74/bt2zAwMMCyZcu4DqdEilYigMJPz+ruVyAQCNC0aVONbwqlKHt7e1SpUkWtc/L5fLRr106tc6qLlZUVxo8fj9TUVLXOGxwcrNGnkSrqjz/+QL169ZCUlIRt27ZxHQ4pI0oitAjDMPKsfdy4cbC1teU4opIpU4kAgAEDBuDAgQNshvRVSUlJePXqlUYtTmVTjRo1sHv3brV+ytuwYYPGdz5VRlxcHEJDQ9U2361bt2Bvb69QYq7pjI2NsXTpUgDAxo0btfqMFn1ESYQWuXTpEp4/fw4TExN5MqGJlKlEAIWr4AcOHKiWtRE5OTnw9fVFly5dVD4Xl6ZNm6a2Dp/h4eGYO3euWubiSv369WFubo7AwECVzyWTydChQwfUqlVL5XNxZdiwYXB0dERmZiY2bNjAdTikDCiJ0BJSqRTz588HAEydOpXTxV3fomwlAig8r+Hff/9lK6QvsrCwQKdOnVQ+jyYIDw9Xy+mA4eHhKp9DEzRo0AANGzZU+TwBAQHw8fHRmHM/VMHAwEC+SHzr1q1ISEjgNiBSapREaImjR48iPDwc5cqVw6xZs7gO56sUOSfiv2xsbNCtWze8f/+erbA+I5PJ4ObmBmtra5XNoUn69++PHTt2qHSO48ePo127djq1+O9LKlasiNOnT6v0kLT8/HwYGBigW7duKptDU/Tp0wctW7aEQCCQdyUmmo+SCC0gEonkHe/mzZuHihUrchzR1xU9e1d2u6RMJpP34FCFwMBAzJs3T2Xjaxoej4fff/8dYrFYJeNLJBL07dtX7Ys4uTR69GikpaWprMIjEon0prcEj8fD2rVrAQC7d+9GVFQUxxGR0qAkQgvs2bMHcXFxqFq1KiZPnsx1ON/ERiUCAGrVqoU3b94gLCyMjbA+Ex8fr7OLKb+Ez+er7NTFe/fu4d27d3r3PU1LS4NQKGR9XKFQiMOHD6N58+asj62pvv/+e/z444+QyWRYsmQJ1+GQUqAkQsPl5eXJS3uLFy+Gubk5xxF9W1El4lsdPEujZ8+eqFWrFuufxjw9PdG1a1dWx9QG5ubmGDFiBF6/fs3quPHx8ahTp47OnKZYFl26dIGXlxfrP6MikQiTJk1idUxtsHr1agCFj8bUsXCVKIeSCA3n4eGB5ORk2NraytvnarqiSgQbSYShoSEuXLiAzMxMpccqIhQKMW7cOJ1eqPY1xsbGrJ8iaWBgoBEdWLkyceJEVo9tZhgGBw8ehKGhIWtjaosWLVpg8ODBAApbhRPNRkmEBktPT8fGjRsBAEuXLtWaTohsViIA4LfffsODBw9Y2/J57tw5vd6LbmlpidzcXDx48ICV8VJSUnD16lU0aNCAlfG0kYmJCaunWF68eFErHl2qyrJly8Dn83H9+nU8efKE63DIV1ASocE2bNiA7OxsODk5YdiwYVyHU2psViKKODg4sDJOWFgY2rdvDxsbG1bG01Zt27ZFixYtWFlkaWhoiLFjxyoflBbj8/n4888/8fjxY1bG++6773Sq50hZNWrUCKNHjwZQ2CpcXxaXaiP9/SnVcPHx8di6dSuAwiZb2vSGwnYlAih8Uy2qyijDwsICpqamLESk3fh8Ph4/foz4+HilxhGLxTh79qxW/XyqiqmpKSs/8zt27ND7JBcoXANmbGyMJ0+ecN45lXwZ/eZrqJUrV0IoFKJ169Za10qZrd0Z/zV//nzExMQofH9YWBgCAgJQvXp1FqPSXj///DNev36N/Px8hcc4d+4cxo8fz2JU2svAwADlypXD6dOnFR4jMTER48ePh5WVFYuRaadatWph4sSJAAqbDnLV3Zd8HSURGujdu3fYt28fAGDNmjVatWVOKpWioKAAgPLnRPyXRCLB8+fPFbqXYRhUrlwZ/fr1YzUmbde0aVOl1tq0adOGxWi0X4MGDdC7d295Il1WkZGRSEpKYjkq7TV37lyUK1cOr169Uio5I6pDSYQGWrp0KaRSKbp06aJ1RzIXHXkNsF+JMDIywg8//IDz58+X+d6srCw8fPhQqxIydahbty42b96s0Kc8Nzc31K5dWwVRaS8ej4egoCCF2tn7+/vDxMREpxuXlVWVKlUwY8YMAIXvi6o6KI0ojsfQihWN8urVK/nhMk+fPtW6g2bi4uLQsGFDGBoaIjs7m/U/2jKZDJmZmahQoUKpn8NLpVKcOHECI0eOZDUWXcEwDN6+fYtGjRqV+p7o6GjY2NhozY4hdbt9+zacnJxKvbaBYRgIhUKIRCKUL19exdFpl5ycHDg4OCAtLQ27d+/Gn3/+yXVI5BNUidAwRfuiBw8erHUJBPD/iyotLS1V8qmfz+cjJiYGd+7cKdM9nTt3Zj0WXVLWQ33Cw8NVckqjrmjSpAmqVq1a6uvj4uJw9uxZSiBKYGlpKW8+6OrqKn9cSjQDJREa5PHjx7h58yb4fL68o522KXoWzPZ6iE81b94cDRo0QGpqaqmu37Bhg063UVYWj8dDnz59sH///lJdf/v2bdja2rL+uEqXVK9eHZ6enqVatCqTyZCYmEiVsq+YMGECatSogcTERHh6enIdDvkEJREagmEYzJ07FwAwZswY2NnZcRyRYj6tRKiSQCAo1QmBoaGh8u8r+TJTU1MMGTLkmw3PpFIpOnbsSM/tS2HGjBl49+7dN68TiURU1fkGU1NTuLq6Aig8FjsrK4vjiEgRSiI0xPXr1/Hs2TMYGxtj0aJFXIejMHVUIgDAyckJN2/eRFpa2levU2UrcV2Tk5ODEydOfPWaoKAgPHz4UG+PDC+r6Ojorx6UxDAMPD090bFjRzVGpZ1GjhwJOzs7ZGVlwd3dnetwyP9QEqEBZDKZfC3EpEmTtLr0XrQ7Qx2l7pEjR3714Kjjx4+jVatWMDAwUHksuqBWrVro3LnzFxOv/Px8SCQS9OjRQ82Raa8+ffpg27ZtX3w9Pz8fU6dOVWNE2svQ0BArV64EAGzatAnJyckcR0QASiI0wsmTJxEUFARLS0vMmTOH63CUUpREqLoSARQ2kjp48CAkEslnr4lEIvTv3x/VqlVTeRz6QiKR6HWTLUXweDyMHz/+i4/e9u/fT9/TMhgwYABatGiBgoICebdPwi1KIjgmFouxbNkyAIXPUCtXrsxxRMpRZyUCACZPnoxr16599vVHjx4hPDyczoUoo1q1auH169cIDw8v9nWRSIT9+/frZatvNpS0aPXu3bv4448/6Ge0DHg8nrwasWvXLqVOsCXsoCSCYwcOHEBkZCSqVKmiE2XNoiSCzb4Z32Jra1vs3+Pj42Ftba2VW2Q1Qd++fWFtbV3sWX5BQQGmTZvGYVTay8LCAr/88gtevXpV7OvW1tYwMzPjJigt1rVrV/zwww/FPoAR7lASwaGCggKsWLECQOH5ELqwZY6LJMLe3h5ubm7yfzcyMqI3ZyUYGBjg8uXLyM7Oln/t33//pbUlSjAxMSm2GPXkyZMwMzOjxmUKKtoCf+TIEYSEhHAbjJ6jn2AObd++HQkJCbCxsdGZJkZcJBFGRkaYNWsW4uPjkZaWhgsXLqB+/fpqm18XjRo1Cnfu3AHDMDh37py8ERJRjKWlJdLT0/Ho0SNkZmaiV69en1XQSOm1bdsWffv2BcMwWLp0Kdfh6DVKIjiSlZWFtWvXAig8hU1XtsxxkUQAhd9PPz8/MAyDP/74Q61z66omTZoAABo3bkxVCBa0a9cOTZo0QXR0NN6+fUtrIZS0fPly8Hg8nD9/Hr6+vlyHo7coieCIm5sbsrKy0KBBA506qY6rJKJSpUpo1KgR3N3dqUTMkoYNG2LQoEG0w4UlfD4fly5dQlBQEFq0aMF1OFrP0dERI0aMAPD/7QKI+tG7LQeSkpKwefNmAED37t0RExODTZs2QSgUYv369QCA9evXIz4+HgcPHoSfnx9u376NCxcuICwsDJ6ensjOzi52bXp6Onbt2oWgoCBcuXIF169fR0BAAPbt24eUlJRi1+bn58PDwwMRERE4c+YMvL294ePjgyNHjiA2NrbYtVKpFG5uboiJicGxY8fw5MkTPHz4EKdOnUJkZCQ2b94MgUAgv6fo+aSfnx/8/Pxw69YtXLx4ESEhIdi5cyeysrKKjZ+ZmQkvLy8EBwfj0qVLuHHjBvz9/bF//34kJSUVu7agoAAeHh6IjIzE6dOncf/+fTx9+hRHjx7Fhw8f4OrqioULF2LdunWQyWTYsGEDPnz4gKNHj+Lp06e4f/8+Tp8+jcjISHh4eKCgoKDY+ElJSdi/fz/8/f1x48YNXLp0CcHBwfDy8kJmZmaxa7OysrBz506EhITg4sWLuHXrFvz8/HDgwAEkJiYWu1YgEGDz5s2IjIzEqVOn8PDhQzx58gTHjh1DTEwM3NzcIJVKi90TGxuLI0eOwMfHB97e3jhz5gwiIiLg4eGB/Pz8YtempKRg3759CAgIwPXr13HlyhUEBQVh165dSE9PL3ZtdnY2PD09ERYWhgsXLuD27dvw8/PDwYMHER8fX+zamJgYtG7dGoGBgTh58iQePnyIx48f4/jx4/K4JRJJsXvi4uJw+PBhPH/+HHfv3sXZs2cRHh6Obdu2ITc3t9i1qamp2Lt3L16/fo1r167h6tWreP36Nfbs2YO0tLRi1+bk5GD79u0ICwvD+fPncefOHfj6+uLQoUOfxS0Wi+Hu7o7379/j+PHjePz4sTzu9+/fw93dHWKx+LPftUOHDsHX1xd37tzB+fPnERYWhu3btyMnJ6fYtWlpadizZw9ev36Nq1ev4tq1a3j9+jX27t2L1NTUYtfm5uZi27ZtCA8Px9mzZ1G9enW8e/cO+/btQ1xcXLFrJRKJ/HetKO6HDx/i5MmTiIqK0pn3iPXr1yMxMREHDhxQ+j1i0qRJAABvb+8y9dMh7KEunhyYMmUKPD090bJlSwwdOhRTp07VmdJmhw4d4O/vj7Nnz6JXr15qnTsmJgahoaGoUaMGmjVrpta5ddHTp09RrVo1vH79GoMHD+Y6HK2Xl5eH/fv3Y+DAgbCxsaFHRCzYtWsXQkNDsXfvXrRo0QJ+fn46816qLagSoWbR0dHYvXs3AGDlypX4448/IBaLOY6KPerqnfFfbm5uqFmzJn7++WcYGBhQLwIl+fv7g8/nw87ODh06dMC5c+e4DknrvXnzBtOmTUPt2rXh7u7+1eOwybeJRCL89ttvWLJkCUxNTeHv74/z589zHZbeoSRCzZYtWwaJRIIff/wRXbp0gZmZGXbs2MF1WKxRV++MT0VFRWHq1KkwNDQEUJjIULtgxTEMAycnJzRu3BgAUK1aNXTt2hUymYzjyLQXwzDF+rzMmzfvswO9SNmcOXMGmZmZqFatGqZPnw4AmD9//jebyBF2URKhRsHBwThy5AiA/9/nzOfz8ddff8HHx4fDyNjDRSUiKiqqWDXHxcUFR44c0akKjzrFx8fj9OnT8sWxfD4fb9++hbe3N8eRaS8PDw/89NNPxb4WFhbGUTTaLyIiAh06dEDNmjUBADNnzkSFChUQGRmJw4cPcxydfqEkQo2WLFkCAOjfvz/atGkj/7qJiYlOnJ/PMIzaKxG3bt1CrVq1Pkta/vzzzxJ7apCvk8lk+PDhA0aNGlXs661atUKtWrWQmprKUWTaq6CgAJMmTSr2rJ7H46Fnz57Yt28fh5Fpr/8e3mVlZYV58+YBKKz20uNM9aEkQk2ePXuGixcvgsfjwdXVtdhrhoaGqFChAs6cOcNRdOwQCATykrc6KhFSqRTff/99iQdLmZiYyNeekNITi8VffFYvk8koMVPA4cOHS/y+mZmZYdiwYfQ9LaOwsDC8evUK1tbWxb4+ceJE2NjYIC4uDrt27eIoOv1DSYQaMAyDhQsXAgBGjBgBBweHz65p0KABevfuLX8coI2KqhCAeioRISEhuH//PoyNjUt8ferUqbh586bK49AVDMNg27ZtaN++fYmvOzg44MqVK0hPT1dzZNrrxYsXGDhw4Bd/HzIyMrT+w4M6MQyDatWqoW/fvp+9ZmZmJn+fXb16dbH3I6I6lESowe3bt3H//n0YGRl9sWEMj8fDq1ev8O7dOzVHx56iBMjCwkLlBz7l5+ejoKAAPXv2/OI1BgYGqF69ukrj0CW5ubmYOXPmV68ZM2aMTjx6U5cKFSp8NaGuU6cO2rVrh6ioKDVGpb2ysrJw586dL76/jBkzBvXr10dqaiq2bNmi5uj0EyURKvZpFWLChAmoW7fuF69t164dUlJSkJiYqK7wWFWU+avjtEqpVFqqJlt16tSh0mYpHTp06JvJn5GREQ4fPkwr4Evh4cOHSEpK+mZVzsDAgM6MKAWZTIbr169j2LBhX7zGyMhI/rjY3d292I4YohqURKjY2bNn4e/vD1NTU8yfP/+b1zs4OKBixYpqiIx96jryWiwWY+/evfLeDl9TqVIl/P777/Rm8g23b9/G2LFjS1VBmjx5Mi5duqSGqLSXUCiEk5MTOnTo8M1ra9WqBT8/P7x9+1YNkWkvHo/3xUdtnxoyZAicnJyQk5ODdevWqSEy/UZJhApJJBL5me4zZswoVQ8CGxsbeHl5IT8/X9XhsU5dSURubi5mzJhR6utjY2NpO9031KlTB+bm5qW+3t7eXoXRaL+kpCQEBASU+vTEgQMHokqVKnQA1Vds2LDhq5XcInw+H6tXrwZQ2Ck5Li5O1aHpNUoiVOjw4cN49+4drKysyvRHb8aMGYiJiVFdYCqiriTixIkTZVpzYWdnB2NjYwQHB6swKu115swZ8Pn8Mn9P3d3dVRiV9srNzcWzZ8/QrVu3Ut/D5/Nx7do1Wgz4BREREZg7d26pr+/Rowfatm0LkUiEFStWqDAyQkmEiggEAvkiygULFsDKyqpM979580brPpWoI4k4ffo0/vzzzzLfZ2dnhwYNGqggIu2WmZmJbt26oWHDhmW6z8jICNOnT8fHjx9VFJn2Mjc3R6dOncp838iRI3Hjxg2t+71Xhzdv3pTpeh6Ph7Vr1wIADhw4gIiICFWERUBJhMp4eXnh48ePqFGjBv76668y3z9gwAB4enqqIDLVUUcS0apVK4UWoVWqVAknT55EQkKCCqLSXgkJCQrvDMjNzUVgYCDLEWk3mUwGNzc31KhRQ6H7W7duzXJE2u/kyZNo27ZtmX/v27dvj549e0Imk2Hx4sUqio5QEqECOTk5WLNmDQBg8eLFpdpF8F88Hg9jxozRqoNoVJ1E7Nq1CxUqVFD4/rFjxyI5OZk+6f1PXFwcoqKi0KJFC4Xur1ixIurXr48HDx6wHJn2CgoKkp+cqIh69ephw4YNLEak3SQSCXr37l2q9WQlKXqUcebMGfj7+7MZGvkfSiJUYPPmzUhLS4OdnR3GjBmj8Dg8Hg979uxhMTLVUmUS8fHjR4wcOVLpnSspKSnUU+N/rK2t0a5dO6XGqFOnDlq1akXNuf7n48ePSrWi5vF4mDt3rlafF8OmBw8eICIiQuHvaZMmTeRbQkuzO46UHSURLEtNTZUvOHN1dZV3llREuXLlMHToUAQFBbEVnkqp8pyIjx8/snJSYrdu3bBz504WItJuBQUF8PT0RKVKlZQax8zMDPfv30dISAhLkWmvXbt2oUuXLkqPI5PJaBEwgMTERNSoUUPhSlmRZcuWwcDAAHfu3KGqmQpQEsGytWvXIjc3F02aNMGgQYOUHs/IyEjlpz+ypegsBkUe33xNQEAAxGJxqbZ3lcZff/2l9w16goKC5O2TldW7d29IpVK9/p4KhUKMGTMGpqamSo9lZGSE1q1b4+zZsyxEpr34fP4Xj7Qvi/r162PcuHEAChe50+NMdmnHXyctERsbix07dgAAVq5cycoffysrK2RmZuLx48dKj6VKt27dkvcAuHz5MmvjMgwDR0dHNG3alLUxTUxM4OXlxdp42ubTbqtsEYlEep1EXLhwgdUDzaytrdG9e3etWhPFprS0NFy6dAl2dnasjLdw4UKYmpri2bNnuHLlCitjkkKURLBoxYoVEIvFcHFxQffu3Vkb18XFBU2bNtXIN5Tk5GSMHTsW/fv3l2f44eHhKCgoYGX8xMREnDx5ktVHJHw+H+PHj9f4xExVNm/ezErZ/VNt2rTBgQMH9HK9SXh4ONq2bYtatWqxNiafz0dYWBgePXrE2pjahM/n4/fff2dtvBo1amDSpEkAgEWLFtEaHhZREsGSN2/e4MCBAwAKH2kos7jqvwwMDPDgwQPEx8ezNqayGIbBoUOH0KxZM5w8eRJ8Ph8TJ05EzZo1kZuby0obbplMhnfv3mH06NEsRFycqampWnp8aJq8vDxMnTqV1Z/PIpMmTYJIJGJ9XE1nYWHBymOM/2rTpg2qV6+O1NRU1sfWZBKJBGfOnGG9n8icOXNQvnx5BAcH49ixY6yOrc8oiWDJ4sWLwTAMevbsqfSK95L06dMHAQEBEAgErI9dVhEREejZsycmTpyIjIwMNG3aFA8fPsSWLVvkB2y5ubkhKytLqXkkEolSC1O/xtDQEGZmZjh37pxKxtdUJ06cUNkfemNjY+zbt08lY2uqiIgI+Pr6wtraWiXj8/l8vfvUfPbsWUyYMIH1cStWrIjZs2cDKFxsqY8JrypQEsGCly9fyv8YrVy5UmXzNGvWjNM2zEKhEGvXrkXr1q3x8OFDmJubY+3atXj8+DFatmwJAPj1119hb2+P9PR0eHh4KDXf1q1b4eLiwkLkJWvUqBF69Ogh35qq654/f47evXurtAIzefJkXLt2TWXjaxKGYVChQgUMHDhQZXPY29vjwoULyMzMVNkcmkaVB25NnjwZlStXxvv37/Uu4VUVSiJYULT/ePjw4XB2dlbZPHXr1oWbmxsnq4ufPHkCFxcXrFq1CiKRCN27d8fLly8xc+bMYtUCQ0NDLF++HEBh85ukpCSF5svKysKsWbPYCP2LeDweAgMDER0drdJ5NEXVqlVhaWmp0jkMDQ1Ru3Ztlc6hKbKzs3H//n2VPBr61B9//KHyOTSFm5ubSn9+LCws5KdXurq6Ii8vT2Vz6QtKIpTk7e2Nu3fvwsDAAEuXLlX5fPPnz0dkZKTK5ymSmZmJKVOmoFu3bnjz5g2qVauGQ4cO4cKFC6hXr16J9/Tr1w+tWrVCXl6ewqfvHTt2TC1vnO3atcPHjx+RmJio8rm49PjxY3z8+BEWFhYqn8vGxkarDklThEwmw5UrVzB06FCVz2VgYIDjx49DKpWqfC4uxcTEYNq0aSqvtv7xxx+oVasWUlNTsW3bNpXOpQ8oiVACwzDyVt9//PEH6tevr5Z5X79+rfI5GIbBmTNn0Lx5c+zfvx9A4bHRAQEBGDp06Ff/wPN4PKxatQoAsG/fvjJ/0r958yZGjhyptvMxmjZtisqVK6tlLi4IhUI0atQIHTt2VMt8lStXxujRo3V6QSCPx0Pnzp3VNt/EiRNx/vx5tc3HhfDwcLWs+TI2NpYfh71u3TpkZGSofE5dRkmEEi5dugRfX1+YmprKkwlV4/F46NWrF/755x+VzRETE4PBgwfjt99+Q2JiIho1aoTbt2/Dy8ur1Ccc/vDDD+jatSvEYrE8oSitBg0aqHXnhLW1NXbs2KERi1ZVISUlBYGBgWoticfGxup058QNGzagZs2aap2zSZMmap1Pne7evYs6deqUuduxooYNGwZHR0fk5ORg48aNaplTV/EYOr5LIVKpFE2bNkVISAjmzJlT5j+UysrKykK5cuVY3QYlkUiwc+dOrFy5Enl5eTAyMsK8efMwd+5cmJiYlHm8ly9fomPHjuDxePD19UXjxo2/ec+5c+fg5OSE7777TpH/C0oJDg4uVYzaJC8vD5cuXcKIESPUPrePjw/Kly8PJycntc+tSmFhYWjYsKHKdg59iUgkwvbt2+U7DHSFVCqFSCQCn89X6H1GUZcvX8bQoUNhZmaGd+/eKdx5Vd9RJUJBR48eRUhICKysrFS+ALAkmZmZOHXqFGvjBQQE4Pvvv8f8+fORl5eHDh06wNfXF0uWLFH4F7tly5YYNGgQGIaRlw+/JiMjAz/88AMnCQQAvH//XueOxDU3N0fXrl05mdvBwYG1Ewc1yfv37zmZ19jYGFOnTkVsbCwn86tKSEgI7t+/r9YEAijcNt+2bVsUFBSodFedrqMkQgEikUi+A2H27NlKd5ZURN26ddGhQwfExMQoNU5eXh4WLFiAjh07IiAgABUqVICnpydu3boFe3t7peN0dXWFgYEBrly5Ah8fn69em5ycrPT/H2X07t1bfmy5LmAYBm5ubgq3UVZWhQoVcPToUYV36GiikydPonnz5mqvQhTJzc1FWFgYJ3OrQkFBAUQiEX7++We1z83j8eTJw969exEVFaX2GHQBJREK2LNnD96/f4+qVati8uTJnMaizCfnGzduoEWLFti6dStkMhmGDBmCgIAAjBs3jrVFjY0aNZKfOLl06dIvxvvx40e8efMGzZs3Z2VeRfD5fIwbN05nDqHx9/fH3LlzOY1h3LhxiIuL04kKj0gkQt++fVV2sFRpVKpUCbVq1dKZbpRisZj1kynL4vvvv0fXrl0hlUrlB+WRsqEkoozy8vLkpflFixbB3Nycs1jq1KmDly9flnnLZ2JiIn777TcMHDgQHz58QJ06dXD+/HkcOXJEJW+QixYtgomJCZ48eYJbt26VeI21tbVaV7t/iUwmU+miVXVKS0vTiPMFsrOzNbLvS1k9ffoUoaGhnH9PbW1t0aZNG60/yVIsFmPfvn2cfnAA/v+AwKNHjyIwMJDTWLQRJRFl5OHhgdTUVNSqVUveXpZLAwYMKHW5uugPZPPmzXHmzBnw+XxMnz4d/v7+6Nmzp8pirFWrFv7++28AhY83/vvmJxAIsGPHDlSoUEFlMZSWpaUlBg4cqPVvJl5eXmrb0vktP/74I7Zv3851GEpJSkpC5cqV0apVK65DgZmZGe7cuYM3b95wHYpSCgoKWGtHr4wWLVpgwIABACA/iIqUHu3OKIP09HTUrVsXubm5OHDgACcr3kty+PBhDBw48KunEb558wZTpkzBkydPAADNmzeHp6en2j4FpKWlwdHREdnZ2Th06FCxQ3qePn2K9u3bqyWO0sjKykJCQgIra0K4UFBQoPaV7t8iFAohk8lgZmbGdSgKSU1NRWZmpkYtFPXz80OTJk1gbGzMdSgK8fT05PxxcJHw8HC0aNECMpkMjx8/RocOHbgOSWtQJaIM1q9fj9zcXDg6OmLYsGFchyM3evRo3Lp1q8TnzkKhEKtXr0abNm3w5MkTmJubY8OGDXj48KFay4iVK1fGzJkzAfx/y3SgcE2Hpq1BsLKyQnJysjzh0jZXr17VuBM4DQ0NtfYUy4yMDJw9e1ajEgig8HdHKBRyHYZCzp8/j4kTJ3Idhtx3332H3377DQCwcOFCnVjDoy5UiSilhIQE1K9fHwKBAKdPn0afPn24DqmYt2/fws7Ortjz2sePH2Py5MnyQ39+/vlnbNmyBXXr1uUkxtzcXDg5OSE5ORnbtm3DhAkT4O7ujlmzZqntdMrSkslkyM7ORrly5Thbia+I8PBwGBoaokGDBlyH8pns7Gy8fv0anTp14jqUMklLS0PFihU17mcUADZv3oxp06Zp1c8oUPh+1bBhQ67DKCYuLg6Ojo4Qi8W4fv26Sh/x6hLN+63QUCtWrIBAIEDbtm3Ru3dvrsP5jJ2dnbxPRUZGBiZNmoSffvoJERERqF69Oo4cOYKzZ89ylkAAQLly5eQne65duxbJycmYPn26Rr458/l8PHnyROM+0X+LpaWlWk/7LAtzc3ONWPdSFlKpFOfPn9fIn1EAmDZtmtadtLp7926NPGb+07Vbc+fO1fqFq+pClYhSuH79Ovr06QOZTIYNGzagSZMmxcpdZflnVV0LFK523rx5M96+fYvs7GwAhT09Vq1axclZFiURiURo0qQJYmJiMGDAAOzZs0flnSWVcfbsWfTu3RumpqZch/JNkZGRCAgIwC+//MJ1KF8UGhqK8PBwlbbPZtOxY8fw66+/ch3GV23fvh1Tp07lOoxSSUxMhKWlpVoawSkiNTUV3333HfLz83HixAmNemytqSiJKIW6deviw4cPXIdRJt999x08PT01coHQ8ePHMW7cOJiYmMDV1RVTp07V2HLshw8fUKtWLY39JFqEYRgkJibC2tqa8y2IX8MwDHJzcwFAo5PHIjExMZxW70pDLBbj5s2bGveItSTPnj2DtbX1FzsAa4I1a9Zg9erVaNCgAd68eaOx702agpKIb3j9+jWaNWsGoHB/tpmZmfxN+tM367L885de/5Qic7x+/RoCgQBWVlaIiYnRqNX5n5JKpejdu7f8wJymTZvC09MTLVu25Diykm3YsAHz5s3T6D/OOTk5uHnzJoYMGcJ1KN/k4+ODChUqwMHBgetQvsrd3R3Tpk3Tit0PYWFhGv/9fP36NfLy8jRqJ1ZJcnJy4OjoiNTUVOzZswcTJkzgOiSNRknEN/Tp0wdXr17FoEGDcPToUa7D+aKjR49i/PjxqFChAp4/f67Rn56fPHkCsViMmJgYLFy4EBkZGeDz+Zg0aRJcXV017pk+wzCIiIjgrKfHt8hkMvz777/yk0G1wbVr19CyZUtUr16d61BK9P79e9SoUUMrHmMBhZ1aL168iPHjx3MdyheJRCIIhUKtqEBt374d8+bNg42NDSIjI7V2a7I6aOZfGQ3x+PFjXL16FXw+H66urlyH80VCoVB+6tqcOXOQkJCgscfiCoVC1K9fHz/88APGjBmDV69eYdiwYZDJZNixYweaN2+Oq1evch3mZ0JDQ7kO4Yt4PB66devGdRhl0rJlS41cXFckKipKq7ZPVq1aFSNGjEBycjLXoZQoISEBx48f14oEAgAmTJiAmjVrIj4+Hp6enlyHo9EoifgChmGwcOFCAIXnMDRq1IjjiL5s//79+PDhA6ytrfH333+jbdu2qFGjBtLT07kO7TPp6enFGghVq1YNBw8exKVLl1CvXj3ExcVhyJAh+PXXX5GQkMBhpP+Px+Ohe/fuOHDgANehlGjDhg1a18a4evXq2LZtm0b+ob579y6sra1hZWXFdShlEh8fj+joaK7D+AzDMIiOjsaYMWO4DqXUTE1NsXTpUgCF5wMVLVQnn6Mk4guuX7+Ox48fw8jISKOPQs3NzcX69esBFB7ZWtTLQyqVyg900hT5+fm4e/cuunTp8tlrP/30E16+fIlZs2bBwMAA58+fR/PmzbF3716N2GplYWGBX375ReN6QISEhGD27NkavV7jS2bNmqVxRzdLpVJ06NBB4w6WKo2GDRtCJBJpXNVMLBZDKpVyHUaZjRw5Eg0bNkRaWhrc3d25DkdjURJRAplMhnnz5gEAJk2ahFq1anEc0Zft2LEDKSkpaNCgQbFM38nJCZcuXUJmZiZ3wf2HqanpVw9wMTc3x5o1a/DkyRO0bNkSWVlZmDZtGrp27aoRb4xpaWk4f/4812EUExcXp5UJRJGPHz9q1OmAb968wd27dzV2UfK3ODs7a1wCtG3bNo3p41IWhoaG8maLmzZt0thHRVyjJKIEJ0+eREhICMzNzTFnzhyuw/mi1NRUbNmyBQCwbNkyGBkZFXt97NixGrM9iWEYuLm5oUqVKt+8tmnTpnjw4AHc3d1Rrlw5PHv2DC4uLvIDv7hSr149tG7dWmNKxidPnkTjxo015r+xInr27Ilt27ZxHQaAwp4jOTk5GnmYXGlZWVnh8OHDGvMHLzc3VyOabClqwIABaNKkCfLz87FmzRquw9FIlET8h1gsxqJFiwAAs2fPLtUfPa64u7sjOzsbTZs2LXFrn5GREQ4dOqQRpcQXL17IqzulYWBggMmTJ8Pf3x99+vSBWCzG+vXr0bp1a04XjRoYGGjErheRSIR+/frBxsaG61CUwufz8eeff2rE2gipVCp/HKjNxo8fj+joaI2o8Bw6dEgjfl8UxePxsHbtWgDAzp07ERMTw3FEmkd7/+uqyIEDBxAdHY3KlStr9ClwcXFx2LVrFwBg5cqVX/xFnTx5Mi5fvqzO0EqUk5OjUNm9du3aOHXqFI4fPw5ra2tERkaiZ8+e+PPPP5GWlqaCSL8dz7Nnz/D+/Xu1z/0pX19fBAUFafWjjCJisRhHjhzhNAaJRILdu3ejSZMmnMbBloKCAs7X79y9exe//fYbDAwMOI1DWV27dkWnTp0gkUg0epceVyiJ+ERBQYH8GdiCBQs0ejvS2rVrIRQK0bFjR/z0009fvZbrnSVeXl5wcXFR+H4ej4cBAwbg1atX+Ouvv8Dj8XDkyBE0a9YMx48fV/snrsGDB3PaAyIpKQkWFhZo06YNZzGwqXz58ujTpw9evXrFWQw5OTmYMWMGZ/OzrXPnzti6dSunMdSqVUtjj7cuq1WrVgEADh8+rBHrszQJJRGf2L59+/+1d+dhUZZ7A8e/AwMqoGaJS6ng9mqkpOWbmuba0Y6ldbK0NC3XytJSczspLpnnuKSWmbmgqWAumZbpq5UG6ZUocFzAHRfCMRQ9sqToMDPP+8cw4wAzwAyzBb/Pdc1VwPPc9z2Pyv2b373xxx9/UL9+fa/epezs2bOsXbsWMGYhSvo02qRJEz755BN3NK2I27dvM2zYMKf8MqlevTqLFy9m79695h3lhg4dSu/evblw4YITWls6vr6+/PDDD+btm93N39//L7f8sCRVqlTxaIezcePGv/wn5sLGjBnDrVu3PFL3tm3bUBSl3DzTtm3b0rt3bxRFMQ93CyMJIvJlZWWZl0pOnTrVq2dnz5w5E4PBwHPPPUf79u1LvN7f358xY8Zw5coVN7SuoJ9++snp9bZv356DBw8yc+ZMKlWqxN69e3n88cdZsGCB25a1Dho0iJ07d7o9C5KZmcnGjRtp1KiRW+t1terVq6PRaIiLi3N73Vu3bvXqDw2O8vHxYc2aNW6vNysri65du9K8eXO31+1K06dPR6VS8d1333H48GFPN8dryLbX+aZOncrHH39Ms2bNSExM9NoIOjExkY4dO6JSqYiPj+eRRx4p1X3Xr18nISGh2CWWznb69Gl8fHxcOpySkpLC6NGjiYmJAYxL3D7//HO3pPovXLhAw4YN3Tov4dq1awQHB5eLuRCFKYrCzZs3qVatmltXnFy4cKHcBWUmWVlZHDt2jE6dOrmtzjNnzpCTk0ObNm3cVqe7jBgxgqioKLp168bevXs93RyvIJkIjGPMpqWS06dP99oAAjBP7Hn11VdLHUAA1KxZkwYNGrB//35XNa2IGjVquDzt3qRJE3bt2sWqVat44IEHSEpKokuXLowbN87lu8w1bNjQnL1yB71ezw8//FCmACI9PZ3hw4eTnp7uxJY5h0ql4uDBg25dnrhixQqqVavmtvrcLTAwkODgYLfVd+XKFc6ePVsuAwgwfthUq9Xs27ePn3/+2dPN8QoSRACzZ8/m9u3btGrVihdeeMHTzbEpJiaGvXv34ufnZ96S1R6hoaE89thjbknBX7hwgZiYGLccsKRSqRg4cCBHjhxhwIABKIrCsmXLaN26Nd9//71L6504caLb5mNs2LCBoUOHlqmMuXPnEh0dzbx585zUKud69tlnOXDggFv2A/njjz949dVXvXoZd1mp1WoMBoNL/x1Yql27ttef0lkWISEh5qGvf/7zn16xjNbTKnwQcenSJZYtWwYYz5H31jSxoijmwGHYsGGEhobaXUZAQAB79+4tcHaFKyiKgr+/P/369XNpPYUFBwcTGRnJzp07adSoEVeuXKF///70798fjUbjkjr1ej1JSUkuKduSoih06dKlTGVoNBpWr1pFY2D1qlUueyZl1aFDB7fMSdJoNFy/ft3l9XhaWFgYnTt3Jicnx6X13LlzhyVLlnj1wWrOMGnSJAICAoiPj/e6HWw9ocIHEdOnT0ev19OpUyerZzp4ix07dpCQkEBgYCCTJ092uJw+ffqg1WpdOgHx1q1bHDp0yGMBWbdu3UhISGDChAmo1Wq+//57WrduzZdffun0jbf8/f1p06YN27dvd2q5hS1cuJBatWqVqYwFCxYQqCjEAAEGg8dW7JTkoYceYu7cuS79lHf8+HFyc3Np2LChy+rwFiqVilOnTrn8QLukpKRytUzWltq1a5v3EJo2bZpXbObnSRV6YmVycjLh4eEoikJsbKzXrrvX6/W0adOG06dPM3ny5DJveBIXF8cjjzzikn0wFEVhzZo1ZU67O0tycjKjRo0iPj4egCeeeIKlS5fSokULp9Wh1+u5ffs2AQEBLplPc/78eerWrVum3RQ1Gg1hzZsTodPxITAbmO3nx4lTp3jooYec1lZnURSFkydP2jXvxx5arRatVktQUJBLyvdG33//PW3btnXJEKOiKOzdu/cvdyS9o7KysmjevDmZmZmsWbOGN954w9NN8pgKHUT06dOHHTt20KdPHzZt2uTp5ti0fv16Ro4cyf3338/JkyedMllx8eLFjB492umdnqIoZGRklPlTszPp9XpWrlxJREQEOTk5qNVqunfvzsCBA6lTp45T6jhx4gQ6nY5HH33UKeVZOnr0KE2aNClTh7d06VJid+wg1WCgGpAFhPr40Ll3b955550yt/HatWucOHGCzp07O2WbY0VROHjwoEvG12/evElcXBx///vfnV62M8XExPDwww87rdPPzMykatWqLgl0N2/ezEsvveTVW1xnZmYSHx9Pt27dnPIMNm/ezKpVq6hZsyZpaWlUrlzZCa3866mwQURcXBzt27dHpVKRkJBAWFiYp5tk1d27d2nZsiVpaWnMmTOHsWPHOqVcrVaLTqdz+lkBc+fOZeLEiV45t0Sj0TBu3Di3TTLzJj7ALMDyUPvZwHTA8wetC/HXtmjRogoxlGNNhQ0iOnfuzK+//sprr73GypUrPd0cmz7//HMmTJjAgw8+SHJyMlWqVHFa2Z9++qlTT9g7fvw4zZs3x9/f32llukKnTp2Ij4+natWqTk3l3717F19fX6ftcWD5T7OsyzoNmZmkAZaLGbOABoBvjRpl+rR79+5d81kiISEhTv07ajAYUKlUqFQqdDod165do1atWg4/Y0VRuHPnjlPb6ApXrlwhOzubgIAAGjRo4NSy79y549RPzZmZmVSrVs2rsxB5eXmcP38eMJ5/46zdUTMzM0lPT+e+++4jLS2tQg2Pmfx1zxAug59++olff/0VPz8/pk6d6unm2JSTk8PcuXMB+PDDD53+i2/UqFHs3r3baRtQXbt2zWszOiZxcXHEx8fj4+NDbGwsDz/8sFPL//LLLxkwYIBT9h44ffo0KSkpPPfcc4AxGJg6dSqzZ88u9TCMeS4EBQMIgOrABOCjnBwe7tKFhQsXWi23pHpfeuklLl68yPPPP8/GjRvte5MlOH/+PMeOHePFF19k7NixfPnllwwcOJCFCxc6VN78+fOZMGGCU9vobElJSbRt2xaAPXv2OH3PhT179tCzZ0+nlBUfH0+9evWoW7euU8pzlSFDhnD+/Hm6d+/ODz/84LRy8/LyePTRR7l48SKLFi1yaOn9X533ho4uoigKU6ZMAYy7j4WEhHi4RbYtWbKE69ev07RpUwYPHuz08tVqtdM+iW/evJlmzZq5dadBeymKwocfGhP6gwcPdnoAATht++Q7d+5w/fp1cwABju3xYFqR0Q94HSi8xdRoQKXXs23bNpvlFlfvgQMH2LlzJ76+vsyaNavU7Sqtxo0b06pVK+Li4sq0PDU9PZ3XBw9m4MCBTm+js0VERKAoCn379nXJpk1PP/00ixcvdkpZwcHBXn1QIcCxY8fMwa3pIC1n8fPzY8aMGQDMmzfPIycLe1qFCyK2bt1KYmIiAQEBTJo0ydPNsSkjI8P8Dz0iIsIlnbNKpaJ27dpERkaWqZy7d+/Sp08f6tWr56SWucauXbv47bffqFy5sssyUL6+vkRHR2MwlG2mgcFgKDCB9ujRo6xYvpwQSt+JmvaFGK/X8xmwDphT6Jo/Ab2iGDvnlSuLlFvc3hKKopif45AhQ1y2vbmfnx+rV68u0/LUuXPnsnnLFubPn++SNjrL/v372b17N2q12tw5OZuvry9vv/02ubm5ZSrnt99+IzU11etT+BEREQD069eP1q1bO738l156ifDwcP7880/mzCn8L6z8q1BBhE6nM6ebxowZ41UrCAqbP38+OTk5tG7dmhdffNFl9dSqVYsBAwaUadOdhIQEjh496pWTKU0s/+zfffddly5rfOutt/j2228dvl+v17N06VJatmxp/t6IESMwKArdKH0naspC9AVWAI3z/2sZJvwbqArEAJX1ej766CPzttjp6en06tXLZue9Y8cODh06REBAgEtPNvTx8WFjdDTj9XrqAeP0eruyERqNhsiVK2kMrImM9NpNtiwzZcOGDaNJkyYuqys3N7dMQ09arZbGjRu79UwOR8TExPDjjz+iVqvLvDTeFh8fH2bOnAkYV0F5698vV6lQQcS6des4ffo0NWrU8OqZtL///jsrVqwAjEd9u3rCUmpqKikpKQ7dm56eTqVKlWjXrp2TW+Vc0dHRnDp1ivvvv5/x48e7tC6VSkV4eLjD92dmZhZoY2JiIieTk2kMbACGl6ITtcxCLAUCMQYKAcBc0zUYg4rxQD3gA0Vh/bp15qGLadOmkXL2rNXOW6fTmT/hjR49usxj4sWd6bFgwQKCVCpG5389GvuyEQsWLCAICgRClvV5y3ki3333HfHx8WXeUK407rvvPnr27MmRI0ccuv/69escO3bMqz84WGbKhg8f7tJD1nr27EmHDh24e/euOaCoKCpMEHHnzh1zJDphwgSXHwxVFnPmzOHu3bt06tSJ7t27u7y+5s2bo9frOX36tN33VqlSxeu3ub19+7Z5vH7ixIncd999Lq8zJCTE4cl/W7duLRA4vvXWW/c6QYyrKkrqRBctWoRWpyMOWMa9QGFc/tcDgB5AFTB3zn0BLIY2NkRHE2Txc8vOe926dZw5c4YHHnjAKcuObc27sAyGTBNDq3MvoDl69GixAYApC1E4EIqIiDDXN336dKKjo102fFAalkHZe++957T9S4oTGBjo0L+F3NxcYmNj6dGjh/Mb5UTbtm0jMTGRoKAg8zw4V1GpVObfMZGRkZw7d86l9XmTChNELFu2jMuXL1OrVi3eeustTzfHpjNnzrB+/XrAmIVwV6QfFhZm9xbA2dnZREVFef3WwV988QVXrlyhfv36vPnmm26ps1KlSrzzzjtcvnzZrvu2bNnCkCFDzF+bshATuRcErKXkbESDBg1o9cgj7MKYhSgQCACbgNMYV2eYOucv8v8/BvDX6UBRmGjxc8vO29ThTpo0qcwBeXHzLkxDMqML3WMKaEaMGGEOAKwFE3PmzCkSCFXW6/k6KsocLEWvW0djIHr9+iLP05SlSEpKcmm2Yu3atZw7d47g4GCnLrsuTvXq1blw4QKHDx+2677KlSt79REBYFw1YfrQ+N5777ll6PrJJ5+kR48eGAwGr17152wVIojIyckxR4kRERFevUZ8xowZGAwGevfubV7m5Q41atRg3bp1ZGRklPqeW7duMWrUKBe2quxu3LjBggULAOOzdeeuctnZ2XZnd5544gn8/PzMX5uyEIWDgJKyEWPGjKFF/iSywoGAaYGjZXBhGtr4AGOwEgYF6jUZDVTW6cjIyCAkJISRI0fa9f6ssXWmh7UshIkpoDmZnEwIxl1dC2cyNBoNUevWFcliNDcYCMivr7Jejy/3hjoKrzAxZUiGDh1q18oYe4ZIbt26xezZswFjUObOo8m7detGaGgoOp2uVNcrisK8efPccjpvWXz11VekpKS4NSiDe6s/Nm/ezH/+8x+31etJFSKIWLhwIZmZmYSGhvL66697ujk2JSQksH37dlQqlUdSqyNGjCAtLa1UBx8ZDAZ+/PFHrx4TBeME1aysLFq2bEn//v3dWndwcDB16tThwIEDpbp+1apVBQJcyyxEgWwAJWcjNBoNG6KjqYb1QCAIeMSi3H9zL6jQAEcoGHyYVMc4d8IHY6BS1tM2LQOFwvMubGUhCr+PtgAGQ5FMxqxZswgwGArcb/neTPNAFECFMbjaEBVlvt/UthAwz0kpzYTO9PR0evbsWeqgY+nSpaSnp1OvXj0SEhLcOjfDtGNvaSdWHzlyhIkTJ7q4VWVz69Yt8yqJyZMnu3UJanh4OC+//DKAeZJseVfug4jr16+bl3XNmjXLq/cxMKXfBg4c6LFNm27evFmqTyXr1q1zyd4VzpSammo+5n327NkuOTOgJI0aNeLxxx8vccnnlStXePnllwukXQtnIUxKk42YMWMGKAofYHuTqSMYO1XLCZbVKBhQWGPqvB2djGtp1qxZGHQ6lgPHuTdMMXPmTJtZiMLv49v89sRw73loNBo2btjAB/lnhZgUfm+m9zLXom5TNsIUxHS1Un5xIiIiOHv2LNUpOei4fv26ubxmzZqxceNGc+DhrgmfvXr14pdffuHu3bslXpuRkeH1Hxw+//xz0tPTCQ0NZfjw4W6vf8aMGfj6+rJ7925iY2PdXr+7lfsg4uOPP+bWrVuEh4fTt29fTzfHpn379rFv3z78/f09Op7WvXt3lixZUuw1iqLQo0cPr/9l8tFHH6HVauncuTN/+9vfPNKGgIAA9uzZU+JEq6tXr/Lf//7X/LW1LIRJSdmI4rIQJqZAZC5FsxCWAYU1ps7bkU2fLGk0GqLWr0cHpAGDuDdMER0dXWwWovD7CIMCmYwJEyag6HT0s6zPynszPcsVGPfMMGUjEhMTWb1qFSP0ejZwL3NR0vJSjUZjnm9xB+NwSXFBx7x588jOziYsLIz9sbEFsh2ObC5WWGkDkS5dupS4Xf3y5cvp2LGjw21xB8ugbMaMGR7Zgr9Ro0bmeU0TJ0506ZH23qBcBxFpaWl88cUXgHuWSjpKURTzzOzhw4d7fBfNUaNGFbsRzaJFi7x+Rcbx48fZsGEDYMxCeDLgeeGFF8jOziYvL8/qz5OSksjKyiowQdVWFsKkuGxEcVkIE1PnuRz7shAF6ndg0ydLM2bMwJA/NNIYSAYSyF8lYjAUm4WwfB+WWRXTxMnt27aRB1iGw7bem2VAZXpfw4YNI1BRyMTKxNRi3vfMmTPN8y0CgIcNBptBR2pqKsuXLweMHY/lvBBTJsbRHTpNShuI1K1bl3nz5tns8HJzc3n99dedduaEq8ybN4+cnBxatWplHlbwhClTpuDv78/hw4edus22NyrXB3ANGzaM1atX07FjR68ev9++fTuvvvoqgYGBnDx50uObYOl0Or744gvGjBlT5GcpKSnUrVvX63+ZvPDCC+zZs4e+ffsSFRXl6eZw+PBhwsLCrO7up9VqycvLMz/TxMREOnXsWOTUzcJm57/GAJ/5+XHi1CkA/qdpU6opCqnYDiLAGIDUBioDvwM5QCMgooR6C9SfX6+9m3dpNBr+p2lTUBSCgBMY52iEAl2BT4F+JbTfJBvYDLwPLAIeA44BDYHLwPn864p7b6ZneR5Yg/F00zeBSCv32HrfGo2G5k2bMkNR+NCizMo+Pgx4880iS36HDRvGhg0bePLJJ0k4fJgInc5833QfH6qrVBzX6wn39eWVkSPtXjJsOjelvk7H5VL8ORkMBpKTk63ucbJ9+3YeffRRr16JlZqaSnh4OFqtlh07dvD00097tD1Tp07lk08+4eGHHyY5OdlrP8SWVfl8VxiXSn711VeAe5dK2kun05knUbprKVJJ1Go1b7zxhtUJgenp6aWeye0psbGx7NmzB7Va7TUbvzzxxBOsWLECvV5f4PsZGRnG1L1FUPb2228Xm4UwGY0xANgG3M3LY/HixcbZ4SVkIUz+BPTcW+Y5H9BinJswohSvJIt67WXKQsC9oYIJGLMRKqAN8D2ws2ZNjrVqZX4lPPII64CfgF1AFHAB4+TKBhizEce4t1TV2pCNNYWzEUHArzbusZWNMGUhCq+ksZaNOH78OF9//TVgzAJYDt1YZmIc2aHTxNaqF1tUKpXVJcnnzp2jRYsWXh1AgPH3vFarpWvXrm7ZX6ck48ePp2rVqpw6dcr8Z10eldsgYurUqRgMBp555hnat2/v6ebYFB0dbd60x51LkUpSpUqVIhvR/PLLLwQFBXn1Rl2Kopi3YB42bBiNGzf2cIvuGT16dIFhIkVROHnyZIF9ITQaDadOnLA6F6Kw6hiXZF5UqXg0LIz69euTkpKCgdIFAj0o2EmGYOyMUzEGCKbXUYyd9S6L/99ZsyYXW7Xif1u1on79+nY9B9OcDR8osnw1CNgHxGP89P/frCy+/uYbfj14kF8PHqTdU08R6OvLbuA6MDX/2jhgLMYMAtxbqmraXOsLSp7nYTk34k3glI17LPfLsFzJ8XVUlNWVNEcoOjfCdMjWs88+y47vviswdGPar8PRHTpN7bG16sUWlUpFx44dWbVqVYHvBwUFeX3mMSkpydxRf/TRR17xobFGjRrmnWenTZuGVqv1cItco1wOZyQmJppPvzt06FCZtiB2pTt37tCyZUsuX77Mv//9b68KIgCSk5O5ePEivXv3Rq/Xo9PpUKlUHpmsVFrffPMNgwYNIigoiBMnTnhFZsfSkiVLGD3a2D3odDoSExML7AfywQcfsHTpUrtT+e+++y7z58/ns88+45tSfOrRarUknzzJTOwYuuBeut/RYQwwLiWOiorCB4oM2czGGAgcApoCDS1S+eZjzXU60oFo4BIFn9NTGAOotPzvZ2HMUNwBMih5eKch8Fr+9ZssyrF6rUXbRo4cyXfr1xe53lRmM+BI/jNLSUnhmWeeQa1W069fP/5v0yYu5QcRGqwPu9g7dDR27Fg2rVxpLrdwe4uTm5uLj48PlSpV4uLFixw6dIhXXnmlxDo96R//+Ae7d+/2muFLk1u3bhEWFsa1a9dYunSp1++r44hyGUT07NmTH3/8kX79+rF27VpPN8emzz77jEmTJvHQQw+RnJzs1o2QSkNRFP788098fHzQaDScOXOG3r17e7pZNuXl5dGqVSsuXLjA1KlTvXKdtlar5eeff6ZXr17MmzevyJr70gYBhb306qtW57DY4miw0gLjSgjLwMUehedCWOt0G2CcG3GMgp3nggUL2LRyJb/p9TxK0Y62uA54FnARKKn7NQVLBozBTIlzUvz82PvLL3R56inzXAhbZZrmRhw+fJjExERee+01Nm/caJ4LAcasg7XgyJ4gwDLYciQQOXfuHKdOnaJ3795cvXqV2rVre8Une1v2799Pjx49UKvVHDlyxKUHlzli2bJljBs3jrp165KSkkJAQICnm+RU5SqIUBSF3bt306tXL3x9fYmLi/PacbycnBzatGnDjRs3WLx4Ma+99pqnm2RVXFwcQUFB3Lx5k6eeesrTzSlWZGQkkyZNIjg4mPj4eK88olhRFM6cOUOdOnWoWrWqR/auAOMvtu1btgDGeRlZN28WuUbBeMy7QVHwU6vxU6upWqMGtYODAXjh5Zd5++237ar3/ffft5mFMDFlI3oDlTAGLJYdrq0sRHEdcAjGAOLJEtpnCpb8KWXmwteX4MaNST97tvisBcZsRKKvL3l6PYGBgTz//POlykKYmIKAhP/8p9gDzyZPnszWNWvM5RZu70tDh/Kvf/2rmHdmnECdl5fHyZMnvXppvKIoPPPMMyQmJjJkyBCvPOpdq9XSrl07fv/9d2bOnMm0adO8Oiizm1KO5OTkKBh/98lLXvLy0pcPKNVAyQJFsfLKzP85oKjzX6Z7nwPFF5TZhe65DIq/le+bXh/l16sGpYqNlx8oAfn//buNcmyVO7EU11UCJajQc7Bs77ug1Cjlcynp+Zb0HDz9d6Aiv9LS0jzdVTpVucpE3Lhxg5o1a3q6GUKIYhSXhTAxZSMs9/lUY1xNolB0+edvGD/J/47tbEADjMtYreV+TOU2Bc5ZKd8WU+ZiFLC0mOsssxGHMb6vatwbzikpC2Fi7bkUZlmutXY0yG+38IzTp0/TrFkzTzfDacpVEKEoCrdv3yY7O9uth9g4Kicnh6CgoPKV2vKwv8qf/V9Fbm4uarW6wKFgjho/fjwrly+3OheiMFNn16BFC+Li4szf/3zJErZu2lTg2tJOEjUNBySdOMGDDz5Y4GemctOvXSPbytCONYqikHvnDg8C72FcHVKc2cBMQAf4qFTMsphDYWsopjDTkET/ESPMB8uZXLlyhRZhYUXmQlhrh63n8FeUl5eHTqfz6oMV4d7vpoCAgHL1O79cBRFCCO+k0WgIbdAAg8FQYhbCxLxSIz7evNrKmtGjRxO9bFmROQCFmTrg10aN4rPPPrOr/da8//77fPrpp3ZnLQKrVOFWbq75vlv537drpYyfH+cvXiwwQdJTz0FUbBJECCFcztThgv1DBa1bt7Z5rLJGo6FRaGiJn75NbHXAjli0cCGboqPtvq9OvXqkW2zqlJqWRnpGhv07dL7/PosWLQI8+xxExea9R1oKIcqNkAYNuK9aNXJzcvjOjvsC1WoaFLOZ1fz589HqdObNtUqSjXGXzQULFpg7YEeNHTeOsePGlakMuBeMpNpxT1so8Fw8+RxExSZBhBDC5ZzV4RYW0qABbR97rEwdsKc549mUh+cg/ppkOEMIIYQQDim3Z2cIIYQQwrUkiBBCCCGEQySIEEIIIYRDJIgQQgghhEMkiBBCCCGEQySIEEIIIYRDJIgQQgghhEMkiBBCCCGEQ1weRKhUKrp06eLqaoQQQgjhZmUKIgYPHoxKpaJOnTrodDpntUkIIYQQdvBUf+xwEJGdnc3WrVtRqVRcvXqVnTt3OrNdQgghhCgFT/bHDgcRX3/9Nbdv32b8+PGoVCoiIyOd2S4hhBBClIIn+2OHg4jIyEj8/f2ZMmUKHTp0YNeuXfzxxx82r09LS6N///488MADBAYG0qVLF3777TdHqxdCCCEEpe+PTXMUMzIyGDp0KLVq1aJKlSq0a9eOmJgYh+p2KIhISkoiPj6eZ599lvvvv5/Bgwej1+tZu3at1etv3rxJhw4duHTpEiNHjqRv374cPHiQrl27OtxwIYQQoqKztz/OzMykQ4cOHD9+nIEDB/Liiy+SkJBAz549SU5Otr8BigPee+89BVC+/fZbRVEUJTMzU6lcubLStGnTItcCCqAMGjRIMRgM5u/HxMQoKpVKadKkiaLX6x1phhBCCFGhOdIfjxo1qkC/u2rVKgVQ3nzzTbvrV+UXXGparZYHH3wQg8FAeno6/v7+ALzyyits2rSJ2NhYOnXqZL5epVLh6+vLxYsXqV+/foGynnvuOXbu3Mn+/fvp2LGj/RGQEEIIUUE50h8HBgaSnp5OUFCQ+fs6nY4qVaoQHh5OYmKiXW2wezhj+/bt3Lhxg/79+5sbDMblJQCrV68uck9ISEiRAALgqaeeAuDo0aP2NkMIIYSo0Bzpj5s2bVoggABQq9XUrl2bzMxMu9tgdxBhatSgQYMKfL9nz57UqVOHLVu2kJ2dXeBntWrVslpW7dq1AcjKyrK3GUIIIUSF5kh/XL16datlqdVq9Hq93W2wK4hIS0vjp59+AqBDhw6oVCrzS61Wk56ezu3bt9m4cWOB+65du2a1vKtXrwK235QQQgghinK0P3Y2tT0Xr1mzBoPBQMeOHWnWrFmRn2u1WtavX09kZCQjR440fz81NZW0tLQiQxr79+8HoFWrVg40XQghhKiYHO2Pna3UEysVRaFRo0akpqZy/vx5GjZsaPW6xx57jCNHjpCUlESLFi1QqVSAMd2ydu1a89exsbF07dqVxo0bc+bMGXx85CwwIYQQoiRl6Y87d+5sdWuF0NBQAC5dumRXW0rdc+/du5dLly7RpUsXmw0GGDJkCECBHbPCw8OJiYmhXbt2TJkyhddff50ePXrg5+fHypUrJYAQQgghSqks/bGzlbr3NjVi6NChxV43YMAA/P39iYqKQqvVAlCjRg0OHDhASEgIy5cvZ8uWLbRr1459+/bJCZ9CCCGEHcrSHzub3ftECCGEEEJAGY8CF0IIIUTFJUGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEcIkGEEEIIIRwiQYQQQgghHCJBhBBCCCEc8v8CrVXs2vm5pwAAAABJRU5ErkJggg==\n",
"text/plain": [
"
"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Plag_tern_points = pt.tern_points_fspar(fspar_comps=Plags)\n",
"\n",
"fig, tax = pt.plot_fspar_classification(major_grid=True, ticks=False)\n",
"\n",
"# Here we plot the data ontop, with plag as red triangles\n",
"tax.scatter(\n",
" Plag_tern_points,\n",
" edgecolor=\"k\",\n",
" marker=\"^\",\n",
" facecolor=\"red\",\n",
" label='Plag',\n",
" s=90\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "b719fdfc-fa9d-4210-9b42-d696fd17da70",
"metadata": {},
"outputs": [],
"source": [
"Plag_components=pt.calculate_cat_fractions_plagioclase(plag_comps=Plags)\n",
"RealPlags=Plags.loc[Plag_components['An_Plag']>0.05].reset_index(drop=True)"
]
},
{
"cell_type": "markdown",
"id": "55b30af5-91f7-478b-9385-0814b4737a9b",
"metadata": {},
"source": [
"## Example 1a - All possible Liq-Plag matches, An-Ab equilibrium\n",
"- You can specify to use the An-Ab equilibrium test of Putirka (2008) as in this example"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "890ee427-77d6-4c51-97a7-e552676ca94f",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Considering N=28 Fspar & N=24 Liqs, which is a total of N=672 Liq-Fspar pairs, be patient if this is >>1 million!\n",
"Applying filter to only average those that pass the An-Ab eq test of Putirka, 2008\n",
"Done!!! I found a total of N=327 Fspar-Liq matches using the specified filter. N=28 Fspar out of the N=28 Fspar that you input matched to 1 or more liquids\n"
]
}
],
"source": [
"MM_dict=pt.calculate_fspar_liq_temp_matching(liq_comps=Liqs, plag_comps=RealPlags, \n",
" equationT=\"T_Put2008_eq24a\", P=5, \n",
" Ab_An_P2008=True, H2O_Liq=2)\n",
"Av_Matches=MM_dict['Av_PTs']\n",
"All_Matches=MM_dict['All_PTs']"
]
},
{
"cell_type": "markdown",
"id": "aa21261f-0f16-42cd-a4f8-96f33ffd8007",
"metadata": {},
"source": [
"## Example 1b - All possible Liq-Plag matches, no filter\n",
"- This could be used if you want to develop your own equilibrium filters"
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "e6f3f0db-56b2-4d3f-a04a-d6d5ea9f0a36",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Considering N=28 Fspar & N=24 Liqs, which is a total of N=672 Liq-Fspar pairs, be patient if this is >>1 million!\n",
"We are returning all pairs, if you want to use the Ab-An equilibrium test of Putirka (2008), enter Ab_An_P2008=True\n",
"Done!!! I found a total of N=672 Fspar-Liq matches using the specified filter. N=28 Fspar out of the N=28 Fspar that you input matched to 1 or more liquids\n"
]
}
],
"source": [
"MM_dict=pt.calculate_fspar_liq_temp_matching(liq_comps=Liqs, plag_comps=RealPlags, \n",
" equationT=\"T_Put2008_eq24a\", P=5, \n",
" Ab_An_P2008=False, H2O_Liq=2)\n",
"Av_Matches=MM_dict['Av_PTs']\n",
"All_Matches=MM_dict['All_PTs']"
]
},
{
"cell_type": "markdown",
"id": "4a04260c-c993-4b1e-b770-82abb15ab18c",
"metadata": {},
"source": [
"## Example 1c - Iterating T and H2O for all possible Plag-Liq matches"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "43ef6cf7-b62f-4faf-b925-e5577c78d66f",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Considering N=28 Fspar & N=24 Liqs, which is a total of N=672 Liq- Fspar pairs, be patient if this is >>1 million!\n"
]
},
{
"name": "stderr",
"output_type": "stream",
"text": [
"100%|██████████████████████████████████████████████████████████████████████████████████| 20/20 [00:02<00:00, 9.58it/s]\n"
]
},
{
"name": "stdout",
"output_type": "stream",
"text": [
"Applying filter to only average those that pass the An-Ab eq test of Putirka, 2008\n",
"Done!!! I found a total of N=327 Fspar-Liq matches using the specified filter. N=28 Fspar out of the N=28 Fspar that you input matched to 1 or more liquids\n"
]
}
],
"source": [
"T_H_Iter_dict=pt.calculate_fspar_liq_temp_hygr_matching(liq_comps=Liqs, plag_comps=RealPlags,\n",
" equationT=\"T_Put2008_eq24a\", equationH=\"H_Waters2015\", P=5)\n",
"T_H_Calc_Av=T_H_Iter_dict.get('Av_HTs')\n",
"T_H_Calc_All=T_H_Iter_dict.get('All_HTs')\n",
"T_H_Evol=T_H_Iter_dict.get('T_H_Evolution')"
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "a322a364-017a-4a2d-8536-f5623f4a7278",
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"Text(0.5, 0, 'H$_2$O Calc (wt%)')"
]
},
"execution_count": 9,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAh8AAAG0CAYAAACSbkVhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAA9hAAAPYQGoP6dpAAAlGklEQVR4nO3de1TUdf7H8dcoMpgLY2qgKEilpmlqaaVlgV285L1tszQvu9W6K1lJW+qvLO2Gtl2oSMuOmZ3U/CM1i7Iswct6yRutluUlWNmUyEvgrZHg8/tjD3OaAJXhOx8YfT7OmXOc73xnvu+PA/I83xkclzHGCAAAwJI6NT0AAAA4txAfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsCqspgf4vdLSUu3bt0+RkZFyuVw1PQ4AADgDxhgdOXJEsbGxqlPn1Oc2al187Nu3T3FxcTU9BgAACEBeXp5atGhxyn1qXXxERkZK+t/wUVFRNTwNAAA4E0VFRYqLi/P9HD+VWhcfZS+1REVFER8AAISYM3nLBG84BQAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAqirHx6pVqzRgwADFxsbK5XJpyZIlvtuKi4s1YcIEXXbZZWrQoIFiY2M1cuRI7du3z8mZAQBACKtyfBw7dkydOnVSenp6uduOHz+uLVu2aPLkydqyZYsWLVqknTt3auDAgY4MCwAAQp/LGGMCvrPLpcWLF2vw4MGV7rNx40ZdddVV+s9//qP4+PjTPmZRUZE8Ho8KCwv5YDkAAEJEVX5+B/1TbQsLC+VyudSwYcMKb/d6vfJ6vb7rRUVFwR4JAADUoKDGxy+//KKJEydq2LBhlVZQamqqpk6dGswxgLNWwsSMmh6hynKn9avpEQDUsKD9tktxcbHuuOMOlZaWasaMGZXuN2nSJBUWFvoueXl5wRoJAADUAkE581FcXKzbb79dOTk5WrFixSlf+3G73XK73cEYAwAA1EKOx0dZeOzatUuZmZlq3Lix04cAAAAhrMrxcfToUe3evdt3PScnR9nZ2WrUqJFiY2N12223acuWLfroo49UUlKi/Px8SVKjRo0UHh7u3OQAACAkVTk+Nm3apJ49e/qup6SkSJJGjRqlKVOmaOnSpZKkzp07+90vMzNTSUlJgU8KAADOClWOj6SkJJ3qvwapxn8bAgAAzgF8tgsAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVVY6PVatWacCAAYqNjZXL5dKSJUv8bjfGaMqUKYqNjVX9+vWVlJSkr7/+2ql5AQBAiKtyfBw7dkydOnVSenp6hbc/99xzevHFF5Wenq6NGzeqadOmuvnmm3XkyJFqDwsAAEJfWFXv0LdvX/Xt27fC24wxSktL06OPPqpbb71VkjR37lzFxMRo/vz5GjNmTPWmBQAAIc/R93zk5OQoPz9fvXr18m1zu91KTEzU2rVrK7yP1+tVUVGR3wUAAJy9HI2P/Px8SVJMTIzf9piYGN9tv5eamiqPx+O7xMXFOTkSAACoZYLy2y4ul8vvujGm3LYykyZNUmFhoe+Sl5cXjJEAAEAtUeX3fJxK06ZNJf3vDEizZs182wsKCsqdDSnjdrvldrudHAMAANRijp75uPDCC9W0aVMtX77ct+3kyZNauXKlrrnmGicPBQAAQlSVz3wcPXpUu3fv9l3PyclRdna2GjVqpPj4eD344IN69tln1bp1a7Vu3VrPPvuszjvvPA0bNszRwQEAQGiqcnxs2rRJPXv29F1PSUmRJI0aNUpvv/22HnnkEZ04cUJjx47V4cOHdfXVV+uzzz5TZGSkc1MDAICQ5TLGmJoe4reKiork8XhUWFioqKiomh4HqNUSJmbU9AhVljutX02PACAIqvLzm892AQAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwKq+kBgNoiYWJGTY8AAOcEznwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABY5Xh8/Prrr3rsscd04YUXqn79+rrooov05JNPqrS01OlDAQCAEBTm9ANOnz5dr7/+uubOnav27dtr06ZN+vOf/yyPx6MHHnjA6cMBAIAQ43h8rFu3ToMGDVK/fv0kSQkJCVqwYIE2bdrk9KEAAEAIcvxllx49euiLL77Qzp07JUlfffWV1qxZo1tuucXpQwEAgBDk+JmPCRMmqLCwUG3btlXdunVVUlKiZ555RnfeeWeF+3u9Xnm9Xt/1oqIip0cCAAC1iONnPhYuXKh3331X8+fP15YtWzR37lw9//zzmjt3boX7p6amyuPx+C5xcXFOjwQAAGoRlzHGOPmAcXFxmjhxopKTk33bnn76ab377rv69ttvy+1f0ZmPuLg4FRYWKioqysnRgFNKmJhR0yOcE3Kn9avpEQAEQVFRkTwezxn9/Hb8ZZfjx4+rTh3/Eyp169at9Fdt3W633G6302MAAIBayvH4GDBggJ555hnFx8erffv22rp1q1588UX95S9/cfpQAAAgBDkeH6+++qomT56ssWPHqqCgQLGxsRozZowef/xxpw8FAABCkOPxERkZqbS0NKWlpTn90AAA4CzAZ7sAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYRXwAAACriA8AAGAV8QEAAKwiPgAAgFXEBwAAsIr4AAAAVhEfAADAKuIDAABYFZT4+OGHH3TXXXepcePGOu+889S5c2dt3rw5GIcCAAAhJszpBzx8+LCuvfZa9ezZU5988omio6O1Z88eNWzY0OlDAQCAEOR4fEyfPl1xcXGaM2eOb1tCQoLThwEAACHK8Zddli5dqq5du+pPf/qToqOjdfnll+vNN9+sdH+v16uioiK/CwAAOHs5Hh/ff/+9Zs6cqdatW+vTTz/V3/72N91///165513Ktw/NTVVHo/Hd4mLi3N6JAAAUIu4jDHGyQcMDw9X165dtXbtWt+2+++/Xxs3btS6devK7e/1euX1en3Xi4qKFBcXp8LCQkVFRTk5GnBKCRMzanqEc0LutH41PQKAICgqKpLH4zmjn9+On/lo1qyZLr30Ur9t7dq10969eyvc3+12Kyoqyu8CAADOXo7Hx7XXXqvvvvvOb9vOnTvVsmVLpw8FAABCkOPxMX78eK1fv17PPvusdu/erfnz52vWrFlKTk52+lAAACAEOR4fV155pRYvXqwFCxaoQ4cOeuqpp5SWlqbhw4c7fSgAABCCHP9/PiSpf//+6t+/fzAeGgAAhDg+2wUAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMCqoMdHamqqXC6XHnzwwWAfCgAAhICgxsfGjRs1a9YsdezYMZiHAQAAISRo8XH06FENHz5cb775ps4///xgHQYAAISYoMVHcnKy+vXrp5tuuumU+3m9XhUVFfldAADA2SssGA/63nvvafPmzdq0adNp901NTdXUqVODMQYAAKiFHD/zkZeXpwceeEDz5s1TRETEafefNGmSCgsLfZe8vDynRwIAALWI42c+Nm/erIKCAnXp0sW3raSkRKtWrVJ6erq8Xq/q1q3ru83tdsvtdjs9BgAAqKUcj48bb7xR27Zt89v25z//WW3bttWECRP8wgMAAJx7HI+PyMhIdejQwW9bgwYN1Lhx43LbAQDAuYf/4RQAAFgVlN92+b2srCwbhwEAACGAMx8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwyvH4SE1N1ZVXXqnIyEhFR0dr8ODB+u6775w+DAAACFGOx8fKlSuVnJys9evXa/ny5fr111/Vq1cvHTt2zOlDAQCAEBTm9AMuW7bM7/qcOXMUHR2tzZs36/rrr3f6cAAAIMQE/T0fhYWFkqRGjRoF+1AAACAEOH7m47eMMUpJSVGPHj3UoUOHCvfxer3yer2+60VFRcEcCQAA1LCgxsd9992nf//731qzZk2l+6Smpmrq1KnBHANALZIwMaOmRwhI7rR+NT3COSEUvz742qi6oL3sMm7cOC1dulSZmZlq0aJFpftNmjRJhYWFvkteXl6wRgIAALWA42c+jDEaN26cFi9erKysLF144YWn3N/tdsvtdjs9BgAAqKUcj4/k5GTNnz9fH3zwgSIjI5Wfny9J8ng8ql+/vtOHAwAAIcbxl11mzpypwsJCJSUlqVmzZr7LwoULnT4UAAAIQUF52QUAAKAyfLYLAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArCI+AACAVcQHAACwivgAAABWER8AAMAq4gMAAFhFfAAAAKuIDwAAYBXxAQAArAqr6QFsS5iYUdMjVFnutH41PUKVheLfMwAEIhT/vavpnyuc+QAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALAqaPExY8YMXXjhhYqIiFCXLl20evXqYB0KAACEkKDEx8KFC/Xggw/q0Ucf1datW3Xdddepb9++2rt3bzAOBwAAQkhQ4uPFF1/U3XffrXvuuUft2rVTWlqa4uLiNHPmzGAcDgAAhJAwpx/w5MmT2rx5syZOnOi3vVevXlq7dm25/b1er7xer+96YWGhJKmoqMjp0SRJpd7jQXncYArW30UwheLfM3Aqofh9GIr4t8OOYHw9lz2mMea0+zoeHwcOHFBJSYliYmL8tsfExCg/P7/c/qmpqZo6dWq57XFxcU6PFrI8aTU9AQC+D3E2CebX85EjR+TxeE65j+PxUcblcvldN8aU2yZJkyZNUkpKiu96aWmpDh06pMaNG1e4f21SVFSkuLg45eXlKSoqqqbHCSrWenZirWcn1nr2qs3rNcboyJEjio2NPe2+jsdHkyZNVLdu3XJnOQoKCsqdDZEkt9stt9vtt61hw4ZOjxVUUVFRte6LIFhY69mJtZ6dWOvZq7au93RnPMo4/obT8PBwdenSRcuXL/fbvnz5cl1zzTVOHw4AAISYoLzskpKSohEjRqhr167q3r27Zs2apb179+pvf/tbMA4HAABCSFDiY+jQoTp48KCefPJJ7d+/Xx06dNDHH3+sli1bBuNwNcbtduuJJ54o97LR2Yi1np1Y69mJtZ69zpb1usyZ/E4MAACAQ/hsFwAAYBXxAQAArCI+AACAVcQHAACwiviogsOHD2vEiBHyeDzyeDwaMWKEfv7551Pe5+jRo7rvvvvUokUL1a9fX+3atQuJD9gLZK0ul6vCyz//+U87QwcokLVK0o4dOzRw4EB5PB5FRkaqW7duIfHJzYGsd/To0eWe127dutkZuBoCfW7LjBkzRi6XS2lpaUGb0SmBrHXKlClq27atGjRooPPPP1833XSTNmzYYGfgaqjqWouLizVhwgRddtllatCggWJjYzVy5Ejt27fP3tABCuR5XbRokXr37q0mTZrI5XIpOzvbyqxVQXxUwbBhw5Sdna1ly5Zp2bJlys7O1ogRI055n/Hjx2vZsmV69913tWPHDo0fP17jxo3TBx98YGnqwASy1v379/td3nrrLblcLv3xj3+0NHVgAlnrnj171KNHD7Vt21ZZWVn66quvNHnyZEVERFiaOnCBrFeS+vTp4/f8fvzxxxamrZ5A1ypJS5Ys0YYNG87ov4quDQJZa5s2bZSenq5t27ZpzZo1SkhIUK9evfTTTz9ZmjowVV3r8ePHtWXLFk2ePFlbtmzRokWLtHPnTg0cONDi1IEJ5Hk9duyYrr32Wk2bNs3SlAEwOCPffPONkWTWr1/v27Zu3TojyXz77beV3q99+/bmySef9Nt2xRVXmMceeyxos1ZXoGv9vUGDBpkbbrghGCM6JtC1Dh061Nx11102RnRUoOsdNWqUGTRokIUJnVOdr+P//ve/pnnz5mb79u2mZcuW5qWXXgrytNXj1PdsYWGhkWQ+//zzYIzpCKfW+uWXXxpJ5j//+U8wxnREddeak5NjJJmtW7cGccrAcObjDK1bt04ej0dXX321b1u3bt3k8Xi0du3aSu/Xo0cPLV26VD/88IOMMcrMzNTOnTvVu3dvG2MHJNC1/taPP/6ojIwM3X333cEa0xGBrLW0tFQZGRlq06aNevfurejoaF199dVasmSJpakDV53nNisrS9HR0WrTpo3uvfdeFRQUBHvcagl0raWlpRoxYoQefvhhtW/f3sao1ebE9+zJkyc1a9YseTwederUKVijVpsTa5WkwsJCuVyuWv1ZYk6ttTYiPs5Qfn6+oqOjy22Pjo4u9yF6v/XKK6/o0ksvVYsWLRQeHq4+ffpoxowZ6tGjRzDHrZZA1/pbc+fOVWRkpG699Vanx3NUIGstKCjQ0aNHNW3aNPXp00efffaZhgwZoltvvVUrV64M9sjVEuhz27dvX82bN08rVqzQCy+8oI0bN+qGG26Q1+sN5rjVEuhap0+frrCwMN1///3BHM9R1fme/eijj/SHP/xBEREReumll7R8+XI1adIkWKNWmxP/Pv3yyy+aOHGihg0bVis/nK2ME2utrc75+JgyZUqlb5Qsu2zatEnS/95Q+XvGmAq3l3nllVe0fv16LV26VJs3b9YLL7ygsWPH6vPPPw/amioT7LX+1ltvvaXhw4fX2HsggrnW0tJSSdKgQYM0fvx4de7cWRMnTlT//v31+uuvB29RpxDs53bo0KHq16+fOnTooAEDBuiTTz7Rzp07lZGREbQ1VSaYa928ebNefvllvf3222f8tR5MNr5ne/bsqezsbK1du1Z9+vTR7bffXiNntWz9+1RcXKw77rhDpaWlmjFjhuPrOBM2/y2urYLy2S6h5L777tMdd9xxyn0SEhL073//Wz/++GO523766SfFxMRUeL8TJ07o//7v/7R48WL169dPktSxY0dlZ2fr+eef10033VT9BVRBMNf6W6tXr9Z3332nhQsXBjxrdQVzrU2aNFFYWJguvfRSv+3t2rXTmjVrAh+6Gmw9t2WaNWumli1bateuXVWetbqCudbVq1eroKBA8fHxvm0lJSV66KGHlJaWptzc3GrNXlU2ntcGDRqoVatWatWqlbp166bWrVtr9uzZmjRpUrVmryobay0uLtbtt9+unJwcrVixosbOetj+fq2VauzdJiGm7I0/GzZs8G1bv379Kd/4U/bmrY8//thv+1//+ldz8803B3Xe6ghkrb81atQo06VLl2CO6JhA19q9e/dybzgdPHiwufPOO4M2qxOq+9yWOXDggHG73Wbu3LnBGNMRgaz1wIEDZtu2bX6X2NhYM2HChCr9/djm1PNqjDEXX3yxeeKJJxye0DmBrvXkyZNm8ODBpn379qagoMDGqNVW3ee1Nr/hlPiogj59+piOHTuadevWmXXr1pnLLrvM9O/f32+fSy65xCxatMh3PTEx0bRv395kZmaa77//3syZM8dERESYGTNm2B6/SgJZqzH/C67zzjvPzJw50+a41RLIWhctWmTq1atnZs2aZXbt2mVeffVVU7duXbN69Wrb41dZVdd75MgR89BDD5m1a9eanJwck5mZabp3726aN29uioqKamIJZyzQr+PfCoXfdjGm6ms9evSomTRpklm3bp3Jzc01mzdvNnfffbdxu91m+/btNbGEM1bVtRYXF5uBAweaFi1amOzsbLN//37fxev11sQSzlggX8MHDx40W7duNRkZGUaSee+998zWrVvN/v37bY9fKeKjCg4ePGiGDx9uIiMjTWRkpBk+fLg5fPiw3z6SzJw5c3zX9+/fb0aPHm1iY2NNRESEueSSS8wLL7xgSktL7Q5fRYGs1Rhj3njjDVO/fn3z888/2xu2mgJd6+zZs02rVq1MRESE6dSpk1myZIm9oauhqus9fvy46dWrl7ngggtMvXr1THx8vBk1apTZu3ev/eGrKNDn9rdCJT6qutYTJ06YIUOGmNjYWBMeHm6aNWtmBg4caL788kv7w1dRVddadgagoktmZqb1+asikK/hOXPmVLjW2nRGy2WMMcF8WQcAAOC3zvnfdgEAAHYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwCsSEpK0oMPPlitxzh48KCio6Otf7psRW677Ta9+OKLNT0GEJKID6CWGz16tAYPHlxue1ZWllwul37++eczepy8vDzdfffdio2NVXh4uFq2bKkHHnhABw8ePO198/PzNW7cOF100UVyu92Ki4vTgAED9MUXX1RxNdWTmpqqAQMGKCEhIaD7VxRA8+bNU1xcnBo1aqSHH37Y77bc3Fy1adNGRUVF5R7r8ccf1zPPPFPhbQBOjfgAzgHff/+9unbtqp07d2rBggXavXu3Xn/9dX3xxRfq3r27Dh06VOl9c3Nz1aVLF61YsULPPfectm3bpmXLlqlnz55KTk62toYTJ05o9uzZuueeexx7zAMHDuiee+7R888/r08//VRz585VRkaG7/a///3vmjZtmqKiosrdt2PHjkpISNC8efMcmwc4VxAfwDkgOTlZ4eHh+uyzz5SYmKj4+Hj17dtXn3/+uX744Qc9+uijld537Nixcrlc+vLLL3XbbbepTZs2at++vVJSUrR+/XrffsuWLVOPHj3UsGFDNW7cWP3799eePXsqfdzS0lJNnz5drVq1ktvtVnx8vJ555plK9//kk08UFham7t27+7Z9+OGHatiwoUpLSyVJ2dnZcrlcfmcwxowZozvvvFOjR4/WypUr9fLLL8vlcsnlcun777+Xx+PR0KFDdeWVV6pnz5765ptvJEnz589XeHi4br311kpnGjhwoBYsWFDp7QAqRnwAZ7lDhw7p008/1dixY1W/fn2/25o2barhw4dr4cKFqugDrg8dOqRly5YpOTlZDRo0KHd7w4YNfX8+duyYUlJStHHjRn3xxReqU6eOhgwZ4guD35s0aZKmT5+uyZMn65tvvtH8+fMVExNT6TpWrVqlrl27+m27/vrrdeTIEW3dulWStHLlSjVp0kQrV6707ZOVlaXExES9/PLL6t69u+69917t379f+/fv18UXX6zjx49r69atOnTokDZu3KiOHTvq0KFDevzxx5Wenl7pPJJ01VVX6csvv5TX6z3lfgD8hdX0AABO76OPPtIf/vAHv20lJSW+P+fl5WnEiBEqKChQWFiYJk+erD/96U+SpF27dskYo3bt2lX42O3atdPhw4f1008/KTo62u+23bt3yxijtm3bnnbGP/7xj37XZ8+erejoaH3zzTfq0KGD321HjhzRyy+/rPT0dI0aNUqSdPHFF6tHjx6VPn5ubq5iY2P9tnk8HnXu3FlZWVnq0qWLsrKyNH78eE2dOlVHjhzRsWPHtHPnTiUlJcnj8Sg8PFznnXeemjZt6nuMuXPnauTIkTpx4oRGjhyp3r176y9/+YvGjRunnJwcDRw4UMXFxZoyZYpuu+02v+M3b95cXq9X+fn5atmy5Wn/jgD8D/EBhICePXtq5syZfts2bNigu+66S5IUFhamtLQ0de7cWQUFBbriiit0yy23VHi24vfKzni4XK4q3fZ7e/bs0eTJk7V+/XodOHDAd8Zj79695eJjx44d8nq9uvHGG0/7uGVOnDihiIiIctuTkpKUlZWllJQUrV69Wk8//bTef/99rVmzRj///LNiYmJOGU9DhgzRkCFDfNezsrK0bds2paenq1WrVlqwYIGaNm2qq666Stdff71foJWdSTp+/PgZrwMAL7sAIaFBgwZq1aqV36V58+a+25s1a6bOnTtLkqKjo9WoUSPfm0hbtWoll8vley/D73377bc6//zz1aRJk3K3tW7dWi6XSzt27DjtjAMGDNDBgwf15ptvasOGDdqwYYMk6eTJk+X2/f3LP2eiSZMmOnz4cLntSUlJWr16tb766ivVqVNHl156qRITE7Vy5UrfSy5nyuv1auzYsXrjjTe0e/du/frrr0pMTNQll1yiNm3a+NZUpuzv+IILLqjyeoBzGfEBnGU2bdqk0tJSxcXFSZIaN26sm2++WTNmzNCJEyf89s3Pz9e8efM0dOjQCs9uNGrUSL1799Zrr72mY8eOlbu97Nd8Dx48qB07duixxx7TjTfe6HsppzKtW7dW/fr1q/SrupdffnmFAVX2vo+0tDQlJibK5XIpMTFRWVlZ5eIjPDzc7+Wq33vqqafUt29fXXHFFSopKdGvv/7qu624uLjcfbdv364WLVpUGG4AKkd8AGeRgwcPauTIkZo1a5bf9vT0dHm9XvXu3VurVq1SXl6eli1bpptvvlnNmzc/5W+ZzJgxQyUlJbrqqqv0/vvva9euXdqxY4deeeUV32+enH/++WrcuLFmzZql3bt3a8WKFUpJSan0MSMiIjRhwgQ98sgjeuedd7Rnzx6tX79es2fPrvQ+vXv31tdff10uasre9/Huu+8qKSlJ0v+CZMuWLb73e5RJSEjQhg0blJub6/fSkCR9/fXXWrhwoZ588klJUtu2bVWnTh3Nnj1bGRkZ+vbbb3XllVf6HXv16tXq1atXpTMDqIQBUKuNGjXKDBo0qNz2zMxMI8kcPnzYGGPML7/8Yq677jrzzjvvVPg4ubm5ZvTo0aZp06amXr16Ji4uzowbN84cOHDgtDPs27fPJCcnm5YtW5rw8HDTvHlzM3DgQJOZmenbZ/ny5aZdu3bG7Xabjh07mqysLCPJLF682BhjTGJionnggQd8+5eUlJinn37atGzZ0tSrV8/Ex8ebZ5999pRzdOvWzbz++uvltj/00ENGktm+fbtvW6dOncwFF1xgSktLfdu+++47061bN1O/fn0jyeTk5BhjjCktLTXXXHON+fDDD/0e98MPPzTx8fEmJibGvPnmm363nThxwkRFRZl169adcmYA5bmMqeD36wCEFGOMhg0bpksuuURTpkyp6XGC5uOPP9Y//vEPbd++XXXq1OyJ29dee00ffPCBPvvssxqdAwhFvOwCnAX+9a9/aeHChVqyZIk6d+6szp07a9u2bTU9luNuueUWjRkzRj/88ENNj6J69erp1VdfrekxgJDEmQ8AAGAVZz4AAIBVxAcAALCK+AAAAFYRHwAAwCriAwAAWEV8AAAAq4gPAABgFfEBAACsIj4AAIBVxAcAALCK+AAAAFYRHwAAwKr/B3/LQx9ECIPOAAAAAElFTkSuQmCC\n",
"text/plain": [
"
"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"Kspar_tern_points = pt.tern_points_fspar(fspar_comps=Kspars)\n",
"fig, tax = pt.plot_fspar_classification(major_grid=True, ticks=False)\n",
"\n",
"tax.scatter(\n",
" Kspar_tern_points,\n",
" edgecolor=\"k\",\n",
" marker=\"o\",\n",
" facecolor=\"cyan\",\n",
" s=60, label='Kspar'\n",
")"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "53b2cdb3-a9c7-47bf-b3c4-4c0252426cfd",
"metadata": {},
"outputs": [],
"source": [
"Kspar_components=pt.calculate_cat_fractions_kspar(kspar_comps=Kspars)\n",
"RealKspars=Kspars.loc[Kspar_components['An_Kspar']<0.05].reset_index(drop=True)"
]
},
{
"cell_type": "markdown",
"id": "840f2ea2-3ff7-4cf6-93e1-17030c1979ca",
"metadata": {},
"source": [
"### 2a - Lets calculate temperature for all liq-kspar matches\n",
"- Currently no equilibrium tests exist for Kspar-Liquid that we are aware off. If you know of ones I can implement, please reach out!"
]
},
{
"cell_type": "code",
"execution_count": 15,
"id": "16d46a40-34a4-48f6-b45b-c57694258b71",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Considering N=15 Fspar & N=15 Liqs, which is a total of N=225 Liq-Fspar pairs, be patient if this is >>1 million!\n",
"Sorry, no equilibrium tests implemented for Kspar-Liquid. Weve returned all possible pairs, you will have to filter them yourselves\n",
"Done!!! I found a total of N=225 Fspar-Liq matches using the specified filter. N=15 Fspar out of the N=15 Fspar that you input matched to 1 or more liquids\n"
]
},
{
"data": {
"text/html": [
"