{ "cells": [ { "cell_type": "markdown", "id": "sticky-wayne", "metadata": {}, "source": [ "Import required libraries" ] }, { "cell_type": "code", "execution_count": 2, "id": "particular-bikini", "metadata": {}, "outputs": [], "source": [ "import numpy as np\n", "import hydrangea as hy\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "id": "divine-gentleman", "metadata": {}, "source": [ "For illustration, let's look at simulation CE-28. Construct a simulation object to handle file paths etc." ] }, { "cell_type": "code", "execution_count": 4, "id": "previous-variance", "metadata": {}, "outputs": [], "source": [ "sim = hy.Simulation(28)" ] }, { "cell_type": "markdown", "id": "indirect-introduction", "metadata": {}, "source": [ "Construct the z = 0.1 (snapshot 26) snapshot and subfind files " ] }, { "cell_type": "code", "execution_count": 8, "id": "sorted-thermal", "metadata": {}, "outputs": [], "source": [ "snapshot_file = sim.get_snapshot_file(26)\n", "subfind_file = sim.get_subfind_file(26)" ] }, { "cell_type": "markdown", "id": "acceptable-nation", "metadata": {}, "source": [ "Set up a loader to extract the `Subhalo` data in the Subfind catalogue" ] }, { "cell_type": "code", "execution_count": 5, "id": "05260e93", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Prepared reading from 'SUBHALO'...\n" ] } ], "source": [ "sub = hy.SplitFile(subfind_file, 'Subhalo')" ] }, { "cell_type": "markdown", "id": "formed-parking", "metadata": {}, "source": [ "Get e.g. the stellar mass of all subhaloes (warning: might take a bit of time). Units are solar mass." ] }, { "cell_type": "code", "execution_count": 14, "id": "freelance-microphone", "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/software/rhel7/lib64/python3.6/site-packages/ipykernel_launcher.py:2: RuntimeWarning: divide by zero encountered in log10\n", " \n" ] }, { "data": { "text/plain": [ "(array([2120., 1572., 1164., 821., 644., 430., 258., 106., 29.,\n", " 16.]),\n", " array([ 9. , 9.3, 9.6, 9.9, 10.2, 10.5, 10.8, 11.1, 11.4, 11.7, 12. ],\n", " dtype=float32),\n", " )" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAD4CAYAAAAAczaOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAULUlEQVR4nO3df4xd5X3n8fdnScJu00ZAPGWJbWKSdaKSaNeQEWGbJmLFFhxSBVKpWdAqOAmqgwpSolaqnEZaolRIpG0aLdoukVOsgJRA2aUsVkMKDtsNW2mdMFADNj+KTYyw5djTgqApFVvId/+4z2xvzIznztzrOzM+75d0dc/9nuec8zy6+DOH5557T6oKSVI3/LOl7oAkaXwMfUnqEENfkjrE0JekDjH0JalD3rDUHZjPqlWrat26dUvdDUlaMR566KG/qaqJ2dYt+9Bft24dU1NTS90NSVoxkjw71zqndySpQwx9SeoQQ1+SOsTQl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDlv03coexbsu3l+S4+2/4yJIcV5Lm45m+JHWIoS9JHWLoS1KHGPqS1CGGviR1iKEvSR1i6EtSh8wb+knWJvmLJI8n2ZPks61+WpIdSZ5uz6e2epLcmGRvkkeTnNu3r02t/dNJNh2/YUmSZjPImf6rwG9V1dnA+cA1Sc4GtgD3V9V64P72GuDDwPr22AzcBL0/EsB1wPuB84DrZv5QSJLGY97Qr6pDVfVwW/474AlgNXApcEtrdgtwWVu+FLi1enYCpyQ5A7gY2FFVz1fVC8AOYOMoByNJOrYFzeknWQecA3wfOL2qDrVVPwJOb8urgef6NjvQanPVZzvO5iRTSaamp6cX0kVJ0jEMHPpJfha4E/hcVb3Uv66qCqhRdaqqtlbVZFVNTkxMjGq3ktR5A4V+kjfSC/xvVtWftvLhNm1Dez7S6geBtX2br2m1ueqSpDEZ5OqdADcDT1TVH/at2g7MXIGzCbi7r35lu4rnfODFNg10L3BRklPbB7gXtZokaUwG+WnlDwCfAB5LsqvVfge4AbgjyVXAs8DH27p7gEuAvcDLwKcAqur5JL8LPNjafamqnh/FICRJg5k39KvqL4HMsfrCWdoXcM0c+9oGbFtIByVJo+M3ciWpQwx9SeoQQ1+SOsTQl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDDH1J6hBDX5I6xNCXpA4x9CWpQwx9SeoQQ1+SOsTQl6QOGeR2iduSHEmyu6/2J0l2tcf+mTtqJVmX5B/61n2tb5v3JXksyd4kN7bbMEqSxmiQ2yV+A/gvwK0zhar6DzPLSb4CvNjXfl9VbZhlPzcBvw58n94tFTcC31lwjyVJizbvmX5VPQDMei/bdrb+ceC2Y+0jyRnAW6pqZ7ud4q3AZQvurSRpKMPO6X8QOFxVT/fVzkryV0m+l+SDrbYaONDX5kCrzSrJ5iRTSaamp6eH7KIkacawoX8FP32Wfwg4s6rOAX4T+FaStyx0p1W1taomq2pyYmJiyC5KkmYMMqc/qyRvAH4VeN9MrapeAV5pyw8l2Qe8CzgIrOnbfE2rSZLGaJgz/X8PPFlV/3/aJslEkpPa8juA9cAzVXUIeCnJ+e1zgCuBu4c4tiRpEQa5ZPM24P8A705yIMlVbdXlvP4D3A8Bj7ZLOP87cHVVzXwI/BvAHwN7gX145Y4kjd280ztVdcUc9U/OUrsTuHOO9lPAexfYP0nSCPmNXEnqEENfkjrE0JekDjH0JalDDH1J6hBDX5I6xNCXpA4x9CWpQwx9SeoQQ1+SOsTQl6QOWfRPK2tu67Z8e8mOvf+GjyzZsSUtf57pS1KHGPqS1CGGviR1yCA3UdmW5EiS3X21LyY5mGRXe1zSt+7zSfYmeSrJxX31ja22N8mW0Q9FkjSfQc70vwFsnKX+1ara0B73ACQ5m94dtd7TtvmvSU5qt1D8I+DDwNnAFa2tJGmMBrlz1gNJ1g24v0uB29sN0n+YZC9wXlu3t6qeAUhye2v7+MK7LElarGHm9K9N8mib/jm11VYDz/W1OdBqc9VnlWRzkqkkU9PT00N0UZLUb7GhfxPwTmADcAj4yqg6BFBVW6tqsqomJyYmRrlrSeq0RX05q6oOzywn+TrwZ+3lQWBtX9M1rcYx6pKkMVnUmX6SM/pefgyYubJnO3B5kpOTnAWsB34APAisT3JWkjfR+7B3++K7LUlajHnP9JPcBlwArEpyALgOuCDJBqCA/cBnAKpqT5I76H1A+ypwTVW91vZzLXAvcBKwrar2jHowkqRjG+TqnStmKd98jPbXA9fPUr8HuGdBvZMkjZTfyJWkDjH0JalDDH1J6hBDX5I6xNCXpA4x9CWpQwx9SeoQQ1+SOsTQl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDDH1J6hBDX5I6ZN7QT7ItyZEku/tqv5/kySSPJrkrySmtvi7JPyTZ1R5f69vmfUkeS7I3yY1JclxGJEma0yBn+t8ANh5V2wG8t6r+NfDXwOf71u2rqg3tcXVf/Sbg1+ndN3f9LPuUJB1n84Z+VT0APH9U7b6qerW93AmsOdY+2o3U31JVO6uqgFuByxbVY0nSoo1iTv/TwHf6Xp+V5K+SfC/JB1ttNXCgr82BVptVks1JppJMTU9Pj6CLkiQYMvSTfAF4FfhmKx0Czqyqc4DfBL6V5C0L3W9Vba2qyaqanJiYGKaLkqQ+b1jshkk+CfwKcGGbsqGqXgFeacsPJdkHvAs4yE9PAa1pNUnSGC3qTD/JRuC3gY9W1ct99YkkJ7Xld9D7wPaZqjoEvJTk/HbVzpXA3UP3XpK0IPOe6Se5DbgAWJXkAHAdvat1TgZ2tCsvd7YrdT4EfCnJPwI/Aa6uqpkPgX+D3pVA/4LeZwD9nwNIksZg3tCvqitmKd88R9s7gTvnWDcFvHdBvZMkjZTfyJWkDjH0JalDFn31jpandVu+vSTH3X/DR5bkuJIWxjN9SeoQQ1+SOsTQl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDDH1J6hBDX5I6xNCXpA4x9CWpQwYK/STbkhxJsruvdlqSHUmebs+ntnqS3Jhkb5JHk5zbt82m1v7pJJtGPxxJ0rEMeqb/DWDjUbUtwP1VtR64v70G+DC92ySuBzYDN0HvjwS9u269HzgPuG7mD4UkaTwGCv2qegB4/qjypcAtbfkW4LK++q3VsxM4JckZwMXAjqp6vqpeAHbw+j8kkqTjaJg5/dPbDc8BfgSc3pZXA8/1tTvQanPVXyfJ5iRTSaamp6eH6KIkqd9IPsitqgJqFPtq+9taVZNVNTkxMTGq3UpS5w0T+ofbtA3t+UirHwTW9rVb02pz1SVJYzJM6G8HZq7A2QTc3Ve/sl3Fcz7wYpsGuhe4KMmp7QPci1pNkjQmA90jN8ltwAXAqiQH6F2FcwNwR5KrgGeBj7fm9wCXAHuBl4FPAVTV80l+F3iwtftSVR394bAk6TgaKPSr6oo5Vl04S9sCrpljP9uAbQP3TpI0Un4jV5I6xNCXpA4x9CWpQwx9SeoQQ1+SOsTQl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDDH1J6pCBfnBNms+6Ld9esmPvv+EjS3ZsaaXxTF+SOsTQl6QOWXToJ3l3kl19j5eSfC7JF5Mc7Ktf0rfN55PsTfJUkotHMwRJ0qAWPadfVU8BGwCSnETvfrd30btT1ler6g/62yc5G7gceA/wNuC7Sd5VVa8ttg+SpIUZ1fTOhcC+qnr2GG0uBW6vqleq6of0bqd43oiOL0kawKhC/3Lgtr7X1yZ5NMm2dhN0gNXAc31tDrSaJGlMhg79JG8CPgr8t1a6CXgnvamfQ8BXFrHPzUmmkkxNT08P20VJUjOKM/0PAw9X1WGAqjpcVa9V1U+Ar/NPUzgHgbV9261ptdepqq1VNVlVkxMTEyPooiQJRhP6V9A3tZPkjL51HwN2t+XtwOVJTk5yFrAe+MEIji9JGtBQ38hN8mbgl4HP9JV/L8kGoID9M+uqak+SO4DHgVeBa7xyR5LGa6jQr6q/B956VO0Tx2h/PXD9MMeUJC2e38iVpA4x9CWpQwx9SeoQQ1+SOsTQl6QO8SYqWvGW6gYu3rxFK5Fn+pLUIYa+JHWIoS9JHWLoS1KHGPqS1CGGviR1iKEvSR1i6EtShxj6ktQhhr4kdcgoboy+P8ljSXYlmWq105LsSPJ0ez611ZPkxiR7kzya5Nxhjy9JGtyozvT/XVVtqKrJ9noLcH9VrQfub6+hdxP19e2xGbhpRMeXJA3geE3vXArc0pZvAS7rq99aPTuBU466kbok6TgaRegXcF+Sh5JsbrXTq+pQW/4RcHpbXg0817ftgVb7KUk2J5lKMjU9PT2CLkqSYDQ/rfxLVXUwyc8DO5I82b+yqipJLWSHVbUV2AowOTm5oG0lSXMb+ky/qg625yPAXcB5wOGZaZv2fKQ1Pwis7dt8TatJksZgqNBP8uYkPzezDFwE7Aa2A5tas03A3W15O3Blu4rnfODFvmkgSdJxNuz0zunAXUlm9vWtqvrzJA8CdyS5CngW+Hhrfw9wCbAXeBn41JDHlyQtwFChX1XPAP9mlvrfAhfOUi/gmmGOKUlaPL+RK0kd4o3RpUXyhuxaiTzTl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDDH1J6hBDX5I6xNCXpA4x9CWpQwx9SeoQQ1+SOsTQl6QOWXToJ1mb5C+SPJ5kT5LPtvoXkxxMsqs9Lunb5vNJ9iZ5KsnFoxiAJGlww/y08qvAb1XVw+2WiQ8l2dHWfbWq/qC/cZKzgcuB9wBvA76b5F1V9doQfZAkLcCiz/Sr6lBVPdyW/w54Alh9jE0uBW6vqleq6of0bpl43mKPL0lauJHcRCXJOuAc4PvAB4Brk1wJTNH7v4EX6P1B2Nm32QHm+CORZDOwGeDMM88cRRelE8ZS3bwFvIHLiWDoD3KT/CxwJ/C5qnoJuAl4J7ABOAR8ZaH7rKqtVTVZVZMTExPDdlGS1AwV+kneSC/wv1lVfwpQVYer6rWq+gnwdf5pCucgsLZv8zWtJkkak2Gu3glwM/BEVf1hX/2MvmYfA3a35e3A5UlOTnIWsB74wWKPL0lauGHm9D8AfAJ4LMmuVvsd4IokG4AC9gOfAaiqPUnuAB6nd+XPNV65I0njtejQr6q/BDLLqnuOsc31wPWLPaYkaTh+I1eSOsTQl6QOMfQlqUMMfUnqEENfkjrE0JekDjH0JalDDH1J6pCR/MqmpG5Yql/49Nc9R8czfUnqEENfkjrE0JekDjH0JalDDH1J6hBDX5I6xNCXpA4Z+3X6STYC/xk4Cfjjqrph3H2QtLIs1fcD4MT7jsBYQz/JScAfAb8MHAAeTLK9qh4fZz8kaVAn2hfSxj29cx6wt6qeqar/C9wOXDrmPkhSZ417emc18Fzf6wPA+49ulGQzsLm9/HGSpxZ5vFXA3yxy2+XmRBnLiTIOcCzL0YkyDvLlocby9rlWLMvf3qmqrcDWYfeTZKqqJkfQpSV3oozlRBkHOJbl6EQZBxy/sYx7eucgsLbv9ZpWkySNwbhD/0FgfZKzkrwJuBzYPuY+SFJnjXV6p6peTXItcC+9Sza3VdWe43jIoaeIlpETZSwnyjjAsSxHJ8o44DiNJVV1PPYrSVqG/EauJHWIoS9JHbLiQz/JZ5PsTrInyedmWZ8kNybZm+TRJOcuQTcHMsBYLkjyYpJd7fGflqCbs0qyLcmRJLv7aqcl2ZHk6fZ86hzbbmptnk6yaXy9nt2QY3mt7/1Z8osU5hjLr7X/xn6SZM5LApNsTPJU+7ezZTw9nrMvw4xjf5LH2nsyNZ4ez22Osfx+kidbRt2V5JQ5th3+PamqFfsA3gvsBn6G3ofS3wX+1VFtLgG+AwQ4H/j+Uvd7iLFcAPzZUvd1jv5/CDgX2N1X+z1gS1veAnx5lu1OA55pz6e25VNX4ljauh8v9XsxwFh+AXg38L+AyTm2OwnYB7wDeBPwCHD2ShtHa7cfWLXU78U8Y7kIeENb/vIc/1ZG8p6s9DP9X6AX4i9X1avA94BfParNpcCt1bMTOCXJGePu6AAGGcuyVVUPAM8fVb4UuKUt3wJcNsumFwM7qur5qnoB2AFsPF79HMQQY1l2ZhtLVT1RVfN9y31Z/WTKEONYduYYy33t3z3ATnrfYTraSN6TlR76u4EPJnlrkp+hd1a/9qg2s/30w+ox9W8hBhkLwL9N8kiS7yR5z3i7uGCnV9Whtvwj4PRZ2qyU92eQsQD88yRTSXYmuWw8XTsuVsr7MogC7kvyUPuJl+Xu0/RmJ442kvdkWf4Mw6Cq6okkXwbuA/4e2AW8tqSdWqQBx/Iw8Paq+nGSS4D/AawfZz8Xq6oqyQlxffA8Y3l7VR1M8g7gfyZ5rKr2jbN/ep1fau/JzwM7kjzZzraXnSRfAF4Fvnm8jrHSz/Spqpur6n1V9SHgBeCvj2qyYn76Yb6xVNVLVfXjtnwP8MYkq5agq4M6PDOV1p6PzNJmpbw/g4yFqjrYnp+hN9d8zrg6OGIr5X2ZV997cgS4i940ybKT5JPArwD/sdok/lFG8p6s+NBvf71Jcia9OfBvHdVkO3Blu4rnfODFvv9NX1bmG0uSf5kkbfk8eu/f3467nwuwHZi5GmcTcPcsbe4FLkpyarsi5qJWW27mHUsbw8lteRXwAWCl3ivihPjJlCRvTvJzM8v0/vvafeytxi+9m0v9NvDRqnp5jmajeU+W+pPsEXwS/r/p/cN6BLiw1a4Grm7LoXfjln3AYxzjU/6lfgwwlmuBPW39TuAXl7rPfX2/DTgE/CO9ucargLcC9wNP07sa6bTWdpLeXdNmtv00sLc9PrVSxwL8Yvtv7JH2fNUyHcvH2vIrwGHg3tb2bcA9fdteQu//NvcBX1iJ46B3pcsj7bFnqcdxjLHspTdfv6s9vna83hN/hkGSOmTFT+9IkgZn6EtShxj6ktQhhr4kdYihL0kdYuhLUocY+pLUIf8PKOvzg88aiqUAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "mstar = sub.MassType[:, 4]\n", "plt.hist(np.log10(mstar), range=[9.0, 12.0])" ] }, { "cell_type": "markdown", "id": "talented-pennsylvania", "metadata": {}, "source": [ "If we are only interested in the BCG (subhalo index 0), we can take a shortcut to only extract the properties of this one subhalo" ] }, { "cell_type": "code", "execution_count": 6, "id": "c2749e76", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Prepared reading from 'SUBHALO'...\n" ] } ], "source": [ "sub_bcg = hy.SplitFile(subfind_file, 'Subhalo', read_index=0)" ] }, { "cell_type": "code", "execution_count": 18, "id": "3b793eaa", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "The BCG has a stellar mass of 5.81e+12 M_Sun.\n" ] } ], "source": [ "print(f\"The BCG has a stellar mass of {sub_bcg.MassType[4]:.2e} M_Sun.\")" ] }, { "cell_type": "markdown", "id": "previous-boutique", "metadata": {}, "source": [ "Load stellar particles within 50 kpc from the BCG centre.\n", "Arguments here are: snapshot file (constructed above) -- particle type code (4: stars) -- centre of the region we are interested in (BCG centre of potential) -- radius of the selection sphere in Mpc -- strictly exclude particles outside the specified region" ] }, { "cell_type": "code", "execution_count": 25, "id": "earned-union", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Checking 1 cells...\n", "Reading 'Coordinates' took 1.859 sec.\n", "Region setup took 2.760 sec.\n", "Selection region contains 1 cells, 16 segments, 7155514 particles, 16 files\n", "Exact selection region contains 857773 particles.\n" ] } ], "source": [ "stars = hy.ReadRegion(snapshot_file, 4, sub_bcg.CentreOfPotential, 0.05, exact=True)" ] }, { "cell_type": "markdown", "id": "italic-hebrew", "metadata": {}, "source": [ "As an illustration, plot the stellar birth expansion factor vs. metallicity (absolute, not relative to Solar):" ] }, { "cell_type": "code", "execution_count": 29, "id": "controlled-accordance", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 29, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAD5CAYAAAAuneICAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAB8OklEQVR4nO29eXzcZ33g/37mHml0jyTrsmTZ8m3HV5wYQ8jhHARCAoEWFpayaWGXbXpl2wLbXehmlx/Llg3QJptC25RyQxLI6cSJEzuH4yS25VO2Zd2WNLpGI82hub/z/P6YwyN5JI1k67D9vF8vJTPfa57v+Dufz/N8TiGlRKFQKBTXHrqFHoBCoVAoFgalABQKheIaRSkAhUKhuEZRCkChUCiuUZQCUCgUimsUpQAUCoXiGsWQzUFCiLuAHwB64J+llP97wv6bgO8DG4HPSCmfSmy/Bfhe2qGrE/ufEUL8GPgw4E7s+6KU8thU47Db7bKuri6bISsUCoUiwZEjR5xSytKJ26dVAEIIPfAYcDvQAxwSQjwnpTyddth54IvAX6afK6XcB2xKXKcYaAVeSTvkr5LKIhvq6uo4fPhwtocrFAqFAhBCdGXans0KYDvQKqVsT1zoV8C9QEoBSCk7E/tiU1znU8BLUkp/lmNWKBQKxRySjQ+gCuhOe9+T2DZTPgP8csK2bwkhTgghvieEMGc6SQjxZSHEYSHE4aGhoVl8rEKhUCgyMS9OYCFEBbAB2JO2+evEfQLXA8XAVzOdK6X8kZRym5RyW2npRSYshUKhUMySbBRAL1CT9r46sW0m/B7wOyllJLlBStkn44SAfyVualIoFArFPJGNAjgENAghlgkhTMRNOc/N8HM+ywTzT2JVgBBCAPcBp2Z4TYVCoVBcAtMqACllFHiQuPnmDPAbKWWTEOJhIcTHAYQQ1wsheoBPAz8UQjQlzxdC1BFfQbwx4dI/F0KcBE4CduB/XYb7USgUCkWWiCupHPS2bdukCgNVKBSKmSGEOCKl3DZxu8oEViiuYgJhjcOdLgJhbaGHoliEKAWgUFyFJAV/Y9cIj+9vo8nhnv4kxTWHUgAKxVVGIKzx9JFuHn29FZB85eblrKssWOhhKRYhWdUCUigUVw5NDje7T/VTkmuizm6jqsi60ENSLFLUCkChuMpYV1nAyvI89jcP8k9vtSv7v2JSlAJQKK5CKvLN5FuMvN0yRGPXyEIPR7FIUQpAobjKaHK4Odju4o615eRZjMCVE+qtmF+UD0ChuMpYV1nAg7euoN5uo93pUw5gxaQoBaBQXGVYTXq21RUDUGwrXuDRKBYzygSkUFxluHxhfnawE5cvvNBDUSxylAJQKK4ynjnWy9/taeaZYzMt2qu41lAKQKG4yqgrsbKkwEJdiYr/V0yNUgAKxVVEIKxhNhj4xsfWsmN5qaoFpJgSpQAUiquIxq4RHtnbDAisJj1NDreqBaSYFKUAFIornHSnbyiq4QtGCUXjM/51lQWqFpBiUlQYqEJxhbP7pIPv7W0hrMVAgsWgp3N4TJl9FNOiFIBCcYVz94ZKwtEYES3G263DbF9WzJstTjZWFxKKaDzy6jkeun0lOxtKF3qoikWGMgEpFFc4xTYTG2sKOdA6zOryPL64cxkP3rIiYfYRICDxH4ViHEoBKBRXIBOje9ZVFnDH2nJaBn30uQNsqyvGatKzpbaIr39kDVtqixZ4xIrFiFIACsUVSGOXi2/vPkNjlwuIl3+4f2sND966YpzDN1kWwmrSL9RQFYuYrBSAEOIuIUSzEKJVCPG1DPtvEkI0CiGiQohPTdinCSGOJf6eS9u+TAjxXuKavxZCmC79dhSKa4W4accTjKYigAJhjRM9bg62DSkHsCIrpnUCCyH0wGPA7UAPcEgI8ZyU8nTaYeeBLwJ/meESASnlpgzbvwN8T0r5KyHEPwJ/CDw+s+ErFNceceEueWjXKpr7PTy2vw2AsbDG9149R3WRlf99vylVEE6hmIxsVgDbgVYpZbuUMgz8Crg3/QApZaeU8gQQy+ZDhRACuBV4KrHp34D7sh20QnEt0+Rw88SBTsxGHfdtruYvdjVwy+pyeobHEMCOZcUq7l+RFdkogCqgO+19T2JbtliEEIeFEO8KIe5LbCsBRqWU0VleU6G4ZklP7rKa9Cyz57KnqZ/2YT+/f30Nt6xZstBDVFwhzEceQK2UslcIUQ+8LoQ4CWSdly6E+DLwZYClS5fO0RAViiuLUETjYNswnU4fz5/oAyT3b6lmmT2XJw50km81KBOQYlqyUQC9QE3a++rEtqyQUvYm/t8uhNgPbAaeBgqFEIbEKmDSa0opfwT8CGDbtm2qt53imqexy8U3n2vC7Q9jNRlYU5HH3RsqybPo8QQj7Fxhp95uW+hhKq4AslEAh4AGIcQy4kL6M8C/y+biQogiwC+lDAkh7MBO4P9IKaUQYh/wKeI+hT8Anp3NDSgU1xqeYJQBTxBvSMPgjzDsC6HT6Wgd8IKAfKuRjdUFqhuYYlqm9QEkZugPAnuAM8BvpJRNQoiHhRAfBxBCXC+E6AE+DfxQCNGUOH0NcFgIcRzYB/zvtOihrwIPCSFaifsE/uVy3phCcbXSNxrAnwjzjElYWZ7H8e4RNCn51NYaHtq16iInsCoLrchEVj4AKeVuYPeEbd9Ie32IuBln4nnvABsmuWY78QgjhUIxA26oL6G8wEzfaIhCi57VFXm0Oce4Y90SvrCjLmPSV7Is9FduXq58A4oUqhicQnGF8dThbvpGQ9TZrRRajCwrsXHvpgpuWDa5YFdloRWZUApAobjC2LS0iIJjvfS6AnTEApzp81Bss7DvrBOzwcCta8ouOidZEkKhSEfVAlIorgDSbfgmgw69iNv/AYQQhKIarrEQnU7fwg5UcUWhVgAKxRVAug2/bzSANxglKsGea+Qj65aQl2PAG4xy5/qKhR6q4gpCKQCF4gog3YZfUWDl2WO9HO/14ByL8NwJB/Y8M4U5JvrcAaqKrAs9XMUVglIACsUVgicQ4ScHOwDoc/tT28fCGteX5PLvbqhTTl7FjFAKQKG4AmjscvHfnjnJgDeERQ8y0eFLD+RbjHQOjWE2CFX3XzEjlAJQKBYxgbBGk8ONJxglFI2hxWAsBiAxJTzBNosBs1GHavuomCkqCkihWMQknb997gBFOSb0aTI+rEnMRh29IwHqy/JU20fFjFErAIViEZN0/tbbbVQUWPnOy6fpHA6m9nvDMWqKLfz5rpXK/KOYMWoFoFAsUpLmn3WVBQQiGj96s3Wc8E8SjUpG/OEFGKHiSketABSKRUbK7h8I89i+Nh66fSUvnHBwtNsz7jiTDnLNem5eVabKPytmhVIACsUckT6Dz9Y8EwhrPH2kmxdO9oGUhLUYIIhq47ut5pl1WAw6rqsupHc0QLvTp8o/K2aMMgEpFHNEY5eLb714hl++15VVGeak8N/TNABS8m7HCFajDqcviMsfGfdj9YViuPxR2ob93LG2XMX/K2aFWgEoFHOGwB+O8vyJPjbWFE5bjK3J4Wbv2UHuXFeO2aDnbL+XQW+Yb798ln53aNyxEjDpBfdvruL+rTXKAayYFWoFoFDMEWsq8vn965fy5ZuWEYrEUquAyZqzrKss4MFbVnDL6nLeaR9ESsmHltvJM2eepwkhVei/4pJQCkChmAMCYY0nj3Tzi3c7OdA6zD+93U6Tww1ciO1Pvk8/50iXi288c5Jnj/bjDmq8dLqP88NjGT/DH4Fnj/VddB2FIluUCUihuMy4fGEe39/KK6f76BkN4g462FZnp6IgXqRtYnOWpLP4RPcoP3ithXA0hi5R7nnYF8FkEAgkcsLnrFli46MbK1UEkGLWqBWAQnGZefLIeX58oINBTwgtBm5/lP3Ng7xwopfDnS4AttUVp+z2TQ43j+5rxR0Is3ZJPrkGiCakvQRC0YuFP8CoP8K+s4O0qx4AilmiFIBCcRkJhDU6nX6iEgIJKR6OgVEv6Bj08ei+1nEmm0BYwxOIsqw4h+ePOTjd56EkzzqtaV8v4Pe21fDQ7StVBJBi1igTkEJxmQiENX75/nneah1EQGrWLgCTQc+pPi+f3FyVcghbTXoau0b4zstnMOl1rCjPI6hJBtyBjDP+dKxGHdfVFLGzoXRub0pxVZPVCkAIcZcQolkI0SqE+FqG/TcJIRqFEFEhxKfStm8SQhwUQjQJIU4IIX4/bd+PhRAdQohjib9Nl+WOFIoFosnh5tfvd9E7EiI9bUsP+EJRttcWsWpJPk8c6EhbBUj0OkEsJnGMBIhJjbHo9J9lNRloGfBklV+gUEzGtApACKEHHgM+AqwFPiuEWDvhsPPAF4FfTNjuB74gpVwH3AV8XwhRmLb/r6SUmxJ/x2Z1BwrFIiAQ1hjyBnH6QhfN3qNAMBLjrXYnIPncDUtTq4A1FQXYLAZO93vpHvHT74lk9Xm+YJQXjqsIIMWlkY0JaDvQKqVsBxBC/Aq4FzidPEBK2ZnYNy5fXUp5Lu21QwgxCJQCo5c6cIVisZDM4P3nt9oZ9k8+fe8fGeN/vXCaT19fw4HW4UQNf3B640le7lB2s3mbSRCNSawmXSqySKGYDdmYgKqA7rT3PYltM0IIsR0wAW1pm7+VMA19TwhhnuS8LwshDgshDg8NDc30YxWKOSeZwVuaZ5nyOG8YOob9+MNRHti5jHq7DU8gzCc3V1NXMvW5SQRgsxgx6HQc7xllz6m+y3AHimuVeYkCEkJUAD8F/oOUMrlK+DqwGrgeKAa+mulcKeWPpJTbpJTbSkuVw0uxeEhm9FYUWLmpwc7ta8oosk5dkkFqMfacHMBs1NHu9PGD11r52XudeEOxKc9LnQ+sLLNRmmcmz2qionB+VgCTZS8rrmyyUQC9QE3a++rEtqwQQuQDLwJ/I6V8N7ldStkn44SAfyVualIorhiSGb0vnOjlxwc6+eFb7YwEphaQUQnlhWaOdMUVx8bKfIa8YZy+7Ov563WCz2xfSk2RlXyL8VJvIysmy15WXNlk4wM4BDQIIZYRF/yfAf5dNhcXQpiA3wE/kVI+NWFfhZSyTwghgPuAUzMZuEKx0KyrLOCBnXW8cKKPvtEAxizqsZkM0Nzv5VDHCEa9jjyrEb0Abbq4zwR6oN8dYn1lPltri+YtB2Bi9rLi6mDaFYCUMgo8COwBzgC/kVI2CSEeFkJ8HEAIcb0Qogf4NPBDIURT4vTfA24Cvpgh3PPnQoiTwEnADvyvy3ljiiuXK8HckCzfEIpKDrY7iQFT+H8vnBeFmJTcsbYMpzfIT9/tIpyd9QcADagpySUUlTPqM3CpWE36cdnLiqsDIWWWU49FwLZt2+Thw4cXehiKOeZwp4vH97fxlZuXT1tCeb5JCf6Ixj++0Y7NYuBIlwunN0y2clwAZTYjIU3iC0aJSbI+tzxPT4ktBy0m+cbH1rGzwT7LO1FcSwghjkgpt03crkpBKBYdi9nckLSFg2B5WS5vtwxh0MmsBTjEHbkxBP5QlFKbCf0MfoUGnZHqQitaTOIJRhb9SkmxuFEKQKGYAXG7/zI8wQiOET+BsIbDnV3yFoARqC/JweULgxD0ecNEZqA91lXmM+AJodcJ+kYDUzpmrwRTmmJhUbWAFIuO5Cx7MZmAAmGNxq4RQBKKxvjWi6fpGQ3O+DpRYNgXLxWhZev5TZBv0bOmMp9PbMnHbNABklVL8iddKWX6HmfTp1hx9aJWAIpFRbw6ZpidK+yLqs59Y9cIf/t8E9989hSvnu4nomnMJgBTEs/4nannTQfIqMYTb3cw7AuRbzXy8/e6MRt1kwryTKY0Fc6pSEetABSLiiaHm+/vbcXlDxPRNL6wY9k4AbdQM9hQVCMQjhKNxXijxUlY08je8HPpxICxKOxaaafEZsYTiPLAzrop/STJyJ100pWCWg0o1ApAsaioKLCSY9Iz7A3y1JGei2aqCzGDDYQ1Wga9GHSC0lwz/lAUn38mbt/LQ55Zz8bqAn70RhuP7W/BbNTPWHCnh3Oq1YBCrQAUi4Kkjf2Zo90c7R7hA/Ul3Le5miOdLjyBKDuWl2A16RckQuhg2zD/9k4nvlCUTldg3j43SZ4B/DGwWY28dmaQezZVsqp8ctt/tszku1SrhasTtQJQLAqaHG4e2dvM3rODhKKScFRj96k+vr/3HN9+6UxqlroQCUmdTh8jYxFisYXJmbHnW7EadCzJt3DbmjLu21TNzgb7JX8HM/ku1Wrh6kStABSLgnq7jS01hXQOxfvbHu92E47FsBh03HtdxTibdb3dRrvTN2+z0fs2V9Pt8vOL97rm/LMyMeQNsaoin+YBD55ghEKrkfu31syrElzMuRmK2aMUgGJRcKbPwy/eP89YOIYAcswGrECB2cDmpUUA/PK9Lp4/0cc911VwoHV4XsJEA2GNY90jGPSCLAt2prCZYAY13i5CJ6DIamQ0EGFluY071i2hZ2SMPacHWF2RP68hspkcyoorH6UAFAuOyxfm5SYHaHEJKwGdDv7L7aupLLSwpbaYJoeb50/04QtFiGixaSNgLhf7mwf52tPHcQdnnkx1KcIfAAnD/gglVgMGETdFnezx8MktVfM2E1e2/6sb5QNQTMtcZ5TuPungt40OxtIu3+8Js6fJQbxyTtwE8Zd3rOSz22s52O7KOgLmUsd+oNU5K+F/ObAa4vc+HIjy0/d7eOpID4PeAHV227wJ46vd9n+tZ0srBaCYlrkWApuXFiXE/Hje73DxF78+yv7mQawmPWsqCohENXbUl1BRYM3qhzvbsQfCGq+fGSQYjsx4mayDjPczU0JpmcJGATkmAzazEbPhclw9O6522//VruCmQ5mAFNMy10LgZ+924csgyAORGN5QmN819mA26Hnr3CBPHe2hssBKrknPvuahjGUOKgqs7Gnqp67Eyqaa4lmNvcnh5n+92ETXsJ+Zzg1l4m+2mHRxm3soolGWa6bfE6K6yEogovHJLVVsqb1gi58LE83Ea17Ntv+rXcFNh1oBKKZlrkIvXb4wP3yjjfc6nBn3yxgsLbbS7fLz3545we6TDlaU2rhtdSklNjMP7FyWsczBP73VzndfOcs3nj3NmT73rMa+rrKAsjxz1o1axo17wvuZzrJMBh2+kIbVbOC2NaXcvq6Mfm+I0UCUHLNx3L3MxQz2WpoVpyu4a9EUpBSAYsH4t3c6+M5LZ2l3XpxcZdGDxaQjEI5SVWRBAiOBCKcdbp480ssjrzYDjBOGydnc5poCck0G9HpJ84Avqx91ui04ENY42DZE57Bv1vdmSnudRZ+YcfjCMaSEEX+Un7/Xw4bKAv745hXct6mCklwTB1qGUvc03Qx2Njbuide8Fuzk15LSS0cpgGucQFjjQMsQB1qc8/4DP9EzOmkdfaOAYCSG0xfhUOcIn9hUxZ/e1sAnt9SQZzEQt7KPn2snZ3Mmgx6TQUdDWT57Tw/w9JHurHwFj+5r5ekj3RxsG+K/PnOSfs/sq/2EgZxLWDAlvxcBLC/L48blJbQMjfH/9rfyyKvnsk6Mm41gmzgrbuxyXfXC8Vo1BSkfwDWMyxfm8f2tvNfhwmgQfP0ja+Ylrr6xa4RQVMPtD120XwCFFh2FOWaWl+fxfscw3pCGP6zxJ7etpLHLxe1rl2A2iHG28PTrv9c+zMhYmNMON9fVFPHSqX6AKZOn6u02Gkpt7Dk9wAdX2Bmbprl7Nvgv8RICsOcbeKdtiA8st/PHNzcAErNBN05QTeUHSPYtDkViBMLajExhSeXxwM5lV71wnOjruFbCX9UK4BrmF+938pN3u1hVlsNDu1ZdtOR3+cJZL/2zNRMcbBvm6789wf/34mkae7wX7ZeA2Wjgw6vK+B8fX8937t/IF3bU8u931PH0kW7+6e0O8q0GdjaUZvxhNnaNcKRrhBWlufiCUTqHfDSU5/F0Y0+inn9mzvS5ebfVic2sR0iJP7Jw5o5Cs2Dn8hLMBkFxroXfHXXwrRfP0un0sWO5/aJ7n2qWbzXpMRv1PHGgY8Yz+OSseEtt0TXXD/haMQkpBXAN09zvJRSNEdTkuNoyyYd/90lH1j+Cg21O/uZ3JznYltmhm6Rl0MuQL4TZoCPXeHE4o0HAjuXFHO0eZU9TPzevKuerd63hjeZBfvV+NxX5Zt5td/L6mYGMyiYUjTLoDXLS4cUT0ojJGN3DY0S0qWNzHKNBzg152XNqgJ+/d35Wzt/LxWhIokU1LEYdvkCEDyy3YzHqeOaYI6M5azrzxUxt+sn9wDUn+JNcKyahrBSAEOIuIUSzEKJVCPG1DPtvEkI0CiGiQohPTdj3B0KIlsTfH6Rt3yqEOJm45t8LIeYvuFkBwLutAwC80dw/ThgkH/67N1Rm/SNoGfTSMxLgtMM9pXAx6XSEozGa+n2MRS6WslEJBp2OYETj+eO9NDncHGxz8i8HOuj3BHj1zCD/b38733n5LE0O90XO285hP95gNCXqmwf9vNY8RE1xzqQmo5dO9vH9V88S0uK29wH3/Ff8nMigL8RYWEOTgs9cX8P/+Ph67ttUyd6zgzR2ucZ9x5P5ASYT5NPNbpscbh59vTUr38nVykIUHVwIplUAQgg98BjwEWAt8FkhxNoJh50Hvgj8YsK5xcA3gRuA7cA3hRBFid2PA18CGhJ/d836LhSzYsgfF5OeUNx0khQqyYe/2GbK+kewJN+Kxajj1dP9/I9nm2jscl10TNugj++9enba2bWUkq/etZq/vGM16yoL4qsGb1wg5hh1fHJzJV+9K74vXZg1OdzsOTWAZUKilEEHFoMuozA72DbMV586Tq/nQt2G4PyX+k9RUWhghT0H91gATYMV9hwAttQW8dkbannwlhWAyGplNpmgn2p2GwhrhCIxbmqws/fs4FVvArnWyWYFsB1olVK2SynDwK+Ae9MPkFJ2SilPwEVBHXcCr0opXVLKEeBV4C4hRAWQL6V8V0opgZ8A913ivSgugVA0OmubZyCsse9sP8NjEY73eun3BknmwqbP0P/Hc6fwhKeWriU5Bs4N+sm3GlNmqYayPAqsRpCSW1eX89Dtq8m3xhsyJoVZvd2GJxBFCMnIBAduNAYvN/Wz+6Tjos873OXEE7pwfJ55Ya2iI54oPaNBfFEdCGgZ8vHY/tZx/y5rKvKzWplNJuinmt02Odw8caCDVUvyefCWFRedm62vZzaho9dCuOliI5unvQroTnvfk9iWDZOdW5V4Pe01hRBfFkIcFkIcHhoayvJjFTOlczgw6wJrTQ43bzRf+Le5b1MFW2qLUvuSimXA45/yOiW5BjZWF3HLKjueQJRAWMPlC9M57Oev71zF729fypduWs6ZPjff3n2Gxi4XVpOeeruNx/e38oO9zbjGwlgzNOu15xoxT1gF9I4EeP54/7jjvDMt+XmZCcagpsSKFo2BhBvri1MO+uR32e70TSrA04XoTMwYyfPq7TYe2FkHkDECJlvn6GycqPPteFUK5wpwAkspfySl3Cal3FZaWrrQw7lq2X2yD7Mx/mOf6Y+iosCKzXJB6h47f+EHnD4L9U8z+9cJHe+2O/nloW7+4fV4rPvukw4e3dfC3rMDvN3i5IUTvYSiMrHAELh8Yb7+9HF+9l4X9WW5RLQY/gzh+92jIf7PnmZ+crCDAy1DuHxh/vdLZ+geWXh7/0Rubijlw6vLMBh0LC/LT62EsnFMzlaIpiuXqaKGJvYUnuxZyfa4yc6ZD66VSJ+pyEYB9AI1ae+rE9uyYbJzexOvZ3NNxRwQCEcvsqlPffyFH/W+swP0jgZT+46cH035ANLjq7tHL477T1JgEjx873ru31pDeb6Fe66rot5uo7LQwp1ry3m3dYh251iiKUuMT2yqos6eyz+8do7XmwcJR2Kc6/PSMZxZoOebdQhi/PpQNw+/cJpnjvXQ5w6SIRBpwXGMBCjOMVKaa2J5qW1GUTmzFaLp5011jfR/z6ePdPO9V8/xyCvNuBK1rzONNdtnai4dr5mU0LUS6TMV2SiAQ0CDEGKZEMIEfAZ4Lsvr7wHuEEIUJZy/dwB7pJR9gEcIcWMi+ucLwLOzGL/iMjEwOpbqtjXdjyIQ1vjJwQ6+/tsTHGxzsnlpEcW2CyuAGOAJji+AMN2PvyDXzEc2VPA3H13LN+9Zx2e3L6Xd6ePn73VjNhoYC0liEurtuYCOJw/38NdPHuNX73ehS3zmucHJSzd4QzG8wSi5iSJr/rDG+sp8MgQiLTgvNQ3wVGMvWixGvsWQlQC91NDNdOGbjSBucrjZ0xSPInu6sSflX8k01nQ/zUKZXDKN61qJ9JmKaRWAlDIKPEhcmJ8BfiOlbBJCPCyE+DiAEOJ6IUQP8Gngh0KIpsS5LuB/Elcih4CHE9sA/jPwz0Ar0Aa8dFnvTDEz9Ppp7ctJ9jcP8tjrrXQ4xzjeM8Kj+1pw+i7YXXTEI37ggmAqyjFNcrU4W5cWXCTAklms6yryqS6yANDvDvLSyV4GfAEOtLsIaBCKm8uJTGFhkkAgCqf7vDh9YV483svw2OQrksvJTO2sMeLjTYauzqXpZ7asqyzgznXlANy/pZq7N1Smtk8ca1LQtjt9M15dXs7xXuuz/UxkVQpCSrkb2D1h2zfSXh9ivEkn/bgngCcybD8MrJ/JYBVzh9cfyfoH8trpPjwhDQG83TLMiZ7RcftjQCQW//EmBVNxTgbPbBrdoyG+vfsMD92+kp0NpamCbPuaB9lzqo9BX3xFcWZgjJahMbbVFjHkCc84YSvXrCeiSRzuIGcHxmZ28iyZiVtZB5h1sKIij8/fWAtk144xm6JwMy1tMPGcie+TQv/uDZUU20zTjjVbIXyhBEUdZqP+spRjuNrLWs+WRe8EVswPMpa96UCTccO5BBrPj2aceRt18Uer3m5j54oSTvVOPuvTA5/YVA0CQlHJ4U4XrzT18xe/Oc4v3+tOCf8k0Ri0DYzOKlvXE9TQYjHcQe2SavbPFTHAYNCzqszGj9/pyphPkYmpzBmBsMbTR7p5dF/rJUXlpL8PhDWeOdbLk0d6eOZYT1az9WxNLklFMTHfQUXtXH6UAlAAkDeD0pUWwwXRmUmICqA830ogrLH7pINnjjpomcI+L4CKQjOf2FQFxHh8fxsvnuzDHYgSnURKJ5PYZsNilx/esMbbbcMM+0LxiKdpmE4wNjnc7D0zyK7VZdPOvtOvNXHGnv6+yeHm2aM9OEYD/Py9rinrLM107ElFsaW2aNznN3aN8O2Xzsz6sxYjC63UlAJQABCLZWeo6B0J8MKJgSmPkcD7HU6+9WITT7zdwdoK26SCHGBpSQ6gY1/zEGaDga/cvJwdy4qv2VK1hVYDq5cUUJRrytj+caLQmM7+v66ygAdvXTFlNdQkjV2ucTkW6TP29Pf1dhtrK/OJSYg/Otkp5JmM/eIVg0x8zGJcu82OhQ5FVQpgEbGQswGDPjtx+92XT+OdJp4f4FDnCL8+1E33aIDyfOuUxw54/PSN+vnKzctZU5FPKBJj9ynHjBupXA2YdKDTCbqGx7j3usqM9YsmCo3pbOvZmF7ifSGcnHK4CUdjvHpmIBXamdyf/my2O30c7hrBH46yuiKfNRXZOVcn1hmaSUObLbXFfP3uNRm/kyuVhXZOKwWwiFjI2cDystyLtmVSSOcGs3OcdroCICHPrMcxMnVnrbEIPH/CwbrKAp48cp4//9VhDnVdm8k5kRiMhaJIJKuW5GcU2nMhNJocbr77SjNPHemhutjKc8cd40pnpJtfkvWCHrp9JZ/aWoVjNMCx7tFxz8pkk5l1lQXsWlOWqjOUSTkFwhqvnxnkibfaOdg2nPpNXI1hmwt9T0oBLBJcvjBHukbYvqyIertt3j//dO8or5/pH9cZLJNCimnZrU6sBtBkPLv35aapS3jYcwz81Z2raexy8dOD5xkaW+RG+jlEEneEL8m3sKYif9y+2Vb3zIZ1lQXcc10FOQY96yryefDmFakon9TIEuaXZL2gfIuJZXYbSHj9bD/fevGCfX6qMVUWWripoTTjc550WH/7pdN8/7UWOp0+Fb45h1yrZtZFxzPHevn+3mZsZiMNZXncuqZ8Xj/fE4bvvNyMzWJIdQbLNNNsH5y6nk+SQBRyTToKc4y0DoWnPFYKgTsQ5rTDSyh6LRp+xlOeZ8ITiHKmz83OhgvlT9I7dJmNulSZhRPdo3zuhtqshWQynLPebqPd6UuFWX52ey0mvY5XTg9wx9rycbPSLbXFPHT7KkCmkgU9gQjPHnNQVWjheNcoYXmh58Jkq5Qmh5vH9rWBgFXltovCPJscbl45PcCN9SVUF1q5b3N1KsR0qnuZLFRVMTVqBbBIKMk1ohOCUX+Yfc2Zm53MJTaT4Kt3rRrXGSyZjJUM+2sb9JFt6pQOKLQa+Q8fWDblcSY95FuMvNkyzG+P9jLsm30f3quB8jwDFqOeSCxGsqJqkgvhkTI1u9590sFj+9twjPpTAi+bqKBMDX+sJj33b63hjrXl7D0zvu9AvLOYjicOdKaSBc2GeN+GbleArfXFfPWu1Sn7/GSmjXWVBTx0+0oe2rWKTGWt6+02VpXn0TnsZ2uiJPlUTBWqqpgepQAWCd0jAcbCMSIxeOPc0Lw/wCad5NY1S8Z1BoMLP6jGrhH+/JdHs75eDOh1h/jXA+1THreqPI/b15Tz8slejPqJIu/aQgAmvQEtJrl7fQXN/Z5xjtgL4ZHFqdn13Rsq+YtdDePMNVMJwbj9XuOBncsyNvxJKoEHb72470B6SYcDLU5CUY2v3rWGT2yton1oDMfo9IX1rCY9OxtKE9ViJZ+7oZZQRBvnXG4Z9HHn2vKsSlFnClV9YOcyPIHwOHNmpu9B5RQoBbBo8IcumD4G3KF59wPEZGbRm/yBNfWOcrLPM+PrTldt80yflxdPOej3RmgeGONa/jlKQKeDAqsRTyDCY4lZOkxe5rnYZuLzO+rGzZTr7TZuWVVKUY6Jnx3sHKdE4vb7TsxG3aQNfyaLw09uP9Pn5m+fP8W3dp+hZdBDXUkuN68qZe+Z7BvIJMfhGPXzxIHOcUpmspDViYotk7knuVJ5bF8bj+xtnrLrmVopKB/AoiAQ1sZlfIY0yVstQ9y7Odu2CzP/vIkPfkiTKeGSTvJH/xe/PjarzwpPk66bY9Yz4JmfmjyLnVyT4LwryG2rbXx+Rx0GvY4dy+0caHHS3O/hzRYnD966YtqSBmf6PPzuWC+tQz6eP95HWJNsrC64qNLnhedOsKW2KKMNfVtdcUr5JLeFojHGwhrhiMZvDveQbzFyz3UVfOlD9Vn7IdJXE6sr8rM6b+JsPynEv3Lz8nHfSdLMBCLrPsnXKmoFsAg42ObkUNfouG27T13cvSobslnaJmOx0wlGp67Yuar84jDRbJguZScQ1ha8C9diYSwskUBpnoURf5h25xgH25w8sreZ50/0sWvN+Ezeyf+t4xE7Nywr5o9vWUEkqvGD11p4+ki8N1N6meZHXj03bqacaWZ8sG2Yv/ndSQ62Daf6LluNegpzTPzetmruua6CN885MRt1WTteJ7YdBTjQMsS/vN2eqjI78T7Tx55uysrU8WxnQ+lF5sxMn3+pjuIr3ZSkVgALjMsX5hfvdV1U12bX6iWzut5ks6J01lUWcF11AfvPXQjPNAgmNTsFwhrvtA/PajzZMBJY2C5ciwGzHiIaNJTn8vkb6xjwBNm5ooRbVpfHQy2RbKktHufoffpIN6+cHuDmlaWsWpKfmsVvqS3mj29ZTudwgJJcI4/uO8+NdUXsPTvI6or81HORPlNOlmrOVA680+mjzx2k0+nDbNDx7LFePrC8mA8st2M2GAhFo9y00k693ZZYVYyMG282kTlJZTTgCTLij9A5PDZu38TicEkT0lduXp51x7O5iA7K5ve2mFEKYIF55lgP+5svjpO3WWb3kCbtv8kfY6aH3mrS85sjnePOk8QdcMW2ix/ig23D+MNzk34/VQnnKxnBzAoWCCDHJHAHo7zd6uSV0/0gYWN1ITsb7Bcdn6zvs6o8j+dP9PHiqb5U+K7VpMcxGuSxfa1srS3CMRqg2LaEVRX545R8Ulk0Odyc6fPwxIGOiwRZIKxRZ8/lj29ZTp3dRigaJRjRONXrYZk9j18fakGLSfQ6waryfMxGHY/sbQYJX787Pp6JQjL9uUzeS73dxkO3r8QTjNLnDnDfpgvFhZPmmlAklrrOTE04mQS1yxdm90nHuGqmM+VKNyUpBbDA5JkNGevknHfNrFVh8kcVimjsax5idSKJ6AevtVBbbKWmOIdPb12aetD73OPj7W9MxP1n4q2WqWv/KC5mpuoyqIEFcPlCuP1hPrqhkroSa+rfZKIyTzpL6+02zvS5Sdq7k8fdsjqeR2I26DmbcN7/7lgvy+y2lEJx+cI8vr+Vs/1e/tOH6yeN2//5e+e5ZVUpP3+viwd2LuOrd62mc3iMklwjUU2ytMhK90iAUDTKltrSRIinzFhELnnNpDAORTS+s6eZrTWF/MltKzMK4nRfBEAoEv//TGbcmQT17pMOvre3BYDP76jL+lqZxnalohTAAvNWizPj9qXFF+rnZLuEfnRfKzvqS9i5ooR6u40Rf5jRsTCHO1zodIJck2HSB73XHeCX75/nznVL6HMHxn3W/rNTZ/IqLg+VhWb6PWG8oSjNA042Vl8wb0ycwaYLxaRZBOJtGveeHeRLH6xndUU+nkCYEpuZQW8QTyBCKHrBVr37pIOnG3vYucLOmoqCjMI3k7O2yeHmQOsw25cVE9FidLn8GA06zAZDwv4eVzDpzuOJTtqkMD7YNky/O8hTQz3UFOewakk+E81d6XQ4ffH7gxn1CsgkqNP7GVwqV2oCmvK+LTBtg96M28PRC7aRbELW6u02GkptvHDMwa8PdXOse4Tv7z3HSYcHTcb4+IYlUz7og54AP9h7jh8f6Ljos8LyynRwXUmY9YJdayr4/W01fPEDyy6arU5mamjsGuF/Pt/Ed14+y8G2oVTZ52SymNlg4P4tVZzq9RCMJHuNxbl7QyX3b6lm1B+h3Tl1vaZ0hZPMPF5abMVs1PPJLVV842PrErH948eWqXxzugPWbBCU55n5YEMpFYVWvvvKWR5+4XTGPgjpZa1D0ViqaulsyRRCO1uu1LBStQJYYEKRzKUP9p4Z4N7NcTtoNhUTf/FeF7tP9bG9rphTDg8tg16GfCEMIl7/ft+5Ib796ckf9O31dmrtNj6zfSkj/vC4zxpyT13KQXHpbKgq4MWTDkpsFj56Xfii2erEGWzK5BeNMuQL8av3u7Dnmnjw1hVp/3ZxT8TdGyoZC0d56kgP6XM+q0nPzavKuHlVaeqcpF38ltXldDrHaO73sPfsIOsq8vnKzStSmcd/fPPyRLYyrK0syOinmFi+OZPNfU1FAdtqizjcNcL54TFuW13Os8d7CUVjk5q91lUWxJWKSN7jwnOl+gKUAlhgOoaDGbefmdC8Orn0zrTEbOwa4V/f6WB4LMKwLwhI/CGNinxTyr8QCEczxvkncbiDjAYi3L6m/KIfcxbVnxWXyLkBLxEpuWPtkqyESJPDzQ9ea2F5qY2dy0t48dQAEU3jSJeLI10jfHprDSD57p5mrq8r5q3WIQxCjOsv0Ng1wiN7m3lo16pUtM7j+1t5+mgPXcN+GrtHiMWgusjK00d7qC3JSa0iKwtz+IfXWwBJy6AXs0HHmooCzvS5CUUlZoNgTcWFKKNkc6B0m3sgrPHPb7fyq8PdWAyC3xzu5veuryHfasRs0NPY5eK7e5q5Z1Mln91eO8EXIHlo16qLVh1TMZdmmivVF6AUwAIzmXFlxB8ZJ7CnimIosZkw6QQ6YNQfISbh6aPdOEYvJFhFtBgH24a4dU3m8FJvKIzTJzna7UqFEwJXbHzzlYIALHr41JZqqotzuG9zdVbCaV1lAavL83j6aA//8abl/NeaIsbCUf7+tXh+h1EnQAh84Sg/f/88Y2GNpUXWCXX7x8/QmxxuzvZ7uX9zNZ/ZvhR7XnyWvmvNEm5YVszdGyqxmvSsrsinosDKPRsr8EfiK4vdx/vYXFvEgVYniHh9p6/fvQYQKSUz0ebe2OXiJ++cJxCJEY4ChCjOiXeGq7PnsqepH184yvPH+9hYXZgS/k8f6WbP6QHuXFs+IwUwlyGbV7UPQAhxlxCiWQjRKoT4Wob9ZiHErxP73xNC1CW2f04IcSztLyaE2JTYtz9xzeS+sst5Y1c6esE42+nEDM7DnS6eOdbDI3vP8V6Hi09srcZq1NE84KXAaiAalcTSQlGCUTjefSGFPh0DEAhp+IJRnj3mGGfHfKWpf07v81pHAiENWobGqLPnzCiR6os7l3H/5mo+trGS+7fW0FCWzx99qI5PbqmixGbm+eMO7t5QQV1xDgDBcCQRMRRnYoOVdZUF/NmuBh66YxUj/jCvnB7glaYBRvxhPr+jDqtJn+otvO/sAK+dHeSN5iGimmRtVT5vnBtkxB/i7vUVPHT7StZVFhCKxp+rUFRLKY9ARONnBzvxBCPUFFkosOjRJIwGNBrPD7P37CA/PtDB/uZBPrWlmnuuq0yFrzY53Ow5PUCh1cgrpwdmZHOfSzPNleoDmFYBCCH0wGPAR4C1wGeFEGsnHPaHwIiUcgXwPeA7AFLKn0spN0kpNwH/HuiQUh5LO+9zyf1SysFLvpsrgKTwdvnCHGiZPLpmNBAl3WGXbgZq7Brh8f1t1JXk8tENS9jfPMjBVieRmAQhaDzvpt873m4fAxwj8eSaiQ9plPiPrzzfxC0ry/AEwikl8cv3uy7LfSsmJ0Y8G/w7L09euyYT8fIQQzT3e2hyuPnJu534QzG6hv0c7Y5PHhyjAc72e7HoBRoCT/CCzyndGZvJ3v7HN6/gjnXleALR1P69ZwbZUV/CqD+M0xvgUNco511+XGMRhrwhnGNR9p8dZE1FQcLJq8NmNgCSnxzs4Ou/Pcmjr7fwf14+y5vnnHz0ukpWlsVDlvUCqgpz2LW6jCaHh1XleTSU53Og1ZlyUq+rLODOteWM+iPckaFgXJL5ztC9Un0A2awAtgOtUsp2KWUY+BVw74Rj7gX+LfH6KeA2IcRE78xnE+de0yR7rj5zrIf/74XTkx8YgyUF8VDQZLu+nxzs4FsvnCYUjfLAzjqGffGKh53DAY72eIjGJNEpmu9GpUi14ZuIBM6PhPjl+138w+sXZjLtgzMvAKeYGQYBH1hewn9OJDtlI7QCYY03W4biGbrDY/FOW6vLONvvpb40l/fbXZTaTDhGgmjE20wGIjH63IFxwjH5urHLdVFpaLNBx1NHevn2S6dp7BpJOWEB/v61FpoTvSE0CW+eG8CXcBYdd3h4fH/ruDaOZoOBp470JCqGxijIMfFu+zD//FYnR87HlVWtPYfyfCtjYY0Pr7LTMuTDbBAXFaS7f2sNf7arYcoex5lm5HM5S1/ozl6zJRsFUAV0p73vSWzLeIyUMgq4gZIJx/w+8MsJ2/41Yf757xkUBgBCiC8LIQ4LIQ4PDV0N8eiCaEzSMuBhJDh5EbQw8ON34qWUD7Y5+dpvT/DrQ934IxpmgwFPMMo3nj1F+/CFhLGYnNynAMRXCm3TlXQQ3LOxIlUaYPAa7s41XxiAntEgw2MRnjjQMa2ACoQ1fvl+F4faXexYXkxFgQUgJRg31xThcAfY1zzImT43eWY91cVWVpXncsOyknGN35P5I80DXh7YuSz17x4Ia4SikkA4SiwGoWiUJoebigIrnU4f2oTAgNXleSlhYgTO9ntTbRyTpqB7N1XyRx9cRlGOmQ8uL2Z9VT6+UIRSmxE94A1E+LtXmnl0XwtGnY4v3FhHc7+Xervtokz29JXLZK0nsw2lXWzM5+plXvIAhBA3AH4p5am0zZ+TUm4APpT4+/eZzpVS/khKuU1Kua20tDTTIYue5Az+9TMDhKJRNlTm89vGXhzuqZufFOeaAWgZ9OL0Bskz6dlYnY8nGOGd1iGC0ZmF53iCUZ480s3+5skze+tLcrhzfQVn+jx868UpViiKy0ZQgmssREmuMSsB1eRw8/wxB+5ghKZeD4/vb081WQfo9wQIRWLodQK9DqKa5NzAGMe6PTx5pAcQqRDKdZUF3NRQyvPH48UHj3WP8De/O8n+5gFOO0aJSfjE5krMBgOP72/jxwc6eO5Y70UTjRO9XpJPo9kk+LNdDan7aOxy8Z2Xz/La6UH8YY1/O9jFcyf6CUdj6AR4Q1EMBsHGqgLsuSY+tNzOnesrcIz6UyWxJwrFyVYuSTLNyK+UWfp8+hOyiQLqBWrS3lcntmU6pkcIYQAKgPSp5meYMPuXUvYm/u8VQvyCuKnpJzMa/RVCY9cI33zuFO5AhOJcI5FojFAWyr26MO68s5mMaBLODHhpHvTxbvsIf3nHKhq7RmjqnzqBJx17joF3Wodw+iZfeTT1efnHN1r5wPIShrwq/n8+KLQItJhkeCycyv2YinWVBfzlnas42j3Kb4/0sKTAwkun+uka9nO6z0O9PYedK0p4p9UJOkEgGsNIfHVYmW+mzp7Llpoi6uy5WE16KgrMOH1hPMEwfe4g/e4g73W4ONA6jDsQASFoHvDwuRuWUluSy6un+xhzjQ9fTn+cP7SydFyUTSgqiWgxyuw53Le5CnueGZBsqy3mvMtP+9AYH2qwc//WatoGx/jdsR5eOOFg15pyWod8lNhMHGwb5rH9LTy0axU7G+zj2mNOpzSnitBZjNE787lSyWYFcAhoEEIsE0KYiAvz5yYc8xzwB4nXnwJel1JKACGEDvg90uz/QgiDEMKeeG0EPgac4irF6QvR7w7i8YcJR2PjzDZT8Yv348089p7pJ6xJDECuSY/JIBj2hfBHZrZEHB6LUllgpal3dNJjosSLv53s8TDozZyjoLg8COKOz4JcEzevKuPOdRVZLf+tJj1rKgoY8AQJRmO09HuwWQy82+Yk16znt409BCJRcs0GhrzxVWaEuLP5paZ+fnygg6eP9rDvbHwleN7lZ3gsxJvnnFQUWPniB+oIhqN8+UP1/JfbV2LUCx7b14ZjNMC5AQ/dacK/0Hyx5XYsqCVWvUMcaHECknBU8laLk/c6XHxhRx3rKwt5prGHjqF4hdEBd4C/fa6Jvaf7cYwE+dnBTvae6edg2zA/erOdTqdvXMhqUkhuqS2adlY/1Yw63Ry2WJjPlcq0KwApZVQI8SCwB9ADT0gpm4QQDwOHpZTPAf8C/FQI0Qq4iCuJJDcB3VLK9N6AZmBPQvjrgb3AP12WO1qE7G8exBfWKM01IET8hrMR3Ue7PfzdnrPxWRjg18Dvj+IORvm/L5/JahWRThToGRljbUU+h89nXl5eV5XPzavK2N88OG0zF8WlIYk7UHNMRnpGAnQ6x+hzB7KKVX/mWC9PHjqPFIICs57GzmHyLCYKrAZ0QvBO20jGgnQ5Bj1nB+Kx/sl4/IayPIpzzOxvHuBQp4t8i4Fj3fECc5/aVkNFgZVck4FbVpfzzWdPogFmA9QX22gezLQClfEkrlfO4QtF+NTWGj6zvZrfNfbiD0X4l7fa+NWhboQAk1FPMKJxsi9+nUFvmBVluRj1OrpdASwGHfdcV8V9m6rYWFOY8lFMrC80FVPPqMWiyiieb7LyAUgpd0spV0opl0spv5XY9o2E8EdKGZRSflpKuUJKuT1d2Esp90spb5xwvTEp5VYp5UYp5Top5Z9JefUWnGlJVGMcGovS7gzMqO3h74720Ng9IWwzBmNaXKDPFG9YTmlbPD88xs/e66J7WEX/zAd2m5HPb69FJyBZQXO65b/LF6Zj0IvRoCMUlZiMBsbCGjo97D09QDASQwJG3cU/8NP9bv7TTfXcvKosNcPcsdzOAx+sw2LSMzwW5EyfB71OIAQ8vr+NPneAz++oo88doDkRFRaKwmgwgiHDJLVnNIgnGGWF3UrfaIAf7m8hEpV87sZaXjs7xC8PdTPkC1FoNRLRYiQXOzkmHRJwjYXZXhc3D316Ww33baqKZxhHYhzrzlxfKFsmrrC21Bbx9Y+smVFC2XwwX45gVQxuHmgf8c/63OAUYZ2zJTCF5hgJagyPRXAp68+8kGM0ENJiPHT7KtZUXFzuI5Pz8/H9rTx/wsHKsnxKcowIBLUluXQ6/YwENCKJRyYSg4lhAu5AvHroxGijuhIb25YWEY5CICrJM+tYWZbHAzvrUsmHQ94Q5TYrRr3AYhBIKZEZHk8pYc8pBy82xUNDXQGN3zb2UJJrpizfxOaaAoqtRs71ewlEJFYDNJTlEotJVpXnUllo5aaVpfzHm5azzJ7Lmb4Lncs6h8fGmYKyIdkB7+kj3akcmiaHe1Ha/5PMlyNYKYA5JhDWEDFVTEdxMTpgSb6Zfc1DmI06zvS5+fbuMxxsG0oJ/YkVNZscbpoHvOxYUUpYi/GRDRWYjTrcgQh5Zv20hgyTbnyJ53gkzQiP7W/hYLuLmJQUWvS4Ahrff62FUw43jV3xY/7vnmaOnh8lpsVt+uX5ZnS6iz8xmOhJEQjHsBp0bKrK5z99eDk/eaeTPU2DvHnOSUNZXqrGVEGOCceIn3BUctNKO7+/rYZVS/JpHvDy3VeaCUUlD92+kod2reK+TdXjspezYV1lAbvWlLH37CAgUyusxZy9m97YaS5RtYDmmMauEfyzsdUormr0wMolNrzBKHetrxhX4bJzOMCB1u54w5RoNFFKIf4Qrass4E9va8ATiPKD186hxWKsqcjjrXODFOWYGQsFpizeF4xdcDIe7oyHUH7uhhpuWVlG25AHi9EIMcmvjvQQiUn+4bUW8ixGvnHPWq5fVkir80Krxu6RAOG0FWqJ1cBwIEq+xUChxcD1dTnYzHpO9Xo40+/llGMUgFA0BgJ0AvQ6wYAnjCRuhT/Y6mJf1MlYKMrLTf34wxpmg2Bnw4UQ8Exd66YimTyW7GeQnO1PLK+yWFYDycJ5yRaeM73fmaAUwBzj9AVn3B1KcfWjAecGfFQWWlm1JD/RnjFuj66329hYXZBSCjaLAbMh/lNNCu/ekQBWo47fHu1FSklYg9HgzLrIJQWgJxDl14fOM+gLs7TYytaaQvQColGNcAwCkRBtg2NsWVrCrw71IoGyPBM2kwGTQU+pzUSu2YAWixEe8KLTCT53Qy0lNjP//ZlTeAJh8qyj3L1xCUfPu7m+tph+dwC9ALvNhMMdD0uWwJl+L3qdoHN4jD+5dWWqqmjS8Xs5hXN6Bc+kMlwMvX1TfQ/WlM15KKhSAHPMiydVMTVFZlaU5fKfb17Bltqi1Ay03m6j3elLCbukUpgoCPac6uPo+RGknH257qQAPNAyhMmgpzjXyLbaYopyjBTlmrAYdFhNes4P+znaPUosFncuWw2Chz++jjyLkVBU8vN3O3kt0de6JNfIfZurWWbP5f3OYUb88VySo90e8sxGPtRgJxKVHO12Yzbq8Yei1BVZcfkCFNqs5Jp09HtCvNvh4mMbq9jZYE8J5wd2LsNs1M1KETR2uXjk1XM8dPvKcauJJIspSzi978Fcr0aUAphj3u+4GspXKOaC1sExhsciWE16DrQ4eWRvMx/dUMn+5niz96/cvIJimynjjLTEZkavEwQjEpsRfFMnlV9E0r+QbBH5l3eu5L2OEY50DOMcC6PFJCuX5PPQ7Sv5xzfa2H92EINeoBfxxK732l1sry/mmaO9HDt/IedzY1U+S/LNPPzCaTZW5aPTQSwWN/ccPj/C4a4RluRbuL6uiJZBXyJiKIDFpKN7JIBOB4UWI+FoLDW2C03htVnN0gNhjVMOD55gJG5+StuebvZZ6Jl/+njmayxKAcwxowHlAFZkpiTHQF1JsvezJBaDiKax3J7L0409VBRY2FhTmBJQ6QKifchHMCqJMXPhn+wj8eTh7oQC0nHvdZXsOztALAaxWAxPMMKb5wYxCNi1dgmhaIzDHcNoEtYtyWPT0iK++VwTTl8YQyKUxKiDQW+Et1uG6HCO8YF6O1UFVjpdAQoterxBjfVVeTSU2XiteYiRhHMsBvjDMSx6yLHosdtM6BC0DHgJRWOcHx6joTyfTTWFs5qlNzncvHyqH50QmNPiVueyP8BsWIjxKAWgUCwQt61dwo7lpYkQT8Fta0p58nAPH91QwX/8cLzl4qP7Wnnwlrg54Okj3ew+2U9VoZkDbcPj+j3MhMf3x4u/ba4t5q2WQWIxyYleNz0jATQJm2vyuWllKXtOD7Dv3BBHu0cpy7fw4dWl9IwE+Nt71jPgCZJj1LHcnos3FMGo1zHsCbKxOp/TDjdhTfJ22yB5ZiMV+RrhaAy9XuAci+A5P8pYKC78dVwIVY1KMOj0tAyMsaW2kKeO9OCPaAz7wtQUW/n2JzemmsLMxCewrrKAP7l1OZ3DAdZU5I/bvljMPjB/kT/pqDBQhWKBSNbiaXK4+eGbbbzZMkS3y8+LJ/rINek52OZi1+qyVMjiS6f6cQfCPN3owOEOYZ3l9O1Er5tV5Xl86UP1PPzxDXxicxWtgz5WlucigPWVeQSjkkBIwyAkOgGbqgs40ePmlMPDU4e7CUU1NCkIaRqBSIxgKEKt3YbNbOBUrxeATqefGBJPUGMsFKW60MKda8opyjGnkr/MBkG+SUehWVCaZ2IwUX+qotCMzWLgpgY7VUUWPrW1OiWoJ4ZvpvfYyJQ8ZTXpybeaxvUVSG5fTMXh2p0+9jUPjRvjXKNWAHOIy6eKqSkyYzOAzWTgibfa+fCqMnLNet5rc2Ix6bljfRljYY0v7FjKppriVBlmm1nPu+0Xqm5OldA3FRurCmgZ8tHnDrCltogXTvTSOxpgW20hJoOf9ztGaRkcIwbENMi1GLh1TTk1Jbn864F2njvZh9Wkx6gTGAw6BtwhPJpkyO8jGNWwmnQEIjF0OoHT50evg7EwtA8HGD3uIN9iwG4z4fSFCUQl+XkGqgtzOJrIeNeJeIe6P7m1gU01RSmnOMSjderttnEz96RCuGVVKfuahy4yoQTCGqGINq7c9WII95zIQoSlKgUwhzxzrGehh6BYpPii8P3XWghGNN5rH+ZErxujXodep+Ngq4vWIR8P3rICx2iQ3Sf7KbQaONY9wuUoz7RjuZ2bV5WmVhZvtzoJRSVObxiTTtDhHCPXpMMbjmE2wIoyG6uW5LOppoh9Zwc41DnC+x0uvMEonmCY1UvyOO8aYyQQo2M4QKFVR77egD+i4QnECKTlClQVmPFFYmgxDYtBEI5K/CGNjuGxlGIryzNzpGuUm1cFuHVNeSoOPlOoZrpwX1ORn4r1T5LsIbz37CAP3rKCdqdvUdn9IemQdwEi1Y97vsJSlQloDukfnVlctuLawSRgeWm8/MLxnlGGfSFCWtz5eqI3HiLZM+rnldMDlNhMvNU2TESTWA2XXrTsJ+920pFI6Kq32zDr49dsHRrDZDJgMerxh2MYdRCOwtstw+xp6mf3SQcguLG+hEA4ikEvsBgNDHrD+MIXhHwgHGMkEKXAauS66sJxnx3WIEevwxvU2LS0kOJcE3q9IBKVqdmoFotx7+YqKgst48w5mWz2yfDOUFQbFz6bJBVTv7qMeruNUCSWKm+xWGhyuPnuK+d4+IWmVMb3fPknlAKYQ356sHOhh6BYpIQlHOp0cbLPQ47ZQG1xLpEYVBdayLPoqcgzsf/sECaDYEtNIR9YXsyHV5ZC5sZ5M6Kh1MaeREP1dqcPT6KsbHm+mUqbiYpCK4h4h7mohByTHrNe8G/vdPLBhhL+8IP1/Jc7V/O5G2spzDWwpMCMKU2SRLR4Vm/Sebu9Np5YBtA86KNpwEc0BkiJxahDhyAYjZFrEeSY9Hx8UyUVBRb+4bXWcWWaM9vs49U8Wwa9qbLO6T6BUCTGlz60jPu31tDu9PHEgQ5ApGoBLQbq7TY2VRdgMepJ1jiaL/+EUgBzyCJ5vhSLFCnj2bY3LCvmgw12DLp4H+ixcIymvniM/JvnnDz+ZhtvtwzTeH6UBnvOJX9uy6CPDy4v4URP3Ldw2+oyjDoYDUQ42e+jecCHlKSijLzBKG+3DTPgDdHk8PKvBzrodI6BjJtvWgbHqCiwogNqiy3cta6cD60ooTTfxNkBH9fVFLKsJIeqfOM4gdM7EsAbjGLUQ65JsL6yCJAEo5LdJ/pxjoUJRWMpgd47EuCHb7TywzdaU/61ZKJcQ1leqqxz0iew+6SDJw50YDbqU60pv3LzckAuqhpA7U4fHcN+7r2uEhApxeTyhfnZwc459SUqH8AcouS/YipigCek8XrzICU5JrREwpSQMmUPv666gLpiK2+1DtPlClxy1fpym4EHb13Bu+3DfG/PWV484eDW1aXUl+ZSnmfm/Q5XvLd0mgJYVWpjWUkuJblGWgfHKMgx8i8HOjDodGgxyXJ7Dl0uPwjoHQniDkSpKrBy++pyghGNnx7sJKiBIdELI0Z85tkzGkISbyAO4PQGsBoNrC7PQ6/TcbTLhdmgT5l51lbk8eSRXnQCck0GPr+jLjVTDoQ18q3GlMkkWexuok8AYE3F4gr/TGb+hiKxhMLSsa2umN0nHXxvbwsAn99RNyefrRSAQrFAmHTxpf5YMMpgoiBac7+XcAxyTTrGwjG6hv0093nQZFwlXKoPWItJttUVc6TLRUSLceT8KE0ON8GopGckSK7FwMhYlBjxbmWahI6RMVoO+viz2xqoL7XxozfbcPnCbKwpZMDjp2XIj6ZpLMkzMeQNMxaKYjLq2Ht2kHAkRjAxE4rKuBJI3oTVwLhCic1DAUw6eLttGE8gwvZlJaypyOdMnwcElBdYqS6ycmN9caqZTZLkDH9iJm16IbXJEq0WuhBcuhJLlroAUvc48V4vJ8oEpFAsEOFYolw4YNHFLRiV+VYsejAgEcCgL0xEk4zNMNt3MgLRGG2DPhrK8thZX4wAcoxxoWcx6ripoYyiHCM2ky61AvAFNfLMBrqH/VQUWLCaDORZjRTnGrHbLEQ1jRyzgSFfmIgEnU7QOuDB6Q3hC4+PVa0qMmPQxVcBEh0WPRh1ApMOTHrB8jIbpx1ubGY9TQ4Pu086qLPn8olNVXxiczX/4+Pr+Mj6yoyCOingG7tcGfMB0stgH2hxcqBlKCX8pzMJzUeDlol2/2Kbic/vqKPYZpqzz1QKQKFYQCIxCGmSquIcdALah32ENPCEZWq2H7iMrTnHwpLvvtLMD15rQa/XkWPS40okFPhDGk5fmFgsRjAaS32+Dlhmt7H7VD9vnhukrthKRJPsPT1Ez2iQmIQRf7wMdE2xGb2I2/GLrHpCaSGgNpMOo95AvT03fl/R+OpgZbmNPIuesCZpG/TR7wnROjTGusp89p4dZM+pPn53tJdOp++iZjbpgvmCjV9kFOhJAdvu9PHdV87y8AunaexyZRVxs5h7B1wKygQ0RyyWCAPF4kYCy8tsjI6F0SQEInJOy4fnWfRsry3i/54b4nREw5g2BYzKGG+3OlPjSlJRaOGv7lzFD99s4+VTA7gD0ZQdv6rQhC8Yxh2EYX+UYX8Ugw7MBj0j/gu/AR3gC8focY1dFBzRNuhNmYl0Orhn4xJWlOXxseuq6HMH8ASiKQfvxGY2Sbt5eqVQ11iYykIz77Y7qbfbxs2gk3kDd62v4OVTfYDIqhDcYisbcblQCmCOSA9fUygmo77EijcQIRCJS8C5Lh24qbqArpEAgYhGQ3kenUPe1L7EEFiSb6Dfc8F0o8kYTx46z3sdLryBaKp5Swzod4fH1SQy6WBdVQGbqwt5+mgPYyENox70Oj2+sAZCsHVpPp2uAMO+MDGg0GpgeCxKRIJegsWk52C7i611xWkOXgP1dluqZPbukw72nhnkCzuWcsuqUkJRjScOdPCVm5dzts/DbxsdABRa42aUZAG8ykILP3+vmwd21rG+sgCQBMLatLb/xVAtdC7IygQkhLhLCNEshGgVQnwtw36zEOLXif3vCSHqEtvrhBABIcSxxN8/pp2zVQhxMnHO3wtxGQKcFxHn+uevnofiyuX8cICe0SCjgSgCqC2yzumsbNATZnNNAbvWlvHvb6ilNM+S2qdPSAO9MIyLNlpRmscrp+Mz/9oSK0U5BorMYNYLNDl+tSBjcNrhYSQQxmKKV/YMRIkLfyAclXQOjzEyFibPEhe6/b4olsRSRGfQc6zbTUOZLVUULSl8z/R5+NaLZ/j718/xyukBdq0pw2wwxFtqGnSpGfrdGyr509tW8Ke3NaQcqMmIms7hAF+5eTlbaosJRaM8/MJpDrYtnpLt89UMPsm0z5oQQg88BtwO9ACHhBDPSSlPpx32h8CIlHKFEOIzwHeA30/sa5NSbspw6ceBLwHvAbuBu4CXZnsjC0FyVnH3hsqLHDU/ertlgUaluJJId5FKoGtkbrPHozLGE2930OcJ0tznZXjsQox5Q2kuVpORFXYLzxzvJ5JYjmgyhsmowxiCygILnmCUkBRE0nwTukTimE4P0ZikscvFkCd80YomBgQjMaIS3EENqx4CGkSlpNiqZ2NNMdtqizjcNUK70zehHaLEH45ytHOEe7dUscyey5qK/JTgT87irSY9//HDK4ALAvWW1eUA436rncMB+t1BOocDCx4JlGS+S0JnswLYDrRKKdullGHgV8C9E465F/i3xOungNummtELISqAfCnlu1JKCfwEuG+mg19okrOKeIr8eAbdlylsQ6G4jGypzqe+zMaoP0IkFiOqxVKz/bMDY5wfHuOV0wMp4Q/Q2u/l5oZSbl1TRofTz4g/AvKCuUrHhZyBkBZ/3T0Sj/HXZ5AC3rQWZoHkRFfCddWFvNcxTO/oGF/60LJUUbTkjHhLbTHf+Ng6/vojq6krsfHI3nOc6fNMmTGbFKh97sBFETX3barir+5cxX2bqhaNk3e+fQ3ZKIAqoDvtfU9iW8ZjpJRR4rkdJYl9y4QQR4UQbwghPpR2fHqltEzXBEAI8WUhxGEhxOGhocWzVAO4ZXU5n9xclZpdpKNcwIpsmc9QvKeP9uMNRgHJgDtEYW5cIC4pMJFv1TPsj+ANjXdDR4H3z49wpGsEhCQmwWDQUWSJjzx9lm/Wx0tHlOYnehjPYDJtNRuIaDHebXdhToSmPn2km0dfb02Vbuhw+lhTUYDZIBK2p6ld5pkEalKpJJvFtzt94yqMzrcZJp35LlE9189eH7BUSrkZeAj4hRAif5pzxiGl/JGUcpuUcltp6cW9PBeSTqePxvMjdE6o360igBQzYT57xtnMgpXlNm5qKMWeZ0KHpMCiRyDIMerGNWhJMjwWxe2P4PRFGAvFKLYaGA1ojAQvHnlIg7GwRigco67ESmhCyeocAxRm0AqhqKTLOYZOCLbVFhGKxNjfPMiP3+mkpsiScvwmV9xbaov5+t1r2FI7tZlkYuP39Lj/g21OHnmlmb9/rYV2py8leBd6NTCfCigbBdAL1KS9r05sy3iMEMIAFADDUsqQlHIYQEp5BGgDViaOr57mmouWQFjjpZMOfvhGGx1DHt7vGP+P9UqTagSvWJyU5Fl56nA373UMs6O+mPKCHO7ZVIndZibXbJxUGYUjGnoBQkiklNQUmqYsS6FPtF+0mMaLmFAUAsGLBVsMODfgI9eko8Pp5//tb+VA6xA9IwEOtDjZfdLBjuX21Io7XbCnJ3VBZgHa2OVKFYtbV1nAAzvreLNliKeOdFNfaiMU0VLHL3TI53wqoGwUwCGgQQixTAhhAj4DPDfhmOeAP0i8/hTwupRSCiFKE05khBD1QAPQLqXsAzxCiBsTvoIvAM9ehvuZF5ocbh5+vok3W4dxBWL84v2ucf9YTx85v4CjUygmZ2QsjHMsgjcUY9gX4a/vXM222hLK88z0j/oznqMnbre3GuIx8yNBjZ7R8JTGlw+vLMVm0hOIjFcpZgOEJjkxImEspHHk/CgmPYSiGgVWAwU5Jn59qJs3zg1xdsDLjw90pAqkNXaN8PALTXz3lXOp32CmjmHNAz5iMh7AajXpMRv1tA+N8amtNdzUUMITBzpTxy90p7D5VEDTRgFJKaNCiAeBPcSfhSeklE1CiIeBw1LK54B/AX4qhGgFXMSVBMBNwMNCiAhxJf+fpJTJAPn/DPwYsBKP/rliIoDq7bZUyBxATMpx/1ie2bZqUijmmLqiHAZ9QRyeMJqUnHK4+fE7HfS5Q5OeYzPriEmBL6zhc8cF73TJaq+eGWQsHBu3otCL8bV/MmHQxbOjw1HJy6cGQUrODngpyTXhD0ewmfQ83dhDbUkO92+toXnAg8Wg556NFanf4EQB2tg1wvPHe7lr/RJCUY0DLU7WVOTzp7c1pI7Jt5oWTZLXfOYcZBVyLKXcTTxUM33bN9JeB4FPZzjvaeDpSa55GFg/k8EuFs70uel3Xwif06JaqvPQ3Rsq6RlROQCKxYWeeGDC4FgQd1DDohec6fPQPuQjpE3thQhrMWqKc8kzGzjR7Sab+LZgJHaROUmTsCTPSL7VzLnB+G8krTYccEFBFNtMfKC+mNP9Hm5cVozFZOCFE/04vUHuXLeEuzdU0uRw82aLk09vq+b+rTXjwkCTheEqCqy8erqfWLz4EN95+Sw5Jj1/89G1VBRYeeSVZr64c9msms1fCosl7FTVApoFnmB0XGs+q8nIQ78+yndePsszR3tw+ufTradQTI8pMdXz+MOsLM/FYop35YpJ2FJTyLJi66TnRmLxstT5VhOGLCRGsVWP1ZT5QKc3QuvghQlSppVEPHRUcm7QRywmsZgMdA37yTPrGR6LoBdxBbGusoAHb1kxTvgnSdr8//nNNp474WD7smKWluQSk5K71sdXCz8+0MGP3+nkxwc6gPm1vS+0ozmJKgUxC9oGveMe3D5vmAGvE7vNRF0ie1GhyIaJM+C5ImmV9EWg3xMix2TAaoT8HCMH21ysrcyjw5U5Cc2sF3QO+znUNZrVZ7kCGpPIf7IxjmoSDrYNo9cJ3AGN1872c8faShpKbXS5/Fy/rGTSc5Mz61BUgoBNS4sw6HV8cecyOp0+8i1G1lfmYzXp+cz2pfR7gnxm+1Jgfm3vC+1oTqJWALMgkqE6ow5YWpJDeb55/gekuGKZD+GfjkkPm2oKWFVuo6HMRpdzDCEkx3ounokmZ/tjEZm18E+i18V/ExMjPgssF0ROnkmQ3H1xRJFgSb6VYDTGiC/Cq6f7OeVwU1VoxW4zp0yu6TkChztdNHaN8Oi+VjqdPh7atYo8i4HDXS72NPWzpqIgFToaCGscbHPiDUYZ8cfNuZfi/J1p6OZCO5qTKAUwCzL9I2vA4a5RvvPSmfkfkEKRJVajnvfbR9h3bphDXSPkmgyU51tSheCSFFr1GC+hOpfVoIs3ldGNFzHuRO6AAEpsZkpsRuBiRegLRVlakoPNpOe66kJcY2HeaBnino2VbKktorFrhCcP91BTbGHIG+KX73Xx6OutgOSmBjvPH3cAErNBhz+s8fzx3oti/V85PUBDmY2KAmsqlNTlC88qBv9ymnQWWx6AYgInekcv2pZ8gLucygGsWFyk/8j9kRiVBRbyTHqiWgx/RGNtRf5FgsAd0C6UaZgFrmAsXgI6ktkfJoFOV5BB38UuZR1QU2Sl0GIgGpNsqSvkvk1VCAkVhVYCYY2Xm/oYGgtysM3F9/ee45njDhrKbAkTrEiVj95SW8yf74pH+wx5Q+Ni/e9YW07LoI99Zwd4ZG8zj7x6jt0nHbMS5JfTpJMpjHWuFILyAcyQQFijZ2Rs0v1ThdMpFAuB1QhLCnNwDPsJS0nLkI9PbKrkzdZhekeDvNQ0cFHEzuUwTc02FCIGtA75CUSiBKMx+kYCjEUkw74QT7zdQduQj98d7UUv4Ia6IpaX5oEQvHZmAHueiZdP9XPX+orUXQz7wjx33MGhTheleebUKuD+rTWsrsin3m5jmd0GSNZUFGTsIzwd6W0dLzWS6OIw1nhP5IduX8nOhstbDUEpgBnS5HAz5AlPuj8830ZdhWIaxiIQjUoCMTAbBA1lNoa8YQZGgwBE5yBozSjiiV2zwQAU55noS4Rav3JmEC0GkZjkUNcIo2MhttQU8G7HCK+fddLuDHDvpkp0ifqTOiEw6gRPHOgkFI3hGA2wotTGZ7cvjTeMSaviO7EXr9Wkn1CBdGZcDmF9cR7AhRXN5UaZgGZIvd1GTdHkIXMKxWJDL6DfEyTPoiMSlZx2eNnbPJRVRM6sP1MXLxE9G4wGwZA3nAq1FrEYOiGRifctTj+alBj0Am8wytqKPCoKrNyxrpwl+RbuWFdORaGVB3bW0TLo5Z/e7uBo9yinEyaVx/e3psw9cLlDMi+fsE6uJtZU5PP1j6xhS23RJV9zImoFMEOOdY/SOjy3NdsVisuFSQ9SxmfF3oQDNtfIuCbzmQrAXSpRLd7eMTaLVUC+VU/YF8+10QuwWIwY9ToabGZO9nqQQN9oECmhJBF6/c9vt+MLXlBpNrOBr9+9hqXFOeSb9SAE2+tLaOwa4Z22YVaW29ix3J5qEfnAzmWXxX6/pbaIr39kzbTXyiYRLKmYHthZl6qOerlRK4AZcrpvdKGHoFBkTViLJ3IForHUj12n04+r059rvvxiIMrsTUsD3ihlibr9moSiHBNL8izcvLos1UWsdyTAzStL+c79G/j01hoe2rWKP9+1kk9uruKPPljP2iV5DHlDmA16qotyuOe6Km5YVgJIRgNhTvS4eaN5kCaHmycOdAIyFU56KaRnIU91relWHemKabIm95cDpQBmzFXVuVJxDREDlpVY2VhVkDKnAHhD00vqmTz1eqDMdmm/E0/owhJlRVkuf3JbA0adjsIcExur8vn4pir+28fWceuaJRTbTOxssFOaZ+Zgu4vXmwd56mgv//eVZjqHx7h3czWdzjHanT621Bbz2e01FOWa8EeiqT4Al1PIZmNSmi5qKKmYzEYdW2qL5ixpTJmAZkh5Wg9VheJKY9gbZjQQmrHJZyaWnBgw6o+fMZtMZwFsqCzg3c5RBHC4cxTQ0TPixxuI0O0KMBqIsKepj1Xl+WypLUrNvHetKePXh7opyTWxuaaQX75/noduX8mDt65ImVs21xTz7DEHzx7rY3NNETsbSsc5gi+VbEJCpyv4ln6NuSwOpxTADHmzZXChh6BQzBpPWENMHsR2WZBAsuvjbAKBJJBrMWJM+BBCkShvtwxx/9YadDFwBTxEtRjPNvZiMPbx9Y+sGRfaucyeCwhOJZK9+twhPpJoDg9xO/1nt9fy/PFeQCyKwmwTxzBfFUGVCWiGtA+qRC/Flc2VEKn8+tkhIjEwGQRrqwrYVltIocVAj9efaEwjqCqy8kcfrB8307aa9OxsKGVng51Pb61J9fxNT6aymvR8dvtS/uaja1NZxd9+6QyNXSOXZezTXS9TYlcms9F8ZAQrBTBDdMoHoFDMOZK4cNJiknfbR9h/bphH97Uy7IuSY9SzrbaQt1qdnHf5aexycaDFeZGgLLaZUo3gJwrl8bV4ZFb9hWcy+lgMmvs9GYV3JmGfyWw0HxVDlQlohvT5VAio4spmviqQXiox4lFMEB+vxajDF44RiGp0Dfspz7cAkkdePQcCvv6RNanewXdvqEzV/FlXWUAoGsUTiHK0O96Pak1FPse6XbQMjrG02MpDt6+ctr9wtmypLeb+LVXsPTvIxprCi0w5mYR9JpPPfFQMVQpghhiEWgEormyuBOGfxKKHtZX5tDvHGA1oGHSgxaCp18OdGyr42MYq1lcWAIJ1lQU8faSb7+1tAWB1RT6Pvt7KrjVllNjMhCIazxzp5bUzg2yrLWbvmX66RwKU2cz83aevu2z2//QyE5mEd7b2/ctZXmIylAloBgTCGsFwNv2QFArF5UAIaB0awxvQMOrBpBcIYGV5Hmf6PDT3u9lSW4zZGBdld2+o5C92NXDL6nJCkRg7lpfw5JEe3jw3gMmg476tVdy1voL3O1wU5hiJxsAXiuD0hXnirXZePzM4pc09ENYuakKfictZ7nkuTUFKAWRgMudLk8PNqKr1plDMG6Eo+IIaeh3YTHr0Oh1GPSwptDLkDdEyOMbTR7r5wd4Wnj7SjdWk5/M76uhzB3jiQAdGvSAY0Xi3w0UoGmN5qQ2jXodOB9tqi8g16ijKNXPs/Ajff62F77x8ZkpB2+RwpyqHzlc3r7k0BSkTUAYmpmBXFFjZd3aAzUsvfy0OhUIxOcl8hXAMtKBGvkXHiooCaoosVBRWsLTYyq8P97C8LJenG3tZZrexs8GeEpr1dht1JTZePzvAyZ5R+twh3mwZ4v4t1Qz5QkRiMZaV5LJpaRFLCi00lE1dCXRdZQEP7VoFyMsikLMJQVV5APNIPAU7xgM76wDBo/tayTMbONDm5HqlABSKeUMnwGLQEYjEkMTLQoyFY5x0uDnlcFNbkssyu40v3FhLy6CX4+dHyeTh6HT6ONY9yvW1RXx4ZSkmveCW1eX84xutCCFoGfLxk4OdfP3uNdMK2niYqf2S7y3VujKi8cSBTr5y8/KLPns+8hOyMgEJIe4SQjQLIVqFEF/LsN8shPh1Yv97Qoi6xPbbhRBHhBAnE/+/Ne2c/YlrHkv8lV22u7oE4inYHZiNerbUFrFrdRnDY2Hu3VhJjlnpS4VivojJeAObZNiFPddAmc2CTkKOSU+eycBTh3t4s2WIA63D3L+1OhXJ09g1wrdePM3fv97CM8cdeAIRjvW4eelkH//4ZhuPvn6O9qExfm9bNZ/aWs0d68qpv8R+3knTcaauYhPNykkrA4hJzTvzEQY6rQIQQuiBx4CPAGuBzwoh1k447A+BESnlCuB7wHcS253APVLKDcAfAD+dcN7npJSbEn+LIsV2or2tsjCHXavL+NDKMlqHVBKY4urhSolnS5qBhBD0uIOEZbzQXFO/G6cvxNGuUe5YW879W2vSZsoSX0hjf/MApTYzH9tYyV3ryul1jeEYCfJKUx/1pTlUF+awv3mI3zb2cqbPc0njTArsTF3FJgrzpJzZUls0qbN4sYSBbgdapZTtAEKIXwH3AqfTjrkX+NvE66eAR4UQQkp5NO2YJsAqhDBLKRetKzU99OqX75/n5+91ImPwgeUlfGhFKSd7L+0hUSgWAwbAYgLfHJeFuFyYdFBoNeAJRNAJKLOZ8IZCrKss4L7NVdht5nHHb6kt5rPX1/DEgU4Otjs52+9BJwR/emsDo8Eozf1e3m13YTN72VCVz3vtLkLR6TskTGWWSfc7JDuNJcM3Jwrz6ez681WeIhsTUBXQnfa+J7Et4zFSyijgBkomHHM/0DhB+P9rwvzz34XIHGAvhPiyEOKwEOLw0NBQFsO9NJJhXr98r4tnj/USikg8oSjPnejjeM/lSRVXKBaaKFeO8Ie4Zb9lKIAEYlLQ6QoQigiOnB/haPcoj+xtvqj0Qp09l/u3VnLvdVV8oL6YUX8Ylz/Mp7dV88mt1dy9voKbVtoJRWLEvQxxETRVCYapyjwkhXqxzcS2umLanb7UrH+mYaHzYf6BeXICCyHWETcL3ZG2+XNSyl4hRB7wNPDvgZ9MPFdK+SPgRwDbtm2b8xyWJoeb777SjC8U4VNba2gos+ENRni7ZYhXTg/M9ccrFIoJJGepJn08CUxKSX2plc9vX0bXiJ91FXkcbHPiCYY50DJEKCrpdPp4/kQf3mCEcDTGf7ljFdfVFFFZaOEfXmtlwBskGIkRlTEikRglNjOdw/7UzPvx/W2pGfv4mfjUZSPSZ+7ps/7JZvTJ7fV2G+1OX2p/vd3GLatKL9kvke13OxW9QE3a++rEtozHCCEMQAEwnHhfDfwO+IKUsi15gpSyN/F/L/AL4qamBWddZQH3XFdBjtEAUhKKxjh2fpR32l24g3NXlEmhUGQmBlgNOix6sOeZiAH+cIxwLMaBVifvtDnxhaK8fmaAbz7bxDefPcmzxx3cs7GC62sLGfKF6HcH+PyOOnYsL+WeTZUIIfAEIwRCGjfWl/DBFXaePdpDY9fIOME9cSa+pbaYr9+9ZtKyEenHp8/6J5vRT+Y3aHf62Nc8RLtzbv2O2awADgENQohlxAX9Z4B/N+GY54g7eQ8CnwJel1JKIUQh8CLwNSnlgeTBCSVRKKV0CiGMwMeAvZd6M5eDeKXAWkx6HU839jLsCzHiD3PXunKebHQs9PAUimuKIosOXyiGJ1FfWpMRrEbBjcuK8YcjdLv8OL0BvCGNvWcG0esEeRYj915XyWdvqOVg2zCN3W6WluSm7PH3barGH45ypHOELpefG5fb2X3CgcsfIRSNjrPPz9R2P5njdrrtSb9Bcv98OIAhCwUgpYwKIR4E9hBv9vOElLJJCPEwcFhK+RzwL8BPhRCtgIu4kgB4EFgBfEMI8Y3EtjuAMWBPQvjriQv/f7qM9zVrXL4wzxzroSTXzF3ryynOMXO8e4SSCU4mhUIx96R3KzPp46GhkZjkzRYnN1FCic1EfUkO73eOIES8F3C9PRenL8y/vN3G8lIb3/jYOkCm6gIts9vY1zxETJM8sHMZdyd6BTx1pJuJRpGZJmFNdvx015m4f776AQgp59ysftnYtm2bPHz48Jx+xs8OdvJ3e5rJtxqx28xsqy3inVYn3aN+PMoEpFDMK2a9ICYlBgFVJTmcd/rR6QRaTFJgMVJkM/HQ7avocwcoyTXzclMf+885iWoxQFJqM/MHO+poKM+nZcDD7lP93L62LG7CF4LlpTZAcqDVycleD3995yp2NpRe8rizjeI53OlK+RsyCfzLFQ0khDgipdw2cbvKbEojENaoLLTwn29ZztLiHM67/Dx7zMFYKIJXCX+FYt7QATkmwSc3VfO7Y714wzFaB/0sLbZi1IFBr8cXDnNjfQk3LCuh2GbicKcLtz/KB+qLiURjnBvwEtIkvznSQ2GOkY9uqCSYqAga1OIrC7NRhzcYJRCK8pnttdTZbfzsYCd3b6ikONGYfjakO5KnEuzJnsTppp50oT/ddS4VVQwujf3NA3zzudMsLc7hhmV2+txBRv0R+j1BKvNn/zAoFIqZkWcBi9FIrydINBYXVHoBPa4AAsF5l59QBA62DfPMsV4CYY11lQXcvWEJTl+YPk+Q29eVU11g5fe2VfPQrlXcuW4Jn9xSzealBcQkmPQ6PrDcTnm+hfu31fCVm1ew7+wA39vbwu6TF/x9U2X4TkY2Td8f399Gu9N3UXhok8PNo6+38vSR7owK4nKiVgBpHGgdpmckXkUwPvvvZUWZjVBUY9CnykArFPNFKCrwh8McPT/CjfVFHDk/gicYn7V3DfuxmnTcuKyIY90eftfYQyRh8inOMVNqM9I6FKY838Jf31XJltp4Da+nj3TzwrFe+r0hcs16blpZypduWk6fO5AysST9AXen9RBOCutbVpWyr3loyro96eGcs3EWJ/ftWlOWCjsfn+F8eVEKII2dK+y8cMLBqV4P19eWsKI0j353kHp7Lo1dows9PIXimiEYjfsmXf4o7UNjjAUvOIMjEmQkxog/QiASodaaww/faMUX0sgx6ghpknyLgdfODHFjfbxw29NHunnl9ACVxTm0OcfYUltEp8ufEv5Jk0uyjWQ6U2X4JgVzNkoinamcvMmGMgB7zwyyuiJ/zhzCygmcRiCs8Y9vtPCzd8+ztiKf9zpcSMCoF/jDsWnPVyiuVhZDG0mDLl4ddH1VIW2DPqIxjZGARr7FgD8cJddkoNBqpKHcxrqqfK6rLmbH8hIOtg3xrd1n+UB9MauX5POD11v4qztWUVmYQygapXM4wJstQzx4y4qsyjN4AhEe29fKQ7evZGdDKYGwlsgMlqypKBiX0JXp/GwdupezHMRkTmDlA0jDatJj1OnwBKOc7HZh0kF1gYVtNYULPTSFYkGZb+FflqunPDduoNALyDPpuWFZEaV5FmoKzdgsBjZUFbKkwMRtq0upLcnhix+o49bVZQx4Qrx8agCzQWA16ekcDtDvDvJuxwgtQ2OEIjFc/ggdTh//8Forzx7toaHUNm3WbXKW3zk8lqgaIVLbnzjQAYhJhX/6+dmWd7icXcUmQ5mAJmA06NEJ0KTAH43R4QrQ6VKN4BWK+WRwLO5ozTXpyLMYMOh09I4EGPFHePXMIFpMEoxo/Ied9XxsYyUvHO+lZ9TPuf4xqotzcIwGSAro+zbFS5fVlVhZtaSAFaW5VBbm8JN3O7lrwxL63CGaB7y0O30U26a329fbbWysLkzZ7ysKrNTbc/EEwzx1oHdSH0Gyz8hUpSHmG6UA0giENZbkm7l5VRlLiyw8e8zBiD/u/I0oC5BCMe/8u+1L2VJbTNuQl8MdI/SMOtFicdEe1MKApM8d4LdHe+lzB/nwylJGAxHu3VQFSAJhjWKbiQc+uAyI/8ZXV+RTlGOiriQXf1jjULuTpXYbFQXWKceSXik4nX1nB/jt0V4qCq2ThnSmN36xmvTTxv/PF0oBpNHkcPMPe1vocQf5jzfVs7G6kH3NQ2gLbfxUKBYYgy5eg3++Odw5wqg/woF2JyadjnyLkUA4Sr09LrwjWowhb4h7rqvAqNeza005B9ucVBZaeOJAJ2ajPuXkrbfb2H3Swd6zgywrzuHJIz3kWfUEwzHOjwbZd3bgIgdwJibG5qdHDk3MHWhyuPn711qoL7XxuRtqLir1kMmhPJ9c8z6A9BjfUCTGhpoCNBnjRI+btqExriAfuUIxa6YTBNHY/DaQMQgw6+BEr5tnjvUy7AuzqjyP9ZV5BKISlz/Ch1eV8eyxPh55tZl9Z4dYX5nPG+cGefJIDy2DPh7YuSwl/B/d18rj+1vZc3qAXavL2LS0iKJcIx9aXkpVgYUPrrBzy+ryrMY2MYQzGTlUbDNdVEp6XWUBq8rzePZYDwfbXKlrJFcT6SWjF4JrXgGkV+P7p7fa8QaigKB92Mu22iJyzQtnn1Mo5hoTccGezeR+JnMh6yVIFpsBohKqiqxUF1nZUV/CDXVFnBvyUpBjwiBgSYGFzUuLsJn13L2+go1VBXiCUZ5t7KXPHWD38T7MRh1WU3wFsGt1GWf7vdy5tpy7N1Qy7AtRarPwkQ1L+OwNtQQiGn3u7Hx9kzlnA2GNp4908/evtfD0ke6UEtixvJh7N1Vztt87TtAHwhqeQJidK+xzXvZ5Mq55E9C6ygIe2FlHKBrjppWl/K6xh7I8Mx9aUYaQ8bZyCsXVSrY9YWYSBiqA0CWYi4ROB8To94RYW5FP96gffyjK8FiE9ZX5FNtMtA6O0e8JcM91lbQPjbH7VB81xTncu6WKZ485uHdT5bgKnvdvrUnF8O8+6WD/uSHu31LFjuXxuj+TxffPhCaHm71nBllVnpeK3wf4ybvn2VFfwtJi6zhB3+Rw89i+NhCwsbpgnAN6vpzE17wCsJr0mI16njjQyae2VrGppoBgJMaRzmEc7uCCxz4rFHNNrgHGpumGOJPfgcUgUolcs8GbyLnxR2Ic7xmlMMdEdVEOQz43NrOR//3JjXQOj1FXksNj+9qIajHu31zNfZursZr0qQiddMGZnLUf7nSx98wgdyRWAknfAMCZPg9PHOiYlWM2HuWj8aUP1bOmIj8VDhoIazSU2th9woFBr2NjdWFK0K+rLOCh21cC4qKM4LmuAZREJYJBIpHDxQsn+vjt0R5yTXoCYY3AJTzECsWVglEXNwWNJWbtBhE3wcwWvYBcsx5fUMtoWtID2a6rl5VYCUVjfOmD9bzf5eKudUv4UEMZZ/o8iR6+cc+E2aBnS23RtLPl5G89FI3ROTzGm+ec7FpTxr7mIR7YWZc4SmR9remKth3udPHovlZ21Jdg1Avu21SdVZG5y70CUIlgUxAIa+xvHuJg2xDRqGQspKmZv+KaIRK7IPyBi6LeSm1GjCJ7J7AmwRPUKMjRZxQwkwn/9OsbdXFFInTgDUR5t32YlgEfj77eypOHz/PwC038w2ut5FuNmA06vrvnLL98vytjobaJjtkO5xj/8FoLvzvSS57FwI7ldr5y83K21BYnrAEdWTll0xO7kqZkTyDKgZahcU7gB29ZwfrKAg60Dmfd4Ws+ksBAKQBcvjDffO4UTx7pZuWSAtZX5xHSJDbzNW8dUyxSLjUap8ByQajoiYd4pjNx8uMPx7hvSyVrl0zvqDSkDc4TyLwCmIyJwshi1NE/GqTOnst5lx93MBJv3i4EOSYD96Ts/AJ/ROP5430ZBXe6oG5yuNl7dpC7NiyhusTKgVYnB9ucKWE7k05c6ccmTcmP7W/hkVfP0djl4nBnPOpnW10xW2qL5qXD10y55k1AT7zdziOvnKMwx0AkqlFVlENjt+eyfoZCsVioLrSg00lGfBGkjBGMxiOABPGZe6FFz+iE3hdGASU2IyP+COkxEVYjBOagSK5egEkXnwVvWVqEJ6ShEwJvMML9W+K2/jN9HkCmevM2drmYzHSTbk4BUolZ//hGO6uX5PGVm1fMqPZ/puvV220c6x6hZdBLQ5kNs8Ewa3/CZJ+lGsLMASW58X/43tEQEhj2x4W/AZjGL6ZQzDlmXTyiRhCfXUcueb4mGfSGCWa4kNWgw52h8ZEQ8VVAUvjrRVxZzIXwB9AJQUiTRMIagUiM0UCYXpefT29byn2bq2l3+ghFtXEF2abq4pWewZsuuP9sl35WgjXd3g+kqoD+7mgvCNhaW8y6ygLMRt0lz/jn2hl8TSuA3pEAT7zVji+soSO+9E0mvCjhr1gMJMMpJZdD+EPPaIjiHD1WA4wENAwi3mc3BuQYdWgyxkQzutEg0OsgRw9hCUuLrLQPj4+Zt+jhcjXNi8QkZh0Y9Tre7XDFVyZWA3abmSePnOflk/3ctWHJuIJsMPlsOVM5hm11xZMK1Klm3clon2SSGcBXbl5ORYGVsBajriT3spl5JtYPmguuSR9A0in08HNNHHd4gfGJMFeOUUxxNWHWz322rVUvGPFrxGISkx4qC63YEj4BbzBKWIvPCm3G+EisBkE4KhFCoAkd0Rj0jgbGCQ49EJtg7J/NfRSY49eCuOKLSbh7Qznbagv59NYanj3ey28O9+CPaDSU5fHHNzdwyjHKSycdHGhx0tjluiirNpmc9ejrrYBIlV94/cwgT7zVjst3cSZEY9cI337pDI1dIxc5kBu7XHz3lXOc6h1NlICO2/j73AEOtA6TbzViNelnXPkzExeqjMZfZ9OJbKZkpQCEEHcJIZqFEK1CiK9l2G8WQvw6sf89IURd2r6vJ7Y3CyHuzPaal5tAWOP1MwN875Vm/ujf3uMP/vkge84MzPXHKhRZoQfCGlQUGLHqLzhTk4JUn3iRb579nM2gg13ryqkosGA0GghrcH4kgBaDPLMutcLQAF9E0lCWy6615ZTmmfnwylIKzTpMeghp4ydMGvGVAcQFyuaaAjbXzHzGGojEzU1JNlYXYDMaONvnobFrhIgWY0tNPp/cUsWmmiJaBj1879UW/ueLp3lkbzNJAZ8+W25yuHnl9AANZTbWJBqrtDt9fOflM3z/tfGtHy8gE7NAeZEgD0VjOH0hnjzcwyN7m1PbJzqP6+02bllVekkZvslrgpizchHTmoCEEHrgMeB2oAc4JIR4Tkp5Ou2wPwRGpJQrhBCfAb4D/L4QYi3wGWAdUAnsFUKsTJwz3TUvK41dI/z3Z0/hGFXJXYrFR45Jhzccw+GOUJJrJBaTeIPRVDx+MjTTM4sU23hebXy2/58+vILHXm9h39kB6kusCB0IqaPWnsuJbhe+kEZ1kYXWoQBSi9E66AMEg94Q/qiccnlsAJbarWyqLuTN1sEZjzPZc8lqENgsRpy+MAOeEEIITvd5yLcaONbjpXUoQK7JQCQaQ6eDDzeU8rGNVRkdwPV2G6vK8zjbf6Hc87rKAr561xpOO0YZC2u4fGGKbSZ6RwL8+EAH922u4hObq1hTETcBpQt2s0GP3WbmrvXlrK+8UBJ6YoevdqePfc1DrK7In7LE9FSk+y4uhz8hE9n4ALYDrVLKdgAhxK+Ae4F0YX0v8LeJ108BjwohRGL7r6SUIaBDCNGauB5ZXPOy4fKF+dXhLnpHg3NxeYVixqSXVtADenFBshZYDJx3BVLC35www0w1cTFmcBAnP0OI+GfIWIz32ofpGvFjMRu5rqaQvWcH+eCKUnzBKBuri3i71UkwEv+stuEAayvykEjebXellJBZH7fPaxJCiTrpJn3889udAbqGu4hJsBp15Jn14/pppyeBZcpAthoEu9aWs2pJHstL8zg/PMYLJ/rYVFPAB1aUYjYIOof97D0zyBd21PJfP7ImYxVOiK/6d590cHbAy53ryscJ61vXlOEY9fO9vS3kmvR8fkcdPz7QwY/f6aRnxM+gN8Qyey47G0rHCfYttUX8zUfXTOs8nkk4aSYm+iHmKhs4GwVQBXSnve8BbpjsGCllVAjhBkoS29+dcG5V4vV01wRACPFl4MsAS5cuzWK4F7P7pIPdJ/tnda5CcbnQERfK9WW5tA6OAXGzzIO3LOef3uogqRL8kRiFOSa0mIZO6Pj965dy2uHm7bZ4LfyJikBwQfjr0rYl7Ucrymx0ufxIIaiz29heW8xTwz1YjHoqC6x8/LpKSvPMvNs+zFutTroTE6WyPBN5FiOblxay93Q/7mAUTZN8cEUpTl+YtZV5vNc+zPZlJbQO+jh8fpSaQivFuUaO93qoLcmle3gsrhw0MOkF9jwzg54gkRhUFubS5hxLmZM2V+Vz18YKXj7Vj2M0yI31dm5eVcbWuuJxAjcQ1jKWe5hIMub/zrXlGRurT2wA/8Wd8Z4B66sK+Mm7nWTyZGQrjC9VaM9XKYhFHwUkpfwR8COI5wHM5hp3b6jEMRpg75kBeobH8Kv6blccScEmgXyLDmQMdyj7842A0HNRhEtpjg5vMBafKRM3Q9QUmukajV+8wARCr6cyz4I3rKE3xAWQzx9Bk5Bn1SEQBMNgNsJd6yt59XQ/eVYjNQVWBrwhbGYD19UUYTXryDEZ2LVmCc8d6+VMn4evfmQNlYVWqgpz+fWh82yuLWTL0mLMBj0Qw2wwsKW2iEBY45ljPQgJu0/1c/eGJYyFo0S0GJEoDHkDlOZZMOqT8Wzx1N0ck55da5bwRvMgdXYbO5aXsKmmkOVlNm5ZXZ5qim416am32zDqBU5viC6Xnz+7bSWBiEa93cZH1lfgCUbpcwe4c10Ffe4A9XZbquaNYzTA9/ee4893rcRi1KdMKe91DFOSa6LfE2RpcS5mg45hX5jXzw7wRx+q581zgxw9P0p5nok/v2MNxbkm1lcWAnLS2W+2wjWZhTuZopjYAL6qyMrffGwtgbBGVZF1QZO2LnUFkS3TJoIJIXYAfyulvDPx/usAUspvpx2zJ3HMQSGEAegHSoGvpR+bPC5x2pTXzMRcN4VXKBSKq5FLqQV0CGgQQiwTQpiIO3Wfm3DMc8AfJF5/CnhdxjXLc8BnElFCy4AG4P0sr6lQKBSKOWRaE1DCpv8gsIe4D+cJKWWTEOJh4LCU8jngX4CfJpy8LuICncRxvyHu3I0Cfyyl1AAyXfPy355CoVAoJuOarwWkUCgUVzuqHLRCoVAoxqEUgEKhUFyjKAWgUCgU1yhKASgUCsU1yhXlBBZCDAFdCz2ORYIdcC70IBYp6ruZGvX9TM7V+t3USikvappwRSkAxQWEEIczefUV6ruZDvX9TM619t0oE5BCoVBcoygFoFAoFNcoSgFcufxooQewiFHfzdSo72dyrqnvRvkAFAqF4hpFrQAUCoXiGkUpAIVCobhGUQpgkSOEuEsI0SyEaBVCfC3D/oeEEKeFECeEEK8JIWoXYpwLwXTfTdpx9wshpBDimgnvy+a7EUL8XuLZaRJC/GK+x7iQZPG7WiqE2CeEOJr4bd29EOOcc6SU6m+R/hEvld0G1AMm4DiwdsIxtwA5iddfAX690ONeLN9N4rg84E3irUm3LfS4F8t3Q7w3x1GgKPG+bKHHvci+nx8BX0m8Xgt0LvS45+JPrQAWN9uBVillu5QyDPwKuDf9ACnlPimlP/H2XaB6nse4UEz73ST4n8B3gOB8Dm6Byea7+RLwmJRyBEBKOTjPY1xIsvl+JJCfeF0AOOZxfPOGUgCLmyqgO+19T2LbZPwh8NKcjmjxMO13I4TYAtRIKV+cz4EtArJ5blYCK4UQB4QQ7woh7pq30S082Xw/fwt8XgjRA+wG/mR+hja/LPqm8IrsEEJ8HtgGfHihx7IYEELogEeALy7wUBYrBuJmoJuJrxrfFEJskFKOLuSgFhGfBX4spfy/ib7oPxVCrJdSxhZ6YJcTtQJY3PQCNWnvqxPbxiGE2AX8DfBxKWVonsa20Ez33eQB64H9QohO4EbguWvEEZzNc9MDPCeljEgpO4BzxBXCtUA2388fAr8BkFIeBCzEC8VdVSgFsLg5BDQIIZYJIUzEey0/l36AEGIz8EPiwv9asuNO+d1IKd1SSruUsk5KWUfcP/JxKeW10FN02ucGeIb47B8hhJ24Sah9Hse4kGTz/ZwHbgMQQqwhrgCG5nWU84BSAIsYKWUUeBDYA5wBfiOlbBJCPCyE+HjisL8DbMCTQohjQoiJD/JVSZbfzTVJlt/NHmBYCHEa2Af8lZRyeGFGPL9k+f38F+BLQojjwC+BL8pESNDVhCoFoVAoFNcoagWgUCgU1yhKASgUCsU1ilIACoVCcY2iFIBCoVBcoygFoFAoFNcoSgEoFArFNYpSAAqFQnGN8v8DMGTw37YMVc0AAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "plt.scatter(stars.StellarFormationTime, stars.SmoothedMetallicity, s=1, alpha=0.5)" ] }, { "cell_type": "code", "execution_count": null, "id": "empty-denmark", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.6.8" } }, "nbformat": 4, "nbformat_minor": 5 }