{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Atmospheric models"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The DRIFT-PHOENIX model spectra are added to the database. The grid is then interpolated and a spectrum for a given set of parameter values and spectral resolution is computed. The spectrum is then plotted together with several filter profiles."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Getting started"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import species\n",
    "from IPython.display import Image"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Initiating species v0.1.1... [DONE]\n",
      "Creating species_config.ini... [DONE]\n",
      "Database: /Users/tomasstolker/applications/species/docs/tutorials/species_database.hdf5\n",
      "Data folder: /Users/tomasstolker/applications/species/docs/tutorials/data\n",
      "Working folder: /Users/tomasstolker/applications/species/docs/tutorials\n",
      "Creating species_database.hdf5... [DONE]\n",
      "Creating data folder... [DONE]\n"
     ]
    },
    {
     "data": {
      "text/plain": [
       "<species.core.setup.SpeciesInit at 0x13095f710>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "species.SpeciesInit()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Adding model spectra to the database"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "database = species.Database()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Downloading DRIFT-PHOENIX model spectra (151 MB)... [DONE]\n",
      "Unpacking DRIFT-PHOENIX model spectra... [DONE]\n",
      "Adding DRIFT-PHOENIX model spectra... [DONE]                                    \n"
     ]
    }
   ],
   "source": [
    "database.add_model(model='drift-phoenix',\n",
    "                   wavel_range=(0.5, 10.),\n",
    "                   teff_range=None,\n",
    "                   spec_res=None)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Interpolating the model grid"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "read_model = species.ReadModel(model='drift-phoenix',\n",
    "                               wavel_range=(1., 5.))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "model_param = {'teff':1510., 'logg':4.1, 'feh':0.1}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "modelbox = read_model.get_model(model_param=model_param,\n",
    "                                spec_res=200.)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "filters = ['MKO/NSFCam.J', 'MKO/NSFCam.H', 'MKO/NSFCam.K', 'MKO/NSFCam.Lp', 'MKO/NSFCam.Mp']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Adding filter: MKO/NSFCam.J... [DONE]\n",
      "Adding filter: MKO/NSFCam.H... [DONE]\n",
      "Adding filter: MKO/NSFCam.K... [DONE]\n",
      "Adding filter: MKO/NSFCam.Lp... [DONE]\n",
      "Adding filter: MKO/NSFCam.Mp... [DONE]\n",
      "Plotting spectrum: model_spectrum.png... [DONE]\n"
     ]
    }
   ],
   "source": [
    "species.plot_spectrum(boxes=[modelbox, ],\n",
    "                      filters=filters,\n",
    "                      offset=(-0.08, -0.06),\n",
    "                      xlim=(1., 5.),\n",
    "                      ylim=(0., 1.1e5),\n",
    "                      legend='upper right',\n",
    "                      output='model_spectrum.png')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkEAAAGiCAYAAADgPBIcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdeVhUZfsH8O8M27ANO6IgorjvG4iaimYuaamllpnlkq+0WBnaZv2ksjLN1Hbf1wxfTS2z0txzwUwtd01BXNgE2WFggJlhluf3By+T4ww4y5k5Z+D+XBdXzZlznnOf48C551lFjDEGQgghhJBmRsx3AIQQQgghfKAkiBBCCCHNEiVBhBBCCGmWKAkihBBCSLNESRAhhBBCmiVKggghhBDSLFESRAghhJBmiZIgQgghhDRLdkuCsrOz7VU0IYQQQojNXLkohDGGoqIiqFQq/bZp06bhxIkTXBRPCCGEEMI5m5OgjRs3Yv78+ZDL5bhzBQ6RSGRr0YQQQgghdmNzc9j777+P3bt3o6amBjqdTv8zdOhQLuIjhBBCCLELka0LqI4aNQoHDhww2q5Wq+Hm5mZL0YQQQgghdmNzEvTZZ5+hY8eOGD16tMH2SZMm4eeff7YpOC5QsxwhhBDSdNmSxticBLVt2xb5+fnw8vKCv7+/PqDCwkLU1NTYUjQnKAkihBBCmi5b0hibO0b7+/sjOTnZYBtjDAsWLLC6zNOnT2PatGl46623MHPmzAb3++677/DJJ59AJBJh2LBh+PjjjxtMemzM9ZyWSCSia2+G6Nrp2psbuvbmd+1cVHLYnAStWLECw4YNM9r+1VdfWVXezz//jG3btsHPz6/R/S5fvozExERcvnwZgYGBuP/++/Hll1/i+eeft+q8hBBCCGlebB4dNnLkSCiVSmzbtg0rV67Ejz/+CIVCgbi4OKvKi4mJwebNm+Hr69voft988w0efPBBBAcHQywWY/bs2fj666+tOichhBBCmh+ba4KuX7+OkSNHoqSkBMHBwfr//vbbb+jYsaPF5UVERJi13+nTp/Hwww/rX3ft2hVXrlyBQqGAp6enxeclhBBCSPNic03QggULsGTJEsjlcmRnZ0Mul+P//u//8PLLL3MRX4MKCwsNmsz8/f3BGENJSYldz+tslixZwncIvKFrb57o2psnunZiDZuToKqqKsyePRticV1RYrEYc+bMccjIMFOdohrqHCYSiQx+kpKS7BydMDSX6zSFrr15omtvnujam76kpCSD5zgXbE6CamtroVAoDLZVVVVBqVTaWnSjQkNDIZPJ9K9lMhlEIhFCQkJM7s8YM/hpLh8aQgghpClISkoyeI5zweYk6JFHHkFcXByWL1+O5ORkfPTRRxg4cCAeffRRLuJrUExMDNLT0/WvU1NT0a1bN+oPRHh39y9nVVUVVqxYgc6dO8PLywvnz583eF+lUkEkEmHHjh2ODJOQZoUxhurq6mY5lJw0zOaO0QsXLoS7uzu+/vprZGdnIzIyEgkJCXjxxRe5iE+vpKQEkyZNwo4dOxAYGIhnnnkGDzzwAEpLSxEQEIDk5GQkJCRwek7S/JSWlgIAgoKCLD72jz/+wJAhQxp8f+rUqQgODkbfvn0B1C0+3LFjRwwYMACDBg3Cjh07MGHCBOsCJ4QDWq0WGo0GaWlpyMrKQm5uLiorK9GzZ0+MHz/eYXGo1WooFArMnTsXERER0Gg0eOyxxxAXF6fvemGJY8eO6deznDdvHo0kJv9gdlJSUmLVcWfOnGHDhg1jfn5+rFOnTmzSpEmMMcZycnJYcHAwu3Xrln7fTZs2sb59+7L+/fuzV155hel0OqPyADA7XqbZbt68yQCwH374ge9QSANWrFih/7wsW7aswf1++eUXJpFIWEBAANu7dy+7ffu2/rjt27czpVLJioqK2PTp09lbb73F0tPTDY7X6XTs3XffZV26dGEA2Ny5c5lOp2OjR4+29yWSJkomk7G1a9cyjUZj8bFZWVlsx44dbPDgwfrPMQA2b948NmDAAINt3bt3Z9OmTWPr1q1rtMyamhqmUqn08RQWFrKUlBS2YMEC5u/vz7Zv387++OMPdvz4cbZmzRq2aNEi5u3tzeLj4xkAFhwczIKCgtgjjzzCXn31VYM4Bg0axH766SemVqvveW1VVVXsueeeYzExMez8+fNMo9GwkJAQduHCBYvvU1Mlk8nYF198wZ599lm2bNkydvPmTb5DMhsXz3e7ZQfDhw+3V9EWEUISpNPp9L+8AMz65SWOlZ2dzQICAlhlZSXTarXMx8eHvf/++wb76HQ69sorr+g/U97e3vr/f/nll22O4cEHHzSZyBNyt5KSErZ79242fvx4gyQlMjKSKRQKs8rIyclho0ePZgBYp06d2NixY9nGjRtZcXGx0b46nY6pVCq2adMmtnDhQtazZ0/WunVrtnfvXqbT6ZhOp2O5ublMpVKxDz74gIWEhOjLrY9tyJAhLC4ujslkMjZ48GAWGhrKJk+ezN588002d+5cNmnSJLZp0yaWmZnJamtrTcaclZXF3nvvPRYZGcmkUimrqqoyud+BAwfYY489xgCwsWPHspycHP17mZmZ7LHHHjPrHjU1Op2O/fDDD2z79u3syJEjLDExkY0cOZL98MMP7MaNGyw5OZnFx8ez119/3eTn4G6lpaW8Ps94S4ImTpzIKisrGWOMiUQiJhaLDX7qtwmBEJKg6dOn62OYPHkyGzlypMH7eXl5fITFGGNs8+bNLCUlxWDb+vXr2RdffMGUSiVPUTmWUqlk48ePZ3///bd+W2VlJWvbti1bvnw50+l0rLy8nI0ZM4aFh4ez69ev2yWO6dOns7KyMruUTZqO/Px8BoCNHj2aTZw4kV26dEmfiLz33nv3rHFOS0tjPXv2ZGFhYWzNmjUGn3tzaTQalpiYqP/76uLiwgIDA/Wvc3Jy2Llz51hRURErKirSPy+4olQq2fbt21l4eDgDwE6ePMmUSiU7e/YsW7NmDYuKimJffvkl27hxo8njJ0+ebHVrhTPbv38/mzVrFvvmm2/Yhx9+yA4ePGj0xUun07Gff/6ZDRkyhH377bdG71+9epXNnDmTDRgwgE2dOpUNHz6cbd682ZGXocdbErR161b9jYmLi2NZWVkGP5mZmSwuLs6mwLjiyCRIq9Wa3N6yZUtWVFTEGKt7uNb/kWCsrprYEfHV1tayGTNmsNWrV+u3bdiwgc2cOZPNnj2bLV26lGm1WvbOO++wt956i61bt47NmzfPZFkajYZVVFTYPWZ7Ki0tZc899xxbtWoVA8DWrl1rtI9Wq2Wurq76z1BsbCyrrq62W0yJiYksNTXVbuWTpmHOnDns4MGDDb5/7do1fWKSnZ2t315QUMDWrFnDunXrxn744QdOfod1Oh2rqKhgOTk5+qYvR9YMqNVqdvToUfbQQw+xmJgYdt9997GlS5feM8FZtWoV27t3r4OiFI4ZM2aw/Px8s/bVaDTs/fffZ927d2cPPPAAe/DBB9nYsWPZ+PHj2e+//67PARQKxT0/k/YiiOawkydPWrTd0RyVBG3cuJEB0Ffj3r59m+3evZutWrWKde/e3WDfTz75hAFgxcXFbM+ePQyA3ZtB3nvvPfbf//6Xvf7662zlypXs999/Z5MmTWIajYbpdDr2xRdfsBEjRrBPP/1UH8u0adPYtWvXDMrJz89n8fHxbNy4ceyll15iR48eZf/5z3/Yhx9+yGpqaux6DVzR6XRs7NixLDExkT311FPso48+avD+63Q6VlZWxuRyud3j+vjjj9nhw4ftfh7ivJRKJRszZsw999NqtWzhwoXMy8vLoLnsq6++sqrfUFNz/Phx9s477/AdhsOZ89m5m1wuZ7dv32YlJSUNfnZkMhkbMWJEgxUB9iKIJKigoIAdPXqUabVaplQq2XvvvceWLFnSYFutozkqCQLAXnrpJX0b+Z1/eOprge504sQJNn78ePbOO++wyMhIzquLGav78M6ePZtt2LCBPfTQQ/oq8w8++IDNmzePlZeXN3p8Wloae/LJJ/WvNRoNGzt2rL624vDhw2zJkiVs8+bNbO3atezRRx91+C+BNVJSUthrr73GdxhGkpOT2ffff893GETAUlJS2NKlS83eX6fTsVOnTrFr166Z/DvUXNXU1LCJEyfyHYZDVVZWsilTptit/Pfff9/htWtcPN85GSIfGBiIQYMGYenSpdi6dSvCw8Nx8+ZNbNy40dbiBUun0+mHaiqVSoSEhGDBggV49913UV1djT59+uDcuXMNHj9w4EBcv34du3btwhNPPAGVSnXPRWMt9euvvyI6Ohq5ublYu3atfobNN954w6zjO3fuDA8PD1y4cAG9e/fG559/jvHjx6NLly4AgOHDh2P48OH6/YuLi7Fv3z48+OCDnF4H17Zt24YXXniB7zCM+Pv7o6ioiO8wiICdO3cOMTExZu8vEoks2r+58PT0hFKpBGOMs5mHhe7q1avo2rWr3cqfPXs2EhMTMWbMGLudwx5snizx9u3bWLNmDVxdXZGcnIzdu3cjJSUFV69e5SI+wRo8eDAOHDiAffv2Yfr06Vi4cCHatGmDrVu3YvLkydi3b989y/j111/x5JNPws/Pzy4zbF+4cAEPP/ww3nzzTbRs2dKqMt5991288sor2Lp1Kw4fPoxnn322wX0TEhLwxRdf6Ccj02q1OH/+vOAmJ7t+/To6d+7MdxhG/P39DWZBJ+Ruly9fRvfu3fkOo0mIiIhAXl4e32E4TFpamv4LrD2EhYWhtLQUWq3WbuewB5uTII1GA6BuMqrWrVvrV44PCAiwtWhBu3TpEs6cOYOxY8dCo9Hg+eefBwB8/fXXGDVqFEJDQ+9ZRocOHbBx40ZIJBKjpUe4cP36dbRv396mMlq1aoVVq1YhNzcXW7ZsafRbU1BQEPr06YN9+/aBMYYFCxZg4cKF2L59u00xcCkvLw9hYWF8h2ESJUHkXvLz863+QkMMdevWDampqXyH4TDZ2dmIioqy6zn69euHs2fP2vUcXLM5CQoICMCzzz6L5557DnPmzAEAnD17Vp8cNVUxMTFYvHgxli1bhh07dsDb2xsAMGPGDOzfv9+isiQSiV1qgmprayGRSGwup1evXli4cCG8vLzuue+iRYuwfPlyzJ8/H+Hh4fjpp5+wdetWm2Pgyl9//YXBgwfzHYZJlASRxuh0Ok4XjmzuunXrhitXrvAdhsPk5uYiPDzcrucYM2aMWa0gQmJzErR+/Xr4+/tjxowZmDNnDgoKCvD555/j1Vdf5SI+wfLx8cG4ceMwe/Zsm8uyVxLEBz8/P2zatAkTJkzAq6++Cj8/P9TU1PAdlt6VK1fQo0cPvsMwiZIg0pjc3FxERkbyHUaT0bVrV6SlpfEdhsMUFBTYvRZ8wIABOH36tF3PwTWbO0YHBgbiww8/1L8OCwvDt99+a2uxTmHXrl2clGOvJIivb4zh4eEG3zjqm/uEsLjt5cuX8dJLL/Edhkm+vr6orKzkOwwiUHl5eYiIiOA7jCajVatWuH37Nt9hOIxGo4Grq82P/Ea5u7tDp9MZDBwSOpuj3Lt3L2bPng2lUom0tDRER0cjKioKJ0+e5CI+weIywbBHEqTVagXzIWzXrh0yMzP5DgMAIJfLIZVK+Q7DJLFYLLhO5EQ4qD8Qt6hZ0T6E9PfeHDY/JVetWoX58+dDIpHg7bffxgMPPIAlS5YgMTGRi/iaBU9PT847RstkMvj7+3NaprWio6Nx8+ZNvsNAbW0t3N3d+Q6DEKtQEsQ9sVjsdKOZrKHVauHi4uKQczlbXyubkyDGGPr06QO5XI6UlBSsWrUKs2bNgpubGxfxNQv2qAkqLy8XzAi99u3b48aNG3yHgWvXrulHLwoVfTslDaEkiHstW7ZEfn4+32HYXXl5OQIDAx1yru7du+Py5csOORcXbE6CKioqUFtbi++++w6jR4/W9/twVNbZFFAS5BhpaWl2nSyMEHuiJIh7kZGRuHXrFt9h2F1JSQmCgoIcci5nqwmyuZfUlClT0KJFC9TW1uLYsWOorKzE4sWL7T4Uj09c97eRSCQoKyvjrDwAKCsrE0wS1Lp1a0H8ocnLy6OJ5ojTKi4uRnBwMN9hNCmRkZHIycnBwIED+Q7FrkpLSx322QkICEB5eblDzsUFm5OgRYsW4eGHH4a3tzciIiKgUCgwefJku85MyTeVSgUPDw/OyvPw8EBtbS1n5QHCqglydXWFTqfjOwzcvn0bo0aN4jsMQqyi0+mohp1jkZGROHXqFN9h2J0ja4KAumZ9Z1mShJPqjE6dOumHbnp6emLYsGFYsWIFF0ULklKp5GQSwnru7u5QqVSclQcIKwkSitu3bwu+OUEkEgkiYSSkOWjdujVycnL4DsPuHFkTBAChoaEoLi522PlsYVVN0E8//YRJkyZBJBI1OFngvn37mmwixHVNkLu7u11qgnr16sVpmc5OSCPmGuLl5YWamhr4+PjwHQoREI1GQ7VAdhAeHo7c3Fy+w7C7kpISREdHO+x89f1AzVk+im9W1QStX78eVVVVAIAjR46gTZs2Rj9c1pQIDTWHWYfvOXCcoXrW29sb1dXVfIdBBKa4uNgpHijOxsPDA2q1mu8w7M7RNUFCmRbFHFbVBN05U/LTTz+NJUuWGO3D9wPPnpylJkhISVD9w53PGg6hJ0BA3X2qqqpCixYt+A6FCAiNDCO2cHSfoOjoaM5WVLA3m/sEJSUlmdzelL+1cD3pXnNIgvgeMeAsSbmPjw/VBBEjlATZj5eXV5P/nSstLXVoEtS+ffumXRN0J51Ohx07duDatWsGnXuTk5Px3HPP2Vq8INkjCeK6YzTXtVW28vf3R0VFBVq3bs3L+YW8XMadKAkiplASZD/1U3h07tyZ71DsRqPROHQCY76/9FrC5iRo2rRpSE1NRa9evQxusrWT/ymVSiQkJODq1avQaDT44IMPTA5rTk5OxurVqw06ui5duhT33XefVee1hFqtFnxNkNBIpVJeFwctLS112IyptvDx8YFcLuc7DCIwBQUFNMeVnTSHJIg0zOYkKD09HRcvXjSaPLChZrJ7SUpKAmMMf/75J65du4a4uDikpaWZ7COxevVqxMfHW3UeW3BdE2SPjtFC4+fnx2sSVFZW5tDqYGv5+PjoBx0QUq+wsLBJdzHgk1Amc21qpFIpKioq4Ofnx3cojbK5T1C3bt1Mbu/Tp4/FZel0Oqxbtw5z5swBAHTs2BF9+vTBpk2bbIqRa87QJ0ho6n8h+FJWVuYUNUG+vr5UE0SMFBYWUmd5O2nqSRBf/SGdpV+QzTVBH330EZ555hn07t3boGlq2bJlmDBhgkVlZWRkoLS01KBasmvXrjhz5ozJ/deuXYslS5ZAJBJh6tSpDuuDJPQkSIhzikilUhQVFfF2fmdqDuPzPhFhqqmpgbe3N99hNElNPQmqrKzkpTamfph83759HX5uS9icBK1YsQK7du3C9evX4er6T3EFBQUWl1VYWAgABv9g/v7+SE1NNdq3RYsWGDVqFGbOnIni4mLEx8dDq9Vi/vz5VlyFZbhOgtzc3DjtGC2TyQQ1MgyoS4L4XERVJpPpZzUXMmoOI8SxWrRooX/2NEV8fQGMjo7G8ePHHX5eS9ncHPbbb78hJycHx44dw5EjR/Q/U6dOtbrMu+dzMVWdN3bsWMyaNQsikQihoaGYP38+vvzyy0bLvPPH2j5LAPdJkIuLC6dLJQhteDzAf8doZ5gtGqDmMEIcTSwWO80UGtbgqz9k/azRXEpKSjJ4jnPB5iSoc+fOJmeHnjVrlsVl1Xf8k8lk+m0ymcysDoGRkZHIzs5u8H3GmMGPkJIggNuJ/ISYBPn5+fHaJ8hZkiCqCSJ3q62tdejwZtK08FUTFBYWZlWLUGPqB07V/3DB5iRo2LBhmD17Nnbu3Inff/9d//Pss89aXFZ0dDQCAwORnp6u35aamoqYmBijfZcvX27wurCwEK1atbL8AqxgjySIS2VlZYJLgviuCRJiYmiKr68vJUHEQFFREY0Ms7OmPDUFXzVB9SvJC53NfYJef/11hIWF4ciRIwbbrWljFYvFmDt3LtavX48hQ4bg+vXruHDhAr777jukpaXhhRdewIEDB+Di4oK9e/diwIABGDZsGGpqarB27VrMmDHD1ssxi9CTICE+8PlOgpypJqip/jEm1qHh8fYXERGB3NxcdOnShe9QOFdaWurQxVPv5OHhIbiJe+9mcxI0ePBgHDx40Gj7xIkTrSovKSkJCQkJiIuLg0ajwZYtWxAWFoasrCxcvXoVarUaLi4uSExMxLvvvgudToeqqircf//9eOONN2y9HLM4QxLUpk0bvsMw4OHhYfUEmlyorq6Gl5cXb+c3V/0q8oTUKyoqouHxdtaUk6CysjLExsbycu6oqChkZWWhU6dOvJzfHDYnQaYSIAD45ZdfrCpPIpEgOTnZaHtcXBzy8vL0r8ePH4/x48dbdQ5bCT0JKisrE9ywRCEsXiqEGO7FGWIkjkU1QfbXunVr5Obm8h2GXfA5R1q7du2QkZEh6CTI5j5Bp06dwrvvvova2lrk5eVh5MiRiI+PN+jX09Q4QxLkDHPiENOcoR2dOM6tW7d4W3OvuaivCWqK+JwjrX6uICGzOQl655130KNHD7i5uWHx4sXw9vbGsGHDmuziqYDwR2tQEkRI05GVlYW2bdvyHUaT1pSTID77Q9bXBAmZzc1htbW1mDRpElQqFXbt2oUbN27A39/fIQuZ8kXoNUHl5eVO0QnYUahmhTgzWjLD/sLCwpCfn893GHbBGDNa29NR6vsECZnNd6Z+OO8vv/yC++67T//w9fT0tLVowRJ6EqTT6Qxm7xYKsVgMrVbr8PMqFAqn6BRNiCmMMeorZmcuLi68/G1q6upHhwmZzU/K++67D127dkVeXh5+/vlnKBQKfPnll/Dx8eEiPkFSq9WCToKESiqVQi6XO7yWylmGx9cTiUTQ6XS8fXsjwqHVaulz4CCurq5Qq9WC7urgrIScyHOydtjIkSMREBCA2NhYVFZWIjg42Ggyw6ZE6DVBQm3+qZ8riJKgxnl7e6OmpqZJf5Eg5snLy0N4eDjfYTQLbdq0QU5ODm9z6tiDRqPhPYlu0aKFoKd54OTujB49Wj8PgVQqxdNPP43Dhw9zUbQg2SMJ4qqpSKgJEMDfhIlCnDyyMbR0BqlHnaIdxx5rXfGtuLiY9+kVhD5CzOaaIK1Wix07duD69esGbX/JycmYN2+ercULkj2SIHd3d/1EkLaoqqqCVCrlKCpu8ZUEOVtNEC2dQeplZWUhKiqK7zCahfokaPTo0XyHwpn8/Hy0bNmS1xjqR4gNGjSI1zgaYnMS9MQTT+DSpUvo3bu3wUKqfM4ObG/2SoJqa2tNLkZrCSEPj+drEVVnGy1HS2eQellZWU3qoSxk7du3x2+//cZ3GJwqKCgQRBK0b98+XmNojM1JUFpaGi5fvmxUg/Hee+/ZWrRg2SsJ4qIXvZCTIL5qgoqLi9GzZ0+Hn9da1BxG6mVmZlJNkINERUUhMzOT7zA4JYSaIKE3h9ncJ6hDhw4mm3BGjRpla9GCpdFobG62upuHhwdqa2ttLqe4uBjBwcEcRMQ9vpKgoqIihISEOPy81vL19aWaIAKAVpB3pPouCU1Jfn4+wsLCeI0hKCgIxcXFvMbQGJtrgoYOHYqnn34ajzzyiEHn05dffhnnzp2ztXjB4nq4n0QigUKhsLmcgoIC3j/0DZFKpbx0PCwuLnaqJEgqlfLSbEiESahDi5sisVgMjUYjyHnWrCGE5jCRSKRv6RDiavI2/0u//vrrCAsLw++//26wvbCw0NaiBcseI7C4+vZfUFCAuLg4DiLiHl8P95KSEsHWjpkSHByM69ev8x0G4Zk9apxJ4zp16oT09HR069aN71A4IZQvxT169MDly5fRr18/vkMxYnNz2ODBg5GZmWn005Sbw+yBq6YioXzoTfHz8+OlOUyj0TjVBGghISGCrj4mjpGXl4eIiAi+w2hWevfujQsXLvAdBmcUCoXNg2240LdvX5w9e5bvMEyyOQk6ePCgye1r1qyxtehmhcuaIKEmQXz1CXI2lAQRoK5TNM0R5Fh9+/ZtMt04dDqdYJpShXxfOWv4LC0tRXV1tf71tGnTcOLECa6KFxR7fLC4SoIqKyvh6+vLQUTc4yMJEvLkkQ2hJIgAwIkTJwTZfNCU1TeHNQVCqkls06aNYBdStTkJOn/+PKZOnYqMjAyDB45QMlCuMcbs8mCVSqWcrWIs1Hvv4eHBSedvS5SWljpVfyCgbvFhR98nIjynT5/Gyy+/zHcYzYqLiwt0Op2g17oyV3p6Ojp16sR3GADqnkne3t6oqKiAn58f3+EYsLk57PXXX8eBAwcwdOhQ6HQ6KJVKbNu2DW+//TYX8QmOvdYN8/X15aSWRMg1H3z8UXHWeVaE/O9I7E+j0UCpVMLLy4vvUJqdLl264OLFi3yHYbPU1FR07tyZ7zD0HnjgAUFORmlzEsQYQ9u2bfV/tN3d3fHoo4/i/PnzNgcnRNXV1XZZ2JKL5jCZTCa4LPtuYrEYOp3OYedz1mUH3NzcOJk3ijinv/76CwMGDOA7jGbpySefxIYNG/gOw2a///67oJaqGDduHPbs2cN3GEZsToI0Gk1dQWKxvvf3tWvXcOnSJVuLFqSSkhIEBQVxXq6/vz/Ky8ttKsMZOlKGhYU5dPoEZ02CoqKikJ2dzXcYhCf79++n5TJ40rt3b1y+fNmpa2PVajVqamoE9aW4devWyMvLc+iXYHPYnATFxMRg5cqVmDdvHuLi4tC6dWt0794dU6ZM4SI+wbHXqrzh4eHIy8uzqQxnSIIiIyORk5PjsPNduXJFMO3ilmjfvj2uXbvGdxiEB1qtFn/88QdiY2P5DqVZEolEGDBgAI4dO8Z3KFY7efIkBg4cyHcYRuLi4nD06FG+wzBgc8fojz76SP//EREROHnyJDp37ozx48fbWrQgFRYW2iUJkkgkNi86+/fffwv+22ObNm2QkZHhsKp+IaydY43Y2Fjs3r0b48aN4zsU4mB79+7FmDFjaKJEHs2ePRsvvooZmPkAACAASURBVPgihgwZ4pQdpPfv34+JEyfyHYaRuXPnYuHChRg+fDjfoejZXBPUrl07PPfccwDqJk5cuHChzQmQUqnEzJkzERcXh/79++PAgQMN7rty5Ur069cP/fr1w4oVK2w6rzlu3ryJdu3a2aVsW6tfL126ZLRQaFJSkk1lcq1///44ffq0Q861YMECtG7d2iHn4lqfPn1smrRNaP/ujuTM137mzBk89NBDmD17tlXHO/O124rLa2/Xrh169+6N3bt3c1amPd197efOnUPfvn35CaYREREREIvFDm0NuCdmo169etlahJHXXnuNPfXUU4wxxtLT01lAQAArKCgw2m/v3r2sU6dOTKFQMIVCwTp16sR27dplsA8AxsFl6s2ZM4fdvn2bs/LuNGvWLJabm2vVsRqNho0ePdpoO5fXzgWtVmsyTnsAwL7//nuHnMsexo0bx6qqqqw6Vmj/7o7kjNeu1WrZt99+ywCwr776yupynPHaucL1tcvlcta1a1f2xx9/cFquPdx57RkZGWzGjBk8RtO4lJQU9vjjjzOdTmdzWVw8322uCerUqZPJlXdfeeUVq8rT6XRYt24d5syZAwDo2LEj+vTpg02bNhntu3btWjzxxBOQSCSQSCSYPn06vv76a6vOaw6tVovMzEy7zcg8dOhQHDp0yKpjDx06JKiRAA0Ri8UICwuz+0KqMpkMAARZJWyup59+GqtXr+Y7DGJH+fn52LBhAx566CHMmjULp0+fRkJCAt9hEQA+Pj74/fff8fXXX+Nf//oXSktL+Q7pnv7++2/Ex8fjjTfe4DuUBg0bNgxdu3bFggULcPv2bb7Dgeh/2ZTVli5dit9++81oFflly5YhNTXV4vJu3LiBDh06GPS9mT9/PkpKSrBlyxaDfSMiIvDJJ59g6tSpAIDt27fjhRdeMJh0sL4918bLBADs3r0bFy5cwOLFi20uy5SqqiqMGzcOGzZssGhE061bt/DUU09h27ZtRhMDikQiwY1ySE1NxaJFi7B582bORy8olUrs3r0bn3/+OVJSUgR37ZZgjGHWrFlo164dXnvtNYtWYBbiv7ujCO3aGWMoLS1FVlYW8vLycO7cOVy4cAE7d+7EkCFD0Lp1a0yaNAmTJ0+2+VxCu3ZHsue1P/PMM7h8+TKioqLQp08fVFdXY/jw4ejfvz+vM/RXVFTA09MTHh4eSEhIwKVLl7BlyxZERkbyFpM5GGM4dOgQ1q1bB3d3d8yZMwcdOnRAWFgYxGLz62a4eL7bnAQFBASgd+/eRtsvXryIsrIyi8s7fvw47rvvPiiVSv0f/bfffhsnTpwwqiVxc3PDrl279J2BDx06hFGjRkGr1er3ccZObYQQQggxjy1pjFWjw959910AdeuDjRs3zmRTla1VuncnLw1dJCU5hBBCCLGGVUnQzz//jDVr1iAoKMhkAgTA6r459U1gMpnM5P/fvW9934/6/UJCQkyWS1XEzQ9dO117c0PXTtfenHBRCWJVx2h/f38MHToUgYGBDe5z/fp1qwKKjo5GYGCgwUq+qampiImJMdo3JibGrP0IIYQQQu5mVRJkTvY1b948a4qGWCzG3LlzsX79egB1ydSFCxcwffp0pKWl4f7779f3+UlISMCWLVugVCqhVCqxefNmGllBCCGEELNY1RyWkpJi19lMk5KSkJCQgLi4OGg0GmzZsgVhYWHIysrC1atXoVar4eLigjFjxuDKlSsYPHgwgLpZPmmGXUIIIYSYw6okqFevXo3OX8IYw4IFC6wOSiKRIDk52Wh7XFyc0fpaiYmJSExMtPpcTd2SJUv4DoE3dO3NE11780TXTqxh1RD5ESNG4PDhw43us2HDBjz99NNWB8YVLucJIoQQQogwcPF8t3nG6IYIIQEihBBCCGmIVUlQfn4+RowYgb/++ovreAghhBBCHMLmGaOFjprDCCGEkKZH0M1hhBAiRGvXrkV8fDy6d++OiIgIxMfH4/7774dGo7HL+U6ePIlevXrhxx9/NNgeFxeH+Ph4xMfHY8WKFY3uu3nzZrz66qt49tlncerUKZPnaWyfEydOoHfv3pg6dSpKSkqwefNmxMTEYOfOnQ3GnZWVhZEjR1pzyXbx119/wc3NDSUlJUbv2XLf6iUnJyM+Pp6TWFNSUtC3b1/8+9//Nrj3SUlJSEpKwogRI+5ZxkMPPWRwrEwmQ0FBASZOnIjY2FicOXMGV69excKFC22K25x7ZOv95fLecs2q0WGEEOKs5s2bh3nz5mHNmjWora3FokWL7Hq+7Oxs9OnTx2j7mDFjkJSUdM99Kyoq8Mknn+DMmTOQy+WIj4/H2bNnLdpn0KBB6N27N8aPHw+NRoMTJ07g8OHDvC7+aQmFQoFNmzYhPDzc5PvW3rc7zZw50+SoZGv17NkT//rXvwBAf+/rF8nds2dPo8fevHkT0dHRBv9u/v7+AICJEyfizJkz6N+/PwDghRdewMyZM62K0dx7ZOv95frecomSIEJIs3Tp0iVMnTrVqmOVSiUmTpxotH3cuHGYP3++wbbHH38c+/btM3n+5cuXQy6XY9asWWjXrp3Jff/66y906dIFAODr6wuNRoP8/Hy0bNnSon0AICcnB8uWLcOKFSvg6elp1rUyxpCQkIArV65Ap9Nh9erViI2NhVarxYwZM5CXl4eBAwfixx9/xI0bN/TH/fDDD1i1ahV8fHwQEhKC2bNn62uWLLl/ALBs2TIsWrQIv/76q8kYrb1vOp0Oc+bMQVBQENq1a2dwzQsXLkRgYCByc3ORkJCAbdu2Yd26dTh16hSmT5+OGTNm6JMcSyxbtgwJCQkmz9GrVy8AwC+//IIJEyZYXHY9c++vuZ8ba+4vH/fWGpQEEUKapUuXLmHp0qVWHSuRSEwmNpZ488030b9/f+Tn52PMmDG4ePGiyf1KSkrg4+Ojf+3r64uSkhKDB5U5+wDAO++8gw8//NDsBAgAduzYAbVajT/++AMZGRl4/PHHcerUKezYsQNSqRSbN2/G8ePHsXXrVv0xN2/exFdffYWjR48CAFq1aoVPP/1U/74l9+/o0aOIiopCZGSk2TED5t2THTt2wM3NDR9//DEyMjLwww8/AAB2794NmUyGlStXIicnB7NmzcKhQ4eQnZ2N7777DlOnTrX4IZ2cnIw//vgDZ86cQUJCQoPnAIBjx47hpZdeMjoWqFseqmPHjo2ey9z7a+7nxppjHXlvbUFJECFEELaduYXccoXN5UQEeGJK/9aN7qPT6VBUVNTgH/v09HTs3LkTsbGxOHXqFOLi4jBkyBD9+0qlEuPHjzc67qGHHjJ4eDWmvjmjZcuWqK2tRXFxsckFoIODg1FVVaV/LZfLERwcbPE+ALB8+XJ88cUXCAkJwZQpU8yKMz09HYMGDQIAtGvXDuXl5QDqljSqX6txwIABBsspfffdd5g7dy7c3d0BAOHh4QbxWHL/fvvtN/j4+GDZsmWoqKjAp59+iunTp6NTp06Nxm3OPbl27Rratm0LAIiKitJvT01Nxe3bt7Fs2TIwxhAUFAQAWLlyJaKiopCZmdnouU2ZOXMmJk+ejMzMTPj4+DR4jpKSEgQGBsLV1dXoWKAuITpz5kyj5zL3/pr7uTHlXsc68t7agpIgQogg3Ctx4dL169fRvn17/esff/wRKpUK6enpWLRoEY4cOQKNRoPLly9Do9Ggb9++BsdLJBIcPHjQ6vNfvXoVf/75J2bOnAmVSoWampoGF6QeMGAA3nzzTQB1DxoXFxd98paTk4PIyMhG97lTUFAQ9u3bhxEjRkAqlWL06NH3jLVTp07YuXMnnnnmGWRkZOj7prRv3x4HDx7EnDlzcPr0aYMROgqFQp8Ubdy4Ea1bG/7bWnL/7qyt+/rrr/Hiiy/qH7b112+KOfekY8eO2Lt3L4C6fi/1unXrhlu3buH1118HYwzbt28HUJfcrVq1Ci+++CK+//57s+K/W31i0NA5fv31V5MJjCXMvb/mfLasORbg595ag5IgQkizc+nSJfTs2VP/+scff8TMmTMRFhaGnJwcdO7c2WB/b29vq8/13//+F5cuXYJKpUKrVq0waNAgSKVS7Ny5E3l5ecjOzsann34KFxcXk/v6+fnhlVdeQWJiIqqqqvDVV18BqOuYOnr0aKSlpTW4T70///xTX+7w4cOxZ88ejB49GjNnzsTixYvRtm3bBr+BP/zww9i9ezfuu+8+aLVafPbZZwDqOuhu27YNw4YNQ0xMDDw8PPTHPPXUU5gyZQq2bNkCHx8ffY2RLT7++GN9Z9zXXnsNAPTXb8l9u/vafvnlF7zyyisIDQ1FdnY2jhw5gnHjxuHYsWN45513oFAoEB8fj6+++gp79uzBpk2b8MEHH+Dll1/G6tWrG713gOG979y5M7p37w4AJs8BAPv378c333xjdOwDDzwApVKJX3/9FXl5eTh37pxRcm4pcz5bgOnP8L3uLxf31hFoniBCSLO3bds2iEQiyOVyTJw4Ed9//z0qKyshlUpRWVmJV199le8Q7ebnn3/GqVOn8OGHH1p8rFqthpubG44fP44PP/wQu3btMtiuVqsxYcIErFixAt26deM6dN41dO9SUlKQnJxs1YiozZs344knnjB7/6ysLMycORMpKSkWn8vZcfF8pySIEEKasVu3bhk1V5nr0UcfRUlJCVQqFdauXasf3bR48WIcO3YMGo0G//rXv6wewi10Dd27ixcvYuXKlRg2bBjmzJljt/Onp6djy5YtkMvlWLlypd3OI1SUBJmBkiBCCCGk6aEZowkhhBBCrCTIJOj06dNo3779PdtTv/vuO/Tr1w/9+/dHYmIi1fYQQgghxGyCS4J+/vlnrFq1Cn5+fo3ud/nyZSQmJmL//v04deoUzp07hy+//NJBURJCCCHE2QkuCYqJicHmzZvvuabNN998gwcffBDBwcEQi8WYPXs2vv76awdFSQghhBBnJ7gkKCIiwqz9Tp8+bTCXR9euXXHlyhUoFLbPOEsIIYSQpk9wSZC5CgsLDZrM/P39wRhDSUkJj1ERQgghxFk4bRIEwGCtmnoNdY4WiUQGP0lJSXaOjhBCCCFcSUpKMniOc8Fpl80IDQ2FTCbTv5bJZBCJRCYXIARoniBCCCHEmSUlJRlUYHCRCDltTVBMTAzS09P1r1NTU9GtWzd4enryGBUhhBBCnIXTJEElJSUYMmQIysrKAADPPPMM9uzZg9LSUuh0OiQnJyMhIYHnKAkhhBDiLASXBJ09exbx8fG4cOECli1bhkceeQQAoFAocPXqVdTU1AAAunfvjo8//hijRo3CgAED0Lt3bzz33HN8hk4IIYQQJ0JrhxFCCCHE6dDaYYQQQgghVqIkiBAb3SiSG3wTySqpxjMbzkCro9pHQggRMmoOI8RMv18rRueWvpBK3HDiZgmO3yiFh6sYtRod4juFYnD7IGh0DEt3peKBrmFQa3UY3jmU77AJIaRJ4uL5TkkQIWbYcSEPSrUW1SotKpVq9GsTgPvaB0MkEoExhuX709Eh1Ad7Lxfg1dGd0D7UBx/sScPicV1xKrMMEjcx5EoNBrcP5vtSCCGkSeDi+e60kyUSYm8rD6RDpdGhpZ8ENbVaPD+8vcn9RCIRdDqGv/Mq8J+n+uu3d2vlh4/3p0Ot1eFkRimGdAimJIgQQgSEaoIIMWHbmVsI85NgSIcQpOVXokOoD1xdGu5CV1SphIerC/y83PTbdDqGzNJqRIf4QKtjWHPoOl55oKMjwieEkCaPRocRYgc3iqpQVl2LIR3qlmDp0lLaaAIEAKFSiUECBABisQjRIT4AABcxN+vcEEII4Q41hxHyP3KlGl8fvQkRRHhpZAe+wyGEEGJnlAQRAqBYrsKnh67j5ZEdEOTjwXc4hBBCHICSIEIAfHs8E6+P7QxvD/qVIISQ5oL6BJFmT6XRwt1V7JAEiDroE0KIcFASRKx2Oa8CBRVKvsOwGmMMpVUqnMuWoV+bALufz8vdBQq11u7nIYQQYh5KgojVtpzKwfenb/EdhtX+zqvAjG9O4VRmGWKiAu1+Pl+JK+RKjd3PQwghxDyUBBGr5Fco0LGFL1QardM28ZzOKsez8dG4ViiHxM3F7ufzlbihUqEGAMR9cMju5yOEENI4SoKIxUqrVPjvyWw80LUF2oX4IKOkmu+QrFJeXYuHerXC50/0ccj5pBJXVCo1UGm0KK1WOW3ySAghTQUlQcRi649nQipxQyt/T8R3CsHaozfx+7VivsOySJ5Mgf9NNqqfddTefCVukCvVKKmqRVSQN6pU1DRGCCF8oiSIWESuVMPD1QXPxkcDAIJ9PLB8ci8cv1ECnc45ajYqlWok/nAB8Z1CHHpe6f/6BFUq1IgM9ELF/5rGCCGE8IOSIGI2xhj2/J2PcT1bGr3XI8IP6YVyHqKy3K8Xb+OjR3uiXxv7d4a+k6/EDZVKNSoVarQO9IKshpIgQgjhEyVBxGwr9qfjr4wy/XpYd4qNCsTprDIeorJcXrkCbYK8HX7e+tFhlUoN1QQRQogAUBJEzHIpV4Y2QV745LHeJt8PlUqQW65wcFSWU2t1cHfl52Pv5e4CRa3WqCaIMYYiufPOt0QIIc6KkiBilmPXS/Bwr/BG92kd6IU/M0odFJF10gvk6Bzmy8u56ztgyxRqRAX9UxN0KbcCs5NP8xITIYQ0Z5QEEbMo1Vp4ujc+l874Hi2RervSQRFZ58ItGXq19uft/Ax1ncsjArwgU9QCAM7llMPf0523mAghpLmi1SLJPWl1zKxh5AHe7pDV1DogIusVVCgRJpXwdv7TmWUIlXrA090FKrUOQN18RQOjg6BUax0yaSMhhJA6VBNE7in5RBa6mNuEdEeylJYvjFqh17dfwm2ZAjeKqlCr1TlsXiBTlj3aA6cy/+lAzhjD1QI5WvpJnHodNkIIcUaUBJFGVSrVKKlSYVS3MIuOU2t1GLvmmJ2iskxmSTUyiqvxzq9XeOsPVC8y0AvDO4fqX2eUVGNM9zCE+3siTyb8juWEENKUUBJEGqTTMXx55CYm94uAi9i82hPR/44rlqt4bXa6U7CPB0qrVejT2h+P9I3gNRaRSIQPJvXQvz6bXY7+bQLRipIgQghxOEqCSIPO35KhZ4SfyXmBGhLi64GSKhUKK5VoIfWwY3TmmfntKQT7uKO0Sph9lXLLFWgd6IkwPwkKqTmMEEIcipIg0qCj14ox4o6mG3P0jwrAkfQiFFaqEBXsjZpa/tbH0ukYTmWWoXWgFworlZDcY3QbX0QiEdxcxFA7ybIjhBDSVFASREy6UVQFXw9Xi0crdQ6T4laZAoWVSnRv5YecshpoeXq4l1bX4rUxnTFzUBSKq1QI9uG/ZsocG09m4XxOOd9hEEJIk0dJEDFp29lbeGpQG6uOlSvVOJdTjq6tpBj36R+8TaBYJK9rknN1EaO0qhbBPgKci4f9kyBmllRDrlQjvVCOczkyHoMihJDmgZIgYuRaoRxRQd7wcLWu+eihXq1QoVCjTZAX+kb6I5+nvi4VNWr4e9UlPqXVKgR5C6smyMNNDJVGp3/97LBoDFp2GFKJG60rRgghDkBJEDGy/VwuHrVhFFX/qEAkz4pFRIAXNs4ZgMJKfpIgmUINfy83AIAIIrQJ8uIljob4eriiUvlPn6muraQ49MowJMRH8xgVIYQ0HzRjNDGQWVKNiAAvzhYZlbi5GNR22Mu7v6ZiZNdQ+Hq4oZW/BEE+HpDVqPXLUfzy/GCzh/k7CgPg5mIYU+gd0wowZt5M3YQQQqxDNUHEwNZTOZjM8Vw6XD7Gy6tND3UvrVahoEKJ17Zfwtnsuk7FMkWtviZIaAkQAEzoHY4X7+/Q4PuDlh3WzyL9y/k83mrUCCGkqaIkiOj9lVGKrq2k91wo1R4+3p9u1n5DVxyBSqM16jPTQipBQaUSbUO8odTocC6nHEq1TtBrcfl5ujU4Ys3L3QWhvh64WlCp72ietPMKGKNh9IQQwhVKgoje0WvFeLhXK7uV39gD/PMjNxp875kNpwEA7+9OhUqjw4UcGZ5ef8pgH4mbC5RqHQK93FFeXYtHvjzBaQ2Uo8W1C8KrYzrjRlEV/m/HFTwxIBLjerbEpdwKvkMjhJAmg5IgAqBuYkGxSGTXPigzvz2NnNIao+2NzSPEGMPBtCLU1Grwn2OZiAjwhEqjg4+HcXc2EepqV2Q1aohEgEZn/75I9tK7tT8Gtw/G+VsyDIwOQucwKR7o2gJ7/s7nOzRCCGkyKAkiAICrBXJ0bmmfxUXFIhG0OoaSKhVKq1VG71cpNfCVmO6jX12rRd9If5zLrps3JyLAC3IT+4tQ19FYLBZBpdEi3N8TBRXG53I2N4uq0CVMCgDwcHWhjtKEEMIhs5IgsVgMFxcXs36io2l4rzM6dr0Yg6KD7VK2r8QV3ZbsQ782AZArjZfRqFSq4efpZvLYCoUaPSP8cTa7HIsf7ILuraSoVKpNLjZaq9HB3UWE8ppadG/lh/wK51+QtFajQ9sQb/1rTzcXaLTOW8NFCCFCYtYQ+V69emH16tX33I8xhjfffNPmoIhjlVfXolKpRqC3fWZUDg/wxFvjuiImKhA3iqqM3q9QNJIE1ajRoYUPjqYXY859bXEyoxSVCjWGdAiGUq016PhcrdKghdQDJVW16N8mAFml1Xa5Hkf67ZVhBiPbWkg9UChXIdzfk8eoCCGkaTArCQoKCsKwYcPMKtDLS1gT0pF7O5BagMdjIu1W/uhuYQCA2zIF5ErjmZArFWpIJW4m58WpVKoRFeSN/xRmIMTXA17uLiiWq9AmyBuVCrVBEjS2exgCvN2x8+Jt9GsTgHnDnL9W8u6h/S38JCioUFISRAghHDCrOezgwYNmF2jJvqYolUrMnDkTcXFx6N+/Pw4cOGByP7lcjhkzZiAmJgaxsbGYOXMmqqqMaxnIveWVKxARYP+Hqq/EtcHmsBZSDyjUWqP36muJiuUqBPt6wNvDFVUqDfy93CC7Y5g8AzCofTC6tJSiUqFGizsmHWxK/DzdUGkikSSEEGI5wXWMTkpKAmMMf/75JzZv3ozHH38chYWFRvu99957yM7Oxp9//ok///wT2dnZWLp0KQ8RNw2O6HDr7e6K6loTSZBCg5b+nqhWNZwEVddq4evhCh+PukTK39MdshrTycBPzw1G68CmWSMpbSCRJIQQYjmrkqAJEyZwHQcAQKfTYd26dZgzZw4AoGPHjujTpw82bdpktO+VK1cQGxsLFxcXiMVixMbG4vz583aJqylTO7CTrVgsgqmpgiqVarTyk6BaZSpBUsPPyw0/zBsIkUgEHw9X1NRq62qCav6ZPfrOFK6h/kVNga/EDZW0uCohhHCi0T5BH3zwgdE2xhjS0tLsEkxGRgZKS0vRuXNn/bauXbvizJkzRvs++OCD+Oabb1BdXQ2RSITDhw9j7NixdomrKXtxy3nMGtzWYee7Owf6cE8abhRVYUr/1iZriapUGvi4uyK2bSAA6JvD/DzdkJZf6YCIhUUqcaOaIEII4UijNUGrV69GREQEwsPD9T8RERHw8DA91b+t6pu9/Pz89Nv8/f1RVFRktO/zzz+P+Ph4tGvXDlFRUejbty/+7//+zy5xNVXFchVi2wbqEww+pOZXorhKBV+Jq1FzWOrtSvyWWgjxHZ2DfTxcUa3SwM/LrcHmsKZM4iY22XeKEEKI5RpNgnr16oWRI0fi6aefNviJjY21a1B3908xtdzC+++/j4sXLyI7Oxs5OTnIyMhAcnJyo2Xe+ZOUlMRx1M7nxM0Su80N1JC7ex4FerujvKYWXu4uRjVBF3NluHLbsLanvibIx73uv80NTZZICGmukpKSDJ7jXGg0Cfrtt9/QqpXxWlLffPMNJye/W2hoKABAJpPpt8lkMv32O3366aeYN28eJBIJJBIJEhIS8PbbbzdYNmPM4IeSIOBaoRwdW/jwGkOYnwT5MqW+hudOlQo1Ns0ZYLDN28MFGi2r61/kyEAJIYTwqn7gVP0PFwQ1Oiw6OhqBgYFIT/9nRfHU1FTExMQY7VtbWws3t386wLq5uUEulzskzqbgVlkNfDzceK9ZkLi6YOOcAfDycEXNXc1h1bVaDG4fZLAt2NsD2xIGGpVDCREhhBBLWZQEJSYm2isOAHXLc8ydOxfr168HAFy/fh0XLlzA9OnTkZaWhvvvvx9abd2DcuTIkfj+++/1GeHmzZsxfPhwu8bXlOz5Ox+Px7R2+HnvTFYYY8gqrcbA6KAGm7fuTtLEYhFa/W+iwPp36hZ/tVPAhBBCmiyLkqCTJ0/aKw69+uquuLg4TJs2DVu2bEFYWBgqKipw9epVqNV1nWG/+OIL/X4DBgyARqPBv//9b7vH11TIFGoE2GmZDHNtPX0LPcLrOsF7ebigxsTosMbUJ1S1Wh08XF0a3bcp+uH0Lb5DIIQQp2bWshmOJJFITHZwjouLQ15env51aGgovv/+ewdG1nTkyRRo5cf/jMrpBXK8PLIDAMDNRQy11rBRy5zKncJKJUqqVPBwFVTLrt2ptTq8teMyxvVsCW8Pwf0aE0KIU2heTw4CANhzKR+ju4fxcm65Uo1nNpwGUDepob+XYW3UQ5/9YXZZIgC/nM/D2qMZ8HBrXh/l3ZfyMX1AJHLLFXyHQgghTqt5PTkIckproNJoEerLT01QZkm1yeUxAECh1hqsMm9uZ2eVRtusmsNcRCKk5lfioV6tcFtGSRAhhFjLoiSIqyFphD8p14owpb/jO0TX+3RaH/SPCsDRa8VG7xXLVQiVWjYRp0bHoNYySJpRTZCfZ910AuH+nrhdQUkQIYRYi7OO0TU1NTYHQ+yvRK7idYV1qaRuWoOn158yeq+0uhaBFnbWVml0YIxB0oxqgib1iUDiqE4I9vFAsVzFdziEEOK0OPv6PH78eK6KIvYkgBmHy6pr4e9lvMhpRU0t/C1Y/FTi5gKV895b/wAAIABJREFURgulWgeJW/NJgvy83BDo7Q4XsQhrj2ZYPKqOEEJIHYuToEOHDqFnz57w9PSEi4uLfhX3o0eP2iM+wpFbZTUmV2nng1Ktg6eJpKVCodavAG9O06uPxBVypQYqjbZZNYfd6ft5cTiYZry2HiGEkHuz+Mkxf/58LF68GH///TcyMjL0P/ZeT4zYZsnOK9h+LlefZPBpXM8w9I8yXrS1vEYNfy93lFap0PaNPXC9xwyIvh6uqFJqml1N0J16hPsZdCYnhBBiPosnGImIiMBjjz1mtH3Lli2cBETso3OYL1LSizF9QCTfoWBE5xa4eKvCaPvWf8Xh6LViXMyVIdzf854Jm4+HK+RKNZTq5lsTJBKJABqwQAghVrH4yfHII4/g4MGDRtsXLFjASUCEe7WauhmVbxRVoU2QF9/hAAC0OgaXu2p6urSUwsfDFRnF1RjTPQxSz8Zz9FCpBwoqVSiSqxDC05B/QgghzsvimqABAwZg8uTJkMlk8POrW/KAMYbCwkLOgyPcyCipQnSoN14d0wntgvldNb5eTa3W5LIdXu4uuFYoR+KoTvesCerYwhefTeuDW+U1gmjm440AOrsTQogzsjgJmjVrFp599ln06dMHrq51hzPGqCZIwNLyK9Ezwh/RIcJIgIC6ZR/cXIwf3l7uriioUEIqufdHU+LmgvahPmgfKpzrIoQQ4jwsToL8/f2xcOFCo+1fffUVJwER7mWW1ODhXuF8h2Ggoa4s3h4uUKi1RqvHk4bRnSKEEOtY3Cdo6NChuHr1qtH2DRs2cBIQ4R5jxv1v+CaC6WHwXu4uqBLIUH5nwUCzuRNCiDUsrgk6duwYVq1ahQ4dOsDf31+//cKFC1QbJEA6HRPk4KFR3cLQwsQSGd7/G/ZOzCdxE0Olab7TBBBCiLUsToIKCgrw+eefG2xjjCE/P5+zoAh3rhXJ0aGF8PrMDG4fbHJ7sI8HSqtrHRyNc/OVuEGu1FASRAghFrI4CVq8eDGefPJJo+131goR4UhJL8a0WP7nBjJXoJc7kh7uxncYTsX3f/MlhfjW1azV1Grg5W7xrzYhhDQ7FvcJMpUAAcDEiRNtDoZwr1qlcarh42KxCA/3asV3GE7FV+Jq0I9q6PIU6iNECCFmsDgJeuGFF7B+/XpcvHgRWq0WAKDRaBpdYZ7wh0ZZNX11M2f/kwSpNFpUKNQ8RkQIIc7B4iTIxcUF69atw6BBg+Dj44PY2FjMnj0b48aNs0d8xAaMMRo+3QzU9wmqF+rrgSK5iseICCHEOVicBK1ZswYnTpxAZWUlTp06hSlTpmD37t2YNWuWPeIjNtjzdwGCfY1HYJGmxVdS1yfo/pUpOJ9Tjlb+niisVPIdFiGECJ7Vq066uLigR48eWLRoEX766SdkZmZyGRexEWMMz28+h3bB3nyHQuysvk9QXLsg7LtSgK4tpSiqpJogQgi5F06W3h4yZAhSUlK4KIpwpLxGjScGRKJvZADfoRA78/ZwRaVCgxBfD1wvrEJ0qA9k1CeIEELuyeJxtL1790ZsbCxiYmLQv39/9OjRA0ePHkWbNm3sER+xUl65AvEdQ+DpTnPHNHVuLmKUVKnQPtQHxXIVokO8cex6Cd9hEUKI4FlcE7Ro0SJIpVJs3boVI0eOhIeHB0aNGoXw8HB89tlnSElJQUkJ/QHmU1ZJNXLKahAR4MV3KMRBCiqVaCGVwN1VjLbBPkazhB+/UYKlu1L5CY4QQgTK4pqg6dOnY/r06frX2dnZuHDhAs6fP4+DBw9ixYoVyMvL0w+fJ40rqlQiVCrhtMxPD13H4fQi/P7qcE7LJcJVWKlEC6kHfkwYaHJahPM55QD+WWOsoakTGGNo+8YeZC2j0Z6EkKbP5mll27RpgzZt2mDChAn6baWlpbYW2yz8lVGKZ/57BnOHtMOL93fgrNyIQC+8NqYzpBLnmSSR2KZWo0PbYO8Gk5tajQ5eHq7YfCoHVUoN5g2LNrlfYaUKrfwkUKq1tAwHIaTJ46Rj9N2CgoLsUWyTMzv5NDxcxdjzN3frrqk0WriJRU61VAax3a7598Hfy73B9xkAV7EIN4uqodbqGtwvs6QaA6ODDUaXabQ6aHU0AzUhpOmxSxJE7q2kSoWXR3bEmw92QfdwP87KvVVWg8gg6gvU3Li6GP4qqzQ6TPjiOACgokYNqcQN7UK8cau8BrXahhOatPxKDOsUgvwKhX7bN39kYvOpHPsETgghPKIkiCf/PZGFR/tF4JG+EYgM9Gr027klJnx+HN1acZdUEeeUVVKNrJJq6HQMF3Nl6Bnhh04tfKHSNP45K6+pRbdWUhTcMdlisVyFmjvWJiOEkKaCkiAHYozhcl6F/nWgd13zRZhUYvMMv5fzKnCrrAazBrdF+1Afm8oizq+4SoUHe4Th+M0SZJdWo22INzq08MWGWTH3PDZMKkF+hVLfBObmKuYsSSeEECFpNAk6d+4chg4dihEjRuDAgQP67ZMnT7Z7YE1RaXUtpv37z7qHyx0dWMMDPJFXrmjkyHvbduYWfksthFhMq4UR4MeEgXh9bBfsuHAbJVW1CPauWz5FJBI1uJ7cbZkCLf084e3hiuTjWVi+7yqAunmINNQniBDSBDWaBL388stYu3Yt/vOf/+DTTz/Fv//9bwA0+staOWU1eKRvOH6/XgyvOyYxbOXviTyZbUmQt4crCuVK+HjQiB5Sl+z4ebqhlZ8EDDArOT6TXY5+bepmGN88dwDcXMTIKa2xc6SEEMKfRpMgkUiELl26IDo6Gjt37sSJEyfw/vvvNzgMlzTuVlkNpg2IxHd/5qCVv6d+e0SAJ26V2ZYEubuKcausBi04nnOIOL+7f1sbqtM5k1WGji3qmlLbhfhAo2MYuuII1Fqd0eSLhBDSFNyzT1B5ed0ka2KxGMnJycjPz8eJEyfsHlhTU1GjxktbLyAqyBvRod7o09pf/15dc4NtfS4YA27LlAbJFSEebi5QasybuNTf083gC85TA9tg+aM9kWtjUy0hhAhVo0nQt99+a7Tt888/x/79++0WUFN1MVeG2LaBkLi54I2xXdA60HAYOxe1a4/HtEZ3GhlG7hAZ6GWySetWWQ1u39EEy5hhPzWgrpl2Sv8IrJzSy+5xEkIIHxpNgtq1a4eAAONVyIcNG2a3gJqq8ppafPhIjwbf56KB8fHYSFowlRhoF+KNbq2kBtvcxCL8eDYXG//M1m+rUKjh72k8w7hIJIK7q/ju/IiYkHq7Em/98jdOZ5XxHQohxEwWDZFPTEy0VxxNXoVCDT8TDxlC7KlbKz+8MMJwSRappxtUGh3Ud8wZVF6jRoB3w59P6hPUuOzSavx66Tbefbg7DqUVoaJGzXdIZktJL8JH+67i2PVi/dpyhDQXFiVBJ0+etFccTV5FTeNJkC1/eqpVGpuOJ82L1NMVBRUKeHn8s3Rg+f+3d+fhTVZpG8DvNGuXdN+gdIFCWxaBUpaySBcoO6J+KoOgoojiCDMsKijO0FEHHRUXdHBDFkWY0dFxYUTZiuwKIoJsLXShFLrvbZI2yfn+KI2UbmmbNmlz/66r10XenCTPyUuap+c95zmVVXB3bHzbDWqcTm/AhwfTsGR8GBwcJLh/ZDD++8sVa4dllnX7LqJKb8RTE8NRoTPg7b0XmQiRXWGxxA5SbRSQSxt/ux0kaPX+TP1XfQ8j67iQmVxVcpRq9TAYjZiz/kcAQHFlFdydzBup3HkmG9pq8yZb24N1SZcw/9ZeUMhqPt/d3R1xrURr88nE5z9fQf/ubpjQ3x8SiQSTBvhjdB9vvL47hb9PyG4wCbIRjnJpq75YhBDo5qZCObc1IDO5OspRWaXHqSslKNZUAQCKKqrh0cQGrA4SiSlJf/n7C0gvqOiQWG3df3+5gsGB7vUWOkQFe+DnjCIrRdW8XzOLUaKpRkyYT53jQ4I8MHmAP17ccQ46M1cVEnVmTIJshELmgKpm9nW62T+TLiK3TIf5t/ZC4m392yky6mrUKhm8nJWI8FdjaLAnhBAoqqxqMglyVkpRWVWTaAd7OiGDRRRRoqlGal4F4iJ86903rq8fvvst2wpRNS+vTIftp67iwdEhDd7ft5srHhgVgtX/O4cybeeZ20TUGrLmm/zO1od3OzO5tOX7M/3r2GWE+6nhrVa2U1TUFbmq5PB1VWLl1H7YcDANJZpqlGn1UKsa/3XgpJChQmeAWiW3yDYvnYUQAt+fycaJy8VQyhxM836Amk2Q7x8Z0uDjpA4SDA3xwMGUfIzp492q1y2oqElMpRbcCqdKb8Rbe1OwYnJEk2U5eng4YfH4MLzy/QWMCvVCfnkVMgsrIXWQoLevC+4c0sNiMRFZU4uSoKYmRldWVsLJyanR+82l1WqxYMECnD9/Hnq9HqtXr8aECRMabHv48GH85S9/QXV1NfLz8/HYY49h0aJFbY7B0sp1+mavsStkDqhqYRLk6aTAxbxy3BLA2kBkPg8nBYYE1ZS+8HdTmXaMb2prDWel1HTJtWZOUdcfISisqMI/ky4iLtwXz0zpi/PZpXj/QCoWxIQis7ASSrkDfJr4A2Rif3+s2ZkMF5UMg28ojnqz/cl5OHQxH0qZA6qNAjIHCYxCwMdFiYt55XhyYoRFVpaWaKrx+q5kPDymF5wUzf/q93BWYNX0/vgtqwS39HBHwPVCrF+cuIIDKXm4tY9PM89AZPtalAQ1Zdq0adi7d2+bnycxMRFCCBw9ehTJycmIjo7GuXPn4OfnV6ddWloaFi9ejG+//Rbe3t44c+YMNmzY0ObXbw/LPz+FK0UaPDExvNE2CmnDl8P2XcjFgAA3eLvU/2Xr5qRAWl5Fvev6RE1xVEgxfVB3AIC3ixL5ZVXNPsZFKUOFTg9ttQGOCil0Lbx0ayuMRoFXd16ASi7FovjejY6GpOaV4+OjGViSEAZXVU0CEuHvigvZZdhyNANnrpZi1fR+Tb6WRCLBsglh2HI0A3vO5SC6lxdG9647KrT3fA7yynR4ekrfBp+jqKIK7/5wCcsnRbSit78/x7+PZ6Jcq8efx/WBh7P5qwClDhIMuimBu3NID7yw/SwGBbqb3hugZpXcp8evIO96Ug2JBHHhPogMql9rjshWtDgJ2rNnD5YsWYKUlBRUVdX88hRCWKTisdFoxPr16/HFF18AAMLCwhAZGYktW7bUq1H02muv4cEHH4S3d80vlf79+2PNmjVtjqE9BLg7wqeBJOZGCpkDqg31R4s+PJiGh8b0RFx4/XkHbo5ypOVXwKsFv9SIbuTtosCpKyXNllhwVspQUaVHWn4Fgr2ckJxT3qLXqdDpcSAlD5MGdGt9sBbw3ZlsxEX4olpvxK6zOZjQ37/O/dUGI/59LBP55TqsnNIXsptWdM4YHIBfM4sxY3B3qOTNFyaVSCS47/ols38mXYS7kxz9r1d1v1aiwS+Xi7FsQuN/HHk4K9DL2xmnr5Tglh6/j/hqqgzYfuoqRoZ6oYfH7yPwQghsPJSO4soqSCQSCCGgUkhxV1QP+Kott6/gY7Gh+GfSRTw9uSZ5K6qowppdF/DQ6J7o5VOz/5zeYMQXJ7JwOqsE90UHc89Jskktnhi9aNEirFy5EqdPn0ZqaqrpZ/jw4W0OJjU1FQUFBYiI+P2vnn79+uH48eP12u7ZswdVVVWYMmUKRo8ejSeeeAJarbbNMbQHlVza7MTlxuYE+bmqkFPScL/cHeVIza9o0V92RDfyViuRX65rtl3NSJAByTll6OOrrlfhPDWvHF//erXRx6flV+CN3SltjLZ5hy7m48TlxldlncwsxrAQT4wM9cLx66u3jEaBd/Zdwupvz+H57WcxNMQDi8eH1UuAag0KdIda1fLLU4/FhOKrk1eRWVgJg1HgvR9S8Xhc72Yfd0dkAL48mWW6LYTAmp0XcEsPN2w6lI6MG1bqfXb8CiK6qbF0QjiWJIRh6YRw/DG2t0UTIADwclEiNswX/0y6iNNXSrBm1wUsnxRhSoAAQCZ1wD3DAhHmp8abe1I4p5RsUotHgnr06IGZM2fWO75t27Y2B5OTkwMAcHP7/S8ed3d3nD17tl7b9PR0vPvuuzhw4ADc3NwwY8YMLF68GO+++26b47AGuVTS4CUGf9ff52wAQHJOGTydFfB2UcLVUYbCCl2T9YeImqJWylCqbb68gpNCigqdHleKKjFpgH+9kaNfLhdjf0oebrt+me1mRZVV8OyAZP1ASj4UMgdkFFTg9sEBdUYfzlwtQbifGkDNCE2IlzNS88px6GI+bu3jjQHtPLfOwUGCJyaE48Ud5yAEcO+IILNGk2RSB0QFe+DH1AKM6OWFT368jEkD/BHh74rlkyPwwvaz+PP4MFwpqkSZTo9RoS2fiN0aI0O94OuqxJmrpfjrtP6mOkk3i+7lBWeFDG/sTsHi8X04IkQ2pcXfnnfeeSd2795d7/iSJUssEhBQfzPRhv6C0Ol0uPfee+Ht7Q25XI7Fixdj06ZNMDayG7tEIqnzk5iYaLF4LaHmclj92KUOkjqTql/flYyjqQUAAJmDA1jTjNqi9rPW3NeSi1KGcp0eVQYBpaz+F3dumQ5+ro2PNhRWVHXIiKVcKsHShDB0c3PEizvO1/nsfH3yqmkuFADMGNwda/ekoFSrb/cEqJZC5oBV0/sj8bb+CLuekJljUn9//O/0NWw+nA4nhRRDQzwB1IwgPzExHGv3pOBYehEeHBXSTpE3LNTHBbcN6t5oAlTrlh5uGBvmg38msSI1tV5iYmKd73FLaPFI0IgRI3DXXXehuLjYNGIjhDCN4rSFr2/NvJfi4uIG/30jDw+POpOle/ToAZ1Oh/z8/AbbW+uDV6HTw8mMTU0VUgdUVjVSnOyGk92apfRETanU6Zsts+B8fWK04fofGTf/+tFU6Zsc1SiqqOrQuWvRvbzgopThzT0pmD0iCJ/8eBkjennW+bJ2Vsrw4p0DoWzmC9wWODhIsHxSBC7llWNgj7oTldUqeaeoExYV7AG9wYgPDqTikbGh1g6HOqHExMQ6AxiWSIRanAQ9+OCDeOyxxxAZGQmZrObhQgiLjASFhobC09MTFy5cMCUyZ8+exZQpU+q1HTx4MHJzc0238/LyoFAo4OXl1eY42mrDwTSE+rpgTG9vnLhcBF8z6vjUJDdNLzte9dVvcHOU19mc8cRfEtocL9m35NxyDO/p2WQbR7kU6QWV6NutZvRCIXOAttpg1uUcACjR6Nt9A+Gb4xkQ4AYHiQRf/3oV948MhlcDixMczfgDxVY4K2X1EqDOZkQvLxRrqrHzTHa9SelE1tDiP4Hc3d3xxBNPYNy4cYiJiUFMTAxiY2PxzjvvtD0YBwfMnz/ftNQ9JSUFJ0+exOzZs3Hu3DmMGzcOBkPNaMn8+fOxbds2VFbWVK7dsGED5syZA6nU+r/U8st12PpjBjILK3Hfhz+ZNSmxscthN45fbT6SASelFBU3jBh1xDwL6tpW3zEAo3o3PY/EwUGCr09mYUTPmj8yvF0UZk2ormUUNXvn6W/4P376SknrAm5EYQOjTf26u+LhW3s1mACRdUzs74/frpbiMquOkw1ocRI0duxYnD9/vt7xzZs3WySg2jpB0dHRmDVrFrZt2wZ/f3+UlJTg/PnzqK6uGQWZOXMm5syZg+joaIwePRpSqRRvvPGGRWJoK5nUAdG9vLD3fM1Ilb9b80mQXOrQ4MRoCWo2V62d26CU/V6jhVfWyRJ6eDjBRdn8oLBSLkWYX83qnwh/V5y7VgYhBN774ZIpybjxsnNumRaf//z7buquKhkKKqpw7lopAGD62wfNjvGV7+v/zrnxdZZ+ehKFFR0z+ZrablF8b2w4lIbcUttc0Uv2o8VJ0IEDBzBs2DAMGTIE8fHxpp9///vfFglIpVJh06ZNOHr0KI4fP26qFh0dHY2srCyoVL8nFE8//TROnTqFQ4cOYcOGDVCrzZ9o2N5iwnxwOqsEGx8chl7ezs22V8kbToKAmkmpxZpqyBx+n8RqNApYsJo+UbN+fna86Rp8uL8aF7JLkZZfgU+PZ+K2Qd3hrJDWmdd2MCUfey/8fsna1VGOAyn5ePqL0wAApczBtB/ZjQrKdcgt1aLaYMRPaYUwGgX+mXSp0biSs8txpUiDgooqeLkwCeoM5FIHLJ8UgXX7LuHUlWJrh0N2rMVzgrKzs/H222/XOSaEwLVr1ywWVGf28nfn4aKSwc9VhbNXS7F4fJ8mtyOopZRJoWtkF3lXlRxXizVwVsogQc0IkE5vbHCVDlF7uXESokpeMyJ55mop1s2OgoezAh7OChRWVMH5+qhSal4FQm/4A8BVJccvmUXwdlFCCAG1So5yrb7eFg7//SULBy/mY8n4MCR+fQafPDwCrtf3NbtarIGPWlmnLERWcSUCPZxQWKFDoEfnnjNjTxwVUvx1Wj/858QV/O/0NfTxVWPyAH/T/x+ijtDi/20rV67EnDlz6h13d+cvHwBYt+8SnpkSAWelDJlFlWYPz6vkDtBWNzwS5OooQ3aJFi5KmekSWM0kUNtf1UJdlwQ1ic6kATUTXD2cFCiqrEKg5w17CN6QOLk6ypGeXwkftQKaagO8nBWorDLUm9BcqtVjYIAbckq16O6uQkGFDgEeTtBUGbBo2y94YkI4RoZ6odpghINEguwSHcL8XJCWX4n48Lrb65Btc3CQ4J6hgQBqaqC9nXQREf5q3DaoO+sJUYdo8bdoQwkQANx+++1tDqazE0LgwdEhuHdEMACgsspg1lwLoPYv67ojQT9nFMFFKYNaJcfqb8/B20UBbbUBEgBavfkrc4jai8FoNI3KeDrLca1Eiy9OXMHR1AKE+avhcuPGq44ypOVXwN1JgTKtHr6uSlRWGTB7/Y84n11qek4JAEgkyCnTwddVhYLyKvTv7orsUi2CPZ2QWVgzofaRj45j60+XYRACId7OOHu1BK6OHEXorML81Fg+KQI+Lkq88L9zKGjBxHui1uJQggWV6/To7uZoSnzWzR5i9l8zNUvk6051vlqsQWy4D9QqGVLzK+DnqkK5Tl9zOaza2Cnqm1DX5e6kQNEN5Ro8nBTYn5yHTYfTcTAlH9MHdkOwlzN+ziiCp7MCrio5Cq5XOC/T6uHnqoKmWg9npQzXGtgapqBcBz+1CtdKtOjbzRW5pVr4uamQc30y7YAAt5otZYRAkKcTLhdWcvSgCxjV2xt/Ht8HGw+l4+Mj6SjTNl06hKgt+C1qQbllOvjcUBNoyi0t2yyyRFNd5wNfXFkFdycFAtwdAQC+rkpor88b4kgQWduAADecyvp9mbunswJ7z+ciMtAdAjWbKvf0dsbeczkI9nKCm6McpRo9JABKtdXwuz4S5O2sQEF5/Z3sjQLwd1PiXHYp+nZTI7NIAxelDHqjqHcJLdDTCYE3bCRKnZurSo4nJoYjJswXGw+l49XvL+BCdpm1w6IuyKwk6MCBA+0dR5eQV6YzqzBiYzYfScfmw+mm20WV1XB3ksPLRYmHx/SE3/V6QxIA2moj5wSRVQ3v6YmvHh9tuu2qqrkcFuavhuv1DUaDPJ2w+UgGooI94KyUYfqgbnBRypBTooWvWoXKKgO8Gqg5pFbKUKnTw0etRHJ2GSL8XXEyswjd3Ws+A7mlNdt0+LmpUHH9svOHc4d1XOepQwR5OeFP4/pg0bje+Cm9EP9MuoiqRlbRErWGWd+ijz32GE6fPt3esXR62SVa+Daxf1JzhECdlWR6ozDNt3h2Wj8MDHRHdE8vCACVVXo4yjn/gWyHg4MEnz46ErNHBOPRmJptEVRyKRbG9Tbtuv7yXYPg6ihDVrEGvuqakU1HhQya60vra+theasVyCvXwff65TAPJzl+zSxBT++aOkXZpVr4uSpxz9AeWJoQZoXeUkdSyqS4LzoY0wZ2w/Pbz5ouiRK1lVlJkBACn332GWbMmIG5c+di8+bNyMrKau/YOpVqgxF//eo39DSjJlBTJDfuynTTfmcxYT64Z1jNSoriymp4OLfvNgRELdXQ9htPTAyvc1utkiOrWAM/N1W9/fJKtdVwdZTDV10z9yfQ0wnpBRWQSCS4VqJBL5+az1dOqRZ+riooZVIuqbYjwV7OeGZKX2w4mIZj6YXWDoe6ALN+e7z00kuYPn06AKC0tBRJSUn4xz/+gStXriAwMBAJCQmIjY2Fi4tLuwZry/LLdVgxuS+kbaxgqG2kVlAdQuCN3cn4eN6INr0WkTXU1r3yVSvxY2ohlDIHVOmNeH//JYzv6wcvZwX6d3dFdzdHuDnK8dmjowAA+56MMy06yC3T4dY+TW/1QV2To0KKFZMjsPWny7iQXYbZI4I4IZ5azayRoNoECABcXV0xY8YMvP7661iyZAlKS0tx++2328TGpdaUX1ZVZ1J0ayyICa1TByi7kSHf7FItknPK4c39kKgTcnWsWQ3mq1bharEG7k41I5rv70/Fb1dL4emsgLuTAq/NHAwAuKWHGwCYEiB3JzkuF1S0+4asZLskEglmjwhGT29nfHQkw9rhUCfWopm1586dw9q1azF9+nR4enoiPj4ev/32G5566il899137RWjTavSG3H4Uj7yy3XwbmPJ/hWTI7B2Twp+yyrBwMSdGNPHp8F2seG+2LVkbJtHnYisQa2So7iyGgqZAwoqdHB3rCmyOKKXFw5fzG+2wGg3NxWulmj51z9hdG9vqFUy7D2fY+1QqJMyKwl64IEHEBgYiAEDBmDt2rUICAjAhx9+iNzcXBw7dgyrV69GXFxce8dqk7JLtHjqP6eQV6az2MhMXrkOVQYjujWy8eqUW7qhj5/t7JNG1BKuKpmpFERxZTU8nOQI9HDClAHdcOJyUbP7f3Vzc8S1Ek1HhEqdwJ1DeuBYehEKK+qXWSBqjllzgjQaDbRaLe68804sWLAAt956KxQKblQIAHnlWviqlcgr17X5chgAvDsnCpmFlQjzc0GQJ+ueUNfj6ijZgHtmAAAgAElEQVQ3LXPOK9PBy0WJ+WN7QQiBP/3rF3g4Nf27pbu7I/Q3FRYl+7YgJhT/TLqIZ6b0tXYo1MmYlQR9+umnEELgl19+we7du7F27VpIJBKMHj0a48ePR2RkZHvHabNqR4B0eqNFiheO6u2FNd9fwPxbe8GvDcvtiWyVXOqAw0+PAwBcKdIgwKOmGKhEIsHE/n7Nfo581Ep8eUN9IiI3RzlG9vLCd79lm/ayIzKH2WtLJRIJhgwZgiFDhgAAqqqqcPDgQaxevRr79+9HfHw8tm3b1m6B2qq8Mh28LTACVEutlCGzSIOxYQ3PByLqCmonNQ8Jdq+zv9662VFmPZ7V0ulmcRG+eHN3Cnr5OCOM0wXITGbNCVq/fr3p3wUFBfj000/x+OOPY968efj8888BAFKpff5SKqyohqeTol5Nn9aSSCQorKjZLoOoq/vk4Whrh0BdyML43vjXT5nIKuacMTKPWSNBr732GlJSUrB79278+uuvUCqVGDNmDB5//HEkJCRg0KBB7R2nzTIKATdHeYMbQLZWUWUVPJy4/JeIqCWkDhIsnxyOv//vHJYlhMONv0epGWaNBCUnJ2PPnj2YMGECdu7ciaKiInz//fd44okn7DoByinV4s09KfByUaCgQtf8A8xUWF4FL9YAIiJqMaVMimUJ4Xht1wVUG7jPGDXNrJGgyZMn45tvvmnvWDqdwooqzB4RBJVcatr7yBLUKhlcVdwKgIioNdyc5Lh/VAg+OJCKP8b2tnY4ZMMkQlhoMouNqi2oZqluaqoMKNFUw99NhR9TC6A3CvT0dkZBeZWpsm1b6fQGKGX2OceKiMhStv10GRH+akQGeVg7FGoHlvh+b1HFaAK+O3MNz20/AwAo1erhqpKju7ujxRIgAEyAiIgs4A/DAvHVyasWHamnroVJUAul5lWgt0/NRrGlmmruX0REZKMkEgkeiw3Fa7suWOxqAHUtTIJayCgEdHojsoo1KNFUw9WRc3eIiGyVn6sK0wd1xwcHUq0dCtkgJkEtoDcYoZBKkV9ehdEv7cWvV4rrFHojIiLbM7CHOwI9nPC/U9esHQrZGCZBLZBbpoOfqxJr7hkEJ4UU14q1kEn5FhIR2brJt3TDtRINTmYWWzsUsiFNfoPXVopOS0tDXFwc+vTpg+joaPz6668dEpytuVqsQTf3mn2Oti8ag5wyyxVIJCKi9jVvTE98e/oarpWwojTVaDIJ2rp1KwBg2bJl+Mtf/oKUlBSsX78ef/7znzskOFtztUSLAPeaTU19XVXIKWUSRETUWUgkEixNCMO6pEtcMUYAzLwcVlRUhPj4eADAgAED2jUgW3alqBLdr48EuShl0FazGikRUWeikkuxML43XtpxjokQNZ0EHT16FGFhYfj5559RWFgIADAajSgtLe2Q4GyNtsoAJ8XvE6GTnoi1XjBERNQqfq4qPB7fG6u/PYcSTbW1wyEranJpU2VlZb1jGo0G7733XrsFZItKGqkH1NPb2QrREBFRW/mqVVg2IQzrki4iopsa0wd250IXO9SibTPS0tJw+PBhZGdnAwD8/f0xcuRI9OrVq90CbCtLlNUemPg9TiVOxGu7krE0IcxSoRERkQ34NbMY3/52DQHujrg7KhCOClbt7wws8f1uVpGb4uJizJ07F19//TWcnJzg4VGzD0tRURE0Gg1mzJiBjRs3ws3NcltH2JJSrZ7VRomIuqhBge4YFOiO9PwK/OO787grqgcGBHTN7zOqy6yxvwULFiA4OBgpKSkoLy9HZmYmMjMzUV5ejuTkZAQHB+ORRx5p71itokpvhFolQ6lWD5mDxNrhEBFROwnxdsZfp/XD92eyceZqibXDoQ5gVhJ0+fJlvPnmmwgNDa13X2hoKF5//XVkZmZaPLiOtvTTkzAaRZ1Rn4IKHUJ9XHD+Win8XJVWjI6IiNqbg4MES8aH4ctfspBVzHpCXZ1ZSZBOp0N5eXmj95eXl0On01ksKGv54kQWEr85g51nc0zH8sp0CPdT49/HMhET5mvF6IiIqCM4OEiwbEI43vvhEvQGlkLpysyaE3Tfffdh0KBBmDt3LgYNGlRnTtCpU6ewadMm/OlPf2rXQDuCs0KK9IJKRAX/Xjsit1SHAQGu+OhIBvzdVFaMjoiIOopKLsV90cHYdDgdD99qu4t/qG3MSoIWL16MgIAAvP7663j++eeh1+trHiyTISoqCi+++CLuvvvudg20I3i6KJBVVFmnbkR6QQUm39INoT4uVoyMiIg6Wh8/Nfaez0VmYSUCPZ2sHQ61gxYtkQcAvV6PvLw8AIC3tzfk8vr1c2yJuUvotNUG/O2bM5g9Ihh7z+fiT+P64OXvzkMmdeCyeCIiO6WpMmDt3hQsnxRh7VDoJpZYIt/iylAymQzdunVDt27dbD4BaoniymoM7OGOAQFu0BsFdHoD3v3hkrXDIiIiK3JUSOGrViKzsH7xYOr8LFYec/v27ZZ6KqsoqNDBw0lhup1XpoOfq4r1gYiI7Nw9QwPx72OdfwU01WexJOi1116z1FNZRVFFNTyda5KgtXtScCG7DJFB7lwiSURk55yVMrg7yTka1AWZNTG6dgf5ppw8ebLNwdTSarVYsGABzp8/D71ej9WrV2PChAmNtq+urka/fv0we/ZsJCYmtuo1Cyur0NdfDQB4/74o/JReiBE9vXAxt/HSAEREZB/uHRGEt/dexFOcG9SlmJUESaVS+Pv7o3fv3o22SU9Pt1RMSExMhBACR48eRXJyMqKjo3Hu3Dn4+fk12P79999Hbm5um16zpLIK7tcvhwV5OeHDg2mYP7sXvJwVzTySiIi6OieFDN3cVLiUV87Vwl2IWUnQBx98gEcffRQff/xxo2327dtnkYCMRiPWr1+PL774AgAQFhaGyMhIbNmyBcuWLavXvry8HJ9++iluu+22Nr3ujTvFd3NzxG9ZJfByVphmnxMRkX27Z1ggXtuVjKcn97V2KGQhZs0JCgkJwYoVK5CVldVom40bN1okoNTUVBQUFCAi4vchx379+uH48eMNtl+zZg0WLVoEqbRtu/5WGQQUspq3w1Ulg05vZAJEREQmSpkUfXzVOJZeaO1QyELMnhgdFxeHgICARu8PCQmxRDzIyanZsuLGHend3d0bvNyVl5eH/fv346677rLIa9eSSCTY92SsRZ+TiIg6v/8bEoC953ORUVBh7VDIAsxKggoKCsx+wpa0bcrNozANLVV//vnnsXLlSrOf78af5iZQ9/BgdVAiIqpLIpFgaUIYtv50GT9nFFk7HLuSmJhY53vcEsxKglqyJUZbt8/w9a3ZpLS4uNh0rLi42HS8VmpqKtLS0sxauQbUJFE3/rR2FRkREdk3udQBKyZF4Ke0Ql4a60C1i6ZqfyzBrInR165dw8cff2zWiza127w5QkND4enpiQsXLpgSn7Nnz2LKlCl12v3www+4evUqYmNjAQDnz5+HSqXCvn37sHHjRvTs2bNNcRARETVGIpFgQUwvvLjjPHp5O8PLRWntkKgVzNo7LCQkxOyhp8DAQOzfv79NQa1YsQI5OTnYuHEjUlJSEB0djTNnzqCoqAgLFy7Ezp07602Enjt3LkJCQuqN8Ji7t8hru5K5RxgREbVIuU6PV7+/gFXT+3ExTQezxN5hZo0EWbIGkDkSExOxYMECREdHQ6/XY9u2bfD390d6ejrOnz+P6upqUxJUVVWFCRMmmEaCTp06ZVpeT0RE1J5clDLcERmArT9dxuwRwdYOh1qoxbvIdzYcCSIiova2/kAqxob5IMxPbe1Q7IZVdpHvinR6A5QyvhVERNQ6D47uiS1HM6CtNlg7FGoBfvMDKNfq4aI068ogERFRPVIHCR6LDcVLO86jQqe3djhkJl4OA/BbVgkyCysx+ZZuHRUWERF1QQXlOmw+nA6BuvXuqg1GzBwaiBBvZ+sF18VY4nKYxZKgjIwMBAfb3qQwc96kLUczMK6vL7q5OXZUWEREZEeqDUa8s+8Swv3VmNjf39rhdAlWmRM0c+bMeseEEBbfuqKjGI0CO367xgSIiIjajVzqgD+N6wOjUeDtvSkwGLv0RZhOo8VJ0I4dO/Dqq6+abl+9ehXjxo3DiRMnLBpYR7lSpEEfX87mJyKi9jf5lm6Y2N8ff/vmDH65XGSxysfUOi1OghYuXIigoCD83//9Hz788EMMHDgQvXr1MlVu7mwu5ZVj+qDu1g6DiIjsRB8/Nf4yrR+yS7R4fvs5/JpZ3PyDqF20eE5QRkYGfHx8MHv2bHz11Vd46qmn8NJLL7VXfG3W3DXD9QdScVdUD7g7KToyLCIiIggh8Pbeixgb5oNBge7WDqdTscqcoPHjxyMyMhJFRUW4cOEC9Ho9nnvuOdx///2tDsJahBAorKhiAkRERFYhkUiwML43vjyZxRpDVtDiJCgjIwNz585FUlIS+vTpg1dffRXDhw/Ht99+2x7xtau/fnUGP6VxB2AiIrIeiUSCB0aG4F8/XbZ2KHanxUnQrFmz8PTTT9epfzBp0iTMmzfPooF1BLVKhrfujbR2GEREZOdCvJ2RV65Dmbba2qHYlRYnQZs3b27w+NWrV9scTEeTSR24NJ6IiGzCnOhgfPIjR4M6Uov3ioiPj2/w+MmTJ9scTEe6XFAJHa+/EhGRjejm5ghttQFFFVXwcOZc1Y7Q4iTo2rVrWLFihel2YWEhvvnmGyxdutSigbW3N/Yk44GRIdYOg4iIyOT+kSF4b38qVkyOsHYodqHFSdCLL76I22+/vc6xhQsX4r777rNYUO1NCAGF1IHLEYmIyKZ4Oiswoqcntp+6imkDWcOuvbV4TtDNCRAAyOVyXLp0ySIBdYRfr5Sgjx+rRBMRke2Ji/DFxdxyZBZWWjuULq/FI0HPPfdcndtarRY///wzvLy8LBZUezuWVoiZwwOtHQYREVGDHosNxQvbz2FpQhjnB7WjFo8EvfPOO0hLSzP9FBQUICEhAf/617/aI752UabTw1Ult3YYREREDVLKpHhqUjjW7LqAoooqa4fTZbV4JOjRRx9FYmJiO4TSgbhhHRER2Ti1So7lkyLwxu4UDAp0R79uajgrZXC5/nNjvT5qnRbvHdaYJ598Eq+88oolnsqibt5bRAiB13clY+mEcGuGRUREZBYhBM5eK8WlvAqUa/Uo11Ujq0iDB0f3RIi3s7XDsxpL7B1mVhL00EMPNftE3333nU0WTLzxTTIYBU5nleBKUSVn3RMRUadVbTDipR3nMW9MT3R3t8+iv5ZIgsy6HJaUlIS5c+c22UalUrU6iI6y8VAa1u5JwYGnGi74SERE1BnIpQ54cmI4Vn97Dn+Z1g9y6e9TfCur9EjNq4Cfqwo+aqUVo7R9ZiVBCxYswPLly5ts0xmSoBJNNUq1erg5cVI0ERF1biq5FPPG9MTL353H8kkRyC+vwn9+zkSV3oj+AW749vQ1jA3zQXSvzrN6u6OZdTnso48+AgAkJCSgW7du7R6UJd04XPbq9xcwY3B31ggiIqIuIzWvHF+dvApXRzn+b0gA3J1qltQLIbDxUDp6+7pgbJhPo48XQiC7VAsnhQxujp1nkKDD5gT17dsXTz/9NMaPHw+9Xm86HhQU1OoX7ii1b1JaXjl2/JaNx2JDrRwRERFRx3n3h0u4JcANo3t717vvWHohdpzORi8fZ2iqDCjWVEHm4ICJ/f3Rr7urFaI1X4clQXFxcUhKSgIAPPjgg/j+++8xadIkbNiwodUv3FFq36S/fHkayyaEd6osl4iIyBI2HEyDj1qJ6YNqFgUJIbD1p5od6+8dHlRnub222oCvTmahqLIaj47tZbNL8TssCYqPj8fevXtNt29Mimxd7Zv0xq5k/Hl8HytHQ0REZB1JF3Jx+GI+nBQylOv0mDTAH8NCPBttfzKzGDt+u4YnJoTXmXhtKzpsdVhjL3yjlStX4u9//3urA2lv3d1tf+I2ERFRe4kL90VcuK/Z7QcHusPdUY6XdpzHM1P6QupgmyNCbWFWElRSUoIjR46Ysq2bbwM1dYJsOQkK8nSydghERESdSoi3M2YND8KanRfw5MTwOoMg+eU6nL1aCkeFFEODPWz2sllTzLoc5uDgAIlE0uSQk0QigcFgsGhwllB7UvLLtPByYb0EIiKiljp1pRhfnbyKu6J6QG8Q2Hk2G04KGYaGeKBUU439yXl4JCYUAR1YuLHD5gSNHDmyyQ1ShRCYNWsWjhw50upA2osl3iQiIiJ7V6HTY8/5XMgdJIiL8IVKLjXdp6024PVdyZjQ3x9RwR4dEk+HJUH79+/H2LFj29zGGpgEERERtb/aFWfXirUY0csTpRo9knPKYBQCEgB6o8D0Qd3Rt5tllt53WBLUmTEJIiIi6jjaagN+zSyGWiVHmJ8LZNdXlhmMAp/8mIEKnQEP39qzzSvOmASZgUkQERGR7biUV45Nh9LxyNheCGzDoiUmQWZgEkRERGRbdHoD1iVdQh8/F0y9pRskEgmuFFUiu0SLMH81XFXNFzZmEmQGJkFERES26WBKPo6mFsBBAvi7OSLAwxG/ZZVApzdi3pieTe7ywCTIDEyCiIiIOpeCch0+OJCG/t1dMW1gtwZrEDEJMgOTICIios7p8KV8HEzJx5g+3hgVWncDWCZBZmASRERE1PVY4vvd9nZEIyIiIuoATIKIiIjILjEJIiIiIrtkc0mQVqvF3LlzER0djaFDh2Lnzp0Ntvvmm28wYcIExMfHY8iQIVi7dm0HR0pERESdmczaAdwsMTERQggcPXoUycnJiI6Oxrlz5+Dn51en3ZIlS/DZZ58hMjIS2dnZ6N+/P4KDgzFjxgwrRU5ERESdiU2NBBmNRqxfvx7z5s0DAISFhSEyMhJbtmyp1/bxxx9HZGQkAMDf3x9xcXGNjhoRERER3cymkqDU1FQUFBQgIiLCdKxfv344fvx4vbZLliypc1ur1cLHx6fdYyQiIqKuwaaSoJycHACAm5ub6Zi7uztyc3ObfFxpaSmOHTuGBx98sNE2Eomkzk9iYqJFYiYiIqL2l5iYWOd73BJsbk4QgHqda64Q0vLly/HXv/4VwcHBjbZhsUQiIqLOKzExsc4AhiUSIZsaCfL19QUAFBcXm44VFxebjjfk/fffh1wux+OPP97u8REREVHXYVNJUGhoKDw9PXHhwgXTsbNnz2LYsGENtv/vf/+LpKQkvPnmmwCAlJSUDomTiIiIOj+bSoIcHBwwf/58bNiwAUBNUnPy5EnMnj0b586dw7hx42AwGAAA+/fvx9q1a/H222+joqIC5eXleOGFF6wZPhEREXUiNpUEAb/XCYqOjsasWbOwbds2+Pv7o6SkBOfPn0d1dTUAYNasWdi3bx+8vb2hVquhVquRkZFh5eiJiIios+Au8kRERNTpcBd5IiIiolZiEkRERER2iUkQERER2SUmQURERGSXmAQRERGRXWISRERERHaJSRARERHZJSZBREREZJeYBBEREZFdYhJEREREdolJEBEREdklJkFERERkl5gEERERkV1iEkRERER2iUkQERER2SUmQURERGSXmAQRERGRXWISRERERHaJSRARERHZJSZBREREZJeYBBEREZFdYhJEREREdolJEBEREdklJkFERERkl5gEERERkV1iEkRERER2iUkQERER2SUmQURERGSXmAQRERGRXWISRERERHaJSRARERHZJSZBREREZJeYBBEREZFdYhJEREREdolJEBEREdklJkFERERkl5gEERERkV1iEkRERER2iUkQERER2SUmQURERGSXmAQRERGRXbLJJEir1WLu3LmIjo7G0KFDsXPnzkbbrlmzBlFRUYiKisIrr7zSgVF2DomJidYOwWrYd/vEvtsn9p1aQyKEENYO4mYrVqzAtWvXsHnzZiQnJyM6Ohrnzp2Dn59fnXbfffcdFi9ejJMnTwIABg8ejDVr1mDq1KmmNhKJBABgg93sEBKJhH23Q+w7+25v2Hf767slvt9tbiTIaDRi/fr1mDdvHgAgLCwMkZGR2LJlS7227733Hu69916oVCqoVCrMnj0b7777bkeHTERERJ2QzSVBqampKCgoQEREhOlYv379cPz48Xptjx07ZlY7IiIiopvJrB3AzXJycgAAbm5upmPu7u44e/Zsg21vbpebm9vg89YOm9kj9t0+se/2iX23T/bc97awuZGgWjef0Mau+fHEExERUWvY3EiQr68vAKC4uLjBf9/ctri42HS7uLgYPj4+ddrY42QxIiIiap7NjQSFhobC09MTFy5cMB07e/Yshg0bVq/tsGHDzGpHREREdDObS4IcHBwwf/58bNiwAQCQkpKCkydPYvbs2Th37hzGjRsHg8EAAFiwYAG2bdsGrVaLgwcP4vnnn0fv3r2bfP5PPvkEUVFRGDp0KJYtW9YlRoqOHTuG3r17Y9OmTY222bdvHyIiIhAbG2v6+c9//tNxQVpYdXU13njjDcTGxiImJgYjR47Enj17Gm3flc57S/re1c47ALz11lsYN24c4uPjccstt+Ctt95qtG1XOu+A+X3viue9VkpKCuRyOfbt29dom65aP665vm/atAmDBw+uc94PHjzYsUFa2Ny5c+v0JzY2FuXl5Y22b/FnXtggjUYjHnjgATFixAgRFRUlvv/+eyGEEEeOHBHdu3cXGo3G1PbVV18VvXr1Eh4eHiIgIEBs3Lix0ec9ffq08PPzE3l5ecJgMIjY2Fjx9ttvt3d32tUXX3whZs2aJYYMGdJk35OSkpq8v7NJS0sTISEhori4WAghxM6dO4Wzs7O4cuVKvbZd7by3pO9d7bwLIcTAgQNFTk6OEEKI5ORkIZVKxY8//livXVc770KY3/eueN5r/eEPfxBOTk4iKSmpwft37NghwsPDhUajERqNRoSHh4vt27d3bJDtpLm+b9y4sdH7OqsHHnjA7Lat+czb3EgQAKhUKmzatAlHjx7F8ePHMWHCBABAdHQ0srKyoFKpTG2XLVuGH374AYWFhc2OAn344YeYMmUKvL294eDggIceeqjT1xUaNmwYtm7dCrVabe1QOpRarcZzzz1nWh2YkJAAlUqFw4cP12vb1c57S/reFX388cemOYJ9+vSBh4cH0tPT67XraucdML/vXdWxY8fg4uJSb+7njbpq/Thz+m7vWvOZt8kkqKV69OhhVruG6gqdOXMGGo2mvUJrd+b2HQC+/PJLxMXF4dZbb8ULL7wAvV7fjpG1Ly8vL9x3332m20IIVFVVNfgLoqud95b0Heha5x0ABg4caPr3559/DhcXF9MfSjfqaucdML/vQNc77wCwatUqrFq1qsk2XbV+nDl9B2qSwJiYGMTGxmLdunUdEFn7W7BgAW699VZMnjwZSUlJjbZrzWfe5laHtaeG6goJIZCfn4/AwEArRtb+3NzcMHLkSCxbtgxarRbTpk1DUVER1qxZY+3QLOKHH35AcHAwxo4dW+++rn7em+p7Vz3vp0+fxsyZM1FRUYFPP/0U7u7u9dp01fNuTt+74nnfsWMH+vfv3+wffi2pH9dZmNt3Pz8/TJgwAXPnzkVeXh5iY2NhMBiwaNGiDorU8vr27Yv4+HgMGzYMx44dQ1xcHA4ePIjBgwfXa9uqz3yrL9TZoJiYmCavg/fu3Vu89957ptsXL14UAERGRkYHRNe+muv7zb799lvh5OQkjEZj+wXVQTQajRg5cqQ4fvx4g/d35fPeXN9v1pXOuxBCnDhxQvj5+TU4L6Yrn3chmu77zTr7eTcYDCImJkYUFBQIIYQIDg5udO6LTCYzzSMVQojdu3cLBweHjgizXbSk7zdbt26diIiIaMfoOt4999wjFixY0OB9rfnMd4nLYeZqqK6QRCKxy2usQUFBqKysRF5enrVDaRMhBB555BEsWbIEUVFRDbbpqufdnL7frKuc91qRkZGYOnUq1q5dW+++rnreazXV95t19vO+detWTJw4EZ6ens22Nad+XGfSkr7fLCgoCBkZGe0QlfU01afWfObtKglqqK5Q//794ejoaMWoOsbatWuh1WpNt3NycqBQKODl5WXFqNpu2bJlGD58OO6++27odDpcvny5Xpuuet7N6XtXO+8FBQX4/PPP6xxzdnZGRUVFvbZd7by3pO9d7bwfOHAA27dvNy2Rzs7OxuLFizFjxox6bbta/biW9P3ll1+uczsnJwfdu3fvqFDbRUv61KrPvIVHqqzq5ktCeXl5YsyYMaZhxNOnTwt/f3+Rn58vDAaDiI+P7/RLZms11/cHHnhAbN68WQghhF6vF7fffruYO3euNUK1mJdeekksWrRIlJWVibKyMnHmzBmxatUquzjv5va9q533tLQ0MWjQIFFRUSGEECI3N1cEBgaKt956q8uf95b0vaud95vdeEno7NmzIj4+Xuj1eiFEzRL5vn37mpbIR0REdJkl8kI03ffY2Fixb98+IYQQFRUVYvjw4SIxMdFaoVqEn5+fqSxEamqqcHV1NfXREp/5LpEEHT9+XMTExAg3NzcRHh4u7rjjDiGEEJcvXxbe3t4iMzPT1HbLli1iyJAhYujQoWLp0qWd9hp5LXP7fvDgQTF58mQRGxsrhg4dKubPny9KSkqsGXqbXLhwQQCo97Nq1aouf95b0veudt41Go149tlnRXR0tBg7dqwYOHCgWLVqlTAYDF3+vLek713tvNf68ccfRUxMjFAqlWLQoEHib3/7W6P144YMGSKGDBkiXn75ZStGbDnm9P2bb74R8fHxpvO+fPlyodPprBx527zyyiti1KhRYuzYsSIqKkp89NFHpvss8ZmXCNHJS6gSERERtYJdzQkiIiIiqsUkiIiIiOwSkyAiIiKyS0yCiIiIyC4xCSIiIiK7xCSIiIiI7BKTICIiIrJLTIKIiIjILjEJIiIiIrvEJIioE5s2bRqkUikiIiLw/vvvm47fc889uP/++3WCJykAAAnVSURBVE23v/76a/Tr1w/9+vXDDz/80G7x/P3vf0dwcDAefvjhdnsNc/zhD3+Am5sbtmzZYvZjxowZg549eyI2NhZnz55t1etWVlaie/fuyM7ObtXjLeWll15CbGwslEpli94DInvDJIioE9u+fTsGDx6M6dOn45FHHgEAGI1G7N27F3v27DG1u+2223DnnXdi69atiImJabd4Vq5cifvuu6/dnr8hu3fvRu/evesc+9e//oVbbrmlxc81b9487Nu3D/369WtVLE5OTjh9+jT8/f1b9XhLWbFiBfbt2wcfHx+rxkFk65gEEXVyCQkJ2L17t+n2L7/8gtjYWFy9ehVnzpwxHT916hQGDRpkjRDtipeXl7VDICIzMQki6uQSEhLw66+/Ii8vD0DNyMijjz6Knj17YteuXQCAkpISuLu7QyKRAACOHj2KSZMmISEhATExMXjmmWcghEB5eTmGDBkCiUSCmJgYaDQanD59GgMHDkSvXr1w5MgRGI1GrF69GiNGjEBsbCxmzJiBrKysRuMrLi7G3LlzMWrUKIwZMwaLFy+GTqcDACxevBh+fn5YunQp7r//fowcORIjR47E5cuXTY/PysrCpEmTEBUVhcmTJ+P555+Ho6Mjxo0bhyNHjmDZsmXIyspCbGwsYmNjUVhYaHpsWloa7rnnHgwZMgQTJ05ESUmJ2e/ryZMnMWLECMhkMmzZsgWTJk1CWFgYPvnkE+zevRu33347wsLC8Oqrr5oeM23atHqX4Q4fPozY2FjExMRg5MiRWLJkCcrKyrBp0yaEh4dj3LhxePbZZzFu3DgoFAocPHgQOp0OTz75JEaNGoVbb70VM2fORE5ODgDUedzKlSsRGxuLvn371kmEichMlt/4nog6klarFY6OjmLbtm1CCCGmT58uNBqNeOSRR8SUKVOEEEL897//FRs2bDA9ZufOneLEiROm27NmzRJbt24VQghRUVEhXF1dxa5du0z3L1++XBw8eFAIIcQrr7wiIiIiRGlpqRBCiNdee02MGDHC1HblypVi3rx5pttTp04Vc+bMEUajURgMBnHHHXeIp556ynT/7NmzRf/+/UVZWZkQQogHHnhA/PnPfzbdP2bMGLFgwQIhhBB6vV5MmTJFhIaGmu7ftWtXndu1Ro8eLcaPHy+qqqqE0WgUsbGx4vXXX2/0fRw9erR4/vnn6xxLSUkREolEfPTRR0IIIQ4fPixcXV3FG2+8IYQQIiMjQygUCpGXl1fneT7++GPT/U5OTmLfvn1CCCE0Go0YOHCgOH36tBBCiA8++EA4OjqKn3/+WQghxFtvvSVOnTolFi1aJCZNmiSqq6uFEEKsWLFCREdHC6PRaHqck5OT+O2334QQQmzcuFEMGjSoXp8CAgJMsRBRfRwJIurklEolxowZg127dplGWFQqFcaPH4/9+/ejuroae/bswfjx402PGTBgAN5++22MGjUKsbGxOHToEI4ePQqgZl7L7NmzTROtq6urcezYMYwePRoA8O6772Lu3LlQq9UAgAceeAA//vgjTp8+XS+2K1eu4H//+x8WLlwIiUQCBwcHzJkzB+vXr6/TbuLEiXBxcQEADB06FCkpKQCA9PR0HDx40DTRWiqVtmjO0YwZMyCXyyGRSBAVFWV63pYQQuCOO+4AAERGRqK0tBTjxo0DAAQFBcHd3R1paWkNPvajjz5CSEiIaR6WSqXC5s2bERgYaGrTt29fDBkyBACwcOFC9OvXD++//z7mzZsHmUwGAPjjH/+Io0eP4sSJE6bHRUREoH///gDqvmdEZD6ZtQMgorZLSEjA22+/jYMHD5qSlXHjxqGyshJHjhxBZmZmnS/ee++9Fz4+Pti/fz9kMhmeffbZOiua5s+fj+joaOTl5eHAgQOYPn266b7Lly9j8+bN2LFjh+lYcHAwcnJy6k1Grr2stWTJEigUCgCAVquFWq1GeXm5KfHx8PAwPUalUpmSudrLbDdO8PX29jb7fXF3d6/zvMXFxWY/tpZUKjXFWZuU3Pi8crncFO/NLl++DD8/vzrHBg8eXOf2jX0HgJycHOh0Ovj6+pqO1T7H5cuXERUVVe9xN75nRGQ+JkFEXUBCQgKeeuoprFu3Ds888wwAwNPTE5GRkdiwYUOdBAgAjhw5gnfeecf0pV5VVVXn/sjISAwcOBAbN27EkSNH8OGHH5ruCwoKwsKFC/HHP/7RdKy4uBjOzs714goKCgIArF+/vs6Kq9zcXFNi0ZSAgABT+9rnys/Pb/ZxtiIoKAiHDh2qc+zSpUtwc3NrNJnz8/ODUqk0zQECYPp37XtARJbBy2FEXcCgQYNMIzuRkZGm4wkJCfj444+RkJBQp314eDj2798PANDpdKYJ1Dd65JFH8Oabb8LV1RWenp6m4wsWLMDHH3+MyspKAMC1a9cwfPhwaDSaes/Ro0cPTJ06tU4Noz179mDmzJlm9SskJARjxozBBx98AAAwGAzYtm1bnTaurq4oLy8HALz66qv48ssvzXrujnD//fcjPT3dVJupoqICd9xxh+m9a4hUKsXDDz+MDz74AHq9HgCwbt06DB8+3HTZjIgsxNqTkojIMmbNmiXuvvvuOsf27NkjZDKZKCkpqXP8+PHjYvDgwWLo0KHinnvuEbfffrvw9/cXf/vb30xtysrKhFqtFj/88EOdx+r1evHiiy+KESNGiLi4OBEXFycOHDgghBDihRdeEEFBQcLf3188+eSTQgghCgsLxUMPPSRGjhwp4uLixIwZM0RWVpYQomYSta+vrwgJCRHvvfee+OKLL0RYWJhwd3cXjz32mBBCiCtXroiJEyeKyMhIMXXqVPHyyy+L3r17m+Kprq4WCQkJIjo6WowZM0bk5+eLefPmCVdXVxERESG++uorsW7dOlNcN/bxRjdPjE5JSRHDhw8XAERMTIwoKCgQY8eOFQBEdHS0yMjIEFOnThUKhUIMHjxYHD58WEydOtX0up988okQQohDhw6JsWPHirFjx4pRo0aJzz//XAghxGeffWbqa0xMjMjNzTW9tkajEUuXLjX16e677xZXr16t97jHHntMpKam1olTp9OZnocTo4maJhFCCCvnYUREjSooKKhTe2fr1q145513cODAAYu+zpgxYzBp0iQ8++yzFn1ea+rRowdeeuklzJkzx9qhENkkXg4jIpu2fPlyJCUlAahZqfbRRx9h2rRpFn+dbt264fPPP2/Tthm2onbbDJVKZdbcKyJ7xZEgIrJp//nPf/CPf/wDLi4uqKioQExMDFavXg25XG7t0Iiok2MSRERERHaJl8OIiIjILjEJIiIiIrvEJIiIiIjsEpMgIiIisktMgoiIiMguMQkiIiIiu/T/f+X5kyL35bkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Image('model_spectrum.png') "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Extracting a spectrum at a grid point"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "It is also possible to extract a spectrum without interpolating at one of the grid points. Let's first check what parameter values are available for the DRIFT-PHOENIX models."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'teff': array([1000., 1100., 1200., 1300., 1400., 1500., 1600., 1700., 1800.,\n",
       "        1900., 2000., 2100., 2200., 2300., 2400., 2500., 2600., 2700.,\n",
       "        2800., 2900., 3000.], dtype=float32),\n",
       " 'logg': array([3. , 3.5, 4. , 4.5, 5. , 5.5], dtype=float32),\n",
       " 'feh': array([-0.6, -0.3, -0. ,  0.3], dtype=float32)}"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "read_model.get_points()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The `get_data` function is used for extracting an spectrum from the discrete parameter grid."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "model_param = {'teff':1500., 'logg':4., 'feh':0., 'radius':1., 'distance':20.}"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "read_model = species.ReadModel(model='drift-phoenix',\n",
    "                               wavel_range=(1., 5.))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "modelbox = read_model.get_data(model_param)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Plotting spectrum: model_spectrum.png... [DONE]\n"
     ]
    }
   ],
   "source": [
    "species.plot_spectrum(boxes=[modelbox, ],\n",
    "                      filters=filters,\n",
    "                      offset=(-0.08, -0.06),\n",
    "                      xlim=(1., 5.),\n",
    "                      ylim=(0., 1.35e-15),\n",
    "                      legend='upper right',\n",
    "                      output='model_spectrum.png')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkEAAAGiCAYAAADgPBIcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd1hT59sH8G9ABZUliqKi4qi7bpQ60Vq3VVu1tdaqWOtorQNtbbU/cbRaR9W2rrfWYt1brLgH7oETkamy9wo7QJLn/QOJhCQYkpOcE7g/18UlOTl5zn2OCefOM0WMMQZCCCGEkErGjO8ACCGEEEL4QEkQIYQQQiolSoIIIYQQUilREkQIIYSQSomSIEIIIYRUSpQEEUIIIaRSoiSIEEIIIZUSJUGEEEIIqZQMlgRFRkYaqmhCCCGEEL1V4aIQxhiSkpKQn5+v2DZhwgTcvn2bi+IJIYQQQjindxK0Z88ezJkzB1lZWSi5AodIJNK3aEIIIYQQg9G7Oeznn3+Gj48PcnNzIZfLFT99+/blIj5CCCGEEIMQ6buA6qBBg3DhwgWV7YWFhahatao+RRNCCCGEGIzeSdAff/yBli1bYvDgwUrbx4wZgxMnTugVHBeoWY4QQgipuPRJY/ROgpo2bYr4+HjUqFEDdnZ2ioASExORm5urT9GcoCSIEEIIqbj0SWP07hhtZ2cHLy8vpW2MMcyfP1/nMv38/DBhwgQsXboUU6ZM0bjfvn378Ntvv0EkEqFfv35Yv369xqRHz1zPZIlEIjr3SojOnc69sqFzr3znzkUlh95J0Lp169CvXz+V7du2bdOpvBMnTuDIkSOwtbUtc7+AgAB4eHggICAA9vb2eP/997F161Z8/fXXOh2XEEIIIZWL3qPDBg4cCIlEgiNHjmDDhg04evQo8vLy4OrqqlN5Li4u2L9/P6ytrcvc7++//8awYcNQp04dmJmZwd3dHdu3b9fpmIQQQgipfPSuCQoLC8PAgQORkpKCOnXqKP69ePEiWrZsWe7ynJyctNrPz88PH374oeJx27Zt8fz5c+Tl5aF69erlPi4hhBBCKhe9a4Lmz5+PZcuWISsrC5GRkcjKysL//vc/zJs3j4v4NEpMTFRqMrOzswNjDCkpKQY9rqlZtmwZ3yHwhs69cqJzr5zo3Iku9E6CsrOz4e7uDjOzoqLMzMwwbdo0o4wMU9cpSlPnMJFIpPTj6elp4OiEobKcpzp07pUTnXvlROde8Xl6eirdx7mgdxJUUFCAvLw8pW3Z2dmQSCT6Fl2munXrQiwWKx6LxWKIRCI4ODio3Z8xpvRTWd40hBBCSEXg6empdB/ngt5J0EcffQRXV1esXbsWXl5e+PXXX/Hee+/h448/5iI+jVxcXBASEqJ4HBgYiHbt2lF/IMK70h/O7OxsrFu3Dq1bt0aNGjXw+PFjpefz8/MhEong7e1tzDAJqVQYY8jJyamUQ8mJZnp3jF64cCGqVauG7du3IzIyEo0bN8bMmTPx7bffchGfQkpKCsaMGQNvb2/Y29vjyy+/xAcffIDU1FTUqlULXl5emDlzJqfHJJVPamoqAKB27drlfu3NmzfRp08fjc+PHz8ederUQZcuXQAULT7csmVL9OjRAz179oS3tzdGjRqlW+CEcEAmk0EqlSIoKAgRERGIiYlBZmYmOnTogBEjRhgtjsLCQuTl5WH69OlwcnKCVCrFJ598AldXV0XXi/K4ceOGYj3LGTNm0Ehi8gYzkJSUFJ1e9+DBA9avXz9ma2vLWrVqxcaMGcMYYywqKorVqVOHRUdHK/bdu3cv69KlC+vWrRtbsGABk8vlKuUBYAY8Ta29fPmSAWCHDx/mOxSiwbp16xTvlzVr1mjc7+TJk8zS0pLVqlWLnT17lsXFxSled+zYMSaRSFhSUhKbOHEiW7p0KQsJCVF6vVwuZytWrGBt2rRhANj06dOZXC5ngwcPNvQpkgpKLBazHTt2MKlUWu7XRkREMG9vb9arVy/F+xgAmzFjBuvRo4fStvbt27MJEyawnTt3lllmbm4uy8/PV8STmJjIfH192fz585mdnR07duwYu3nzJrt16xbbvHkzW7RoEatZsyZzc3NjAFidOnVY7dq12UcffcS+++47pTh69uzJjh8/zgoLC996btnZ2Wz27NnMxcWFPX78mEmlUubg4MCePHlS7utUUYnFYrZlyxY2a9YstmbNGvby5Uu+Q9IaF/d3g2UH/fv3N1TR5SKEJEgulys+vAC0+vAS44qMjGS1atVimZmZTCaTMSsrK/bzzz8r7SOXy9mCBQsU76maNWsqfp83b57eMQwbNkxtIk9IaSkpKczHx4eNGDFCKUlp3Lgxy8vL06qMqKgoNnjwYAaAtWrVig0dOpTt2bOHJScnq+wrl8tZfn4+27t3L1u4cCHr0KEDa9SoETt79iyTy+VMLpezmJgYlp+fz3755Rfm4OCgKLc4tj59+jBXV1cmFotZr169WN26ddnYsWPZjz/+yKZPn87GjBnD9u7dy8LDw1lBQYHamCMiItjKlStZ48aNmY2NDcvOzla734ULF9gnn3zCALChQ4eyqKgoxXPh4eHsk08+0eoaVTRyuZwdPnyYHTt2jF29epV5eHiwgQMHssOHD7MXL14wLy8v5ubmxhYvXqz2fVBaamoqr/cz3pKg0aNHs8zMTMYYYyKRiJmZmSn9FG8TAiEkQRMnTlTEMHbsWDZw4ECl52NjY/kIizHG2P79+5mvr6/Stl27drEtW7YwiUTCU1TGJZFI2IgRI9izZ88U2zIzM1nTpk3Z2rVrmVwuZ+np6WzIkCGsYcOGLCwszCBxTJw4kaWlpRmkbFJxxMfHMwBs8ODBbPTo0czf31+RiKxcufKtNc5BQUGsQ4cOzNHRkW3evFnpfa8tqVTKPDw8FH9fzc3Nmb29veJxVFQUe/ToEUtKSmJJSUmK+wVXJBIJO3bsGGvYsCEDwO7cucMkEgl7+PAh27x5M3N2dmZbt25le/bsUfv6sWPH6txaYcrOnz/Ppk6dyv7++2+2evVqdunSJZUvXnK5nJ04cYL16dOH/fPPPyrPBwcHsylTprAePXqw8ePHs/79+7P9+/cb8zQUeEuCDh48qLgwrq6uLCIiQuknPDycubq66hUYV4yZBMlkMrXb69evz5KSkhhjRTfX4j8SjBVVExsjvoKCAjZp0iS2adMmxbbdu3ezKVOmMHd3d7Zq1Somk8nY8uXL2dKlS9nOnTvZjBkz1JYllUpZRkaGwWM2pNTUVDZ79my2ceNGBoDt2LFDZR+ZTMaqVKmieA91796d5eTkGCwmDw8PFhgYaLDyScUwbdo0dunSJY3Ph4aGKhKTyMhIxfaEhAS2efNm1q5dO3b48GFOPsNyuZxlZGSwqKgoRdOXMWsGCgsL2bVr19jIkSOZi4sL6927N1u1atVbE5yNGzeys2fPGilK4Zg0aRKLj4/Xal+pVMp+/vln1r59e/bBBx+wYcOGsaFDh7IRI0aw69evK3KAvLy8t74nDUUQzWF37twp13ZjM1YStGfPHgZAUY0bFxfHfHx82MaNG1n79u2V9v3tt98YAJacnMzOnDnDABi8GWTlypXs33//ZYsXL2YbNmxg169fZ2PGjGFSqZTJ5XK2ZcsWNmDAAPb7778rYpkwYQILDQ1VKic+Pp65ubmx4cOHs7lz57Jr166xv/76i61evZrl5uYa9By4IpfL2dChQ5mHhwf74osv2K+//qrx+svlcpaWlsaysrIMHtf69evZlStXDH4cYrokEgkbMmTIW/eTyWRs4cKFrEaNGkrNZdu2bdOp31BFc+vWLbZ8+XK+wzA6bd47pWVlZbG4uDiWkpKi8b0jFovZgAEDNFYEGIogkqCEhAR27do1JpPJmEQiYStXrmTLli3T2FZrbMZKggCwuXPnKtrIS/7hKa4FKun27dtsxIgRbPny5axx48acVxczVvTmdXd3Z7t372YjR45UVJn/8ssvbMaMGSw9Pb3M1wcFBbHPP/9c8VgqlbKhQ4cqaiuuXLnCli1bxvbv38927NjBPv74Y6N/CHTh6+vLvv/+e77DUOHl5cUOHTrEdxhEwHx9fdmqVau03l8ul7P79++z0NBQtX+HKqvc3Fw2evRovsMwqszMTDZu3DiDlf/zzz8bvXaNi/s7J0Pk7e3t0bNnT6xatQoHDx5Ew4YN8fLlS+zZs0ff4gVLLpcrhmpKJBI4ODhg/vz5WLFiBXJyctC5c2c8evRI4+vfe+89hIWF4fTp0/jss8+Qn5//1kVjy+u///5D8+bNERMTgx07dihm2Pzhhx+0en3r1q1hYWGBJ0+eoFOnTvjzzz8xYsQItGnTBgDQv39/9O/fX7F/cnIyzp07h2HDhnF6Hlw7cuQIvvnmG77DUGFnZ4ekpCS+wyAC9ujRI7i4uGi9v0gkKtf+lUX16tUhkUjAGONs5mGhCw4ORtu2bQ1Wvru7Ozw8PDBkyBCDHcMQ9J4sMS4uDps3b0aVKlXg5eUFHx8f+Pr6Ijg4mIv4BKtXr164cOECzp07h4kTJ2LhwoVo0qQJDh48iLFjx+LcuXNvLeO///7D559/DltbW4PMsP3kyRN8+OGH+PHHH1G/fn2dylixYgUWLFiAgwcP4sqVK5g1a5bGfWfOnIktW7YoJiOTyWR4/Pix4CYnCwsLQ+vWrfkOQ4WdnZ3SLOiElBYQEID27dvzHUaF4OTkhNjYWL7DMJqgoCDFF1hDcHR0RGpqKmQymcGOYQh6J0FSqRRA0WRUjRo1UqwcX6tWLX2LFjR/f388ePAAQ4cOhVQqxddffw0A2L59OwYNGoS6deu+tYx33nkHe/bsgaWlpcrSI1wICwtDixYt9CqjQYMG2LhxI2JiYnDgwIEyvzXVrl0bnTt3xrlz58AYw/z587Fw4UIcO3ZMrxi4FBsbC0dHR77DUIuSIPI28fHxOn+hIcratWuHwMBAvsMwmsjISDg7Oxv0GF27dsXDhw8Negyu6Z0E1apVC7NmzcLs2bMxbdo0AMDDhw8VyVFF5eLigiVLlmDNmjXw9vZGzZo1AQCTJk3C+fPny1WWpaWlQWqCCgoKYGlpqXc5HTt2xMKFC1GjRo237rto0SKsXbsWc+bMQcOGDXH8+HEcPHhQ7xi4cu/ePfTq1YvvMNSiJIiURS6Xc7pwZGXXrl07PH/+nO8wjCYmJgYNGzY06DGGDBmiVSuIkOidBO3atQt2dnaYNGkSpk2bhoSEBPz555/47rvvuIhPsKysrDB8+HC4u7vrXZahkiA+2NraYu/evRg1ahS+++472NraIjc3l++wFJ4/f453332X7zDUoiSIlCUmJgaNGzfmO4wKo23btggKCuI7DKNJSEgweC14jx494OfnZ9BjcE3vjtH29vZYvXq14rGjoyP++ecffYs1CadPn+akHEMlQXx9Y2zYsKHSN47i5j4hLG4bEBCAuXPn8h2GWtbW1sjMzOQ7DCJQsbGxcHJy4juMCqNBgwaIi4vjOwyjkUqlqFJF71t+mapVqwa5XK40cEjo9I7y7NmzcHd3h0QiQVBQEJo3bw5nZ2fcuXOHi/gEi8sEwxBJkEwmE8ybsFmzZggPD+c7DABAVlYWbGxs+A5DLTMzM8F1IifCQf2BuEXNioYhpL/32tD7Lrlx40bMmTMHlpaW+Omnn/DBBx9g2bJl8PDw4CK+SqF69eqcd4wWi8Wws7PjtExdNW/eHC9fvuQ7DBQUFKBatWp8h0GITigJ4p6ZmZnJjWbShUwmg7m5uVGOZWp9rfROghhj6Ny5M7KysuDr64uNGzdi6tSpqFq1KhfxVQqGqAlKT08XzAi9Fi1a4MWLF3yHgdDQUMXoRaGib6dEE0qCuFe/fn3Ex8fzHYbBpaenw97e3ijHat++PQICAoxyLC7onQRlZGSgoKAA+/btw+DBgxX9PoyVdVYElAQZR1BQkEEnCyPEkCgJ4l7jxo0RHR3NdxgGl5KSgtq1axvlWKZWE6R3L6lx48ahXr16KCgowI0bN5CZmYklS5YYfCgen7jub2NpaYm0tDTOygOAtLQ0wSRBjRo1EsQfmtjYWJpojpis5ORk1KlTh+8wKpTGjRsjKioK7733Ht+hGFRqaqrR3ju1atVCenq6UY7FBb2ToEWLFuHDDz9EzZo14eTkhLy8PIwdO9agM1PyLT8/HxYWFpyVZ2FhgYKCAs7KA4RVE1SlShXI5XK+w0BcXBwGDRrEdxiE6EQul1MNO8caN26M+/fv8x2GwRmzJggoatY3lSVJOKnOaNWqlWLoZvXq1dGvXz+sW7eOi6IFSSKRcDIJYbFq1aohPz+fs/IAYSVBQhEXFyf45gSRSCSIhJGQyqBRo0aIioriOwyDM2ZNEADUrVsXycnJRjuePnSqCTp+/DjGjBkDkUikcbLAc+fOVdhEiOuaoGrVqhmkJqhjx46clmnqhDRiTpMaNWogNzcXVlZWfIdCBEQqlVItkAE0bNgQMTExfIdhcCkpKWjevLnRjlfcD1Sb5aP4plNN0K5du5CdnQ0AuHr1Kpo0aaLyw2VNidBQc5hu+J4DxxSqZ2vWrImcnBy+wyACk5ycbBI3FFNjYWGBwsJCvsMwOGPXBAllWhRt6FQTVHKm5MmTJ2PZsmUq+/B9wzMkU6kJElISVHxz57OGQ+gJEFB0nbKzs1GvXj2+QyECQiPDiD6M3SeoefPmnK2oYGh69wny9PRUu70if2vhetK9ypAE8T1iwFSScisrK6oJIiooCTKcGjVqVPjPXGpqqlGToBYtWlTsmqCS5HI5vL29ERoaqtS518vLC7Nnz9a3eEEyRBLEdcdormur9GVnZ4eMjAw0atSIl+MLebmMkigJIupQEmQ4xVN4tG7dmu9QDEYqlRp1AmO+v/SWh95J0IQJExAYGIiOHTsqXWRdJ/+TSCSYOXMmgoODIZVK8csvv6gd1uzl5YVNmzYpdXRdtWoVevfurdNxy6OwsFDwNUFCY2Njw+vioKmpqUabMVUfVlZWyMrK4jsMIjAJCQk0x5WBVIYkiGimdxIUEhKCp0+fqkweqKmZ7G08PT3BGMPdu3cRGhoKV1dXBAUFqe0jsWnTJri5uel0HH1wXRNkiI7RQmNra8trEpSWlmbU6mBdWVlZKQYdEFIsMTGxQncx4JNQJnOtaGxsbJCRkQFbW1u+QymT3n2C2rVrp3Z7586dy12WXC7Hzp07MW3aNABAy5Yt0blzZ+zdu1evGLlmCn2ChKb4A8GXtLQ0k6gJsra2ppogoiIxMZE6yxtIRU+C+OoPaSr9gvSuCfr111/x5ZdfolOnTkpNU2vWrMGoUaPKVdarV6+QmpqqVC3Ztm1bPHjwQO3+O3bswLJlyyASiTB+/Hij9UESehIkxDlFbGxskJSUxNvxTak5jM/rRIQpNzcXNWvW5DuMCqmiJ0GZmZm81MYUD5Pv0qWL0Y9dHnonQevWrcPp06cRFhaGKlXeFJeQkFDushITEwFA6T/Mzs4OgYGBKvvWq1cPgwYNwpQpU5CcnAw3NzfIZDLMmTNHh7MoH66ToKpVq3LaMVosFgtqZBhQlATxuYiqWCxWzGouZNQcRohx1atXT3HvqYj4+gLYvHlz3Lp1y+jHLS+9m8MuXryIqKgo3LhxA1evXlX8jB8/XucyS8/noq46b+jQoZg6dSpEIhHq1q2LOXPmYOvWrWWWWfJH1z5LAPdJkLm5OadLJQhteDzAf8doU5gtGqDmMEKMzczMzGSm0NAFX/0hi2eN5pKnp6fSfZwLeidBrVu3Vjs79NSpU8tdVnHHP7FYrNgmFou16hDYuHFjREZGanyeMab0I6QkCOB2Ij8hJkG2tra89gkylSSIaoJIaQUFBUYd3kwqFr5qghwdHXVqESpL8cCp4h8u6J0E9evXD+7u7jh16hSuX7+u+Jk1a1a5y2revDns7e0REhKi2BYYGAgXFxeVfdeuXav0ODExEQ0aNCj/CejAEEkQl9LS0gSXBPFdEyTExFAda2trSoKIkqSkJBoZZmAVeWoKvmqCileSFzq9+wQtXrwYjo6OuHr1qtJ2XdpYzczMMH36dOzatQt9+vRBWFgYnjx5gn379iEoKAjffPMNLly4AHNzc5w9exY9evRAv379kJubix07dmDSpEn6no5WhJ4ECfGGz3cSZEo1QRX1jzHRDQ2PNzwnJyfExMSgTZs2fIfCudTUVKMunlqShYWF4CbuLU3vJKhXr164dOmSyvbRo0frVJ6npydmzpwJV1dXSKVSHDhwAI6OjoiIiEBwcDAKCwthbm4ODw8PrFixAnK5HNnZ2Xj//ffxww8/6Hs6WjGFJKhJkyZ8h6HEwsJC5wk0uZCTk4MaNWrwdnxtFa8iT0ixpKQkGh5vYBU5CUpLS0P37t15ObazszMiIiLQqlUrXo6vDb2TIHUJEACcPHlSp/IsLS3h5eWlst3V1RWxsbGKxyNGjMCIESN0Ooa+hJ4EpaWlCW5YohAWLxVCDG9jCjES46KaIMNr1KgRYmJi+A7DIPicI61Zs2Z49eqVoJMgvfsE3b9/HytWrEBBQQFiY2MxcOBAuLm5KfXrqWhMIQkyhTlxiHqm0I5OjCc6Opq3Nfcqi+KaoIqIzznSiucKEjK9k6Dly5fj3XffRdWqVbFkyRLUrFkT/fr1q7CLpwLCH61BSRAhFUdERASaNm3KdxgVWkVOgvjsD1lcEyRkejeHFRQUYMyYMcjPz8fp06fx4sUL2NnZGWUhU74IvSYoPT3dJDoBGwvVrBBTRktmGJ6joyPi4+P5DsMgGGMqa3saS3GfICHT+8oUD+c9efIkevfurbj5Vq9eXd+iBUvoSZBcLleavVsozMzMIJPJjH7cvLw8k+gUTYg6jDHqK2Zg5ubmvPxtquiKR4cJmd53yt69e6Nt27aIjY3FiRMnkJeXh61bt8LKyoqL+ASpsLBQ0EmQUNnY2CArK8votVSmMjy+mEgkglwu5+3bGxEOmUxG7wMjqVKlCgoLCwXd1cFUCTmR52TtsIEDB6JWrVro3r07MjMzUadOHZXJDCsSodcECbX5p3iuIEqCylazZk3k5uZW6C8SRDuxsbFo2LAh32FUCk2aNEFUVBRvc+oYglQq5T2JrlevnqCneeDk6gwePFgxD4GNjQ0mT56MK1eucFG0IBkiCeKqqUioCRDA34SJQpw8siy0dAYpRp2ijccQa13xLTk5mffpFYQ+QkzvmiCZTAZvb2+EhYUptf15eXlhxowZ+hYvSIZIgqpVq6aYCFIf2dnZsLGx4SgqbvGVBJlaTRAtnUGKRUREwNnZme8wKoXiJGjw4MF8h8KZ+Ph41K9fn9cYikeI9ezZk9c4NNE7Cfrss8/g7++PTp06KS2kyufswIZmqCSooKBA7WK05SHk4fF8LaJqaqPlaOkMUiwiIqJC3ZSFrEWLFrh48SLfYXAqISFBEEnQuXPneI2hLHonQUFBQQgICFCpwVi5cqW+RQuWoZIgLnrRCzkJ4qsmKDk5GR06dDD6cXVFzWGkWHh4ONUEGYmzszPCw8P5DoNTQqgJEnpzmN59gt555x21TTiDBg3St2jBkkqlejdblWZhYYGCggK9y0lOTkadOnU4iIh7fCVBSUlJcHBwMPpxdWVtbU01QQQArSBvTMVdEiqS+Ph4ODo68hpD7dq1kZyczGsMZdG7Jqhv376YPHkyPvroI6XOp/PmzcOjR4/0LV6wuB7uZ2lpiby8PL3LSUhI4P1Nr4mNjQ0vHQ+Tk5NNKgmysbHhpdmQCJNQhxZXRGZmZpBKpYKcZ00XQmgOE4lEipYOIa4mr/f/9OLFi+Ho6Ijr168rbU9MTNS3aMEyxAgsrr79JyQkwNXVlYOIuMfXzT0lJUWwtWPq1KlTB2FhYXyHQXhmiBpnUrZWrVohJCQE7dq14zsUTgjlS/G7776LgIAAdO3ale9QVOjdHNarVy+Eh4er/FTk5jBD4KqpSChvenVsbW15aQ6TSqUmNQGag4ODoKuPiXHExsbCycmJ7zAqlU6dOuHJkyd8h8GZvLw8vQfbcKFLly54+PAh32GopXcSdOnSJbXbN2/erG/RlQqXNUFCTYL46hNkaigJIkBRp2iaI8i4unTpUmG6ccjlcsE0pQr5unLW8JmamoqcnBzF4wkTJuD27dtcFS8ohnhjcZUEZWZmwtramoOIuMdHEiTkySM1oSSIAMDt27cF2XxQkRU3h1UEQqpJbNKkiWAXUtU7CXr8+DHGjx+PV69eKd1whJKBco0xZpAbq42NDWerGAv12ltYWHDS+bs8UlNTTao/EFC0+LCxrxMRHj8/P8ybN4/vMCoVc3NzyOVyQa91pa2QkBC0atWK7zAAFN2TatasiYyMDNja2vIdjhK9m8MWL16MCxcuoG/fvpDL5ZBIJDhy5Ah++uknLuITHEOtG2Ztbc1JLYmQaz74+KNiqvOsCPn/kRieVCqFRCJBjRo1+A6l0mnTpg2ePn3Kdxh6CwwMROvWrfkOQ+GDDz4Q5GSUeidBjDE0bdpU8Ue7WrVq+Pjjj/H48WO9gxOinJwcgyxsyUVzmFgsFlyWXZqZmRnkcrnRjmeqyw5UrVqVk3mjiGm6d+8eevTowXcYldLnn3+O3bt38x2G3q5fvy6opSqGDx+OM2fO8B2GCr2TIKlUWlSQmZmi93doaCj8/f31LVqQUlJSULt2bc7LtbOzQ3p6ul5lmEJHSkdHR6NOn2CqSZCzszMiIyP5DoPw5Pz587RcBk86deqEgIAAk66NLSwsRG5urqC+FDdq1AixsbFG/RKsDb2TIBcXF2zYsAEzZsyAq6srGjVqhPbt22PcuHFcxCc4hlqVt2HDhoiNjdWrDFNIgho3boyoqCijHe/58+eCaRcvjxYtWiA0NJTvMAgPZDIZbt68ie7du/MdSqUkEonQo0cP3Lhxg+9QdHbnzh289957fIehwtXVFdeuXeM7DCV6d4z+9ddfFb87OTnhzp07aN26NUaMGKFv0YKUmJhokCTI0tJS70Vnnz17Jvhvj02aNMGrV6+MVtUvhLVzdNG9e3f4+Phg+PDhfIdCjOzs2bMYMmQITZTII3d3d3z77bfo06ePSXF5epkAACAASURBVHaQPn/+PEaPHs13GCqmT5+OhQsXon///nyHoqB3TVCzZs0we/ZsAEUTJy5cuFDvBEgikWDKlClwdXVFt27dcOHCBY37btiwAV27dkXXrl2xbt06vY6rjZcvX6JZs2YGKVvf6ld/f3+VhUI9PT31KpNr3bp1g5+fn1GONX/+fDRq1Mgox+Ja586d9Zq0TWj/78Zkyuf+4MEDjBw5Eu7u7jq93pTPXV9cnnuzZs3QqVMn+Pj4cFamIZU+90ePHqFLly78BFMGJycnmJmZGbU14K2Ynjp27KhvESq+//579sUXXzDGGAsJCWG1atViCQkJKvudPXuWtWrViuXl5bG8vDzWqlUrdvr0aaV9ADAOTlNh2rRpLC4ujrPySpo6dSqLiYnR6bVSqZQNHjxYZTuX584FmUymNk5DAMAOHTpklGMZwvDhw1l2drZOrxXa/7sxmeK5y2Qy9s8//zAAbNu2bTqXY4rnzhWuzz0rK4u1bduW3bx5k9NyDaHkub969YpNmjSJx2jK5uvryz799FMml8v1LouL+7veNUGtWrVSu/LuggULdCpPLpdj586dmDZtGgCgZcuW6Ny5M/bu3auy744dO/DZZ5/B0tISlpaWmDhxIrZv367TcbUhk8kQHh5usBmZ+/bti8uXL+v02suXLwtqJIAmZmZmcHR0NPhCqmKxGAAEWSWsrcmTJ2PTpk18h0EMKD4+Hrt378bIkSMxdepU+Pn5YebMmXyHRQBYWVnh+vXr2L59O7766iukpqbyHdJbPXv2DG5ubvjhhx/4DkWjfv36oW3btpg/fz7i4uL4Dgei19mUzlatWoWLFy+qrCK/Zs0aBAYGlru8Fy9e4J133lHqezNnzhykpKTgwIEDSvs6OTnht99+w/jx4wEAx44dwzfffKM06WBxe66epwkA8PHxwZMnT7BkyRK9y1InOzsbw4cPx+7du8s1oik6OhpffPEFjhw5ojIxoEgkEtwoh8DAQCxatAj79+/nfPSCRCKBj48P/vzzT/j6+gru3MuDMYapU6eiWbNm+P7778u1ArMQ/9+NRWjnzhhDamoqIiIiEBsbi0ePHuHJkyc4deoU+vTpg0aNGmHMmDEYO3as3scS2rkbkyHP/csvv0RAQACcnZ3RuXNn5OTkoH///ujWrRuvM/RnZGSgevXqsLCwwMyZM+Hv748DBw6gcePGvMWkDcYYLl++jJ07d6JatWqYNm0a3nnnHTg6OsLMTPu6GS7u73onQbVq1UKnTp1Utj99+hRpaWnlLu/WrVvo3bs3JBKJ4o/+Tz/9hNu3b6vUklStWhWnT59WdAa+fPkyBg0aBJlMptjHFDu1EUIIIUQ7+qQxOo0OW7FiBYCi9cGGDx+utqlK3yrd0smLppOkJIcQQgghutApCTpx4gQ2b96M2rVrq02AAOjcN6e4CUwsFqv9vfS+xX0/ivdzcHBQWy5VEVc+dO507pUNnTude2XCRSWITh2j7ezs0LdvX9jb22vcJywsTKeAmjdvDnt7e6WVfAMDA+Hi4qKyr4uLi1b7EUIIIYSUplMSpE32NWPGDF2KhpmZGaZPn45du3YBKEqmnjx5gokTJyIoKAjvv/++os/PzJkzceDAAUgkEkgkEuzfv59GVhBCCCFEKzo1h/n6+hp0NlNPT0/MnDkTrq6ukEqlOHDgABwdHREREYHg4GAUFhbC3NwcQ4YMwfPnz9GrVy8ARbN80gy7hBBCCNGGTklQx44dy5y/hDGG+fPn6xyUpaUlvLy8VLa7urqqrK/l4eEBDw8PnY9V0S1btozvEHhD51450blXTnTuRBc6DZEfMGAArly5UuY+u3fvxuTJk3UOjCtczhNECCGEEGHg4v6u94zRmgghASKEEEII0USnJCg+Ph4DBgzAvXv3uI6HEEIIIcQodEqCgoKCcOXKFfTo0YPreAghHNuxYwfc3NzQvn17ODk5wc3NDe+//z6kUqlBjnfnzh107NgRR48eVdru6uoKNzc3uLm5Yd26dQCA/fv347vvvsOsWbNw//59xb6atpdU1j63b99Gp06dMH78eKSkpGD//v1wcXHBqVOnNMYdERGBgQMH6nranLt37x6qVq2KlJQUlef0uW7FvLy84Obmxkmsvr6+6NKlC/7v//5P6dp7enrC09MTAwYMeGsZI0eOVHqtWCxGQkICRo8eje7du+PBgwcIDg7GwoUL9Ypbn/eWtteXy2sLFC0d1blzZ/Tu3Rt9+vRBenq6XuVp+owW0+f8uY7V4PRaftUEgONV5AkxVZs2bWJr1641+HEOHDjAJk+ezI4cOaK0fdmyZUqPxWIx69q1K2OMsczMTNalS5cyt2vz2pKKY4iPj2dff/01y8zMLDPu8PBw9v7772t3kgaWm5vLvvnmG9akSROWnJys9Jw+1620fv36cRLv1atX2eTJkxWPS///+/j4lPn6Fy9esLlz56p97T///MO+/vprxePw8HCd49bnvVXe68vVtWWMsT59+rDw8HDGGGPTpk1je/fu1as8TZ9RxvQ/f65jLQsX93edRocRQkyPv7+/YrHh8pJIJBg9erTK9uHDh2POnDlK2z799FOcO3dO7fHXrl2LrKwsTJ06FS9evECbNm0AANbW1pBKpYiPj8ezZ8/Ubq9fv76irHv37r11HwCIiorCmjVrsG7dOlSvXl2rc2WMYebMmXj+/Dnkcjk2bdqE7t27QyaTYdKkSYiNjcV7772Ho0eP4sWLF4rXHT58GBs3boSVlRUcHBzg7u6uqFkqz/UDihagXrRoEf777z+V5zSduzbXTS6XY9q0aahduzaaNWumdM4LFy6Evb09YmJiMHPmTBw5cgQ7d+7E/fv3MXHiREyaNAlfffWVVtew9LnMnDlT7TE6duwIADh58iRGjRpV7rKLaXt9tXnf6Hp9DXlto6Oj4ezsjNzcXAQEBKjsW973l6bPqD7nr22sQkNJECGVhL+/P1atWqXTay0tLTX+0dTWjz/+iG7duiE+Ph5DhgzB999/DysrK8Xz1tbWSElJQUpKitrtJf/YarMPACxfvhyrV6/WOgECAG9vbxQWFuLmzZt49eoVPv30U9y/fx/e3t6wsbHB/v37cevWLRw8eFDxmpcvX2Lbtm24du0aAKBBgwb4/fffFc+X5/pdu3YNzs7OGlcC13Tu2lwTb29vVK1aFevXr8erV69w+PBhAICPjw/EYjE2bNiAqKgoTJ06FZcvX0ZkZCT27duH8ePHl/tm5uXlhZs3b+LBgweYOXOmxmMAwI0bNzB37lyV1wJFKwG0bNmyzGNpe331eW+97bWGurZRUVEQi8Xo27ev4otE9+7ddTp/bejz/tImVqGhJIgQnhx5EI2Y9Dy9y3GqVR3jujUqcx+5XI6kpCSVJKFYSEgITp06he7du+P+/ftwdXVFnz59FM9LJBKMGDFC5XUjR45UunmVpVu3bgCA+vXro6CgAObm5sjOzlY8n5WVhTp16qBOnTpqt5ekzT4AsHbtWmzZsgUODg4YN26cVnGGhISgZ8+eAIBmzZop+jSEhYUpluXp0aOH0sz5+/btw/Tp01GtWjUAQMOGDZXiKc/1u3jxIqysrLBmzRpkZGTg999/x8SJE9GqVasyz12baxIaGoqmTZsCAJydnRXbAwMDERcXhzVr1oAxhtq1awMANmzYAGdnZ4SHh2tz6ZRMmTIFY8eORXh4OKysrDQeIyUlBfb29qhSpYrKa4GihOjBgwdlHkvb66vPe+ttrzXUtX369CkmT56MTZs2Yd++fbh165bO568Nfd5fmmL18/ODh4cHvL29UatWrXLFY2iUBBHCk7clLlwKCwtDixYtFI+PHj2K/Px8hISEYNGiRbh69SqkUikCAgIglUrRpUsXpddbWlri0qVLOh8/ODgYd+/exZQpU5Cfn4/c3FwMGTJE0UE6KysL5ubmqF+/PmrUqIEff/xRZTtQ9E2zcePG6NGjh8Z9SqpduzbOnTuHAQMGwMbGBoMHD35rrK1atcKpU6fw5Zdf4tWrV7CzswMAtGjRApcuXcK0adPg5+enNDdJXl6eIinas2cPGjVS/r8tz/UrWVu3fft2fPvtt4qbTVRUlMZzL+u6FWvZsiXOnj0LAIiMjFRsb9euHaKjo7F48WIwxnDs2DEARcndxo0b8e233+LQoUNaxV9acWKg6Rj//fef2ht4eWh7fct637ztvfW262uoa+vv74927doBAMaNG4elS5eCMaaUhOv7+eTi/MuK1cXFBQ0aNMC6deswaNAgXLlyBe7u7jh79iyio6PRtGlTWFlZYcKECXqdgy4oCSKkEvD390eHDh0Uj48ePYopU6bA0dERUVFRaN26tdL+NWvW1PlY//77L/z9/ZGfn48GDRqgZ8+esLGxwalTpxAbG4vIyEj8/vvvsLW1xYIFC+Dh4YHs7Gxs27YNADRuz8jIwODBgxEUFKRxn2J3795VxNC/f3+cOXMGgwcPxpQpU7BkyRI0bdpU4zfwDz/8ED4+PujduzdkMhn++OMPAMDo0aNx5MgR9OvXDy4uLrCwsFC85osvvsC4ceNw4MABWFlZcbKQ8/r165GRkYHffvsN33//PQAozr881630uZ08eRILFixA3bp1ERkZiatXr2L48OG4ceMGli9fjry8PLi5uWHbtm04c+YM9u7di19++QXz5s3Dpk2byrx2gPK1b926Ndq3bw8Aao8BAOfPn8fff/+t8toPPvgAEokE//33H2JjY/Ho0SOV5Ly89Hlvve36Gura+vv7Y968eQCAatWqoXXr1vDz89Orman0Z7Rdu3Z6n782sRbXZjZp0gQA0KZNGzRr1gyDBw/GxIkTFUnQ295jXNJpxmhTQjNGE6LqyJEjEIlEyMrKwujRo3Ho0CFkZmbCxsYGmZmZ+O677/gO0WBOnDiB+/fvY/Xq1eV+bWFhIapWrYpbt25h9erVOH36tNL2wsJCjBo1CuvWrVN8I65INF07X19feHl5qV3u6G3279+Pzz77TOv9IyIiMGXKFPj6+pb7WEKmz/vSFEyYMAELFizA8ePHkZOTg2HDhsHS0hI3b96Eo6MjatasiQkTJpTrOnBxf6ckiBBSqURHR6s0V2nr448/RkpKCvLz87Fjxw7F6KYlS5bgxo0bkEql+OqrrzBlyhQOIxYOTdfu6dOn2LBhA/r164dp06YZ7PghISE4cOAAsrKysGHDBoMdhw/6vC+F7OjRo8jIyEBBQQFmzZql9Ny2bdvAGMPs2bMV28pzHSgJ0gIlQYQQQkjFI+i1wwghhBBChIySIEIIIYRUSpQEEUIIIaRSoiSIEEIIIZUSJUGEEEIIqZQoCSKEEEJIpURJECGlLD7mD7lcechlvlSGpSefQSZnkBTKAAB5BTI+wiOEEMIRSoIIKeXwg2iUnnXi39uROPowBqeexuJ/3gEAgDb/O4ffLoTgUVQ6nsdlKPbdezcSe+9GKn4v/vdmWArCU3KMcg6EEELeTpCTJfr5+WHChAlYunSpxplXb968iTVr1iAvLw9isRhDhgzBypUrYWamnNfRZIlEk7HbbsNjUCu817y2YtvTaDHGbL2Ff6Z2R2KGBIuP++PV6uFwXuyj2KeZQ9G6Wq+SlROaNR+9i/PPE3A1JBkAMLWXM/65FYFxXZ1w5GEMAGDTJ50wunNDQ58aIYRUeFzc3wW3gOqJEydw5MgR2Nralrnf0qVLsWDBAnz44YfIzs5Gt27dULduXcydO9dIkRJT9yAyHQmZeQCA74/6AwAOPYgGAEzedV+x35ITz5ReVzr5Kbb4uPJ+/9yKAABFAgQAZmYi5BZIYSYSoVAmh7VlVf1OghBCiM4E1xzm4uKC/fv3w9rausz9Ro8ejZEjRwIArKysMGLECFy4cMEYIZIKZP6hp4hJz8WhB9GKBKi0ffeiODuemQjotPwiVp8JwrueRe/X449iMGvvQzyMTOPsOIQQQt5OcDVBTk5OWu03b948pccSiQQODg6GCIlUINt8X8LFuRa6Odsrtu153W/HGO69SkOBTI7dd4qOyRjDgsNPAQCN7WugaxP7sl5OCCGEQ4JLgnQhk8lw/vx57N27l+9QiEDMO/gYmz7trHicKSnE/04G4OSTOADA6E4NFM/tuPbKaHGVTria/nBG8Xt6bgGAolhtXjeTlfydEEIItwTXHKaL9evXY/To0ejRo4fGfUQikdKPp6en8QIkBrfN9yWANx3kTj6Jw5NoMW6/SAEAdFt1SZEAFT8vNIcfxIAxhg6eFxCSkAUA6OBJTbyEEAIAnp6eSvdxLph8EnTmzBk8ePAAa9asKXM/xpjSDyVBFcOV4ERcCU7Er+eC8TwuA1/teYjsfCkA4N6rVGy79hKrzwahQCrnOVLt/Pi6E/bgTdd5joQQQoTF09NT6T7OBZNOgu7evYvNmzdjz549MDc3R1hYGN8hEQO58zIV26+9VNomzi2Au9cDuHs9AAAcuB+F+Iw8RLyei0ckAm6EpRi1uUtfB+6/6Zy9715R01kHz/Mmk8QRQogpMZkkKCUlBX369EFaWtEImqCgIMydOxc7d+6EVCpFdnY2VqxYwXOUxFDixHkIjs9EgVQOSaEMmZJCdFpxUWmfvXejEBCbiV/OBAEAfjkTzEeonFlyomhSxkyJFJ/83x2eoyGEkIpHcJMlPnz4EB4eHnjy5AkcHR3Rtm1bHD9+HNHR0ejSpQseP34MJycn9OrVC7dv31Z6bZMmTRAREaG0jSZLNH0x6bnYcCEUJx7Hwrl2DUSk5vIdEi/cWjng9wmdqaM0IYSAm/u74JIgrlESZJqkMjn8YzPQpXEtrPgvELtuhfMdkiD8901vAMC7TmVPJkoIIRUdF/d3k2kOI8KVmp2v+L1AKkdGXqHWr03LKYBczpTKAIAsiRRT//FDUpZEMAlQRwEkHiP/vIkFh5/wHQYhhFQIlAQRncnlDDI5Q9dVl/AwMh1yOYP3k1h4nnqudRkjfr+BuIw8dF11SWm7VM6QnS9F958vcx22zvZNd4Vz7Rp8h6GyuGtJUhl1oCaEEG1REkR09n83XuH/rheNvPrpZAB23QrHNt+XGqsmp//7AP/eiVDaxgCcf56oeHz2WTweR6Vj4s67BopaO9YWqvOIWllU4WxuCn28SMoGUDQ9QFpOAY69XpssKUuCoZtv8BkaIYSYFEqCiE6iUnMRnpyDX88VjcAKjM+EpFCGvEKZYp8Rf9yA82IfrDtftM/FwET8z/tNLZFfRBoYA1aeDlRsm7XvERIyJMjMkxrpTNRrWKu61vtO6elsuEDKsOXqS7xKzobHkafIlBTieVwmMiXaN0USQkhlR0kQ0YnX7QiVBUfXXwhFfIYEcRkS5ORLERCbCaDoZl26mSYhQ4Jx2+8gIVOi2Oa82AdA0fw+CZkSyOT8d2Yf0aG+0mMRgEsL+iptG9Le0YgRFZm97yESMyWYve8RAOBptBg/+wQZPQ5CCDFllAQRnSSWSF5Kux+ehl03lTsz/3n1heL3fKkMrqs19/UpkPGf/BT787MuSo/NzUQozs2WDGsDAHCwtjB2WDjzLAEx6XlIyirqUL716ktFMxkhhBDtUBJEdBKXkVfm8xsuhio93nTpzWzerZaeK/O13x54rHtgHGEM2Dqxi8r2XVNc0LROTQDA9L7NAADNHayMGps6d16lAgASM/PfsichhJBilASRcrkUmIiWS87icZSY71AMrplDUbJjV+PN5ISN7GugqrkZWtWz5iusckvJzsejqHS+wyCEEMGhJIhoLSA2A1/++wAFFXgY9kedGwIAbKpXQRWzopFgniPbqex3bl4fo8ZVHp+WWmLjUWQ6tl59qWFvQgipvCgJIlob8cdNvkMwuHHdGgEADkx3RYu6RbU9o18nRiUJYai8JndfpSE+Iw8/+wS+fWdCCKnEKAkiWpGUGPpeGVQxN+2PxsH70bgUlAQAr2vuhNPZnBBChMK0/9ITo2n9U9mdmYmwbL4cBjljuBaajG/2v+loXjxqTy6A6QcIIYRvqtPiElJCanY+ghOy+A7DaBgY1LV0HfrKVe3+6rbbWFZBpoTfyR4BID5Dgv95BwAoaiIDgBWnA+FgbYHbL1MwuJ0j3FrV5TNEQgjhFdUEkTIFxWdh4s57fIdhVN2d7VW29WhWW+2+6rZ3alyL85h0USCVIzI1FwCQnS/FnZdFw+jnHHgMcW4hcvIrVxMnIYSURkkQ0ej2ixR4HKkcK5Z3amTHWVma1k7j24S/3qzHZmYmQr5UhlfJNMEiIaTyoiSIaOQfm1FpJt/rX6JZiIuBXyXLsLYUXqszYwx+EemY/u8Dlec8Dj9FSnbl+H8nhFRulAQRjdacDeY7BIPq6GRrlONUFeBIszPPEnDgfhREIhHyCmRKtVfHHsWg26pLisd5BTLqSE0IqZCE9xWVCML10OS37iOUDsC6slJTQ1PX2hKuGvr/lMXB2gI+c3pDKmf4/pg/ejS1V3RGFjIRgEl/38O7TrbIzZfhl4/eVdln4G/X8H6buqhnY4mv+7eATM4gQlGTGiGEmDLhfUUlgvDFrvtv3efIzJ5GiMRwGAO+er3+FwPDmW/7oEVdK8wb2FKnsiACGthVBwDMdmtR4jnh1qKEJWXjQWQ6DvlF49CDaHRafkFln1hxHs4/T8C68yEAgN8vh+HfOxHGDZQQQgyAkiCiQtubditHa8XSEqaqTf03a4C1bWCjR0kMIhRdCwHnPBrlFhSNFMvKV1+zV7I1TFIoQ3puIc4/TzBGaIQQYjCUBBEVfddd1XpfE7zfG1x9W0vUtCi7pdncxJJHpT5BoqJFWVf8R8tyEEJMG/UJIkq8n8QiOi1P6/0b2FmWa/+KqmTtz9qxHZTWFlOXKFYzN0OeXNjz9FwJToSVRVUAQGpOAQDgfngaRBAhTpyHWDH9vxNCTBvVBBElPv7xGp979cswAMDKUW9WVbeyqIqmdWoaPC5D46IJqzjvKb24KmPAuw2NMxKNS+5eDzB+h/KK9ON33MH2ay9xNUS54zxjDAGxGcYMjxBC9EZJEFFyITBR43PFo4EmveeMH4a2BlB08zOthp03GDNe/x0Ta/3S2vxDRZNpMgaM+OMmz9EQQkj5UBJEFNJfN3loY0a/5m8eVNAbfHmUlUvVt7XEeJdGRovFmIo7R2uaXHHB4cox4zghxDRplQSZmZnB3Nxcq5/mzZu/vUAiOJ/vvIfOKy/q9NoP2tTjOBrj4WJ2aKDsGrF6Npb4rHtjbg4kMFI5w48nniFTUqi0Pe11Qn38USwfYRFCiFa06hjdsWNHbNq06a37Mcbw448/6h0UMb6bL1LKfF7doqLFRnduiJj0PPg809yfSKiKm8NmuemfvJfuC6TtcyaNAfvvRWFa76ZKm7usvIiINcN5CooQQrSjVRJUu3Zt9OvXT6sCa9SooVdAxPg2XQrV+Fy1KmYokMoxslMDtc8XJxFbJnaBz2IfQ4RnFJZVzfV6fWWdKqBAJgegvkVUUijs0W+EEKJVc9ilS5fevpMO+xJh2HQpTO12p1rV8XmPJmW+loFx1qRk6spzGVgFS5uKa7rypTLkFhRNuLjyNM0jRAgRNuoYXYkxxpCUKdH4/LFZPdGxkS0a2VcvowwoZkqe0tOZ6xANioFhaPv6fIdRIdx63Zz624VQtP3feZ6jIYQQ7eiUBI0aNYrrOAgPHkSmo/svlzU+X8/GEqM6NUTfdxwU264v6q+0z4DWdWFXo2hCPSGvkaVJ9Wr6NYMVM8FT59TSkwEAgB3XXym2mWmoItSmmUxSKKPmNEKIwZXZJ+iXX35R2cYYQ1BQkMECIsYjk2u+c68psZr4vIEtYVm1KF9uXFu5z9cPw9oYJjgT9LZmQSuLKsjWsDZXRbTnbqTa7e2WnceDJQMxdPMN3P3xfcV2SaFM0TdrzNbbCIrPxNP/DYLt6ySbEEK4VmZN0KZNm+Dk5ISGDRsqfpycnGBhYWHQoPz8/NCiRQt4eXmVud++ffvQtWtXdOvWDR4eHiZZE8Gn/feiND5X8ko6WFvA2pJuRPra/GknvkPghaRQBrmcITtfCkmhDDI5w5XgJCRkShS1PUmZErT+6ZziNflUC0QIMYIya4I6duyIgQMHokED5ZFB169fN1hAJ06cwJEjR2BrW/YyAwEBAfDw8EBAQADs7e3x/vvvY+vWrfj6668NFltFc+ppnNrtR2a+h9aO1mqfI+p5f90LNqUSxe5N7XE/PE3xuLLm6APW+2LSe8749VwwatesBgCQv74YrX86h4DlgyF7/bhAKke1KtRVkRBiHGX+tbl48aJKAgQAf//9t8ECcnFxwf79+2FtXfZN+O+//8awYcNQp04dmJmZwd3dHdu3bzdYXBWN9PXQ5tJOft0LLs72OtX8VNJ7PADAuU5NxbIixX7/tDNP0QhLXIYEv54LBvBmIda/brzpOxSSkKX4XSJVrgGqaKPoCCHCIrivXE5OTlrt5+fnh9atWyset23bFs+fP0deHq1srY0WS86q3d60tukvhqqtylozIwShidmK3wPj3iy8yuTAxouhyC1QTobyCqh5jBDCvXIlQR4eHoaKo9wSExOVmszs7OzAGENKivqZj0UikdKPp6enkSIVnpyyOufqMecPTRekjOZP0l5xQno5OBGbL4chodTUDR2XX+AhKkKIkHh6eirdx7mg1YzRxe7cucPJQbmi7iJo6hxNnabfoBsKERIGYM3ZouayBYefKj/3+mNboKH5lhBSeXh6eipVYHCRCAmuOUxbdevWhVgsVjwWi8UQiURwcHAo41Vk+X/PIS1jaLw+KM0sG9UMqReTnqexk76s1JeXFf/RLNSEEO6YbBLk4uKCkJAQxePAwEC0a9cO1atrnt24sovPyMM/tyLUPvfPFBf4ew6CjWW5KgdJGUonPVQZqd7/lZhgsTT/GDH+7/pLxeNdt8IBGPGeDwAAIABJREFUAPvuRSIgNkPTywghRCsmkwSlpKSgT58+SEsrGnL85Zdf4syZM0hNTYVcLoeXlxdmzpzJc5TC9t7qKxqfc2vlABvLqhV3tXM1DJ2TOFhZYOvELuqPTQmRVgqkcuy+rTzpYnRaLq4GJyMyNZenqAghFUW5kiBj9Kt5+PAh3Nzc8OTJE6xZswYfffQRACAvLw/BwcHIzS36w9e+fXusX78egwYNQo8ePdCpUyfMnj3b4PGZqo+33S7zeS6SH7qxKxOJRJwty1FZPYvNQKy4aMTn1eAkAECftVcRGJeBr/c/UuwXGJfJS3yEENPGWcfo3Nxc1KhRQ+Pz2uratSt8fX1Vtjdq1AjJyclK2yZOnIiJEyfqfczK4GFkOt8hCIaZCDBQt6gyUY5YfiX/n6Z6+Sl+j8tQHj027Pcb8F3ohpoWVeBgbdgZ7QkhFQdnzWEjRozgqijCsSvBiUY5jqm1pDnaWPIdAnmL8tQu/nXjFc49T1A8Hr9DWKNZCSHCU+4k6PLly+jQoQOqV68Oc3NzmJubw8zMDNeuXTNEfIQD0Wk0gaQ6xljL692Gtvi6fwsANI+SLi4HlZ3AD//9Bv69EwFAuabtUVS60pIlhBCiTrmHAs2ZMwfLli1D165dUbVq0dIKjDFMmDCB8+AIN5adem6U45hKn6DiMI3RCbyOlQXqWFkoHZdoLywpu8znn8dlIk78pmnsp5MBmOTaBP/ejjBwZISQiqDcSZCTkxM++eQTle0HDhzgJCDCLTkfnV8IMaLt14qG0JdMws1KJLiMMeRL5bCsqtxJXVIoU9lGCKlcyt0c9tFHH+HSpUsq2+fPn89JQIQ7OflSbLn6gu8wCDGKA/ejAACZkkJFLV98Rh7uvEzFF7vuq+zf+qdziBXnITtfChl9WSCkUip3TVCPHj0wduxYiMVixdpdjDEkJhqn8y3RziG/KCw+/syoTVRVzE2j1wtfzXa0dItxdPC8gLFdixZiLp4by8ayCrLzpSiUylGrZjXFvnI5w+x9j+DeyxlureryEi8hhD/lToKmTp2KWbNmoXPnzqhSpejljDGqCRKQtJwCHPSL1vpmX8eq2tt30sLioa01zkhNiDGZlcrHMyVSHH0QjVcpORjb1QkdnOwUz10PTcbk95qolPH1vkfYomGyS0JIxVDuJMjOzg4LFy5U2b5t2zZOAiL6O/MsHo+jxGXuU83cDAUyOeYMaIG577/DyXEtqlD/irJUptm4+XYtNFllm2VVc0gKZfjwz1uIWDNc6bni9fQYY4r/J59n8dhi+FAJITwqd5+gvn37Ijg4WGX77t27OQmI6E+b/g3bJxV9w/UY1ApVzE1m9RSTRs1hxpOYma+yrYq5mSLZufcqFcCbz4pcznD3VSr+522ckZSEEGEod03QjRs3sHHjRrzzzjuws3tTpfzkyROqDRIIbYbED2hdzwiRECIcIrzpD3bySdGq9d8cKFp648jDGHzq0gjxpWaiJoRUbOVOghISEvDnn38qbWOMIT4+nrOgiO5CErL4DoFogeqEjM/M7E1tXPFIsoy8QgDAleAkpOcWQC5nSk1iK08H4qcRbRVlMMaQlJWPejTbOCEVQrmToCVLluDzzz9X2V6yVogYX75UhpdJObgfnsp3KEQDSnz4JYJIUQNUrGQLZXE/urScApwNKFp+4++b4RjVqYGiI7WkUI5+664ieOVQxevyCmSISstFK0drreKIz8iDuUiEupRIEcK7cncGuXv3Lnbt2oWnT59CJpMBAKRSKerVo+YVPsWJJRj2+w14/hf41n3Xj+uo9C+XPndtzHmZhHAhNFG7WtJP/u8ulp4MUDz+8M9buPB6TTKRmsV3QxKzsOjoU43lXSixntnjqHSsPx+KIw9jyhE5IcRQyp0EmZubY+fOnejZsyesrKzQvXt3uLu7Y/jw4W9/MTGY/fcitd63eA6V4n+5ZE4joDSiK8Ovrb4vVbape7u+KLFUR/HzX+15iAKpvOgBK+pIXSCVI18qQxUzkcbBCM9iMvDVnodIzynAuO238deNV3gYKcw1zfKlMr5DIMToyp0Ebd68Gbdv30ZmZibu37+PcePGwcfHB1OnTjVEfEQLj6PS8deNcL7DAEBNPqRiKdlctujoUxzyiwZERc1kLZeeRZcVF2EmUk6CPE89R0RKDoCivkYAECvOg19EeolyhfdJGbD+GsS5BcjJlyK3QKp2n+Qs1VF3hJgyncdGm5ub491338WiRYtw/PhxhIcL4yZcGQXEZfIdgoJcgH/c32bhoJZ8h0B4Ep2Wp/W+kkIZlp16jgKpHD+fCQIA5BTIYGZW9L5/Gi3G02gxvG5HwG29LwBg46VQlXIiUnORKSlKMp5GizH34GPF78WKt5X07QHVbVzKKZCCMWCb70vs1PClyuVn1SWTCDFlnEwQ06dPH/j6+nJRFDFxprgE0zcDuJksklRs55+rXxrIXCRCaGI2Rm25hVFbbmlVVm6BFM6LfTBqyy14v+6sXfK13iU6cDPGIJcznHoap1IOF26EJeNlcjbEuUUj5RhYmU232izK/NPJAKw7rzqfHCFCU+4kqFOnTvjqq6/w119/4fHjx5BKpbh69SqaNFGddp4Y3o5rL/FTiU6cfBNiNb9Q0JWpmMxKr9HxmqRQpvZ3ANh7N0rp8Y5rqv2Virftvx+lcSHknHz1zVZlOfE4Bj8cf6Z4vPFiKA77RSseM6a+r1SxXr9eeesxwpKysOWq6jkRIjTlToIWLVoEGxsbHDx4EAMHDoSFhQUGDRqEhg0b4o8//oCvry9SUlIMESspJTwlB6vPCuvblinkQDP6NuM7BFKBaBoMsLzESM25B58AKBqmr466z3HxtmyJFFkakp2Oyy+UK1YAmH/oKQ7cj8LL5KIO4GYikaIZu/PKi2AoWuIlO1+KhNeTR+65G4mg+KJmd5pQklQk5U6CJk6ciPXr1+Py5ctITU3Fq1evcPz4cXTr1g2XLl3CF198QcPljWT8jjt8h6DCFPoE/TCsDS/HpdFhFZOmZqriCRkBQCovGllW3FFak1l7Hyo9Zowp1SBeCU5Uek6qpmmq5D5leX/DNQBFtT4lB1YUf4QvBSbil9d9n66HJiMqLVerckuWkSUpxOx9D8vemRAe6d0nqEmTJhg1ahQ8PT3h7e2NqKgoJCWV/UEn3NBQC88rE8iBeEOXpmL67aJq5+fSitcyyyssexj62YAEpGS/GYHV9IczAN4k0O5eDxTPzdr7SOm1eQUylX1uvUh56yzyJUetAa/7BImAeYeeQM4YJIUylea8YrtuvkmeCqRySGVypecLpHLceUkTuBLhMsjKmbVr1zZEsaSEF0lZaheJ5Nuk96hvGCH6eBSpnJQUSuV4mZyjst+5EpMwAkDHFapNY6f94/FAh3mJ1p4LUby+9U/ncCNMfReHFacDUfg68dl4KRTL/wuETM4UfYpENG8YEThaPtyETPr7HoCiCdgG/nad52jU6+Bkh34tHfgOQ/ioWoho8NUe5eajg37RuBT0pomr9Ois4jl9CqRypecCYjOQkJFXZu3sMm81gyo07H/3lfoaHZmcIa9Ahm2+L7HnbiRSs/OpRpiYDEqCTEjxt7Fdt4Q9J1NVc3pbqUM3BqKLWLHyXEZ91l5Velyyo3LJ53bfjsCNsBQwANkaOlbvvqM607ymfn3/3IpQ/F6yPJEIKCjZDCYC7oW/qX2itz0RMrpbCUCmpFAxw6w2+/o8izdwRIQQoYoV5+H7o/6Kx6JSzwFFnaYLZXJI5QyBcRn4eOttLD35DAGxGW8tX5vZ5zPyCpUel5wao+QIuNTsfMX8QwC0Oj4hxlRmEvTo0SP07dsXAwYMwIULb9qbx44da/DAKpNbYSlY5fP2hU8BYPmpwDdrGBGTQt0jCFeOPnqzAKv8dcJTUkx6nqLP4IH70ZBIZTj6MAYj/rjJyfELpXJF4nPiUazSF7OS49kOlph/6MRj7o5PCFeqlPXkvHnzsGPHDlSrVg1z585FREQEvvrqK6SmUm9/LsSJ87DnbiR23QxHvpaJzbFHprP6dKt61gjRcuXuyoCawwhXSq5VNvC367CrUVXp+YeR6ZCVesNJCrn78jRjz0PM7t8cwP+zd97xUZN/HP/kruPa0l3aUgpd0JYySksnq6UFZIOoTIEqIrhBXIhKFTdO1J8LGYLg9uf4ORCkKAIyBESZChQsu1CgpS20fX5/XO96I7lL7pJLcve8efXFXfLkyTeX5Mk33+c7gIc+343RGW2N697f2DLF9l5z9Nhz3+/DmywFbCkUubGpBDEMg06d9DlVvvrqK9x88804c+YM9fgXiZ7P2s+8aoCroKGSIQrxBpjWJ0ExBWYpFCkwnXIC9OHtWXGhxu/llfxz/PChur7BmAASAD7fUWH8/DpLdmuqAFGUil2foPPn9eGaGo0GS5cuxYkTJ7Bx40bJBaPoefzrv/DehsP4s0I5RVLVhp+PTV2fQnFLtlmE2ouJxglv0iNna6zyCVEocmHz6bBkyRKrZa+//jpuuOEGyQTyFPjU2Npy+JxZRIba4CoR4KmoIZs2hcKHY+dq7TfioPCFMrw0Jh2jM2NFlIhCcQyb+nxiYiJCQ0OtlhcUFEgmkKdgyARrCyWWxVAlClE+Qvy87TeiyML0AlpPjkLxRAQZNWfPni2VHEbq6upQUlKCvLw8ZGVlmUWlmXLp0iVMmjQJ2dnZyMnJQUlJCaqrqyWXz1lOXKhF2X77ZUUe+HSXJPvvmeS6bN5K8QkykBodKOv+813421OEwVUElUKhuDeClKBNm6S3TJSWloIQgs2bN2PlypUYN24cTp2yLgg4f/58lJeXY/Pmzdi8eTPKy8vx5JNPSi6fs6zacgzP2qn8PuTVX/DxNmmiwLq2DZakX3OUpfwYuG9gitwiUBSKJ+tAeYlhLt9n7dVGLN+sjyIz1DyjUORAUR6jTU1NWLRoET7//HMAQHJyMjIyMrBixQorK9Rff/2FnJwcaLVaAEBOTg527NjhcpmFwgCcxQj/Pl2NDpGtsOeEdE7QclVQp1CUjCf7r+m8tS7f59wv9OU6fj5wBn/8W4XC5Eg8d303l8tBoSgqY/ShQ4dQWVmJ1NRU47K0tDRs27bNqu2QIUPw008/oaamBpcvX8ZPP/2E3NxcV4ormCe/2YNX1x7EEY5wVXsWIop4FKa4vr6ZaWoJpU0VUjwXOacCf9xzCg2NBB9tO2a/MYUiAYpSggzTXsHBLVM2ISEhOH3a2ofmjjvuQGFhIRITExEfH4/MzEw89thjLpPVESxrAFlyrqYe244Ir/isVBTzds0yyC+9KUcGQSgU5eGllfc+Ndyev/59Fn/8W2VczieClkJxFkFKkKsuSstkjGz7feqpp7Br1y6Ul5fj6NGjOHToEJYuXWqzT9O/0tJSkaW2z+lL9TbXV9c34Pq33CcijFo7KBTl4+VM0h9R0I/3Exf9hlkftSRgHPDyz8bPjU2EswgsxXMoLS01e46LgaIcoyMjIwEAVVUtbwNVVVXG5aYsXLgQ06dPh06ng06nw4wZM/Doo49y9k0IMfuTQwnabid52YFTyo9uUzOe7PxKsY0nXxtjs9shWMb0DWerW14O/zlTYywJ8vfplvFw3b7TmPmh3ufTMIZTPA9D4JSY14BorwCXLzuflj0pKQlhYWHYv3+/cdmePXuQnZ1t1fbKlSvw9m65cb29vXHpknLrVB0VOW09RTh03KRw4cE6EPx9tAjUKSdGJunhbzHn890AgJMX6gCYJxp9b8NhY00yCsVZRFOChg0b5nQfGo0G06ZNw+LFiwEABw8exM6dOzFx4kTs3bsXxcXFaGzUR1b1798fH330kVEjXLlyJfr16+e0DFJx14fKj1yTishAX7lFoFBsI6IpqG+y653unUGjYRT3grBqy1EAQN4za1FeWYNd/1bBoKpeqmvAxTo6NUYRB8FK0Nq1a9GtWzf4+flBq9VCq9VCo9Fg/fr1oghkMHfl5eVh/PjxWLVqFaKjo3HhwgXs27cPV6/qCwW+8cYbxna5ubloaGjAO++8I4oMUrBPwrB3pWIYWP18XB+Ca6AoNRI6b7l9HiiexHPXdZVbBEF4aRizqvT2iGjl2pea3RUX8Ma6f4x6Kpu+Ou/LP10qE8V9EGwDveuuuzBv3jz06NHDOB1FCMH48eNFEUin07E6OOfl5aGioqVScWRkJD766CNR9ukKfLw0qG/wjKKBSnqrvKV3AvISw/H89/vtN6Z4LGJOh7UJ9hOxN+mYnB+H9zeVQ6th0CjopnXtDf7r32cBANV1Dfhsuz6J7KotR6Hz1uD2wg4AgGWbyvH4yC4ulYviHgh+RY6NjcXYsWPRoUMHxMXFIS4uDvHx8Vi1apUU8rkNl2Q230a08pF1/3LBMAw0Gk/2+KDwwRMdowd1jgYAaBhh02Gufsn5tFnxOXDqEt5a/w8YMDhzqR7Pf78fZftPC7JiUSiWCFaCRo8ejTVr1lgtnzVrligCuSPHzsnvFP31Xb0R6Ota50clPFiUIAOF4gjx4f7S7qD53tAwDAzWnSFdo6XdpwNcbdTLRqCX9dj5lvG0ZMlWHD5bI5NkFHdAsBKUm5uL6dOnIzw8HImJiUhMTERCQgJ++OEHKeRzC/o8v07W/Y/LbodQfx9sf3SAS/d7XWYsYoJ1Lt0nheIIiknsaULnGFfU+QM0DNBEgMdHdOaVe+W56+Qpb9FECBimxTJk4OHmSDJTlm86guWbjrhELoq6EWwauOmmm3DbbbchIyMDXl76zQkh1BKkYAakRclSH6hzTDDiIwLsZsqmUNTO0K5t8L/dJ8Tt1EG9LDLQ125iVrPdmOyHzy77p0UJF0oEmpoIjrJY1bewZNl/9Mu/AACT8uOlFouicgQrQSEhIbjvvvuslr/55puiCERRP3OHdsLaffpSJ66ejrqldwIWqSCHiJKcxynOX6dSZEd3VCS+x2KwfjEMY0w8J1YWXi4MVidHIOD3m9BEihQhCJ4O69u3L/btsy70uWzZMlEEoqifxNatZNs39QGiuAtSKyQG/Ly16Jeqz8ov9R6fcCKC61JdA2quNHKur7vaiDMCLGAUCuCAEvTLL78gOzsbmZmZKCoqMv6pKVzdldRd5b5pXQVVDPTQ34HChTtdGslRgbzaGUplBPh44aUx3QFIf49I2f+Ytzch+6k11MpKEYTg6bCTJ0/i9ddfN1tGCMGJEyLPh7sBhBBMfm+L3GKoflAY2T0GX+487nQ/VAmicKGEa+OOfkl4Y90/xu+OitSaZzLDtJggqx2Z7rN9mD+rDw4XMwqS8Nb6f+w3lIg//r1g/Mww6h/3KK5BsBI0d+5c3HjjjVbLQ0JCRBHInZj31V+sTnuehisHIzrwUdRKbKh5SLyrFDMzx2iTL/4CMr0zDNA21H6SyNToIEGyUShSI3g6jE0BAoBRo0Y5LYy78dsh+wqQhgG8tQp4DaVQKA4jhfLt6KjAJkpJz3jOmmaG/VgqXTf1iseU/DjO/QxPj7HYsf0foUdcqN02zkJMRDGtUE+hsEGLKknAS6v3I/6h/2H/KdtV7Wf1T8a4nPaSWy+UYOqXCw8+dIoAXOWE7ArYoqO0GgZcidNNj52xWG7Its5mFXptfAZvmXp1CDd+7tMxgvd2jvDC6pYSOY4mql2xuVwscSgKhypBInPh8lUs/OlvXm3v6d8RT1+rrmKLfCkdniZaX3SKi0Lhj9DbheH8It5+/UzylC2fmuv4TnjwZlmLX9LFugas239acB+P/JcWZPUUeClBv/zyi9RyuA3pT6yWWwRFUNIrQW4RKBTeOGIIah3Y4oAsyXSYg9YpU1nimktvMLCv3xBiKKFhu08uQvxt1yd8ZGgn4+ec+DD7HYrA0coavLiaFk+mcMNLCbrtttuwe7d1anKKOUpN0hUXHiDr/qVIJEehyE37sBZHZjGucbGGD9Nu2pk4W4cGsCsplhmjfb1aHgsGmewdHwNgRHoMRnaPwctj083WfXt3HwDmv5erZh/1CRb1O3v4C/oMo1jDSwkihOCTTz7ByJEjUVJSgmXLlqGiokJq2VTHewrMVJwaHYgkGZMXuho3cu3AhNz2covgMYhdO8yQg0cOTF/GDPcDwwBPjrJOVOjrpTHJHN28PYCXxqSjY2Qrs/spgEe0GCEtv+XH0/MBGELxzX/f0hGdWbfv3UF8f6G6q43Ye+IiVv52lFd7pb7MUqSBlxL07LPP4oknnsCXX36JhQsXIiQkBM899xxGjx6Ne+65B9988w2qq6ulllXxVNc38GonJPTUWeR2+BRjPHEnxUYIXI6sFHl4nOPBzUa3WH7FT00dhh2BzcmY65YzKCc6b41ZW/MQef3C0ZmxyGgfatzG8j4O0plnVzGMM6b9ZcWFct67ndqwh8ozjLklSgwOnq7GjBXbAQDr9tn3Dyp8oUzU/VOUDa+rbfjw4cbPQUFBGDlyJF5++WXMmjULFy9exKhRoxAe7tzN7EmsnV1g9l3K9w45n6NKUF6U+k7HRy76QqoslHA9W8L6ksNy3Zi2m2+jdAUDxmzqy/CZWPSx/dEBrNvf2icR+YktzwLTa9hSrJ5J7M+M568Xr0q9Yf8Xa68CAGqu2H9RLa90LKKMok4Eqdx79+7FwoULMXz4cISFhaGoqAh//vknHnjgAXz//fdSyagaXllz0G6b96ZkoU2w/aRiYuEj8luVHHiiMjDZRn4Wivg4XUCV2P4OAKH+8k2RaTUM/Hy0GNwlGt5a9jFBq2Hg7cXlGG1+QFx9dI0NRmSQDoD93/ThIZ1sNxCBpma56xuaAABzPjP3C3pjnXUkrxKVXYp08HpCTpkyBe3atUOXLl2wcOFCtG3bFu+99x5Onz6NrVu34umnn0a/fv2kllX1ZMeHorhTlNVyKe+596fmSNg7P+TMGE3HMwoflHCdiBVAwNbPrP7JAFheikyajstuj4cGsysmXhoNBqRZj11ikt48hdjQKN6A0dhcsv5yc+HVS/UNGPH6BuP6BT+0RI4Nf02/XAnXAsV18CqbUVtbi7q6OowePRozZsxAnz594ONjOxySYs2wbjH2G4mMzst1/kcUijtj+XBke1gKSSAodH982zQ1AT5aDa40NhmX2bIIGywfWg0DnZcGAztHW7Xx9dbg5bHd8cUO9oAYe7Lysa7c3q8Dlm8qNyouYvDk//ZaLTOtMWbK7gr9chF3T1EBvCxBH3/8MU6fPo05c+Zg+/btuOGGGzBq1CgsWLAAO3bskFpGt6HAJG394C7mA83P90tjSfM0066nHS9FHBy5bky3MTw3DaUkhFh1Xh3XnaN/+0Jx7WXXvIG89u3rpTGznnppNXhjQmaLDM3qjWH6S4iSJzQowyCHrzf3Y0ms6LEPfivHyDd+BQC887NpcsWrovRPUQ+8HUYYhkFmZiYeeOABfPnll/j444/Ro0cPPP3004iKisL48eOllFP1TM6PQ3xES76eN2/sYfxMAPhJFDFGdQIKxTW8PsG2gsAZtcWhLEQH86sGb70fwjmeWE4Xb5nbHzpv7rHHoMwZXtKs6oU144gSaZDFMsItjCOfEQCM7C6ONf14VS12HasCADz97T7j8p7P/CRK/xT1wEsJWrRokfFzZWUlPv74Y9xxxx2YOnUqPvvsMwCAVkunXWzxhI2IDIBaMGxBrdMUMUmOss6bJUaeINPpbiF+cPby0ggNGReyb74vX0Jf0oSkAbl3gN5fKbN9iF0L2ujMWEFycHHqInthVb5pTijuAy+foJdeegkHDx7EmjVrsGvXLvj6+qJ379644447MGDAAKSnp9vvxI15Zc0Bs8ysQpFS/5EzT5Ax26yTWgxNXkYRE7bIJl63iUUjpV6WcohlqUTueMw6hN4qgq5ZUm+tBmltghAZpDNLtvju5CxMe3+b2TZakZJnNVk4/tgaY+77ZBdKR3RGK19ej0uKyuB1Vg8cOIC1a9di4MCBWLBgAXr37k0do02wFxr/6Yx8u32Ipar06RiBXw6eFb1fR1GKhatfSmv7jSgegdKUF64XFUetU3IcX75Fzh9fi4AMR46EKxptdGZbfP67cxULGi1+pIQ533K2/XT7v+jfKQqDLPw4a6804t1fDuHu4o5OyUKRF15K0ODBg/H1119LLQtFBLLiwsyVIIUoIRSKUujePgR7Tlx0uh8zx2geeYLsIZ7y4notaFKe7bxW/VIjER/huLVcbL7ceVxQe7ZxtOZKA5ZuPEKVIJXDa7KZKkDsEEJwpaHJZpvVs/o2187hJrF1AGf1ZqFYdiN32QwxcSTZnDsdP0UcAnWOTWtYh8gr89qypUxZlslwpi8hjM9pj9Ro2+OgYX+W8lmer6LUSHx9Z28AME5Rjc+Rts7ek//bg8ss2aaVeQVQhKD+dMIyUrb/DJIf+c5mm+SoQPj72B50f5jZF15akZQgUXpRJn42olgMKG2qg+LerLm3r9N9iHnPFqdGqiaQgBB+97ThgAy/9bBuMejanFixJceR+PI1NRGj4nPsXC1ufX877l61A2erW5yq1fJbU7ihSpATXG3ktgLNH9kZHSL5VW+n1gp+ODPg0MGKYgtHLbEdIgNZl7NFOQlV0B3NXaSWQAICoCNLpB4XbL+14ScSy5JuysZ/KpH++Grj9w1/n8VXu44j68k1ou+LIh9UCXICW8rLpPx4rLm3gHO9kL6EYpqUUW7EGI/Fqn+mpOkLtTyoPAU+t5+tNpbr+J5eW9ekmi6RLm2DWZcPYsk+bQoDoHNMED6angdArzw6cpdKoQTd+N5vuMpRwuP6NzcCcG/Lu6dAlSCBXG1swp/N6dXFvAHE7MtRnwexEWtc6ttRuFIXFeRYojkKhQuuEG+2dbz7NKvYzs4NWdy5cayVaUYWq2d0sI51+VuTerAuB1pkZxgGgTpufz9bx+PjpUXrQNff69vKz1uF2VMYZDZpAAAgAElEQVTUiU0lyJAk8fDhw+jXrx86duyIvLw87Nq1yyXCKZGz1fWYumwrnv1uHy7UKi/FuultObO/MqIWnLV6jMpoC8Cz3rrUZAVQGwwYPDAohVfb/MSW0G/bSj2x8Y3P9tztnxzVlfc2QTovRAexKyRSwHfKnxM+pUHs3Ayb5xTLEgWb/6x5duk5n+/maElRMjaVoJUrVwIAZs+ejUcffRQHDx7EokWLcM8990gmUF1dHUpKSpCXl4esrCysXr2as+3GjRtRXFyMvn37Ii0tDa+99ppkchkgRJ9t9K31/+Cn/adZ29hK+86FFDex1gN9jdzpkN3pWJSEkLpetpQlMaZXne3Dcho9Kz4Mz4xuUZqev76bU/3b48dZzjuGm0IIrN52suLDbG4jVgJFoRgKvV6su4qGxias2nJUFjkozsFrOuz8+fMoKioCAHTpYrv8g7OUlpaCEILNmzdj5cqVGDduHE6dOmXV7vDhw5g5cyY++ugj/Pzzz/jkk09w5MgRSWUDzPNL/O+PE6xtvr27j+RycMGg5S3U3R6iQvymlOT/4yjUGiQdltcHn6uFZebJRmN+chAQ0e9T0/ukTbCf+f6aD+KVsexFW53Zl1A4a6lZ/LDLbs7h3We7sJbjnT+ysyNiCYJhgPs+3oU1e81fiNftO409x53PRUWRHptK0ObNm5GcnIzt27fj3LlzAICmpiZcvCjNyW1qasKiRYswdepUAEBycjIyMjKwYsUKq7YvvfQSbrrpJkRE6Ivvde7cGS+++KIkcpnt98f9dttwzY/bQqyHthIVHz7Pg4hW9uf1HZlWS23DHr2jdJR4Ht0ZV0Rocl2+zii7hBAktm4pzMzHysUwjHGKWW6c/dUNp80wfnppWh5pk/LjneydlwSsv/j3f57Ern+rzJaduliHA6cuuUAmihBsKkGXL1/GgQMHcPHiRYSF6U2StbW1ePvttyUR5tChQ6isrERqaqpxWVpaGrZt22bVdu3atbhy5QqGDBmCXr164b777kNdXZ0kchm40tDEGS1g4O6iDg71LcUYHKCiWjdzBqfab+QAQTYcLikUA1z3n9TGOLOXHwttiO+Q8NiwNPEEYkEqHyM25U9txs+z1fU4du4yPv/9XwDA7R9sN66rumzuM/rTvtNYvOGwS+Wj2EdwdFhAQACys7OlkMU47RUc3BJuGRISgtOnrX1vjhw5grfeegvvv/8+ysrKsGfPHsycOVMSuf46ro8Ge/77fXbbJjnrKCgCDIAQf2/Z5sop0tA+TDllB5QMr6SaMj9uX5+QIUo/lhYsU6Xq1XHiTHlNL0gUpR9LIlr5oE/HCKvlpoe0a95AXn0xDDC0Wxu7iuMjQzsJkJAf+05ewuo9+mfXt7tPGqPGnmt+XgiJIqMRZ67HoRD5kSNHii2HGZY3Nts0SH19PSZMmICIiAh4e3tj5syZWLp0KZqa2BMYMgxj9ldaWspbnqELNwAAFtnR4n28NOjIkTzNVRh+OyX5k/CRxVXTP5/wKGZLcQ/8fdiVIYOicP81KSbL7CPklnJIyRL5JhjZ3fkpLynzWcWFB2D2QHbH89Ro/Tga7MffkpvWJggv2/F1coVl+K5VO4zn/+CpS5iyZAuv7RqbCG5ethW7jlXZb+yhlJaWmj3HxcDmfMnTTz9ttYwQgr1794qyc0siIyMBAFVVVayfTQkNDUVUVEuV4djYWNTX1+Ps2bOs7V2RnG7XYwPhxzHw2kMKJYDagazJthNpokT0WYD1n/ultMa6/WfkFUgleHFYQg0PKMFjgpNjiN0pNxv9v3VjJmas+J1FJMvQfPHHOZeOI83if3eP/eCSDQ/2M7N2+2g1SG8XIpVkvDl/+QoCfPUO2vUNTaisvsJru//uqEDZ/jO4vdAxlwpPoLS01MyAIYYiZNMS9MorryA2NhZt27Y1/sXGxsLXV5rkVElJSQgLC8P+/S3Ox3v27GGdfuvevbvZNNmZM2fg4+OD8PBwq7auwlEFSEy6tg1W3Ly6YWAuHc7tuyC2jupoBl+lEBvaEuUSYxLhs2iKNFPRngYDxvyaEPktxLTvCbnSFvcU8iBw5DZw9a3DgN8xxYb6G6PfchLC0Lmt/QKtrqChieDf87XG7wwDPP71X3a3M9SPpEERrsWmEpSeno7+/ftjypQpZn85OfxDFgUJo9Fg2rRpWLx4MQDg4MGD2LlzJyZOnIi9e/eiuLgYjY2NAIBp06Zh1apVuHz5MgBg8eLFuPHGG6HVyqOImIZmOoJY0WG9OljPsSuFPgoq56F0/FiK7tLBkR98fyc+7Wwpz6ab22oX0Zw3zLSN4XPX2GDjeocEhLk1iM844o6X0TWdo9Ezyf7Y5wpfsC2Hz2HjP5Vmy5b8eoSz/fGqWlxtbDKW/uA6P7TUjjTYVIJ+/PFHxMTEWC1/7733JBPIkCcoLy8P48ePx6pVqxAdHY0LFy5g3759uHpV73E/duxY3HjjjcjLy0OvXr2g1WrxyiuviC7Pr3+f5dVuWh/nnAclecAp7KnpSmksxwuF/RSKs9Z5IqbWBvF9glpoZaOMzYODUtGpjXgWDLVPh8ntsC42lorvh1uPWbW5Zdk27D/ZEjrPNVYVv7Te+LmpiaC+oVE0OT0ZxcVQ63Q6LF261Gp5Xl4eKioqzJbNmTMHc+bMEV0GQggYhsGSXw/j1bUHAQDHzl3mbJ8SpZxcNIY58pToQLMpFTlx5QsMfVmiGHBF7h8DQi87MUUTcpxdYoIF5TGT43aS8rzJkUTVcDjl52qMyw6dqUZia/NI4mUbj9js59CZlu2//fMEftp3Gi+NEScC0JMRFB02e/ZsqeRQDN/uPoFnv9OHNj7+9R7jQ3XyYm4P//bhzocui31rTsqLQ78UawdxT0Fplh+K6+CjCBMQ82kkEafGAOCGHrE2X0Is+wrx9xZF4bD3kL+tMElwcIBLlUmJtS45LU1vrz8EQJ9NuuhFvVVn9V8njeu3lZ9v/mT/925oJGiwk7OOwg9BStCmTZukkkMxVNc1oLLmCuqu6k2NhiKph8/WcG4TKUIVY1cONK4iOSoQ0/pKk2PEHu73a1LEht8Dl7uRrVu2XZg/2oXqX44ctT7w3UpqXxFX+6I4OxRumVts9n1anwTnOhSZsuaak8s3l+PW5dvRRAgamqwV8jFvcz9vpSi34qk4lCfIXdl5rAoHT18CIRCU3lyjsKsxqXUrXmUopKZ1oC/6pUS69N2LvhtR+MCAMXuIOKKo2NINQvy9kdKc68ZgfTDbHx/Lk2CJKAAQGagzS4/A9oLprXXNmF2yZAv+rDAvM1XfoM9lt6Q579y+k5cw6o1fjes/267PPr3l8DkcPHUJO46eBxvKeuqoF6oEmfD1ruP44a9TsphMxbyg7x2QjPwk+VIFSIU7Wsso6kGIRaRzTDBKR5gX8DRsnpMQZu4w66RcpvdFvAhT87b6lxqxRl57maZfdJEvTWWNdY4gg3P0IY7ZhQ9+a6lGv/7AGXzDUqib+j6KB1WCTGhsIjh67jI+/71C0EUW4OsFnZf8OYLcHSEPIXccI6gK6DhBNiK0uCgwSenw9LVdrdab6gZ8r7fYUD88OlSfL2to1zYtfbG0deR8p4oYaSYXYlznpnUT2cqodI+VP6kiH5oIQaNFKY0J724GIebK6Rvr/jZrs/vfC2hq0m9LS3HYRpAS5O55CpaaeOcLOdL7r0nBDVmxTu3bnY0cbn7ZSAP90RzGWFnc5J5qG2ptIbH3E5uubhPCHVEV6u+Nkp7x/AWE9f1OeMgDAAEWCVk7xwQjzF+fZyjE31t0RVkpl6Gvl2Pv64efGYK7izsav0cHKyNili/Hq+qwdOMRfL3ruHHZxn8qQaBXFg3P5AU/7MeqLUdx2wp9AdfRb/6Kq01NeOnH/Vi55ShLzxQD1DFaBLQa5+uYiGFyvq0wyek+1A5j8T+FwgfOkhY8pq0iWvkKCju37FcIGotSIA8NTkXX2GCO1uqD60WbbyFVSxiGMSutUZDcGm1D/BT/0vnLQX1pnKrL+um013/SW3oMkcsAAAYoWFBm/Lrxn0p896c+2owQvZ9b3dUmfNZc4Z6N0xfrsOXwOZGlVxeCbMSHDx/Gxo0bcfKk/oeOjo5Gfn4+EhPliQASk3X7zSvV37x0q0ySOM6Dg1LlFoEDGXysGAbjc9q5fL+iofRRWuWIEV1jcKYWlETRrFQH60c7HTi1WhWwnRcdy5SW0D6VYtXiw6T39ClZrjZPZe1vDtQxWISe/nYvClNa46hJ/jpDRDOgvw4YBtAwwI6j3AVZdxyrwifb/kVOgvpqKooFLyXowoULmDJlCr766iv4+/sjNDQUAHD+/HnU1tZi5MiRWLJkCYKD1ftGctMSc6XnHItDG8U5nLV2OTOGqWkApLgG02uCyyfH18TXT0ieIFsYbwM716RS9GClyOEMpufu4SGdEN6qpVRJbKifWa0vJdFokQvIcC7O1VzBVYt1P+451bJdE8HcL3YjlKskC8UIr+mw6dOnIy4uDgcPHkR1dTWOHTuGY8eOobq6GgcOHEBcXBxuvfVWqWWlUHjDalZ3g8GcYh8+Ci+fkHgC2JxqkiqK1NSHRShqv8Rd8a4ytFsb+JvU5usRF+qCvTpGQ1MT5zpTPyFT3l7/DwDgnzM1xvQtL/ywn7Wt2q8XMeClBB09ehSvvvoqkpKsfU6SkpLw8ssv49gx65ooFIq8KPcWd8QyFeDrpeAjUg7GvDwCtrFl7TBEiY3KaIv0ds5HFRnPPcc+k1oHAJC++jwfXJ0uJMDHC14a1wYt2/Ol3Dq3v4sksWbvCfN8dXwsVs80+w1pNQwM7lB/Hr/A2Z4QIigvnrvB62qrr69HdXU15/rq6mrU19eLJhTF86CzVdwYHtBf3N5LXkFUhr1rSug0T1iAD9oEtThA802waDbtZlZ63rCeXdLIQGHO1gE+Xqqfulpcko1uCnP0bi1CRQBHqaiyr/Qs31zOue6NdXqrkCFKfu4XuzHh3c04WnkZZftPg0CfvHHE6xvEEFeV8PIJmjRpEtLT01FSUoL09HQzn6A//vgDS5cuxd133y2poBSKEJSaWDHYz9tYikUoPl4amvNDROyGyLM0kOqyYhiGVwqSrrHB2HSoklX2n+4rMPNhEkUuF9sefRwMhfdkHv3vn+wrzCrY67+s23caxy/U4b87K7D/5CWM6B5jvKY/+K0cE3PjWIu7ujO8rriZM2fi2WefxXfffYfrr78eBQUFKCgowPXXX4///e9/eOaZZ6gSJBEdItV/MVKnZIor4efv4/xFyTdZIqfiZLo9z5vkiZGdkRARwLpObAUIkLfgqNT06RghtwiSwnbujl+oAwC89OMB7D15Eb+Xn0djE4GGYTD3C70yVfTieiz99bBLZZUT3mr3DTfcgI0bN+Ly5cuoqKhARUUFampqsGnTJtxwww1SykihCMbdE3tS7CPIJ8gJi4fDWzp0ier39tltPR3dq4N7dD+WT83lVCjdAdOirGzD4aEzNXj750NoIgTaZk196xF9zqBd/3L7ELkbgm2PXl5eaNOmDdq0aQNvb28pZFI0/5mYKbcIHovQwVihM2KCCdJ53n0mBramRB1VetgeJnx1GVuKOcMwGJYeg4z29hyv9X24KqLJnd4l/nuHtU+daSHVEH/9fZZp9xyoA9P8QLYsek1NLWPlsea8QwyA+oZGfPBbOb7cWYHKavf1+RVtAvabb74RqytVMW94mtwiKB6xxlGl+vlIzZd3UodoZ2FTQLqbPuwcuLTs5Rlia2eGxUYRrXzRNoSrrAPD8okihO4ckX0GhdgQBeaOfkm7jnFbdn4/eh4X6xrMFzJA1eWrWPDDfry34TCOKTSPkhiIdrZfeuklsbpSNJYD0E29EmSRwxMRMsWlVIXJkWk6b637DcquwN5vber8aetqKekZj9ToQH07Jy4rq2uyWbwQfx98OiPfztbymGQUehtJgjvfZ9X1DZwJgE2nzQznW8MwyH16LRiY1yizZNM/lSJL6np4RYcVFRXZbbNz506nhVEDbmQddjkeNJ5SFAabUsx1Pb46rjvu+bBlPOuXGsnRpxiS6fO5GBQye8p7QkQrfHBLnjg75oE7TYfZ4teHWp5xro6IcxWZ83+028Zwvj/drq83xjD6uphr957G/G/24PPbe+HI2Ro0NDWhQ2Qgxr+7GUeeHYoLtVcR7KfOaXteSpBWq0V0dDQ6dOjA2ebIkSNiyURxM+RwUrbc593FHREfbl1JnOK+mCoU9h5rpm1Hdm9rpgQ5LwfXCutF3PeKvrFWwwgu1uo0HmAOspyG9PXSoL6BO1uzp8BAf/pfX6cv4FpRVYv7PtmF5OhAPH1tV2O79MdX48izQ2WS0jl4KUHvvvsupk+fjuXLl3O2KSsrE0sml7Pn+EW5ReDEnYYfm2HEEuzLtE+lpMZX6jSdu2N57REXBn8LfQcwvUb6pbQ29CKeQBS7RAXp4KVhcOhsjdyiuBTL4amy5goqTabRej37EwDA39fL2N6Qu+z9TUcwOT/e6rPS4TUJGh8fj4ceeggVFRWcbZYsWSKaUK5myMJf5BbBKVKiAuUWwSauHL5NKxLQx4ZnIoeeae9aK0hu7ZBFdMlNOY4JZMKTo7o4vO1dRR0xvFsbp2WgqINqSwdpDn4+cAYfbjkKQoBd/+qj0B778i/jetPPbLy/6QjOK6RIOW9PsH79+qFt27ac6+Pj48WQR/HQ93h1oHqLi6c4Y0iAWAVUxSTFacdqx+W9MS/O4W2jg3UI8XffSuSW18orY7vLI4hC2HOC/6zI4Uq9leza/2w0LjPNaH/WRlj9so1HbK53JbymwyorKxEeHs6rQyFt5abuaiN03sKyrNJHE0VO6PUnDc6qREK2d+cszGrEVDHNSQgDAMSF++P57/d73HSYI+9eUUG+OHVRr9A8/e1ehLfS11qrunwFEa2s667d8cHvinpJ5WUJEpIRWk3Zo/OfWYsjCr/I3WG4dKVRw3BrqT5jtIIGCbViVqtUxuvBOEXLMNTAp3BiQvwwKqMtMtqHIrG1+2aTFhN/nxZbyrnLV/Dc9/uav7HXw/vf7hMukowfvCxBJ06cwPLly3kNJLaqzSuNmiuNeGH1frnFoEiEEtUInbcGF9w375hqcaXOaWsqjupIriU7IQytfNkfg09d2xUf/HbUxRLJy5VG/hFxH245BsB8nP29/LzxM8MABQvKMHtgMlr5eqG4U5RxXd3VRrN9nbxQBw0DRAa5OPIRPJWg2tpaPPbYY7w6bNeunVMCuZozl8znJXPiw7DlyDn4eWtRe7VRJqk8D7EGf6U/RKQockkxR6zpJrGsNoLLvYizW5fx0ph0uUVwGNMwbwrw+e/cwU+WXKi9CsB8zLXMLH303GXsOX4Rwf7eZkrQv+drMXThBmNY/fLNR+DnrcWdRR0dF95BeClB7poDSMMANVfMveFDA/QJnz6ZkY9hr22w2sbVA5TaBkQ2LC2I7cL8cOyctOYQJc05myLGA1rDAE0u0vYy2oeY1SBSK/auB1dcLlaWdKVr7DwZnRkrtwgUGTls4lLSyDYwMcDz3+/HjqNVeHdyFmc/hOgLuGbFhRrv19s/2I7/TOwhusym8FKC3BUGjFlIoGltGcvqwkO6RuPb3SddJps7Yus5o0yVhRLg4x5DBNtUviN+Qs4qS4btCcB60c8f2QX+PvyshVzTOBSKEpjz+W4AwNvrDwEAftxzirXdJ9uOGaeIJ7y7GY+P6IL9Jy9ieHqMS5657lsshQcaBjhSedn4vU2wDmEB1t7sphiGzXSOYnxyoFCjhyBEK7Jq6I96oHoshgHVVfeFvWuNLVEjwK74hwX48I5Y/f3RAbzaUShysOXwOV7t7v/0D+Pnq40ED3+xG78cPIuLdfrptsrqemw7wq8vR/BYJWjbkXO42mg+PBECBDS/hVnVOrQYydxA73AZciRLBNxDObREwwBzBndy2f7c8Td0BsM4sGJqrqDtGOP/THM/zt8V7ljtXKkYCuhSnIPLGrRyy1GzsduQgwgAej33E1776W/JZPLYu+juVTusPOGz4kPRu2MEdN4aqwiOJsL99kahiEXrQNuWSIZhMK1vIoo4inqKjVJ9q8TAvLaYsOM0/Vn4/EbE+L9pxW73/W3djdcnZLIuf/EG9TqFy8G097exLj9XcwUv/XjA+N30HaHuqrQ13ERTgsrLy8XqCnV1dSgpKUFeXh6ysrKwevVqm+2vXr2Kjh07orS0lFf/VxqacPxCndXyjPYhKEyJZI0KCfFz36ypkmPpD0pnqpxmcUm2S/bjzo9pgzWGa2qby4k9TmAhXqMViCo9bgc9pc5h8Btio95E+Vl/4Izx86P//RNXG5tEq/kpWAkaO3as1TJCCK6//npRBAKA0tJSEEKwefNmrFy5EuPGjcOpU+xmNAB45513cPr0ad79/3v+MutyWw/nYH9vXNM5irsBRVHQB47nIjQCr22IzuphZmu664FrUgXK0/y/yQBDiHsrmO4H+zVFX+icY9UW7jxMt33wu9WyKw1NWL65HDX1DaLV/BSsBH333Xd44YUXjN+PHz+O4uJi/P67tcCO0NTUhEWLFmHq1KkAgOTkZGRkZGDFihWs7aurq/Hxxx9jxIgR/PfBcuHe1CseXdoGs7Yf0jWad99i4y7P8rAAakmjuAYpH0yO3o+M2Wc3uak9AHsReATA5HzHa7NRhJH8yHcAgO5P/Chan4KVoDvvvBPt27fHddddh/feew/dunVDYmIiCgsLRRHo0KFDqKysRGpqy9tWWloatm1jn0t88cUXcdddd0Gr5Z+Ejs0pMTEiwBiVQWA+2Hlr9T8THbwco3v7EGg10v52txcmmX2n0WEUvhju6+HpMcZltL4XBeAXgUeHGnUjWAmaPn06hg0bhqamJkybNg233HILFi1ahLVr14oikGHaKzi4xSoTEhLCOt115swZ/Pzzz7ym4hiGMf4lRwehasMHxnXjstvxGvKKOkUiobmejLtYaFzB8qm5ZrlPuH67oV3bOLyP2FBrPw0lKq1qGzDVeJ23yGxbeKkVcy5I8z9bv+2r48yrmdM6VvJgiMBT233rjryy5gCqNnyA8ueGGf/EQLAS1L9/f2RkZOD8+fPYv38/Ghoa8MQTT2Dy5MmiCGTA0qeD7c1+/vz5mDt3Lq/+CCHGv7gHv0FI74nGdZ3aBCEm2M+kMXsfY7LaITU6qLk/XrulNEN/L/t8MiNfbhE8ijam97wFXAo023XMx+poLKIK04g0bkZ2b2v2PTmKhmg7w20WlmKxoBZn1/HKmoMI6T0RcQ9+Y/wTA8FKUHl5OUpKSrBu3Tp07NgRL7zwAnJycvDtt9+KIlBkpD70t6qqJVV/VVWVcbmBQ4cO4fDhwygqKnJ6n1N6xqN/GrfTswpfhhWJMVuuBOOG5Vu1Gh2j48KERR1R2HHUAujMA43v5WaYZqPPTtfy4CBhzuxCoFOn6kZw3vXx48djzpw5ZssGDRpkdGR2lqSkJISFhWH//v1GxWfPnj0YMmSIWbv169fj+PHjRl+kffv2QafToaysDEuWLEFCQoJV3/UNjcaibwZ2sMz5bn2kP3wVkoiMDpYUivSwKc1s9x6XsmPvPmXbTKiivnKasASNFPGwdXrpGK1uBD/ply1bxrr8+PHjTgsDABqNBtOmTcPixYsBAAcPHsTOnTsxceJE7N27F8XFxWhsbMRNN92E7du3o6ysDGVlZRg0aBBKSkpQVlbGqgABwIaDZ/HQZ+Z5CUJZopaC/bxVZ0l47jr3r4bMZ6wxDEhqNFOrT2J1MrhLGxSmtLZa7mi2cSFjRUvGaP79G+iZFCF8I4qkCD2NvTvQc6g0BFuCuKafdu7c6bQwBkpLSzFjxgzk5eWhoaEBq1atQnR0NI4cOYJ9+/bh6tWrxmiwK1euYODAgUZL0B9//IHPP/+ctd8HP9uNs9X1osnpCviOr2Oz20sriMIR8tauZJQmstLkEYPulskRGXmO0x1/W3clMYLDMZ2AV8HbQJ0XLtU1YEx2O2z4+6zI0lGcQbASdOLECTz00EPG7+fOncPXX3+Ne++9VzShdDodli5darU8Ly8PFRUVZst8fHxQVlbGq1+xFKDS4Wn4cpc4li+KuKhR8aGIi8M+GiabaTUM70Kmdrtlc6YWpWeKq/DSOuceYUizokYLtbsjWAl65plnMGrUKLNld955JyZNmiSaUFJgWpfEwL75gxzqq6RXAlWCREbswUFt05kG6BApHs5cArkJYchoz15OQwwITRlNoSgCweqtpQIEAN7e3vjnn39EEUgqFq49aPZ9/sjOor3pUYTB+XDy8IeChx++aBj8bvj+nlxOy75e1uPDmKx26GqRWd6eAq9SfZzCwctjW4qm0sgw9SPYEvTEE0+Yfa+rq8P27dsRHh4umlCuYFJ+vNwiUCxxIMKGsyuFmp1v7ZuIx778S24x3BpjGLrpMhuXgyFDPJ8HWmGKPmK1oqrWOMUB8FN06EuXe2B6LTkzzLQP88fRc+x1LCmuQ7Al6M0338Thw4eNf5WVlRgwYAA+/PBDKeRTLBr6eicqUkxfKfEMTVaZ8q3WaUUu7ujHnjRPaG6hvMQwvD2ph6Btfn6gH4L9vAE0K16KvEIpQnFUEbIs9UORB8GWoOnTp6O0tFQCUZQL24Ng2c056DLvBxmkUSeuMhur9bHy/HXd5BbBI7ifpQK84ZoR8jBjGHPHaZuWpuZ1pu3bBOuQHNWK/w4pstMvpTXW7T9DpzfdDMGWIC4F6P7773dWFlVhr7qwWLjj26LUM1VFqZEYYVF2QOmMyW4ntwhujb0HlzMPNl7WMosmfTq2xvQCaglQE48MS7NapsxJd4oQeD3Jb775Zrttvv/+eyxYsMBpgVxBP5ZEabYYkR7jkgKGH92ah7HvbJZ8P3JgqsxJqQQRAItLsrH01xvqHYIAACAASURBVMM4fLZBuh2pmOz4UGw9cl5uMSTHVO/gc81JqpyL5EdCkY+2IX5IiAiw6xMUHaTDyYt1dvsz1pNj6DUhJ7yUoHXr1qGkpMRmG51OJ4Y8LiGluQgqX6KDdYgLl14Jyk1Ul3O5kokJ4S6OqVRcNQ56yoCbFR+Kb3eftNuOYaS1t5oaiuhUinrReWuRHhtstdzydrJ1jpNat0yBGpIsMix9UFwHLyVoxowZePDBB222UZMS9NBg6YrpUZTBwM7RcovgEEp7RipNHiG8Oi4D3+7+zm47QggCfL0Q1sq6hA4f7EUiihVNRFEGziiyXZrTK3hrGYxIj8E9H4pXaYHiGLx8gtq0aYP3338fJ06c4GxjT0miUAzQt2F5ob+/NX2TW+OpUV0c3t7ubyqxNeiuog7id0qxi6MBHwzDuF3kpVrhpQQ988wzAPRvPEePHjX+UShKwVPHk49uzZNbBFXA9bDy9dIYH0ZSPJQyJcw6bcrsgSku2Q/FvmXP8iq6oUcsSyeiiqQKZig0EICXEhQdHY3JkycjJiYG8+bNQ15enmrD5NmqRyuZ5643D53Op35DNvG06YYEiR32NRp5tEvOgpUC4CP589dLm5rg89t7Sdo/xbWYKsqGj/bGnAU3pFstM1XKPcUi5OPlXP01qeAllelJWrJkCVJSUrB48WLJhJKSpTflCN5GzmvUtOK1hgEiAn2t2gTpvF0pEkUhvDkx06HthCiKA9KiHNqHs2hlUr4oFFsY/L+6xQbDi+MadVSpcfR+pjiHQ6oZ20meO3eu08JQbMP17Fo5LdelcjgL2xih81bmW4KSGdy1jUPbCTGWjUiPcWgfzjJW5LxJtuK/xFC3bCmWJT3jcVdRR15tKcqHYYDkqEAwYEAA3DswWZR+Hb2f1YwSEobyig67cOECNm3aZNSCLb8D+jxBTz31lDRSisRz13WVWwSn4Co87aVVvgJhav7Njg/DsXMV5uvpg8ExHPjdhNRV85XJhJ2f5Pi0b3SQDkcqzWsySZmx3N6bv59PS6ZoD5n5cGsMt49Wo1etHUmca3oLevIlkRUfhgOnqmWVgdfZ27FjB3r37m02ePbqZT7XrfR5zc1zihEd7FgYf4II/gmUlrdxzgerCJdQK51rMnlTlElWXCiuNhGgkhampIiP6XPu4SGpGCfAYtm/U8vUsie+8ylVQ+D1xMjNzbVZIJUQgvHjx4smlBQ4qgAxDODrRas/qwEfrQa3KTQCQUko/YXFWTpFB2LXsSqXHacQyxrFffD11hqt8Nf3iMWn2/8FAGhM3vGeNEm7oPaZCGdR6l3Cy9b93HPPIS4ujvMvPj4ezz33nNSyysKueQPlFsEI25ge4KNGBc3xh1NbW5mgGTrdIAY94kLlFsFhGAZ4VoZitHyvO6ovqRszhdfk4wsmEWDf3t3H+FlDByQA+oLBSr34eSlBffv2FaWNGlFS5JXl7RQbqr7SEGLDFaGhRpQyRnx2W0+5RRD9t+DqT47fnD4X1Y+tUxho8szg8kUzVaaUcD2IkZLCFlyHqIQxT/ketRQjltdLanSgLHI4jzhXPsO4Xyi1EgZET8PVv7kSBn6KY5hOsTrjbE8vAT1KGO+oEqQyTG8efx8vrLuvUC5RHEaMh8AtvROs3y7oyGJG17bWxR7VgqODY68OEaL2JzZKkYMiDZ/MyDf7zmesG+oBofEMwyh2eKZKkIqwHD8ZBogMUk/hWgNiPAhmDWDPzeFxDxkbxxvJkljT3ZnZv+W64HMptAnWoVObIOkEorgtbLmnsuPDAAB/lOp9Sfk8+J8e3RV+3jL7drpg3GRTCGuvNEq/YzvQeGIVoVRNmqJMPE4htAHXG3luYjhyRShFQ6e4PINrOkcjISIAm/6ptDkdFticO0jLcRMq7XqReqhgGPbpw82HKiXes32oEqQyTC9Wpd1IYhPRygcLx2XILYaycfNrgEJREoO6RAPg/4IhduZztcL1eynBp5NOh9ngmdGenddBTOxVXmbj5bHd0ZPDx8OSnIQwB6SiKBE1nktqdaOYYnCg5vOQd/eXWYC7bI0Soq+pEmSD8Tnt5RbBJmobeKWUN6m1/DVo1IIaLhtPeDBQKAYMY+PojLbyCiIRDMN+T8eF+7teGAuoEqQCpvdNlFsEUbGlDD0xsrPrBFEBH92aJ7cIqmTDg/3MrjOl5NR6Z1IW0mKoI7an06cju4U7TwT/NEeQOrs6A+XO3FMliKI65gxOlVsElyGG066ruZcjcs+VxIa2vGGOSI/B4pJsGaVpITpYB28VFDym2MdRa2XrQF8sn5rLui6jfYgTEjmO9I7RjGKtu/RuVAOM2X9uwfQCc+tWRvsQTOA5/RjRyj1Dv3XeGmTFi+kPwxj7NUXqsUiJ16m710ujuBYuHxc+bJ5TzLmuY5Q8CXDv6NdB0v5pxmgKxYKIAHNFJjbUnzPRnYEgP9N09O7D/OYpwBB/H7wxIVP0/u8p1ltmoptzSkld8FMMfUNMnYXqPxQlYeksrYSxLMWk+oCPCy2VWq38NydVgtSAEu4SkdDfYI5d+KZKEOA+b/eT8uNdsh9vL9f8XhoFhL26ki5tg+HnQ7ONUJRJWICPoPYdIiUIMuHIE6SE2o+KU4Lq6upQUlKCvLw8ZGVlYfXq1aztvv76awwcOBBFRUXIzMzEwoULXSyp67G8hOS/fMQhqXVL8T53OSYl4irTc/9OUU73IYasrlKSXx2XgbYhynC8psjHO5N6yC2CYmEAqwfY7YVJ0Fjco3K4Oiju9aW0tBSEEGzevBkHDhxAXl4e9u7di6go84F11qxZ+OSTT5CRkYGTJ0+ic+fOiIuLw8iRI10m687HBrhsXwbcxPghCgyA8AAfHL9QJ7coisbRQo93FTnmJ+DsJUovcYoasLyrBnaO5mzbMykc3l7cNgd3v+bZaoc1EmKlBMmBoixBTU1NWLRoEaZOnQoASE5ORkZGBlasWGHV9o477kBGhj6bcHR0NPr168dpNZKKEH9hZkZKywPZ9Np39I3d10uDnx/oJ4ZYHgVflSjYz9ttphwpFDHhc1vMHdLJ+Pn1CZmciQG9NAwy2oeKJZoVhBBkx9vuX2orsYax9kUkRL/cFDmGG0UpQYcOHUJlZSVSU1tCoNPS0rBt2zartrNmzTL7XldXh9atW0suo1JIU1nRxwXXp6N1c0HPnIRQpEaLI78XDTcGANxd3JFznTOhvI5A9SYKBZjGM7+bzluLFbe0hMwvuzlHdFk+mdETb/OcrpNCH2KLpmtsIlZO4nIMHYqaDjt16hQAIDg42LgsJCQEe/bssbndxYsXsXXrVrzxxhuSykdxnPyklnw3Q7vGWL0BWNKbZ7kMAJiYq+zM3q7Al8XU7qwyMrK7Y9lrlRD2aiBXhSU4KMpnSs94BPtJU/KhIFncl3nD7ZgiU/g9F15aBj42pghdhfwSsGBpgrcX0vvggw/iscceQ1xcnM0+Tf9KS0vFEJUiIXwfpk9dS2u8mRITrLO5Xm5DzbZH+rtsX+Ny2sl+vBT3o1ObIMSowBn+v3f0klsEAOxlM2YPSMGNedzPbDaqNnyA8ueGGf/EQFGWoMjISABAVVUV62c23nnnHXh7e+OOO+6w2bfUuVHkQI3+GnxPw5s3ZsLfIuy4a9tgVNc1SCCVZ6CUW4BPBAjbpR0b6od/z9dKIBGF4p6EW4THazUMGpu4B4JZ/ZPx/V8nJZHFcq8+Xhq7MwKWhPSeiJDeE43fxVCEFGUJSkpKQlhYGPbv329ctmfPHmRns6e8/+KLL7Bu3Tq8+uqrAICDBw+6RE65uCaNO/pATfDR3QJ13lbzxY8OS0NKdCDn9rMGJEOrQsXQ1ahBeRaisNmallCK4kehiEX7MP5FRy0tMPZSOUg5NPC5F+W4XRWlBGk0GkybNg2LFy8GoFdqdu7ciYkTJ2Lv3r0oLi5GY2MjAODnn3/GwoUL8frrr6OmpgbV1dV48skn5RRfMgwXRtfYYJvt1ISzD2ItwyA+IsBs2W2FSapwlN7woPCItlB/jge98vUZwfj5aB3abnh6DOtyNSh9FApfIgUELJhe+7ZuA1tpNKT021HCS4rinhiGPEF5eXkYP348Vq1ahejoaFy4cAH79u3D1atXAQDjx49HWVkZIiIiEBgYiMDAQJSXl8ssvbS445QeXywdf/18tPjmrt4ySeMcpsU9+bJ1LrsfTetWvpg9IBmtWaaYuOv1KLtsBlfkoz2xh3RxD0sphSIFsaH++Hh6vs02xZ0iBY0PmTIVfBUTRfkEAYBOp8PSpUutlufl5aGiosL43fSzJ6H2d1pHH7+75g20WuZJb/hcFi6GYXBnUQcwDIMHPvvDxVK1cFdRB7z209+y7R8ABndtY7UsUOctuGwAhaIUQvy9UXX5qtkyR4c9rYZBtJ2gic4x/GcbGMa6Dhp3W4ZX0lY53vMVZwmicMPXtKl0bIrOcWA6b8emSNTIjIIkQe3tKYOGNzspr5nZA1Ok65wHsaHsvg4D0qLwyNBOrOvYkCJHC4XCRf9O3EE/ALcFmC8MnLP8PjYszWbfNnysrVDqRAZVglSG4Toa0rUNHrVxgVLUwbQ+CVbLHhqcytKSP55kIYsK0k8DLpqSxdlGyO8hdo4WCsUWi6awB/0Y8HbSx5FhnHc2/s+ETNblGobhrWApeUSiSpAKYLuAfL00bmPmV/INIjVSOnIbnv0GJcAwXkn9QqbVaNDNxIk/3QmH/r5UKaFQzGDLvuwsXLoMwwA9O4Szr4TtseT567o5JMuHt+Y5tJ2jUCVIpaj1bd84NSOzHEqgk8pKn/AlISIAn87o6VQfHSJbAQDe5zE9dWNee1bHcEv+d7c6HekpFEfx93be7ddLw+WPaHs7DQ9/IWsliiAvkVvpkgKqBFFcjloVOLEZ1s3akVcK5IgqdCaslmGsfcC+u6cPa9sVU3Mxq38ywnkoQUKcPikUt4CBU6ZfQrjvZYZhbPr5eGnMfVgJIehlw6okF1QJUijTeRbfo1CE4nKdiEPpFaIocVnNOkYFqiI3FIXiDPOGm/h/CniHlPJ9k0u/ahfmh84xQdB56+9L0xxnzo49Uvjs0dFDocwakGxzvbvYUnomCXsz8PXS8A7LVDPdY0PQL8V25Ig9RqS3RZe2yplyMy10Wzo8Dbses057YAtPzpNF8Wxu6mUdQGFZEsOUsVntAADeGg16d+RfjBowV1RsKVEMA1athgGDV8dloGdzEeytc/sbHbSVOAmguDxBFPcmLMAXhBCHp8SeGNnFzMzqrnSNDUZXODd90zrQF7Eh/ooJTTU9byUsgzoXQzmmDd0lMIBCsYVlolgDfIZQPx8t3ryxh8gStcA2tDCM+b1usNRmxYcho30Ifv27smV7iw5oniCKXdT++N88p8ip6QsfLw0vhzuKnhfGpGNgZ/NMykp8GwOAO/qx50d6gyNEd8vDxVKKQ6EoArZEsXq4b2Rn7nG2bT+7zTrQgStKjUuRiQ/3R24Ct+X/7UnWylpeYph535xbOw5VghTG/FFdrJZZWk3G57RX7IPMHgYFSK3yq41Wvl6iTB/G2Mk0Kwb3XyMsPxL1BaJ4AqZBAqH+3pK/CLMpMZ1jrKfVLYuzWq03+2xf6nYs5YQiA3UYLHE5HDqKKIyiVGs/kLhw84uDyzzqDuQlhNlv5AZwFft0BY4MonOG8M+67Cx3FXV02b4oFLWQ2T4Eb0xkt4q6mj4dIxAfHoAmNp+g5gHGco0jVeQZBpJO5wFUCVIl7uQgmhARgP+Y3NgdowJllMZ1vDY+Q7K+7+lvqUSo63rhKpGS0T4UtxUKKylCobgLn9/eCz5aDVr5ut6V19JyPyarHRJaBxgVmz4mztdsjyc+jtGm695hmRrT9y3+WEaVIJXiLrNJPl4aK8Xn6Wu7yiSNaxmd2RYTctrbbyiQMc2RIVyIPYyU9Iy3uV5I7S5b+PlokRARAEBfsJVC8TR6xIXiLQ4FwRJbdb9sYVro1OZzhgAdo1oJ65vn4GPwY3TFc44qQRTFMSFXfMVAibQJ9kO7MOt5cDVg6mdUOqKzzbbZ8fopzol5cc7tk2nxLJC7YCuFIgcMw/CuJyak6DSXcsKlsxisNq+OyzDKZbmO135N9sBuQTLvzDSb9OsTxLGmUyVIhahrcoMdmjVa3fzBGbHCTbKTU52Pj+yMazPaOtUHheIJSDm85iSE2bXoEFjkG7LTpyHdhWU7y++3FyZhzxPXYOH4DAzrJo5fJVWCVApVIihCSYnWKyHOVqYGgAAZ/BJ03lp4aTVWYbMUiqcxrQ//PFt84HqcMMb/9Z8iWvniKUMEs8BHEJs/T3y4fnpbq2HMrEJWAhjlZODv44URIgaWUCVIYbiT0zNFWSwuyQagn1ZSMx/emi+3CBSKrEzrI25ZJdPHTlSQDpPz2aeu5wxOtRm88ssD/QCw6EccY47plB0h1s00LhirqBKkYB4YZO33YK9oHcVz4XIWFsPyw5cBaVGCt/n6TlrdnUJREuEB+oLEgbqWul/eWgbX9YgFAGS0DzEL6jCoKgYfR8cfUeZKTxeW/ERiQ5UgBXN7IXcEjMpf5ikSMHtgCvqyFBiMazY5t3wX7ozdnqcD97uTswT33TWWVnenUExhu4/NcGL8v6U3+1Tae1PM79302GCMzmT3wYsN9UevDgJrkglqDRx+Zoig8jqOQpUgimyMzmyLyEDpMxF7Eu/fnGO3zfSCJMQ3K0LbH+nPq9/0diGsyx8eIizLMxczCmj+HwqFN07MBoSy1NyLj/BHcScLKy7DOOd7KsAxGtBPh5kmt3eV3ytVglTI+Jz2TlcYVwIRrXxBy4BJw7uTs5DZnl1xMSVQ5412YX4O7aNPxwjc2lcc5eWhwdbKFF8FjULxRCx1hBduSDf7zupozIHptJezmIoVqGsJoLDlxjEwLQq+3hpZZjhoFXkVkuaCeVJX4Oulcam/iifB1zdHyEBpyfKpuQ5vy8VUE1N9eCtf0funUNQAIQTXN/vfcLcx/67zlm4sFTJOhLfyNasiDwBBft7w1nJrOE9d2xUnLtTyqjEmNvQJRJGNe4o7YoqdbMMU6bC0wj3JUrzX1TzqYJZbCsXd4FNfkO2eZas/KQS2KC0hfDw93+ggbehm8ZQsdIttsUyzFUWVK+CHKkEKw5Py/3hpNaJUOKc4hpdGf/vHhvgjJyEMN3JkdJ7VP1m0fUpZM41C8SQYBqz3rD4VBv9xNSbYvl+mLQtNkM58QskwpntpGaMvoZeFxZ+tKKq92mJSQZUghdE2hN0/IyUqEBGtrB3aKBRHWTi+O7w1Gtxd3BELru/G2c60IKsQB+bOMUEY0tX8jY/P2y2F4uk8M7orchNsJwW1tJyYfu/SNghxYeZRoVxseLDI7DubImJrOuzVcRms2wT4ellNmYexOGWb7Ztj+SQnS+7YgvoEKRBfL3PddGi3NuiXEonyyhpcqmuQSSqKmllc0hL+GuLnjesyYzGoSxsAQH5SONdmVgzuEo231v/Dq+3/7u5jt824bNvFXikUTyQ21LmaghNzuZUG04S8iyZnQWNhjScECPbzRte25r6nYvjr2KpnZitR8PxRXbB8c7nT+2eDKkEqIEjnjSCdN6KDdchN5P/AolAMFKW2OEqHBvjgruKOrO1sOVSP6m7bipMssKI0ADx7HbcFikKhcGNqffH30aIgxU5uIRN8vTSob2iChmMuKLF1gMv9NfW+SNQxmkKhyIitZIevjMuAl5ZBTw7L0Q8z+0olFoVCsUG7UH8E8QxzZxgGU3snoJ8Apclbo8FbkzJZ1zkTYSqETm2kiYqmliAF8tWdveFDQ8cpEsDXr+z7mexTWZ1jgjkTMtp6i7MX7kuhUPgTGuCDNybolZLtj/SHtxf/54W9+pQFya3RJtjcN1WjYcysyVIQpPM2luVg47t77E+vOwJVghSIodo3hSI2vz3MLwFhajT3W5dlpAcfLBO5USgUx/HWapDT7DjtSD4tW7NO0/qKW5yVL8H+3rh3gHiRqHyh5gYKxYOgKQkoFPdlDkvmdUsYhgEDBlqNRrTkhAz4h7cbKs0bePrarqLI4CgMsWcbUzkGE72bHyaFQqFQKHa50tAEAoKmJr2DtGV0mFAuX2lAZfUVhPh7i1p+gw9iPN8VaQmqq6tDSUkJ8vLykJWVhdWrV3O2ffHFF9GjRw/06NEDCxYscKGU6qC0tFRuEWSDHrtnQo/dM6HHzg8fLw18vbTw89E6rQABgL+PF9qF+btcARILRVqCHnroIZw4cQLLli3DgQMHkJeXh7179yIqytwx6/vvv8fMmTOxc+dOAED37t3x4osvYujQocY2nm4JYhiGHrsHQo+dHrunQY/d847dLS1BTU1NWLRoEaZOnQoASE5ORkZGBlasWGHV9u2338aECROg0+mg0+kwceJEvPXWW64WmUKhUCgUigpRnBJ06NAhVFZWIjW1xcErLS0N27Zts2q7detWXu0oFAqFQqFQLFFciPypU6cAAMHBwcZlISEh2LNnD2tby3anT59m7deTCpNaQo/dM6HH7pnQY/dMPPnYnUFxliADlieUa86PnngKhUKhUCiOoDhLUGRkJACgqqqK9bNl26qqKuP3qqoqtG5tngrcE53FKBQKhUKh2EdxlqCkpCSEhYVh//79xmV79uxBdna2Vdvs7Gxe7SgUCoVCoVAsUZwSpNFoMG3aNCxevBgAcPDgQezcuRMTJ07E3r17UVxcjMbGRgDAjBkzsGrVKtTV1WHDhg2YP38+OnToYLP/Dz74AD169EBWVhZmz57tFpairVu3okOHDli6dClnm7KyMqSmpqKwsND49+mnn7pOSJG5evUqXnnlFRQWFqKgoAD5+flYu3YtZ3t3Ou9Cjt3dzjsAvPbaayguLkZRURG6du2K1157jbOtO513gP+xu+N5N3Dw4EF4e3ujrKyMs4275o+zd+xLly5F9+7dzc77hg0bXCukyJSUlJgdT2FhIaqrqznbC77niQKpra0lU6ZMIbm5uaRHjx7khx9+IIQQsmnTJhITE0Nqa2uNbV944QWSmJhIQkNDSdu2bcmSJUs4+929ezeJiooiZ86cIY2NjaSwsJC8/vrrUh+OpHz++edk/PjxJDMz0+axr1u3zuZ6tXH48GESHx9PqqqqCCGErF69mgQEBJB///3Xqq27nXchx+5u550QQrp160ZOnTpFCCHkwIEDRKvVkt9++82qnbudd0L4H7s7nncD48aNI/7+/mTdunWs67/77juSkpJCamtrSW1tLUlJSSHffPONa4WUCHvHvmTJEs51amXKlCm82zpyzyvOEgQAOp0OS5cuxebNm7Ft2zYMHDgQAJCXl4eKigrodDpj29mzZ2P9+vU4d+6cXSvQe++9hyFDhiAiIgIajQY333yz6vMKZWdnY+XKlQgM9Kyiq4GBgXjiiSeM0YEDBgyATqfDxo0brdq623kXcuzuyPLly40+gh07dkRoaCiOHDli1c7dzjvA/9jdla1bt6JVq1ZWvp+muGv+OD7H7uk4cs8rUgkSSmxsLK92bHmF/vrrL9TW1kolmuTwPXYA+O9//4t+/fqhT58+ePLJJ9HQ0CChZNISHh6OSZMmGb8TQnDlyhXWAcLdzruQYwfc67wDQLdu3YyfP/vsM7Rq1cr4omSKu513gP+xA+533gFg3rx5mDdvns027po/js+xA3olsKCgAIWFhfjPf/7jAsmkZ8aMGejTpw8GDx6MdevWcbZz5J5XXHSYlLDlFSKE4OzZs2jXrp2MkklPcHAw8vPzMXv2bNTV1WHYsGE4f/48XnzxRblFE4X169cjLi4Offv2tVrn7ufd1rG763nfvXs3xo4di5qaGnz88ccICQmxauOu553Psbvjef/uu+/QuXNnuy9+QvLHqQW+xx4VFYWBAweipKQEZ86cQWFhIRobG3HXXXe5SFLx6dSpE4qKipCdnY2tW7eiX79+2LBhA7p3727V1qF73uGJOgVSUFBgcx68Q4cO5O233zZ+//vvvwkAUl5e7gLppMXesVvy7bffEn9/f9LU1CSdUC6itraW5Ofnk23btrGud+fzbu/YLXGn804IIb///juJiopi9Ytx5/NOiO1jt0Tt572xsZEUFBSQyspKQgghcXFxnL4vXl5eRj9SQghZs2YN0Wg0rhBTEoQcuyX/+c9/SGpqqoTSuZ4xY8aQGTNmsK5z5J53i+kwvrDlFWIYxiPnWNu3b4/Lly/jzJkzcoviFIQQ3HrrrZg1axZ69OjB2sZdzzufY7fEXc67gYyMDAwdOhQLFy60Wueu592ArWO3RO3nfeXKlbjmmmsQFhZmty2f/HFqQsixW9K+fXuUl5dLIJV82DomR+55j1KC2PIKde7cGX5+fjJK5RoWLlyIuro64/dTp07Bx8cH4eHhMkrlPLNnz0ZOTg5uuOEG1NfX4+jRo1Zt3PW88zl2dzvvlZWV+Oyzz8yWBQQEoKamxqqtu513Icfubuf9l19+wTfffGMMkT558iRmzpyJkSNHWrV1t/xxQo79+eefN/t+6tQpxMTEuEpUSRByTA7d8yJbqmTFckrozJkzpHfv3kYz4u7du0l0dDQ5e/YsaWxsJEVFRaoPmTVg79inTJlCli1bRgghpKGhgYwaNYqUlJTIIapoPPvss+Suu+4ily5dIpcuXSJ//fUXmTdvnkecd77H7m7n/fDhwyQ9PZ3U1NQQQgg5ffo0adeuHXnttdfc/rwLOXZ3O++WmE4J7dmzhxQVFZGGhgZCiD5EvlOnTsYQ+dTUVLcJkSfE9rEXFhaSsrIyQgghNTU1JCcnh5SWlsolqihERUUZ00IcOnSIBAUFGY9RjHveLZSgbdu2kYKCAhIcHExSUlLItddeSwgh5OjRoyQiIoIcO3bM2HbFihUkMzOTZGVlkXvvvVe1c+QG+B77hg0byODBg0lhYSHJuxWFBQAACvFJREFUysoi06ZNIxcuXJBTdKfYv38/AWD1N2/ePLc/70KO3d3Oe21tLXnkkUdIXl4e6du3L+nWrRuZN28eaWxsdPvzLuTY3e28G/jtt99IQUEB8fX1Jenp6eTxxx/nzB+XmZlJMjMzyfPPPy+jxOLB59i//vprUlRUZDzvDz74IKmvr5dZcudYsGAB6dmzJ+nbty/p0aMHef/9943rxLjnGUJUnkKVQqFQKBQKxQE8yieIQqFQKBQKxQBVgigUCoVCoXgkVAmiUCgUCoXikVAliEKhUCgUikdClSAKhUKhUCgeCVWCKBQKhUKheCRUCaJQKBQKheKRUCWIQqFQKBSKR0KVIAqFQqFQKB4JVYIoFBUzbNgwaLVapKam4p133jEuHzNmDCZPnmz8/tVXXyEtLQ1paWlYv369ZPI89dRTiIuLwy233CLZPvgwbtw4BAcHY8WKFby36d27NxISElBYWIg9e/Y4tN/Lly8jJiYGJ0+edGh7sXj22WdRWFgIX19fQb8BheJpUCWIQlEx33zzDbp3747hw4fj1ltvBQA0NTXhp59+wtq1a43tRowYgdGjR2PlypUoKCiQTJ65c+di0qRJkvXPxpo1a9ChQwezZR9++CG6du0quK+pU6eirKwMaWlpDsni7++P3bt3Izo62qHtxeKhhx5CWVkZWrduLascFIrSoUoQhaJyBgwYgDVr1hi/79ixA4WFhTh+/Dj++usv4/I//vgD6enpcojoUYSHh8stAoVC4QlVgigUlTNgwADs2rULZ86cAaC3jEyfPh0JCQn48ccfAQAXLlxASEgIGIYBAGzevBmDBg3CgAEDUFBQgIcffhiEEFRXVyMzMxMMw6CgoAC1tbXYvXs3unXrhsTERGzatAlNTU14+umnkZubi8LCQowcORIVFRWc8lVVVaGkpAQ9e/ZE7969MXPmTNTX1wMAZs6ciaioKNx7772YPHky8vPzkZ+fj6NHjxq3r6iowKBBg9CjRw8MHjwY8+fPh5+fH4qLi7Fp0ybMnj0bFRUVKCwsRGFhIc6dO2fc9vDhwxgzZgwyMzNxzTXX4MKFC7x/1507dyI3NxdeXl5YsWIFBg0ahOTkZHzwwQdYs2YNRo0aheTkZLzwwgvGbYYNG2Y1Dbdx40YUFhaioKAA+fn5mDVrFi5duoSlS5ciJSUFxcXFeOSRR1BcXAwfHx9s2LAB9fX1uP/++9GzZ0/06dMHY8eOxalTpwDAbLu5c+eisLAQnTp1MlOEKRQKT8QvfE+hUFxJXV0d8fPzI6tWrSKEEDJ8+HBSW1tLbr31VjJkyBBCCCFffPEFWbx4sXGb1atXk99//934ffz48WTlypWEEEJqampIUFAQ+fHHH43rH3zwQbJhwwZCCCELFiwgqamp5OLFi4QQQl566SWSm5trbDt37lwydepU4/ehQ4eSG2+8kTQ1NZHGxkZy7bXXkgceeMC4fuLEiaRz587k/+3dW0hUWxgH8P/kqJPpZFJqaJNQlF1QJ8UmHbygg9IU5oNFJRlZoWJEQRfClyi6EVQkWmmJSvVQRj1EDyKV9+jykL1JpWaFt7C8pWP+z0O4cB9P53g4cvQcv9/TzLrttdbLfOz17T09PT0kybS0NO7fv1/VW61WZmRkkCSHh4e5fv16LlmyRNWXl5drvo+KjIxkfHw8h4aGODIywpiYGF64cOGX+xgZGckTJ05oyhobG6nT6VhSUkKSrK2tpdFo5MWLF0mSzc3NdHFxYUdHh2ac0tJSVe/m5sYnT56QJAcGBhgUFMSGhgaSZEFBAWfPns2XL1+SJC9fvszXr19z3759TExMpMPhIEkePXqUFouFIyMjqp+bmxvfvHlDkiwqKmJwcPC4Nfn5+am5CCHGkztBQvzHubq6wmq1ory8XN1hMRgMiI+PR2VlJRwOByoqKhAfH6/6rF69Grm5uYiIiEBMTAxqampQX18P4Gdey/bt21WitcPhwPPnzxEZGQkAuHLlCnbu3AkPDw8AQFpaGp49e4aGhoZxc2ttbcXDhw+RnZ0NnU6HWbNmITU1FYWFhZp2CQkJcHd3BwCEhYWhsbERANDU1ITq6mqVaO3k5PS3co6SkpLg7OwMnU6H0NBQNe7fQRLJyckAALPZjG/fviEuLg4AYDKZ4Onpiffv3/9h35KSEgQEBKg8LIPBgOLiYixatEi1WbFiBdasWQMAyM7OxsqVK3Ht2jWkp6dDr9cDALKyslBfX49Xr16pfoGBgVi1ahUA7Z4JISZOP9UTEEL8czabDbm5uaiurlbBSlxcHPr7+1FXV4cPHz5ofni3bduGBQsWoLKyEnq9Hjk5OZonmvbs2QOLxYKOjg5UVVVh48aNqq6lpQXFxcV49OiRKlu8eDHa2trGJSOPHmsdOHAALi4uAIDv37/Dw8MDvb29KvCZN2+e6mMwGFQwN3rMNjbBd/78+RPeF09PT8243d3dE+47ysnJSc1zNCgZO66zs7Oa7++1tLTAx8dHUxYSEqL5PnbtANDW1obBwUF4e3urstExWlpaEBoaOq7f2D0TQkycBEFC/A/YbDYcPnwYeXl5OHbsGADAy8sLZrMZN27c0ARAAFBXV4f8/Hz1oz40NKSpN5vNCAoKQlFREerq6nD9+nVVZzKZkJ2djaysLFXW3d2NOXPmjJuXyWQCABQWFmqeuGpvb1eBxZ/x8/NT7UfH6uzs/Mt+04XJZEJNTY2m7O3bt5g7d+4vgzkfHx+4urqqHCAA6vPoHgghJocchwnxPxAcHKzu7JjNZlVus9lQWloKm82mab98+XJUVlYCAAYHB1UC9Vh79+7FpUuXYDQa4eXlpcozMjJQWlqK/v5+AMDnz58RHh6OgYGBcWP4+/vDbrdr3mFUUVGBLVu2TGhdAQEBsFqtKCgoAAD8+PEDt2/f1rQxGo3o7e0FAJw/fx7379+f0Nj/hh07dqCpqUm9m6mvrw/Jyclq7/6Ik5MTdu/ejYKCAgwPDwMA8vLyEB4ero7NhBCTZKqTkoQQk2Pr1q1MSUnRlFVUVFCv1/Pr16+a8hcvXjAkJIRhYWHcvHkzN23aRF9fXx4/fly16enpoYeHB58+farpOzw8zNOnT3Pt2rWMjY1lbGwsq6qqSJInT56kyWSir68vDx06RJL88uULd+3axXXr1jE2NpZJSUn8+PEjyZ9J1N7e3gwICODVq1d57949Llu2jJ6enszMzCRJtra2MiEhgWazmXa7nefOnePSpUvVfBwOB202Gy0WC61WKzs7O5menk6j0cjAwEA+ePCAeXl5al5j1zjW7xOjGxsbGR4eTgCMjo5mV1cXo6KiCIAWi4XNzc202+10cXFhSEgIa2trabfb1XVv3rxJkqypqWFUVBSjoqIYERHBsrIykuSdO3fUWqOjo9ne3q6uPTAwwIMHD6o1paSk8NOnT+P6ZWZm8t27d5p5Dg4OqnEkMVqIP6cjySmOw4QQ4pe6uro07965desW8vPzUVVVNanXsVqtSExMRE5OzqSOO5X8/f1x5swZpKamTvVUhJiW5DhMCDGtHTlyBI8fPwbw80m1kpISbNiwYdKvs3DhQpSVlf2jv82YLkb/NsNgMEwo90qImUruBAkhprW7d+/i7NmzcHd3R19fH6Kjo3Hq1Ck4OztP9dSEEP9xEgQJIYQQYkaS4zAhhBBCzEgSBAkhhBBiRpIgSAghhBAzkgRBQgghhJiRJAgSQgghxIwkQZAQQgghZqTfACWJXwi/kh7TAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "Image('model_spectrum.png')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Calculating synthetic photometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Adding filter: Paranal/NACO.Mp... [DONE]\n"
     ]
    }
   ],
   "source": [
    "read_model = species.ReadModel(model='drift-phoenix',\n",
    "                               filter_name='Paranal/NACO.Mp')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Flux density [W m-2 micron-1] = 3.43e-16\n"
     ]
    }
   ],
   "source": [
    "flux = read_model.get_flux(model_param)\n",
    "print(f'Flux density [W m-2 micron-1] = {flux:.2e}')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Downloading Vega spectrum (270 kB)... [DONE]\n",
      "Adding Vega spectrum... [DONE]\n",
      "Apparent magnitude [mag] = 12.00\n",
      "Absolute magnitude [mag] = 10.50\n"
     ]
    }
   ],
   "source": [
    "app_mag, abs_mag = read_model.get_magnitude(model_param)\n",
    "print(f'Apparent magnitude [mag] = {app_mag:.2f}')\n",
    "print(f'Absolute magnitude [mag] = {abs_mag:.2f}')"
   ]
  }
 ],
 "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.0"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
