{ "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.2.0... [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" ] }, { "data": { "text/plain": [ "" ] }, "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": [ "Unpacking DRIFT-PHOENIX model spectra... [DONE]\n", "Adding DRIFT-PHOENIX model spectra... [DONE] \n" ] } ], "source": [ "database.add_model(model='drift-phoenix',\n", " wavel_range=None,\n", " spec_res=None,\n", " teff_range=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": [ { "data": { "text/plain": [ "{'teff': (1000.0, 3000.0), 'logg': (3.0, 5.5), 'feh': (-0.6, 0.3)}" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "read_model.get_bounds()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [], "source": [ "model_param = {'teff':1510., 'logg':4.1, 'feh':0.1}" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "modelbox = read_model.get_model(model_param=model_param,\n", " spec_res=200.)" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [], "source": [ "filters = ['MKO/NSFCam.J', 'MKO/NSFCam.H', 'MKO/NSFCam.K', 'MKO/NSFCam.Lp', 'MKO/NSFCam.Mp']" ] }, { "cell_type": "code", "execution_count": 10, "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.1e5),\n", " legend='upper right',\n", " output='model_spectrum.png')" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkEAAAGiCAYAAADgPBIcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdd3xT5f4H8E/Sle49oAtaZG9oKSBSEBmCCoooIsoQqQMVC+rV66UqKooIeF3cH2q5IKDIRRRkiFBEhoBQsHRQ6KJ7pjNp1vP7ozYSmrYZJzknzff9evWlOX3ynO85DT3fPlPEGGMghBBCCLEzYr4DIIQQQgjhAyVBhBBCCLFLlAQRQgghxC5REkQIIYQQu0RJECGEEELsEiVBhBBCCLFLlAQRQgghxC5REkQIIYQQu2SxJCg/P99SVRNCCCGEmM2Ri0oYYygvL0dzc7P22Ny5c3Hq1CkuqieEEEII4ZzZSdDWrVuxbNky1NfX4+YdOEQikblVE0IIIYRYjNndYW+//Tb279+PpqYmaDQa7dcdd9zBRXyEEEIIIRYhMncD1cmTJ+Pw4cNtjiuVSjg5OZlTNSGEEEKIxZidBP373/9G7969MWXKFJ3js2bNwp49e8wKjgvULUcIIYR0XeakMWYnQT179kRJSQnc3Nzg4+OjDaisrAxNTU3mVM0JSoIIIYSQrsucNMbsgdE+Pj5ITk7WOcYYw/Lly02u89y5c5g7dy7++c9/YsGCBe2W+/rrr/Hhhx9CJBJh/Pjx+OCDD9pNeszM9WyWSCSia7dDdO107faGrt3+rp2LRg6zk6C1a9di/PjxbY5/9tlnJtW3Z88e7Nq1C97e3h2WS0tLQ2JiItLS0uDn54c777wTn376KZ555hmTzksIIYQQ+2L27LBJkyZBLpdj165dWLduHb777jvIZDLExcWZVF9MTAy2b98OT0/PDst98cUXuPvuuxEQEACxWIxFixbh888/N+mchBBCCLE/ZrcEZWdnY9KkSaisrERAQID2vz///DN69+5tdH1hYWEGlTt37hzuvfde7ev+/fvjypUrkMlkcHV1Nfq8hBBCCLEvZrcELV++HKtWrUJ9fT3y8/NRX1+Pf/3rX3jhhRe4iK9dZWVlOl1mPj4+YIyhsrLSoue1NatWreI7BN7Qtdsnunb7RNdOTGF2EtTQ0IBFixZBLG6pSiwWY/HixVaZGaZvUFR7g8NEIpHOV1JSkoWjEwZ7uU596NrtE127faJr7/qSkpJ0nuNcMDsJUigUkMlkOscaGhogl8vNrbpDQUFBkEql2tdSqRQikQiBgYF6yzPGdL7s5UNDCCGEdAVJSUk6z3EumJ0E3X///YiLi8P777+P5ORkvPfeexg9ejQeeOABLuJrV0xMDLKysrSv09PTMWDAABoPRHh36z/OhoYGrF27Fn379oWbmxsuXryo8/3m5maIRCLs3bvXmmESYlcYY2hsbLTLqeSkfWYPjF6xYgWcnZ3x+eefIz8/HxEREUhISMBzzz3HRXxalZWVmDVrFvbu3Qs/Pz888cQTuOuuu1BVVQVfX18kJycjISGB03MS+1NVVQUA8Pf3N/q9v/32G8aNG9fu9+fMmYOAgAAMHz4cQMvmw71798aoUaMwZswY7N27F/fdd59pgRPCAbVaDZVKhYyMDOTl5aGwsBB1dXUYPHgwZsyYYbU4lEolZDIZlixZgrCwMKhUKjz00EOIi4vTDr0wxokTJ7T7WS5dupRmEpO/MQuprKw06X3nz59n48ePZ97e3qxPnz5s1qxZjDHGCgoKWEBAALtx44a27LZt29jw4cPZyJEj2Ysvvsg0Gk2b+gAwC16mwa5fv84AsG+//ZbvUEg71q5dq/28rFmzpt1y33//PZNIJMzX15cdOHCAFRcXa9+3e/duJpfLWXl5OZs3bx775z//ybKysnTer9Fo2Jtvvsn69evHALAlS5YwjUbDpkyZYulLJF2UVCplmzZtYiqVyuj35uXlsb1797KxY8dqP8cA2NKlS9moUaN0jg0cOJDNnTuXbd68ucM6m5qaWHNzszaesrIylpKSwpYvX858fHzY7t272W+//cZOnjzJNm7cyFauXMnc3d1ZfHw8A8ACAgKYv78/u//++9lLL72kE8eYMWPY//73P6ZUKju9toaGBvb000+zmJgYdvHiRaZSqVhgYCBLTU01+j51VVKplH3yySfsqaeeYmvWrGHXr1/nOySDcfF8t1h2MGHCBEtVbRQhJEEajUb7jxeAQf94iXXl5+czX19fVldXx9RqNfPw8GBvv/22ThmNRsNefPFF7WfK3d1d+/8vvPCC2THcfffdehN5Qm5VWVnJ9u/fz2bMmKGTpERERDCZTGZQHQUFBWzKlCkMAOvTpw+bNm0a27p1K6uoqGhTVqPRsObmZrZt2za2YsUKNnjwYBYeHs4OHDjANBoN02g0rLCwkDU3N7N33nmHBQYGauttjW3cuHEsLi6OSaVSNnbsWBYUFMRmz57NXn31VbZkyRI2c+ZMtm3bNpabm8sUCoXemPPy8thbb73FIiIimJeXF2toaNBb7vDhw+yhhx5iANi0adNYQUGB9nu5ubnsoYceMugedTUajYZ9++23bPfu3ezYsWMsMTGRTZo0iX377bfs2rVrLDk5mcXHx7NXXnlF7+fgVlVVVbw+z3hLgmbOnMnq6uoYY4yJRCImFot1vlqPCYEQkqB58+ZpY5g9ezabNGmSzveLior4CIsxxtj27dtZSkqKzrEvv/ySffLJJ0wul/MUlXXJ5XI2Y8YM9ueff2qP1dXVsZ49e7L333+faTQaVlNTw6ZOncpCQ0NZdna2ReKYN28eq66utkjdpOsoKSlhANiUKVPYzJkz2eXLl7WJyFtvvdVpi3NGRgYbPHgwCwkJYRs3btT53BtKpVKxxMRE7e9XBwcH5ufnp31dUFDALly4wMrLy1l5ebn2ecEVuVzOdu/ezUJDQxkAdvr0aSaXy9kff/zBNm7cyHr06ME+/fRTtnXrVr3vnz17tsm9Fbbs0KFDbOHCheyLL75g7777Ljty5EibP7w0Gg3bs2cPGzduHPvqq6/afD8zM5MtWLCAjRo1is2ZM4dNmDCBbd++3ZqXocVbErRz507tjYmLi2N5eXk6X7m5uSwuLs6swLhizSRIrVbrPd6tWzdWXl7OGGt5uLb+kmCspZnYGvEpFAo2f/58tmHDBu2xLVu2sAULFrBFixax1atXM7Vazd544w32z3/+k23evJktXbpUb10qlYrV1tZaPGZLqqqqYk8//TRbv349A8A2bdrUpoxarWaOjo7az1BsbCxrbGy0WEyJiYksPT3dYvWTrmHx4sXsyJEj7X7/6tWr2sQkPz9fe7y0tJRt3LiRDRgwgH377bec/BvWaDSstraWFRQUaLu+rNkyoFQq2fHjx9k999zDYmJi2O23385Wr17daYKzfv16duDAAStFKRzz589nJSUlBpVVqVTs7bffZgMHDmR33XUXu/vuu9m0adPYjBkz2K+//qrNAWQyWaefSUsRRHfY6dOnjTpubdZKgrZu3coAaJtxi4uL2f79+9n69evZwIEDdcp++OGHDACrqKhgP/30EwNg8W6Qt956i/33v/9lr7zyClu3bh379ddf2axZs5hKpWIajYZ98sknbOLEieyjjz7SxjJ37lx29epVnXpKSkpYfHw8mz59Onv++efZ8ePH2f/93/+xd999lzU1NVn0Grii0WjYtGnTWGJiInvsscfYe++91+7912g0rLq6mtXX11s8rg8++IAdPXrU4uchtksul7OpU6d2Wk6tVrMVK1YwNzc3ne6yzz77zKRxQ13NyZMn2RtvvMF3GFZnyGfnVvX19ay4uJhVVla2+9mRSqVs4sSJ7TYEWIogkqDS0lJ2/PhxplarmVwuZ2+99RZbtWpVu3211matJAgAe/7557V95Df/4mltBbrZqVOn2IwZM9gbb7zBIiIiOG8uZqzlw7to0SK2ZcsWds8992ibzN955x22dOlSVlNT0+H7MzIy2KOPPqp9rVKp2LRp07StFUePHmWrVq1i27dvZ5s2bWIPPPCA1f8RmCIlJYW9/PLLfIfRRnJyMvvmm2/4DoMIWEpKClu9erXB5TUaDTt79iy7evWq3t9D9qqpqYnNnDmT7zCsqq6ujj344IMWq//tt9+2eusaF893TqbI+/n5YcyYMVi9ejV27tyJ0NBQXL9+HVu3bjW3esHSaDTaqZpyuRyBgYFYvnw53nzzTTQ2NmLYsGG4cOFCu+8fPXo0srOzsW/fPjzyyCNobm7udNNYY/3444+Ijo5GYWEhNm3apF1h8x//+IdB7+/bty9cXFyQmpqKoUOH4uOPP8aMGTPQr18/AMCECRMwYcIEbfmKigocPHgQd999N6fXwbVdu3bh2Wef5TuMNnx8fFBeXs53GETALly4gJiYGIPLi0Qio8rbC1dXV8jlcjDGOFt5WOgyMzPRv39/i9W/aNEiJCYmYurUqRY7hyWYvVhicXExNm7cCEdHRyQnJ2P//v1ISUlBZmYmF/EJ1tixY3H48GEcPHgQ8+bNw4oVKxAZGYmdO3di9uzZOHjwYKd1/Pjjj3j00Ufh7e1tkRW2U1NTce+99+LVV19Ft27dTKrjzTffxIsvvoidO3fi6NGjeOqpp9otm5CQgE8++US7GJlarcbFixcFtzhZdnY2+vbty3cYbfj4+Oisgk7IrdLS0jBw4EC+w+gSwsLCUFRUxHcYVpORkaH9A9YSQkJCUFVVBbVabbFzWILZSZBKpQLQshhVeHi4dud4X19fc6sWtMuXL+P8+fOYNm0aVCoVnnnmGQDA559/jilTpiAoKKjTOm677TZs3boVEomkzdYjXMjOzkavXr3MqqN79+5Yv349CgsLsWPHjg7/avL398ewYcNw8OBBMMbwwgsvYMWKFdi9e7dZMXCpqKgIISEhfIehFyVBpDMlJSUm/0FDdA0YMADp6el8h2E1+fn56NGjh0XPMWLECPzxxx8WPQfXzE6CfH198dRTT+Hpp5/G4sWLAQB//PGHNjnqqmJiYvDaa69hzZo12Lt3L9zd3QEA8+fPN6gV6GYSicQiLUEKhQISicTseoYMGYIVK1bAzc2t07IrV67E+++/j2effRZhYWH43//+h507d5odA1d+//13jB07lu8w9KIkiHREo9FwunGkvRswYACuXLnCdxhWU1hYiNDQUIueY+rUqUY///hmdhL05ZdfwsfHB/Pnz8fixYtRWlqKjz/+GC+99BIX8QmWh4cHpk+fjkWLFpldl6WSID54e3tj27ZtmDlzJl566SV4e3tbpJXLVFeuXMGgQYP4DkMvSoJIRwoLCxEREcF3GF1G//79kZGRwXcYVlNaWmrxVvBRo0bh3LlzFj0H18weGO3n54d3331X+zokJARfffWVudXahH379nFSj6WSIL7+YgwNDdX5i8PFxQVyuZyTVilzpaWl4fnnn+c7DL08PT1RV1fHdxhEoIqKihAWFsZ3GF1G9+7dUVxczHcYVqNSqeDoaPYjv0POzs7QaDQ6E4eEzuwoDxw4gEWLFkEulyMjIwPR0dHo0aMHTp8+zUV8gsVlgmGJJEitVgvmQxgVFYXc3Fy+wwAA1NfXw8vLi+8w9BKLxYIbRE6Eg8YDcYu6FS1DSL/vDWH2U3L9+vVYtmwZJBIJXn/9ddx1111YtWoVEhMTuYjPLri6unLeZSSVSuHj48NpnaaKjo7GtWvX+A4DCoUCzs7OfIdBiEkoCeKeWCy2udlMplCr1XBwcLDKuWxtrJXZSRBjDMOGDUN9fT1SUlKwfv16LFy4EE5OTlzEZxcs0RJUU1MjmBl6vXr1EkQSdPXqVe3sRaGiv05JeygJ4l63bt1QUlLCdxgWV1NTAz8/P6uca+DAgUhLS7PKubhgdhJUW1sLhUKBr7/+GlOmTIGrqysAWC3r7ArsIQm6fv0632EgIyPDoouFEWJJlARxLyIiAjdu3OA7DIurrKyEv7+/Vc5lay1BZo+SevDBBxEcHAyFQoETJ06grq4Or732msWn4vGJ6/E2EokE1dXVnNUHANXV1YJJgsLDw1FQUMB3GCgqKqKF5ojNqqioQEBAAN9hdCkREREoKCjA6NGj+Q7Foqqqqqz22fH19UVNTY1VzsUFs5OglStX4t5774W7uzvCwsIgk8kwe/Zsi65Mybfm5ma4uLhwVp+LiwsUCgVn9QHCaglydHSERqPhOwwUFxdj8uTJfIdBiEk0Gg21sHMsIiICZ8+e5TsMi7NmSxDQ0q1vK1uScNKc0adPH+3UTVdXV4wfPx5r167lompB4nq6t7OzM5qbmzmrDxBWEiQUxcXFgu9OEIlEgkgYCbEHQmmltjRrtgQBQFBQECoqKqx2PnOY1BL0v//9D7NmzYJIJGp3scCDBw922USI65YgZ2dni7QEDRkyhNM6bZ2QZsy1x83NDU1NTfDw8OA7FCIgKpWKWoEsIDQ0FIWFhXyHYXGVlZWIjo622vlaJ8MYsn0U30xqCfryyy/R0NAAADh27BgiIyPbfAlhYTxLoe4w0/C9Bo4tNM+6u7ujsbGR7zCIwFRUVNjEA8XWuLi4QKlU8h2GxVm7JSg6OloQk2EMYVJL0M0rJT/++ONYtWpVmzJ8P/AsyVZagoSUBLU+3Pls4RB6AgS03KeGhgYEBwfzHQoREJoZRsxh7TFB0dHRnO2oYGlmjwlKSkrSe7wr/9XC9aJ79pAE8T1jwFaScg8PD2oJIm1QEmQ5bm5uXf7fXFVVlVWTIKEsi2IIs2eHaTQa7N27F1evXtUZ3JucnIynn37a3OoFyRJJENcDo7lurTKXj48PamtrER4ezsv5hbxdxs0oCSL6UBJkOeHh4bhx4wb69u3LdygWo1KprLqAMd9/9BrD7CRo7ty5SE9Px5AhQ3RusqmL/8nlciQkJCAzMxMqlQrvvPOO3mnNycnJ2LBhg85A19WrV+P222836bzGUCqVgm8JEhovLy9eNwetqqqy2oqp5vDw8EB9fT3fYRCBKS0tpTWuLMQekiDSPrOToKysLFy6dKnN4oHtdZN1JikpCYwxnDlzBlevXkVcXBwyMjL0jpHYsGED4uPjTTqPObhuCbLEwGih8fb25jUJqq6utmpzsKk8PDy0kw4IaVVWVtalhxjwqTUJItzy8vJCbW0tvL29+Q6lQ2aPCRowYIDe48OGDTO6Lo1Gg82bN2Px4sUAgN69e2PYsGHYtm2bWTFyzRbGBAlN6z8IvlRXV9tES5Cnpye1BJE2ysrKaLC8hXT1JIiv8ZC2Mi7I7Jag9957D0888QSGDh2q0zW1Zs0a3HfffUbVlZOTg6qqKp1myf79++P8+fN6y2/atAmrVq2CSCTCnDlzrDYGSehJkBDXFPHy8kJ5eTlv57el7jA+7xMRpqamJri7u/MdRpfU1ZOguro6XlpjWqfJDx8+3OrnNobZSdDatWuxb98+ZGdnw9Hx7+pKS0uNrqusrAwAdH5gPj4+SE9Pb1M2ODgYkydPxoIFC1BRUYH4+Hio1WosW7bMhKswDtdJkJOTE6cDo6VSqaBmhgEtSRCfO8lLpVLtquZCRt1hhFhXcHCw9tnTFfH1B2B0dDROnjxp9fMay+zusJ9//hkFBQU4ceIEjh07pv2aM2eOyXXeup6Lvua8adOmYeHChRCJRAgKCsKyZcvw6aefdljnzV+mjlkCuE+CHBwcON0qQWjT4wH+B0bbwmrRAHWHEWJtYrHYZpbQMAVf4yFbV43mUlJSks5znAtmJ0F9+/bVuzr0woULja6rdeCfVCrVHpNKpQYNCIyIiEB+fn6732eM6XwJKQkCuF3IT4hJkLe3N69jgmwlCaKWIHIrhUJh1enNpGvhqyUoJCTEpB6hjrROnGr94oLZSdD48eOxaNEi/PDDD/j111+1X0899ZTRdUVHR8PPzw9ZWVnaY+np6YiJiWlT9v3339d5XVZWhu7duxt/ASawRBLEperqasElQXy3BAkxMdTH09OTkiCio7y8nGaGWVhXXpqCr5ag1p3khc7sMUGvvPIKQkJCcOzYMZ3jpvSxisViLFmyBF9++SXGjRuH7OxspKam4uuvv0ZGRgaeffZZHD58GA4ODjhw4ABGjRqF8ePHo6mpCZs2bcL8+fPNvRyDCD0JEuIDn+8kyJZagrrqL2NiGpoeb3lhYWEoLCxEv379+A6Fc1VVVVbdPPVmLi4uglu491ZmJ0Fjx47FkSNH2hyfOXOmSfUlJSUhISEBcXFxUKlU2LFjB0JCQpCXl4fMzEwolUo4ODggMTERb775JjQaDRoaGnDnnXfiH//4h7mXYxBbSIIiIyP5DkOHi4uLyQtocqGxsRFubm68nd9QrbvIE9KqvLycpsdbWFdOgqqrqxEbG8vLuXv06IG8vDz06dOHl/MbwuwkSF8CBADff/+9SfVJJBIkJye3OR4XF4eioiLt6xkzZmDGjBkmncNcQk+CqqurBTctUQiblwohhs7YQozEuqglyPLCw8NRWFjIdxgWwecaaVFRUcjJyRF0EmT2mKCzZ8/izTffhEKhQFFRESZNmoT4+HidcT1djS0kQbawJg7Rzxb60Yn13Lhxg7c99+xFa0tQV8TnGmmtawUJmdlJ0BtvvIFBgwbByckJr732Gtzd3TF+/Pguu3kqIPzZGpQEEdJ15OXloWfPnnyH0aV15SSIz/GQrS1BQmZ2d5hCocCsWbPQ3NyMffv24dq1a/Dx8bHKRqZ8EXpLUE1NjU0MArYWalkhtoy2zLC8kJAQlJSU8B2GRTDG2uztaS2tY4KEzOw70zqd9/vvv8ftt9+uffi6urqaW7VgCT0J0mg0Oqt3C4VYLIZarbb6eWUymU0MiiZEH8YYjRWzMAcHB15+N3V1rbPDhMzsJ+Xtt9+O/v37o6ioCHv27IFMJsOnn34KDw8PLuITJKVSKegkSKi8vLxQX19v9VYqW5ke30okEkGj0fD21xsRDrVaTZ8DK3F0dIRSqRT0UAdbJeREnpO9wyZNmgRfX1/Exsairq4OAQEBbRYz7EqE3hIk1O6f1rWCKAnqmLu7O5qamrr0HxLEMEVFRQgNDeU7DLsQGRmJgoIC3tbUsQSVSsV7Eh0cHCzoZR44uTtTpkzRrkPg5eWFxx9/HEePHuWiakGyRBLEVVeRUBMggL8FE4W4eGRHaOsM0ooGRVuPJfa64ltFRQXvyysIfYaY2S1BarUae/fuRXZ2tk7fX3JyMpYuXWpu9YJkiSTI2dlZuxCkORoaGuDl5cVRVNziKwmytZYg2jqDtMrLy0OPHj34DsMutCZBU6ZM4TsUzpSUlKBbt268xtA6Q2zMmDG8xtEes5OgRx55BJcvX8bQoUN1NlLlc3VgS7NUEqRQKPRuRmsMIU+P52sTVVubLUdbZ5BWeXl5XeqhLGS9evXCzz//zHcYnCotLRVEEnTw4EFeY+iI2UlQRkYG0tLS2rRgvPXWW+ZWLViWSoK4GEUv5CSIr5agiooKDB482OrnNRV1h5FWubm51BJkJT169EBubi7fYXBKCC1BQu8OM3tM0G233aa3C2fy5MnmVi1YKpXK7G6rW7m4uEChUJhdT0VFBQICAjiIiHt8JUHl5eUIDAy0+nlN5enpSS1BBADtIG9NrUMSupKSkhKEhITwGoO/vz8qKyt5jaEjZrcE3XHHHXj88cdx//336ww+feGFF3DhwgVzqxcsrqf7SSQSyGQys+spLS3l/UPfHi8vL14GHlZUVNhUEuTl5cVLtyERJqFOLe6KxGIxVCqVINdZM4UQusNEIpG2p0OIu8mb/ZN+5ZVXEBISgl9//VXneFlZmblVC5YlZmBx9dd/aWkp4uLiOIiIe3w93CsrKwXbOqZPQEAAsrOz+Q6D8MwSLc6kY3369EFWVhYGDBjAdyicEMofxQMHDkRaWhpGjBjBdyhtmN0dNnbsWOTm5rb56srdYZbAVVeRUD70+nh7e/PSHaZSqWxqAbTAwEBUVFTwHQbhWVFREcLCwvgOw64MHToUqampfIfBGZlMZvZkGy4MHz5csD1DZidBR44c0Xt848aN5lZtV7hsCRJqEsTXmCBbQ0kQAVoGRdMaQdYl5Ie1sTQajWC6UoV8Xznr+KyqqkJjY6P29dy5c3Hq1CmuqhcUS3ywuEqC6urq4OnpyUFE3OMjCRLy4pHtoSSIAMCpU6cE2X3QlbV2h3UFQmpJjIyMFOxGqmYnQRcvXsScOXOQk5Oj88ARSgbKNcaYRR6sXl5enO1iLNR77+Liwsngb2NUVVXZ1HggoGXzYWvfJyI8586dwwsvvMB3GHbFwcEBGo1G0HtdGSorKwt9+vThOwwALc8kNzc31NbWwtvbm+9wdJjdHfbKK6/g8OHDuOOOO6DRaCCXy7Fr1y68/vrrXMQnOJbaN8zT05OTVhIht3zw8UvFVtdZEfLPkVieSqWCXC6Hm5sb36HYnX79+uHSpUt8h2G29PR09O3bl+8wtO666y5BLkZpdhLEGEPPnj21v7SdnZ3xwAMP4OLFi2YHJ0SNjY0W2diSi+4wqVQquCz7VmKxGBqNxmrns9VtB5ycnDhZN4rYpt9//x2jRo3iOwy79Oijj2LLli18h2G2X3/9VVBbVUyfPh0//fQT32G0YXYSpFKpWioSi/HHH38AAK5evYrLly+bW7UgVVZWwt/fn/N6fXx8UFNTY1YdtjCQMiQkxKrLJ9hqEtSjRw/k5+fzHQbhyaFDh2i7DJ4MHToUaWlpNt0aq1Qq0dTUJKg/isPDw1FUVGTVP4INYXYSFBMTg3Xr1mHp0qWIi4tDeHg4Bg4ciAcffJCL+ATHUrvyhoaGoqioyKw6bCEJioiIQEFBgdXOd+XKFcH0ixujV69euHr1Kt9hEB6o1Wr89ttviI2N5TsUuyQSiTBq1CicOHGC71BMdvr0aYwePZrvMNqIi4vD8ePH+Q5Dh9kDo9977z3t/4eFheH06dPo27cvZsyYYW7VglRWVmaRJEgikZi96eyff/4p+L8eIyMjkZOTY7WmfiHsnWOK2NhY7N+/H9OnT+c7FGJlBw4cwNSpU2mhRB4tWrQIzw98EDgAACAASURBVD33HMaNG2eTA6QPHTqEmTNn8h1GG0uWLMGKFSswYcIEvkPRMrslKCoqCk8//TSAloUTV6xYYXYCJJfLsWDBAsTFxWHkyJE4fPhwu2XXrVuHESNGYMSIEVi7dq1Z5zXE9evXERUVZZG6zW1+vXz5cpuNQpOSksyqk2sjR47EuXPnrHKu5cuXIzw83Crn4tqwYcPMWrRNaD93a7Llaz9//jzuueceLFq0yKT32/K1m4vLa4+KisLQoUOxf/9+zuq0pFuv/cKFCxg+fDg/wXQgLCwMYrHYqr0BnWJmGjJkiLlVtPHyyy+zxx57jDHGWFZWFvP19WWlpaVtyh04cID16dOHyWQyJpPJWJ8+fdi+fft0ygBgHFym1uLFi1lxcTFn9d1s4cKFrLCw0KT3qlQqNmXKlDbHubx2LqjVar1xWgIA9s0331jlXJYwffp01tDQYNJ7hfZztyZbvHa1Ws2++uorBoB99tlnJtdji9fOFa6vvb6+nvXv35/99ttvnNZrCTdfe05ODps/fz6P0XQsJSWFPfzww0yj0ZhdFxfPd7Nbgvr06aN3590XX3zRpPo0Gg02b96MxYsXAwB69+6NYcOGYdu2bW3Kbtq0CY888ggkEgkkEgnmzZuHzz//3KTzGkKtViM3N9diKzLfcccd+OWXX0x67y+//CKomQDtEYvFCAkJsfhGqlKpFAAE2SRsqMcffxwbNmzgOwxiQSUlJdiyZQvuueceLFy4EOfOnUNCQgLfYREAHh4e+PXXX/H555/jySefRFVVFd8hderPP/9EfHw8/vGPf/AdSrvGjx+P/v37Y/ny5SguLuY7HIj+yqZMtnr1avz8889tdpFfs2YN0tPTja7v2rVruO2223TG3ixbtgyVlZXYsWOHTtmwsDB8+OGHmDNnDgBg9+7dePbZZ3UWHWztzzXzMgEA+/fvR2pqKl577TWz69KnoaEB06dPx5YtW4ya0XTjxg089thj2LVrV5uFAUUikeBmOaSnp2PlypXYvn0757MX5HI59u/fj48//hgpKSmCu3ZjMMawcOFCREVF4eWXXzZqB2Yh/tytRWjXzhhDVVUV8vLyUFRUhAsXLiA1NRU//PADxo0bh/DwcMyaNQuzZ882+1xCu3ZrsuS1P/HEE/jzzz/Rs2dPDBs2DI2NjZgwYQJGjhzJ6wr9tbW1cHV1hYuLCxISEnD58mXs2LEDERERvMVkCMYYfvnlF2zevBnOzs5YvHgxbrvtNoSEhEAsNrxthovnu9lJkK+vL4YOHdrm+KVLl1BdXW10fSdPnsTtt98OuVyu/aX/+uuv49SpU21aSZycnLBv3z7tYOBffvkFkydPhlqt1paxxUFthBBCCDGMOWmMSbPD3nzzTQAt+4NNnz5db1eVuU26tyYv7V0kJTmEEEIIMYVJSdCePXuwceNG+Pv7602AAJg8Nqe1C0wqler9/1vLto79aC0XGBiot15qIrY/dO107faGrp2u3Z5w0Qhi0sBoHx8f3HHHHfDz82u3THZ2tkkBRUdHw8/PT2cn3/T0dMTExLQpGxMTY1A5QgghhJBbmZQEGZJ9LV261JSqIRaLsWTJEnz55ZcAWpKp1NRUzJs3DxkZGbjzzju1Y34SEhKwY8cOyOVyyOVybN++nWZWEEIIIcQgJnWHpaSkWHQ106SkJCQkJCAuLg4qlQo7duxASEgI8vLykJmZCaVSCQcHB0ydOhVXrlzB2LFjAbSs8kkr7BJCCCHEECYlQUOGDOlw/RLGGJYvX25yUBKJBMnJyW2Ox8XFtdlfKzExEYmJiSafq6tbtWoV3yHwhq7dPtG12ye6dmIKk6bIT5w4EUePHu2wzJYtW/D444+bHBhXuFwniBBCCCHCwMXz3ewVo9sjhASIEEIIIaQ9JiVBJSUlmDhxIn7//Xeu4yGEEEIIsQqzV4wWOuoOI4QQQroeQXeHEUKIEG3atAnx8fEYOHAgwsLCEB8fjzvvvBMqlcoi5zt9+jSGDBmC7777Tud4XFwc4uPjER8fj7Vr13ZYdvv27XjppZfw1FNP4ezZs3rP01GZU6dOYejQoZgzZw4qKyuxfft2xMTE4Icffmg37ry8PEyaNMmUS7aI33//HU5OTqisrGzzPXPuW6vk5GTEx8dzEmtKSgqGDx+O//znPzr3PikpCUlJSZg4cWKnddxzzz0675VKpSgtLcXMmTMRGxuL8+fPIzMzEytWrDArbkPukbn3l8t7yzWTZocRQoitWrp0KZYuXYqNGzdCoVBg5cqVFj1ffn4+hg0b1ub41KlTkZSU1GnZ2tpafPjhhzh//jzq6+sRHx+PP/74w6gyY8aMwdChQzFjxgyoVCqcOnUKR48e5XXzT2PIZDJs27YNoaGher9v6n272YIFC/TOSjbV4MGD8eSTTwKA9t63bpL7008/dfje69evIzo6Wufn5uPjAwCYOXMmzp8/j5EjRwIAnn32WSxYsMCkGA29R+beX67vLZcoCSKE2KXLly9jzpw5Jr1XLpdj5syZbY5Pnz4dy5Yt0zn28MMP4+DBg3rP//7776O+vh4LFy5EVFSU3rK///47+vXrBwDw9PSESqVCSUkJunXrZlQZACgoKMCaNWuwdu1auLq6GnStjDEkJCTgypUr0Gg02LBhA2JjY6FWqzF//nwUFRVh9OjR+O6773Dt2jXt+7799lusX78eHh4eCAwMxKJFi7QtS8bcPwBYs2YNVq5ciR9//FFvjKbeN41Gg8WLF8Pf3x9RUVE617xixQr4+fmhsLAQCQkJ2LVrFzZv3oyzZ89i3rx5mD9/vjbJMcaaNWuQkJCg9xxDhgwBAHz//fe47777jK67laH319DPjSn3l497awpKggghduny5ctYvXq1Se+VSCR6ExtjvPrqqxg5ciRKSkowdepUXLp0SW+5yspKeHh4aF97enqisrJS50FlSBkAeOONN/Duu+8anAABwN69e6FUKvHbb78hJycHDz/8MM6ePYu9e/fCy8sL27dvx8mTJ7Fz507te65fv47PPvsMx48fBwB0794dH330kfb7xty/48ePo0ePHoiIiDA4ZsCwe7J37144OTnhgw8+QE5ODr799lsAwP79+yGVSrFu3ToUFBRg4cKF+OWXX5Cfn4+vv/4ac+bMMfohnZycjN9++w3nz59HQkJCu+cAgBMnTuD5559v816gZXuo3r17d3guQ++voZ8bU95rzXtrDkqCCCGCsOv8DRTWyMyuJ8zXFQ+ODO+wjEajQXl5ebu/7LOysvDDDz8gNjYWZ8+eRVxcHMaNG6f9vlwux4wZM9q875577tF5eHWktTujW7duUCgUqKio0LsBdEBAABoaGrSv6+vrERAQYHQZAHj//ffxySefIDAwEA8++KBBcWZlZWHMmDEAgKioKNTU1ABo2dKoda/GUaNG6Wyn9PXXX2PJkiVwdnYGAISGhurEY8z9+/nnn+Hh4YE1a9agtrYWH330EebNm4c+ffp0GLch9+Tq1avo2bMnAKBHjx7a4+np6SguLsaaNWvAGIO/vz8AYN26dejRowdyc3M7PLc+CxYswOzZs5GbmwsPD492z1FZWQk/Pz84Ojq2eS/QkhCdP3++w3MZen8N/dzo09l7rXlvzUFJECFEEDpLXLiUnZ2NXr16aV9/9913aG5uRlZWFlauXIljx45BpVIhLS0NKpUKw4cP13m/RCLBkSNHTD5/ZmYmzpw5gwULFqC5uRlNTU3tbkg9atQovPrqqwBaHjQODg7a5K2goAAREREdlrmZv78/Dh48iIkTJ8LLywtTpkzpNNY+ffrghx9+wBNPPIGcnBzt2JRevXrhyJEjWLx4Mc6dO6czQ0cmk2mToq1btyI8XPdna8z9u7m17vPPP8dzzz2nfdi2Xr8+htyT3r1748CBAwBaxr20GjBgAG7cuIFXXnkFjDHs3r0bQEtyt379ejz33HP45ptvDIr/Vq2JQXvn+PHHH/UmMMYw9P4a8tky5b0AP/fWFJQEEULszuXLlzF48GDt6++++w4LFixASEgICgoK0LdvX53y7u7uJp/rv//9Ly5fvozm5mZ0794dY8aMgZeXF3744QcUFRUhPz8fH330ERwcHPSW9fb2xosvvojExEQ0NDTgs88+A9AyMHXKlCnIyMhot0yrM2fOaOudMGECfvrpJ0yZMgULFizAa6+9hp49e7b7F/i9996L/fv34/bbb4darca///1vAC0DdHft2oXx48cjJiYGLi4u2vc89thjePDBB7Fjxw54eHhoW4zM8cEHH2gH47788ssAoL1+Y+7brdf2/fff48UXX0RQUBDy8/Nx7NgxTJ8+HSdOnMAbb7wBmUyG+Ph4fPbZZ/jpp5+wbds2vPPOO3jhhRewYcOGDu8doHvv+/bti4EDBwKA3nMAwKFDh/DFF1+0ee9dd90FuVyOH3/8EUVFRbhw4UKb5NxYhny2AP2f4c7uLxf31hponSBCiN3btWsXRCIR6uvrMXPmTHzzzTeoq6uDl5cX6urq8NJLL/EdosXs2bMHZ8+exbvvvmv0e5VKJZycnHDy5Em8++672Ldvn85xpVKJ++67D2vXrsWAAQO4Dp137d27lJQUJCcnmzQjavv27XjkkUcMLp+Xl4cFCxYgJSXF6HPZOi6e75QEEUKIHbtx40ab7ipDPfDAA6isrERzczM2bdqknd302muv4cSJE1CpVHjyySdNnsItdO3du0uXLmHdunUYP348Fi9ebLHzZ2VlYceOHaivr8e6dessdh6hoiTIAJQEEUIIIV0PrRhNCCGEEGIiQSZB586dQ69evTrtT/36668xYsQIjBw5EomJidTaQwghhBCDCS4J2rNnD9avXw9vb+8Oy6WlpSExMRGHDh3C2bNnceHCBXz66adWipIQQgghtk5wSVBMTAy2b9/e6Z42X3zxBe6++24EBARALBZj0aJF+Pzzz60UJSGEEEJsneCSoLCwMIPKnTt3Tmctj/79++PKlSuQycxfcZYQQgghXZ/gkiBDlZWV6XSZ+fj4gDGGyspKHqMihBBCiK2w2SQIgM5eNa3aGxwtEol0vpKSkiwcHSGEEEK4kpSUpPMc54LNbpsRFBQEqVSqfS2VSiESifRuQAjQOkGEEEKILUtKStJpwOAiEbLZlqCYmBhkZWVpX6enp2PAgAFwdXXlMSpCCCGE2AqbSYIqKysxbtw4VFdXAwCeeOIJ/PTTT6iqqoJGo0FycjISEhJ4jpIQQgghtkJwSdAff/yB+Ph4pKamYs2aNbj//vsBADKZDJmZmWhqagIADBw4EB988AEmT56MUaNGYejQoXj66af5DJ0QQgghNoT2DiOEEEKIzaG9wwghhBBCTERJECFmUKg0OJdXrXOsSCpDvVzJU0SEEEIMZbNT5AmxthvVTSiSyuDn7gwA+D2nCvlVTZA4OcBT4oi+IV4olsrwxYlc+Hs4I9DTBXNGhvMcNSGEkPZQEkRIJzQahg8OZ6GbjyuiAtxxtaweGgZM7BeMUB9XqNQavPNTJu4fHoodZwvw+oz+kDg5YM2BTMiVajSrNHByEMHNmf65EUKIkNDAaELasfNsAaoaFRCJgHG9AjEozLvdshcKanCjuglTB4bAxdEBAHC9ogHfXyxCrUwJEYDXZ/SHowP1QBNCCBe4eL5TEkSIHt+evwEfVycMDvNBTkUDxvQKMKmevalFmDowBMVSOQ5dKUXC+GiOIyWEEPvExfOd2ucJucXZ3Goo1RpMHhACAAjxlphc131DQwEAPQPc0aRQcxIfIYQQblASRMhfTl2vxKlrVXBxFOPZib34DocQQoiFURJECFoGPx9JL8frM/pxtjsxIYQQYaNRmoQAOHilFFMHhlACRAghdoSSIEIApN6QIrann0XPIRYBag0N0CeEEKGgJIiYTNqksOlZdyq1Bkq1BmlFtRjQ3cvi5/OSOKFORitJE0KIUFASREz2xo/pOJNT3XlBgdqbWoxX//cn9lxsmcZuab7uTpBSEkQIIYJBSRAxSZ1ciehAdxzNLOM7FJNllzfg5Wl9MSTcR7vAoSX5uDpD2qSw+HkIIYQYhpIgYpJ9l0owdWA39OvmhYsFNXyHYzRpkwI+bk4I8HDBvUO6W+WcPm5OkDZRSxAhhAgFJUHEaJ8cuwZnRzF6BXlg5tBQfH+xCEq1hu+wjHLoSimmDLB8F9jNfNycIZW1tATJlbRwIiGE8I2SIGKUsjo5XBzFmD0iDAAgFouw6PaeeO9AJlQ2lAhdK29AzwB3q57T180JNY1KyJVqvPNThlXPTQghpC1KgohRdpwtwAPDw3SORfq74+7B3XDwSilPURnn0g0pBof5WP28nhIn1MmVyK1sRGGNzOrnJ4QQoouSIGKwr3/PR79uXvB1d27zvWHhPrhSXMdDVMY79NfCiNbmIBZBw4CcikY4OdCijIQQwjdKgohBiqUyNMhV7Y6jEYlEcHIQo1kl7LEuGg2Do7glVr7kVTWih79uV1xFfTNP0RBCiP2iJIgYZOfZAsyLi+ywTFxPP5zLFfZMsavl9egd4snb+UUAmpVq+Hs4o6FZpT3+3I6LNjWmihBCugJKgkin1BoGBsDDpeP9dgeEeiOjRNhdYufyajAy0rLbYxiim7crSqQt44JKa+WolSlRRq1BhBBiVZQEkU5dLKjBsIjOBxJ7uzqhXi7sdXDK6+QI8Zbwdv7WTUa6+7ii6K8k6GBaCRbd3hOF1U28xUUIIfaIkiDSqdPXqzA6KsDo92kEslnov3/J1v4/31udNTWr4O/hgu4+EpTUygEAZfXNGBHpq02KCCGEWAclQaRDtU1KNCnVcHU2cFsJkQiMMag1DG/8eMWywRmgoVmFrWfy0axSo7RWjmAeW4GAlmnyUYHuCPKUoLyuGSW1MoR4SdDN+++kiBBCiHV0PMiD2L1PU67hqfhog8u3PswZgPP5/A+STskqx/3Dw1BQ1YRCqQz9u/E3KBoAogLd0TvYEw5iEdSM4WBaKe4e1A0SJwc0q2hgNCGEWBO1BJF2ldTK0M1bAh+3tusCtad3sAeyyxuQX9UIb1cnKHh+sKcX12HG4G64XtGIq6X1uC2Y3yToniHdEez1d2tUWV2zzmtCCCHWQ0kQade+SyWYYeTmor2DPZFZUof8qiZM7BuE6xUNForOMA5iEXoGuCOvqhGNzSp4SZx4jedm0iYFuvu0TYBu0ABpQgixCkqCiF6VDc2oalQgwMPFqPd5SpzQ0KxCaa0c8X2CeJ0yXytTwtvVCe4ujmhSCG8Rx2AvCabetPikCABjDAnb/qANVgkhxAooCSJtMMaw8Ug2lk3sZeL7W6aC9wxwR34Vf60aRTUyhPm6AgDUGg1EImFtVfHMhF4IuqkrLMDDGb9dq0SkvxuulfPbgkYIIfaAkiDSRkpWBaYMCIF7J4sjtifU1xVltfK/9srib056YU0TQn3cAAB1MhUi/d14i8UQob6uSD6Zh+fv7I10gS86SQghXQElQaSNk9cqMbaXv8nvHxruAzeXv6fUM54SoSLp3y1Bkf5u6MPjdhmGCPVxg5uLI3oHe9Au84QQYgU0RZ7ouHRDikFh3mZ1HfUO9sRjo3sAAIK8JCiv52cGVE2TEj5uLQOhHxvdQ/A7t/cMcMeKyb1b7j3fqzoSQogdoJYgouOntBJMH9TNrDpaZ2QBQP9unhbp2jmfV93mmL4VqluTOWdHseDGBN3K2VGMyJt2l+erBY0QQuwFJUFE62pZPSL93OHowN3Hok+IF7JK6zmpS9qkQE2jAgCw89yNNt9ftuMiJ+cRgjA/NxTWyJBb2YgLBTWobRL2nmyEEGKLKAkiWrsvFGL2iDBO6/RwcURTs8qsOn76swQAcDSzHKdzqqBSa1Bwy1o6CpUGacW1yK1sNOtcQtG/mxfe2peO37IrUCtTYuMv2dQyRAghHKMkiAAAZAo1PF0c4ezI/UfC190Zh6+U4vmdF6FSG7+C9JZTeaiVKXGtvAGltXJUNihQL9dNrKRNCswbFYG9qUUAWvYMczd0vzMBGtDdCx8+NBTzR/fAhD5BmNQvCAfSSvkOixBCuhRKgggA4Hx+NUZE+lmk7gVjekCmVKNngDsaTGgVYgy4UlQLB7EIUpkS5fVy3NpjV9OkRKiPGxxEIpTXy5FWVKsdl2SLRCIRPG5aomBMrwBcyK+h1iBCCOGQQUmQWCyGg4ODQV/R0YZvtkmE40K+FMMjfSxSt0gkwn1DQxHm69amBadVe3uMqdQaDA7zRlpxLVqHNVfUNyPIU3e2WU2TAr5uTlgwtgc+PXYdB9NKMalfMJeXwbtBYd7IpkUUCSGEMwZNkR8yZAg2bNjQaTnGGF599VWzgyLWp1Cr4eJo2e4jT4kj6uT6B/i+ue8KVs8c1OZ4Sa0cvYM9calQil5BHqhpUqK8vhmBt2znIW1SIMLPHZ4SJ0wZEILewR4Qi4U9G8xYcVH+OJhWit48bwJLCCFdhUFJkL+/P8aPH29QhW5uwl6Vl7RVXic3eo8wU3hKHFEn098SdD6vRu/xIqkMob6uOHGtEn1DvHA6pwqV9c0I9nKBWsPg8FeiU9OkxJDwljWBRkebvtCjkAV7SVBeL+c7DEII6TIM6g47cuSIwRUaU1YfuVyOBQsWIC4uDiNHjsThw4f1lquvr8f8+fMRExOD2NhYLFiwAA0N1FVgiiMZ5VbpOvKSOKFeT0uQXKlGUY1M76ahhX/t/3VX/2D07+YFEQClWgMvVyc0Kv5OqFq6w5wtGb4giEUiveshEUIIMZ7gBkYnJSWBMYYzZ85g+/btePjhh1FWVtam3FtvvYX8/HycOXMGZ86cQX5+PlavXs1DxLbvRk0Twv0s34LnKXHUOyboRnUTBod7o/qvNYBuViKVoZu3K+4d0h3ebk7wdnVCrUzZsjN8899JU7NSA4mT7c4GM9RtwZ64Ws7NukuEEGLvTEqC7rvvPq7jAABoNBps3rwZixcvBgD07t0bw4YNw7Zt29qUvXLlCmJjY+Hg4ACxWIzY2FhcvNh1Fsuzll8yyjAkzNsq52qvJSi3shEjInz1JkFKtUZn2n4375ZtONxdHE2aaWbr4nr64Vxu29WyCSGEGK/DMUHvvPNOm2OMMWRkZFgkmJycHFRVVaFv377aY/3798f58+fblL377rvxxRdfoLGxESKRCEePHsW0adMsEldXxRjD6etV+OeM/lY5n0c7LUF5VY0YHumLKj1J0K2CvSVwc3aEu7MDGu0wCQr0dEFFQ+f3iRBCSOc6TII2bNiADz74oM3aJC4ulhlE29rt5e39d8uEj48P0tPT25R95plncP36dURFRYExhlmzZuFf//qXReLqqvKrmtDbijurOzmIodIznqWhWY1If3ek3tA/OPpmIV4SBHm5wN3FUWdMkL0Q+v5nhBBiSzrsDhsyZAgmTZqExx9/XOcrNjbWokHd+ote3wJxb7/9Ni5duoT8/HwUFBQgJycHycnJHdZ581dSUhLHUduek9crMUYgM6n83J1RdUsLh1rD2nwWgr0kmDMyHB4ujmhsbjuQmhBCSNeUlJSk8xznQodJ0M8//4zu3bu3Of7FF19wcvJbBQUFAQCkUqn2mFQq1R6/2UcffYSlS5dCIpFAIpEgISEBr7/+ert1M8Z0vigJAoqlMoT5WndJg1vTWY2m5efhpaerrLxejmAv3UURW3eod7PT7jAAEAE0Q4wQYndaJ061fnFBULPDoqOj4efnh6ysLO2x9PR0xMTEtCmrUCjg5OSkfe3k5IT6epo1Y6h6udLiiyMaIqeyEVGB7nqz+qKaljWC9PG4qTuMMQZ76iUK8ZagtI7WCyKEEHMZlQQlJiZaKg4ALdtzLFmyBF9++SUAIDs7G6mpqZg3bx4yMjJw5513Qq1u6QKZNGkSvvnmG21GuH37dkyYMMGi8XUley4WYdawUKuf99Zc5Y/8aoxsZ8+y1jWC9HF3cdS2BNU3q3T22erqIvzckF/VxHcYhBBi84xKgk6fPm2pOLRam7vi4uIwd+5c7NixAyEhIaitrUVmZiaUypYp1p988om23KhRo6BSqfCf//zH4vF1FcVSuVXWBtLncqEUaUW1AIAb1e0nOkVSGUJ99H/P1ckBTYqWhLi6wT4WSmwV4eeGgupGvsMghBCbJ7g/nyUSid4BznFxcSgqKtK+DgoKwjfffGPFyLqOq2X16BPiwcu5GYBjmRUI9XXFwNCWWYDtDXBrVqrbXQBRLBahtUv4+9QiPBoXaYlwBam7jyuKpXIwxlDZoECgp+W3PCGEkK5IUGOCiHX8nF6GO3ncYb2hWYmyOjnkSjVcnc0bl3S9ogHerk5W2ftMKBzEIjDGcOhKGT4/fp3vcAghxGZREmSHGppV8JI4dV7QApwdRJA4OUCh0rTp7nJ2FKNZZdy09wN/lmBubATXYQqeUsNwPq8a7mYmkYQQYs+MSoK4mpJG+JNdVo9egfx0hQGAp8QJA7p7gaHtwOdADxdU/rVWkEbPGkH6NKvsY8+wW6k1DA+ODOc7DEIIsWmcDYxuaqLZKrbg54wyTOrPX1dYXJQ/RkcHAAAKb9m4NdDTBRX1zQCA9JI69O1kNWu1hkFsT3Pjb/KPaX3RJ8QTdrU2ACGEcIyz7rAZM2ZwVRWxoMZmFbxd+ekKA4A+IZ7a85fXNSPwprE8gZ4uKP9r/ZszOVWIi+p4Neucygb06+ZluWAFrLWVrGWWnH0uGkkIIeYyOgn65ZdfMHjwYLi6usLBwUG7i/vx48ctER/hkFyphkQACyQCgJNYBIVaA7H475aMIE8XVDS0tARJm5Twde942ruDWIwRkb4WjVPoQn1dUSyV8R0GIYTYJKOnyC9btgyrVq3CiBEjtCs2M8Ywd+5czoMj3Nl2Jh+9gz0xMMy788JWEOjpgvP5uhumtu4fplBp4CDuvJvnkdgIu58eHuojYBHPfQAAIABJREFUQWGNDL2CrLcRLiGEdBVGJ0FhYWF46KGH2hzfsWMHJwERyziWWY69qUX4z/yRfIcCAAj2lsDJQbch0tFBDLWG4UpxLQaFdp6sjRbI5q98CvVxw9HMcr7DIIQQm2R0d9j999+PI0eOtDm+fPlyTgIi3GtsVmFED18MCfPptIvJWrp5S9pdKTqtuA6DBNJiJXSBni4or6d9xAghxBRGtwSNGjUKs2fPhlQqhbd3y4OKMYaysjLOgyPcyCytQ78QL0zoG8R3KFpRAR7wGtZ2gDYDUF4nR5Cdd3MZykEsAm0oTwghpjE6CVq4cCGeeuopDBs2DI6OLW9njFFLkIBdKa7D1IEhfIehw9lRjO7t7AsGtL+VBiGEEMIVo5MgHx8frFixos3xzz77jJOACPcq65sR5CnhOwyDOYppIXNCCCGWZ/TT5o477kBmZmab41u2bOEkIGLffFydcFswfyta2yIHkQhq6hMjhBCjGd0SdOLECaxfvx633XYbfHx8tMdTU1OpNUiAlGoNHGyoZWVkD1+EeNtOq5UQtK60TfeNEEKMY3QSVFpaio8//ljnGGMMJSUlnAVFuJNd1mBTLSuDw3w6L0R0dPOWoKRWRkkQIYQYyegk6LXXXsOjjz7a5vjNrUJEOK4U1yK2px/fYRALCvGWIK+yke8wCCHE5hjdT6IvAQKAmTNnmh0M4V5BdRPCfd06L0hsVktLEK0VRAghxjI6CXr22Wfx5Zdf4tKlS1Cr1QAAlUrV4Q7zhB+MsZad1g3YgoLYLm9XJ0hlSgDAgT9LkFZUy3NEhBBiG4xOghwcHLB582aMGTMGHh4eiI2NxaJFizB9+nRLxEfMkFvZiJ4B7nyHQSzs5jWVMkrrcTa3msdoCCHEdhg9Jmjjxo0AALVajfT0dBw8eBBr1qzBwoULOQ+OmOfXqxW4e1A3vsMgVlbTpOA7BEIIsQlGJ0GtHBwcMGjQIAwaNAixsbHa5IgIR1WjAkFeNGOIEEII0YeTBWTGjRuHlJQULqoihJhABECh0sBRLEKQpwvK6migNCGEdMbolqChQ4ciNjYWMTExGDlyJAYNGoTjx48jMjLSEvERExVLZejm3f7eXKRrCfBwxoWCGkT6uyE60AMXC2owdeDfXaGMMdTJVPB2a7tpLSGE2CujW4JWrlwJLy8v7Ny5E5MmTYKLiwsmT56M0NBQ/Pvf/0ZKSgoqKystESsxUEmtDMmn8jA8ktZushch3q44ea0S0YEe6BviiYySep3vVzUq8MR/z4Ex2l6DEEJaiZiZvxXz8/ORmpqKixcvar+Kioq00+f51jpzxp5++X98NBuT+gejT7An7cZuJ9KKavGvvWnY9sQouDk74v2DmXhpal/t98/lVeNYZjkGh3nrtBDdqrpRAUcHEbwk1GJECBE2Lp7vJg+MbhUZGYnIyEjcd9992mNVVVXmVmsX9qYWwdXJAZMHhHBab7NKg74hXpzWSYQtxFuCygYF3Jxb/kmH+7khv6oRkf4tSyTkVDTgkVER+OpkHib1C4ajg/5G4MNXSsEAzI2NsFbohBDCG4vsrOnv72+JarsUuVKN9JI6/J5bDaVaw1m99tTiRf7m5+aM3sGe2tdTB4TgYFqp9nWxVI7u3q4YGu6DnA622CiuleNGdZNFYyWEEKGwne3Fu5grxXWI6+mPx0ZH4qNfsjmrt7BGhjBfGhBtb8RiEV68q7f2ta+7s3YVaQBgf5WJCnRHTkVDh3U5iEVQayiZJoR0fZQE8eRKcS0GdPdCpL87HMQiqDhqDbp4Q4qh4b6c1EVsS//uul2gvYM9kFWqO0C6h787cis7bukZGOpNW28QQuwCJUE8qahv1i5kGO7rhsIamVn1VTU04+S1SmSX1eO2IA8uQiQ27q7+Ifg5vRRKdcv6QQDg7uIImUKlt3xVQzP83JwwOtofp3P+Htf3c3oZ8qtol3pCSNfTYRJ04cIF3HHHHZg4cSIOHz6sPT579myLB9ZVXSyogeaWroboIA9c76SLojOFNTJ8dTIXjIE2TCUAAA8XRzQq1LhR3YRwv867SLPK6tE7xBNeEifU/dWVplJrsPlEDoqk5iXphBAiRB3ODnvhhRewadMmODs74/nnn0deXh6efPJJmv1lhvVHsnHfkO7wd3fWHosOdMfZ3Crc2S/Y5HprZUqIRCJ4u9LUZvK3IWHe+P5ikUGfraul9bhnSHcAwNheAdh6Jh8eLg6YOSwU9XL9rUeEEGLLOmwJEolE6NevH6Kjo/HDDz/g1KlTePvtt2ntGROoNQyHrpRi2sAQbPr1OmYNC9N+z1PihAYzHzJSmRKLb++J6YNpw1Tyt/g+Qfjuj0L0DHTXHvNydUJtk7JN2epGBfw9XAC0JEH+7s74Pacat/cK0LYMkc7RDE1CbEenY4JqampaCorFSE5ORklJCU6dOmXxwLqaa+UN+LOwFrOGhSJ5YWzb7QvMTCxrZUpEB3qguw/NDCN/kzg5YOn4aJ3FD3sGuCPXgDE+dw/qhnfvHwQviRO1BBnh/UNZOJdXzXcYhBADdJgEffXVV22Offzxxzh06JDFAuqqqhsVGNPLHxInB72Jirlta7VNCuoKI3o9PqaHzuueAe7IrTRsDJpIJIKHxBENzX8nQVvP5HMZXpdyuVAKb1cnpBZI+Q6FEGKADpOgqKgo+Pq2nW49fvx4iwXUVdU0KeB30zigW5nbgK5UMzg70mQ/0rlwPzcUVOkOdJYp1JA4O+gtf/O6QXKlGmdyaEygPkq1Brv/KMST46JQJ7et7sNr5Q2U3BK7ZNRTMzEx0VJxdHnVjQr4ubWfBJnj1DXasJYYzslB3GaV8rI6OUL+WrKhIzVNCpTVyi0Vmk3JLqvXSXa2nMrDY2N6QCwWmd2ya02/51Rhb2oRGptVtFo4sTtGJUGnT5+2VBxdnrRJAZ8OkiBTf2nKlWp8/XuBie8m9oqB6QzgLauTI9iQJKhRqdM1Zk8+OXYN5fV/J4A/Xi7Bx0evgTGmXeIiOrBljS5/DxdU1DfzEqehrpbVY8ORq7haVo8X7+qNx0f3oN8lxO5Q/4mVWKq7StqkRG4He0ERok9UgIfOHmKlBiZBUlnL2DN7nAF1JqcK18pvGkvFGCb2DcKu84X476k8nbFXQ8J9cOmGcMcFbTp+HX8W1mLJuCjMH90DIpEIrs4OiApwx5ViWi2c2A9Kgqyks0eGqY+UmiYFVBruNmAl9mFML3+cuv732J7yumYEe7l0+j5pkxI9/N1R24WmzNfKlJAr1Z2WuS3IEzkVun9wxEX5o06uxAMjwuDk8Pev037dPJFeUmeReM0lV6rR0KzCAyPC4O6iu1Tc/cNDsf33Ar0/3yaFCs98fYGzLX4IEQJKgmxcTZMCvYM9UdOk4DsUYkO6ebuitPbvwdH1zSp4uHS4diqAls9b326eKKsTdlePMb77//buOzzKKv0b+HdKkkmvpEAqhBAChBSEUBMIQaqguxYEIYoIrvguZVUUd82qq6591cWGFEHYn32VtVAEAlIElCIJSSANSE9In0mZOe8fMSOTSZkkM5kk8/1cV66LeZ7zzNwnD5O559TTVzuc0n48sxTzIwaiqLKpO0ytEdqV2e+fPBjhvi465W3kMr1xV73Fd78WYNbI1tcTk8ukWJsQgncPXcanp6/qtPh9cDgLt0QMxNfn8noqVCKT61QSZIlN4D1Fgq79fstrGxDl7wpbq9Zn9hC1RQKJzhYuhiyCWl7bgGFejiis7LuDoyuUDSi4YXB3dklNq13KP2WVacc/XbhWgZGDnLXn8sqVBq3J1dF7ukGtwTfn8/HKnjTtIOsGtQbbj2Vj85EsHEgr0ttmpyM/XCxEbmnbA5xT8yv1Ntu9kbuDDR6dGYpAdzs8+79U1NY3IiWvEi721rh5hDdS8ipbbQ0qqe4/iTFZDqMNjK6tNc6sApVKhcTERMTExGDMmDE6e5a1dPToUcTHx2PKlCkICwvDm2++aZQYzMFKJkFjJ//YAU3fzKMCXNsddE3UmrCBTp3usqlr1MDPzQ4FfTgJ+v5CAT4+dUX72NnWCiXVv7ekltfW4/lvU5FdUoONBy4BADSiaamAZlklNRjs8fsq3K0JdLfH5eK2x+vVNarx+OfnMdDFFvdNDMJre9Ox+1went2dggnBHrg1chAqlQ3Yf7GoU/U7nlmGXSdbH+CcXliFYAM3WB4T6IYVsYPx/DcX8dGJHNw91h8AcMvoQfjqrG5rkFoj8MgnZ3/bv7Dp75gQosNuRiJzM1p32Ny5c43yPElJSRBC4Pjx49i5cyfuuusuFBYW6pXLysrC6tWr8X//939ITk7GJ598guzsbKPEYGzf/VqAIQPa/4NpI5d16Q9GeW0Dhvs44v7JQV0NjyzU2CA3/JRl2MrGNlZSqBrUkADwdLLRdgv1Rbmltdr3WlGlCp4txkLt/CkXK6YMwR03+WHkIGfsPJELH5emQePW8qbfQ2ZxNYI6SIJmjPDCjuM52pacokqVzkbJJzLL8IcoX0T4ucDV3hp/mxsGF1tr/HVuGIYMcICrvTVuGT0QxzNLW21RUrfypalK1QBnWyv4u9khraAKAFBT14gPj2Xj41NXsOXHLMwNH2jw78rTUYEn5w7Hn+OHapPAUb7Oeq1Bp3OuY8n4QIwY6IyX96ThgyNZeOG7i/j71ykcQ0S9WqeToP379yM8PBy2traQyWSQyWSQSqU4dOhQt4PRaDTYtGkTli1bBgAICQlBZGQkduzYoVf21Vdfxb333gsPDw8AwIgRI/DKK690OwZTuJBXgfkRg9otY2MlRV1j5/9Y1DVqYCOX6QzKJDKEh4MNymrqDeqGdVRYoVLVAIGmhL1e3be6xut/e29pNAISCeBiZ4Xy2nr8nFuOSD/dBWGV9WrtwqazR/ngzJXrmBw8AAAQ6GGP7NKapnW/2ln8FGj6nf0x2hdbjmYjs7gaGw9exmenr2rPn8q5jjGBv7+2RCLBpKEekN/wXpZIJBgb5IaT2dd1nvtiQSXu33ZSLxE6mFaMqcM88cdoX3x6+gqKqlR4/ttUTB3miQlD3PHwtKGwbWNRzLbYyGXwbDFzcE64D76/8PuX0yOXSjAh2B1jg9ywcKw/5ozyweOzhmPZpCBs+TG7U69H1JM6/cn58MMPY8OGDTh//jwyMzO1P2PHju12MJmZmSgtLUVoaKj2WFhYGE6dOqVXdv/+/aivr8fs2bMxceJE/OUvf4FK1Tu/nRoy1kIhl7WaBNU3avBe8uVWr+nL4zKo96hQNsClgy1XnBTybu0fll+hxK/XTD/1OumrC3qJwa/XKrDuk7MAgNSCSoR6O2FisAeOXCrBxYJKhPo4wkoqQV1j6y2xL/5xNPzd7QA0rQPUPEPMkPf1yEHOcFLI8cUv1/DknOFQWMlQ89s4I7VGY9CXl4ThXtibUqBz7Iufr2FNQgh2/qTb7ZWaX4nhPo6wkkkxytcFb+zPwBOzh8PPzQ6+rnZG21sw0t8VZ640JWYajUCjuunLGAD4utrB27kpaQr2dICNlZTT7qnX6nQS5OvrizvvvBPBwcEICAhAQEAAAgMDsWvXrm4H09zt5ez8+wBEFxcXFBXp94lnZ2fjnXfewYcffoiDBw8iJSUFq1ev7nYMxmboYGcbKynqWukOyyqpQXph6/s8PfrpOQ5Wp27xclbg7NWKDtcIclJYdWsn+VPZ102+3UajWoP0wiocStf9e/H9hQLEDHbDhbwKHLtciglD3DHc2wkX86ug1ghYyaTwd7fDlbJaVNQ2tLsHX9O+azWd2vD49jF+WDdjGOQyKW4e4Y3vLxT81p1m2NgcqVSCEQOdcSi9GACQVlCFIA97hPu6oFLZoF3lua5RDSuZVJuc3TJ6IJ5dMAp21h3P+uuKKH9XnM65jjNXyxHh59JmucXjAvDZ6Wt464cMfH02D1V9bEsR6t86nQTddttt2Ldvn97xNWvWGCUgQP8bVmsf9HV1dbj77rvh4eEBKysrrF69Glu3boWmjTVzJBKJzk9SUpLR4m1PeW0DXFvuGN8KG7kUqgb92NMKq7RN+QBw4GKRdjxDhbIB6YVVxguWLE6Uvwu++zW/wzWCHH9rCerqyuZZJTUorzXth19GUTUWjQvAiczfxzk1j2O6Y4wfvvzlGq7X1sPV3hpSqQRqIbTjXAZ7OCCrpBYX8iranTmlsJKhUtkAmy4ufDrM2xFphVU4mFaM2JABBl+3IHIQcktr8NXZPHz+81XcFuULAFg2KQjvHLqM45ml2gSvp8wY4Y09KQU4nF6CKe3URSqV4G/zwnD/5MEI8rDHh8e4Rxl1TVJSks7nuDF0+p08btw4rFixAu7u7hg8eDAGDx6MoKAgo+ws7+npCQAoL/99pdXy8nLt8Ru5urrCy8tL+9jX1xd1dXUoKWl9Hy0hhM5PTyVBRVV1Bq3EayOXtdocn1VcAz+335uwj1wq0X6YjBvsZtBzE7VlmJcjjlwq6bglyNZKZ58sqUR3YO75q+13dzSqNdp1dUzl7JVyjPZzhr+7HbJ/m/L+/YUC3DzSG1YyKYZ6OUIu/f1PXqi3I8J8mhKeQA87ZJVU40JeJUYMdG71+ZtdLVci4Lfusa7wdFQgo6gKAxw7XpzyRveMD4RUAowY5KxdfV5hJcMz80eiUtmA3efyER2gv+G1qcikEgx0tsX12nooDFiiQ2Elw8hBzqipa9T5YkdkqOaJU80/xtDpdtJ7770XDz74ICIjIyGXN10uhDBKS9CQIUPg5uaGtLQ0beKTkpKC2bNn65WNiIjQ6SYrLi6GtbU13N177ptQW7745Spmj/KBjVyGwkoVPA34Y9fWwGi1RqPzh7ugUoUqVQO8nGxgI5dh/cxQvWuIDCWXSeHvZqc3Q6olR4UceeVKOCqa3vPuDjYoramDp6MCao3Ahi/P46tVk9p/MRN33eaVKzHIxRYLIgbh719fwMyR3jh75fdJCX+I8kVN/e/jmm4Z/fssKUeFFarr1BBCtNsdBgADHGw6nBnWnnnhPtib2rWWpNZmdkmlEswY4Y0ZI7y7HFNX/THaF1evKzsueINZI33w7a/5HU4WIeoJnX4nuri44C9/+Qvi4+MRGxuL2NhYxMXF4e233+5+MFIpli9fjs2bNwMAMjIycObMGSxatAipqamIj4+HWt3UWrJ8+XLs2rVLuz7R5s2bsXjxYshk5l808NjlUuw43jRg0dCNKW3aGBjdcuxBcVUdKlWNqK1Xw85aZrQmQbJc62YM0w5qbYujwgq5ZbXa9ai8HG1QWNG0OF5uWS3qGzWoaKO7S9WghrVcCplUavLp0hKJBPY2cjwxezhc7axx9zh/7TmZVAInhZVO2ZbvH0PeTQsiB2k3Su0KTycFFo0L6PL1vYm9jRzDvB07dc0oX2dcyOudW4qQ5el0EjRlyhRcvHhR7/i2bduMElBzc1dMTAwWLlyIXbt2wdvbGxUVFbh48SIaGpr+0N55551YvHgxYmJiMHHiRMhkMrz++utGiaG7vJ1tUVFbj+q6RuSW1RrU7G3z2/ojLUnQlAc1N/3ZyKWoUjWgStWo/VZO1B1R/h13odhby5BXrtKObwv2dEBGUdN4tPTCKiwc64/Tua2vOXS5uBrBng4Y6KJAXrlpZjQ2J1rNXOysEenvavDCgEDTwGJbAwYRR/i5GNT9Q20L9XZEai/dW40sS6c/RQ8fPozXXnsNQ4cOhYvL7zMCzpw5Y5TWIIVCga1bt+odj4mJwbVr13SOPf7443j88ce7/ZrGVKFsgJNCjlmjfPDy92kIG+hkYH+5fneYqkENGysprGVSVNU1wtFGDieFFapUjahSNeh8qyUyJYlEgrKaerj8lgQFutvjy1+a3o+Xi6uxaGwAth7NxrRQL71r0wurEO7rgtLqeuSU1WinmxuTIWN5OuJqZ43QTrZqUNfMHuWDv/33V9w9LgAjBzrprI1E1JM6nQQVFBTgrbfe0jkmhEB+fr7RgurLsktqEORhj0Eutki6ZYTB19nIZXpT5Iur6jDAwQYSiQQVtQ1oVAv4udmhStWISrYEUQ+rVDVou8OkUgmaR/io6tVwtrPS2zD01T1pWBE7BDmltZgXPhAONnLsTSnE5KHtv05RlUr7/95Q566W64zx6Yo5o3w6XACRjENhJcOGOWE4kVmKN27YFsTVzgr3xAQwKaIe0+lP0Q0bNmDx4sV6x29sFbJkl4urMbqdNTPa0trA6OLqOgxwtIFaI/D1uTzklSsR6eeKkuo6VKka4MiWIOpBCrkMrjfsUSeXSnVmNFrJmh7byGW4kFcBexs5/n3gEuRSCeQyKTwdbVBU1fEmm6/tzcDNI7wQN0x/Vmh9owYNag3sW+x4n1+hgrtD52ZbteTnZvwWKmqbs62V3oDujMIqPL07BaumBuutUk1kCp1Ot1tLgABgwYIF3Q6mP7hcXI0g987PHGlt77CiyqYkyMXOCnsuFKJa1YhAD3vU1qtRpWqEE1uCqAc521rpzJwa7uOIC3mV2qnv4b7O2qnyX53Nw32TghAd4KpNfCRNg9s6fB13e2scyWh9qYv3D2fqrTPzS+51jGhnbR/qO4Z6OWL9rFB8cCQLJ7MN29eOqDvY5mhEzYOXu7Ieio287ZYgFztr+LnZ4e+3jMSoQU3jHirZEkQ9bJCrrc5O6hF+Lth9Nh/+v7WgxAx2x77UIryfnAlfVztYyaSIH+6Fv80L+/1JOujiat7QdLSfC37O1d0vK72wCvbWMp31igDgf+fyO7UpKPVudtZyrJ8VitT8Sp291ohMgUmQEaUXViPEq2sDK23kUr0FxMqq6+Fub4OBzrZYNikIznZW2hkwnB1GPW1+hG6i4emkwImsUgz1bPo/b2stw/pZoRg/xB23R/tqy7XctqG9Rc5+yi7DTYFumDXSG7vP5mvLqjUCO0/kYnFMAAa62OLq9aalMU7nXMdoPxed5Iz6PolEgiXjA+HuYI0Xvr2I1/am4/3kzFZn0BJ1h0FJ0OHDh00dR79wOuc6xgS6dena1gaBNi/rb2st09ubR/nbOkFEPWXyUP2tEQLc7fSmoY8c5NzmjEhfF1tkl9bqHCutrtPui5VeUIVhXo6Qy6SYE+6Dd5MzodYIvPj9RSz+bcBswnAv7EspRKNag6/OXMOcUT5GqiH1NnHDPLF+VijWJIRgdrgPnv8mlYkQGZVBSdCDDz6I8+fPmzqWPi+/QomBzl0fzKdqsW1GR8uCc6FEMrdn5o+EbSeS8YQwL3z3q+6O6P85eQX7UgpRXlsPgd+7k6MDXDHcxwkPffQz/hDlq022vJ0VKKyqw9aj2bhnfIDJt+Og3mGQiy0ejAvGc9+kQlnPRIiMw6AkSAiBTz75BPPnz0diYiK2bdumt2YPNelOYuLlqEB+RdMS9BoNd4en3q+zM7Jc7a1RrqzXJviVqgbUN2qwNiEEr+5N1+vijQ0ZgI2LovS6md1+m6UW7Ml1fSyJt7MCq6YF4+U9adifWmi0/aPIchk0qOSFF17AvHnzAACVlZU4cOAA/vnPf+Lq1avw8/NDQkIC4uLi4ODQ9aXk+zq1RnR7XMKccB9892sBlk4IxKXfVtltTavbaxD1ERG+Ljh7tQIRfi74+OQV3DXWD6721hgX5I5AD/1p6q219CydEAg5W4AskqejAn+dG4aDaUV4ZncqHp8dCiuuK0RdJBFdTKXVajWOHj2KzZs3Y/v27ZDJZKir63gNkJ7W3DJj6m8MBRUq7L9Y2K09gYQQeG1vOtbOGIbtx3MwfbgnfJxt9cptPHgJ8yMGYZCL/jmi3k7VoMZr+9Lh62KLBrXAfZOCzB0S9VGXiqpxOKMY907k/yFLZIzP905NL0pNTcXevXuxd+9eJCcno7a2FhEREXj00UeRkJDQ5SD6sga1BnWNGlwrV2JgN5OSG7vSiipVrSZAAPCnuOBuvQ6ROSmsZBg50BkjBjphcDc2IiUK9nTA9xcKkGeEv79kmQxKgpYuXYoffvgBeXl5CAoKwvTp0/HBBx8gPj4erq4db77Yn/2SW47csloorKTaqcJE1L553dzigqjZfROD8OreNGyYE9ZxYaIWDEqClEolVCoVbrvtNqxcuRKTJ0+GtTX32AGAvHIlCitVsJZJMSVEfwpxZ8llTbvJSznzi4ioQ7bWMowNcsf+1ELED9ffwJeoPQYlQR9//DGEEPjll1+wb98+vPHGG5BIJJg4cSKmT5+OyMhIU8fZa+VXqFClaoS1TGKUXd0D3O2QnF7c6gBRIiLSlxDmhdf3pcPPza7LC9aSZTJ4SL1EIkFUVBQeffRR/Pe//8XHH3+M6OhoPPfcc/Dy8sLChQtNGWevpWxQa1dxNoYJQzyw+ccsDPPiXkhERIZ6eNpQ7PopF0VVKnOHQn2IQZ/emzZt0v67tLQUH3/8MR566CEsW7YMn332GQBAJuPqxcYwwNEGby+KRqg3v80QERlKJpXgkZuH4c39l1BR29DxBUQwcIp8WFgY5s2bh3379uHs2bOwsbHBpEmTkJCQgISEBIwePbonYu0SU0+Rf3VvOqpVjfB1teVUXyIiM6tQNuCVPWl4YMpg+LpyWEF/ZozPd4OSILlcjoiICG3SM2nSpD4zMLonkiAhBKaEDMBNXdw3jIiIjKeuUY2XvkvDH8f4ItSbQwv6qx5LgubNm4evv/66yy9iTqZMghrVGmw8eBlxwwZgqKdjp/ZQIiIi09FoBJ7enYIn5wyHnCtK90s9lgT1ZaZMgoqqVPj+QiHuien6KtFERGQaaQVVOJFViiXjA80dCpmAMT7fmR53QXFV0/YgZTX1cLfvG92CRESWZpi3I8pq6jljjNrEJKiTjmSU4IkvzgMASqvr4cYkiIio11o2KQgfHMkydxjUSzEJ6qQTWaWI8nfFjuM5uFxcDQ8HJkFERL2Vo8IKYwLc8N8z18wdCvVCTII6SSLAwNC9AAAgAElEQVSRYPJQD3x/oQAH04rhZm9j7pCIiKgdCWFeuF5Tj5+yyswdCvUyTII6QaMRkAAYOcgZby+OxqnsMrjYdn+rDCIiMq2lEwKRnF6MrJIac4dCvUi7SVDzStFZWVmYOnUqhg4dipiYGJw9e7ZHgutt8iqUGOiiAAA42MgR4uUIqZQbnRIR9XYSiQSrpw/F1h+zUKXiitLUpN0kaOfOnQCAdevW4a9//SsyMjKwadMm/PnPf+6R4HqbzOIaDB7goH28jCtEExH1GXKZFH+eHoJX96ZDo+nXq8OQgQzqDrt+/TqmTZsGABg5cqRJA+rNskpqEORhr308a5SPGaMhIqLOcrO3xm2Rvtj8I2eMUQdJ0PHjxxESEoLTp0+jrKxpQJlGo0FlZWWPBNfblHJdICKiPm+UrzN8nG2x8eAlKOvV5g6HzEje3sna2lq9Y0qlEu+++67JAurtmleoJCKivmtOuA+ulSvx2r50+Lna4nptA+oa1Rju44S54QPNHR71kHaToJaysrJw9OhRFBQUIDk5Gd7e3hg/fjwGDx5sqviIiIhMYpCLLZ6YPRxXymrh6WQDG7kMz32TilkjfSDjpBeLYFASVF5ejsTERHz11Vews7ODq6srgKaxQkqlEvPnz8eWLVvg7Oxs0mDNSVmvho2cKwoQEfU3fm522n/PGumNb87nY95otgZZAoM+1VeuXImAgABkZGSguroaV65cwZUrV1BdXY309HQEBATggQceMHWsZpWSX4kwHydzh0FERCYU6e+Ks1fKTbLpNvU+Bu0iP2HCBBw9erTbZcyhu7vMFlWp8OJ3afB2UiBxYiA8HLhCNBFRf3b0UgmUDWrED/cydyjUjh7bRb6urg7V1dVtnq+urkZdXV2Xg+gtXt2Thpq6Rp31I9IKqhAf6om9KYVMgIiILMD4Ie44cqkEjWqNuUMhEzNoTNA999yD0aNHIzExEaNHj9YZE3Tu3Dls3boV/+///T+TBtoTvvm1ABcLqrB6egjCBjZ1fRVX1SHK3xWv3jnazNEREVFPkEgkWDI+EM/+LxWPzQyFrbXM3CGRiRjUHQYAn3zyCV577TWcOnUKjY2NAAC5XI7o6GisXbsWt99+u0kD7arONJet+b8zKK2pxz0xAUgIa2oGfefQZdwTEwB7m05NpCMioj6uqEqFN/dfwtqEELhyjbhexxjdYQYnQc0aGxtRXFwMAPDw8ICVVe/eQNTQX1JtfSO2Hc3BHWN88fXZPCRODEJOaQ0+O30Va2cM64lQiYiol6lSNeCl79OQNG8E94rsZXpsTNCN5HI5fHx84OPj0+sToM4oqarHAEcbuNlbo6ymHgDw5Je/ooariRIRWSxHhRUWRA7Cpz9fNXcoZAJGW/hm9+7dxnoqsyiuVmGAo402s9RoBJT1aqQXVpk5MiIiMqcof1dcKatFSXXfnwBEuoyWBL366qvGeiqzKK6qh4fD732+2aU1mDd6IJxs+09rFxERdc3yKYPx/uFMc4dBRmbQaN/mHeTbc+bMmW4H00ylUmHlypW4ePEiGhsb8dxzz2HGjBltlm9oaEBYWBgWLVqEpKSkLr1mSXUdovxdmh5IJDhzpRwTg90xMdi9S89HRET9h5PCCiMHOuNUdhnGBLqZOxwyEoOSIJlMBm9vbwQHB7dZJjs721gxISkpCUIIHD9+HOnp6YiJiUFqaiq8vFpfuOq9995DUVFRt16zUtWgbfXxc7VFcnoxFkQM4kA4IiICAMwZ5YPnv01lEtSPGJQEvf/++1ixYgW2b9/eZpmDBw8aJSCNRoNNmzbh888/BwCEhIQgMjISO3bswLp16/TKV1dX4+OPP8Ytt9zSrddVNWi0e4OFejvhq7N5TICIiEhLKpUgxMsRqfmVGM5tlPoFg8YEBQYGYv369bh27VqbZbZs2WKUgDIzM1FaWorQ0FDtsbCwMJw6darV8q+88goefvhhyGTdX8yqeVD0UC8HxIYM6PbzERFR/zI/YhC+PNP2ZyH1LQYPjJ46dSoGDRrU5vnAwEBjxIPCwkIA0NmR3sXFpdXuruLiYiQnJ+OPf/yjUV67mcJKhvsnDzbqcxIRUd9nLZfCx0mB3NJac4dCRmBQElRaWmrwE3ambHuaW2WatbYY0jPPPIMNGzYY/Hw3/rQcQM2OLyIiMsTtY/zw6ekr5g7D4iQlJel8jhuDQUlQZ7bE6O72GZ6engCA8vJy7bHy8nLt8WaZmZnIysoyaOYa0JRE3fjTMgnq+nqTRERkSext5HB3sOE6cj2sedJU848xGDQwOj8/H9u3bzfoRdvbbd4QQ4YMgZubG9LS0rSJT0pKCmbPnq1T7tChQ8jLy0NcXBwA4OLFi1AoFDh48CC2bNmCoKCgbsVBRETUlkXj/PHM7hT8bd4IyDiJps8yaO+wwMBAg5ue/Pz8kJyc3K2g1q9fj8LCQmzZsgUZGRmIiYnBhQsXcP36daxatQp79uzRGwidmJiIwMBA/W4uA/cWeXVvOtYmhHQrbiIishyp+ZU4mV2GJeMDzR2KRTLG3mEGtQQZcw0gQyQlJWHlypWIiYlBY2Mjdu3aBW9vb2RnZ+PixYtoaGjQJkH19fWYMWOGtiXo3Llz2un1ncE8noiIOmO4jxMOphUjt7QW/u525g6HuqDTu8j3NWwJIiIiU6lrVOPfBy5DAmDeaB8EezqaOySLYZZd5PsjtUaAXbpERNRZNnIZ1iaEYEXsYJzMvo5/H7hktEG7ZHpMggBU1zXCwcagnkEiIiI9dtZyLBzrjylDB+DF79PQqNaYOyQyAD/5AeSVK+HtrDB3GERE1MeN8nWGg0KOp766AA8HG53lVzQagdF+LkgIa30fTOp5RhsTlJOTg4CAAGM8lVEZ0me4+UgWbokYCA8Hm54Ki4iILNCzu1Pw+Ozh2mn1Qgj85+QVLBzrb+bI+h6zjAm688479Y4JIYy+dUVPKq2pYwJEREQmNzvcB/87n699/MUv17A/tQi19Y1mjMpydToJ+vbbb/Hyyy9rH+fl5SE+Ph4///yzUQPrKQ1qDeRSDo0iIiLTi/J3xdkr5VBrBEqr65BVUoMHpgzGT1ll5g7NInX603/VqlXw9/fHH/7wB3zwwQcIDw/H4MGDtSs39zVnr5Qjws/F3GEQEZGFmBvug93n8vDe4UwsnzIYkf4uOHOlvOMLyeg6nQStWLECc+fOhUajwfLly3H//fdj06ZN2L9/vyniM7kTWWUYG+Rm7jCIiMhCRPq7Ijm9BJF+LnBSWMFKJkUDZ5OZRaeToOnTpyMyMhLXr19HWloaGhsb8fTTT2PJkiWmiM/klPVq2HN6PBER9aBnF4zEzJE+2sfu9jYorqozY0SWqdNJUE5ODhITE3HgwAEMHToUL7/8MsaOHYtvvvnGFPGZVEVtA2ytZR0XJCIiMqKWnz0Tgz1w9HKJmaKxXJ1OghYuXIjHH39cZ0PVmTNnYtmyZUYNrCf8++AlLBrHaYlERGReIV4OyCisNncYFqfTSdC2bdtaPZ6Xl9ftYHpSo1oDWysZXOyszR0KERFZOIlEAgHBLTd6WKcHw0ybNq3V42fOnOl2MD0pvbAaw7y50R0REfUOQR4OyCqpweABDuYOxWJ0OgnKz8/H+vXrtY/Lysrw9ddfY+3atUYNzNR+uXIdccM8zR0GERERAGBSsAe++OUaHoxjEtRTOp0EPf/881iwYIHOsVWrVuGee+4xWlA94ep1JQa52Jo7DCIiIgCAt7MCQR52eGZ3CuaNHsg17HpAp8cEtUyAAMDKygqXL182SkA9oVLVwF3jiYio15k50gdPzB6O7JIavJ+c2W7ZjMKqHoqq/+r0BqpPP/20zmOVSoXTp09DIpHgu+++M2pwxtDaBmtf/nINw32cOCaIiIh6rf2phSivbcAfon31zp2/WoGkry/gswcnmCGy3sEsG6i+/fbbyMrK0v6UlpYiISEB//nPf7ocRE/LKKpCiBf7XImIqPeKH+6FerUGhzOK9c598cs1zA33QVGlygyR9R+d7hNasWIFkpKSTBBKz5FJJDrrHBEREfVGC8f6498HLsHOWobogKYtnpLTizFhiDt83WzxU3YZ5oYPNHOUfVenW4LaSoAeeeSR7sbSI4oqVRjgaGPuMIiIiAzyp7gh+DmnHPtSCqHRCBxIK0L8cE+EeDoivUB3XFB9o4bbb3SCQS1B9913X4dlvvvuO7z00kvdDsjUzl+rwIhBzuYOg4iIyCASiQTLpwzGJ6euYN0nZ7E4JgASiQQSCaBpMRzmf+fz8OOlUiyOCeDsMgMYlAQdOHAAiYmJ7ZZRKBTGiMfkjmeW4tGZoeYOg4iIqFNuH+OH6ABXncUUXeysUFZTDzf7pt0PLhZU4aU/huODI1nIK1di9iiftp6OYGAStHLlSjz22GPtlukLSdD/nczFpKEDYCXrdC8gERGR2bVcTXpskBtOZpfh5hHeuFJWC19XO0gkEtw/eTB2n8vDf89cw/yIQWaKtvczKBvw8fHBhx9+iPz8/DbLdJQk9QZXypSIDRlg7jCIiIiMIszHCRfyKgEAX5/Lwy03DJKeGz4QeeUqpLeznlB5bT0uFVUhJa8SF/IqUF3XaPKYexODkqDnn38eQNNc/NzcXO1PX5JXroSXc+9vrSIiIjKUXCaFRiOg0QhUKhvhbGelc/6BKYPx0fEc1LSS3FSpGvDS92lIya9CblktrpQp8dHxHLy6Jw3//O4iSqv7/wBrg7rDvL29sWTJEgDAvffei++//x4zZ87E5s2bTRqcMW07lo0l4wPNHQYREZFROSrk2JNSgAlD3PXOyaQSrJo2FG/sz8D6WaE6y8O8l5yJtQkhcHfQnzFdpWrA6/syMDfcB5H+riaN35wMWjF62rRp+OGHH7SPp06digMHDpg0MGMxxoqSREREvdXpnDI8981FfLxiPGTS1tfA+yX3OvakFGL19KGwkctw9FIJSmvqMW9022sMCSGw7Wg2XO2te+W4IrOsGH3jC99ow4YNXQ6CiIiIumbUIBcsjvFvMwECgEh/V9w91h/P7k5FTmkN9l8swtzw9meOSSQSJE4MwtXrSuRXKPXO51co+3wDg0HdYRUVFTh27Ji2si0fA03rBP3jH/8wTZRERETUKmu5FLdG6u8v1pKfmx0enx2KV/ek44Epgw3eOSFxQiDe/OES1s/6fXmZ/53LR1pBJRo0ArZWMswe5YNgz763HZVB3WFSqRQSiaTdjE8ikUCtVhs1OGNgdxgREVH3fHr6KoI9HRDh54Kc0hp8+Use/jx9KABAWa/G2wcvYcYIb4xsYzHiKlUD3j2UiXvGB8DLyTiTlIzx+W5QS9C4cePa3SBVCIGFCxd2OQgiIiLqvW6LHISnd6cg1NsRmw5n4a9zw7TnbK1lWJMQghe/T4PCSopgT0eda4UQ+Ne+DCydEIhPT1+FjVyKxTEBUFjJeroaegxqCUpOTsaUKVO6XcYc2BJERETUfT9lleGtA5fwjwUj4edmp3derRF4/ptULJ0QqHP+oxM5GObliDGBTRvA5pTWYOvRbDw1b0S34jHG57tBSVBfxiSIiIjIOPIrlPBxtm3zfH2jBu8lX0ajRuCOMX4orqrDL7nXkTgxSKfcf37KxZhAt26NI2ISZAAmQURERD2rUtWAT05dxfWaeqybEaI3CLumrhHvJWdiTUKIzvH8CiW8nRQGDdrusTFBRERERIZyUlhh2aSgNs/b28ghhICyXg1b66axQVev1+KtHy7B1lqGReP89cYWmQJbgoiIiKjHpRVU4dzVctw+xg9CCDy9OwWPzQyFTCrBzhO5qK5rxJ/ihrTZKmS2xRKJiIiIumOYtyMyiqoBALvP5SNhuBcUVjJYyaRYOiEQ44Lc8NnP10waA5MgIiIiMosRA53w46USpOZXYkKwh865MYFuuFJWi8JKlclen0kQERERmcXMkd54eU8aHpgyuNXzK2IH491DmSYb0sIxQURERGQ2jWoN5LK222SOXi5BcVWd3iauHBNEREREfVp7CRAATBjigaLKOqg1xm/MYEsQERER9Tn9siVIpVIhMTERMTExGDNmDPbs2dNqua+//hozZszAtGnTEBUVhTfeeKOHIyUiIqK+rNctlpiUlAQhBI4fP4709HTExMQgNTUVXl5eOuXWrFmDTz75BJGRkSgoKMCIESMQEBCA+fPnmylyIiIi6kt6VUuQRqPBpk2bsGzZMgBASEgIIiMjsWPHDr2yDz30ECIjIwEA3t7emDp1aputRkREREQt9aokKDMzE6WlpQgNDdUeCwsLw6lTp/TKrlmzRuexSqXCgAEDTB4jERER9Q+9KgkqLCwEADg7O2uPubi4oKioqN3rKisrcfLkSdx7771tlpFIJDo/SUlJRomZiIiITC8pKUnnc9wYet2YIAB6leto5Pdjjz2Gv/3tbwgICGizDGeHERER9V1JSUk6DRjGSIR6VUuQp6cnAKC8vFx7rLy8XHu8Ne+99x6srKzw0EMPmTw+IiIi6j96VRI0ZMgQuLm5IS0tTXssJSUFN910U6vlv/jiCxw4cAD/+te/AAAZGRk9EicRERH1fb0qCZJKpVi+fDk2b94MoCmpOXPmDBYtWoTU1FTEx8dDrVYDAJKTk/HGG2/grbfeQk1NDaqrq/Hss8+aM3wiIiLqQ3pVEgT8vk5QTEwMFi5ciF27dsHb2xsVFRW4ePEiGhoaAAALFy7EwYMH4eHhAUdHRzg6OiInJ8fM0RMREVFfwW0ziIiIqM/pl9tmEBEREfUEJkFERERkkZgEERERkUViEkREREQWiUkQERERWSQmQURERGSRmAQRERGRRWISRERERBaJSRARERFZJCZBREREZJGYBBEREZFFYhJEREREFolJEBEREVkkJkFERERkkZgEERERkUViEkREREQWiUkQERERWSQmQURERGSRmAQRERGRRWISRERERBaJSRARERFZJCZBREREZJGYBBEREZFFYhJEREREFolJEBEREVkkJkFERERkkZgEERERkUViEkREREQWiUkQERERWSQmQURERGSRmAQRERGRRWISRERERBaJSRARERFZJCZBREREZJGYBBEREZFFYhJEREREFolJEBEREVkkJkFERERkkZgEERERkUViEkREREQWiUkQERERWaRemQSpVCokJiYiJiYGY8aMwZ49e9os+8orryA6OhrR0dF46aWXejDKviEpKcncIZgN626ZWHfLxLpTV0iEEMLcQbS0fv165OfnY9u2bUhPT0dMTAxSU1Ph5eWlU+67777D6tWrcebMGQBAREQEXnnlFcyZM0dbRiKRAAB6YTV7hEQiYd0tEOvOulsa1t3y6m6Mz/de1xKk0WiwadMmLFu2DAAQEhKCyMhI7NixQ6/su+++i7vvvhsKhQIKhQKLFi3CO++809MhExERUR/U65KgzMxMlJaWIjQ0VHssLCwMp06d0it78uRJg8oRERERtSQ3dwAtFRYWAgCcnZ21x1xcXJCSktJq2ZblioqKWn3e5mYzS8S6WybW3TKx7pbJkuveHb2uJahZyxvaVp8fbzwRERF1Ra9rCfL09AQAlJeXt/rvlmXLy8u1j8vLyzFgwACdMpY4WIyIiIg61utagoYMGQI3NzekpaVpj6WkpOCmm27SK3vTTTcZVI6IiIiopV6XBEmlUixfvhybN28GAGRkZODMmTNYtGgRUlNTER8fD7VaDQBYuXIldu3aBZVKhSNHjuCZZ55BcHBwu8//0UcfITo6GmPGjMG6dev6RUvRyZMnERwcjK1bt7ZZ5uDBgwgNDUVcXJz259NPP+25II2soaEBr7/+OuLi4hAbG4vx48dj//79bZbvT/e9M3Xvb/cdAN58803Ex8dj2rRpGDVqFN588802y/an+w4YXvf+eN+bZWRkwMrKCgcPHmyzTH9dP66jum/duhURERE69/3IkSM9G6SRJSYm6tQnLi4O1dXVbZbv9Hte9EJKpVIsXbpUjBs3TkRHR4vvv/9eCCHEsWPHxMCBA4VSqdSWffnll8XgwYOFq6urGDRokNiyZUubz3v+/Hnh5eUliouLhVqtFnFxceKtt94ydXVM6vPPPxcLFy4UUVFR7db9wIED7Z7va7KyskRgYKAoLy8XQgixZ88eYW9vL65evapXtr/d987Uvb/ddyGECA8PF4WFhUIIIdLT04VMJhMnTpzQK9ff7rsQhte9P973ZnfddZews7MTBw4caPX8t99+K4YNGyaUSqVQKpVi2LBhYvfu3T0bpIl0VPctW7a0ea6vWrp0qcFlu/Ke73UtQQCgUCiwdetWHD9+HKdOncKMGTMAADExMbh27RoUCoW27Lp163Do0CGUlZV12Ar0wQcfYPbs2fDw8IBUKsV9993X59cVuummm7Bz5044OjqaO5Qe5ejoiKefflo7OzAhIQEKhQJHjx7VK9vf7ntn6t4fbd++XTtGcOjQoXB1dUV2drZeuf523wHD695fnTx5Eg4ODnpjP2/UX9ePM6Tulq4r7/lemQR1lq+vr0HlWltX6MKFC1AqlaYKzeQMrTsAfPnll5g6dSomT56MZ599Fo2NjSaMzLTc3d1xzz33aB8LIVBfX9/qH4j+dt87U3egf913AAgPD9f++7PPPoODg4P2i9KN+tt9BwyvO9D/7jsAPPXUU3jqqafaLdNf148zpO5AUxIYGxuLuLg4bNy4sQciM72VK1di8uTJmDVrFg4cONBmua6853vd7DBTam1dISEESkpK4OfnZ8bITM/Z2Rnjx4/HunXroFKpMHfuXFy/fh2vvPKKuUMzikOHDiEgIABTpkzRO9ff73t7de+v9/38+fO48847UVNTg48//hguLi56ZfrrfTek7v3xvn/77bcYMWJEh1/8OrN+XF9haN29vLwwY8YMJCYmori4GHFxcVCr1Xj44Yd7KFLjGz58OKZNm4abbroJJ0+exNSpU3HkyBFERETole3Se77LHXW9UGxsbLv94MHBweLdd9/VPr506ZIAIHJycnogOtPqqO4tffPNN8LOzk5oNBrTBdVDlEqlGD9+vDh16lSr5/vzfe+o7i31p/suhBA///yz8PLyanVcTH++70K0X/eW+vp9V6vVIjY2VpSWlgohhAgICGhz7ItcLteOIxVCiH379gmpVNoTYZpEZ+re0saNG0VoaKgJo+t5d9xxh1i5cmWr57rynu8X3WGGam1dIYlEYpF9rP7+/qitrUVxcbG5Q+kWIQQeeOABrFmzBtHR0a2W6a/33ZC6t9Rf7nuzyMhIzJkzB2+88Ybeuf5635u1V/eW+vp937lzJ26++Wa4ubl1WNaQ9eP6ks7UvSV/f3/k5OSYICrzaa9OXXnPW1QS1Nq6QiNGjICtra0Zo+oZb7zxBlQqlfZxYWEhrK2t4e7ubsaoum/dunUYO3Ysbr/9dtTV1SE3N1evTH+974bUvb/d99LSUnz22Wc6x+zt7VFTU6NXtr/d987Uvb/d98OHD2P37t3aKdIFBQVYvXo15s+fr1e2v60f15m6v/jiizqPCwsLMXDgwJ4K1SQ6U6cuveeN3FJlVi27hIqLi8WkSZO0zYjnz58X3t7eoqSkRKjVajFt2rQ+P2W2WUd1X7p0qdi2bZsQQojGxkaxYMECkZiYaI5QjeaFF14QDz/8sKiqqhJVVVXiwoUL4qmnnrKI+25o3fvbfc/KyhKjR48WNTU1QgghioqKhJ+fn3jzzTf7/X3vTN37231v6cYuoZSUFDFt2jTR2NgohGiaIj98+HDtFPnQ0NB+M0VeiPbrHhcXJw4ePCiEEKKmpkaMHTtWJCUlmStUo/Dy8tIuC5GZmSmcnJy0dTTGe75fJEGnTp0SsbGxwtnZWQwbNkzceuutQgghcnNzhYeHh7hy5Yq27I4dO0RUVJQYM2aMWLt2bZ/tI29maN2PHDkiZs2aJeLi4sSYMWPE8uXLRUVFhTlD75a0tDQBQO/nqaee6vf3vTN172/3XalUiieffFLExMSIKVOmiPDwcPHUU08JtVrd7+97Z+re3+57sxMnTojY2FhhY2MjRo8eLf7+97+3uX5cVFSUiIqKEi+++KIZIzYeQ+r+9ddfi2nTpmnv+2OPPSbq6urMHHn3vPTSS2LChAliypQpIjo6Wnz44Yfac8Z4z0uE6ONLqBIRERF1gUWNCSIiIiJqxiSIiIiILBKTICIiIrJITIKIiIjIIjEJIiIiIovEJIiIiIgsEpMgIiIiskhMgoiIiMgiMQkiIiIii8QkiKgPmzt3LmQyGUJDQ/Hee+9pj99xxx1YsmSJ9vFXX32FsLAwhIWF4dChQyaL5x//+AcCAgJw//33m+w1DHHXXXfB2dkZO3bsMPiaSZMmISgoCHFxcUhJSenS69bW1mLgwIEoKCjo0vXG8sILLyAuLg42Njad+h0QWRomQUR92O7duxEREYF58+bhgQceAABoNBr88MMP2L9/v7bcLbfcgttuuw07d+5EbGysyeLZsGED7rnnHpM9f2v27duH4OBgnWP/+c9/MGrUqE4/17Jly3Dw4EGEhYV1KRY7OzucP38e3t7eXbreWNavX4+DBw9iwIABZo2DqLdjEkTUxyUkJGDfvn3ax7/88gvi4uKQl5eHCxcuaI+fO3cOo0ePNkeIFsXd3d3cIRCRgZgEEfVxCQkJOHv2LIqLiwE0tYysWLECQUFB2Lt3LwCgoqICLi4ukEgkAIDjx49j5syZSEhIQGxsLJ544gkIIVBdXY2oqChIJBLExsZCqVTi/PnzCA8Px+DBg3Hs2DFoNBo899xzGDduHOLi4jB//nxcu3atzfjKy8uRmJiICRMmYNKkSVi9ejXq6uoAAKtXr4aXlxfWrl2LJUuWYPz48Rg/fjxyc3O111+7dg0zZ85EdHQ0Zs2ahWeeeQa2traIj4/HsWPHsG7dOly7dg1xcWW3DuIAAAfvSURBVHGIi4tDWVmZ9tqsrCzccccdiIqKws0334yKigqDf69nzpzBuHHjIJfLsWPHDsycORMhISH46KOPsG/fPixYsAAhISF4+eWXtdfMnTtXrxvu6NGjiIuLQ2xsLMaPH481a9agqqoKW7duxbBhwxAfH48nn3wS8fHxsLa2xpEjR1BXV4dHHnkEEyZMwOTJk3HnnXeisLAQAHSu27BhA+Li4jB8+HCdRJiIDGT8je+JqCepVCpha2srdu3aJYQQYt68eUKpVIoHHnhAzJ49WwghxBdffCE2b96svWbPnj3i559/1j5euHCh2LlzpxBCiJqaGuHk5CT27t2rPf/YY4+JI0eOCCGEeOmll0RoaKiorKwUQgjx6quvinHjxmnLbtiwQSxbtkz7eM6cOWLx4sVCo9EItVotbr31VvHoo49qzy9atEiMGDFCVFVVCSGEWLp0qfjzn/+sPT9p0iSxcuVKIYQQjY2NYvbs2WLIkCHa83v37tV53GzixIli+vTpor6+Xmg0GhEXFydee+21Nn+PEydOFM8884zOsYyMDCGRSMSHH34ohBDi6NGjwsnJSbz++utCCCFycnKEtbW1KC4u1nme7du3a8/b2dmJgwcPCiGEUCqVIjw8XJw/f14IIcT7778vbG1txenTp4UQQrz55pvi3Llz4uGHHxYzZ84UDQ0NQggh1q9fL2JiYoRGo9FeZ2dnJ3799VchhBBbtmwRo0eP1qvToEGDtLEQkT62BBH1cTY2Npg0aRL27t2rbWFRKBSYPn06kpOT0dDQgP3792P69Onaa0aOHIm33noLEyZMQFxcHH788UccP34cQNO4lkWLFmkHWjc0NODkyZOYOHEiAOCdd95BYmIiHB0dAQBLly7FiRMncP78eb3Yrl69iv/9739YtWoVJBIJpFIpFi9ejE2bNumUu/nmm+Hg4AAAGDNmDDIyMgAA2dnZOHLkiHagtUwm69SYo/nz58PKygoSiQTR0dHa5+0MIQRuvfVWAEBkZCQqKysRHx8PAPD394eLiwuysrJavfbDDz9EYGCgdhyWQqHAtm3b4Ofnpy0zfPhwREVFAQBWrVqFsLAwvPfee1i2bBnkcjkA4E9/+hOOHz+On3/+WXtdaGgoRowYAUD3d0ZEhpObOwAi6r6EhAS89dZbOHLkiDZZiY+PR21tLY4dO4YrV67ofPDefffdGDBgAJKTkyGXy/Hkk0/qzGhavnw5YmJiUFxcjMOHD2PevHnac7m5udi2bRu+/fZb7bGAgAAUFhbqDUZu7tZas2YNrK2tAQAqlQqOjo6orq7WJj6urq7aaxQKhTaZa+5mu3GAr4eHh8G/FxcXF53nLS8vN/jaZjKZTBtnc1Jy4/NaWVlp420pNzcXXl5eOsciIiJ0Ht9YdwAoLCxEXV0dPD09tceanyM3NxfR0dF61934OyMiwzEJIuoHEhIS8Oijj2Ljxo144oknAABubm6IjIzE5s2bdRIgADh27Bjefvtt7Yd6fX29zvnIyEiEh4djy5YtOHbsGD744APtOX9/f6xatQp/+tOftMfKy8thb2+vF5e/vz8AYNOmTTozroqKirSJRXsGDRqkLd/8XCUlJR1e11v4+/vjxx9/1Dl2+fJlODs7t5nMeXl5wcbGRjsGCID2382/AyIyDnaHEfUDo0eP1rbsREZGao8nJCRg+/btSEhI0Ck/bNgwJCcnAwDq6uq0A6hv9MADD+Bf//oXnJyc4Obmpj2+cuVKbN++HbW1tQCA/Px8jB07FkqlUu85fH19MWfOHJ01jPbv348777zToHoFBgZi0qRJeP/99wEAarUau3bt0inj5OSE6upqAMDLL7+ML7/80qDn7glLlixBdna2dm2mmpoa3HrrrdrfXWtkMhnuv/9+vP/++2hsbAQAbNy4EWPHjtV2mxGRkZh7UBIRGcfChQvF7bffrnNs//79Qi6Xi4qKCp3jp06dEhEREWLMmDHijjvuEAsWLBDe3t7i73//u7ZMVVWVcHR0FIcOHdK5trGxUTz//PNi3LhxYurUqWLq1Kni8OHDQgghnn32WeHv7y+8vb3FI488IoQQoqysTNx3331i/PjxYurUqWL+/Pni2rVrQoimQdSenp4iMDBQvPvuu+Lzzz8XISEhwsXFRTz44INCCCGuXr0qbr75ZhEZGSnmzJkjXnzxRREcHKyNp6GhQSQkJIiYmBgxadIkUVJSIpYtWyacnJxEaGio+O9//ys2btyojevGOt6o5cDojIwMMXbsWAFAxMbGitLSUjFlyhQBQMTExIicnBwxZ84cYW1tLSIiIsTRo0fFnDlztK/70UcfCSGE+PHHH8WUKVPElClTxIQJE8Rnn30mhBDik08+0dY1NjZWFBUVaV9bqVSKtWvXaut0++23i7y8PL3rHnzwQZGZmakTZ11dnfZ5ODCaqH0SIYQwcx5GRNSm0tJSnbV3du7cibfffhuHDx826utMmjQJM2fOxJNPPmnU5zUnX19fvPDCC1i8eLG5QyHqldgdRkS92mOPPYYDBw4AaJqp9uGHH2Lu3LlGfx0fHx989tln3do2o7do3jZDoVAYNPaKyFKxJYiIerVPP/0U//znP+Hg4ICamhrExsbiueeeg5WVlblDI6I+jkkQERERWSR2hxEREZFFYhJEREREFolJEBEREVkkJkFERERkkZgEERERkUViEkREREQW6f8DQF815qoH4w0AAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "execution_count": 11, "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": 12, "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.]),\n", " 'logg': array([3. , 3.5, 4. , 4.5, 5. , 5.5]),\n", " 'feh': array([-0.6, -0.3, -0. , 0.3])}" ] }, "execution_count": 12, "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": 13, "metadata": {}, "outputs": [], "source": [ "model_param = {'teff':1500., 'logg':4., 'feh':0., 'radius':1., 'distance':20.}" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [], "source": [ "read_model = species.ReadModel(model='drift-phoenix',\n", " wavel_range=(1., 5.))" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [], "source": [ "modelbox = read_model.get_data(model_param)" ] }, { "cell_type": "code", "execution_count": 16, "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": 17, "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/+yKZPn85Gjx7N9u7dy8LDw1lBQYHamCMiItjKlStZ48aNmY2NDcvOzla734ULF9gnn3zCALChQ4eyqKgoxXPh4eHsk08+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/f388ePAAQ4cOhVQqxddffw0A2L59OwYPHoy6deu+tYx33nkHe/bsgaWlpcrSI1wICwtDixYt9CqjQYMG2LhxI2JiYnDgwIEyvzXVrl0bnTt3xrlz58AYw7x587Bw4UIcO3ZMrxi4FBsbC0dHR77DUIuSIPI28fHxOn+hIcratWuHwMBAvsMwmsjISDg7Oxv0GF27dsXDhw8Negyu6Z0E1apVC7NmzcLs2bMxbdo0AMDDhw8VyVFF5eLigiVLlmDNmjXw9vZGzZo1AQCTJk3SqhaoJEtLS4PUBBUUFMDS0lLvcjp27IiFCxeiRo0ab9130aJFWLt2Lb755hs4OTnh+PHjOHjwoN4xcOXevXvo1asX32GoRUkQKYtcLud04cjKrl27dnj+/DnfYRhNTEwMGjZsaNBjDBkypNz3P77pnQTt2rULdnZ2mDRpEqZNm4aEhAT8+eef+O6777iIT7CsrKwwfPhwuLu7612WoZIgPtja2mLv3r0YPXo0vvvuO9ja2hqklktXz58/x7vvvst3GGpREkTKEhMTg8aNG/MdRoXRtm1bBAUF8R2G0SQkJBi8FrxHjx7w8/Mz6DG4pnfHaHt7e6xevVrx2NHREf/884++xZqE06dPc1KOoZIgvr4xNmzYUOkbh4WFBSQSCSe1UvoKCAjA3Llz+Q5DLWtra2RmZvIdBhGo2NhYODk58R1GhdGgQQPExcXxHYbRSKVSVKmi9y2/TNWqVYNcLlcaOCR0ekd59uxZuLu7QyKRICgoCM2bN4ezszPu3LnDRXyCxWWCYYgkSCaTCeZN2KxZM4SHh/MdBgAgKysLNjY2fIehlpmZmeA6kRPhoP5A3KJmRcMQ0t97beh9l9y4cSPmzJkDS0tL/PTTT/jggw+wbNkyeHh4cBFfpVC9enXOm4zEYjHs7Ow4LVNXzZs3x4sXL/gOAwUFBahWrRrfYRCiE0qCuGdmZmZyo5l0IZPJYG5ubpRjmVpfK72TIMYYOnfujKysLPj6+mLjxo2YOnUqqlatykV8lYIhaoLS09MFM0KvRYsWgkiCQkNDFaMXhYq+nRJNKAniXv369REfH893GAaXnp4Oe3t7oxyrffv2CAgIMMqxuKB3EpSRkYGCggLs27cPgwcPRvXq1QHAaFlnRVAZkqCXL1/yHQaCgoIMOlkYIYZESRD3GjdujOjoaL7DMLiUlBTUrl3bKMcytZogvXtJjRs3DvXq1UNBQQFu3LiBzMxMLFmyxOBD8fjEdX8bS0tLpKWlcVYeAKSlpQkmCWrUqBGioqL4DgOxsbE00RwxWcnJyahTpw7fYVQojRs3RlRUFN577z2+QzGo1NRUo713atWqhfT0dKMciwt6J0GLFi3Chx9+iJo1a8LJyQl5eXkYO3asQWem5Ft+fj4sLCw4K8/CwgIFBQWclQcIqyaoSpUqkMvlfIeBuLg4DBo0iO8wCNGJXC6nGnaONW7cGPfv3+c7DIMzZk0QUNSsbypLknBSndGqVSvF0M3q1aujX79+WLduHRdFCxLXw72rVauG/Px8zsoDhJUECUVcXJzgmxNEIpEgEkZCKgOh1FIbmjFrggCgbt26SE5ONtrx9KFTTdDx48cxZswYiEQijZMFnjt3rsImQlzXBFWrVs0gNUEdO3bktExTJ6QRc5rUqFEDubm5sLKy4jsUIiBSqZRqgQygYcOGiImJ4TsMg0tJSUHz5s2NdrziwTDaLB/FN51qgnbt2oXs7GwAwNWrV9GkSROVHyFMjGco1BymG77nwDGF6tmaNWsiJyeH7zCIwCQnJ5vEDcXUWFhYoLCwkO8wDM7YNUHNmzcXxGAYbehUE1RypuTJkydj2bJlKvvwfcMzJFOpCRJSElR8c+ezhkPoCRBQdJ2ys7NRr149vkMhAkIjw4g+jN0nqHnz5pytqGBoevcJ8vT0VLu9In9r4XrSvcqQBPE9YsBUknIrKyuqCSIqKAkynBo1alT4z1xqaqpRkyChTIuiDb1Hh8nlcnh7eyM0NFSpc6+Xlxdmz56tb/GCZIgkiOuO0VzXVunLzs4OGRkZaNSoES/HF/JyGSVREkTUoSTIcBo1aoTo6Gi0bt2a71AMRiqVGnUCY76/9JaH3knQhAkTEBgYiI4dOypdZF0n/5NIJJg5cyaCg4MhlUrxyy+/qB3W7OXlhU2bNil1dF21ahV69+6t03HLo7CwUPA1QUJjY2PD6+KgqampRpsxVR9WVlbIysriOwwiMAkJCTTHlYFUhiSIaKZ3EhQSEoKnT5+qTB6oqZnsbTw9PcEYw927dxEaGgpXV1cEBQWp7SOxadMmuLm56XQcfXBdE2SIjtFCY2try2sSlJaWZtTqYF1ZWVkpBh0QUiwxMbFCdzHgU3ESRLhlY2ODjIwM2Nra8h1KmfTuE9SuXTu12zt37lzusuRyOXbu3Ilp06YBAFq2bInOnTtj7969esXINVPoEyQ0xR8IvqSlpZlETZC1tTXVBBEViYmJ1FneQCp6EsRXf0hT6Rekd03Qr7/+ii+//BKdOnVSappas2YNRo0aVa6yXr16hdTUVKVqybZt2+LBgwdq99+xYweWLVsGkUiE8ePHG60PktCTICHOKWJjY4OkpCTejm9KzWF8XiciTLm5uahZsybfYVRIFT0JyszM5KU2pniYfJcuXYx+7PLQOwlat24dTp8+jbCwMFSp8qa4hISEcpeVmJgIAEr/YXZ2dggMDFTZt169ehg0aBCmTJmC5ORkuLm5QSaTYc6cOTqcRflwnQRVrVqV047RYrFYUCPDgKIkiM+V5MVisWJWcyGj5jBCjKtevXqKe09FxNcXwObNm+PWrVtGP2556d0cdvHiRURFReHGjRu4evWq4mf8+PE6l1l6Phd11XlDhw7F1KlTIRKJULduXcyZMwdbt24ts8ySP7r2WQK4T4LMzc05XSpBaMPjAf47RpvCbNEANYcRYmxmZmYmM4WGLvjqD1k8azSXPD09le7jXNA7CWrdurXa2aGnTp1a7rKKO/6JxWLFNrFYrFWHwMaNGyMyMlLj84wxpR8hJUEAtxP5CTEJsrW15bVPkKkkQVQTREorKCgw6vBmUrHwVRPk6OioU4tQWYoHThX/cEHvJKhfv35wd3fHqVOncP36dcXPrFmzyl1W8+bNYW9vj5CQEMW2wMBAuLi4qOy7du1apceJiYlo0KBB+U9AB4ZIgriUlpYmuCSI75ogISaG6lhbW1MSRJQkJSXRyDADq8hTU/BVE1S8krzQ6d0naPHixXB0dMTVq1eVtuvSxmpmZobp06dj165d6NOnD8LCwvDkyRPs27cPQUFB+Oabb3DhwgWYm5vj7Nmz6NGjB/r164fc3Fzs2LEDkyZN0vd0tCL0JEiIN3y+kyBTqgmqqH+MiW5oeLzhOTk5ISYmBm3atOE7FM6lpqYadfHUkiwsLAQ3cW9peidBvXr1wqVLl1S2jx49WqfyPD09MXPmTLi6ukIqleLAgQNwdHREREQEgoODUVhYCHNzc3h4eGDFihWQy+XIzs7G+++/jx9++EHf09GKKSRBTZo04TsMJRYWFjpPoMmFnJwc1KhRg7fja6t4FXlCiiUlJdHweAOryElQWloaunfvzsuxnZ2dERERgVatWvFyfG3onQSpS4AA4OTJkzqVZ2lpCS8vL5Xtrq6uiI2NVTweMWIERowYodMx9CX0JCgtLU1wwxKFsHipEGJ4G1OIkRgX1QQZXqNGjRATE8N3GAbB5xxpzZo1w6tXrwSdBOndJ+j+/ftYsWIFCgoKEBsbi4EDB8LNzU2pX09FYwpJkCnMiUPUM4V2dGI80dHRvK25V1kU1wRVRHzOkVY8V5CQ6Z0ELV++HO+++y6qVq2KJUuWoGbNmujXr1+FXTwVEP5oDUqCCKk4IiIi0LRpU77DqNAqchLEZ3/I4pogIdO7OaygoABjxoxBfn4+Tp8+jRcvXsDOzs4oC5nyReg1Qenp6SbRCdhYqGaFmDJaMsPwHB0dER8fz3cYBsEYU1nb01iK+wQJmd5Xpng478mTJ9G7d2/Fzbd69er6Fi1YQk+C5HK50uzdQmFmZgaZTGb04+bl5ZlEp2hC1GGMUV8xAzM3N+flb1NFVzw6TMj0vlP27t0bbdu2RWxsLE6cOIG8vDxs3boVVlZWXMQnSIWFhYJOgoTKxsYGWVlZRq+lMpXh8cVEIhHkcjlv396IcMhkMnofGEmVKlVQWFgo6K4OpkrIiTwna4cNHDgQtWrVQvfu3ZGZmYk6deqoTGZYkQi9JkiozT/FcwVRElS2mjVrIjc3t0J/kSDaiY2NRcOGDfkOo1Jo0qQJoqKieJtTxxCkUinvSXS9evUEPc0DJ1dn8ODBinkIbGxsMHnyZFy5coWLogXJEEkQV01FQk2AAP4mTBTi5JFloaUzSDHqFG08hljrim/Jycm8T68g9BFietcEyWQyeHt7IywsTKntz8vLCzNmzNC3eEEyRBJUrVo1xUSQ+sjOzoaNjQ1HUXGLryTI1GqCaOkMUiwiIgLOzs58h1EpFCdBgwcP5jsUzsTHx6N+/fq8xlA8Qqxnz568xqGJ3knQZ599Bn9/f3Tq1ElpIVU+Zwc2NEMlQQUFBWoXoy0PIQ+P52sRVVMbLUdLZ5BiERERFeqmLGQtWrTAxYsX+Q6DUwkJCYJIgs6dO8drDGXROwkKCgpCQECASg3GypUr9S1asAyVBHHRi17ISRBfNUHJycno0KGD0Y+rK2oOI8XCw8OpJshInJ2dER4ezncYnBJCTZDQm8P07hP0zjvvqG3CGTRokL5FC5ZUKtW72ao0CwsLFBQU6F1OcnIy6tSpw0FE3OMrCUpKSoKDg4PRj6sra2trqgkiAGgFeWMq7pJQkcTHx8PR0ZHXGGrXro2UlBReYyiL3jVBffv2xeTJk/HRRx8pdT6dN28eHj16pG/xgsX1cD9LS0vk5eXpXU5CQgLvb3pNbGxseOl4mJycbFJJkI2NDS/NhkSYhDq0uCIyMzODVCoV5DxruhBCc5hIJFK0dAhxNXm9/6cXL14MR0dHXL9+XWl7YmKivkULliFGYHH17T8hIQGurq4cRMQ9vm7uKSkpgq0dU6dOnToICwvjOwzCM0PUOJOytWrVCiEhIWjXrh3foXBCKF+K27dvj4CAAHTt2pXvUFTo3RzWq1cvhIeHq/xU5OYwQ+CqqUgob3p1bG1teWkOk0qlJjUBmoODA5KTk/kOg/AsNjYWTk5OfIdRqXTq1AlPnjzhOwzO5OXl6T3YhgtdunQRbMuQ3knQpUuX1G7fvHmzvkVXKlzWBAk1CeKrT5CpoSSIAEWdommOIOMS8s26vORyuWCaUoV8XTlr+ExNTUVOTo7i8YQJE3D79m2uihcUQ7yxuEqCMjMzYW1tzUFE3OMjCRLy5JGaUBJEAOD27duCbD6oyIqbwyoCIdUkNmnSRLALqeqdBD1+/Bjjx4/Hq1evlG44QslAucYYM8iN1cbGhrNVjIV67S0sLDjp/F0eqampJtUfCChafNjY14kIj5+fH+bNm8d3GJWKubk55HK5oNe60lZISAhatWrFdxgAiu5JNWrUQEZGBmxtbfkOR4nezWGLFy/GhQsX0LdvX8jlckgkEhw5cgQ//fQTF/EJjqHWDbO2tuaklkTINR98/FEx1XlWhPz/SAxPKpVCIpGgRo0afIdS6bRp0wZPnz7lOwy9BQYGonXr1nyHofDBBx8IcjJKvZMgxhiaNm2q+KNdrVo1fPzxx3j8+LHewQlRTk6OQRa25KI5TCwWCy7LLs3MzAxyudxoxzPVZQeqVq3KybxRxDTdu3cPPXr04DuMSunzzz/H7t27+Q5Db9evXxfUUhXDhw/HmTNn+A5Dhd5JkFQqLSrIzAwPHz4EAISGhsLf31/fogUpJSUFtWvX5rxcOzs7pKen61WGKXSkdHR0NOr0CaaaBDk7OyMyMpLvMAhPzp8/T8tl8KRTp04ICAgw6drYwsJC5ObmCupLcaNGjRAbG2vUL8Ha0DsJcnFxwYYNGzBjxgy4urqiUaNGaN++PcaNG8dFfIJjqFV5GzZsiNjYWL3KMIUkqHHjxoiKijLa8Z4/fy6YdvHyaNGiBUJDQ/kOg/BAJpPh5s2b6N69O9+hVEoikQg9evTAjRs3+A5FZ3fu3MF7773HdxgqXF1dce3aNb7DUKJ3x+hff/1V8buTkxPu3LmD1q1bY8SIEfoWLUiJiYkGSYIsLS31XnT22bNngv/22KRJE7x69cpoVf1CWDtHF927d4ePjw+GDx/OdyjEyM6ePYshQ4bQRIk8cnd3x7fffos+fcAZBHgAACAASURBVPqYZAfp8+fPY/To0XyHoWL69OlYuHAh+vfvz3coCnrXBDVr1gyzZ88GUDRx4sKFC/VOgCQSCaZMmQJXV1d069YNFy5c0Ljvhg0b0LVrV3Tt2hXr1q3T67jaePnyJZo1a2aQsvWtfvX391dZKNTT01OvMrnWrVs3+Pn5GeVY8+fPR6NGjYxyLK517txZr0nbhPb/bkymfO4PHjzAyJEj4e7urtPrTfnc9cXluTdr1gydOnWCj48PZ2UaUulzf/ToEbp06cJPMGVwcnKCmZmZUVsD3orpqWPHjvoWoeL7779nX3zxBWOMsZCQEFarVi2WkJCgst/Zs2dZq1atWF5eHsvLy2OtWrVip0+fVtoHAOPgNBWmTZvG4uLiOCuvpKlTp7KYmBidXiuVStngwYNVtnN57lyQyWRq4zQEAOzQoUNGOZYhDB8+nGVnZ+v0WqH9vxuTKZ67TCZj//zzDwPAtm3bpnM5pnjuXOH63LOysljbtm3ZzZs3OS3XEEqe+6tXr9ikSZN4jKZsvr6+7NNPP2VyuVzvsri4v+tdE9SqVSu1K+8uWLBAp/Lkcjl27tyJadOmAQBatmyJzp07Y+/evSr77tixA5999hksLS1haWmJiRMnYvv27TodVxsymQzh4eEGm5G5b9++uHz5sk6vvXz5sqBGAmhiZmYGR0dHgy+kKhaLAUCQVcLamjx5MjZt2sR3GMSA4uPjsXv3bowcORJTp06Fn58fZs6cyXdYBICVlRWuX7+O7du346uvvkJqairfIb3Vs2fP4Obmhh9++IHvUDTq168f2rZti/nz5yMuLo7vcCB6nU3pbNWqVbh48aLKKvJr1qxBYGBguct78eIF3nnnHaW+N3PmzEFKSgoOHDigtK+TkxN+++03jB8/HgBw7NgxfPPNN0qTDha35+p5mgAAHx8fPHnyBEuWLNG7LHWys7MxfPhw7N69u1wjmqKjo/HFF1/gyJEjKhMDikQiwY1yCAwMxKJFi7B//37ORy9IJBL4+Pjgzz//hK+vr+DOvTwYY5g6dSqaNWuG77//vlwrMAvx/91YhHbujDGkpqYiIiICsbGxePToEZ48eYJTp06hT58+aNSoEcaMGYOxY8fqfSyhnbsxGfLcv/zySzx79gxNmzZF586dkZOTg/79+6Nbt268ztCfkZGB6tWrw8LCAjNnzoS/vz8OHDiAxo0b8xaTNhhjuHz5Mnbu3Ilq1aph2rRpeOedd+Do6AgzM+3rZri4v+udBNWqVQudOnVS2f706VOkpaWVu7xbt26hd+/ekEgkij/6P/30E27fvq1SS1K1alWcPn1a0Rn48uXLGDRoEGQymWIfU+zURgghhBDt6JPG6DQ6bMWKFQCK1gcbPny42qYqfat0Sycvmk6SkhxCCCGE6EKnJOjEiRPYvHkzateurTYBAqBz35ziJjCxWKz299L7Fvf9KN7PwcFBbblURVz50LnTuVc2dO507pUJF5UgOnWMtrOzQ9++fWFvb69xn7CwMJ0Cat68Oezt7ZVW8g0MDISLi4vKvi4uLlrtRwghhBBSmk5JkDbZ14wZM3QpGmZmZpg+fTp27doFoCiZevLkCSZOnIigoCC8//77ij4/M2fOxIEDByCRSCCRSLB//34aWUEIIYQQrejUHObr62vQ2Uw9PT0xc+ZMuLq6QiqV4sCBA3B0dERERASCg4NRWFgIc3NzDBkyBM+fP0evXr0AFM3ySTPsEkIIIUQbOiVBHTt2LHP+EsYY5s+fr3NQlpaW8PLyUtnu6uqqsr6Wh4cHPDw8dD5WRbds2TK+Q+ANnXvlROdeOdG5E13oNER+wIABuHLlSpn77N69G5MnT9Y5MK5wOU8QIYQQQoSBi/u73jNGayKEBIgQQgghRBOdkqD4+HgMGDAA9+7d4zoeQgghhBCj0CkJCgoKwpUrV9CjRw+u4yGEcGzHjh1wc3ND+/bt4eTkBDc3N7z//vuQSqUGOd6dO3fQsWNHHD16VGm7q6sr3Nzc4ObmhnXr1gEA9u/fj++++w6zZs3C/fv3Fftq2l5SWfvcvn0bnTp1wvjx45GSkoL9+/fDxcUFp06d0hh3REQEBg4cqOtpc+7evXuoWrUqUlJSVJ7T57oV8/LygpubGyex+vr6okuXLvi///s/pWvv6ekJT09PDBgw4K1ljBw5Uum1YrEYCQkJGD16NLp3744HDx4gODgYCxcu1Ctufd5b2l5fLq8tULR0VOfOndG7d2/06dMH6enpepWn6TNaTJ/z5zpWg9Nr+VUTAI5XkSfEVG3atImtXbvW4Mc5cOAAmzx5Mjty5IjS9mXLlik9FovFrGvXrowxxjIzM1mXLl3K3K7Na0sqjiE+Pp59/fXXLDMzs8y4w8PD2fvvv6/dSRpYbm4u++abb1iTJk1YcnKy0nP6XLfS+vXrx0m8V69eZZMnT1Y8Lv3/7+PjU+brX7x4webOnav2tf/88w/7+uuvFY/Dw8N1jluf91Z5ry9X15Yxxvr06cPCw8MZY4xNmzaN7d27V6/yNH1GGdP//LmOtSxc3N91Gh1GCDE9/v7+isWGy0sikWD06NEq24cPH445c+Yobfv0009x7tw5tcdfu3YtsrKyMHXqVLx48QJt2rQBAFhbW0MqlSI+Ph7Pnj1Tu71+/fqKsu7du/fWfQAgKioKa9aswbp161C9enWtzpUxhpkzZ+L58+eQy+XYtGkTunfvDplMhkmTJiE2Nhbvvfcejh49ihcvXihed/jwYWzcuBFWVlZwcHCAu7u7omapPNcPKFqAetGiRfjvv/9UntN07tpcN7lcjmnTpqF27dpo1qyZ0jkvXLgQ9vb2iImJwcyZM3HkyBHs3LkT9+/fx8SJEzFp0iR89dVXWl3D0ucyc+ZMtcfo2LEjAODkyZMYNWpUucsupu311eZ9o+v1NeS1jY6OhrOzM3JzcxEQEKCyb3nfX5o+o/qcv7axCg0lQYRUEv7+/li1apVOr7W0tNT4R1NbP/74I7p164b4+HgMGTIE33//PaysrBTPW1tbIyUlBSkpKWq3l/xjq80+ALB8+XKsXr1a6wQIALy9vVFYWIibN2/i1atX+PTTT3H//n14e3vDxsYG+/fvx61bt3Dw4EHFa16+fIlt27bh2rVrAIAGDRrg999/Vzxfnut37do1ODs7a1wJXNO5a3NNvL29UbVqVaxfvx6vXr3C4cOHAQA+Pj4Qi8XYsGEDoqKiMHXqVFy+fBmRkZHYt28fxo8fX+6bmZeXF27evIkHDx5g5syZGo8BADdu3MDcuXNVXgsUrQTQsmXLMo+l7fXV5731ttca6tpGRUVBLBajb9++ii8S3bt31+n8taHP+0ubWIWGkiBCeHLkQTRi0vP0LsepVnWM69aozH3kcjmSkpJUkoRiISEhOHXqFLp374779+/D1dUVffr0UTwvkUgwYsQIldeNHDlS6eZVlm7dugEA6tevj4KCApibmyM7O1vxfFZWFurUqYM6deqo3V6SNvsAwNq1a7FlyxY4ODhg3LhxWsUZEhKCnj17AgCaNWum6NMQFhamWJanR48eSjPn79u3D9OnT0e1atUAAA0bNlSKpzzX7+LFi7CyssKaNWuQkZGB33//HRMnTkSrVq3KPHdtrkloaCiaNm0KAHB2dlZsDwwMRFxcHNasWQPGGGrXrg0A2LBhA5ydnREeHq7NpVMyZcoUjB07FuHh4bCystJ4jJSUFNjb26NKlSoqrwWKEqIHDx6UeSxtr68+7623vdZQ1/bp06eYPHkyNm3ahH379uHWrVs6n7829Hl/aYrVz88PHh4e8Pb2Rq1atcoVj6FREkQIT96WuHApLCwMLVq0UDw+evQo8vPzERISgkWLFuHq1auQSqUICAiAVCpFly5dlF5vaWmJS5cu6Xz84OBg3L17F1OmTEF+fj5yc3MxZMgQRQfprKwsmJubo379+qhRowZ+/PFHle1A0TfNxo0bo0ePHhr3Kal27do4d+4cBgwYABsbGwwePPitsbZq1QqnTp3Cl19+iVevXsHOzg4A0KJFC1y6dAnTpk2Dn5+f0twkeXl5iqRoz549aNRI+f+2PNevZG3d9u3b8e233ypuNlFRURrPvazrVqxly5Y4e/YsACAyMlKxvV27doiOjsbixYvBGMOxY8cAFCV3GzduxLfffotDhw5pFX9pxYmBpmP8999/am/g5aHt9S3rffO299bbrq+hrq2/vz/atWsHABg3bhyWLl0KxphSEq7v55OL8y8rVhcXFzRo0ADr1q3DoEGDcOXKFbi7u+Ps2bOIjo5G06ZNYWVlhQkTJuh1DrqgJIiQSsDf3x8dOnRQPD569CimTJkCR0dHREVFoXXr1kr716xZU+dj/fvvv/D390d+fj4aNGiAnj17wsbGBqdOnUJsbCwiIyPx+++/w9bWFgsWLICHhweys7Oxbds2ANC4PSMjA4MHD0ZQUJDGfYrdvXtXEUP//v1x5swZDB48GFOmTMGSJUvQtGlTjd/AP/zwQ/j4+KB3796QyWT4448/AACjR4/GkSNH0K9fP7i4uMDCwkLxmi+++ALjxo3DgQMHYGVlxclCzuvXr0dGRgZ+++03fP/99wCgOP/yXLfS53by5EksWLAAdevWRWRkJK5evYrhw4fjxo0bWL58OfLy8uDm5oZt27bhzJkz2Lt3L3755RfMmzcPmzZtKvPaAcrXvnXr1mjfvj0AqD0GAJw/fx5///23yms/+OADSCQS/Pfff4iNjcWjR49UkvPy0ue99bbra6hr6+/vj3nz5gEAqlWrhtatW8PPz0+vZqbSn9F27drpff7axFpcm9mkSRMAQJs2bdCsWTMMHjwYEydOVCRBb3uPcUmnGaNNCc0YTYiqI0eOQCQSISsrC6NHj8ahQ4eQmZkJGxsbZGZm4rvvvuM7RIM5ceIE7t+/j9WrV5f7tYWFhahatSpu3bqF1atX4/Tp00rbCwsLMWrUKKxbt07xjbgi0XTtfH194eXlpXa5o7fZv38/PvvsM633j4iIwJQpU+Dr61vuYwmZPu9LUzBhwgQsWLAAx48fR05ODoYNGwZLS0vcvHkTjo6OqFmzJiZMmFCu68DF/Z2SIEJIpRIdHa3SXKWtjz/+GCkpKcjPz8eOHTsUo5uWLFmCGzduQCqV4quvvsKUKVM4jFg4NF27p0+fYsOGDejXrx+mTZtmsOOHhITgwIEDyMrKwoYNGwx2HD7o874UsqNHjyIjIwMFBQWYNWuW0nPbtm0DYwyzZ89WbCvPdaAkSAuUBBFCCCEVj6DXDiOEEEIIETJKggghhBBSKVESRAghhJBKiZIgQgghhFRKlAQRQgghpFKiJIgQQgghlRIlQYSUsviYP+Ry5SGX+VIZlp58BpmcQVIoAwDkFcj4CI8QQghHKAkipJTDD6JRetaJf29H4ujDGJx6Gov/eQcAANr87xx+uxCCR1HpeB6Xodh3791I7L0bqfi9+N+bYSkIT8kxyjkQQgh5O0FOlujn54cJEyZg6dKlGmdevXnzJtasWYO8vDyIxWIMGTIEK1euhJmZcl5HkyUSTcZuuw2PQa3wXvPaim1Po8UYs/UW/pnaHYkZEiw+7o9Xq4fDebGPYp9mDkXrar1KVk5o1nz0Ls4/T8DVkGQAwNRezvjnVgTGdXXCkYcxAIBNn3TC6M4NDX1qhBBS4XFxfxfcAqonTpzAkSNHYGtrW+Z+S5cuxYIFC/Dhhx8iOzsb3bp1Q926dTF37lwjRUpM3YPIdCRk5gEAvj/qDwA49CAaADB5133FfktOPFN6Xenkp9ji48r7/XMrAgAUCRAAmJmJkFsghZlIhEKZHNaWVfU7CUIIIToTXHOYi4sL9u/fD2tr6zL3Gz16NEaOHAkAsLKywogRI3DhwgVjhEgqkPmHniImPReHHkQrEqDS9t2L4ux4ZiKg0/KLWH0mCO96Fr1fjz+Kway9D/EwMo2z4xBCCHk7wdUEOTk5abXfvHnzlB5LJBI4ODgYIiRSgWzzfQkX51ro5myv2Lbndb8dY7j3Kg0FMjl23yk6JmMMCw4/BQA0tq+Brk3sy3o5IYQQDgkuCdKFTCbD+fPnsXfvXr5DIQIx7+BjbPq0s+JxpqQQ/zsZgJNP4gAAozs1UDy349oro8VVOuFq+sMZxe/puQUAimK1ed1MVvJ3Qggh3BJcc5gu1q9fj9GjR6NHjx4a9xGJREo/np6exguQGNw235cA3nSQO/kkDk+ixbj9IgUA0G3VJUUCVPy80Bx+EAPGGDp4XkBIQhYAoIMnNfESQggAeHp6Kt3HuWDySdCZM2fw4MEDrFmzpsz9GGNKP5QEVQxXghNxJTgRv54LxvO4DHy15yGy86UAgHuvUrHt2kusPhuEAqmc50i18+PrTtiDN13nORJCCBEWT09Ppfs4F0w6Cbp79y42b96MPXv2wNzcHGFhYXyHRAzkzstUbL/2UmmbOLcA7l4P4O71AABw4H4U4jPyEPF6Lh6RCLgRlmLU5i59Hbj/pnP2vntFTWcdPM+bTBJHCCGmxGSSoJSUFPTp0wdpaUUjaIKCgjB37lzs3LkTUqkU2dnZWLFiBc9REkOJE+chOD4TBVI5JIUyZEoK0WnFRaV99t6NQkBsJn45EwQA+OVMMB+hcmbJiaJJGTMlUnzyf3d4joYQQioewU2W+PDhQ3h4eODJkydwdHRE27Ztcfz4cURHR6NLly54/PgxnJyc0KtXL9y+fVvptU2aNEFERITSNpos0fTFpOdiw4VQnHgcC+faNRCRmst3SLxwa+WA3yd0po7ShBACbu7vgkuCuEZJkGmSyuTwj81Al8a1sOK/QOy6Fc53SILw3ze9AQDvOpU9mSghhFR0XNzfTaY5jAhXana+4vcCqRwZeYVavzYtpwByOVMqAwCyJFJM/ccPSVkSwSRAHQWQeIz88yYWHH7CdxiEEFIhUBJEdCaXM8jkDF1XXcLDyHTI5QzeT2Lheeq51mWM+P0G4jLy0HXVJaXtUjlDdr4U3X++zHXYOts33RXOtWvwHYbK4q4lSWXUgZoQQrRFSRDR2f/deIX/u1408uqnkwHYdSsc23xfaqyanP7vA/x7J0JpGwNw/nmi4vHZZ/F4HJWOiTvvGihq7VhbqM4jamVRhbO5KfTxIikbQNH0AGk5BTj2em2ypCwJhm6+wWdohBBiUigJIjqJSs1FeHIOfj1XNAIrMD4TkkIZ8gplin1G/HEDzot9sO580T4XAxPxP+83tUR+EWlgDFh5OlCxbda+R0jIkCAzT2qkM1GvYa3qWu87paez4QIpw5arL/EqORseR54iU1KI53GZyJRo3xRJCCGVHSVBRCdetyNUFhxdfyEU8RkSxGVIkJMvRUBsJoCim3XpZpqEDAnGbb+DhEyJYpvzYh8ARfP7JGRKIJPz35l9RIf6So9FAC4t6Ku0bUh7RyNGVGT2vodIzJRg9r5HAICn0WL87BNk9DgIIcSUURJEdJJYInkp7X54GnbdVO7M/OfVF4rf86UyuK7W3NenQMZ/8lPsz8+6KD02NxOhODdbMqwNAMDB2sLYYeHMswTEpOchKauoQ/nWqy8VzWSEEEK0Q0kQ0UlcRl6Zz2+4GKr0eNOlN7N5t1p6rszXfnvgse6BcYQxYOvELirbd01xQdM6NQEA0/s2AwA0d7Ayamzq3HmVCgBIzMx/y56EEEKKURJEyuVSYCJaLjmLx1FivkMxuGYORcmOXY03kxM2sq+BquZmaFXPmq+wyi0lOx+PotL5DoMQQgSHkiCitYDYDHz57wMUVOBh2B91bggAsKleBVXMikaCeY5sp7LfuXl9jBpXeXxaaomNR5Hp2Hr1pYa9CSGk8qIkiGhtxB83+Q7B4MZ1awQAODDdFS3qFtX2jH6dGJUkhKHymtx9lYb4jDz87BP49p0JIaQSoySIaEVSYuh7ZVDF3LQ/GgfvR+NSUBIAvK65E05nc0IIEQrT/ktPjKb1T2V3ZibCsvlyGOSM4VpoMr7Z/6ajefGoPbkAph8ghBC+qU6LS0gJqdn5CE7I4jsMo2FgUNfSdegrV7X7q9tuY1kFmRJ+J3sEgPgMCf7nHQCgqIkMAFacDoSDtQVuv0zB4HaOcGtVl88QCSGEV1QTRMoUFJ+FiTvv8R2GUXV3tlfZ1qNZbbX7qtveqXEtzmPSRYFUjsjUXABAdr4Ud14WDaOfc+AxxLmFyMmvXE2chBBSGiVBRKPbL1LgcaRyrFjeqZEdZ2VpWjuNbxP+erMem5mZCPlSGV4l0wSLhJDKi5IgopF/bEalmXyvf4lmIS4GfpUsw9pSeK3OjDH4RaRj+r8PVJ7zOPwUKdmV4/+dEFK5URJENFpzNpjvEAyqo5OtUY5TVYAjzc48S8CB+1EQiUTIK5Ap1V4dexSDbqsuKR7nFcioIzUhpEIS3ldUIgjXQ5Pfuo9QOgDrykpNDU1da0u4auj/UxYHawv4zOkNqZzh+2P+6NHUXtEZWchEACb9fQ/vOtkiN1+GXz56V2Wfgb9dw/tt6qKejSW+7t8CMjmDCEVNaoQQYsqE9xWVCMIXu+6/dZ8jM3saIRLDYQz46vX6XwwMZ77tgxZ1rTBvYEudyoIIaGBXHQAw261FieeEW4sSlpSNB5HpOOQXjUMPotFp+QWVfWLFeTj/PAHrzocAAH6/HIZ/70QYN1BCCDEASoKICm1v2q0crRVLS5iqNvXfrAHWtoGNHiUxiFB0LQSc82iUW1A0UiwrX33NXsnWMEmhDOm5hTj/PMEYoRFCiMFQEkRU9F13Vet9TfB+b3D1bS1R06LslmZzE0selfoEiYoWZV3xHy3LQQgxbdQniCjxfhKL6LQ8rfdvYGdZrv0rqpK1P2vHdlBaW0xdoljN3Ax5cmHP03MlOBFWFlUBAKk5BQCA++FpEEGEOHEeYsX0/04IMW1UE0SU+PjHa3zu1S/DAAArR71ZVd3Koiqa1qlp8LgMjYsmrOK8p/TiqowB7zY0zkg0Lrl7PcD4Hcor0o/fcQfbr73E1RDljvOMMQTEZhgzPEII0RslQUTJhcBEjc8Vjwaa9J4zfhjaGkDRzc+0GnbeYMx4/XdMrPVLa/MPFU2myRgw4o+bPEdDCCHlQ0kQUUh/3eShjRn9mr95UEFv8OVRVi5V39YS410aGS0WYyruHK1pcsUFhyvHjOOEENOkVRJkZmYGc3NzrX6aN2/+9gKJ4Hy+8x46r7yo02s/aFOP42iMh4vZoYGya8Tq2Vjis+6NuTmQwEjlDD+eeIZMSaHS9rTXCfXxR7F8hEUIIVrRqmN0x44dsWnTprfuxxjDjz/+qHdQxPhuvkgp83l1i4oWG925IWLS8+DzTHN/IqEqbg6b5aZ/8l66L5C2z5k0Buy/F4VpvZsqbe6y8iIi1gznKShCCNGOVklQ7dq10a9fP60KrFGjhl4BEePbdClU43PVqpihQCrHyE4N1D5fnERsmdgFPot9DBGeUVhWNdfr9ZV1qoACmRyA+hZRSaGwR78RQohWzWGXLl16+0467EuEYdOlMLXbnWpVx+c9mpT5WgbGWZOSqSvPZWAVLG0qrunKl8qQW1A04eLK0zSPECFE2KhjdCXGGENSpkTj88dm9UTHRrZoZF+9jDKgmCl5Sk9nrkM0KAaGoe3r8x1GhXDrdXPqbxdC0fZ/53mOhhBCtKNTEjRq1Ciu4yA8eBCZju6/XNb4fD0bS4zq1BB933FQbLu+qL/SPgNa14VdjaIJ9YS8RpYm1avp1wxWzARPnVNLTwYAAHZcf6XYZqahilCbZjJJoYya0wghBldmn6BffvlFZRtjDEFBQQYLiBiPTK75zr2mxGri8wa2hGXVony5cW3lPl8/DGtjmOBM0NuaBa0sqiBbw9pcFdGeu5Fqt7dbdh4PlgzE0M03cPfH9xXbJYUyRd+sMVtvIyg+E0//Nwi2r5NsQgjhWpk1QZs2bYKTkxMaNmyo+HFycoKFhYVBg/Lz80OLFi3g5eVV5n779u1D165d0a1bN3h4eJhkTQSf9t+L0vhcySvpYG0Ba0u6Eelr86ed+A6BF5JCGeRyhux8KSSFMsjkDFeCk5CQKVHU9iRlStD6p3OK1+RTLRAhxAjKrAnq2LEjBg4ciAYNlEcGXb9+3WABnThxAkeOHIGtbdnLDAQEBMDDwwMBAQGwt7fH+++/j61bt+Lrr782WGwVzamncWq3H5n5Hlo7Wqt9jqjn/XUv2JRKFLs3tcf98DTF48qaow9Y74tJ7znj13PBqF2zGgBA/vpitP7pHAKWD4bs9eMCqRzVqlBXRUKIcZT51+bixYsqCRAA/P333wYLyMXFBfv374e1ddk34b///hvDhg1DnTp1YGZmBnd3d2zfvt1gcVU00tdDm0s7+XUvuDjb61TzU0nv8QAA5zo1FcuKFPv90848RSMscRkS/HouGMCbhVj/uvGm71BIQpbid4lUuQaooo2iI4QIi+C+cjk5OWm1n5+fH1q3bq143LZtWzx//hx5ebSytTZaLDmrdnvT2qa/GKq2KmvNjBCEJmYrfg+Me7PwKpMDGy+GIrdAORnKK6DmMUII98qVBHl4eBgqjnJLTExUajKzs7MDYwwpKepnPhaJREo/np6eRopUeHLK6pyrx5w/NF2QMpo/SXvFCenl4ERsvhyGhFJTN3RcfoGHqAghQuLp6al0H+eCVjNGF7tz5w4nB+WKuougqXM0dZp+g24oREgYgDVni5rLFhx+qvzc649tgYbmW0JI5eHp6alUgcFFIiS45jBt1a1bF2KxWPFYLBZDJBLBwcGhjFeR5f89h7SMofH6oDSzbFQzpF5Mep7GTvqyUl9eVvxHs1ATQrhjskmQi4sLQkJCFI8DAwPRrl07VK+ueXbjyi4+Iw//3IpQ+9w/U1zg7zkINpblqhwkZSid9FBlpHr/V2KCxdL8Y8T4v+svFY93rhdhdwAAIABJREFU3QoHAOy7F4mA2AxNLyOEEK2YTBKUkpKCPn36IC2taMjxl19+iTNnziA1NRVyuRxeXl6YOXMmz1EK23urr2h8zq2VA2wsq1bc1c7VMHRO4mBlga0Tu6g/NiVEWimQyrH7tvKki9FpubganIzI1FyeoiKEVBTlSoKM0a/m4cOHcHNzw5MnT7BmzRp89NFHAIC8vDwEBwcjN7foD1/79u2xfv16DBo0CD169ECnTp0we/Zsg8dnqj7edrvM57lIfujGrkwkEnG2LEdl9Sw2A7HiohGfV4OTAAB91l5FYFwGvt7/SLFfYFwmL/ERQkwbZx2jc3NzUaNGDY3Pa6tr167w9fVV2d6oUSMkJycrbZs4cSImTpyo9zErg4eR6XyHIBhmIsBA3aLKRDli+ZX8f5rq5af4PS5DefTYsN9vwHehG2paVIGDtWFntCeEVBycNYeNGDGCq6IIx64EJxrlOKbWkuZoY8l3COQtylO7+NeNVzj3PEHxePwOYY1mJYQIT7mToMuXL6NDhw6oXr06zM3NYW5uDjMzM1y7ds0Q8REORKfRBJLqGGMtr3cb2uLr/i0A0DxKurgcVHYCP/z3G/j3TgQA5Zq2R1HpSkuWEEKIOuUeCjRnzhwsW7YMXbt2RdWqRUsrMMYwYcIEzoMj3Fh26rlRjmMqfYKKwzRGJ/A6VhaoY2WhdFyivbCk7DKffx6XiTjxm6axn04GYJJrE/x7O8LAkRFCKoJyJ0FOTk745JNPVLYfOHCAk4AIt+R8dH4hxIi2XysaQl8yCTcrkeAyxpAvlcOyqnIndUmhTGUbIaRyKXdz2EcffYRLly6pbJ8/fz4nARHu5ORLseXqC77DIMQoDtyPAgBkSgoVtXzxGXm48zIVX+y6r7J/65/OIVach+x8KWT0ZYGQSqncNUE9evTA2LFjIRaLFWt3McaQmGiczrdEO4f8orD4+DOjNlFVMTeNXi98NdvR0i3G0cHzAsZ2LVqIuXhuLBvLKsjOl6JQKketmtUU+8rlDLP3PYJ7L2e4tarLS7yEEP6UOwmaOnUqZs2ahc6dO6NKlaKXM8aoJkhA0nIKcNAvWuubfR2ram/fSQuLh7bWOCM1IcZkViofz5RIcfRBNF6l5GBsVyd0cLJTPHc9NBmT32uiUsbX+x5hi4bJLgkhFUO5kyA7OzssXLhQZfu2bds4CYjo78yzeDyOEpe5TzVzMxTI5JgzoAXmvv8OJ8e1qEL9K8pSmWbj5tu10GSVbZZVzSEplOHDP28hYs1wpeeK19NjjCn+n3yexWOL4UMlhPCo3H2C+vbti+DgYJXtu3fv5iQgoj9t+jdsn1T0DddjUCtUMTeZ1VNMGjWHGU9iZr7KtirmZopk596rVABvPityOcPdV6n4n7dxRlISQoSh3DVBN27cwMaNG/HOO+/Azu5NlfKTJ0+oNkggtBkSP6B1PSNEQohwiPCmP9jJJ0Wr1n9zoGjpjSMPY/CpSyPEl5qJmhBSsZU7CUpISMCff/6ptI0xhvj4eM6CIroLScjiOwSiBaoTMj4zsze1ccUjyTLyCgEAV4KTkJ5bALmcKTWJrTwdiJ9GtFWUwRhDUlY+6tFs44RUCOVOgpYsWYLPP/9cZXvJWiFifPlSGV4m5eB+eCrfoRANKPHhlwgiRQ1QsZItlMX96NJyCnA2oGj5jb9vhmNUpwaKjtSSQjn6rbuK4JVDFa/LK5AhKi0XrRyttYojPiMP5iIR6lIiRQjvyt0Z5O7du9i1axeePn0KmUwGAJBKpahXj5pX+BQnlmDY7zfg+V/gW/ddP66j0r9c+ty1MedlEsKF0ETtakk/+b+7WHoyQPH4wz9v4cLrNclEahbfDUnMwqKjTzWWd6HEemaPo9Kx/nwojjyMKUfkhBBDKXcSZG5ujp07d6Jnz56wsrJC9+7d4e7ujuHDh7/9xcRg9t+L1Hrf4jlUiv/lkjmNgNKIrgy/tvq+VNmm7u36osRSHcXPf7XnIQqk8qIHrKgjdYFUjnypDFXMRBoHIzyLycBXex4iPacA47bfxl83XuFhpDDXNMuXyvgOgRCjK3cStHnzZty+fRuZmZm4f/8+xo0bBx8fH0ydOtUQ8REtPI5Kx183wvkOAwA1+ZCKpWRz2aKjT3HILxoQFTWTtVx6Fl1WXISZSDkJ8jz1HBEpOQCK+hoBQKw4D34R6SXKFd4nZcD6axDnFiAnX4rcAqnafZKzVEfdEWLKdB4bbW5ujnfffReLFi3C8ePHER4ujJtwZRQQl8l3CApyAf5xf5uFg1ryHQLhSXRantb7SgplWHbqOQqkcvx8JggAkFMgg5lZ0fv+abQYT6PF8LodAbf1vgCAjZdCVcqJSM1FpqQoyXgaLcbcg48Vvxcr3lbStwdUt3Epp0AKxoBtvi+xU8OXKpefVZdMIsSUcTJBTJ8+feDr68tFUcTEmeISTN8M4GaySFKxnX+ufmkgc5EIoYnZGLXlFkZtuaVVWbkFUjgv9sGoLbfg/bqzdsnXepfowM0Yg1zOcOppnEo5XLgRloyXydkQ5xaNlGNgZTbdarMo808nA7DuvOp8coQITbmToE6dOuGrr77CX3/9hcePH0MqleLq1ato0kR12nlieDuuvcRPJTpx8k2I1fxCQVemYjIrvUbHa5JCmdrfAWDv3SilxzuuqfZXKt62/36UxoWQc/LVN1uV5cTjGPxw/Jni8caLoTjsF614zJj6vlLFev165a3HCEvKwparqudEiNCUOwlatGgRbGxscPDgQQwcOBAWFhYYNGgQGjZsiD/++AO+vr5ISUkxRKyklPCUHKw+K6xvW6aQA83o24zvEEgFomkwwPISIzXnHnwCoGiYvjrqPsfF27IlUmRpSHY6Lr9QrlgBYP6hpzhwPwovk4s6gJuJRIpm7M4rL4KhaImX7HwpEl5PHrnnbiSC4oua3WlCSVKRlDsJmjhxItavX4/Lly8jNTUVr169wvHjx9GtWzdcunQJX3zxBQ2XN5LxO+7wHYIKU+gT9MOwNrwcl0aHVUyamqmKJ2QEAKm8aGTZ1ZCkMsuatfeh0mPGmFIN4pXgRKXnpGqapkruU5b3N1wDUFTrU3JgRfFH+FJgIn553ffpemgyotJytSq3ZBlZkkLM3vew7J0J4ZHefYKaNGmCUaNGwdPTE97e3oiKikJSUtkfdMINDbXwvDKBHIg3dGkqpt8uqnZ+Lq14LbPcgrKHoZ8NSEBK9psRWE1/OAPgTQLt7vVA8dysvY+UXpv3uuyS+9x6kfLWWeRLjloDXvcJEgHzDj2BnDFICmUqzXnFdt18kzwVSOWQyuRKzxdI5bjzkiZwJcJlkJUza9eubYhiSQkvkrLULhLJt0nvUd8wQvTxKFI5KSmUyvEyOUdlv3MlJmEEgI4rVJvGTvvH44EO8xKtPReieH3rn87hRpj6Lg4rTgei8HXis/FSKJb/FwiZnCn6FIlo3jAicLR8uAmZ9Pc9AEUTsA387TrP0ajXwckO/Vo68B2G8FG1ENHgqz3KzUcH/aJxKehNE1fp0VnFc/oUSOVKzwXEZiAhI6/M2tll3moGVWjY/+4r9TU6MjlDXoEM23xfYs/dSKRm51ONMDEZlASZkOJvY7tuCXtOpqrm9LZSh24MRBexYuW5jPqsvar0uGRH5ZLP7b4dgRthKWAAsjV0rN59R3WmeU39+v65FaH4vWR5IhFQULIZTATcC39T+0RveyJkdLcSgExJoWKGWW329XkWb+CICCFCFSvOw/dH/RWPRaWeA4o6TRfK5JDKGQLjMvDx1ttYevIZAmIz3lq+NrPPZ+QVKj0uOTVGyRFwqdn5ivmHAGh1fEKMqcwk6NGjR+jbty8GDBiACxfetDePHTvW4IFVJrfCUrDK5+0LnwLA8lOBb9YwIiaFukcQrhx99GYBVvnrhKekmPQ8RZ/BA/ejIZHKcPRhDEb8cZOT4xdK5YrE58SjWKUvZiXHsx0sMf/QicfcHZ8QrlQp68l58+Zhx44dqFatGubOnYuIiAh89dVXSE2l3v5ciBPnYc/dSOy6GY58LRObY49MZ/XpVvWsEaLlyt2VATWHEa6UXKts4G/XYVejqtLzDyPT8f/snXd81PT/x1+567i2dJe2lEIXtKWM0tLJamkB2SAqU6CKCG4QF6JSxY0TtyJDENz+HF8HghRFQIaAKFOBgmUXCrS0hbaf3x/Xu95I7pK75JLcfZ48+uAu+eSTdy7JJ++8P+/RaHHB1V0V7+Vp+vLtuL1fEgDgoS92Y3RGW+O6Dza2TLG93xw99twP+/AWSwFbCkVubCpBDMOgUyd9TpWvv/4aN998M86cOUM9/kWi57P2M68a4CpoqGSIQrwBpvVJUEyBWQpFCkynnAB9eHtWXKjxe3kl/xw/fKiubzAmgASAL3ZUGD+/zpLdmipAFKVi1yfo/Hl9uKZGo8HSpUtx4sQJbNy4UXLBKHoe/+ZvvL/hMP6qUE6RVLXh52NT16dQ3JJtFqH2YqJxwpv0yNkaq3xCFIpc2Hw6LFmyxGrZ66+/jhtuuEEygTwFPjW2thw+ZxaRoTa4SgR4KmrIpk2h8OHYuVr7jTgofKEML41Jx+jMWBElolAcw6Y+n5iYiNDQUKvlBQUFkgnkKRgywdpCiWUxVIlClI8QP2/7jSiyML2A1pOjUDwRQUbN2bNnSyWHkbq6OpSUlCAvLw9ZWVlmUWmmXLp0CZMmTUJ2djZycnJQUlKC6upqyeVzlhMXalFmp34QADzw2S5J9t8zyXXZvJXiE2QgNTpQ1v3nu/C3pwiDqwgqhUJxbwQpQZs2SW+ZKC0tBSEEmzdvxsqVKzFu3DicOmVdEHD+/PkoLy/H5s2bsXnzZpSXl+PJJ5+UXD5nWbXlGJ61U/l9yKu/4pNt0kSBdW0bLEm/5ihL+TFw38AUuUWgKBRP1oHyEsNcvs/aq41YvlkfRVZrp54ahSIlivIYbWpqwqJFi/DFF18AAJKTk5GRkYEVK1ZYWaH+/vtv5OTkQKvVAgBycnKwY8cOl8ssFAbgLEb4z+lqdIhshT0npHOClquCOoWiZDzZf03nrXX5Pud+qS/X8cuBM/jzvyoUJkfiueu7uVwOCkVRGaMPHTqEyspKpKamGpelpaVh27ZtVm2HDBmCn3/+GTU1Nbh8+TJ+/vln5ObmulJcwTz57R68uvYgjnCEq9qzEFHEozDF9fXNTFNLKG2qkOK5yDkV+NOeU2hoJPh42zH7jSkUCVCUEmSY9goObpmyCQkJwenT1j40d9xxBwoLC5GYmIj4+HhkZmbisccec5msjmBZA8iSczX12HZEeMVnpaKYt2uWQX7pTTkyCEKhKA8vrbz3qeH2/O2fs/jzvyrjcj4RtBSKswhSglx1UVomY2Tb71NPPYVdu3ahvLwcR48exaFDh7B06VKbfZr+lZaWiiy1fU5fqre5vrq+Ade/7T4RYdTaQaEoHy9nkv6Ign68n7jod8z6uCUB44CXfzF+bmwinEVgKZ5DaWmp2XNcDBTlGB0ZGQkAqKpqeRuoqqoyLjdl4cKFmD59OnQ6HXQ6HWbMmIFHH32Us29CiNmfHErQdjvJyw6cUn50m5rxZOdXim08+doYm90OwTKmbzhb3fJy+O+ZGmNJkH9Ot4yH6/adxsyP9D6fhjGc4nkYAqfEvAZEewW4fNn5tOxJSUkICwvD/v37jcv27NmD7Oxsq7ZXrlyBt3fLjevt7Y1Ll5Rbp+qoyGnrKcKh4yaFCw/WgeDvo0WgTjkxMkkPf4c5X+wGAJy8UAfAPNHo+xsOG2uSUSjOIpoSNGzYMKf70Gg0mDZtGhYvXgwAOHjwIHbu3ImJEydi7969KC4uRmOjPrKqf//++Pjjj40a4cqVK9GvXz+nZZCKuz5SfuSaVEQG+sotAoViGxFNQX2TXe907wwaDaO4F4RVW44CAPKeWYvyyhrs+q8KBlX1Ul0DLtbRqTGKOAhWgtauXYtu3brBz88PWq0WWq0WGo0G69evF0Ugg7krLy8P48ePx6pVqxAdHY0LFy5g3759uHpVXyjwjTfeMLbLzc1FQ0MD3n33XVFkkIJ9Eoa9KxXDwOrn4/oQXANFqZHQecvt80DxJJ67rqvcIgjCS8OYVaW3R0Qr177U7K64gDfW/WvUU9n01Xlf/eVSmSjug2Ab6F133YV58+ahR48exukoQgjGjx8vikA6nY7VwTkvLw8VFS2ViiMjI/Hxxx+Lsk9X4OOlQX2DZxQNVNJb5S29E5CXGI7nf9hvvzHFYxFzOqxNsJ+IvUnH5Pw4fLCpHFoNg0ZBN61rb/Df/jkLAKiua8Dn2/VJZFdtOQqdtwa3F3YAACzbVI7HR3ZxqVwU90DwK3JsbCzGjh2LDh06IC4uDnFxcYiPj8eqVaukkM9tuCSz+TailY+s+5cLhmGg0XiyxweFD57oGD2oczQAQMMImw5z9UvOZ82Kz4FTl/D2+n/BgMGZS/V4/of9KNt/WpAVi0KxRLASNHr0aKxZs8Zq+axZs0QRyB05dk5+p+hv7uqNQF/XOj8q4cGiBBkoFEeID/eXdgfN94aGYWCw7gzpGi3tPh3gaqNeNgK9rMfOt4ynJUu24vDZGpkko7gDgpWg3NxcTJ8+HeHh4UhMTERiYiISEhLw448/SiGfW9Dn+XWy7n9cdjuE+vtg+6MDXLrf6zJjEROsc+k+KRRHUExiTxM6x7iizh+gYYAmAjw+ojOv3CvPXSdPeYsmQsAwLZYhAw83R5KZsnzTESzfdMQlclHUjWDTwE033YTbbrsNGRkZ8PLSb04IoZYgBTMgLUqW+kCdY4IRHxFgN1M2haJ2hnZtg//tPiFupw7qZZGBvnYTs5rtxmQ/fHbZPy1KuFAi0NREcJTFqr6FJcv+o1/9DQCYlB8vtVgUlSNYCQoJCcF9991ntfytt94SRSCK+pk7tBPW7tOXOnH1dNQtvROwSAU5RJTkPE5x/jqVIju6oyLxPRaD9YthGGPiObGy8HJhsDo5AgG/34QmUqQIQfB0WN++fbFvn3Whz2XLlokiEEX9JLZuJdu+qQ8QxV2QWiEx4OetRb9UfVZ+qff4hBMRXJfqGlBzpZFzfd3VRpwRYAGjUAAHlKBff/0V2dnZyMzMRFFRkfFPTeHqrqTuKvdN6yqoYqCH/g4ULtzp0kiOCuTVzlAqI8DHCy+N6Q5A+ntEyv7HvLMJ2U+toVZWiiAET4edPHkSr7/+utkyQghOnBB5PtwNIIRg8vtb5BZD9YPCyO4x+Grncaf7oUoQhQslXBt39EvCG+v+NX53VKTWPJMZpsUEWe3IdJ/tw/xZfXC4mFGQhLfX/2u/oUT8+d8F42eGUf+4R3ENgpWguXPn4sYbb7RaHhISIopA7sS8r/9mddrzNFw5GNGBj6JWYkPNQ+JdpZiZOUabfPEXkOmdYYC2ofaTRKZGBwmSjUKRGsHTYWwKEACMGjXKaWHcjd8P2VeANAzgrVXAayiFQnEYKZRvR0cFNlFKesZz1jQz7MdS6bqpVzym5Mdx7md4eozFju3/CD3iQu22cRZiIopphXoKhQ1aVEkCXlq9H/EP/Q/7T9muaj+rfzLG5bSX3HqhBFO/XHjwoVME4ConZFfAFh2l1TDgSpxueuyMxXJDtnU2q9Br4zN4y9SrQ7jxc5+OEby3c4QXV7eUyHE0Ue2KzeViiUNROFQJEpkLl69i4c//8Gp7T/+OePpadRVb5Evp8DTR+qJTXBQKf4TeLgznF/H262eSp2z51FzHd8KDN8ta/JIu1jVg3f7Tgvt45P9oQVZPgZcS9Ouvv0oth9uQ/sRquUVQBCW9EuQWgULhjSOGoNaBLQ7IkkyHOWidMpUlrrn0BgP7+g0hhhIatvvkIsTfdn3CR4Z2Mn7OiQ+z36EIHK2sMbMMUSiW8FKCbrvtNuzebZ2anGKOUpN0xYUHyLp/KRLJUShy0z6sxZFZjGtcrOHDtJt2Js7WoQHsSoplxmhfr5bHgkEme8fHABiRHoOR3WPw8th0s3Xf3d0HgPnv5arZR32CRf3OHv6SPsMo1vBSgggh+PTTTzFy5EiUlJRg2bJlqKiokFo21fG+AjMVp0YHIknG5IWuxo1cOzAht73cIngMYtcOM+TgkQPTlzHD/cAwwJOjrBMV+nppTDJHN28P4KUx6egY2crsfgrgES1GSMtv+cn0fACGUHzz37d0RGfW7Xt3EN9fqO5qI/aeuIiVvx/l1V6pL7MUaeClBD377LN44okn8NVXX2HhwoUICQnBc889h9GjR+Oee+7Bt99+i+rqaqllVTzV9Q282gkJPXUWuR0+xRhP3EmxEQKXIytFHh7neHCz0S2WX/FTU4dhR2BzMua65QzKic5bY9bWPERev3B0Ziwy2ocat7G8j4N05tlVDOOMaX9ZcaGc926nNuyh8gxjbokSg4OnqzFjxXYAwLp99v2DCl8oE3X/FGXD62obPny48XNQUBBGjhyJl19+GbNmzcLFixcxatQohIc7dzN7EmtnF5h9l/K9Q87nqBKUF6W+0/GRi76QKgslXM+WsL7ksFw3pu3m2yhdwYAxm/oyfCYWfWx/dADr9rf2SUR+YsuzwPQathSrZxL7M+P568WrUm/Y/8XaqwCAmiv2X1TLKx2LKKOoE0Eq9969e7Fw4UIMHz4cYWFhKCoqwl9//YUHHngAP/zwg1QyqoZX1hy02+b9KVloE2w/qZhY+Ij8ViUHnqgMTLaRn4UiPk4XUCW2vwNAqL98U2RaDQM/Hy0Gd4mGt5Z9TNBqGHh7cTlGmx8QVx9dY4MRGaQDYP83fXhIJ9sNRKCpWe76hiYAwJzPzf2C3lhnHcmrRGWXIh28npBTpkxBu3bt0KVLFyxcuBBt27bF+++/j9OnT2Pr1q14+umn0a9fP6llVT3Z8aEo7hRltVzKe+6DqTkS9s4POTNG0/GMwgclXCdiBRCw9TOrfzIAlpcik6bjstvjocHsiomXRoMBadZjl5ikN08hNjSKN2A0Npesv9xcePVSfQNGvL7BuH7Bjy2RY8Nf0y9XwrVAcR28ymbU1tairq4Oo0ePxowZM9CnTx/4+NgOh6RYM6xbjP1GIqPzcp3/EYXizlg+HNkelkISCArdH982TU2Aj1aDK41NxmW2LMIGy4dWw0DnpcHAztFWbXy9NXh5bHd8uYM9IMaerHysK7f364Dlm8qNiosYPPm/vVbLTGuMmbK7Qr9cxN1TVAAvS9Ann3yC06dPY86cOdi+fTtuuOEGjBo1CgsWLMCOHTukltFtKDBJWz+4i/lA88v90ljSPM2062nHSxEHR64b020Mz01DKQkhVp1Xx3Xn6N++UFx72TVvIK99+3ppzKynXloN3piQ2SJDs3pjmP4SouQJDcowyOHrzf1YEit67MPfyzHyjd8AAO/+Yppc8aoo/VPUA2+HEYZhkJmZiQceeABfffUVPvnkE/To0QNPP/00oqKiMH78eCnlVD2T8+MQH9GSr+etG3sYPxMAfhJFjFGdgEJxDa9PsK0gcEZtcSgL0cH8qsFb74dwjieW08Vb5vaHzpt77DEoc4aXNKt6Yc04okQaZLGMcAvjyGcEACO7i2NNP15Vi13HqgAAT3+3z7i85zM/i9I/RT3wUoIWLVpk/FxZWYlPPvkEd9xxB6ZOnYrPP/8cAKDV0mkXWzxhIyIDoBYMW1DrNEVMkqOs82aJkSfIdLpbiB+cvbw0QkPGheyb78uX0Jc0IWlA7h2g91fKbB9i14I2OjNWkBxcnLrIXliVb5oTivvAyyfopZdewsGDB7FmzRrs2rULvr6+6N27N+644w4MGDAA6enp9jtxY15Zc8AsM6tQpNR/5MwTZMw266QWQ5OXUcSELbKJ121i0Uipl6UcYlkqkTsesw6ht4qga5bUW6tBWpsgRAbpzJItvjc5C9M+2Ga2jVak5FlNFo4/tsaY+z7dhdIRndHKl9fjkqIyeJ3VAwcOYO3atRg4cCAWLFiA3r17U8doE+yFxn82I99uH2KpKn06RuDXg2dF79dRlGLh6pfS2n4jikegNOWF60XFUeuUHMeXb5Hzx9ciIMORI+GKRhud2RZf/OFcxYJGix8pYc53nG0/2/4f+neKwiALP87aK41479dDuLu4o1OyUOSFlxI0ePBgfPPNN1LLQhGBrLgwcyVIIUoIhaIUurcPwZ4TF53ux8wxmkeeIHuIp7y4XgualGc7r1W/1EjERzhuLRebr3YeF9SebRytudKApRuPUCVI5fCabKYKEDuEEFxpaLLZZvWsvs21c7hJbB3AWb1ZKJbdyF02Q0wcSTbnTsdPEYdAnWPTGtYh8sq8tmwpU5ZlMpzpSwjjc9ojNdr2OGjYn6V8luerKDUS39zZGwCMU1Tjc6Sts/fk//bgMku2aWVeARQhqD+dsIyU7T+D5Ee+t9kmOSoQ/j62B90fZ/aFl1YkJUiUXpSJn40oFgNKm+qguDdr7u3rdB9i3rPFqZGqCSQghN89bTggw289rFsMujYnVmzJcSS+fE1NxKj4HDtXi1s/2I67V+3A2eoWp2q1/NYUbqgS5ARXG7mtQPNHdkaHSH7V26m1gh/ODDh0sKLYwlFLbIfIQNblbFFOQhV0R3MXqSWQgADoyBKpxwXbb234icSypJuy8d9KpD++2vh9wz9n8fWu48h6co3o+6LIB1WCnMCW8jIpPx5r7i3gXC+kL6GYJmWUGzHGY7Hqnylp+kItDypPgc/tZ6uN5Tq+p9fWNammS6RL22DW5YNYsk+bwgDoHBOEj6fnAdArj47cpVIoQTe+/zuucpTwuP6tjQDc2/LuKVAlSCBXG5vwV3N6dTFvADH7ctTnQWzEGpf6dhSu1EUFOZZojkLhgivEm20d7z7NKrazc0MWd24ca2WakcXqGR2sY13+9qQerMuBFtkZhkGgjtvfz9bx+Hhp0TrQ9ff6tvI7Bn0YAAAgAElEQVTzVmH2FHViUwkyJEk8fPgw+vXrh44dOyIvLw+7du1yiXBK5Gx1PaYu24pnv9+HC7XKS7FuelvO7K+MqAVnrR6jMtoC8Ky3LjVZAdQGAwYPDErh1TY/sSX027ZST2x847M9d/snR3XlvU2QzgvRQewKiRTwnfLnhE9pEDs3w+Y5xbJEweY/a55des4XuzlaUpSMTSVo5cqVAIDZs2fj0UcfxcGDB7Fo0SLcc889kglUV1eHkpIS5OXlISsrC6tXr+Zsu3HjRhQXF6Nv375IS0vDa6+9JplcBgjRZxt9e/2/+Hn/adY2ttK+cyHFTaz1QF8jdzpkdzoWJSGkrpctZUmM6VVn+7CcRs+KD8Mzo1uUpuev7+ZU//b4aZbzjuGmEAKrt52s+DCb24iVQFEohkKvF+uuoqGxCau2HJVFDopz8JoOO3/+PIqKigAAXbrYLv/gLKWlpSCEYPPmzVi5ciXGjRuHU6dOWbU7fPgwZs6ciY8//hi//PILPv30Uxw5ckRS2QDz/BL/+/MEa5vv7u4juRxcMGh5C3W3h6gQvykl+f84CrUGSYfl9cHnamGZebLRmJ8cBET0+9T0PmkT7Ge+v+aDeGUse9FWZ/YlFM5aahY/7LKbc3j32S6s5Xjnj+zsiFiCYBjgvk92Yc1e8xfidftOY89x53NRUaTHphK0efNmJCcnY/v27Th37hwAoKmpCRcvSnNym5qasGjRIkydOhUAkJycjIyMDKxYscKq7UsvvYSbbroJERH64nudO3fGiy++KIlcZvv9ab/dNlzz47YQ66GtRMWHz/MgopX9eX1HptVS27BH7ygdJZ5Hd8YVEZpcl68zyi4hBImtWwoz87FyMQxjnGKWG2d/dcNpM4yfXpqWR9qk/Hgne+clAesv/sNfJ7HrvyqzZacu1uHAqUsukIkiBJtK0OXLl3HgwAFcvHgRYWF6k2RtbS3eeecdSYQ5dOgQKisrkZqaalyWlpaGbdu2WbVdu3Ytrly5giFDhqBXr1647777UFdXJ4lcBq40NHFGCxi4u6iDQ31LMQYHqKjWzZzBqfYbOUCQDYdLCsUA1/0ntTHO7OXHQhviOyQ8NixNPIFYkMrHiE35U5vx82x1PY6du4wv/vgPAHD7h9uN66oum/uM/rzvNBZvOOxS+Sj2ERwdFhAQgOzsbClkMU57BQe3hFuGhITg9Glr35sjR47g7bffxgcffICysjLs2bMHM2fOlESuv4/ro8Ge/2Gf3bZJzjoKigADIMTfW7a5coo0tA9TTtkBJcMrqabMj9vXJ2SI0o+lBctUqXp1nDhTXtMLEkXpx5KIVj7o0zHCarnpIe2aN5BXXwwDDO3Wxq7i+MjQTgIk5Me+k5eweo/+2fXd7pPGqLHnmp8XQqLIaMSZ63EoRH7kyJFiy2GG5Y3NNg1SX1+PCRMmICIiAt7e3pg5cyaWLl2Kpib2BIYMw5j9lZaW8pZn6MINAIBFdrR4Hy8NOnIkT3MVht9OSf4kfGRx1fTPpzyK2VLcA38fdmXIoCjcf02KyTL7CLmlHFKyRL4JRnZ3fspLynxWceEBmD2Q3fE8NVo/jgb78bfkprUJwst2fJ1cYRm+a9UO4/k/eOoSpizZwmu7xiaCm5dtxa5jVfYbeyilpaVmz3ExsDlf8vTTT1stI4Rg7969ouzcksjISABAVVUV62dTQkNDERXVUmU4NjYW9fX1OHv2LGt7VySn2/XYQPhxDLz2kEIJoHYga7LtRJooEX0WYP3nfimtsW7/GXkFUgleHJZQwwNK8Jjg5Bhid8rNRv9v35iJGSv+YBHJMjRf/HHOpeNIs/jf32M/uGTDg/3MrN0+Wg3S24VIJRlvzl++ggBfvYN2fUMTKquv8Nru/3ZUoGz/Gdxe6JhLhSdQWlpqZsAQQxGyaQl65ZVXEBsbi7Zt2xr/YmNj4esrTXKqpKQkhIWFYf/+FufjPXv2sE6/de/e3Wya7MyZM/Dx8UF4eLhVW1fhqAIkJl3bBituXt0wMJcO5/ZdEFtHdTSDr1KIDW2JcokxifBZNEWaqWhPgwFjfk2I/BZi2veEXGmLewp5EDhyG7j61mHA75hiQ/2N0W85CWHo3NZ+gVZX0NBE8N/5WuN3hgEe/+Zvu9sZ6kfSoAjXYlMJSk9PR//+/TFlyhSzv5wc/iGLgoTRaDBt2jQsXrwYAHDw4EHs3LkTEydOxN69e1FcXIzGxkYAwLRp07Bq1SpcvnwZALB48WLceOON0GrlUURMQzMdQazosF4drOfYlUIfBZXzUDp+LEV36eDID76/E592tpRn081ttYtozhtm2sbwuWtssHG9QwLC3BrEZxxxx8voms7R6Jlkf+xzhS/YlsPnsPHfSrNlS347wtn+eFUtrjY2GUt/cJ0fWmpHGmwqQT/99BNiYmKslr///vuSCWTIE5SXl4fx48dj1apViI6OxoULF7Bv3z5cvar3uB87dixuvPFG5OXloVevXtBqtXjllVdEl+e3f87yajetj3POg5I84BT21HSlNJbjhcJ+CsVZ6zwRU2uD+D5BLbSyUcbmwUGp6NRGPAuG2qfD5HZYFxtLxfejrces2tyybBv2n2wJnecaq4pfWm/83NREUN/QKJqcnoziYqh1Oh2WLl1qtTwvLw8VFRVmy+bMmYM5c+aILgMhBAzDYMlvh/Hq2oMAgGPnLnO2T4lSTi4awxx5SnSg2ZSKnLjyBYa+LFEMuCL3jwGhl52Yogk5zi4xwYLymMlxO0l53uRIomo4nPJzNcZlh85UI7G1eSTxso1HbPZz6EzL9t/9dQI/7zuNl8aIEwHoyQiKDps9e7ZUciiG73afwLPf60MbH/9mj/GhOnkxt4d/+3DnQ5fFvjUn5cWhX4q1g7inoDTLD8V18FGECYj5NJKIU2MAcEOPWJsvIZZ9hfh7i6Jw2HvI31aYJDg4wKXKpMRal5yWpnfWHwKgzyZd9KLeqrP675PG9dvKzzd/sv97NzQSNNjJWUfhhyAlaNOmTVLJoRiq6xpQWXMFdVf1pkZDkdTDZ2s4t4kUoYqxKwcaV5EcFYhpfaXJMWIP9/s1KWLD74HL3cjWLdsuzB/tQvUvR45aH/huJbWviKt9UZwdCrfMLTb7Pq1PgnMdikxZc83J5ZvLcevy7WgiBA1N1gr5mHe4n7dSlFvxVBzKE+Su7DxWhYOnL4EQCEpvrlHY1ZjUuhWvMhRS0zrQF/1SIl367kXfjSh8YMCYPUQcUVRs6QYh/t5Iac51Y7A+mO2Pj+VJsEQUAIgM1JmlR2B7wfTWumbMLlmyBX9VmJeZqm/Q57Jb0px3bt/JSxj1xm/G9Z9v12ef3nL4HA6euoQdR8+DDWU9ddQLVYJM+GbXcfz49ylZTKZiXtD3DkhGfpJ8qQKkwh2tZRT1IMQi0jkmGKUjzAt4GjbPSQgzd5h1Ui7T+yJehKl5W/1LjVgjr71M0y+6yJemssY6R5DBOfoQx+zCh7+3VKNff+AMvmUp1E19H8WDKkEmNDYRHD13GV/8USHoIgvw9YLOS/4cQe6OkIeQO44RVAV0nCAbEVpcFJikdHj62q5W6011A77XW2yoHx4dqs+XNbRrm5a+WNo6cr5TRYw0kwsxrnPTuolsZVS6x8qfVJEPTYSg0aKUxoT3NoMQc+X0jXX/mLXZ/d8FNDXpt6WlOGwjSAly9zwFS02884Uc6f3XpOCGrFin9u3ORg43v2ykgf5oDmOsLG5yT7UNtbaQ2PuJTVe3CeGOqAr190ZJz3j+AsL6fic85AGAAIuErJ1jghHmr88zFOLvLbqirJTL0NfLsff1w88Mwd3FHY3fo4OVETHLlxMX6rB04xF8s+u4cdnGfytBoFcWDc/kBT/ux6otR3HbCn0B19Fv/YarTU146af9WLnlKEvPFAPUMVoEtBrn65iIYXK+rTDJ6T7UDmPxP4XCB86SFjymrSJa+QoKO7fsVwgai1IgDw1ORdfYYI7W6oPrRZtvIVVLGIYxK61RkNwabUP8FP/S+etBfWmc883Taa//rLf0GCKXAQAMULCgzPh147+V+P4vfbQZIXo/t7qrTfi8ucI9G6cv1mHL4XMiS68uBNmIDx8+jI0bN+LkSf0PHR0djfz8fCQmyhMBJCbr9ptXqr956VaZJHGcBwelyi0CBzL4WDEMxue0c/l+RUPpo7TKESO6xuBMLSiJolmpDtaPdjpwarUqYDsvOpYpLaF9KsWqxYdJ7+tTshiixvY3B+oYLEJPf7cXhSmtcdQkf50hohnQXwcMA2gYYMdR7oKsO45V4dNt/yEnQX01FcWClxJ04cIFTJkyBV9//TX8/f0RGhoKADh//jxqa2sxcuRILFmyBMHB6n0juWmJudJzjsWhjeIczlq7nBnD1DQAUlyD6TXB5ZPja+LrJyRPkC2Mt4Gda1IperBS5HAG03P38JBOCG/VUqokNtTPrNaXkrDMBWQ4F+dqruCqxbqf9pwyfm5sIpj75W6EcpVkoRjhNR02ffp0xMXF4eDBg6iursaxY8dw7NgxVFdX48CBA4iLi8Ott94qtawUCm9YzepuMJhT7MNH4eUTEk8Am1NNUkWRmvqwCEXtl7gr3lWGdmsDf5PafD3iQl2wV8doaGriXGfqJ2TKO+v/BQD8e6bGmL7lhR/3s7ZV+/UiBryUoKNHj+LVV19FUpK1z0lSUhJefvllHDtmXROFQpEX5d7ijlimAny9FHxEysGYl0fANrasHYYosVEZbZHezvmoIuO559hnUusAANJXn+eDq9OFBPh4wUvj2qBle76UW+f2d5Ek1uw9YZ6vjo/F6plmvyGthoHBHeqv4xc42xNCBOXFczd4XW319fWorq7mXF9dXY36+nrRhKJ4HnS2ihvDA/rL23vJK4jKsHdNCZ3mCQvwQZugFgdovgkWzabdzErPG9azSxoZKMzZOsDHS/VTV4tLstFNYY7erUWoCOAoFVX2lZ7lm8s5172xTm8VMkTJz/1yNya8txlHKy+jbP9pEOiTN454fYMY4qoSXj5BkyZNQnp6OkpKSpCenm7mE/Tnn39i6dKluPvuuyUVlEIRglITKwb7eRtLsQjFx0tDc36IiN0QeZYGUl1WDMPwSkHSNTYYmw5Vssr+830FZj5MosjlYtujj4Oh8J7Mo//3F/sKswr2+i/r9p3G8Qt1+L+dFdh/8hJGdI8xXtMf/l6OiblxrMVd3RleV9zMmTPx7LPP4vvvv8f111+PgoICFBQU4Prrr8f//vc/PPPMM1QJkogOkeq/GKlTMsWV8PP3cf6i5JsskVNxMt2e503yxMjOSIgIYF0ntgIEyFtwVGr6dIyQWwRJYTt3xy/UAQBe+ukA9p68iD/Kz6OxiUDDMJj7pV6ZKnpxPZb+dtilssoJb7X7hhtuwMaNG3H58mVUVFSgoqICNTU12LRpE2644QYpZaRQBOPuiT0p9hHkE+SExcPhLR26RPV7+/y2no7u1cE9uh/Lp+ZyKpTugGlRVrbh8NCZGrzzyyE0EQJts6a+9Yg+Z9Cu/7h9iNwNwbZHLy8vtGnTBm3atIG3t7cUMimaNydmyi2CxyJ0MFbojJhggnSed5+Jga0pUUeVHraHCV9dxpZizjAMhqXHIKO9PcdrfR+uimhyp3eJ/7vD2qfOtJBqiL/+Psu0ew7UgWl+IFsWvaamlrHyWHPeIQZAfUMjPvy9HF/trEBltfv6/Io2Afvtt9+K1ZWqmDc8TW4RFI9Y46hS/Xyk5qs7qUO0s7ApIN1NH3YOXFr28gyxtTPDYqOIVr5oG8JV1oFh+UQRQneOyD6DQmyIAnNHv6Rdx7gtO38cPY+LdQ3mCxmg6vJVLPhxP97fcBjHFJpHSQxEO9svvfSSWF0pGssB6KZeCbLI4YkImeJSqsLkyDSdt9b9BmVXYO+3NnX+tHW1lPSMR2p0oL6dE5eV1TXZLF6Ivw8+m5FvZ2t5TDIKvY0kwZ3vs+r6Bs4EwKbTZobzrWEY5D69FgzMa5RZsunfSpEldT28osOKiorsttm5c6fTwqgBN7IOuxwPGk8pCoNNKea6Hl8d1x33fNQynvVLjeToUwzJ9PlcDAqZPeU9IaIVPrwlT5wd88CdpsNs8dtDLc84V0fEuYrM+T/ZbWM4359t19cbYxh9Xcy1e09j/rd78MXtvXDkbA0amprQITIQ49/bjCPPDsWF2qsI9lPntD0vJUir1SI6OhodOnTgbHPkyBGxZKK4GXI4KVvu8+7ijogPt64kTnFfTBUKe48107Yju7c1U4Kcl4NrhfUi7ntF31irYQQXa3UaDzAHWU5D+nppUN/Ana3ZU2CgP/2vr9MXcK2oqsV9n+5CcnQgnr62q7Fd+uOrceTZoTJJ6Ry8lKD33nsP06dPx/LlyznblJWViSWTy9lz/KLcInDiTsOPzTBiCfZl2qdSUuMrdZrO3bG89ogLg7+FvgOYXiP9UlobehFPIIpdooJ08NIwOHS2Rm5RXIrl8FRZcwWVJtNovZ79GQDg7+tlbG/IXfbBpiOYnB9v9Vnp8JoEjY+Px0MPPYSKigrONkuWLBFNKFczZOGvcovgFClRgXKLYBNXDt+mFQnoY8MzkUPPtHetFSS3dsgiuuSmHMcEMuHJUV0c3vauoo4Y3q2N0zJQ1EG1pYM0B78cOIOPthwFIcCu//RRaI999bdxvelnNj7YdATnFVKknLcnWL9+/dC2bVvO9fHx8WLIo3joe7w6UL3FxVOcMSRArAKqYpLitGO14/LemBfn8LbRwTqE+LtvJXLLa+WVsd3lEUQh7DnBf1bkcKXeSnbtmxuNy0wz2p+1EVa/bOMRm+tdCa/psMrKSoSHh/PqUEhbuam72gidt7Asq/TRRJETev1Jg7MqkZDt3TkLsxoxVUxzEsIAAHHh/nj+h/0eNx3myLtXVJAvTl3UKzRPf7cX4a30tdaqLl9BRCvrumt3fPiHol5SeVmChGSEVlP26Pxn1uKIwi9ydxguXWnUMNxaqs8YraBBQq2Y1SqV8XowTtEyDDXwKZyYED+MymiLjPahSGztvtmkxcTfp8WWcu7yFTz3w77mb+z18P63+4SLJOMHL0vQiRMnsHz5cl4Dia1q80qj5kojXli9X24xKBKhRDVC563BBffNO6ZaXKlz2pqKozqSa8lOCEMrX/bH4FPXdsWHvx91sUTycqWRf0TcR1uOATAfZ/8oP2/8zDBAwYIyzB6YjFa+XijuFGVcV3e10WxfJy/UQcMAkUEujnwETyWotrYWjz32GK8O27Vr55RArubMJfN5yZz4MGw5cg5+3lrUXm2USSrPQ6zBX+kPESmKXFLMEWu6SSyrjeByL+Ls1mW8NCZdbhEcxjTMmwJ88Qd38JMlF2qvAjAfcy0zSx89dxl7jl9EsL+3mRL03/laDF24wRhWv3zzEfh5a3FnUUfHhXcQXkqQu+YA0jBAzRVzb/jQAH3Cp09n5GPYaxustnH1AKW2AZENSwtiuzA/HDsnrTlESXPOpojxgNYwQJOLtL2M9iFmNYjUir3rwRWXi5UlXekaO09GZ8bKLQJFRg6buJQ0sg1MDPD8D/ux42gV3pucxdkPIfoCrllxocb79fYPt+PNiT1El9kUXkqQu8KAMQsJNK0tY1ldeEjXaHy3+6TLZHNHbD1nlKmyUAJ83GOIYJvKd8RPyFllybA9AVgv+vkju8Dfh5+1kGsah0JRAnO+2A0AeGf9IQDAT3tOsbb7dNsx4xTxhPc24/ERXbD/5EUMT49xyTPXfYul8EDDAEcqLxu/twnWISzA2pvdFMOwmc5RjE8OFGr0EIRoRVYN/VEPVI/FMKC66r6wd62xJWoE2BX/sAAf3hGrfzw6gFc7CkUOthw+x6vd/Z/9afx8tZHg4S9349eDZ3GxTj/dVlldj21H+PXlCB6rBG07cg5XG82HJ0KAgOa3MKtahxYjmRvoHS5DjmSJgHsoh5ZoGGDO4E4u2587/obOYBgHVkzNFbQdY/yfae7H+bvCHaudKxVDAV2Kc3BZg1ZuOWo2dhtyEAFAr+d+xms//yOZTB57F929aoeVJ3xWfCh6d4yAzltjFcHRRLjf3igUsWgdaNsSyTAMpvVNRBFHUU+xUapvlRiY1xYTdpymPwuf34gY/zet2O2+v6278fqETNblL96gXqdwOZj2wTbW5edqruClnw4Yv5u+I9RdlbaGm2hKUHl5uVhdoa6uDiUlJcjLy0NWVhZWr15ts/3Vq1fRsWNHlJaW8ur/SkMTjl+os1qe0T4EhSmRrFEhIX7umzVVciz9QelMldMsLsl2yX7c+TFtsMZwTW1zObHHCSzEa7QCUaXH7aCn1DkMfkNs1JsoP+sPnDF+fvT//sLVxibRan4KVoLGjh1rtYwQguuvv14UgQCgtLQUhBBs3rwZK1euxLhx43DqFLsZDQDeffddnD59mnf//52/zLrc1sM52N8b13SO4m5AURT0geO5CI3Aaxuis3qY2ZrueuCaVIHyNP9vMsAQ4t4KpvvBfk3RFzrnWLWFOw/TbR/+YbXsSkMTlm8uR019g2g1PwUrQd9//z1eeOEF4/fjx4+juLgYf/xhLbAjNDU1YdGiRZg6dSoAIDk5GRkZGVixYgVr++rqanzyyScYMWIE/32wXLg39YpHl7bBrO2HdI3m3bfYuMuzPCyAWtIorkHKB5Oj9yNj9tlNbmoPwF4EHgEwOd/x2mwUYSQ/8j0AoPsTP4nWp2Al6M4770T79u1x3XXX4f3330e3bt2QmJiIwsJCUQQ6dOgQKisrkZra8raVlpaGbdvY5xJffPFF3HXXXdBq+SehY3NKTIwIMEZlEJgPdt5a/c9EBy/H6N4+BFqNtL/d7YVJZt9pdBiFL4b7enh6jHEZre9FAfhF4NGhRt0IVoKmT5+OYcOGoampCdOmTcMtt9yCRYsWYe3ataIIZJj2Cg5uscqEhISwTnedOXMGv/zyC6+pOIZhjH/J0UGo2vChcd247Ha8hryiTpFIaK4n4y4WGlewfGquWe4Trt9uaNc2Du8jNtTaT0OJSqvaBkw1XuctMtsWXmrFnAvS/M/Wb/vqOPNq5rSOlTwYIvDUdt+6I6+sOYCqDR+i/Llhxj8xEKwE9e/fHxkZGTh//jz279+PhoYGPPHEE5g8ebIoAhmw9Olge7OfP38+5s6dy6s/QojxL+7BbxHSe6JxXac2QYgJ9jNpzN7HmKx2SI0Oau6P124pzdDfyz6fzsiXWwSPoo3pPW8BlwLNdh3zsToai6jCNCKNm5Hd25p9T46iIdrOcJuFpVgsqMXZdbyy5iBCek9E3IPfGv/EQLASVF5ejpKSEqxbtw4dO3bECy+8gJycHHz33XeiCBQZqQ/9rapqSdVfVVVlXG7g0KFDOHz4MIqKipze55Se8eifxu30rMKXYUVizJYrwbhh+VatRsfouDBhUUcUdhy1ADrzQON7uRmm2eiz07U8OEiYM7sQ6NSpuhGcd338+PGYM2eO2bJBgwYZHZmdJSkpCWFhYdi/f79R8dmzZw+GDBli1m79+vU4fvy40Rdp37590Ol0KCsrw5IlS5CQkGDVd31Do7Hom4EdLHO+Wx/pD1+FJCKjgyWFIj1sSjPbvcel7Ni7T9k2E6qor5wmLEEjRTxsnV46RqsbwU/6ZcuWsS4/fvy408IAgEajwbRp07B48WIAwMGDB7Fz505MnDgRe/fuRXFxMRobG3HTTTdh+/btKCsrQ1lZGQYNGoSSkhKUlZWxKkAAsOHgWTz0uXleglCWqKVgP2/VWRKeu879qyHzGWsMA5IazdTqk1idDO7SBoUpra2WO5ptXMhY0ZIxmn//BnomRQjfiCIpQk9j7w70HCoNwZYgrumnnTt3Oi2MgdLSUsyYMQN5eXloaGjAqlWrEB0djSNHjmDfvn24evWqMRrsypUrGDhwoNES9Oeff+KLL75g7ffBz3fjbHW9aHK6Ar7j69js9tIKonCEvLUrGaWJrDR5xKC7ZXJERp7jdMff1l1JjOBwTCfgVfA2UOeFS3UNGJPdDhv+OSuydBRnEKwEnThxAg899JDx+7lz5/DNN9/g3nvvFU0onU6HpUuXWi3Py8tDRUWF2TIfHx+UlZXx6lcsBah0eBq+2iWO5YsiLmpUfCji4rCPhslmWg3Du5Cp3W7ZnKlF6ZniKry0zrlHGNKsqNFC7e4IVoKeeeYZjBo1ymzZnXfeiUmTJokmlBSY1iUxsG/+IIf6KumVQJUgkRF7cFDbdKYBOkSKhzOXQG5CGDLas5fTEANCU0ZTKIpAsHprqQABgLe3N/79919RBJKKhWsPmn2fP7KzaG96FGFwPpw8/KHg4YcvGga/G76/J5fTsq+X9fgwJqsdulpklrenwKtUH6dw8PLYlqKpNDJM/Qi2BD3xxBNm3+vq6rB9+3aEh4eLJpQrmJQfL7cIFEsciLDh7EqhZudb+ybisa/+llsMt8YYhm66zMblYMgQz+eBVpiij1itqKo1TnEA/BQd+tLlHpheS84MM+3D/HH0HHsdS4rrEGwJeuutt3D48GHjX2VlJQYMGICPPvpICvkUi4a+3omKFNNXSjxDk1WmfKt1WpGLO/qxJ80TmlsoLzEM70zqIWibXx7oh2A/bwDNipcir1CKUBxVhCxL/VDkQbAlaPr06SgtLZVAFOXC9iBYdnMOusz7UQZp1ImrzMZqfaw8f103uUXwCO5nqQBvuGaEPMwYxtxx2qalqXmdafs2wTokR7Xiv0OK7PRLaY11+8/Q6U03Q7AliEsBuv/++52VRVXYqy4sFu74tij1TFVRaiRGWJQdUDpjstvJLYJbY+/B5cyDjZe1zKJJn46tMb2AWgLUxCPD0qyWKXPSnSIEXk/ym2++2W6bH374AQsWLHBaIFfQjyVRmi1GpMe4pIDhx7fmYey7myXfjxyYKnNSKkEEwOKSbOInWbgAACAASURBVCz97TAOn22QbkcqJjs+FFuPnJdbDMkx1Tv4XHOSKuci+ZFQ5KNtiB8SIgLs+gRFB+lw8mKd3f6M9eQYek3ICS8laN26dSgpKbHZRqfTiSGPS0hpLoLKl+hgHeLCpVeCchPV5VyuZGJCuItjKhVXjYOeMuBmxYfiu90n7bZjGGntraaGIjqVol503lqkxwZbLbe8nWyd46TWLVOghiSLDEsfFNfBSwmaMWMGHnzwQZtt1KQEPTRYumJ6FGUwsHO03CI4hNKekUqTRwivjsvAd7u/t9uOEIIAXy+EtbIuocMHe5GIYkUTUZSBM4psl+b0Ct5aBiPSY3DPR+JVWqA4Bi+foDZt2uCDDz7AiRMnONvYU5IoFAP0bVhe6O9vTd/k1nhqVBeHt7f7m0psDbqrqIP4nVLs4mjAB8Mwbhd5qVZ4KUHPPPMMAP0bz9GjR41/FIpS8NTx5ONb8+QWQRVwPax8vTTGh5EUD6VMCbNOmzJ7YIpL9kOxb9mzvIpu6BHL0omoIqmCGQoNBOClBEVHR2Py5MmIiYnBvHnzkJeXp9owebbq0UrmuevNQ6fzqd+QTTxtuiFBYod9jUYe7ZKzYKUA+Ej+/PXSpib44vZekvZPcS2mirLho70xZ8EN6VbLTJVyT7EI+Xg5V39NKnhJZXqSlixZgpSUFCxevFgyoaRk6U05greR8xo1rXitYYCIQF+rNkE6b1eKRFEIb03MdGg7IYrigLQoh/bhLFqZlC8KxRYG/69uscHw4rhGHVVqHL2fKc7hkGrGdpLnzp3rtDAU23A9u1ZOy3WpHM7CNkbovJX5lqBkBndt49B2QoxlI9JjHNqHs4wVOW+SrfgvMdQtW4plSc943FXUkVdbivJhGCA5KhAMGBAA9w5MFqVfR+9nNaOEhKG8osMuXLiATZs2GbVgy++APk/QU089JY2UIvHcdV3lFsEpuApPe2mVr0CYmn+z48Nw7FyF+Xr6YHAMB343IXXVfGUyYecnOT7tGx2kw5FK85pMUmYst/fm7+fTkinaQ2Y+3BrD7aPV6FVrRxLnmt6CnnxJZMWH4cCpalll4HX2duzYgd69e5sNnr16mc91K31ec/OcYkQHOxbGnyCCfwKl5W2c88EqwiXUSueaTN4UZZIVF4qrTQSopIUpKeJj+px7eEgqxgmwWPbv1DK17InvfErVEHg9MXJzc20WSCWEYPz48aIJJQWOKkAMA/h60erPasBHq8FtCo1AUBJKf2Fxlk7Rgdh1rMplxynEskZxH3y9tUYr/PU9YvHZ9v8AABqTd7wnTdIuqH0mwlmUepfwsnU/99xziIuL4/yLj4/Hc889J7WssrBr3kC5RTDCNqYH+KhRQXP84dTWViZohk43iEGPuFC5RXAYhgGelaEYLd/rjupL6sZM4TX5+IJJBNh3d/cxftbQAQmAvmCwUi9+XkpQ3759RWmjRpQUeWV5O8WGqq80hNhwRWioEaWMEZ/f1lNuEUT/Lbj6k+M3p89F9WPrFAaaPDO4fNFMlSklXA9ipKSwBdchKmHMU75HLcWI5fWSGh0oixzOI86VzzDuF0qthAHR03D1b66EgZ/iGKZTrM4429NLQI8SxjuqBKkM05vH38cL6+4rlEsUhxHjIXBL7wTrtws6spjRta11sUe14Ojg2KtDhKj9iY1S5KBIw6cz8s2+8xnrhnpAaDzDMIodnqkSpCIsx0+GASKD1FO41oAYD4JZA9hzc3jcQ8bG8UayJNZ0d2b2b7ku+FwKbYJ16NQmSDqBKG4LW+6p7PgwAMCfpXpfUj4P/qdHd4Wft8y+nS4YN9kUwtorjdLv2A40nlhFKFWTpigTj1MIbcD1Rp6bGI5cEUrR0Ckuz+CaztFIiAjApn8rbU6HBTbnDtJy3IRKu16kHioYhn36cPOhSon3bB+qBKkM04tVaTeS2ES08sHCcRlyi6Fs3PwaoFCUxKAu0QD4v2CInflcrXD9Xkrw6aTTYTZ4ZrRn53UQE3uVl9l4eWx39OTw8bAkJyHMAakoSkSN55Ja3SimGByo+Tzk3f1lFuAuW6OE6GuqBNlgfE57uUWwidoGXinlTWotfw0ataCGy8YTHgwUigHD2Dg6o628gkgEw7Df03Hh/q4XxgKqBKmA6X0T5RZBVGwpQ0+M7Ow6QVTAx7fmyS2CKtnwYD+z60wpObXenZSFtBjqiO3p9OnIbuHOE8E/zRGkzq7OQLkz91QJoqiOOYNT5RbBZYjhtOtq7uWI3HMlsaEtb5gj0mOwuCRbRmlaiA7WwVsFBY8p9nHUWtk60BfLp+ayrstoH+KERI4jvWM0o1jrLr0b1QBj9p9bML3A3LqV0T4EE3hOP0a0cs/Qb523BlnxYvrDMMZ+TZF6LFLideru9dIoroXLx4UPm+cUc67rGCVPAtw7+nWQtH+aMZpCsSAiwFyRiQ3150x0ZyDIzzQdvfswv3kKMMTfB29MyBS9/3uK9ZaZ6OacUlIX/BRD3xBTZ6H6D0VJWDpLK2EsSzGpPuDjQkulViv/zUmVIDWghLtEJPQ3mGMXvqkSBLjP2/2k/HiX7MfbyzW/l0YBYa+upEvbYPj50GwjFGUSFuAjqH2HSAmCTDjyBCmh9qPilKC6ujqUlJQgLy8PWVlZWL16NWu7b775BgMHDkRRUREyMzOxcOFCF0vqeiwvIfkvH3FIat1SvM9djkmJuMr03L9TlNN9iCGrq5TkV8dloG2IMhyvKfLx7qQecougWBjA6gF2e2ESNBb3qByuDop7fSktLQUhBJs3b8aBAweQl5eHvXv3IirKfGCdNWsWPv30U2RkZODkyZPo3Lkz4uLiMHLkSJfJuvOxAS7blwE3MX6IAgMgPMAHxy/UyS2KonG00ONdRY75CTh7idJLnKIGLO+qgZ2jOdv2TAqHtxe3zcHdr3m22mGNhFgpQXKgKEtQU1MTFi1ahKlTpwIAkpOTkZGRgRUrVli1veOOO5CRoc8mHB0djX79+nFajaQixF+YmZHS8kA2vfYdfWP39dLglwf6iSGWR8FXJQr283abKUcKRUz43BZzh3Qyfn59QiZnYkAvDYOM9qFiiWYFIQTZ8bb7l9pKrGGsfREJ0S83RY7hRlFK0KFDh1BZWYnU1JYQ6LS0NGzbts2q7axZs8y+19XVoXXr1pLLqBTSVFb0ccH16WjdXNAzJyEUqdHiyO9Fw40BAHcXd+Rc50woryNQvYlCAabxzO+m89ZixS0tIfPLbs4RXZZPZ/TEOzyn66TQh9ii6RqbiJWTuBxDh6Kmw06dOgUACA4ONi4LCQnBnj17bG538eJFbN26FW+88Yak8lEcJz+pJd/N0K4xVm8AlvTmWS4DACbmKjuztyvwZTG1O6uMjOzuWPZaJYS9GshVYQkOivKZ0jMewX7SlHwoSBb3Zd5wO6bIFH7PhZeWgY+NKUJXIb8ELFia4O2F9D744IN47LHHEBcXZ7NP07/S0lIxRKVICN+H6VPX0hpvpsQE62yul9tQs+2R/i7b17icdrIfL8X96NQmCDEqcIb/vzt6yS0CAPayGbMHpODGPO5nNhtVGz5E+XPDjH9ioChLUGRkJACgqqqK9TMb7777Lry9vXHHHXfY7Fvq3ChyoEZ/Db6n4a0bM+FvEXbctW0wqusaJJDKM1DKLcAnAoTt0o4N9cN/52slkIhCcU/CLcLjtRoGjU3cA8Gs/sn44e+TkshiuVcfL43dGQFLQnpPREjvicbvYihCirIEJSUlISwsDPv37zcu27NnD7Kz2VPef/nll1i3bh1effVVAMDBgwddIqdcXJPGHX2gJvjoboE6b6v54keHpSElOpBz+1kDkqFVoWLoatSgPAtR2GxNSyhF8aNQxKJ9GP+io5YWGHupHKQcGvjci3LcropSgjQaDaZNm4bFixcD0Cs1O3fuxMSJE7F3714UFxejsbERAPDLL79g4cKFeP3111FTU4Pq6mo8+eSTcoovGYYLo2tssM12asLZB7GWYRAfEWC27LbCJFU4Sm94UHhEW6g/x4Ne+fqMYPx8tA5tNzw9hnW5GpQ+CoUvkQICFkyvfVu3ga00GlL67SjhJUVxTwxDnqC8vDyMHz8eq1atQnR0NC5cuIB9+/bh6tWrAIDx48ejrKwMERERCAwMRGBgIMrLy2WWXlrccUqPL5aOv34+Wnx7V2+ZpHEO0+KefNk6l92PpnUrX8wekIzWLFNM3PV6lF02gyvy0Z7YQ7q4h6WUQpGC2FB/fDI932ab4k6RgsaHTJkKvoqJonyCAECn02Hp0qVWy/Py8lBRUWH8bvrZk1D7O62jj99d8wZaLfOkN3wuCxfDMLizqAMYhsEDn//pYqlauKuoA177+R/Z9g8Ag7u2sVoWqPMWXDaAQlEKIf7eqLp81WyZo8OeVsMg2k7QROcY/rMNDGNdB427LcMraasc7/mKswRRuOFr2lQ6NkXnODCdt2NTJGpkRkGSoPb2lEHDm52U18zsgSnSdc6D2FB2X4cBaVF4ZGgn1nVsSJGjhULhon8n7qAfgNsCzBcGzll+HxuWZrNvGz7WVih1IoMqQSrDcB0N6doGj9q4QCnqYFqfBKtlDw1OZWnJH0+ykEUF6acBF03J4mwj5PcQO0cLhWKLRVPYg34MeDvp48gwzjsbvzkhk3W5hmF4K1hKHpGoEqQC2C4gXy+N25j5lXyDSI2UjtyGZ79BCTCMV1K/kGk1GnQzceJPd8Khvy9VSigUM9iyLzsLly7DMEDPDuHsK2F7LHn+um4OyfLRrXkObecoVAlSKWp92zdOzcgshxLopLLSJ3xJiAjAZzN6OtVHh8hWAIAPeExP3ZjXntUx3JL/3a1OR3oKxVH8vZ13+/XScPkj2t5Ow8NfyFqJIshL5Fa6pIAqQRSXo1YFTmyGdbN25JUCOaIKnQmrZRhrH7Dv7+nD2nbF1FzM6p+McB5KkBCnTwrFLWDglOmXEO57mWEYm34+XhpzH1ZCCHrZsCrJBVWCFMp0nsX3KBShuFwn4lB6hShKXFazjlGBqsgNRaE4w7zhJv6fAt4hpXzf5NKv2oX5oXNMEHTe+vvSNMeZs2OPFD57dPRQKLMGJNtc7y62lJ5Jwt4MfL00vMMy1Uz32BD0S7EdOWKPEelt0aWtcqbcTAvdlg5Pw67HrNMe2MKT82RRPJubelkHUFiWxDBlbFY7AIC3RoPeHfkXowbMFRVbShTDgFWrYcDg1XEZ6NlcBHvr3P5GB20lTgIoLk8Qxb0JC/AFIcThKbEnRnYxM7O6K11jg9EVzk3ftA70RWyIv2JCU03PWwnLoM7FUI5pQ3cJDKBQbGGZKNYAnyHUz0eLt27sIbJELbANLQxjfq8bLLVZ8WHIaB+C3/6pbNneogOaJ4hiF7U//jfPKXJq+sLHS8PL4Y6i54Ux6RjY2TyTshLfxgDgjn7s+ZHe4AjR3fJwsZTiUCiKgC1RrB7uG9mZe5xt289vsw504IpS41Jk4sP9kZvAbfl/Z5K1spaXGGbeN+fWjkOVIIUxf1QXq2WWVpPxOe0V+yCzh0EBUqv8aqOVr5co04cxdjLNisH91wjLj0R9gSiegGmQQKi/t+QvwmxKTOcY62l1y+KsVuvNPtuXuh1LOaHIQB0GS1wOh44iCqMo1doPJC7c/OLgMo+6A3kJYfYbuQFcxT5dgSOD6Jwh/LMuO8tdRR1dti8KRS1ktg/BGxPZraKupk/HCMSHB6CJzSeoeYCxXONIFXmGgaTTeQBVglSJOzmIJkQE4E2TG7tjVKCM0riO18ZnSNb3Pf0tlQh1XS9cJVIy2ofitkJhJUUoFHfhi9t7wUerQStf17vyWlrux2S1Q0LrAKNi08fE+Zrt8cTHMdp03bssU2P6vsUfy6gSpFLcZTbJx0tjpfg8fW1XmaRxLaMz22JCTnv7DQUypjkyhAuxh5GSnvE21wup3WULPx8tEiICAOgLtlIonkaPuFC8zaEgWGKr7pctTAud2nzOEKBjVCthffMcfAx+jK54zlEliKI4JuSKrxgokTbBfmgXZj0PrgZM/YxKR3S22TY7Xj/FOTEvzrl9Mi2eBXIXbKVQ5IBhGN71xIQUneZSTrh0FoPV5tVxGUa5LNfx2q/JHtgtSOadmWaTfn2CONZ0qgSpEHVNbrBDs0armz85I1a4SXZyqvPxkZ1xbUZbp/qgUDwBKYfXnIQwuxYdAot8Q3b6NKS7sGxn+f32wiTseeIaLByfgWHdxPGrpEqQSqFKBEUoKdF6JcTZytQAECCDX4LOWwsvrcYqbJZC8TSm9eGfZ4sPXI8Txvi//lNEK188ZYhgFvgIYvPniQ/XT29rNYyZVchKAKOcDPx9vDBCxMASqgQpDHdyeqYoi8Ul2QD000pq5qNb8+UWgUKRlWl9xC2rZPrYiQrSYXI++9T1nMGpNoNXfn2gHwAW/YhjzDGdsiPEupnGBWMVVYIUzAODrP0e7BWto3guXM7CYlh++DIgLUrwNt/cSau7UyhKIjxAX5A4UNdS98tby+C6HrEAgIz2IWZBHQZVxeDj6Pgjylzp6cKSn0hsqBKkYG4v5I6AUfnLPEUCZg9MQV+WAoNxzSbnlu/CnbHb83Tgfm9yluC+u8bS6u4Uiils97EZToz/t/Rmn0p7f4r5vZseG4zRmew+eLGh/ujVQWBNMkGtgcPPDBFUXsdRqBJEkY3RmW0RGSh9JmJP4oObc+y2mV6QhPhmRWj7I/159ZveLoR1+cNDhGV55mJGAc3/Q6HwxonZgFCWmnvxEf4o7mRhxWUY53xPBThGA/rpMNPk9q7ye6VKkAoZn9Pe6QrjSiCilS9oGTBpeG9yFjLbsysupgTqvNEuzM+hffTpGIFb+4qjvDw02FqZ4qugUSieiKWO8MIN6WbfWR2NOTCd9nIWU7ECdS0BFLbcOAamRcHXWyPLDAetIq9C0lwwT+oKfL00LvVX8ST4+uYIGSgtWT411+FtuZhqYqoPb+Urev8UihoghOD6Zv8b7jbm33Xe0o2lQsaJ8Fa+ZlXkASDIzxveWm4N56lru+LEhVpeNcbEhj6BKLJxT3FHTLGTbZgiHZZWuCdZive6mkcdzHJLobgbfOoLst2zbPUnhcAWpSWET6bnGx2kDd0snpKFbrEtlmm2oqhyBfxQJUhheFL+Hy+tRpQK5xTH8NLob//YEH/kJIThRo6MzrP6J4u2TylrplEongTDgPWe1afC4D+uxgTb98u0ZaEJ0plPKBnGdC8tY/Ql9LKw+LMVRbVXW0wqqBKkMNqGsPtnpEQFIqKVtUMbheIoC8d3h7dGg7uLO2LB9d0425kWZBXiwNw5JghDupq/8fF5u6VQPJ1nRndFboLtpKCWlhPT713aBiEuzDwqlIsNDxaZfWdTRGxNh706LoN1mwBfL6sp8zAWp2yzfXMsn+RkyR1bUJ8gBeLrZa6bDu3WBv1SIlFeWYNLdQ0ySUVRM4tLWsJfQ/y8cV1mLAZ1aQMAyE8K59rMisFdovH2+n95tf3f3X3sthmXbbvYK4XiicSGOldTcGIut9JgmpB30eQsaCys8YQAwX7e6NrW3PdUDH8dW/XMbCUKnj+qC5ZvLnd6/2xQJUgFBOm8EaTzRnSwDrmJ/B9YFIqBotQWR+nQAB/cVdyRtZ0th+pR3W1bcZIFVpQGgGev47ZAUSgUbkytL/4+WhSk2MktZIKvlwb1DU3QcMwFJbYOcLm/pt4XiTpGUygUGbGV7PCVcRnw0jLoyWE5+nFmX6nEolAoNmgX6o8gnmHuDMNgau8E9BOgNHlrNHh7UibrOmciTIXQqY00UdHUEqRAvr6zN3xo6DhFAvj6lf0wk30qq3NMMGdCRltvcfbCfSkUCn9CA3zwxgS9UrL9kf7w9uL/vLBXn7IguTXaBJv7pmo0jJk1WQqCdN7GshxsfH+P/el1R6BKkAIxVPumUMTm94f5JSBMjeZ+67KM9OCDZSI3CoXiON5aDXKaHacdyadla9ZpWl9xi7PyJdjfG/cOEC8SlS/U3ECheBA0JQGF4r7MYcm8bgnDMGDAQKvRiJackAH/8HZDpXkDT1/bVRQZHIUh9mxjKsdgonfzw6RQKBQKxS5XGppAQNDUpHeQtowOE8rlKw2orL6CEH9vUctv8EGM57siLUF1dXUoKSlBXl4esrKysHr1as62L774Inr06IEePXpgwYIFLpRSHZSWlsotgmzQY/dM6LF7JvTY+eHjpYGvlxZ+PlqnFSAA8PfxQrswf5crQGKhSEvQQw89hBMnTmDZsmU4cOAA8vLysHfvXkRFmTtm/fDDD5g5cyZ27twJAOjevTtefPFFDB061NjG0y1BDMPQY/dA6LHTY/c06LF73rG7pSWoqakJixYtwtSpUwEAycnJyMjIwIoVK6zavvPOO5gwYQJ0Oh10Oh0mTpyIt99+29UiUygUCoVCUSGKU4IOHTqEyspKpKa2OHilpaVh27ZtVm23bt3Kqx2FQqFQKBSKJYoLkT916hQAIDg42LgsJCQEe/bsYW1r2e706dOs/XpSYVJL6LF7JvTYPRN67J6JJx+7MyjOEmTA8oRyzfnRE0+hUCgUCsURFGcJioyMBABUVVWxfrZsW1VVZfxeVVWF1q3NU4F7orMYhUKhUCgU+yjOEpSUlISwsDDs37/fuGzPnj3Izs62apudnc2rHYVCoVAoFIolilOCNBoNpk2bhsWLFwMADh48iJ07d2LixInYu3cviouL0djYCACYMWMGVq1ahbq6OmzYsAHz589Hhw4dbPb/4YcfokePHsjKysLs2bPdwlK0detWdOjQAUuXLuVsU1ZWhtTUVBQWFhr/PvvsM9cJKTJXr17FK6+8gsLCQhQUFCA/Px9r167lbO9O513IsbvbeQeA1157DcXFxSgqKkLXrl3x2muvcbZ1p/MO8D92dzzvBg4ePAhvb2+UlZVxtnHX/HH2jn3p0qXo3r272XnfsGGDa4UUmZKSErPjKSwsRHV1NWd7wfc8USC1tbVkypQpJDc3l/To0YP8+OOPhBBCNm3aRGJiYkhtba2x7QsvvEASExNJaGgoadu2LVmyZAlnv7t37yZRUVHkzJkzpLGxkRQWFpLXX39d6sORlC+++IKMHz+eZGZm2jz2devW2VyvNg4fPkzi4+NJVVUVIYSQ1atXk4CAAPLff/9ZtXW38y7k2N3tvBNCSLdu3cipU6cIIYQcOHCAaLVa8vvvv1u1c7fzTgj/Y3fH825g3LhxxN/fn6xbt451/ffff09SUlJIbW0tqa2tJSkpKeTbb791rZASYe/YlyxZwrlOrUyZMoV3W0fuecVZggBAp9Nh6dKl2Lx5M7Zt24aBAwcCAPLy8lBRUQGdTmdsO3v2bKxfvx7nzp2zawV6//33MWTIEERERECj0eDmm29WfV6h7OxsrFy5EoGBnlV0NTAwEE888YQxOnDAgAHQ6XTYuHGjVVt3O+9Cjt0dWb58udFHsGPHjggNDcWRI0es2rnbeQf4H7u7snXrVrRq1crK99MUd80fx+fYPR1H7nlFKkFCiY2N5dWOLa/Q33//jdraWqlEkxy+xw4A//d//4d+/fqhT58+ePLJJ9HQ0CChZNISHh6OSZMmGb8TQnDlyhXWAcLdzruQYwfc67wDQLdu3YyfP//8c7Rq1cr4omSKu513gP+xA+533gFg3rx5mDdvns027po/js+xA3olsKCgAIWFhXjzzTddIJn0zJgxA3369MHgwYOxbt06znaO3POKiw6TEra8QoQQnD17Fu3atZNRMukJDg5Gfn4+Zs+ejbq6OgwbNgznz5/Hiy++KLdoorB+/XrExcWhb9++Vuvc/bzbOnZ3Pe+7d+/G2LFjUVNTg08++QQhISFWbdz1vPM5dnc8799//z06d+5s98VPSP44tcD32KOiojBw4ECUlJTgzJkzKCwsRGNjI+666y4XSSo+nTp1QlFREbKzs7F161b069cPGzZsQPfu3a3aOnTPOzxRp0AKCgpszoN36NCBvPPOO8bv//zzDwFAysvLXSCdtNg7dku+++474u/vT5qamqQTykXU1taS/Px8sm3bNtb17nze7R27Je503gkh5I8//iBRUVGsfjHufN4JsX3slqj9vDc2NpKCggJSWVlJCCEkLi6O0/fFy8vL6EdKCCFr1qwhGo3GFWJKgpBjt+TNN98kqampEkrnesaMGUNmzJjBus6Re94tpsP4wpZXiGEYj5xjbd++PS5fvowzZ87ILYpTEEJw6623YtasWejRowdrG3c973yO3RJ3Oe8GMjIyMHToUCxcuNBqnbuedwO2jt0StZ/3lStX4pprrkFYWJjdtnzyx6kJIcduSfv27VFeXi6BVPJh65gcuec9SgliyyvUuXNn+Pn5ySiVa1i4cCHq6uqM30+dOgUfHx+Eh4fLKJXzzJ49Gzk5ObjhhhtQX1+Po0ePWrVx1/PO59jd7bxXVlbi888/N1sWEBCAmpoaq7budt6FHLu7nfdff/0V3377rTFE+uTJk5g5cyZGjhxp1dbd8scJOfbnn3/e7PupU6cQExPjKlElQcgxOXTPi2ypkhXLKaEzZ86Q3r17G82Iu3fvJtHR0eTs2bOksbGRFBUVqT5k1oC9Y58yZQpZtmwZIYSQhoYGMmrUKFJSUiKHqKLx7LPPkrvuuotcunSJXLp0ifz9999k3rx5HnHe+R67u533w4cPk/T0dFJTU0MIIeT06dOkXbt25LXXXnP78y7k2N3tvFtiOiW0Z88eUlRURBoaGggh+hD5Tp06GUPkU1NT3SZEnhDbx15YWEjKysoIIYTU1NSQnJwcUlpaKpeoohAVFWVMC3Ho0CESFBRkPEYx7nm3UIK2bdtGCgoKSHBwMElJSSHXXnstIYSQo0ePkoiICHLs2DFj2xUrVpDMzEySlZVF7r33XtXOkRvge+wbNmwggwcPJoWFiQi5ZwAACvBJREFUhSQrK4tMmzaNXLhwQU7RnWL//v0EgNXfvHnz3P68Czl2dzvvtbW15JFHHiF5eXmkb9++pFu3bmTevHmksbHR7c+7kGN3t/Nu4PfffycFBQXE19eXpKenk8cff5wzf1xmZibJzMwkzz//vIwSiwefY//mm29IUVGR8bw/+OCDpL6+XmbJnWPBggWkZ8+epG/fvqRHjx7kgw8+MK4T455nCFF5ClUKhUKhUCgUB/AonyAKhUKhUCgUA1QJolAoFAqF4pFQJYhCoVAoFIpHQpUgCoVCoVAoHglVgigUCoVCoXgkVAmiUCgUCoXikVAliEKhUCgUikdClSAKhUKhUCgeCVWCKBQKhUKheCRUCaJQVMywYcOg1WqRmpqKd99917h8zJgxmDx5svH7119/jbS0NKSlpWH9+vWSyfPUU08hLi4Ot9xyi2T74MO4ceMQHByMFStW8N6md+/eSEhIQGFhIfbs2ePQfi9fvoyYmBicPHnSoe3F4tlnn0VhYSF8fX0F/QYUiqdBlSAKRcV8++236N69O4YPH45bb70VANDU1ISff/4Za9euNbYbMWIERo8ejZUrV6KgoEAyeebOnYtJkyZJ1j8ba9asQYcOHcyWffTRR+jatavgvqZOnYqysjKkpaU5JIu/vz92796N6Ohoh7YXi4ceeghlZWVo3bq1rHJQKEqHKkEUisoZMGAA1qxZY/y+Y8cOFBYW4vjx4/j777+Ny//880+kp6fLIaJHER4eLrcIFAqFJ1QJolBUzoABA7Br1y6cOXMGgN4yMn36dCQkJOCnn34CAFy4cAEhISFgGAYAsHnzZgwaNAgDBgxAQUEBHn74YRBCUF1djczMTDAMg4KCAtTW1mL37t3o1q0bEhMTsWnTJjQ1NeHpp59Gbm4uCgsLMXLkSFRUVHDKV1VVhZKSEvTs2RO9e/fGzJkzUV9fDwCYOXMmoqKicO+992Ly5MnIz89Hfn4+jh49aty+oqICgwYNQo8ePTB48GDMnz8ffn5+KC4uxqZNmzB79mxUVFSgsLAQhYWFOHfunHHbw4cPY8yYMcjMzMQ111yDCxcu8P5dd+7cidzcXHh5eWHFihUYNGgQkpOT8eGHH2LNmjUYNWoUkpOT8cILLxi3GTZsmNU03MaNG1FYWIiCggLk5+dj1qxZuHTpEpYuXYqUlBQUFxfjkUceQXFxMXx8fLBhwwbU19fj/vvvR8+ePdGnTx+MHTsWp06dAgCz7ebOnYvCwkJ06tTJTBGmUCg8Eb/wPYVCcSV1dXXEz8+PrFq1ihBCyPDhw0ltbS259dZbyZAhQwghhHz55Zdk8eLFxm1Wr15N/vjjD+P38ePHk5UrVxJCCKmpqSFBQUHkp59+Mq5/8MEHyYYNGwghhCxYsICkpqaSixcvEkIIeemll0hubq6x7dy5c8nUqVON34cOHUpuvPFG0tTURBobG8m1/9/evYVEtYVxAP9P3ibTyaTU0CahKMtQJ8UmHbygg9IU5oNFKRlZoWJEQRbhSxTdCCoSrbREpXooox6iB5HKe3R5yN6kUrPCW1je0jH/5yFcuI+nczwcOXaO3+9pZt32WutlPvb69p7kZObm5qr61NRUBgYGsq+vjySZnp7O/fv3q3qLxcLMzEyS5OjoKDds2MBly5ap+srKSs33cZGRkYyPj+fIyAjHxsYYExPD8+fP/3QfIyMjefz4cU1Zc3MzdTody8rKSJL19fU0GAy8cOECSbK1tZXOzs7s6urSjFNeXq7qXV1d+fjxY5Lk0NAQg4KC2NTURJIsKiri3Llz+eLFC5LkpUuX+OrVK+7bt4+JiYm02+0kySNHjtBsNnNsbEz1c3V15evXr0mSJSUlDA4OnrQmX19fNRchxGRyJ0iI/zgXFxdYLBZUVlaqOyx6vR7x8fGorq6G3W5HVVUV4uPjVZ81a9YgPz8fERERiImJQV1dHRobGwH8yGtJTU1VidZ2ux3Pnj1DZGQkAODy5cvYuXMn3N3dAQDp6el4+vQpmpqaJs2tvb0dDx48QE5ODnQ6HebMmYO0tDQUFxdr2iUkJMDNzQ0AEBYWhubmZgBAS0sLamtrVaK1g4PD38o5SkpKgpOTE3Q6HUJDQ9W4fwdJJCcnAwBMJhO+fv2KuLg4AIDRaISHhwfevXv3h33Lysrg7++v8rD0ej1KS0uxZMkS1WbVqlVYu3YtACAnJwerV6/G1atXkZGRAUdHRwBAdnY2Ghsb8fLlS9UvICAAgYGBALR7JoSYOseZnoAQ4p+zWq3Iz89HbW2tClbi4uIwODiIhoYGvH//XvPDu337dixatAjV1dVwdHREXl6e5ommPXv2wGw2o6urCzU1Ndi0aZOqa2trQ2lpKR4+fKjKli5dio6OjknJyOPHWgcOHICzszMA4Nu3b3B3d0d/f78KfBYsWKD66PV6FcyNH7NNTPBduHDhlPfFw8NDM25vb++U+45zcHBQ8xwPSiaO6+TkpOb7e21tbfD29taUhYSEaL5PXDsAdHR0YHh4GF5eXqpsfIy2tjaEhoZO6jdxz4QQUydBkBD/A1arFbm5uSgoKMDRo0cBAJ6enjCZTLh+/bomAAKAhoYGFBYWqh/1kZERTb3JZEJQUBBKSkrQ0NCAa9euqTqj0YicnBxkZ2erst7eXsybN2/SvIxGIwCguLhY88RVZ2enCiz+jK+vr2o/PlZ3d/df9vtVGI1G1NXVacrevHmD+fPn/zSY8/b2houLi8oBAqA+j++BEGJ6yHGYEP8DwcHB6s6OyWRS5VarFeXl5bBarZr2K1euRHV1NQBgeHhYJVBPtHfvXly8eBEGgwGenp6qPDMzE+Xl5RgcHAQAfPr0CeHh4RgaGpo0hp+fH2w2m+YdRlVVVdi6deuU1uXv7w+LxYKioiIAwPfv33Hr1i1NG4PBgP7+fgDAuXPncO/evSmN/W/YsWMHWlpa1LuZBgYGkJycrPbujzg4OGD37t0oKirC6OgoAKCgoADh4eHq2EwIMU1mOilJCDE9tm3bxpSUFE1ZVVUVHR0d+eXLF0358+fPGRISwrCwMG7ZsoWbN2+mj48Pjx07ptr09fXR3d2dT5480fQdHR3lqVOnuG7dOsbGxjI2NpY1NTUkyRMnTtBoNNLHx4eHDh0iSX7+/Jm7du3i+vXrGRsby6SkJH748IHkjyRqLy8v+vv788qVK7x79y5XrFhBDw8PZmVlkSTb29uZkJBAk8lEm83Gs2fPcvny5Wo+drudVquVZrOZFouF3d3dzMjIoMFgYEBAAO/fv8+CggI1r4lrnOj3idHNzc0MDw8nAEZHR7Onp4dRUVEEQLPZzNbWVtpsNjo7OzMkJIT19fW02Wzqujdu3CBJ1tXVMSoqilFRUYyIiGBFRQVJ8vbt22qt0dHR7OzsVNceGhriwYMH1ZpSUlL48ePHSf2ysrL49u1bzTyHh4fVOJIYLcSf05HkDMdhQgjxUz09PZp379y8eROFhYWoqamZ1utYLBYkJiYiLy9vWsedSX5+fjh9+jTS0tJmeipC/JLkOEwI8Us7fPgwHj16BODHk2plZWXYuHHjtF9n8eLFqKio+Ed/m/GrGP/bDL1eP6XcKyFmK7kTJIT4pd25cwdnzpyBm5sbBgYGEB0djZMnT8LJyWmmpyaE+I+TIEgIIYQQs5IchwkhhBBiVpIgSAghhBCzkgRBQgghhJiVJAgSQgghxKwkQZAQQgghZiUJgoQQQggxK/0Geu5fAoh9lqoAAAAASUVORK5CYII=\n", "text/plain": [ "" ] }, "execution_count": 17, "metadata": {}, "output_type": "execute_result" } ], "source": [ "Image('model_spectrum.png')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Calculating synthetic photometry" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "read_model = species.ReadModel(model='drift-phoenix',\n", " filter_name='Paranal/NACO.Mp')" ] }, { "cell_type": "code", "execution_count": 19, "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[0]:.2e}')" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "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 }