{ "cells": [ { "cell_type": "markdown", "id": "VSC-69c2a767", "metadata": { "language": "markdown" }, "source": [ "# Transmission Spectrum Basics\n", "\n", "The atmospheric components from [Temperature, pressure, and chemistry](https://escience-taurex.github.io/taurex3/examples/02_atmosphere_components.html) are now ready to be turned into a spectrum. This notebook assembles the first complete forward model in the series: a transmission geometry, where starlight grazes the planet's limb and the wavelength-dependent opacity of the atmosphere stamps its signature onto the light. [Emission spectrum basics](https://escience-taurex.github.io/taurex3/examples/04_emission_basics.html) works through the same steps for emission, and [Cloud models](https://escience-taurex.github.io/taurex3/examples/05_clouds.html) and [Contribution analysis](https://escience-taurex.github.io/taurex3/examples/07_contribution_analysis.html) extend this setup further.\n", "\n", "In transmission, the effective area blocked by the planet at each wavelength depends on how high in the atmosphere the gas becomes optically thick. Wavelengths that overlap with strong molecular absorption bands probe higher, cooler layers and show a larger transit depth.\n", "\n", "This notebook builds a `TransmissionModel` step by step and then shows how adding radiative contributions — molecular absorption first, then CIA and Rayleigh scattering — reshapes the spectrum. More information about forward models is [here](../user/taurex/models.rst), planet parameters are [here](../user/taurex/planet.rst), and star parameters are [here](../user/taurex/star.rst)." ] }, { "cell_type": "markdown", "id": "63ece2a7", "metadata": {}, "source": [ "## Data Note\n", "\n", "This notebook uses the opacity and CIA files set up in [Setup and opacity data](https://escience-taurex.github.io/taurex3/examples/01_setup_and_data.html). TauREx provides the software to work with these datasets; the files themselves are third-party products from [ExoMol](https://www.exomol.com) and [HITRAN](https://hitran.org/)." ] }, { "cell_type": "code", "execution_count": 1, "id": "VSC-c765f003", "metadata": { "language": "python" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Planet: 0.74 M_jup, 1.38 R_jup\n", "Star: 6117 K, 1.16 R_sun\n" ] } ], "source": [ "from _shared import build_base_components\n", "from taurex.constants import RSOL\n", "\n", "context = build_base_components(download=False)\n", "iso_t = context['iso_t']\n", "press = context['press']\n", "chemistry = context['chemistry']\n", "planet = context['planet']\n", "star = context['star']\n", "\n", "print(f'Planet: {planet.mass:.2f} M_jup, {planet.radius:.2f} R_jup')\n", "print(f'Star: {star.temperature:.0f} K, {star.radius / RSOL:.2f} R_sun')\n" ] }, { "cell_type": "markdown", "id": "VSC-72da9cfa", "metadata": { "language": "markdown" }, "source": [ "## Building the Model\n", "\n", "Assembling a `TransmissionModel` connects the atmosphere from [Temperature, pressure, and chemistry](https://escience-taurex.github.io/taurex3/examples/02_atmosphere_components.html) with a planet and a host star. Each physical process — molecular absorption, CIA, Rayleigh scattering, clouds — is added as an independent **contribution** object, which keeps the code composable and easy to extend.\n", "\n", "The model must be explicitly `build()`-ed before running. Calling `build()` again after adding or removing a contribution refreshes the radiative-transfer kernel. More information about forward models is [here](../user/taurex/models.rst), planet parameters are [here](../user/taurex/planet.rst), and star parameters are [here](../user/taurex/star.rst)." ] }, { "cell_type": "code", "execution_count": 2, "id": "VSC-8bfee465", "metadata": { "language": "python" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Contributions: ['Absorption']\n" ] } ], "source": [ "from taurex.model import TransmissionModel\n", "from taurex.contributions import AbsorptionContribution\n", "\n", "tm = TransmissionModel(\n", " planet=planet,\n", " temperature_profile=iso_t,\n", " chemistry=chemistry,\n", " star=star,\n", " pressure_profile=press,\n", ")\n", "tm.add_contribution(AbsorptionContribution())\n", "tm.build()\n", "\n", "print('Contributions:', [c.name for c in tm.contribution_list])" ] }, { "cell_type": "code", "execution_count": 3, "id": "VSC-08952604", "metadata": { "language": "python" }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Computed 76744 spectral points.\n", "Wavelength range: 0.300 to 50.002 um\n" ] } ], "source": [ "wngrid, rprs, tau, _ = tm.model()\n", "wlgrid = 10000 / wngrid[::-1]\n", "abs_only_rprs = rprs[::-1]\n", "\n", "print(f'Computed {len(wngrid)} spectral points.')\n", "print(f'Wavelength range: {wlgrid.min():.3f} to {wlgrid.max():.3f} um')" ] }, { "cell_type": "code", "execution_count": 4, "id": "VSC-047d132a", "metadata": { "language": "python" }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAGOCAYAAAAHL7LOAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAgz9JREFUeJzt3XtczfcfB/DXqXS/UelOLiG3IgpDLhGaLXNnbsMwjLXZMDQzw4Yfm/vmuglzmZnrkuvIrdwpREIqia66nu/vj9bh6FSnOtd6PR+PHo/O93y+n8/7nM45vc/n+7mIBEEQQERERERVho66AyAiIiIi1WICSERERFTFMAEkIiIiqmKYABIRERFVMUwAiYiIiKoYJoBEREREVQwTQCIiIqIqhgkgERERURXDBJCIiIioimECSFQJxcTEQCQSYdOmTQqvWyQS4ZtvvlF4vVS1jBw5EqampuoOQ6FOnDgBkUiEEydOqDuUClPmZwhpBiaAVOmJRCK5firDhzapx/fff4+9e/eqOwxSkVWrVjExIq2np+4AiJTtt99+k7q9ZcsWhISEFDnu5uamyrCUqnbt2nj16hWqVaum8LpfvXoFPT1+dLzp+++/R79+/RAQEKDuUEgFVq1aBWtra4wcOVLqeMeOHfHq1Svo6+urJzCiMuCnOFV6H374odTtc+fOISQkpMjxt2VmZsLY2FiZoSmNSCSCoaGhUupWVr1VRUZGBkxMTNQdBv0nKysL+vr60NGp+AUxHR0dvj9Ia/ASMBGATp06oWnTpggPD0fHjh1hbGyMmTNnAgD++usv+Pv7w8HBAQYGBqhXrx7mzZuH/Px8mXXcunULnTt3hrGxMRwdHfHDDz8Uae/nn39GkyZNYGxsjOrVq6NVq1YIDg6W3P/NN99AJBLhzp07+PDDD2FhYQEbGxvMnj0bgiDg0aNHeP/992Fubg47OzssWbJEqn5Z43fi4+MxatQoODk5wcDAAPb29nj//fcRExMjKXPp0iX4+fnB2toaRkZGqFOnDj766COpumWNAbx8+TJ69uwJc3NzmJqaomvXrjh37pxUmU2bNkEkEuHMmTMIDAyEjY0NTExM0KdPHzx79qzUv5E88bu4uODdd9/FP//8Aw8PDxgaGqJx48bYs2dPkfpevnyJqVOnwtnZGQYGBqhfvz4WLVoEsVgsVU4sFmP58uVo1qwZDA0NYWNjgx49euDSpUuS5yMjIwObN2+WDCco7Bkq/DveunULQ4YMQfXq1dG+fXsABa+XTp06FYlr5MiRcHFxkdwu/FsuXrwYK1euRN26dWFsbIzu3bvj0aNHEAQB8+bNg5OTE4yMjPD+++8jOTm51OdTWU6fPo3+/fujVq1aMDAwgLOzMz777DO8evVKZvn79+/Dz88PJiYmcHBwwLfffgtBEKTKbN++HZ6enjAzM4O5uTmaNWuG5cuXF6mnf//+qFGjBoyNjdGmTRscOHBAqkzhGL3t27dj1qxZcHR0hLGxMVJTUyWvz1OnTmHcuHGwsrKCubk5hg8fjhcvXkjqcHFxwc2bN3Hy5EnJ37vw71jcGMCdO3fC09MTRkZGsLa2xocffognT55IlSkcE/nkyRMEBATA1NQUNjY2+OKLL4p81hRn1apVaNKkCQwMDODg4ICJEyfi5cuXUmXK8jn1po0bN0IkEuHy5ctF7vv++++hq6tb5DGRZmMPINF/nj9/jp49e2LQoEH48MMPYWtrC6AgcTE1NUVgYCBMTU1x7NgxzJkzB6mpqfjxxx+l6njx4gV69OiBDz74AAMGDMCuXbvw1VdfoVmzZujZsycA4JdffsGnn36Kfv36YcqUKcjKysK1a9dw/vx5DBkyRKq+gQMHws3NDQsXLsSBAwfw3XffoUaNGli7di26dOmCRYsWYevWrfjiiy/QunVrdOzYsdjH17dvX9y8eROTJ0+Gi4sLEhMTERISgtjYWMnt7t27w8bGBtOnT4elpSViYmJkJk9vunnzJjp06ABzc3N8+eWXqFatGtauXYtOnTrh5MmT8Pb2lio/efJkVK9eHUFBQYiJicGyZcswadIk7Nixo8R2Sou/0N27dzFw4ECMHz8eI0aMwMaNG9G/f38cPnwY3bp1A1DQu+vj44MnT55g3LhxqFWrFs6ePYsZM2bg6dOnWLZsmaS+0aNHY9OmTejZsyfGjBmDvLw8nD59GufOnUOrVq3w22+/YcyYMfDy8sLHH38MAKhXr55U7P3794erqyu+//77IsmNvLZu3YqcnBxMnjwZycnJ+OGHHzBgwAB06dIFJ06cwFdffYV79+7h559/xhdffIENGzaUq52K2rlzJzIzMzFhwgRYWVnhwoUL+Pnnn/H48WPs3LlTqmx+fj569OiBNm3a4IcffsDhw4cRFBSEvLw8fPvttwCAkJAQDB48GF27dsWiRYsAALdv38aZM2cwZcoUAEBCQgLatWuHzMxMfPrpp7CyssLmzZvx3nvvYdeuXejTp49Uu/PmzYO+vj6++OILZGdnS12ynTRpEiwtLfHNN98gKioKq1evxsOHDyXJ3bJlyzB58mSYmpri66+/BgDJZ4UsmzZtwqhRo9C6dWssWLAACQkJWL58Oc6cOYPLly/D0tJS6vnw8/ODt7c3Fi9ejKNHj2LJkiWoV68eJkyYUOLz/s0332Du3Lnw9fXFhAkTJLFfvHgRZ86ckRoOIs/n1Nv69euHiRMnYuvWrWjRooXUfVu3bkWnTp3g6OhYYoykYQSiKmbixInC2y99Hx8fAYCwZs2aIuUzMzOLHBs3bpxgbGwsZGVlFaljy5YtkmPZ2dmCnZ2d0LdvX8mx999/X2jSpEmJMQYFBQkAhI8//lhyLC8vT3BychJEIpGwcOFCyfEXL14IRkZGwogRIyTHHjx4IAAQNm7cKCkDQPjxxx+LbfPPP/8UAAgXL14sMTYAQlBQkOR2QECAoK+vL0RHR0uOxcXFCWZmZkLHjh0lxzZu3CgAEHx9fQWxWCw5/tlnnwm6urrCy5cvi21TnvgFQRBq164tABB2794tOZaSkiLY29sLLVq0kBybN2+eYGJiIty5c0fq/OnTpwu6urpCbGysIAiCcOzYMQGA8OmnnxZp683HYGJiIvX8Fyr8Ow4ePLjIfT4+PoKPj0+R4yNGjBBq164tuV34t7SxsZF6jmbMmCEAENzd3YXc3FzJ8cGDBwv6+vpSr01VkvV+WbBggSASiYSHDx9Kjo0YMUIAIEyePFlyTCwWC/7+/oK+vr7w7NkzQRAEYcqUKYK5ubmQl5dXbJtTp04VAAinT5+WHEtLSxPq1KkjuLi4CPn5+YIgCMLx48cFAELdunWLxFn4+vT09BRycnIkx3/44QcBgPDXX39JjjVp0kTm366w/uPHjwuCIAg5OTlCzZo1haZNmwqvXr2SlNu/f78AQJgzZ06R5+Pbb7+VqrNFixaCp6dnsY9dEAQhMTFR0NfXF7p37y55rIIgCCtWrBAACBs2bJAck/dz6u3PEEEoeG05ODhItREREVGkHGkHXgIm+o+BgQFGjRpV5LiRkZHk97S0NCQlJaFDhw7IzMxEZGSkVFlTU1OpsYX6+vrw8vLC/fv3JccsLS3x+PFjXLx4sdSYxowZI/ldV1cXrVq1giAIGD16tFR9DRs2lGpD1mPQ19fHiRMnpC5nvamwJ2L//v3Izc0tNTagoMfin3/+QUBAAOrWrSs5bm9vjyFDhuDff/9Famqq1Dkff/wxRCKR5HaHDh2Qn5+Phw8fVij+Qg4ODlI9PoWX8S5fvoz4+HgABb1UHTp0QPXq1ZGUlCT58fX1RX5+Pk6dOgUA2L17N0QiEYKCgoq08+ZjKM348ePlLluc/v37w8LCQnK7sGf1ww8/lJqU4+3tjZycHLVdjnvz/ZKRkYGkpCS0a9cOgiDIvHw4adIkye8ikQiTJk1CTk4Ojh49CqDgdZmRkYGQkJBi2zx48CC8vLwkl9eBgvfixx9/jJiYGNy6dUuq/IgRI6TifNPHH38s1Vs2YcIE6Onp4eDBg6U88qIuXbqExMREfPLJJ1JjA/39/dGoUaMil6iBoq+VDh06lPjeBoCjR48iJycHU6dOlRrLOHbsWJibmxdpR57PKVmGDx+OuLg4HD9+XHJs69atMDIyQt++fUs8lzQPE0Ci/zg6OsqcvXfz5k306dMHFhYWMDc3h42NjeTDMyUlRaqsk5NTkcSgevXqUknLV199BVNTU3h5ecHV1RUTJ07EmTNnZMZUq1YtqdsWFhYwNDSEtbV1keMlJUYGBgZYtGgRDh06BFtbW3Ts2BE//PCDJCECAB8fH/Tt2xdz586FtbU13n//fWzcuBHZ2dnF1vvs2TNkZmaiYcOGRe5zc3ODWCzGo0ePSnxM1atXB4AKx1+ofv36Rf4GDRo0AADJeMG7d+/i8OHDsLGxkfrx9fUFACQmJgIAoqOj4eDggBo1ahQbmzzq1KlTofMB2a8FAHB2dpZ5vKTn89WrV4iPjy/Xz9uv+bfFxsZi5MiRqFGjhmQcm4+PD4Ci7xcdHR2pLw5A0b/VJ598ggYNGqBnz55wcnLCRx99hMOHD0ud8/Dhw2Jfg4X3v6mkv4erq6vUbVNTU9jb20uNNZVXYbuyYmvUqFGRuArHmL7p7c+PsrSjr6+PunXrFmlHns8pWbp16wZ7e3ts3boVQMH42G3btuH999+HmZlZieeS5mECSPQfWT0CL1++hI+PD65evYpvv/0Wf//9N0JCQiRjkd6eMKCrqyuzbuGNcV9ubm6IiorC9u3b0b59e+zevRvt27eX2cskqz552pBl6tSpuHPnDhYsWABDQ0PMnj0bbm5ukl4ZkUiEXbt2ISwsDJMmTcKTJ0/w0UcfwdPTE+np6SXWXRbKir8sxGIxunXrhpCQEJk/iu7NkPXaKq4HsbgB/8U9b+V5Pnfs2AF7e/ty/RSOuysu9m7duuHAgQP46quvsHfvXoSEhEgmI739fpFHzZo1ceXKFezbtw/vvfcejh8/jp49e2LEiBFlrqtQcb1/6lbc31JV7ZT2HtTV1cWQIUOwe/duZGVl4fjx44iLiyt1RQXSTJwEQlSCEydO4Pnz59izZ4/UBIsHDx5UqF4TExMMHDgQAwcORE5ODj744APMnz8fM2bMUOoyEvXq1cPnn3+Ozz//HHfv3oWHhweWLFmC33//XVKmTZs2aNOmDebPn4/g4GAMHToU27dvl7ocXcjGxgbGxsaIiooqcl9kZCR0dHSK9FApO/579+5BEASpBOvOnTsAIJksUq9ePaSnp0t6/Epq78iRI0hOTi6xF7Asl4MLVa9eXeYlt5IuhSuKn59fiZdUS+Lg4FDsfdevX8edO3ewefNmDB8+XHK8uLbEYjHu378v6fUDiv6tgIKerN69e6N3794Qi8X45JNPsHbtWsyePRv169dH7dq1i30NAgXrYsrr7t276Ny5s+R2eno6nj59il69ekmOyfv3Lmw3KioKXbp0kbovKiqqTHHJ286bPao5OTl48OBBqa/zshg+fDiWLFmCv//+G4cOHYKNjQ38/PwUVj+pDnsAiUpQ+E35zW/GOTk5WLVqVbnrfP78udRtfX19NG7cGIIgyD32rqwyMzORlZUldaxevXowMzOTXOJ98eJFkR4ADw8PACj2MrCuri66d++Ov/76S+oSWUJCAoKDg9G+fXuYm5urJP5CcXFx+PPPPyW3U1NTsWXLFnh4eMDOzg4AMGDAAISFheHIkSNF2nr58iXy8vIAFMw8FgQBc+fOLVLuzefKxMSkyHIbpalXrx4iIyOllsC5evVqscMBFMne3h6+vr7l+mncuHGx9cp6vwiCUGTJljetWLFCquyKFStQrVo1dO3aFUDR94uOjg6aN28O4PXrslevXrhw4QLCwsIk5TIyMrBu3Tq4uLiUGPPb1q1bJ/U+XL16NfLy8qRmx8r7927VqhVq1qyJNWvWSL1ODx06hNu3b8Pf31/uuEri6+sLfX19/PTTT1LP/fr165GSkqKwdgCgefPmaN68OX799Vfs3r0bgwYN4sLwWop/NaIStGvXDtWrV8eIESPw6aefQiQS4bfffiv3Uh4A0L17d9jZ2eGdd96Bra0tbt++jRUrVsDf319p42ju3LmDrl27YsCAAWjcuDH09PTw559/IiEhAYMGDQIAbN68GatWrUKfPn1Qr149pKWl4ZdffoG5ublU78fbvvvuO4SEhKB9+/b45JNPoKenh7Vr1yI7O7vUtcUUGX+hBg0aYPTo0bh48SJsbW2xYcMGJCQkYOPGjZIy06ZNw759+/Duu+9i5MiR8PT0REZGBq5fv45du3YhJiYG1tbW6Ny5M4YNG4affvoJd+/eRY8ePSAWi3H69Gl07txZMoHB09MTR48exdKlS+Hg4IA6deoUWf7mbR999BGWLl0KPz8/jB49GomJiVizZg2aNGlSZOKMtmjUqBHq1auHL774Ak+ePIG5uTl2795d7NgyQ0NDHD58GCNGjIC3tzcOHTqEAwcOYObMmZKxcGPGjEFycjK6dOkCJycnPHz4ED///DM8PDwkY/ymT5+Obdu2oWfPnvj0009Ro0YNbN68GQ8ePMDu3bvLtMhzTk6O5LUWFRWFVatWoX379njvvfckZTw9PbF69Wp89913qF+/PmrWrFmkhw8AqlWrhkWLFmHUqFHw8fHB4MGDJcvAuLi44LPPPivL01ssGxsbzJgxA3PnzkWPHj3w3nvvSWJv3bq1wi/RDh8+HF988QWAogvtkxZR9bRjInUrbhmY4pZmOXPmjNCmTRvByMhIcHBwEL788kvhyJEjUss9lFTH28t6rF27VujYsaNgZWUlGBgYCPXq1ROmTZsmpKSkSMoULh9SuBTGm3WZmJgUaePttt9ewiEpKUmYOHGi0KhRI8HExESwsLAQvL29hT/++ENyTkREhDB48GChVq1agoGBgVCzZk3h3XffFS5duiTVFt5aBqbwXD8/P8HU1FQwNjYWOnfuLJw9e1aqTOEyG28vM/P20hmyyBO/IBQsA+Pv7y8cOXJEaN68uWBgYCA0atRI2LlzZ5E609LShBkzZgj169cX9PX1BWtra6Fdu3bC4sWLpZYBycvLE3788UehUaNGgr6+vmBjYyP07NlTCA8Pl5SJjIwUOnbsKBgZGQkAJEvCFPd3LPT7778LdevWFfT19QUPDw/hyJEjxS4D8/YSOIXP29uPrbjnWVVu3bol+Pr6CqampoK1tbUwduxY4erVq0WWCil8LUdHRwvdu3cXjI2NBVtbWyEoKEhqmZFdu3YJ3bt3F2rWrCno6+sLtWrVEsaNGyc8ffpUqt3o6GihX79+gqWlpWBoaCh4eXkJ+/fvlypT3HMmCK+ft5MnTwoff/yxUL16dcHU1FQYOnSo8Pz5c6my8fHxgr+/v2BmZiYAkCwJU9xreceOHUKLFi0EAwMDoUaNGsLQoUOFx48fS5Up7r1d+BqSx4oVK4RGjRoJ1apVE2xtbYUJEyYIL168kCoj7+eUrGVgCj19+lTQ1dUVGjRoIFdcpJlEglCBrgwiIg3i4uKCpk2bYv/+/eoOhbRM4YLNFy9eRKtWrdQdjkZLSkqCvb095syZg9mzZ6s7HConjgEkIiIiuW3atAn5+fkYNmyYukOhCuAYQCIiIirVsWPHcOvWLcyfPx8BAQFSM7VJ+zABJCIiolJ9++23OHv2LN555x38/PPP6g6HKkgrLgGvXLkSLi4uMDQ0hLe3Ny5cuFBi+Z07d6JRo0YwNDREs2bNpLbwyc3NlWx6bWJiAgcHB8n2Nm+aP38+2rVrB2NjY6nNuolIc8XExHD8H5XLyJEjIQgCx/+V4MSJE8jJycHx48fh6Oio7nCogjQ+AdyxYwcCAwMRFBSEiIgIuLu7w8/PT7JV09vOnj2LwYMHY/To0bh8+TICAgIQEBCAGzduAChYTywiIgKzZ89GREQE9uzZg6ioKKkp/kDBUgD9+/fHhAkTlP4YiYiIiFRJ42cBe3t7o3Xr1pLFQsViMZydnTF58mRMnz69SPmBAwciIyNDqhegTZs28PDwwJo1a2S2cfHiRXh5eeHhw4dF9tvctGkTpk6dWuZFXomIiIg0lUaPAczJyUF4eDhmzJghOaajowNfX1+pFd/fFBYWhsDAQKljfn5+2Lt3b7HtpKSkQCQSKfxSr1gsRlxcHMzMzMq1VRQRERGRvARBQFpaGhwcHEpdAF2jE8CkpCTk5+fD1tZW6ritra1kj8e3xcfHyywfHx8vs3xWVha++uorDB48uMJbVmVnZ0tt9/PkyZMybUFEREREVFGPHj2Ck5NTiWU0OgFUttzcXAwYMACCIGD16tUVrm/BggUy9wy9ceOG0rb4ospFEARkZGTAxMSEvcZEKsD3HFUmaWlpaNq0qVw5h0YngNbW1tDV1UVCQoLU8YSEBMmm7m+zs7OTq3xh8vfw4UMcO3ZMIRvWz5gxQ+ryc2pqKpydneHk5KSQ+qnyEwQBKSkpsLCw4D8jIhXge44qk8J9xOV5LWt0Aqivrw9PT0+EhoYiICAAQMG4utDQUMkm7G9r27YtQkNDMXXqVMmxkJAQtG3bVnK7MPm7e/cujh8/DisrK4XEa2BgAAMDgyLHRSIRP1hIboWvF75miFSD7zmqLMryGtboBBAAAgMDMWLECLRq1QpeXl5YtmwZMjIyMGrUKADA8OHD4ejoiAULFgAApkyZAh8fHyxZsgT+/v7Yvn07Ll26hHXr1gEoSP769euHiIgI7N+/H/n5+ZLxgTVq1IC+vj4AIDY2FsnJyYiNjUV+fj6uXLkCAKhfvz5MTU1V/CwQERERKY7GJ4ADBw7Es2fPMGfOHMTHx8PDwwOHDx+WTPSIjY2VmunSrl07BAcHY9asWZg5cyZcXV2xd+9eNG3aFEDBxIx9+/YBADw8PKTaOn78ODp16gQAmDNnDjZv3iy5r0WLFkXKEBEREWkjjV8HUJulpqbCwsICKSkpJY4BzM/PR25urgojI01VOIVfWUsHVatWDbq6ugqvl0hbcQwgVSby5h2AFvQAVmaCICA+Pp6LTJMUsViM58+fK61+S0tL2NnZ8Z8dEVEVxgRQjQqTv5o1a8LY2Jj/kAmCICA/Px+6uroKfz0IgoDMzEzJNor29vYKrZ+IiLQHE0A1yc/PlyR/ipqFTNpPmQkgABgZGQEAEhMTUbNmTV4OJiKqokreJ4SUpnDMn7GxsZojoaqm8DXHcadERFUXE0A142VfUjW+5oiIiAkgERERKUROnhjRz9LVHQbJgQkgKdyJEycgEom0anbzyJEjJbvNaDoXFxcsW7ZM3WEQkQYSiwU8S8tWS9uCIKDfmrPouuQk1v/7QOXt34pLxZBfzmHdqWiVt62NmABSuYSFhUFXVxf+/v7qDqVMYmJiIBKJJDu7FFq+fDk2bdqklpiIiBRBEAQM+uUcvL4/ip2XHqm8/ccvXuHa4xQAwLz9t1Te/oC1YTgb/RzfH4xEQmqWytvXNkwAqVzWr1+PyZMn49SpU4iLi1N3OMjJyanQ+RYWFrC0tFRMMEREahDzPBMXHiRDEIBpu66pvH2xmveVSM/Ok/yelK6eXlBtwgSQyiw9PR07duzAhAkT4O/vX2zP2ZkzZ9C8eXMYGhqiTZs2uHHjhuS+hw8fonfv3qhevTpMTEzQpEkTHDx4UHL/yZMn4eXlBQMDA9jb22P69OnIy3v95u7UqRMmTZqEqVOnwtraGn5+fgAKJjisXr0aPXv2hJGREerWrYtdu3ZJzqtTpw6Agq39RCKRZFu/ty8BZ2dn49NPP0XNmjVhaGiI9u3b4+LFi5L7Cy9zh4aGolWrVjA2Nka7du0QFRVV4nN3/fp1dOnSBUZGRrCyssLHH3+M9PTX42VGjRqFvn37YvHixbC3t4eVlRUmTpxY7Izdjz76CO+++67UsdzcXNSsWRPr168vMRYiqlxy88XqDoG0CBNAKrM//vgDjRo1QsOGDfHhhx9iw4YNkLWj4LRp07BkyRJcvHgRNjY26N27tySRmThxIrKzs3Hq1Clcv34dixYtgqmpKYCC/Zp79eqF1q1b4+rVq1i9ejXWr1+P7777Tqr+zZs3Q19fH2fOnMGaNWskx2fPno2+ffvi6tWrGDp0KAYNGoTbt28DAC5cuAAAOHr0KJ4+fYo9e/bIfIxffvkldu/ejc2bNyMiIgL169eHn58fkpOTpcp9/fXXWLJkCS5dugQ9PT189NFHxT5vGRkZ8PPzQ/Xq1XHx4kXs3LkTR48exaRJk6TKnThxAtHR0Th+/Dg2b96MTZs2FZtkjxkzBocPH8bTp08lx/bv34/MzEwMHDiw2FiIiKhq40LQGqT3z/+qZfCujZkB/p7cXu7y69evx4cffggA6NGjB1JSUnDy5ElJb1qhoKAgdOvWDUBBsubk5IQ///wTAwYMQGxsLPr27YtmzZoBAOrWrSs5b9WqVXB2dsaKFSsgEonQqFEjxMXF4auvvsKcOXOgo1PwvcXV1RU//PBDkfj69++PMWPGAADmzZuHkJAQ/Pzzz1i1ahVsbGwAAFZWVrCzs5P5+DIyMrB69Wps2rQJPXv2BAD88ssvCAkJwfr16zFt2jRJ2fnz58PHxwcAMH36dPj7+yMrKwuGhoZF6g0ODkZWVha2bNkCExMTAMCKFSvQu3dvLFq0CLa2tgCA6tWrY8WKFdDT00OjRo3g7++P0NBQjB07tkid7dq1Q8OGDfHbb7/hyy+/BABs3LgR/fv3lyTURFQ1hUU/x8JDtxHQwhGj3qmj7nBIwzAB1CDP0rIRr+EDV6OionDhwgX8+eefAAA9PT0MHDgQ69evL5IAtm3bVvJ7jRo10LBhQ0lP3KeffooJEybgn3/+ga+vL/r27YvmzZsDAG7fvo22bdtKrVf3zjvvID09HY8fP0atWrUAAJ6enjJjfLPdwttvT/ooSXR0NHJzc/HOO+9IjlWrVg1eXl6S+AsVxgy83lotMTFREuObbt++DXd3d0nyV/i4xGIxoqKiJAlg48aNpXbosLe3x/Xr14uNd8yYMVi3bh2+/PJLJCQk4NChQzh27Jjcj5eIKqfBv5wDAFx9nIIP29RGNV1e9KPXmABqEBszA41vd/369cjLy4ODg4PkmCAIMDAwwIoVK2BhYSFXPWPGjIGfnx8OHDiAf/75BwsWLMCSJUswefJkuWN5M5FSl2rVqkl+L0xYxeKKjcN5s87Cekuqc/jw4Zg+fTrCwsJw9uxZ1KlTBx06dKhQDERUuaRl5UEQBFiZquf/DGkeJoAapCyXYdUhLy8PW7ZswZIlS9C9e3ep+wICArBt2zaMHz9ecuzcuXOSnrAXL17gzp07cHNzk9zv7OyM8ePHY/z48ZgxYwZ++eUXTJ48GW5ubti9ezcEQZAkVWfOnIGZmRmcnJxKjfPcuXMYPny41O0WLVoAAPT19QEU7MVcnHr16knGFtauXRtAwcSKixcvYurUqaW2Xxw3Nzds2rQJGRkZkuT1zJkz0NHRQcOGDctdr5WVFQICArBx40aEhYVh1KhR5a6LiCqndgtDIRYDuye0QzMn+b6oazMRuONRadgfTHLbv38/Xrx4gdGjR6Np06ZSP3379i0y6/Tbb79FaGgobty4gZEjR8La2loy03bq1Kk4cuQIHjx4gIiICBw/flySHH7yySd49OgRJk+ejMjISPz1118ICgpCYGCgZPxfSXbu3IkNGzbgzp07CAoKwoULFyQTLWrWrAkjIyMcPnwYCQkJSElJKXK+iYkJJkyYgGnTpuHw4cO4desWxo4di8zMTIwePbrcz9/QoUNhaGiIESNG4MaNGzh+/DgmT56MYcOGSS7/lteYMWOwefNm3L59GyNGjKhQXURU+WTlipGTL8b438PVHYrEvqtx+F/IHaRlcV9ydWACSHJbv349fH19ZV7m7du3Ly5duoRr116vPbVw4UJMmTIFnp6eiI+Px99//y3VAzdx4kS4ubmhR48eaNCgAVatWgUAcHR0xMGDB3HhwgW4u7tj/PjxGD16NGbNmiVXnHPnzsX27dvRvHlzbNmyBdu2bUPjxo0BFIxZ/Omnn7B27Vo4ODjg/fffl1nHwoUL0bdvXwwbNgwtW7bEvXv3cOTIEVSvXr1Mz9mbjI2NceTIESQnJ6N169bo168funbtihUrVpS7zkK+vr6wt7eHn5+f1OV5Iqo65OnzSs6o2JqpihIVn4ZPt13G8tC7WHgoUt3hVEkiQdb6HaQQqampsLCwQEpKCszNzaXuy8rKwoMHD1CnTh2ZM0apfEQiEf7880+t2dbtbYIgID8/H7q6ulKTYEqTnp4OR0dHbNy4ER988EGJZfnaI3pNEASkpKTAwsKiTO85TXQ3IQ3d/neqxDJG1XRxe14PpbQf+zwTHX88Lrkds7D4naLW//tAareQksrKy2X6AcnvBz/tgMYO5iWUrpxKyjvexjGARFpMLBYjKSkJS5YsgaWlJd577z11h0RERFqACSCRFouNjUWdOnXg5OSETZs2QU+Pb2kiIiod/1tQpVLVRjS4uLhUucdMRFQaLb+arxKcBEJERFRFCFDOF8aM7Dz87+gductz32L1Yw+gmrH3hlSNrzkiUrQl/9zBn5efyF1+3an7ZW5DLBaw/eIjCBAwuHUt6Oiwm68imACqSeFuD5mZmTAyMlJzNFSVZGZmAii64wgRaTd1XvbccOZBmcqXZzma/defYuafBdtimujrIaCFY5nroNeYAKqJrq4uLC0tkZiYCKBgjThtX4KAKq68y8DIW3dmZiYSExNhaWkptd8wEZGm2/RGkrnpbAwTwApiAqhGdnZ2ACBJAomAgqVd5NnxpLwsLS0lrz0i0lx3EtLgXN0YRvr8skaKxwRQjUQiEezt7VGzZk3k5nIrHCropUtLS4OZmZlSeoSrVavGnj8iLfDbuYeYvfcG6tmYIOQzH4WOd3uRkYPqJvoVqiNfLCD1VW6F61EWXlArHRNADaCrq8t/ygSgIAHMzs6GoaEhhwQQVWGz994AAEQ/y8DNuFQ0cyq6BWd5ZOWK0WJeCL7o3gCTuriWqw5BEBCw8gxuPU3FyiEt0KOpvUJiK7VdlbRSdXAZGCIiIg2Wr4SZ+4v/kX/JlrddeJCM609SkC8WMP73CIXFJBYLGLP5InyXnsT9Z+klluX344pjAkhERERyy8pTzhp+f119gqO3E3EvMR0TFJhYZubkKayuyoQJIBERESlc7PNM/H01Dlm5+XKVf5T8SvJ7VEJakfvL0xG68vg9NA06gu8P3pZ5/6+n78Nr/lHsvPSo7JVrOSaAREREVCpBEBAW/RxPX74qvTAA3/+dxORtlzH4l3PIyxcjMycPaVnln/D4Zv53OfYlPv/jaqnJ5Y9HoiAWil94+rsDt5GYlo1pu66VOy5txUkgREREVKqd4Y/xZRkSpZz/LhVfjn2JUZsu4vTdJABA8FhvhcSzO+IxnGsYYapvA4XUBwDHIxPx27mHGN2+Dt6pb62wejURewCJiIg0mKZs31iW5O9thckfAAz55bwiwgEAnLzzTOZxEco3S2TUpos4FpmIob8qLkZNxQSQiIhIS+Xll39CRtzLV/h4yyUs+SdKgRGRtmACSEREpIWGb7iA+l8fQtsFoQCAe4klL53ytqk7ruCfWwn4+dg9nL//XBkhVsiNJynSB8rQE8plYkrHBJCIiEjLpGbl4tR/lz+fpmQhKj6tzGvyXXiQLPk9Mr7orFt5vczMKfe5JZn7902p21cfpxRTksqDCSAREZEGe5qSVeRSb3qW9Np2spZNUZXZf90svZCSaMjwSK2kNQngypUr4eLiAkNDQ3h7e+PChQsllt+5cycaNWoEQ0NDNGvWDAcPHpTcl5ubi6+++grNmjWDiYkJHBwcMHz4cMTFxUnVkZycjKFDh8Lc3ByWlpYYPXo00tPL1sVORERUEZ9sjUCfVWdLnAxS1su/ivT31bjSC2mglFflX5KmMtCKBHDHjh0IDAxEUFAQIiIi4O7uDj8/PyQmJsosf/bsWQwePBijR4/G5cuXERAQgICAANy4UbC3YmZmJiIiIjB79mxERERgz549iIqKwnvvvSdVz9ChQ3Hz5k2EhIRg//79OHXqFD7++GOlP14iIqI3XX+SUqHLtJVVcWP95BkCOHGr4nYbKfQiIwfTdl7FT6F3FV63omnFOoBLly7F2LFjMWrUKADAmjVrcODAAWzYsAHTp08vUn758uXo0aMHpk2bBgCYN28eQkJCsGLFCqxZswYWFhYICQmROmfFihXw8vJCbGwsatWqhdu3b+Pw4cO4ePEiWrVqBQD4+eef0atXLyxevBgODg5KftRERFQZpWblYsO/D9DQ1gw9m9nLfV6+uPJe7xSr4bH9ey+p9EL/ycrNxydbI5CVm4/VQz1hYVxNct/Ze0mYtusaujexRVJ6jqRH1MPZEh0b2Cg8bkXR+AQwJycH4eHhmDFjhuSYjo4OfH19ERYWJvOcsLAwBAYGSh3z8/PD3r17i20nJSUFIpEIlpaWkjosLS0lyR8A+Pr6QkdHB+fPn0efPn2K1JGdnY3s7GzJ7dTUVAAFazhpyjpOpNkKXyt8vRCphjrec9/tv4U/Lj0GABz73Ad1rE3kOu9YZCI+/u0SAjwcMdjLWeq+isYvVvPnjvvcf2D5RlIFADfjUkuPSZD92GX9TUurq6Tyq0/cw7HIgquO8w/ewqK+zSX3DflvzcCNZ2Kkzr8Uk4zEtCwcvZWAKb6uaGRnXvJjUYCy/A01PgFMSkpCfn4+bG1tpY7b2toiMjJS5jnx8fEyy8fHx8ssn5WVha+++gqDBw+Gubm5pI6aNWtKldPT00ONGjWKrWfBggWYO3dukeMpKSn8h05yEQRBMs5UxHUMiJROHe+5wuQPAE7eeoIa7rYllH5tacgdAMCqE9F4lZUldV/WW7fLKuvVK6SkqG+WbVp2HtKypSe2ZObklxpTXn6ezDJp6elIMZSeOFNaXW/fv+3sPQRfisOYdk4Ij3m9TM6+K3Ew0RPQrk51tK5tUWx9yWmZ+OnYPQAFvY3/ftamxPYVobDjSR4anwAqW25uLgYMGABBELB69eoK1TVjxgypnsfU1FQ4OzvDwsJCklgSlaTwi4KFhQUTQCIVUPd7zsjICBYWFsjMyUPk0zR4OFtCR6f0OO4nSy+9YmBgWKE4DP+LQ9OUFpOurp7MMmamprCwMCtTXW/fP/PvgoR76u5IdG74+lJuVp4YG889wcZzTxA5z6/Y+rLFr6dZpGfnq+T5LctrWOMTQGtra+jq6iIhIUHqeEJCAuzs7GSeY2dnJ1f5wuTv4cOHOHbsmFSSZmdnV2SSSV5eHpKTk4tt18DAAAYGBkWOi0Qi/jMnuRW+XviaIVINdb7nCpsc/Mt5XHucgomd62GaXyO5zyvudpnjgGZedZAnJlllZP09S6urPI8/LStf7rKqeH7L0obGzwLW19eHp6cnQkNDJcfEYjFCQ0PRtm1bmee0bdtWqjwAhISESJUvTP7u3r2Lo0ePwsrKqkgdL1++RHh4uOTYsWPHIBaL4e2tmI2siYiIMnPyce2/RY5XHo+W65y3l31R5Cij3HwxktKztWLoUrGzgDUvl9U4Gt8DCACBgYEYMWIEWrVqBS8vLyxbtgwZGRmSWcHDhw+Ho6MjFixYAACYMmUKfHx8sGTJEvj7+2P79u24dOkS1q1bB6Ag+evXrx8iIiKwf/9+5OfnS8b11ahRA/r6+nBzc0OPHj0wduxYrFmzBrm5uZg0aRIGDRrEGcBERKQw5UmzHr94JXVbXMFkTUDB5fDhGy7g9N2C2bEftHREf09n6OvpwLN29QrVr+00PxUuO61IAAcOHIhnz55hzpw5iI+Ph4eHBw4fPiyZ6BEbGwsdndedme3atUNwcDBmzZqFmTNnwtXVFXv37kXTpk0BAE+ePMG+ffsAAB4eHlJtHT9+HJ06dQIAbN26FZMmTULXrl2ho6ODvn374qefflL+AyYiIiqDVSfk6zksybn7yZLkDwD2RDzBnognAIDxPvUqXL8yaEEnpcbSigQQACZNmoRJkybJvO/EiRNFjvXv3x/9+/eXWd7FxUWuru0aNWogODi4THESERHJ6/dzsehVhrUAlamknTHWnKx4glnVZOWJSy+kRho/BpCIiEhb5YsFRD9LL7bT4fqTFPx4JErqWP2ZB2WWJe2i6VvkMQEkIiJSktGbL6LrkpP439HitwbbEvZQ6naeGnbF4KXUqocJIBERkZKciHoGABq/N+y3+2/hm3031R2GAil2GrAiEuR3fz6NH4/I3sBCHZgAEhEREeJTK7abSGVW3CxroQzzg288ScXK49F4lJypqLAqhAkgERGRCiSmMcHSVm/Ojq6o5xk5pRdSASaAREREKjB8/QV1h0AkwQSQiIhIBSLj09QdQpU2ZvMlhdcZnZih8DpVRWvWASQiIqKq53hUInKKWVPvyqOXuPEkBU0dLUqt5+jtBNxJSIOpgR4mbI1QSGyPXpR9PJ+mbLHHBJCIiIg00vn7zzFq48USy7z787+4/W0PqWPF7QWcmZOPuX/fxNVHLxUUofbiJWAiIiLSSCuO35Or3MNk+S/FhkU/L284lQoTQCIiIiIV0YwLwEwAiYiIqJLZE/EY2Xn56g5Do3EMIBEREWmk8s6XWHk8GgZ6ukWOK3Z/EGhOd145sAeQiIhIRQ7feKruELRKVm75e/GWhtwpU/k7CVVrmR4mgEREpFLHIhPg/9NpBJ+PVXcoFZaXL3t5kuKM/10xy49UFXli1XWxdf/fKZW0oyGrwDABJCIi1fpo0yXcjEvFzD+vF7lv+dG7GLb+PO4/S1dDZPITiwUMXncOXt+H4nLsC3WHU2k9lnOdPXmTquKWhyk3hV9TVh0mgEREpBFuP03F/47ewem7SVK7NjxPz8a5+88hVmFvUGmO3k5A2P3nSM7IweBfzqk7nEorKV2x++b+evoBNOhlpFacBEJERBrh4fPXa7ndT8rA0VsJ2BwWg9N3kwAAXRvVxOoPPaGvp/6+i+SM14lJVq7sy8CasuMDvbbvapxC60tMzSrHWZrxumACSEREGmnMFum9W0MjEzHnrxtY2Le5miKST1R8GvZEPIaNmYG6QyEl23bhUZnPyS7mC4OqMQEkIiKtsf3iI41OAB8kZcBvmWomE5D6lWdM4dYLsWhX31rxwZSR+vvRiYio0ol+lo6/r8ZVucV4B6/jeMCqpDxX+VMycxUfSDmwB5CIiBTqZlwK/H/6FwDgaGmEM9O7yHVeZRgyF1+uMWFUUep67eSXY0aJoCFjANkDSERECjVt5zXJ709evlJjJMqj8OVEqEJuP01VS7vanPAzASQiIoXKfWtx5M//uIph688jMU32P8vMnDwcuv4UyZmKXfKDqo7Pd15Vdwhah5eAiYhIqXZHPAYATN99HRtGti5y//Td1xW+PIc6PE/PVncIRHJjDyAREZXoZlwK3l95BosOR8pVvrgRTsciE4uWFYRKkfzl5YvR66fT6g6DtEBaVp66QwDAHkAiIipGWlYuUrPyMHjdOaRm5eHqo5fo29IR9WualXheSQsgH4+STgK3VoL9gAHg6uOXSEhlDyCVLj5FM8YNMgEkIqIiMrLz0PGH43j5KldqhmVCanapCWD0s4xi7xu18aLU7Vl7b1QoTnURvbUJbGWYwUyqkZimGV8UeAmYiAjAy8wcbt31ht/OPcSLzFyZiU1aVi72XY3DMw35R6YJHr+onLOdqfJiAkhEVd7OS4/Qcl4IRm26WHrhKiIrt/gFnKftvIZPt13G4F/Uv+hxRnYeMrIrNqYqMj4dq07cK3aWskxvLQMzdceVCsVApGpMAImoSsnJE+NiTLLUUiXTdl2DWABORD3D0xT25JQkTyzg8M14AMC9xHS1xJCZU5DwPXn5Ct7fh8L7+9ByrzcoFgsYtOkqfjxyB+N/C1dkmEQajQkgEVUpn+24gv5rwvD5H7LXDcvLF5CckVOuFf612bd/38KANWGIflZyUjfut0sqiqh4jeccgSAICPrrBtKz85CenYc55RxLmPPGF4GI2JcKipBI8zEBJKIq5cD1pwBQ7NIj/9xKQOv5R9H7538hriJJ4KWYZGw48wAXYpIxdnPJCV5WrvQizzfjUvDwefGTPpTlflIGnme8Xji6tMRV0bgRCGk7JoBERG+Yt/8W8sUCbj1NxZnoJHWHU6KcPDEmBkdg9KaLSHlV/g3mP1x/XvL7/aSyJXP+P/2LLktOIk7FW769ysnHo+RMye2Y55kllCaitzEBJKIq680EQpY3e7sEQUBMUoZGzRRe/+8DHLj2FKGRiRi4Nqzc9bzdqwcA1x6nyH1+vljAkn/ulLv98nj353+RlK76reNe5eRj24VYXi4mrccEkIiqjLeTt57LT8s91m/u37fQafEJfLHzmjJCK5fwhy8kv0fGpyms3vP3n8vctaMkQrH7f2ifQ9efov+aszh6K6HIff87egcz9lzHtguVYwFrqrqYABJRlRH81j/t9Oy8Ei+dvpkwbjobA+D1vrbqlJSejdN3nyE9u/yXfUsycJ36l3dRFVlj+SZsjcDFmBcYs0V6PGROnhjrTt1XTWBESqYVCeDKlSvh4uICQ0NDeHt748KFCyWW37lzJxo1agRDQ0M0a9YMBw8elLp/z5496N69O6ysrCASiXDlypUidURHR6NPnz6wsbGBubk5BgwYgISEot8GiUh7fP2ndu468aZv/76FVt8dxbD1F3DufrK6w5F4np6DG0/kv2ysbQRBQJ9VZ9QdBpHCaHwCuGPHDgQGBiIoKAgRERFwd3eHn58fEhNlX544e/YsBg8ejNGjR+Py5csICAhAQEAAbtx4/cGfkZGB9u3bY9GiRTLryMjIQPfu3SESiXDs2DGcOXMGOTk56N27N8TiomNliIhUZcOZB8XeF1vGiRD5YgGTgiMqGhIA4OSdZ3j353/xdzGzq7Xd3cR03IxLVXcYRAqj8Qng0qVLMXbsWIwaNQqNGzfGmjVrYGxsjA0bNsgsv3z5cvTo0QPTpk2Dm5sb5s2bh5YtW2LFihWSMsOGDcOcOXPg6+srs44zZ84gJiYGmzZtQrNmzdCsWTNs3rwZly5dwrFjx5TyOIlI8xReAM7OK35XDE3S8cfjZSq/J+Ix9l97qtAYJm+7rND6NMWsStB7TPQmjU4Ac3JyEB4eLpWo6ejowNfXF2Fhsme8hYWFFUns/Pz8ii0vS3Z2NkQiEQwMDCTHDA0NoaOjg3///beMj4KINFlJs3rH/RaO/mvOYtj6koedaJKybIt2V007eWib6GfpuBCjOZfbiRRBT90BlCQpKQn5+fmwtbWVOm5ra4vIyEiZ58THx8ssHx8fL3e7bdq0gYmJCb766it8//33EAQB06dPR35+Pp4+Lf7bcnZ2NrKzX2+OnppacLlAEASNWjqCNFfha4WvF9URl/JcX4x5UeSYJv99fgq9i+k9G8lVVpMfR3mU5/G8fY6s28/Ts0GkSMp675WlXo1OANXFxsYGO3fuxIQJE/DTTz9BR0cHgwcPRsuWLaGjU3yn6YIFCzB37twix1NSUirdBy0phyAISE8v6JURibjXgCrM3Vv2ZV1SUjR3ssPV2GS543vzC2tlsPLobZgb6uG9ZjXlPicrV/ry/tvPXUpKCh7Es/ePFEtZnyGFHU/y0OgE0NraGrq6ukVm3yYkJMDOzk7mOXZ2dmUqX5zu3bsjOjoaSUlJ0NPTg6WlJezs7FC3bt1iz5kxYwYCAwMlt1NTU+Hs7AwLCwuYm5uXqX2qmgq/KFhYWDABVJH9N5+V+RwLCwslRKIYzzJy5Y7vzWEulcHi0IIJMvXsq6NdPWu5zjHIkb5k/vZzZ2FhgZWnK+e4RlIfZX2GlOX/hkYngPr6+vD09ERoaCgCAgIAAGKxGKGhoZg0aZLMc9q2bYvQ0FBMnTpVciwkJARt27YtVwzW1gUfIseOHUNiYiLee++9YssaGBjI/EAViUT8Z05yK3y98DWjuTT5bxP9LEPu+DT5cVTEnog4vFPfRq6ybz8Hsm7n5HH1B1IsZb33Kk0CCACBgYEYMWIEWrVqBS8vLyxbtgwZGRkYNWoUAGD48OFwdHTEggULAABTpkyBj48PlixZAn9/f2zfvh2XLl3CunXrJHUmJycjNjYWcXEFyxVERUUBKOg9LOwp3LhxI9zc3GBjY4OwsDBMmTIFn332GRo2bKjKh09ERGr08HkG4lOz1B0GkcJpfAI4cOBAPHv2DHPmzEF8fDw8PDxw+PBhyUSP2NhYqXF57dq1Q3BwMGbNmoWZM2fC1dUVe/fuRdOmTSVl9u3bJ0kgAWDQoEEAgKCgIHzzzTcACpLCGTNmIDk5GS4uLvj666/x2WefqeARExGRypTSYeLz4wmVhEGkahqfAALApEmTir3ke+LEiSLH+vfvj/79+xdb38iRIzFy5MgS21y4cCEWLlxYljCJSENl5uTh4y3h6g6DtMC6U9HqDoFIJTR6HUAiIkX4KfQe/r2XpO4wVEbeVQfKunOIthBQ/lUXvj8oe4kxosqGCSARVXrXn7xUaH1f/3m9yPIhmuTCA/mWLTl8U/71UbXJo+RM5IsFfLPvJgL/uIK0rFx1h0SkcZgAEhGV0dbzsVh1QnMvFc4/eFvdIajVxZgX+OPSI2w6G4M9EU+w+EiUukMi0jhMAImIyuHwDcXuoatI1x5r7kLVqnI2+rnk981hD7Hy+D0kcjYvkQQTQCKq9JSxEQ8399EuPx6JwrjfORGIqBATQCKicribmK7uEKiMLse+xNd/XseZKjQhiKg4TACJiMrJZfoBjNhwAU9TXqk7FHpLccv7bT0fi6G/npfc5i4fVFUxASQiqoCTd56h3+owdYdB5ZAvFtBj2Wl1h0GkFkwAiYgq6MlLzesBvBRT8lIw+WIOYjx6OwFxKZwYQlUTE0AiqvSq4oSNfmtk90qmZeXCZfoB1Jt5UMURqVZiWsmJ3b3ENIz7jZNCqOpiAkhEpCXk3eGjJF/uuqaASDTfufsl94C+t+KMiiIh0kxMAImItERE7MsK13HoRuXc/aOsMnM0dycXIlVgAkhElV5F9obVJKkV3NIsN58zXomoABNAIqr0SrscWBUIggDv70PVHQYRaQgmgESkFg+fZ+DHI5G48aRybFsWGZ+q7hBKdP1JCpIzctQdBhFpCCaARKQWA9aGYeXxaLz787/qDkUheiw7XeFLtMqSmy9G4B9X1R0GEWkQJoBEpBYJqdnqDkHhFh+JUncIMn2z7ybuces6InoDE0AiIgW5/yxD3SFIycjOA1Cw/RkR0ZuYABIRVVJL/rmj7hCISEMxASQiqqQ2nHmg7hCISEMxASQiIiKqYpgAEpHa/XHxEfqvOYuw6OfqDkWzVY71rIlIAzABJCKVS0zNkrr95e5ruBjzAoN/Oafw3Spuxalufb5XudxejIi0AxNAIioiMTULJ+88Q56Stg47dTep2Pu85h/F1UcvFdZW7xWqW2cw/OEL7Il4rLwGRGU/5afQu4qPg4i0HhNAIpKSmy+G/8//YsSGC1h76r7K23+RmYsP159XWH35YtVeNw384yruJqSptM2SLA3hTGAiKqrMCeCrV6/w5MmTIsdv3rypkICISL2uP0nBs7SCRZp/VNPCxmlZeWppV1FGbryo7hCIiEpUpgRw165dcHV1hb+/P5o3b47z519/Sx82bJjCgyOiykcQBPx5WYmXSTXAk5ev1B0CEVGJypQAfvfddwgPD8eVK1ewceNGjB49GsHBwQAKPtSJSLucvPMMa09GIz07D4mpWfj7ahxe5Sh3IsOpu0k4c081s31z8pQzhpGISNvplaVwbm4ubG1tAQCenp44deoU+vTpg3v37kEkKsfoZCJSm7iXrzBiwwUAQGxyJk5EPcOTl6/gYmWstDYzc/Iw568bSqv/bev/5ULIRESylKkHsGbNmrh27Zrkdo0aNRASEoLbt29LHScizXf+weteuK3nYyWXLWOeZyqlvYzsPDSecwQP5az/+uOUCl9ZWHQ4skLnV8SBa0+xK/wxHiRp1v7ARERAGXsAf/vtN+jpSZ+ir6+Pbdu2YdKkSQoNjIgqlyG/lm1mb+HyLQ8W9NLKKwwTgyMAANV0Rbj1bQ+8yMxBUloOGjuYqzkyIqIy9gA6OTnBzs6uyPHc3FzUqlULUVFRSE5OVlhwRKQ8IjkXlbv+OEUh7ZV3bT9VjRdUltx8AdHP0tH5xxPo9dNp/HMzXt0hERGVfx3AtLQ0rF69Gj4+PjA3N4eLiwvc3NxgY2OD2rVrY+zYsbh4kUshEGkqeTvVvtyt3uEdD5O1/xJq8PlYZPw3uebj38LLXxHn2hGRgpQrAVy6dClcXFywceNG+Pr6Yu/evbhy5Qru3LmDsLAwBAUFIS8vD927d0ePHj1w9y5XoifSVoq4+FrVZ+OeLmHnEyIidSjTGMBCFy9exKlTp9CkSROZ93t5eeGjjz7CmjVrsHHjRpw+fRqurq4VCpSIFEuV4+qWh1bt3SgUNhFE+4ZCEpGGKlcCuG3bNrnKGRgYYPz48eVpgoiU7Px9+cbWKSJPXHk8utznyjtWkYiI5Me9gImqqK3nY1XSzsk7zyp0/uMXylmWhoioKlNKAjh48GCkpqYCAPbt24edO3cqoxkiUoGbcakIf1j+2f3hMRVbGWDVifL3HlY6nARCRAqilATw5s2bMDc3x61btzBz5kycOHECkydPLnd9K1euhIuLCwwNDeHt7Y0LFy6UWH7nzp1o1KgRDA0N0axZMxw8eFDq/j179qB79+6wsrKCSCTClStXitQRHx+PYcOGwc7ODiYmJmjZsiV2795d7sdApM36rg4r96LMzFmIiDSPUhLAatWqQRAEbNy4ETNmzMDKlStx9uzZctW1Y8cOBAYGIigoCBEREXB3d4efnx8SExNllj979iwGDx6M0aNH4/LlywgICEBAQABu3Hi9/VRGRgbat2+PRYsWFdvu8OHDERUVhX379uH69ev44IMPMGDAAFy+fLlcj4OIVCP84Qt8sOoMlh/V3NUHsnLLud8yh0MSkYIoJQGcMGECWrZsiV27diEgIABAQdJVHkuXLsXYsWMxatQoNG7cGGvWrIGxsTE2bNggs/zy5cvRo0cPTJs2DW5ubpg3bx5atmyJFStWSMoMGzYMc+bMga+vb7Htnj17FpMnT4aXlxfq1q2LWbNmwdLSEuHhFVjDi6gKquBubmXWd/VZRMS+xP+Oau7M48VHotQdAhFVceWaBVyaMWPGoF+/ftDT04OJiQnu3buHNm3alLmenJwchIeHY8aMGZJjOjo68PX1RVhYmMxzwsLCEBgYKHXMz88Pe/fuLVPb7dq1w44dO+Dv7w9LS0v88ccfyMrKQqdOnYo9Jzs7G9nZ2ZLbheMgBUGo8J6mVDUUvlY08fUiFgvQKcdXRrECHosmPh8V8eu/D/C1v1uZzknPzkNmdp6SIiIiVVLWZ1pZ6q1QAvjw4UNERUWhefPmRbaIs7S0RFxcHExNTVG/fn1s2rSpzPUnJSUhPz8ftra2UsdtbW0RGSl7k/f4+HiZ5ePjy7b90h9//IGBAwfCysoKenp6MDY2xp9//on69esXe86CBQswd+7cIsdTUiq+qT1VDYIgID09HYBq1+mTR0pqCnTKEdObX4rK3XaKYraj0yRleUyXH6di1O/XlRgNEamSsj7TCjue5FHuBHDbtm0YPnw48vPzYWhoiLVr12LYsGGIjY1FcHAw9uzZg4iICOTlaec31tmzZ+Ply5c4evQorK2tsXfvXgwYMACnT59Gs2bNZJ4zY8YMqd7H1NRUODs7w8LCAubm3ACeSlf4RcHCwkLjEkALcwvo6JQ9JgMDg4q3bWFR4To0zdwjMTAz1MOcdxtDt5TnddTCMyqKiohUQVmfaWX5v1HuBHDevHmYPHkyRo8ejZkzZ2LChAm4c+cOFi5ciHr16qFr165Sl27Lw9raGrq6ukhISJA6npCQUKTHsZCdnV2ZyssSHR2NFStW4MaNG5LdTtzd3XH69GmsXLkSa9askXmegYGBzH92IpFI4/6Zk+YqfL1o3GumnDEpou/7ZlwqmjpWriTwz8tPAABu9uYY7FVLzdEQkSop6/O9LPWWexJIdHQ0pkyZgiZNmmDlypXIzMzEmTNncO3aNURGRmLlypXo06dPeasHAOjr68PT0xOhoaGSY2KxGKGhoWjbtq3Mc9q2bStVHgBCQkKKLS9LZmbBwrM6bw140tXVhVhctfc0JVKHd3/+t9IOozgn544sRESKVO4ewNzcXBgZGQEAnJycYGhoiMWLF8PNrWwDm0sTGBiIESNGoFWrVvDy8sKyZcuQkZGBUaNGAShYrsXR0RELFiwAAEyZMgU+Pj5YsmQJ/P39sX37dly6dAnr1q2T1JmcnIzY2FjExcUBAKKiCmbk2dnZwc7ODo0aNUL9+vUxbtw4LF68GFZWVti7dy9CQkKwf/9+hT4+Im1RkICVowdQQXnbqhPRmNi5+DG42qqS5rVEpOEqtAxMcHCwZDKGrq4uqlevrpCg3jRw4EAsXrwYc+bMgYeHB65cuYLDhw9LJnrExsbi6dOnkvLt2rVDcHAw1q1bB3d3d+zatQt79+5F06ZNJWX27duHFi1awN/fHwAwaNAgtGjRQnJpt1q1ajh48CBsbGzQu3dvNG/eHFu2bMHmzZvRq1cvhT9GIirdj3IsnXLw+tNSy2iafVfj1B0CEVVBIqGc11V8fHxw5coVpKeno3r16khJScHEiRPRrl07NG3aFA0aNICenlJWmdEaqampsLCwQEpKCieBkFwEQUBKSopKJoG4TD9QpvL35veEnm7ZvzO2nBeC5IycMp8nS8xC/2Lvu/EkBe/+/K9C2lG1kh4XUPa/FRFpttLe8+VVlryj3BnayZMnAQB3795FeHg4IiIiEBERgS1btuDly5fQ19dHgwYNcO3atfI2QUQapLxXKhWV/JXmyM2yLfVERFSVVbiLztXVFa6urhg0aJDk2IMHD3Dp0iVum0ZERESkgco9BnDOnDnFbotWp04d9O/fH99//325AyMizcLJCsoT/vCFQuopxzKNRFRFlTsBfPz4MXr27AknJydMmDABhw4dQk6Oai71EJF2+OvKE3WHoBX6rj4LsbjiGfbBKR0UEA0RVQXlTgA3bNiA+Ph4bNu2DWZmZpg6dSqsra3Rt29fbNmyBcnJyYqMk4jU7LdzD8u8Zt2U7VeUE0wl9OB5BvLyi64zKutYcRrZmWPV0JZobM9JZ0RUsgotA6Ojo4MOHTrghx9+QFRUFM6fPw9vb2+sXbsW9vb26NixIxYvXownT9gLQKTt5u2/hUHrzuFpyit1hyKTtl/97LrkJOp/fQg5edIJ3+c7r8p1/rvN7QEAvZrZ4+CUDrj4ta/CYyTV2j+5Pcb51FV3GKQEmrCwfYUSwLe5ubnhyy+/xJkzZ/D48WOMHDkSp0+fxrZt2xTZDBGp0b93k+Qq9/B5hpIjeW3v5Sf46dg9lbWnTHsvS39h/utK6esELh/kgZ8GtZA6ZmNW8T2YSb2aOlpgfMd66g6DKimlLNQXERGBJUuW4MWLF2jWrBn69u2rjGaISIOdv6+aYSDxKVmYuuOKStpShfTsvDKf08TBAjqcAVIpadqW4KQYgqD+v61CewALDRw4EP7+/pg/fz4aNGiADz74AP/8848ymiIiFTsg524bQrlXDiyerEkl95PSFd6OOhX+U1DEpJD/DXSvcB2kXiYGVXtDhcpK/ReAlZQAVq9eHUOGDEGLFi0wevRoHD16FNOmTVNGU0SkYieinslVThlDXKZsvwKX6Qdw5t7ry9AirR/9J00EYGJwBFrPP4qwaPkm3RTXkxDg4ai4wEgtqunq4I9xbTGynYu6Q6FKRikJYL169bBkyRLJIEdLS0tlNENEanLhQemXd5X5DXfor+clv6v7MoqiLQm5gwPXnuJ5Rg4G/3KuQnUpeztBUp7b3/aQ/O5Vpwa+ea+JGqMhRat0k0AKZWdnY/Xq1ahVqxZ69OiBpk2bwtfXl7OBiSqJAWvDSi2jqs+3qPg01TSkImlZZR8DaMbLhJWOkb5ukWOc2FN5qD/9U9IkkD179gAAMjIycP36dVy7dg3Xrl3D4MGDERcXh3v3KsdsPSIqnjLGAMoStO+mStrRZDXNDdUdApVBxwY2OHXn9VAKo2q6eJWbX+p5G0e2Rv81YXKVJSqNUr82mpiYoE2bNmjTpo0ymyEiDaQBVzgIgLuzJa4+eqnuMOgNG0a0Qv2vDwEAPu/WAOM71cPz9By0WRAKoGBZH1maOlrgSlA3NJx1WFWhkpJowudjhS4BP3z4EP/88w/i4+Nl3h8XV/r6VURUOSWmZas7BAJQq4axukMok/E+lXvdu+WDPKCnq4M73/XEyWmdMLmrK6rp6sDOwhBHpnbElo+80Lu5Q7HnG+jpYv2IViqMmJRBVVdISlLuBHDbtm2oX78+evTogbp16+K3334DAMTGxmLhwoXw9vZGrVq1FBYoEWmXW3Gp6g6BALjZm6k7BLl92qU+Pu/eQN1hKNX7/83M1tfTQW0rE6n7GtqZoWMDm1LXdOzSqCam+TVUWoxUNZQ7AZw3bx4mT56M69evo1u3bpgwYQJmz56NevXqYdOmTWjVqhV27typyFiJSIscvZ2g1PrzFbBOXlUwrE1tdYcgt8+6NUA1Xdn/llq7VEcjO+1JZpVJJBJhYuf66g6DKkCrLwFHR0djypQpaNKkCVauXInMzEycOXMG165dQ2RkJFauXIk+ffooMlYiApCXL8aXu64i8I8reJWTj9/PPcQ7C49h56VHctdxMy6lwnEcviHfgtDK8v3B22ptX1uYGVZTdwhyK1y2Zs8n7Yrct3JIS3jWrq7qkCqsg6u15PfBXs4KrVuvhJ7C4DHeCm2LFEtHA5ZoKvckkNzcXBgZGQEAnJycYGhoiMWLF8PNzU1hwRFRUd8duI0/Lj0GAEQ8fIGY55kAgGm7rqF/K/n+wVxRwKSA8b9HYOkAd3zQ0qnCdZXH+n8fYPa7jdXSNilXy1rSiZ6XSw3UNDdEy1rVsfV8rJqiKruaZgZYPqgFzAz18PjFK9SxNin9pDJ4sxNpqHctNHeywL3EdAzyqoV6NqYI+awjuv3vlKTM8kEe2Hw2BvcS0zGhU30sOhyp0HhIfvp6SlmFr0wqNAs4ODgYPXr0QKNGjaCrq4vq1bXv2xmRttl0Nkbye2Hypy6Bf1xVWwKoaYLHegMCMOSNRapJfp/5Fj/275f/Jj180NIRn++8qqqQKuzUl51hWK1gPT9FJ39v09URYWBr6XH3rrZm6OfphL2Xn+C7gKZ438NRMgYRAKobV8P0PdeVGhdprnKnoB06dEBQUBCaNGkCa2trZGVlYfny5fjjjz9w69Yt5OWVfTFTIqqYH49EIi9frO4wqiTvOlZoV9+69IJq0FcLkvRPuxY/ps3CqOAytkgkwokvOqkoote2jW2Dxf3dMd6nHhrYmpZa3sxAD5HzekiSP2VZ+EEzye/D28oe67m4vzuuf+OHQV5FJ2UO8qqFNR96YqqvKy5+7au0OEkzlbsH8OTJkwCAu3fvIjw8HBEREYiIiMCWLVvw8uVL6Ovro0GDBrh27ZrCgiWikq08Hg07c0MMa+ui7lBUIlbNPaBv0v1vPNa37zfBnL9eL0796/BWGLPlkrrCAgDM6d0YuyMeqzWGkvTzdCqybd3CD5rhf0fv4JNO0omhi7UJYhb6F/w+/YDSYopZ6A+xWCgyI/erHg1RZ8ZBmed0bmiDUe/UgZu9udKTPwD4oKUTTA30YGVqgPo1i58gI2tXkUI9mtqhR1M7AIBPAxuc/G+BaqfqRnj84pViAyaNUuGFoF1dXeHq6opBgwZJjj148ACXLl3C5cuXK1o9EZXR39eeVpkEsOOPx9UdQhHD27pIJYAetSzVF8x/CnvQNNWP/ZoXOTbIq5bMXitVcK1Z0MsnazkWkUiEBwt64VlaNmqaG+JWXCp6/XQaberWwLKBLWBhrLrnWldHhJ7N7BVWX2uX6pIEsH19a2y/KP/EMtI+5UoAY2NjS1zjr06dOqhTpw769+8PAHjy5AkcHR2LLU9EqqXIJQjy8sXQK2bpjqri3Iyuxd5nZsh9ekvyfZ9mRXr/5OXrZquU5YZkzUJ+k0gkkmy/19jBXNIjqe3GdqyL2/FpyM4VY5pfQyaAlVy5PrVbt26NcePG4eLFi8WWSUlJwS+//IKmTZti9+7d5Q6QiBRPUGAGOOSXqj3pwcxQD3YW0nvx1v1vwL+bvTkM9HRhZaKvjtC0QkWWRilpGZTyOhroo1VL5yiSgZ4uVg5piV9HtIL5W73G73tI704yiesQar1yfTW9desW5s+fj27dusHQ0BCenp5wcHCAoaEhXrx4gVu3buHmzZto2bIlfvjhB/Tq1UvRcRNRMTJzSp+Atf+a4tbwuxCTjKT0bFibGiisTm2y5kPPIsdCAn3wICkd9WwKLiVemuVb7LgxVWlsb45bTzVvd5by9v4BgEE1xfc8169Z+iSPqkBPR4TmTha49jgFo9vXweQu9WFnbogD15/C180WEzvXx4rj99QdJlVAud49VlZWWLp0KZ4+fYoVK1bA1dUVSUlJuHv3LgBg6NChCA8PR1hYGJM/orcosvdNlhtPUiEuZZeMa48rvhD0m8SasKy9Grg7WeAdGTN/dXVEqF/TTJLciEQi3J3fU9XhSdn8kZda23/TP591xOj2dbB/cvsK1fNVj0ZSt2f5u+HBgl4w52X3ChOJRPh9jDe2jvHG9J6NYGmsjxm93PDvV13wzXtNYKSvi4mdy7Zv83cBTZUS68lpnTC/j3Lqrswq9C4xMjJCv3790K9fP0XFQ1SpLTociZ2XHuO7gKbo0dQO9xLTse1CLHq7O8DD2VJh7dyIS0Fzp+LrU/hG5G9Vdy8xTbH1q5ClcTW8zMwttdzWMd4yk7/iVNPVwdWg7riTkAbPWtWx4cwDfHdAdbuZ2Jippoe2rrUJ7idlFHu/UTVdNLA1U8gi3g6WRtg/uT1GbbqI5o4WGNHOBSKRCFeDumPr+VhsCYvBnYT0CrdTVZkbVivxNT6psyuMqumiprkh3qlvjRrG+nCbc7jY8oO9aiFo303kiwVM7lIfPx+T7kF0rmEEPR0dPCjh9fM2Fytj1LYyQW0rE9xNSJdaJ5VKVrVHbhOpUHZePlafiEZSejbG/x4OAOiz8gzW//sAASvPyFXHxZhkucptu/CoxL1ylb2N7sSt2rsCQIOaZhjZzqXUcmVJ/gpZGFVDa5ca0NER4aN36pQjOs13NNAHs/zd0L2xLY4Gdixyv4mBYpdHaepogYtf+2L9yNaSfYRFIhE+bFMb/3zmU6a6lDGmsDIz0tfFpC6uGNDKGY6WRjDS18X5mV3xZY+GOPhpB6lxg592qQ9dHRGiv++FmIX++Lx7Q7jZm0vV911AM2z/uI3k9tphRYdXvG3JAA/J735N7Cr+oKoQ9pMTqYhYxvrMadllWzD9wTP5vhlvuxCLFs6WGNBa9gD7nDzFLhb9dj4Z/UyLe11EQFDvxhjY2hmPX7zC2P/W8Ns2tg2epWdj1p/X8YVfwwo3I2uJEW33RfcG0NERYUyHuhjToa66wwFQcLm5+xvboZWkag5kUCxbc0PJ2o1LB3igb0sn5IsFtHct+oVp5ZAWmPB7BKIS0rBySEt0dLWGSCTC0UAfPE/PhledGsW2M82vIepam0jtD922npXkd//m9jDQ1cGey08U+OgqFyaARFqkLJduv9x9rdgEkIrXopYlRCIR3OzN0cjODL8ObwVdXZHkn8t77g6l1KC53OzNcVtJE0H2TnxH5jCGee83wew31kX8vk+zImWUqYGtmUoWjqaidHVE6NjAptj769qY4shnRXuJ69c0lUzG+bqXGzaeeYCvejaCcw1jzNt/C35N7DDeR/b4wwcLeuHh80zUtjKGSCRCQAtHDN9wQTEPqJJhAkikIhWY7Kjxrjx6WWkuv7y5vIVIJIJvY1ultbV34jtyX/5XhK6NaiotAaxezALIw9q6YLBXLVx/koK0rDx0kNETpCqFu1t0b2yLnwa3QNzLV4iIfYEvdhbsWDVWQ3ot6bWxHetibMfXf5c/P3mnxPIikQgub+y73MHVGr8MbwVBENDKpQZazgtRWqzahgkgEVXYuN/C8WBBrwot6aEpVLkGnCIn/shjQqd6Slu6o6RlgPR0ddCiVvVi71eV3RPaISz6Obq61YRhNV3UtTFFHWsTJKekIzVXB5+UcVYraT6RSIRu/32Je3sFhh/7Nce0XVV3u1pOAiHSImlZZRszqEqLDkdJfq8EeWClZGKgvO/8yqxbUWzNDRHQwrFIkj+gpT0+794Axvqa/xio/EQiEf6e1B5edWog9HMf9G9VtYfIMAEkUpG3k6JXOfllrmPhoUgFRaN4a05GS37PzedweiLSPM2cLPDHuLaSRdrPTu8iuS9sRhfELPRH6Odlmz2urZgAEqmICNIZYEnrZRUnT9nrt1RQl8UncDwqUd1haJWxHVS7HEzbulalFyqjlUNaKrxOIlVwsDTC7W974Pa3PWBvYQQAqGdjihk9G5VypvZjAkikIZS9Q4gq3E/KwKiNxe8Rrul6NVP9RJZPu7qqtD2/Joqf1OLf3F7hdRKpipG+Loz0pdenHOdTD190b6CmiFRDKxLAlStXwsXFBYaGhvD29saFCyVP6d65cycaNWoEQ0NDNGvWDAcPSu/BuWfPHnTv3h1WVlYQiUS4cuWK1P0xMTEQiUQyf3bu3Knoh0dVRGnj4rr/7xSeprxSTTAkk4kaxoCpctIJAHzYprZC67M21VdofUSaYlIXV1z/pru6w1AajU8Ad+zYgcDAQAQFBSEiIgLu7u7w8/NDYqLsy0xnz57F4MGDMXr0aFy+fBkBAQEICAjAjRs3JGUyMjLQvn17LFq0SGYdzs7OePr0qdTP3LlzYWpqip491bufJ1VedxPTMX33dXWHUaVVhckrero6mNlLMZe3fuzXHH9Nqth+vkSazMywGnZ83AZj2ksP1bgypxuWDfRQT1AKovEJ4NKlSzF27FiMGjUKjRs3xpo1a2BsbIwNGzbILL98+XL06NED06ZNg5ubG+bNm4eWLVtixYoVkjLDhg3DnDlz4OvrK7MOXV1d2NnZSf38+eefGDBgAExNTZXyOKnykye3uP4kRelxkOb5sV9zlbY3sHWtCtexdIA7+v+3BRhRZeZd1wqz3m2MAA8H6Ovp4H8D3WFprI929aXH0w7xrvj7SpU0es57Tk4OwsPDMWPGDMkxHR0d+Pr6IiwsTOY5YWFhCAwMlDrm5+eHvXv3ljuO8PBwXLlyBStXrix3HUTyrJGnimGA2Xlln32sbSLn9YCBng5ORD3D3itP8NeVOHWHVKIPWjqpdD0yC6OyXXae5tcQE3zqQUdHhOhn6Xj6Mgvt6il+MgmRJls2qAUW5ubDsFrBeMG3J/bNe78pgs/HqiO0ctHoBDApKQn5+fmwtZUetGxra4vISNnLYcTHx8ssHx8fX+441q9fDzc3N7Rr167EctnZ2cjOzpbcTk0tWHFfEIRKMcCfKkbe14AiXy+y6tGmD6jy2DSyNQz0Ci5udGpog04NbTDBpx56LD8t1/nqeK9WdFvg8sTsXN0Ij16UPub0wYJeUu3UtTZB3f92WqgMn2uF77fK8FhI+Qz0dCSvlbdfMzoiYPf4tvh851XEPM8ssR5lvd7KUq9GJ4Ca4NWrVwgODsbs2bNLLbtgwQLMnTu3yPGUlBR+uBAS0rJLLSMWBKSnpwOQr8ewNCkp0peUBUHA3L9vVbheTeZhp1/kcdsZAePeccbaM49KPLdTPfMi52qD8sS87+MWaLHobJHjpga6SM8u6CX+vncDrXw+ykJQ8HuOqg4hTyz53c684HOnnqUO9o5tAY+FxW/xOLN7XaW9rwo7nuSh0QmgtbU1dHV1kZCQIHU8ISEBdnayl2uws7MrU/nS7Nq1C5mZmRg+fHipZWfMmCF1+Tk1NRXOzs6wsLCAubl5udqnymPzpbtylTM1NYWFhUWRf0bZuWW/dGthYSF1+8jN8veEazrn6kaYF9C0yGMuNK1X0xITwO8CmqCnh2JnyJbF+x4O5b5UXdxjLo+rc7ojN1+MhLRs1KphrLB6NVXhl3NZ7zmi0qwb5onQyESM61gXFhav9yDeNtYbg385X6S8Z+3qGNWxAfR0lTMFoyyvYY1OAPX19eHp6YnQ0FAEBAQAAMRiMUJDQzFp0iSZ57Rt2xahoaGYOnWq5FhISAjatm1brhjWr1+P9957DzY2NqWWNTAwgIFB0f0wC5eQoapOjjGAgNSyQ29adfJ+mVvMFwuSD5rUrFyM/z2izHVog2vfdId5Kcup6OmKsHtCW/RdLXv88IdtXJQQmfyWDvAodwJY3s+Xfz7riO7/OyW5PcvfDbq6OtDV1UFtK43+96BQxb3niErTvYkdujcp2sHUtp41No5qjX9uJuDa45e4GVfQM7d2mCeq6ekWKa8olSYBBIDAwECMGDECrVq1gpeXF5YtW4aMjAyMGjUKADB8+HA4OjpiwYIFAIApU6bAx8cHS5Ysgb+/P7Zv345Lly5h3bp1kjqTk5MRGxuLuLiCD9uoqII9TAtn/Ba6d+8eTp06VWQdQSJ1+ClUvh7EN+WJBRR+1iw4eFvBEWmGy7O7lZr8FfKsXQMxC/0BAOEPk4tNBtVBV0cEo2q6eFXGnt5JneuXu80GtmZ4sKAX/rj0CKmv8jC8rUu56yIiaZ0b1kTnhjXxMjMHv597CHdnS1ibFu0kUheNXwZm4MCBWLx4MebMmQMPDw9cuXIFhw8flkz0iI2NxdOnTyXl27Vrh+DgYKxbtw7u7u7YtWsX9u7di6ZNm0rK7Nu3Dy1atIC/f8E/gkGDBqFFixZYs2aNVNsbNmyAk5MTunevvAtBkoZR8FDR2OTXA5Fvxck/NkTd+ns64dv3mwAAapjo45/POqK5k/RlTnNDPTxY0AvVTcq3ELFn7RoVjlPRVn1Y9i3VJnUpfwIIFPQYDGxdC2M71oW+nsb/SyDSOpbG+pjUxRUdXEu/kqhKIoGzE5QmNTUVFhYWSElJ4RhAwrKjd7DsaMm9eKYGevj3M2+Z45Fcph8oV7vnZ3aFrbkh3lvxL6491vwB/X9Pao9mThYQiwWE3X+OejamsLMwRG6+GBN+D8ftp2kY7OWMId61UaOcyV8h/59O42ZcKvp5OmFxf3cFPYKKKevfubBHk8pHEASkpKRwDCBVCmXJOzT+EjBRZfFYjiU30rPzFD5jfMWxe5gX0BR5+drxXc/qv63FdHREeKe+teR4NV0d/DqitULb+n20Ny7GJKNjA835Zv73pPboveJfdYdBRJUc+/uJNMzJe8kKre+3cw8REfsCt55qxyXgsi5SXBHVTfTRvYmdZGFXTdDEgVcLiEj5mAASqYi8F5dik7MU3vYHq4qu96aJpnR1hYlB1b4woaMjwqEpHdQdBhFVckwAiVRE3uFF2nGhVjk+69ZA3SFoBDd79gISkXIxASTSMEIVTQE9a1dXdwga5cqcbqWW2fKRlwoiIaLKiAkgkYapqvPyF/Vtpu4QNIqlsT6OTO2IDq7WxZbRpMkrRKRdmAASqYhIzlGAy088REKq4scBarr6Nc3UHYLGaWhnhi0feSF4jDf2TnxH3eEQUSXCBJBIiTJz8vDLqfsIvZ0g9xhAAPhy1zXlBVVBV+cofmH0jzvWVXidlYVIJEK7+tbwcLZEp4ave/w8nC3VFxQRab2qPd2OSMmW/HMH6/99AABoU1f+nSdO3U1SVkgV8n2fZrAwrobQz33QdclJhdVb28pYYXVVZj8PboHd4Y+Rky/GBy2d1B0OEWkxJoBESlSY/AHAufuKXd9PlT7v1gCTutSX7JRQz8ZUofUbadA6fJrMzLAaRr5TR91hEFElwEvARBoq9HaCukMAAAR4OGByV9ci22RdmNkV/T0r3gtVq4Yxers7VLgeIiKSHxNAIg01evMlnI1W/6XgOb2byDxe09wQP/Z3x8+DW1So/qOBPqimy48iIiJV4qcukQYb8st5tbYf/X0v1DDRL7FMb3cHPFjQq1z1z/J3g74eP4aIiFSNn7xEJJODhSF0deSbuiwSiXBvfs8yt/ERx7MREakFE0AiKqJzQxsEj21TpnP0dHUQs9Afs/zd5CofMbsbdORMMImISLGYABKV0/XHKRj32yUcvP5Uqe3EvXyl1Ppl+ea9JnCxNinXuWM61C11h4pZ/m6lXlomIiLlYQJIVE69V/yLIzcT8MnWCJyNTkKTOYex9mS0wtv5cL3qxwFamxpU6PwtH3nhwtddZd634INmGN2el36JiNSJ6wASKUDhZI0FhyLxTn1rNHW0UFjd959lKKwueZkYVPyjoaaZIWzMDPAsLRsAcHOun0LqJSKiiuOnMZGC3YxLUWgCCABisaDQ+kpybobsnrvy+GNcW/x+7iF6NbNj8kdEpEH4iUykBY5HJaqknbPTu8DOwlBh9dWxNsHsdxsrrD4iIlIMjgEkUrCcPLHC6xy9+ZLC63zbH+PawsHSSOntEBGR+jEBJFKw7w7cVncIREREJWICSPQfQRDwRAFLrmQroQdQFRyrs/ePiKiqYAJI9J9JwZfxzsJjWHb0jrpDUQtHXv4lIqoymAAS/efAfws6Lzt6Fy7TD+CDVWeQm1++3rzI+FRFhqZ0ipz5S0REmo8JIFExImJfYuu5h+U6d+r2K4oNRomm+roqdOYvERFpPiaARCVYFnq3XOe9yMxRcCTKM6FTPXWHQEREKsYEkEgJBNWt21xhBnq66g6BiIhUjAkgUQm0KZErj4ru+UtERNqJCSBRCdKz82QeT0rPLvG8xLSS71elBR80K3JMX1cHHVytsW2stxoiIiIideNWcEQlyBcLEAQBIpFI6vi6U/dLPVdQY/fhb6O9MP/AbfTzdMJgr1poX98aq07cg1N1YzR3skDbulbQ0+X3PyKiqooJIFEp6sw4iPBZvrB643JpRjE9g29Kk6OMMpgZ6qGDqw0OT7WRHHOuYYwFHzRXSzxERKR52AVAJIeZf14v8znq6AB0d7bEjo/bqr5hIiLSKuwBJJLDkZsJUrefp5e+zMvnf1xVVjgy1bMxwV8T31Fpm0REpJ3YA0hUDhk5pV/ePXo7odQyivRR+zoqbY+IiLQXE0Cicrj+JEXdIUiZ2asRBrZyVncYRESkJXgJmKgcXmbmqjsEKR935G4eREQkP/YAEhEREVUxWpEArly5Ei4uLjA0NIS3tzcuXLhQYvmdO3eiUaNGMDQ0RLNmzXDw4EGp+/fs2YPu3bvDysoKIpEIV65ckVlPWFgYunTpAhMTE5ibm6Njx4549eqVoh4WkUIc/LSDukMgIiIto/EJ4I4dOxAYGIigoCBERETA3d0dfn5+SExMlFn+7NmzGDx4MEaPHo3Lly8jICAAAQEBuHHjhqRMRkYG2rdvj0WLFhXbblhYGHr06IHu3bvjwoULuHjxIiZNmgQdHY1/yqgK+bqXGxo7mKs7DCIi0jIiQZ3bFcjB29sbrVu3xooVKwAAYrEYzs7OmDx5MqZPn16k/MCBA5GRkYH9+/dLjrVp0wYeHh5Ys2aNVNmYmBjUqVMHly9fhoeHh9R9bdq0Qbdu3TBv3rxyx56amgoLCwukpKTA3Jz/pDWdy/QDJd4fs9Bf7rKqsri/O/p5Oqk7DCKtJQgCUlJSYGFhUWTHHyJtU5a8Q6MngeTk5CA8PBwzZsyQHNPR0YGvry/CwsJknhMWFobAwECpY35+fti7d6/c7SYmJuL8+fMYOnQo2rVrh+joaDRq1Ajz589H+/btiz0vOzsb2dmv94BNTU0FUPABo+F5dpUnz99HE/+G3RrX1Mi4iLRF4ecz30dUGZTldazRCWBSUhLy8/Nha2srddzW1haRkZEyz4mPj5dZPj4+Xu52798v2Of1m2++weLFi+Hh4YEtW7aga9euuHHjBlxdXWWet2DBAsydO7fI8ZSUFH64aDh5/j4pKZqz9MtP/dzgYmUEITsTKdmllyci2QRBQHp6OgCwB5C0XmHHkzw0OgFUF7FYDAAYN24cRo0aBQBo0aIFQkNDsWHDBixYsEDmeTNmzJDqfUxNTYWzszMsLCx4CViDicUCfjgSVWo5E1Mz6OmqZgxog5qmCGjhiOFta6PpN/9I3dfIzgy9PbnoM5EiFH754yVgqgzK8hrW6ATQ2toaurq6SEiQ3lEhISEBdnZ2Ms+xs7MrU3lZ7O3tAQCNGzeWOu7m5obY2NhizzMwMICBgUGR4yKRiB8sGmzvlSdYe+p+qeUEKPfvGLPQHxnZubj7+Bnc69pL2opZ6A9BEPDvvSSIIEJ7V2ulxUBUFRV+RvNzmrRdWV7DGj2lVV9fH56enggNDZUcE4vFCA0NRdu2sje8b9u2rVR5AAgJCSm2vCwuLi5wcHBAVJR0r9CdO3dQu3btMjwC0gY7wx/JVW5LWAwGrAlD+MMXSovFWF8PLlbGRY6LRCJ0cLVh8kdERAqh0T2AABAYGIgRI0agVatW8PLywrJly5CRkSG5NDt8+HA4OjpKLstOmTIFPj4+WLJkCfz9/bF9+3ZcunQJ69atk9SZnJyM2NhYxMXFAYAk0bOzs4OdnR1EIhGmTZuGoKAguLu7w8PDA5s3b0ZkZCR27dql4meAlOleYjrO3U+Wq+x3B24DAPquPqvMkIiIiJRO4xPAgQMH4tmzZ5gzZw7i4+Ph4eGBw4cPSyZ6xMbGSq3N165dOwQHB2PWrFmYOXMmXF1dsXfvXjRt2lRSZt++fZIEEgAGDRoEAAgKCsI333wDAJg6dSqysrLw2WefITk5Ge7u7ggJCUG9etxyq7KIfZ6Jbv87qe4wAAANbc3UHQIREVUhGr8OoDbjOoCa7ZOt4Th4Xf7Z4coU4OGAZYNacE0yIhXje44qk7LkHRo9BpBIXk9evsLhG/HIzsuX+xx+9SEioqpK4y8BE5XkxpMUfHfglmQc33ifepjes5Gaoyo79jwQEZEqsQeQtNr7K89ITeJYczJa7nM1qQfQwdJQ3SEQEVEVwgSQtFq+WP4sLjEtCxv+fYAHSRlKjKh8Jnaur+4QiIioCuElYKoyxm4Jx9VHL/HDkUhEzuuJqIQ0tcbTp4UjzA31MKC1M4z1+VYkIiLV4X8d0lq5+eIylb/66CUAICu34Dx19wS2qGWJ4W1d1BoDERFVTUwASWtlZss34/dlZg52hT9WcjRlN6CVs7pDICKiKooJIFV6n/9xFaGRieoOQ8qtb/1gWE1X3WEQEVEVxUkgVOlpWvLXq5kdx/wREZFaMQEkUrGg3k3UHQIREVVxTACJVMjFyhi25lzzj4iI1IsJIFVKt+JSsSUsBimvctUdisRUX1ecmNZZ3WEQERFxEghVPouPRGHF8XsAgAsPkmWW+S0sRoURAaPb18FU3wYqbZOIiKg47AGkSqcw+QOA/deeyiwz+6+bqgoHPZvaaeX+xEREVHmxB5BIyVZ/6KnuEIiIiKSwB5BIiQI8HNQdAhERURFMAImUaNmgFuoOgYiIqAgmgERERERVDMcAEimBXxNbjPOpp+4wiIiIZGICSKRgdaxNsHZYK3WHQUREVCxeAiatdemh7DX+1GH10JaS35cMcFdjJERERKVjDyBprbPRz9UdAgBg48jW6NyoJo5/0QlZuflwszdXd0hEREQlYgJIWkssCOoOARdmdkXN//b2rWNtouZoiIiI5MNLwKS1ElOz1dr++TeSPyIiIm3CBJC0kiAIOHBd9jZvqvBxx7qwZfJHRERaigkgaaXDN+LV1raDhSFm9nJTW/tEREQVxQSQtNKtp6lqa/vM9C5qa5uIiEgRmACSVlLnBBCRSKS2tomIiBSBCSBpFUEQsOzoHaw8Hq3UdnaObyv5fUz7OjAzLJgw71WnhlLbJSIiUgUuA0Na5ejtRCw7elfp7VibGuD0l51x7XEKurrVxPC2Ljh5JxG9mtkrvW0iIiJlYwJIapcvFpCckQMbM4Mi973MzMHkbZdhVE0XPw9pgWORiSqJqXBNP+caxgCAWlbGGNbWRSVtExERKRsTQFIrQRBQb+ZBye2783uimu7rkQnfHbiN03eTAACD1p3D5diXqg6RiIio0uEYQFKryPg0qduuXx/C8/TXCzyHvbHdG5M/IiIixWACSGqVkycucmzW3hsACnoHn7x8peqQMMufa/wREVHlxgSQ1Coxreh2bof+W+Q5Qk09fqPeqaOWdomIiFSFYwBJrX48EinzeLOgI9DVVf16e8sHeUBXh+v8ERFR5cYEkNQqIztf5vG07DwVRwLs+LgN1/kjIqIqQSsuAa9cuRIuLi4wNDSEt7c3Lly4UGL5nTt3olGjRjA0NESzZs1w8OBBqfv37NmD7t27w8rKCiKRCFeuXClSR6dOnSASiaR+xo8fr8iHVak9fJ6BzWdjkJRe9BKvpjE10EPMQn9417XiLh9ERFQlaHwCuGPHDgQGBiIoKAgRERFwd3eHn58fEhNlrwd39uxZDB48GKNHj8bly5cREBCAgIAA3LhxQ1ImIyMD7du3x6JFi0pse+zYsXj69Knk54cfflDoY6vMei0/jaB9N/HJ1ogSy73MzFFRRMCGka2kbu+b9A7uf98LN+b6qSwGIiIiTaDxl4CXLl2KsWPHYtSoUQCANWvW4MCBA9iwYQOmT59epPzy5cvRo0cPTJs2DQAwb948hISEYMWKFVizZg0AYNiwYQCAmJiYEts2NjaGnZ2dAh9N5ScIAi4/eomMnIJLuxceJJdYvrCcKnRpZIvDUzvg6qOX6OBqAwdLI5W1TUREpEk0OgHMyclBeHg4ZsyYITmmo6MDX19fhIWFyTwnLCwMgYGBUsf8/Pywd+/eMre/detW/P7777Czs0Pv3r0xe/ZsGBsbF1s+Ozsb2dmvL3mmpqYCKEiKBEEoc/va6M/LTxD4x1WpY+ExyWhZuzryxQJ0RJBcZlXlc/LnJ+0gCAIa2pqhoa2ZytuXV+FrRRNjI6qM+J6jyqQsr2ONTgCTkpKQn58PW1tbqeO2traIjJQ9ezQ+Pl5m+fj4+DK1PWTIENSuXRsODg64du0avvrqK0RFRWHPnj3FnrNgwQLMnTu3yPGUlJRK8eESmZCOF5l58HaxgE4xY+XeTv4AoO8a2cm6KvzUzw2talnAWF+ElJQUtcUhL0EQkJ6eDgAcj0ikAnzPUWVS2PEkD41OANXp448/lvzerFkz2Nvbo2vXroiOjka9evVknjNjxgyp3sfU1FQ4OzvDwsIC5ubmCo/x3P3nyM1/vZDy2znmmzdf5eTjWVo2nGsYSX/ICYAAATl5YiSmZSPs/nPoikRo7GCO+JQsnLqbBFtzA5y7X/KlXE2w4+M2+Gr3NeTmCzj+hY/UlnLaovCLgoWFBf8ZEakA33NUmZTlNazRCaC1tTV0dXWRkJAgdTwhIaHYsXl2dnZlKi8vb29vAMC9e/eKTQANDAxgYGBQ5HjhLGJFG/97BFJe5Sq8XgD4+9pTye8PkjKU0oai7BrfFq1cCpZvOTGts5qjqbg3Z54TkfLxPUeVRaVJAPX19eHp6YnQ0FAEBAQAAMRiMUJDQzFp0iSZ57Rt2xahoaGYOnWq5FhISAjatm1boVgKl4qxt7evUD1UftP8GiIjOw/NHC3Qrp41zI0KXr780CYiIiobjU4AASAwMBAjRoxAq1at4OXlhWXLliEjI0MyK3j48OFwdHTEggULAABTpkyBj48PlixZAn9/f2zfvh2XLl3CunXrJHUmJycjNjYWcXFxAICoqCgABb2HdnZ2iI6ORnBwMHr16gUrKytcu3YNn332GTp27IjmzZur+Bko3scd6yI7961ZtG8lQ4W3ouLTcPNpCno0sYORvp7UfSIR8PjFK5x/8ByPklW/9+4Q71qY0tUVtuaGyM0XI18swEBPh4kdERGRkmh8Ajhw4EA8e/YMc+bMQXx8PDw8PHD48GHJRI/Y2Fjo6Lwe69WuXTsEBwdj1qxZmDlzJlxdXbF37140bdpUUmbfvn2SBBIABg0aBAAICgrCN998A319fRw9elSSbDo7O6Nv376YNWuWih61fCZ2rq/uEBSumq4OqumqOwoiIqLKTSRUhumpGio1NRUWFhZISUlRyiQQqnwEQUBKSgoHpBOpCN9zVJmUJe/QvmmSRERERFQhTACJiIiIqhgmgERERERVDBNAIiIioiqGCSARERFRFcMEkIiIiKiKYQJIREREVMUwASQiIiKqYpgAEhEREVUxGr8VnDYr3GQlNTVVzZGQthAEAampqRCJRNyVgEgF+J6jyqQw35BnkzcmgEqUlpYGAHB2dlZzJERERFRVpKWlwcLCosQy3AtYicRiMeLi4mBmZqa0b5atW7fGxYsXlVK3KttTRL0VqaM855blHHnLpqamwtnZGY8ePary+0er+rVdFnzfKaYOTXjf8T0nje875benzPedIAhIS0uDg4MDdHRKHuXHHkAl0tHRgZOTk1Lb0NXVVemHlrLaU0S9FamjPOeW5Zyy1m9ubl7l/xmp+rVdFnzfKaYOTXrf8T1XgO875ben7PddaT1/hTgJRMtNnDixUrSniHorUkd5zi3LOar+O1UGmvyc8X2nmDr4vtM8mvyc8X2n2Dp4CZhIg6SmpsLCwgIpKSka+y2cqDLhe46qKvYAEmkQAwMDBAUFwcDAQN2hEFUJfM9RVcUeQCIiIqIqhj2ARERERFUME0AiIiKiKoYJIBEREVEVwwSQiIiIqIphAkhERERUxTABJNIS+/fvR8OGDeHq6opff/1V3eEQVQl9+vRB9erV0a9fP3WHQqRQXAaGSAvk5eWhcePGOH78OCwsLODp6YmzZ8/CyspK3aERVWonTpxAWloaNm/ejF27dqk7HCKFYQ8gkRa4cOECmjRpAkdHR5iamqJnz574559/1B0WUaXXqVMnmJmZqTsMIoVjAkikAqdOnULv3r3h4OAAkUiEvXv3FimzcuVKuLi4wNDQEN7e3rhw4YLkvri4ODg6OkpuOzo64smTJ6oInUhrVfR9R1SZMQEkUoGMjAy4u7tj5cqVMu/fsWMHAgMDERQUhIiICLi7u8PPzw+JiYkqjpSo8uD7jqh4TACJVKBnz5747rvv0KdPH5n3L126FGPHjsWoUaPQuHFjrFmzBsbGxtiwYQMAwMHBQarH78mTJ3BwcFBJ7ETaqqLvO6LKjAkgkZrl5OQgPDwcvr6+kmM6Ojrw9fVFWFgYAMDLyws3btzAkydPkJ6ejkOHDsHPz09dIRNpPXned0SVmZ66AyCq6pKSkpCfnw9bW1up47a2toiMjAQA6OnpYcmSJejcuTPEYjG+/PJLzgAmqgB53ncA4Ovri6tXryIjIwNOTk7YuXMn2rZtq+pwiRSOCSCRlnjvvffw3nvvqTsMoirl6NGj6g6BSCl4CZhIzaytraGrq4uEhASp4wkJCbCzs1NTVESVG993VNUxASRSM319fXh6eiI0NFRyTCwWIzQ0lJeaiJSE7zuq6ngJmEgF0tPTce/ePcntBw8e4MqVK6hRowZq1aqFwMBAjBgxAq1atYKXlxeWLVuGjIwMjBo1So1RE2k3vu+Iiset4IhU4MSJE+jcuXOR4yNGjMCmTZsAACtWrMCPP/6I+Ph4eHh44KeffoK3t7eKIyWqPPi+IyoeE0AiIiKiKoZjAImIiIiqGCaARERERFUME0AiIiKiKoYJIBEREVEVwwSQiIiIqIphAkhERERUxTABJCIiIqpimAASERERVTFMAImIKuibb76Bh4eHusOQEIlE2Lt3b5nPi4qKgp2dHdLS0hQf1BuSkpJQs2ZNPH78WKntEFHxmAASkVZYs2YNzMzMkJeXJzmWnp6OatWqoVOnTlJlT5w4AZFIhOjoaBVHqVqKTjxnzJiByZMnw8zMTGF1ymJtbY3hw4cjKChIqe0QUfGYABKRVujcuTPS09Nx6dIlybHTp0/Dzs4O58+fR1ZWluT48ePHUatWLdSrV08doWql2NhY7N+/HyNHjlRJe6NGjcLWrVuRnJyskvaISBoTQCLSCg0bNoS9vT1OnDghOXbixAm8//77qFOnDs6dOyd1vHPnzgCA3377Da1atYKZmRns7OwwZMgQJCYmAgDEYjGcnJywevVqqbYuX74MHR0dPHz4EADw8uVLjBkzBjY2NjA3N0eXLl1w9erVEuP99ddf4ebmBkNDQzRq1AirVq2S3BcTEwORSIQ9e/agc+fOMDY2hru7O8LCwqTq+OWXX+Ds7AxjY2P06dMHS5cuhaWlJQBg06ZNmDt3Lq5evQqRSASRSIRNmzZJzk1KSkKfPn1gbGwMV1dX7Nu3r8R4//jjD7i7u8PR0VFyTFYP47Jly+Di4iK5PXLkSAQEBOD777+Hra0tLC0t8e233yIvLw/Tpk1DjRo14OTkhI0bN0rV06RJEzg4OODPP/8sMS4iUg4mgESkNTp37ozjx49Lbh8/fhydOnWCj4+P5PirV69w/vx5SQKYm5uLefPm4erVq9i7dy9iYmIkvVw6OjoYPHgwgoODpdrZunUr3nnnHdSuXRsA0L9/fyQmJuLQoUMIDw9Hy5Yt0bVr12J7r7Zu3Yo5c+Zg/vz5uH37Nr7//nvMnj0bmzdvlir39ddf44svvsCVK1fQoEEDDB48WHKJ+8yZMxg/fjymTJmCK1euoFu3bpg/f77k3IEDB+Lzzz9HkyZN8PTpUzx9+hQDBw6U3D937lwMGDAA165dQ69evTB06NASe9tOnz6NVq1alfj8F+fYsWOIi4vDqVOnsHTpUgQFBeHdd99F9erVcf78eYwfPx7jxo0rMubPy8sLp0+fLlebRFRBAhGRlvjll18EExMTITc3V0hNTRX09PSExMREITg4WOjYsaMgCIIQGhoqABAePnwos46LFy8KAIS0tDRBEATh8uXLgkgkkpTPz88XHB0dhdWrVwuCIAinT58WzM3NhaysLKl66tWrJ6xdu1YQBEEICgoS3N3dpe4LDg6WKj9v3jyhbdu2giAIwoMHDwQAwq+//iq5/+bNmwIA4fbt24IgCMLAgQMFf39/qTqGDh0qWFhYSG6/3W4hAMKsWbMkt9PT0wUAwqFDh2Q+J4IgCO7u7sK3334rdUxW/f/73/+E2rVrS26PGDFCqF27tpCfny851rBhQ6FDhw6S23l5eYKJiYmwbds2qbo+++wzoVOnTsXGRETKwx5AItIanTp1QkZGBi5evIjTp0+jQYMGsLGxgY+Pj2Qc4IkTJ1C3bl3UqlULABAeHo7evXujVq1aMDMzg4+PD4CCMW8A4OHhATc3N0kv4MmTJ5GYmIj+/fsDAK5evYr09HRYWVnB1NRU8vPgwQOZk0wyMjIQHR2N0aNHS5X/7rvvipRv3ry55Hd7e3sAkFyejoqKgpeXl1T5t2+X5M26TUxMYG5uLqlbllevXsHQ0FDu+t/UpEkT6Oi8/ndia2uLZs2aSW7r6urCysqqSPtGRkbIzMwsV5tEVDF66g6AiEhe9evXh5OTE44fP44XL15IkjkHBwc4Ozvj7NmzOH78OLp06QKgIBnz8/ODn58ftm7dChsbG8TGxsLPzw85OTmSeocOHYrg4GBMnz4dwcHB6NGjB6ysrAAUzDR+e+xhocLxeG9KT08HUDB+z9vbW+o+XV1dqdvVqlWT/C4SiQAUjEtUhDfrLqy/pLqtra3x4sULqWM6OjoQBEHqWG5urlxtydN+cnIybGxsin8QRKQ0TACJSKt07twZJ06cwIsXLzBt2jTJ8Y4dO+LQoUO4cOECJkyYAACIjIzE8+fPsXDhQjg7OwOA1CziQkOGDMGsWbMQHh6OXbt2Yc2aNZL7WrZsifj4eOjp6UlNfiiOra0tHBwccP/+fQwdOrTcj7Nhw4a4ePGi1LG3b+vr6yM/P7/cbbypRYsWuHXrltQxGxsbxMfHQxAESYJ65coVhbQHADdu3CiyhA8RqQYvARORVuncuTP+/fdfXLlyRdIDCAA+Pj5Yu3YtcnJyJBNAatWqBX19ffz888+4f/8+9u3bh3nz5hWp08XFBe3atcPo0aORn5+P9957T3Kfr68v2rZti4CAAPzzzz+IiYnB2bNn8fXXX8tMJoGCCRgLFizATz/9hDt37uD69evYuHEjli5dKvfjnDx5Mg4ePIilS5fi7t27WLt2LQ4dOiRJxArjfvDgAa5cuYKkpCRkZ2fLXf/b/Pz8EBYWJpVQdurUCc+ePcMPP/yA6OhorFy5EocOHSp3G2/KzMxEeHg4unfvrpD6iKhsmAASkVbp3LkzXr16hfr168PW1lZy3MfHB2lpaZLlYoCCHqxNmzZh586daNy4MRYuXIjFixfLrHfo0KG4evUq+vTpAyMjI8lxkUiEgwcPomPHjhg1ahQaNGiAQYMG4eHDh1Ltv2nMmDH49ddfsXHjRjRr1gw+Pj7YtGkT6tSpI/fjfOedd7BmzRosXboU7u7uOHz4MD777DOpcXp9+/ZFjx490LlzZ9jY2GDbtm1y1/+2nj17Qk9PD0ePHpUcc3Nzw6pVq7By5Uq4u7vjwoUL+OKLL8rdxpv++usv1KpVCx06dFBIfURUNiLh7QEeRESkkcaOHYvIyEilLZ2ycuVK7Nu3D0eOHFFK/W9q06YNPv30UwwZMkTpbRFRURwDSESkoRYvXoxu3brBxMQEhw4dwubNm6UWlFa0cePG4eXLl0hLS1PqdnBJSUn44IMPMHjwYKW1QUQlYw8gEZGGGjBgAE6cOIG0tDTUrVsXkydPxvjx49UdFhFVAkwAiYiIiKoYTgIhIiIiqmKYABIRERFVMUwAiYiIiKoYJoBEREREVQwTQCIiIqIqhgkgERERURXDBJCIiIioimECSERERFTFMAEkIiIiqmL+Dyha6OGDKZMXAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "import matplotlib.pyplot as plt\n", "\n", "plt.figure(figsize=(7, 4))\n", "plt.plot(wlgrid, abs_only_rprs, lw=2, label='Absorption only')\n", "plt.xscale('log')\n", "plt.xlabel('Wavelength (um)')\n", "plt.ylabel('$(R_p/R_s)^2$')\n", "plt.title('Transmission spectrum — absorption only')\n", "plt.legend()\n", "plt.grid(alpha=0.2)" ] }, { "cell_type": "markdown", "id": "5337a83e", "metadata": {}, "source": [ "## CIA and Rayleigh Scattering\n", "\n", "With absorption-only features in place, two more fundamental processes fill in the spectral baseline. Collision-Induced Absorption (CIA) from H₂–H₂ and H₂–He pairs adds a smooth continuum at long wavelengths, while Rayleigh scattering raises the transit depth in the blue. Each is added as a separate contribution and the model is rebuilt to register it.\n", "\n", "[Cloud models](https://escience-taurex.github.io/taurex3/examples/05_clouds.html) extends this model with cloud opacities. More information about available contributions is [here](../user/taurex/models.rst)." ] }, { "cell_type": "markdown", "id": "f35b863e", "metadata": {}, "source": [ "### Composite transmission regions\n", "\n", "TauREx now includes the former `taurex-multimodel` plugin directly in core. If you need a transmission spectrum built from multiple 1D regions, use `model_type = multi_transit` in a parameter file and provide one regional parameter file per atmosphere in `parfiles`. The optional `fractions` list sets the weighting of each region, and the last fraction can be inferred automatically by providing only `N-1` values." ] }, { "cell_type": "code", "execution_count": null, "id": "2eea8fb9", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Contributions: ['Absorption', 'CIA', 'Rayleigh']\n" ] } ], "source": [ "from taurex.cache import CIACache\n", "from taurex.contributions import CIAContribution, RayleighContribution\n", "from _shared import CIA_DIR, TMP_DIR\n", "\n", "clean_cia_dir = TMP_DIR / 'cia_clean'\n", "clean_cia_dir.mkdir(parents=True, exist_ok=True)\n", "\n", "for source_path in CIA_DIR.glob('*.cia'):\n", " cleaned = [\n", " f\"H2-H2 {line.split(maxsplit=3)[-1]}\"\n", " if line.lstrip().startswith('eq-H2 -- eq-H2')\n", " else line\n", " for line in source_path.read_text().splitlines(keepends=True)\n", " ]\n", " (clean_cia_dir / source_path.name).write_text(''.join(cleaned))\n", "\n", "cia_cache = CIACache()\n", "cia_cache.cia_dict.clear()\n", "cia_cache.set_cia_path(str(clean_cia_dir))\n", "\n", "# Keep the cell rerunnable by starting from absorption only.\n", "tm.contribution_list = tm.contribution_list[:1]\n", "\n", "spectra = {}\n", "for contribution, key in [\n", " (CIAContribution(cia_pairs=['H2-H2', 'H2-He']), 'cia'),\n", " (RayleighContribution(), 'full'),\n", "]:\n", " tm.add_contribution(contribution)\n", " tm.build()\n", " spectra[key] = tm.model()[1][::-1]\n", "\n", "rprs_cia = spectra['cia']\n", "rprs_full = spectra['full']\n", "\n", "print('Contributions:', [c.name for c in tm.contribution_list])" ] }, { "cell_type": "code", "execution_count": 6, "id": "ca2f9bc7", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoAAAAGOCAYAAAAHL7LOAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAj2hJREFUeJzs3XlczNv/B/DXzNTMKISkRGQJIUVELgpR1pubLa5ws15c5Npd2beLaynrpWzZQnxdXGS9ZKlEWbILqSRNm7aZz++Pfg1jppqZZqvez8fjPm7z+ZzPOe8ZM9O78zkLi2EYBoQQQgghpMJgazsAQgghhBCiWZQAEkIIIYRUMJQAEkIIIYRUMJQAEkIIIYRUMJQAEkIIIYRUMJQAEkIIIYRUMJQAEkIIIYRUMJQAEkIIIYRUMJQAEkIIIYRUMJQAElKOZWRkYMyYMTAzMwOLxcK0adMAAImJiRg4cCCMjY3BYrGwYcMGrcapiKKekzpYWlpi1KhRJZYLDAwEi8XC69evxcecnZ3h7OysttjKi+9fpytXroDFYiE4OFgj7Y8aNQqWlpYaaYsQXUIJICFlTGGyUdR/t27dEpddsWIFAgMDMXHiROzbtw8jRowAAEyfPh3//vsv5s6di3379sHNzU3lca5YsQIhISFqqVfWcyLaFR8fj0WLFiEqKkrboUjR5dgI0RY9bQdACFHOkiVL0KBBA6njjRs3Fv986dIldOjQAb6+vhJlLl26hB9//BG///672uJbsWIFBg4cCHd3d5XWW9Rz0jXnz5/XdggaFR8fj8WLF8PS0hJ2dnZyX6eJ16m42Hbu3AmRSKT2GAjRNZQAElJG9erVC23bti22TFJSEpo3by7zeLVq1dQUmXoV9Zx0DZfL1XYIOi0rKwsGBgZaf5309fW12j4h2kK3gAkphwrHUb169Qr//POP+PZw4e1jhmHg7+8vPl4oNTUV06ZNg4WFBXg8Hho3bozVq1dL9ZCIRCJs3LgRNjY24PP5MDExgZubG8LDwwEALBYLmZmZ2LNnj7iNksbSJSUlwdvbG6ampuDz+bC1tcWePXtKfE7fjrv7XkBAALp164ZatWqBx+OhefPm2Lp1q1Q5hmGwbNky1K1bFwYGBujatSsePnwos86HDx+iW7duqFSpEurWrYtly5bJ7EEqamzbkSNHsHz5ctStWxd8Ph/du3fH8+fPpa739/dHw4YNUalSJTg4OOD69esKjSvcv38/HBwcYGBggOrVq6NLly5SvW1btmxBixYtwOPxYG5ujkmTJiE1NVXqebRs2RKPHj1C165dYWBggDp16mDNmjUSz61du3YAgNGjR0u8376tIyIiAl26dIGBgQHmzZsn83UqJBQKMW/ePJiZmcHQ0BD9+/fH27dvJcoUNUbz2zpLik3WGMDMzEzMmDFD/Dlo2rQp1q5dC4ZhJMqxWCxMnjwZISEhaNmyJXg8Hlq0aIFz585JlEtPT8e0adNgaWkJHo+HWrVqoUePHoiMjJSKnRBNoR5AQsoogUCA5ORkiWMsFgvGxsawtrbGvn37MH36dNStWxczZswAALRu3Vo8bq5Hjx7w8vISX5uVlQUnJye8f/8e48ePR7169XDz5k3MnTsXHz58kJgo4u3tjcDAQPTq1QtjxoxBfn4+rl+/jlu3bqFt27bYt28fxowZAwcHB4wbNw4A0KhRoyKfy5cvX+Ds7Iznz59j8uTJaNCgAY4ePYpRo0YhNTUVU6dOLfI5mZiYFFnv1q1b0aJFC/Tv3x96enr43//+h19//RUikQiTJk0Sl1u4cCGWLVuG3r17o3fv3oiMjETPnj2Rm5srUV9CQgK6du2K/Px8zJkzB4aGhtixYwcqVapUwr/WV6tWrQKbzcbvv/8OgUCANWvWYPjw4bh9+7ZE3JMnT0bnzp0xffp0vH79Gu7u7qhevTrq1q1bYhuLFy/GokWL0LFjRyxZsgRcLhe3b9/GpUuX0LNnTwDAokWLsHjxYri4uGDixImIjY3F1q1bcffuXdy4cUOiZ+zz589wc3PDTz/9hMGDByM4OBizZ8+GjY0NevXqBWtrayxZsgQLFy7EuHHj0LlzZwBAx44dxXV8+vQJvXr1wtChQ/Hzzz/D1NS02OewfPlysFgszJ49G0lJSdiwYQNcXFwQFRWl0OstT2zfYhgG/fv3x+XLl+Ht7Q07Ozv8+++/mDlzJt6/f4+//vpLovx///2H48eP49dff0WVKlWwadMmeHh4IC4uDsbGxgCACRMmIDg4GJMnT0bz5s3x6dMn/Pfff3j8+DHatGkj93MhRKUYQkiZEhAQwACQ+R+Px5MoW79+faZPnz5SdQBgJk2aJHFs6dKljKGhIfP06VOJ43PmzGE4HA4TFxfHMAzDXLp0iQHA/Pbbb1L1ikQi8c+GhobMyJEj5XpOGzZsYAAw+/fvFx/Lzc1lHB0dmcqVKzNpaWklPidZsrKypI65uroyDRs2FD9OSkpiuFwu06dPH4n4582bxwCQeA7Tpk1jADC3b9+WuN7IyIgBwLx69Up83MnJiXFychI/vnz5MgOAsba2ZnJycsTHN27cyABgoqOjGYZhmJycHMbY2Jhp164dk5eXJy4XGBjIAJCoU5Znz54xbDabGTBgACMUCiXOFT6/wufcs2dPiTJ+fn4MAGb37t0SzwMAs3fvXvGxnJwcxszMjPHw8BAfu3v3LgOACQgIkIqpsI5t27bJPCfrdapTp47Ev/uRI0cYAMzGjRvFx+rXry/zPfZ9ncXFNnLkSKZ+/frixyEhIQwAZtmyZRLlBg4cyLBYLOb58+fiYwAYLpcrcez+/fsMAGbz5s3iY0ZGRlKfN0K0jW4BE1JG+fv748KFCxL/nT17Vun6jh49is6dO6N69epITk4W/+fi4gKhUIhr164BAI4dOwYWiyVzEsa3t5MVcebMGZiZmcHT01N8TF9fH7/99hsyMjJw9epVper9tqeosMfUyckJL1++hEAgAABcvHgRubm5mDJlikT8spaXOXPmDDp06AAHBwfxMRMTEwwfPlzumEaPHi0x7q2wR+rly5cAgPDwcHz69Aljx46Fnt7XmzTDhw9H9erVS6w/JCQEIpEICxcuBJst+RVf+PwKn/O0adMkyowdOxZVq1bFP//8I3Fd5cqV8fPPP4sfc7lcODg4iGOWB4/Hw+jRo+Uu7+XlhSpVqogfDxw4ELVr18aZM2fkrkMZZ86cAYfDwW+//SZxfMaMGWAYRuoz5uLiItG73apVK1StWlXitalWrRpu376N+Ph4tcZOiCLoFjAhZZSDg0OJk0AU8ezZMzx48KDIW6pJSUkAgBcvXsDc3Bw1atRQWdtv3ryBlZWVVMJibW0tPq+MGzduwNfXF2FhYcjKypI4JxAIYGRkJK7byspK4ryJiYlUwvXmzRu0b99eqp2mTZvKHVO9evUkHhe28fnzZ3EbgORsbgDQ09OTa726Fy9egM1mFztRprCN7+Pmcrlo2LCh1Otdt25dqeS+evXqePDgQYnxFKpTp45CEz6+//dgsVho3LhxsWM+VeHNmzcwNzeXSD6Bot+L3/97AgWvTeG/JwCsWbMGI0eOhIWFBezt7dG7d294eXmhYcOGangGhMiHEkBCCICCiR09evTArFmzZJ5v0qSJhiMqnRcvXqB79+5o1qwZ1q9fDwsLC3C5XJw5cwZ//fWX1pb+4HA4Mo8z300w0CWqiFmRcXvyKqrHWSgUFhmzqsnz2gwePBidO3fGiRMncP78efz5559YvXo1jh8/jl69emkkTkK+R7eACSEACiZpZGRkwMXFReZ/hT0djRo1Qnx8PFJSUoqtT5HbwfXr18ezZ8+kkrInT56Izyvqf//7H3JycnDq1CmMHz8evXv3houLi1QiUlj3s2fPJI5//PhRohfn2zi/Fxsbq3B8RSmM5/uZwfn5+XL1fjVq1AgikQiPHj0qsY3v487NzcWrV6+Uer2Vvf1flO9fZ4Zh8Pz5c4le0OrVq0vNWgake+kUfS/Gx8cjPT1d4nhp3osAULt2bfz6668ICQnBq1evYGxsjOXLlytVFyGqQAkgIQRAQS9FWFgY/v33X6lzqampyM/PBwB4eHiAYRgsXrxYqty3vR6GhoYyfznL0rt3byQkJODw4cPiY/n5+di8eTMqV64MJycnBZ/N156Zb2MSCAQICAiQKOfi4gJ9fX1s3rxZoqys7fF69+6NW7du4c6dO+JjHz9+xIEDBxSOryht27aFsbExdu7cKX7NAeDAgQNSCaks7u7uYLPZWLJkiVRCXfj8XFxcwOVysWnTJonnvGvXLggEAvTp00fhuA0NDQFA7n/zkuzdu1ciCQsODsaHDx8keswaNWqEW7duSczWPn36tNRyMYrE1rt3bwiFQvj5+Ukc/+uvv8BisRTusRMKheLxpoVq1aoFc3Nz5OTkKFQXIapEt4AJKaPOnj0r7pX4VseOHZUaWzRz5kycOnUKffv2xahRo2Bvb4/MzExER0cjODgYr1+/Rs2aNdG1a1eMGDECmzZtwrNnz+Dm5gaRSITr16+ja9eumDx5MgDA3t4eFy9exPr162Fubo4GDRrIHD8HAOPGjcP27dsxatQoREREwNLSEsHBwbhx4wY2bNggNR5LHj179gSXy0W/fv0wfvx4ZGRkYOfOnahVqxY+fPggLmdiYoLff/8dK1euRN++fdG7d2/cu3cPZ8+eRc2aNSXqnDVrlnjrvKlTp4qXgalfv75C4+GKw+VysWjRIkyZMgXdunXD4MGD8fr1awQGBqJRo0Yl9mY1btwY8+fPx9KlS9G5c2f89NNP4PF4uHv3LszNzbFy5UqYmJhg7ty5WLx4Mdzc3NC/f3/ExsZiy5YtaNeuncSED3k1atQI1apVw7Zt21ClShUYGhqiffv2MnerkUeNGjXQqVMnjB49GomJidiwYQMaN26MsWPHisuMGTMGwcHBcHNzw+DBg/HixQvs379faskhRWLr168funbtivnz5+P169ewtbXF+fPncfLkSUybNq3Y5YxkSU9PR926dTFw4EDY2tqicuXKuHjxIu7evYt169Yp9doQohJamn1MCFFSccvA4LulLhRZBoZhGCY9PZ2ZO3cu07hxY4bL5TI1a9ZkOnbsyKxdu5bJzc0Vl8vPz2f+/PNPplmzZgyXy2VMTEyYXr16MREREeIyT548Ybp06cJUqlRJajkVWRITE5nRo0czNWvWZLhcLmNjYyNz2Q5FloE5deoU06pVK4bP5zOWlpbM6tWrmd27d0st2SIUCpnFixcztWvXZipVqsQ4OzszMTExMpcZefDgAePk5MTw+XymTp06zNKlS5ldu3bJvQzM0aNHJep79eqVzCVKNm3axNSvX5/h8XiMg4MDc+PGDcbe3p5xc3OT67nv3r2bad26NcPj8Zjq1aszTk5OzIULFyTK+Pn5Mc2aNWP09fUZU1NTZuLEicznz58lyjg5OTEtWrSQqv/75VMYhmFOnjzJNG/enNHT05N4TkXVUXhO1ut08OBBZu7cuUytWrWYSpUqMX369GHevHkjdf26deuYOnXqMDwej/nhhx+Y8PBwqTqLi03W80hPT2emT5/OmJubM/r6+oyVlRXz559/SiwTxDBFf46+fd/k5OQwM2fOZGxtbZkqVaowhoaGjK2tLbNlyxaZrwchmsJiGB0eeUwIIQRAwSQdExMT/PTTT9i5c6e2wyGElHE0BpAQQnRMdna21AzbvXv3IiUlRe6t4AghpDjUA0gIITrmypUrmD59OgYNGgRjY2NERkZi165dsLa2RkREhELr6RFCiCw0CYQQQnSMpaUlLCwssGnTJqSkpKBGjRrw8vLCqlWrKPkjhKhEmbgF7O/vD0tLS/D5fLRv315iCQZZjh49imbNmoHP58PGxkZi66C8vDzxJuaGhoYwNzeHl5eX1BY9y5cvR8eOHWFgYIBq1aqp42kRQohMlpaWOHXqFBISEpCbm4uEhATs3r0btWrV0nZohJByQucTwMOHD8PHxwe+vr6IjIyEra0tXF1dxdtSfe/mzZvw9PSEt7c37t27B3d3d7i7uyMmJgYAkJWVhcjISPzxxx+IjIzE8ePHERsbi/79+0vUk5ubi0GDBmHixIlqf46EEEIIIZqk82MA27dvj3bt2okX5RSJRLCwsMCUKVMwZ84cqfJDhgxBZmYmTp8+LT7WoUMH2NnZYdu2bTLbuHv3LhwcHPDmzRupfR0DAwMxbdo0lS1uSgghhBCibTo9BjA3NxcRERGYO3eu+BibzYaLiwvCwsJkXhMWFgYfHx+JY66urggJCSmyHYFAABaLpfJbvSKRCPHx8ahSpYrKt0kihBBCCPkWwzBIT0+Hubk52Ozib/LqdAKYnJwMoVAIU1NTieOmpqYyd0AAgISEBJnlExISZJbPzs7G7Nmz4enpiapVq5Yq3pycHImtfd6/f4/mzZuXqk5CCCGEEEW8ffsWdevWLbaMTieA6paXl4fBgweDYRhs3bq11PWtXLlS5v6oMTExSm1lRSoehmGQmZkJQ0ND6jUmRAPoM0fKk/T0dLRs2VKunEOnE8CaNWuCw+EgMTFR4nhiYiLMzMxkXmNmZiZX+cLk782bN7h06VKpe/8AYO7cuRK3n9PS0mBhYYG6deuqpH5S/jEMA4FAACMjI/plRIgG0GeOlCdpaWkAINd7WacTQC6XC3t7e4SGhsLd3R1Awbi60NBQ8Ybz33N0dERoaCimTZsmPnbhwgU4OjqKHxcmf8+ePcPly5dhbGysknh5PB54PJ7UcRaLRV8sRG6F7xd6zxCiGfSZI+WFIu9hnU4AAcDHxwcjR45E27Zt4eDggA0bNiAzMxOjR48GAHh5eaFOnTpYuXIlAGDq1KlwcnLCunXr0KdPHxw6dAjh4eHYsWMHgILkb+DAgYiMjMTp06chFArF4wNr1KghXmQ1Li4OKSkpiIuLg1AoRFRUFACgcePGqFy5soZfBUIIIYQQ1dH5BHDIkCH4+PEjFi5ciISEBNjZ2eHcuXPiiR5xcXESM106duyIoKAgLFiwAPPmzYOVlRVCQkLQsmVLAAUTM06dOgUAsLOzk2jr8uXL4n02Fy5ciD179ojPtW7dWqoMIYQQQkhZpPPrAJZlaWlpMDIygkAgKHYMoFAoRF5engYjI7qqcAq/LiwdxOVyS1xGgJCyjsYAkvJE3rwDKAM9gOUZwzBISEigRaaJBJFIhE+fPmk7DLDZbDRo0ID2niWEkHKIEkAtKkz+atWqBQMDA/rrk4BhGAiFQnA4HK2+HwoXMf/w4QPq1atH701CCClnKAHUEqFQKE7+VDULmZR9upIAAoCJiQni4+ORn58PfX19rcZCCCFEtWiAj5YUjvkzMDDQciSEyFZ461coFGo5EkIIIapGCaCWabuXh5Ci0HuTEELKL0oACSGEEKISmVnpuHn/rLbDIHKgBJCo3JUrV8BiscrU7OZRo0aJd5vRdZaWltiwYYO2wyCE6KD8/Dw8ffNAK22LhEIMP/ADxkfNwrJ9IzTe/sXbRzF0R2ss3uup8bbLIkoAiVLCwsLA4XDQp08fbYeikNevX4PFYol3dim0ceNGBAYGaiUmQghRBZFQiGG722Hg5WFYf1j2dqnqFPXsJl5wC5YWPiyK0nj78x8twkNePoKZGDx+GaHx9ssaSgCJUnbt2oUpU6bg2rVriI+P13Y4yM3NLdX1RkZGqFatmmqCIYQQLbjzKBSPeUIwLBYCsq9qvH2hsHTfw6WV9c3C9e+SnmsxkrKBEkCisIyMDBw+fBgTJ05Enz59iuw5u3HjBlq1agU+n48OHTogJiZGfO7Nmzfo168fqlevDkNDQ7Ro0QJnzpwRn7969SocHBzA4/FQu3ZtzJkzB/n5+eLzzs7OmDx5MqZNm4aaNWvC1dUVQMHEha1bt6JXr16oVKkSGjZsiODgYPF1DRo0AFCwtR+LxRJv6/f9LeCcnBz89ttvqFWrFvh8Pjp16oS7d++Kzxfe5g4NDUXbtm1hYGCAjh07IjY2ttjXLjo6Gt26dUOlSpVgbGyMcePGISMjQ3x+9OjR8PDwwNq1a1G7dm0YGxtj0qRJRe4U88svv6Bv374Sx/Ly8lCrVi3s2rWr2FgIIeVLTu4XbYdAyhBKAInCjhw5gmbNmqFp06b4+eefsXv3bsjaUXDmzJlYt24d7t69CxMTE/Tr10+cyEyaNAk5OTm4du0aoqOjsXr1alSuXBlAwX7NvXv3Rrt27XD//n1s3boVu3btwrJlyyTq37NnD7hcLm7cuIFt27aJj//xxx/w8PDA/fv3MXz4cAwdOhSPHz8GANy5cwcAcPHiRXz48AHHjx+X+RxnzZqFY8eOYc+ePYiMjETjxo3h6uqKlJQUiXLz58/HunXrEB4eDj09Pfzyyy9Fvm6ZmZlwdXVF9erVcffuXRw9ehQXL17E5MmSt2quXLmCFy9e4PLly9izZw8CAwOLTLLHjBmDc+fO4cOHD+Jjp0+fRlZWFoYMGVJkLIQQQio2Wghah/Tb/B8+pudovF2TKjz8b0onucvv2rULP//8MwDAzc0NAoEAV69eFfemFfL19UWPHj0AFCRrdevWxYkTJzB48GDExcXBw8MDNjY2AICGDRuKr9uyZQssLCzg5+cHFouFZs2aIT4+HrNnz8bChQvF+9NaWVlhzZo1UvENGjQIY8aMAQAsXboUFy5cwObNm7FlyxaYmJgAAIyNjWFmZibz+WVmZmLr1q0IDAxEr169AAA7d+7EhQsXsGvXLsycOVNcdvny5XBycgIAzJkzB3369EF2djb4fL5UvUFBQcjOzsbevXthaGgIAPDz80O/fv2wevVqmJqaAgCqV68OPz8/6OnpoVmzZujTpw9CQ0MxduxYqTo7duyIpk2bYt++fZg1axYAICAgAIMGDRIn1ISQiunIhU04+OpvtDVojfk/79F2OETHUAKoQz6m5yAhLVvbYRQrNjYWd+7cwYkTJwAAenp6GDJkCHbt2iWVADo6Oop/rlGjBpo2bSruifvtt98wceJEnD9/Hi4uLvDw8ECrVq0AAI8fP4ajo6PEOnQ//PADMjIy8O7dO9SrVw8AYG9vLzPGb9stfPz9pI/ivHjxAnl5efjhhx/Ex/T19eHg4CCOv1BhzABQu3ZtAEBSUpI4xm89fvwYtra24uSv8HmJRCLExsaKE8DmzZuDw+FI1BsdHV1kvGPGjMGOHTswa9YsJCYm4uzZs7h06ZLcz5cQUj4tjd8J8IDnwkhMz86EAd+w5ItIhUEJoA4xqcLT+XZ37dqF/Px8mJubi48xDAMejwc/Pz8YGRnJVc+YMWPg6uqKf/75B+fPn8fKlSuxbt06TJkyRe5Yvk2ktOXbLdIKE1aRSKSyOgvrLa5OLy8vzJkzB2FhYbh58yYaNGiAzp07lyoGQkj5kvjpLRiREA0tWmg7FKIjKAHUIYrchtWG/Px87N27F+vWrUPPnj0lzrm7u+PgwYOYMGGC+NitW7fEPWGfP3/G06dPYW1tLT5vYWGBCRMmYMKECZg7dy527tyJKVOmwNraGseOHQPDMOKk6saNG6hSpQrq1q1bYpy3bt2Cl5eXxOPWrVsDkG97s0aNGonHFtavXx9AwcSKu3fvYtq0aSW2XxRra2sEBgYiMzNTnLzeuHEDbDYbTZs2VbpeY2NjuLu7IyAgAGFhYRg9erTSdRFCyqdB5wdCxAJWNJkLt47DtR2O2rFYNMWhJPQKEbmdPn0anz9/hre3N1q2bCnxn4eHh9Ss0yVLliA0NBQxMTEYNWoUatasKZ5pO23aNPz777949eoVIiMjcfnyZXFy+Ouvv+Lt27eYMmUKnjx5gpMnT8LX1xc+Pj7i8X/FOXr0KHbv3o2nT5/C19cXd+7cEU+0qFWrFipVqoRz584hMTERAoFA6npDQ0NMnDgRM2fOxLlz5/Do0SOMHTsWWVlZ8Pb2Vvr1Gz58OPh8PkaOHImYmBhcvnwZU6ZMwYgRI8S3f5U1ZswY7NmzB48fP8bIkSNLVRchpPzJYbOQx2Jh7aMV2g5FzP/4TMwNcEdC8ltth1IhUQJI5LZr1y64uLjIvM3r4eGB8PBwPHjwdQX6VatWYerUqbC3t0dCQgL+97//SfTATZo0CdbW1nBzc0OTJk2wZcsWAECdOnVw5swZ3LlzB7a2tpgwYQK8vb2xYMECueJcvHgxDh06hFatWmHv3r04ePAgmjdvDqBgzOKmTZuwfft2mJub48cff5RZx6pVq+Dh4YERI0agTZs2eP78Of79919Ur15dodfsWwYGBvj333+RkpKCdu3aYeDAgejevTv8/PyUrrOQi4sLateuDVdXV4nb84SQioMtx/7dAo5u7PF9JfwEtqWfw2n2Cyw74VXyBUTlWIys9TuISqSlpcHIyAgCgQBVq1aVOJednY1Xr16hQYMGMmeMEuWwWCycOHGizGzr9j2GYSAUCsHhcCQmwZQkIyMDderUQUBAAH766SeVxELvUVIRMAwDgUAAIyMjhT5zuuh65Cn8Gj2/2DJ8EYO7o2OKLaOs8IdXMDr86zju6JFFT15btm+ExG4hxZWVl80eG/HPfzVbCJf2g0pdZ1lTXN7xPRoDSEgZJhKJkJycjHXr1qFatWro37+/tkMihBBSBlACSEgZFhcXhwYNGqBu3boIDAyEnh59pAkhhJSMfluQcqWijWiwtLSscM+ZEEJKIs94yIqOJoEQQgghFYS6/lz8+DkeO24UP/7wW/ki2fubE82hBJAQQgghpbLi2CiE8dLkLn9RqPikj/z8PKw5OA5rgsYiP58SyNKiBJAQQggpB7S5+PFF/Q8KlRdwFI91+6l52Jcbhn15t7D1xCyFryeSKAEkhBBCiM67lHxe/PPlz6FajKR8oASQEEII0UFXI07is+CjtsMg5RQlgIQQQoiOWRXkjckxCzDsaFeVj3d7E/+01HXk5uaopB51KeuLemsCJYCEEEKIjjmQdwcA8E6fhYt3jqis3hw2C30veGBBoIfSdYiEQgwObIcfz/+EHSfln/lLdAslgETjEhISMGXKFDRs2BA8Hg8WFhbo168fQkO/jumwtLTEhg0bpK5duXIlOBwO/vzzTw1GTAgh2iNUw5IpJ1nK994dv7INL3gMhCwWNqeeUllM+fl58NrRDr13tkTYg3PFlmWBevhKixJAUirOzs4IDAyUu/zr169hb2+PS5cu4c8//0R0dDTOnTuHrl27YtKkSSVev3v3bsyaNQu7d+8uRdSEEEKUlZ2boZZ6t5yYiXu8bLzlsrD89kyV1fspNUFldZUnlAASjfr111/BYrFw584deHh4oEmTJmjRogV8fHxw69atYq+9evUqvnz5giVLliAtLQ03b97UUNSEEEIUFf7wCrYcn4XU9GS5yidmxIl/fsOVPs8osYz1wj2D0C3EBTP+dpN5fum+n+G8qwXWH56scN1lHSWARGNSUlJw7tw5TJo0CYaGhlLnq1WrVuz1u3btgqenJ/T19eHp6Yldu3apKVJCCCHfEwmFOHJhExIEr+QqP/7OZGxNPwvvg92RnZOFT6kJSEh+q5JYnvJEmLCjEwQZKcWWO4EnELFYOK//Xub5I6L7+KTHRkD2VZXEVZbQXsC6ZLsTkJGk+XYr1wLGq//N//z5czAMg2bNmil8bVpaGoKDgxEWFgYA+Pnnn9G5c2ds3LgRlStXVnWohBBCvrMh+DcEZF+Tu3wuu2Cc3lOeCL8EdkI0v2As4wJzb5XEc4MnwOqj3lgx+oRK6gOAPf8sx8V3J9CnwTAM7emjsnp1ESWAuiQjCUiP13YUxVqxYgVWrFghfvzlyxfcunULkyd/7T5/9OgR6tWrJ3Utwyi/C+XBgwfRqFEj2NraAgDs7OxQv359HD58GN7eqvkyIYQQXSQSibQdAgAolPx9rzD5A4Bl8aq7exOb91LmcWV3RVmbfAjgA1EfAjAUlAASTalcS+fbnTBhAgYPHix+PHz4cHh4eOCnn34SHzM3N5d5rZWVFVgsFp48eaJwiLt27cLDhw+hp/f1LSsSibB7925KAAkhFVZ2Thb4PAOlrn3w7BbWX/oNZvpmWPWL6mbzkrKBEkBdooHbsKVVo0YN1KhRQ/y4UqVKqFWrFho3bizXta6urvD398dvv/0mNQ4wNTVV5jjA6OhohIeH48qVKxJtp6SkwNnZGU+ePFHqtjIhhJRlw7fb4wE/Fyb5IlzyfohHcWEKXb/s8gQ85gsBvELbUH8M7F7ySgyadD7sIHo6eoofK3IPSZv7IpcV9AoRjfL394dQKISDgwOOHTuGZ8+e4fHjx9i0aRMcHR1lXrNr1y44ODigS5cuaNmypfi/Ll26oF27djQZhBBS4cR/fIMH/FwAwEc9Nq6En1B4Tb7HPKH452cJd5WO5e2HZ0pfW5yd0aslHj/nKT+MiEijBJBoVMOGDREZGYmuXbtixowZaNmyJXr06IHQ0FBs3bpVqnxubi72798PDw/Zq9Z7eHhg7969yMtT/UKphBCiC+I/PUN2TpbEseTUDxKPH8Zpb1msxf/z0lrbyiwNQwqUmQTQ398flpaW4PP5aN++Pe7cuVNs+aNHj6JZs2bg8/mwsbHBmTNnxOfy8vIwe/Zs2NjYwNDQEObm5vDy8kJ8vOQEjJSUFAwfPhxVq1ZFtWrV4O3tjYwM9SyAWVZduXIFo0aNUuia2rVrw8/PD69fv0ZOTg7evXuHkydPwtnZWVzm9evXmDZtGrhcLpKTkzFzpuxFQWfNmoXExETo6+uX4lkQQoju8ks7jaF7OkAkFBZZ5p1Ae/vy3uaVzd+L75JeazsErSoTCeDhw4fh4+MDX19fREZGwtbWFq6urkhKkr1kys2bN+Hp6Qlvb2/cu3cP7u7ucHd3R0xMDAAgKysLkZGR+OOPPxAZGYnjx48jNjYW/fv3l6hn+PDhePjwIS5cuIDTp0/j2rVrGDdunNqfLyGEEPKtFzwGVyJUt9xJeVHUlnBsVslbxc0M+anEMop6E/8UE3Z0xvwA1detamUiAVy/fj3Gjh2L0aNHo3nz5ti2bRsMDAyK3A5s48aNcHNzw8yZM2FtbY2lS5eiTZs28PPzAwAYGRnhwoULGDx4MJo2bYoOHTrAz88PERERiIsrWIn88ePHOHfuHP7++2+0b98enTp1wubNm3Ho0CGpnkJCCCFEXvEf32BewABsD5mn0HV5wlw1RaR9+fmaH8YTw5O/zdT0ZIzY3hZDd7TG2wTJpWcOnf8LLn+3wLS/e2Dp6ZG4wUvFKfYzHDi3RtUhq5TOzwLOzc1FREQE5s6dKz7GZrPh4uIiXhT4e2FhYfDxkVy/x9XVFSEhIUW2IxAIwGKxxLNQw8LCUK1aNbRt21ZcxsXFBWw2G7dv38aAAQOk6sjJyUFOTo74cVpaGoCC9e++XwOv8LGsc4QApVs3UZXt03uUlGeF729NvseXnvgZ//FSAcFztHrgjA42PeS67vqTY1jzaDnsYYmhP8yWOFfa+EWMdr9zOu2zQxUhC2B/7bl7o58nV0yyyhT3e1feer59vDrYG1H8gt/vK06NxJaxX9dEXP5hN6DPRiIk9xyOib+BDUd+w/3PtzGizSx0bav+XkFF/g11PgFMTk6GUCiEqampxHFTU9Mi15NLSEiQWT4hQfaG0NnZ2Zg9ezY8PT1RtWpVcR21akmuj6enp4caNWoUWc/KlSuxePFiqeMCgUDqHyU3NxcikQhCoRDCYsZ1kIpHVxZ9FQqFEIlESE9Pl/jDhpDyhGEY8dhulhy3DVXhP16q+OdrD47Cup6DXNedZD8F2GycRRyE11cC3+yXKxKJAI7yMQmFQggEAuUrKKVMNhuZ392T/MJmlxgTA0ZmmS9fsqWOl1TX9+c3HJ6Ga2mX4VqzH97mvgP4Bcfv6n/G73/3Rmvz7uj7w5gi6/uSn4VdXy4DfOBxzEKcs+pebPuqUNjxJA+dTwDVLS8vD4MHDwbDMDJnoSpi7ty5Ej2PaWlpsLCwgJGRkTixLJSdnY1Pnz6Bw+GAwynFp5aUS7rwnuBwOGCz2ahSpQr4fL62wyFELQr/ODcyMtJYAvgtPT09GBkZ4VNqAm7c/we9f/CCnl7Jk9oSRJJj4Nns0o3o4nA4MDIyKlUd6iBPTLLKVKrElzpeUl3fn9+dcwngAc/TT8EWPPHxHDYL57nvcD55D/rrFb0RQTayxT9nstkaeX0VeQ/rfAJYs2ZNcDgcJCYmShxPTEyEmZmZzGvMzMzkKl+Y/L158waXLl2SSNLMzMykJpnk5+cjJSWlyHZ5PB54PJ7UcRaLJfWPUvhY1jlScX3bU6zt9wW9R0lFUfge19b7nBGJMCbYFS+4DG68DMGaX06XeI06ItXFz7k8MckqU9zv3dK09b2ET2/lLquJ11eRNnR+EgiXy4W9vT1CQ0PFx0QiEUJDQ4tcONjR0VGiPABcuHBBonxh8vfs2TNcvHgRxsbGUnWkpqYiIiJCfOzSpUsQiURo3769Kp4aIYQQgk9piXjBLfjj7yznjVzXJHCyJR6rcj28rOxMPI+LKXbZGV1R1Cxg2gmkZDrfAwgAPj4+GDlyJNq2bQsHBwds2LABmZmZGD16NADAy8sLderUwcqVKwEAU6dOhZOTE9atW4c+ffrg0KFDCA8Px44dOwAUJH8DBw5EZGQkTp8+DaFQKB7XV6NGDXC5XFhbW8PNzQ1jx47Ftm3bkJeXh8mTJ2Po0KFF7nVLCCGEKEokUjx5S9SXTHAYprRjhxmIhEL8vKsdov9/dmzHHCN0bzAYPK4BfnQqeqxbxVD+JsKViQRwyJAh+PjxIxYuXIiEhATY2dnh3Llz4okecXFxEuMfOnbsiKCgICxYsADz5s2DlZUVQkJC0LJlSwDA+/fvcepUwZY5dnZ2Em1dvnxZvCDxgQMHMHnyZHTv3h1sNhseHh7YtGmT+p8wIYQQooCzenGlriP4kr84+QOAmzwBbsbvBAD89/y4TmYMtBOI8nTwn1O2yZMnY/LkyTLPXblyRerYoEGDMGjQIJnlLS0t5ZoqXaNGDQQFBSkUJyHKsLS0xLRp0zB16lS5yi9atAghISGIioqSuw1nZ2fY2dlhw4YNygVJCFG562k3MLzUvXeqkZr1schz5/TkH+tGCuQjX9shFItukhONS0hIwJQpU9CwYUPweDxYWFigX79+EuM2LS0tZSYqK1euBIfDwZ9//qmRWL8dHF61alW0a9cOJ0+e1Ejbxfn999+lxrkSQnRPbm4Obt4/W+R4uhc8BitOjJQ41jqwpSZCI2qm61vkUQJISsXZ2RmBgYFyl3/9+jXs7e1x6dIl/Pnnn4iOjsa5c+fQtWtXTJo0qcTrd+/ejVmzZhW5C0xJbSszCysgIAAfPnxAeHg4fvjhBwwcOBDR0dEK16NKlStXlpq4RAjRPb8EdMT4qFmYv7foRYAvcyV73vK1MBuXFnuveCgBJBr166+/gsVi4c6dO/Dw8ECTJk3QokUL+Pj44NatW8Vee/XqVXz58gVLlixBWloabt68qZGYq1WrBjMzMzRp0gRLly5Ffn4+Ll++LD5/7tw5dOrUCdWqVYOxsTH69u2LFy9eiM9369ZNavjCx48fweVyi+zFS01NxZgxY2BiYoKqVauiW7duuH//vvj8okWLJMav5ufn47fffhPHMHv2bIwcORLu7u4S9YpEIsyaNQs1atSAmZkZFi1apPwLQwgp0X1+wfZtp9kvSyipXYdF93Aw+bi2w1AdFc8CVkV6/ONOG8za3VcFNakGJYBEY1JSUnDu3DlMmjQJhoaGUucLt+Eryq5du+Dp6Ql9fX14enpi165daopUtvz8fHGbXO7XJfgzMzPh4+OD8PBwhIaGgs1mY8CAAeIdPcaMGYOgoCCJ3TT279+POnXqoFu3bjLbGjx4MJKSknD27FlERESgTZs26N69O1JSUmSWX716NQ4cOICAgADcuHEDaWlpMrc+3LNnDwwNDXH79m2sWbMGS5YswYULF5R9SQgh5UiyHqUERWFYRaWA8o/ffMktWOYn8tFV1QRVSmVmEkhFMOT0ECR/SdZ4uzUr1cThvofV3s7z58/BMAyaNWum8LVpaWkIDg4W7//8888/o3Pnzti4cSMqV66s6lAleHp6gsPh4MuXLxCJRLC0tMTgwYPF5z08PCTK7969GyYmJnj06BFatmyJn376CZMnT8bJkyfF1wUGBmLUqFEyb0n/999/uHPnDpKSksQLi69duxYhISEIDg7GuHHjpK7ZvHkz5s6dK96j2s/PD2fOnJEq16pVK/j6+gIArKys4Ofnh9DQUPToId9epIQQ5T19E4Um9e20HQZRwrezo0vrQ8orAE4qq09ZlADqkOQvyUjKSiq5oBatWLECK1asED/+8uULbt26JXGL89GjR6hXr57UtaUZY3Lw4EE0atQItra2AAqW76lfvz4OHz4Mb++it+Jp0aIF3rx5I9H+twlj586dcfbs2WLb/uuvv+Di4oKXL19i+vTp2LRpE2rUqCE+/+zZMyxcuBC3b99GcnKyuOcvLi4OLVu2BJ/Px4gRI7B7924MHjwYkZGRiImJES9F9L0HDx4gIyNDaozfly9fJG4tFxIIBEhMTISDw9f9RDkcDuzt7aX2FW7VqpXE49q1a0vteEMIUQ+f8yNweqx2xw8TUogSQB1Ss1JNnW93woQJEr1fw4cPh4eHB3766esA56IWyraysgKLxcKTJ08UjnHXrl14+PAh9PS+vmVFIhF2795dbAJ45swZ5OUV/OX2/v17ODs7SyydUqlSpRLbNjMzQ+PGjdG4cWMEBASgd+/eePToEWrVqgUA6NevH+rXr4+dO3fC3NwcIpEILVu2RG5urriOMWPGwM7ODu/evUNAQAC6deuG+vXry2wvMzMTtWvXlrm8UUm3yUuiry+5xyiLxZJKEgkh6vGGW3IZoj5e29th7/i7Kq0z6uk1ldanSZQA6hBN3IYtrRo1akj0flWqVAm1atVC48aN5brW1dUV/v7++O2336TGAaampspMcKKjoxEeHo4rV65ItJ2SkgJnZ2c8efKkyNvK3yZZhcmjPLEWxcHBAfb29li+fDk2btyIT58+ITY2Fjt37kTnzp0BFNzC/Z6NjQ3atm2LnTt3IigoCH5+fkW20bp1ayQkJEBPTw+WlpYlxmRkZARTU1PcvXsXXbp0AQAIhUJERkZKLXROCCFlzd4zK5Cbny3z3DMeg/NhB9HT0bPEeu7xs3E14iSqVzXFosvjAF7pY3uXEqvwNbryRzeN+CQa5e/vD6FQCAcHBxw7dgzPnj3D48ePsWnTpiL3dt61axccHBzQpUsXtGzZUvxfly5d0K5dO41PBpk2bRq2b9+O9+/fo3r16jA2NsaOHTvw/PlzXLp0CT4+PjKvGzNmDFatWgWGYcRj9WTp3r07HB0d4e7ujvPnz+P169e4efMm5s+fj/DwcJnXTJkyBStXrsTJkycRGxuLqVOn4vPnzzq5uTshhMgrONQff348iI2fTxRZZsbTFfgskFxKh13Ed1/Gl89YeXkinvFo2RtKAIlGNWzYEJGRkejatStmzJiBli1bokePHggNDcXWrVulyufm5mL//v1SEy0KeXh4YO/eveLbvJrg5uaGBg0aYPny5WCz2Th06BAiIiLQsmVLTJ8+vchFqj09PaGnpwdPT0/w+fwi62exWPjnn3/QpUsXjB49Gk2aNMHQoUPx5s0b8faH35s9ezY8PT3h5eUFR0dHVK5cGa6ursW2Qwghuu74c/n+wFfkVuwjruZ+X+gyugVMSkXWOLWS1K5dG35+fsXeBn39+rX45+TkomdGz5o1C7NmzZKrXXm3APyWrPIsFguPHz8WP3ZxccGjR49KvC45ORnZ2dkyxywWPt/C66pUqYJNmzYVuff0okWLJNbw09PTw+bNm7F582YABbcYrK2tJcZryvq3krVUDCGEEPXRlf2LKQEkRM3y8vLw6dMnLFiwAB06dECbNm1U3sabN29w/vx5ODk5IScnB35+fnj16hWGDRum8rYIIUTXhUT4w6ax7GFFpAAlgISo2Y0bN9C1a1c0adIEwcHBammDzWYjMDAQv//+OxiGQcuWLXHx4kVYW1urpT1CCNEERslhzGc5b6B/5BepgW4slY98043ePGVQAkiImjk7O6t9n00LCwvcuHFDrW0QQkpvx8n52g6hTMmDUOlrT7GfKVT+asRJpdsqi2gSCCGEEI0K/GcpftxpgzUHpXe1KWuyc7IUKr85VfYC8EQ2kQZ72CbHLNBIO+ruEJAXJYBapitvBEK+R+9Noi7rko/gJRfYlxsmdW5+wAAM29EGYQ/OaSEy+eXn52HojtZwOdAO/7u2W9vhlFvJevL1AMr7faX6pbHK7lJblABqSeGODFlZiv31SIimFO5kwuFwtBwJqSgu3QnGKfZzRPPysPjO7+LjL98+xNGLfsjP153lOwLPLMNDXj4EHDYWv1iv7XDKrc8c1aYpx2K2QETrowKgMYBaw+FwUK1aNfE+rAYGBrRoLwHDMBAKheBwOFp9P4hEInz8+BEGBgYS2+8Rok7PP9wT//xen4Xd/1uCix9CEM0rSPxO7QrEthGXYWhQRVshiqVmJYp/zmHL/qyKhMqPXyPqcZufodL6Ur4kKN4JyOjGTiD0za5FZmZmACBOAgkBCpIvNlv7nfNsNhv16tWjP0yI1vyVclRiu64ofg5m7u+LLeOuai8oOVwJP4HTUTtgxDPWdihEzW7lPwH0Ffu+/pKbqaZoFEMJoBaxWCzUrl0btWrV0uhOFkR3MQyD9PR0VKlSReuJF5fL1YlElJBvXeelaDuEYt2KPo8pDxcC+gBE77QdDtFBF14dxlBM13YYlADqAg6HQ+OsCICCBDAnJwd8Pl/rCSAhpXHz/llEvbiMEa7zUMWwmrbD0Zi5d6YDevSHEylaFitX2yEAoEkghBBCVOzCrcMYHzULW9PPwv3QD3JfVx5mnidT8qcVIi29d4RK/Z2uG+9zeqcSQghRqc3Ry8Q/J5XThIhVhpf/KI9i3lzXSrtlOeEvu5ETQgjRSfnf9XBM2NEJw7a3wdM3UTLLf0pNwPaQeRBkfdRAdKQ82vXlsrZDKHNoDCAhhBC1usETAACWnhuDfePDpc7PPjxA5ctzaMPLtw+1HQIhcqMeQEIIIcW6cOswftrRCjN39SlVPVH8HKljIqGwXCR/2TlZ+OX8YG2HQcqALyzdWB+SEkBCCCEyJSS/xYOnN+H7aAme8Ric04vDf/dOl3gdU8zwuL1nVkg8/vPw+NKGqRP+DduPT2V4PBjRnFS2biSAdAuYEEKIlI+f4zHwlBvS2Cww32zH9e7j0xKvfadfdAb458eDEo/3591WPkgt+n6ZpvIwg5lohq78oaAbURBCiJa9/fCMtu76hv8pHwg4bDAy1qNMSH4L/+Mz8fTNAy1EppveJj/RdgiEKIQSQEJIhbf+8GT0/XcAvP520HYoOiNXJD1er9DcYwOxLf0cpp0fpsGIZPv4OR4fP8eXqo5r907Ad+/gImcpyyaZGO/IuliqGAjRNEoACSEVSmZWOk5c3oas7K/7cQZkX4WIxcJ9fi5inpfNW5Kaki/MQzg/CwDwlqudtfA+pSYAAB48vYm+IT3QN6QHHjy9qVRd+fl5mB+3HifwBHP/HanKMAnRaZQAEkIqlCn7emBhnD9+2+Mi83xubjZevX+C3Nyie8DKo+l/98SgHXa4ef9sseU2JOzTUERFcz7ZAyKhEGuvTEUWm40sNht/Xv5NqboyvqSJf37KE6kqREJ0HiWAhJAK5S6/oOevqKVHzkT8jQEXBmLgnrbIz8/TZGhaE3J5Oy7qf8ATnhBL7s4stmwOW7LX78Ktw7j7MFSd4cl0K+Y8MvF1T9UkTraGI6CdQEjZRgkgIYR847AoCkIWC6+4wJHQDdoOp1iZWenw3t4BI7a3xbuk10rXs+z1ZvHP74uZwSuLT+wyjL07FQ+e3VK6fWVkZKXio16++HG8gnETUtFRAkgIqbAinxS/f2h2bpb4Z5FQiDvRF3VqpvCao2Nwh5+JKH4OJpzsp3Q93/fqAcDb3LdyXy9ksbDl8u9Kt6+MGU9X4DNH87/CPgs+4s+D4/FCEKPxtglRJUoACSEVxvfJ28SwCXKP9fPZ7QrvyOn4dZeTOkJTyqvsl+Kf33BVV29wqL/MXTuKU55WwdseMg+Ddthh9/+WSJ1bdnQE9ubexDXeJy1ERojq0ELQhJAKY+3hCRKPs9hsvP/4qsjyIubrpIBQbiKAr/vaatPzuBjceXwWOaz8kgsrYfG7bWqpVxexZaxz6Cf4H8ADnqQcxS9YKD6emZWO8/rvNRkeIWpTJnoA/f39YWlpCT6fj/bt2+POnTvFlj969CiaNWsGPp8PGxsbnDlzRuL88ePH0bNnTxgbG4PFYiEqKkqqjhcvXmDAgAEwMTFB1apVMXjwYCQmJqryaRFCNGxfnmbHqanD9L97YsBlT6xM2ItHPPUkgMrIYL7g37AgbYehNiKhEMMOdNR2GISojM4ngIcPH4aPjw98fX0RGRkJW1tbuLq6IikpSWb5mzdvwtPTE97e3rh37x7c3d3h7u6OmJiv4zUyMzPRqVMnrF69WmYdmZmZ6NmzJ1gsFi5duoQbN24gNzcX/fr1g0hEywQQQrTnov6HIs+FP7yiUF25uTnw3uFYyogKPODn4venK7Hl+CyV1KdrrkedxksV3mYnRNt0PgFcv349xo4di9GjR6N58+bYtm0bDAwMsHv3bpnlN27cCDc3N8ycORPW1tZYunQp2rRpAz8/P3GZESNGYOHChXBxkb0O2I0bN/D69WsEBgbCxsYGNjY22LNnD8LDw3Hp0iW1PE9CiC4qGNmWnpmq3TDkNDp8ikLl/U7MwB2e7OVwlLU1vfh1BMsqvwhfbYdAiErpdAKYm5uLiIgIiUSNzWbDxcUFYWFhMq8JCwuTSuxcXV2LLC9LTk4OWCwWeDye+BifzwebzcZ///2n4LMghOgyRlT0rN6Nn09g0A47eB9w1lxApaTItmjv05+rMZLy4+b9s3jC053Z34Sogk5PAklOToZQKISpqanEcVNTUzx5Invj7YSEBJnlExIS5G63Q4cOMDQ0xOzZs7FixQowDIM5c+ZAKBTiw4eib7/k5OQgJ+frzLm0tIIV5hmGAcOUpzlyRF0K3yv0ftEcIVP8sA5Zv/h1+d9n9XFv/PnLmZILlkPK/Lt8f42sxwmfXpcmLEKkqOs7RJF6dToB1BYTExMcPXoUEydOxKZNm8Bms+Hp6Yk2bdqAzS6603TlypVYvHix1HGBQKDTvzCI7mAYBhkZBbfkWDJmJxLV+/PseICv2DUCgfZnAhflfX6i3PGJSkh+y5ole4fDQN8IY3uvlPuatAzJW+Dfv3YCgQBvPz5VSXyEFFLXd0hhx5M8dDoBrFmzJjgcjtTs28TERJiZmcm8xszMTKHyRenZsydevHiB5ORk6OnpoVq1ajAzM0PDhg2LvGbu3Lnw8fERP05LS4OFhQWMjIxQtWpVhdonFVPhHwpGRkaUAGpIGD9d4WuMjIzUEIlqpHJy5Y6PzdLpUUAKC0Y0kAdY3AnAkB7T5LuILTmT+vvXzsjICP+kXwBopxGiQur6DlHk94ZOJ4BcLhf29vYIDQ2Fu7s7AEAkEiE0NBSTJ0+WeY2joyNCQ0Mxbdo08bELFy7A0VG5mW41a9YEAFy6dAlJSUno379/kWV5PJ7EuMFCLBaLfpkTuRW+X+g9o7t0+d/mnT69d668PoahrOlylf3+tZL1OI/FgPb+Jaqkrs9ouUkAAcDHxwcjR45E27Zt4eDggA0bNiAzMxOjR48GAHh5eaFOnTpYubKgy3/q1KlwcnLCunXr0KdPHxw6dAjh4eHYsWOHuM6UlBTExcUhPr5gsHRsbCyAgt7Dwp7CgIAAWFtbw8TEBGFhYZg6dSqmT5+Opk2bavLpE0II0aK7D0ORrFe+ekoJAcpAAjhkyBB8/PgRCxcuREJCAuzs7HDu3DnxRI+4uDiJcXkdO3ZEUFAQFixYgHnz5sHKygohISFo2bKluMypU6fECSQADB06FADg6+uLRYsWAShICufOnYuUlBRYWlpi/vz5mD5dvr8oCSGElA2sEm6D/xI+TTOBEKJhLIZmJ6hNWloajIyMIBAIaAwgkQvDMBAIBDQGUMU+pSZg8pHeiOHlqaS+6JHRKqlHGTZ7bEosI298M/52K5dbm/2QY4Rt4+Rbsist4zN+ONZF/HggqyWCmZhiriCk9NT1HaJI3kH92oSQcm/VsV9UlvyVBSKhfGvWpQg/qzkS7ShNrwYlf6SioASQEFLuvRMWvX6nMibv7IrU9GSV1qlKx69sk6tcOD9LzZFoRwoykJubg6k7XTB+RyckJL/VdkiE6BxKAAkhREFXuclYHeyt7TCKtP/ldm2HoFVPeEJsPPYbLnETcZMnwMoTo7QdEiE6hxJAQghRwoP8l9oOoUgvuDS0Ozb9gfjnS9wk/BE4CLGv7mkxIkJ0CyWAhJByTx3pEKOWWom6hLCeYN6FUdoOgxCdQQkgIYQo4S2XZmmXNU95Ikze2RWHzq/XdiiEaB0lgIQQoiSbPTYYvt0eMc9vazsU8r0iOmivcpOx/EOA+HHmF8W3AiSkPKAEkBBCSuEBPxfTrvyi7TCIEnJzczD8VC9th0GIVlACSAghpZSor3tfpSGXi58JnJubo6FIdNfec8vxkbZ5IxUUvfMJIaQc+iPOT+bxhOS3sNljA/uDbTUckWals7KLPf/fvdPY+PmEhqIhRPdQAkgIIWWEvDt8FGfusYEqiET3PeLlF3veJ2q2hiIhRDdRAkgIIWXE//4LKLlQCcrr7h+K+sKmX3+kYqNPACGk3CsvK/YJMpNKdX1WdqaKIiGElHWUABJCyr2SbgdWBCKhED2DHLQdBiFER1ACSAjRirsPQzFrd1+cDzuo7VBU4vLdY9oOoVj/3gqCgENf+YSQAvRtQAjRipm3fsNZzhvMeLpC26GoxG+PFiH+4xtthyFTVnYmNj9are0wCCE6hBJAQohWfCqH66+tChmt7RBkmrWvL21dRwiRUP6+gQkhREs+MqnaDkHCx8/xAAq2PyOEkG9RAkgIIeXUimOjtB0CIURHUQJICCHl1EX9D9oOgRCioygBJIQQQgipYCgBJIRo3frDv2LQDjscubBJ26HoNIYpL0taE0K0jRJAQojGxb66J/E4IPs6nvCEWBq/U+W7VVy8fVSl9RUnD6Xfq5cQQjSBEkBCiJTYV/ew/+xqZOeoZ9/Ys3d3F3muZ5ADzvy3V2Vt/f54scrqKkksT4SNR6eqrX4WS/GlXOYH/KSGSAghZR0lgIQQCVnZmRh36WesTtqP5QdHarx9AYeNpU9Vt2ixUImkqTT+zrqE65GnNNpmcU6xn2k7BEKIDlI4Afzy5Qvev38vdfzhw4cqCYgQol0X7xxCyv8v0hzCeqKVGDLK+JZli+7N1XYIhBBSLIW+ZYODg2FlZYU+ffqgVatWuH37tvjciBEjVB4cIaT8EQmFuJ36n7bDUKukcrjLCSGkfFHoW2rZsmWIiIhAVFQUAgIC4O3tjaCgIAA0O42Qsmj/2dVYvHcoEj+9R+yre9hyfBaysgVqbTPo/FrE8PLV2kahzKx0jbRDCCFljZ4ihfPy8mBqagoAsLe3x7Vr1zBgwAA8f/5cqcHJhBDtefDsFlYn7QcAfDz+M56wkpCoz4Z5CgPoq+fz/Ck1AXvf71Nb/d9bc3SMRtohhJCyRqEewFq1auHBgwfixzVq1MCFCxfw+PFjieOEEN13JeqQ+Oer3GQk6hd8HcSrKTn7+Dkezid74IOc9Z+7eQAiYemWVTmOR6W6vjS2npiLDUem4Fb0ea3FQAghRVEoAdy3bx9q1aolcYzL5eLgwYO4evWqSgMjhJQv44+6KVR+5rNVsN1vV+okUFu2pJ3Gri9XMDHCB1nZmXj6JkqjaxISQkhxFEoA69atCzMzM6njeXl5qFevHmJjY5GSkqKy4Agh6sOC/D1xqvCMp9w44UMXN6ikfW3JZ7FwO/ocRlwajulPluDvU77aDokQQpRfBzA9PR1bt26Fk5MTqlatCktLS1hbW8PExAT169fH2LFjcffuXVXGSghRJTnH7fo/XKXmQIr37pN2lqJRpRNR/shiF3zdbvx8XOl6aLIdIURVlEoA169fD0tLSwQEBMDFxQUhISGIiorC06dPERYWBl9fX+Tn56Nnz55wc3PDs2e0ECkhuobFku/jr4oRgRV9Nm4sK0nbIRBCiASFZgEXunv3Lq5du4YWLVrIPO/g4IBffvkF27ZtQ0BAAK5fvw4rK6tSBUoIUS1NztxffMgT4GisOZ2jqok1tNoCIURVlEoADx48KFc5Ho+HCRMmKNMEIUTNHn8KB7gll1NFynGW80b5iynpIYQQlaPl6gmpoK5ykzXSzv6zpdvXNzlLeutJQgghpaOWBNDT0xNpaWkAgFOnTuHoUVr6gJCy6iUXOHllh9LXP/xwo1Ttn9WLK9X15QlNAiGEqIpaEsCHDx+iatWqePToEebNm4crV65gypQpStfn7+8PS0tL8Pl8tG/fHnfu3Cm2/NGjR9GsWTPw+XzY2NjgzJkzEuePHz+Onj17wtjYGCwWC1FRUVJ1JCQkYMSIETAzM4OhoSHatGmDY8eOKf0cCCnLFrzZrPR6fJSyEEKI7lFLAqivrw+GYRAQEIC5c+fC398fN2/eVKquw4cPw8fHB76+voiMjIStrS1cXV2RlCR7Vt3Nmzfh6ekJb29v3Lt3D+7u7nB3d0dMTIy4TGZmJjp16oTVq4u+NeXl5YXY2FicOnUK0dHR+OmnnzB48GDcu3dPqedBCNGMk1f/hscOW8wPGKDtUIqUmq7c7XeaBEIIURW1JIATJ05EmzZtEBwcDHd3dwAFSZcy1q9fj7Fjx2L06NFo3rw5tm3bBgMDA+zevVtm+Y0bN8LNzQ0zZ86EtbU1li5dijZt2sDPz09cZsSIEVi4cCFcXFyKbPfmzZuYMmUKHBwc0LBhQyxYsADVqlVDRESEUs+DkIqK0XAf4ILXG/GUJ8Ip9nONtquIJYeHazsEQkgFp9Qs4JKMGTMGAwcOhJ6eHgwNDfH8+XN06NBB4Xpyc3MRERGBuXPnio+x2Wy4uLggLCxM5jVhYWHw8fGROObq6oqQkBCF2u7YsSMOHz6MPn36oFq1ajhy5Aiys7Ph7Oxc5DU5OTnIyckRPy4cB8kwDI3dIXIpfK/o4vtFKBKCxVb8b0ZVPBddfD1K44J+vMLPKfHTe2TlVOz1FAkpL9T1naZIvaVKAN+8eYPY2Fi0atVKaou4atWqIT4+HpUrV0bjxo0RGBiocP3JyckQCoUwNTWVOG5qaoonT2TvDpCQkCCzfEJCgkJtHzlyBEOGDIGxsTH09PRgYGCAEydOoHHjxkVes3LlSixevFjquEAgKHe/wIh6MAyDjIwMALp3u0+QJoAeR1/h61Tx3hcIBKWuQ9co8pzO3dqL5Yk71RgNIUST1PWdVtjxJA+lE8CDBw/Cy8sLQqEQfD4f27dvx4gRIxAXF4egoCAcP34ckZGRyM/PV7YJrfrjjz+QmpqKixcvombNmggJCcHgwYNx/fp12NjYyLxm7ty5Er2PaWlpsLCwgJGREapWraqp0EkZVpgsGRkZ6VwCaFTVCHp6iieAqngaRkZGpa9Ex8w60gd8Fh9/ev0DLpdXbFlK/ggpX9T1nabI7w2lE8ClS5diypQp8Pb2xrx58zBx4kQ8ffoUq1atQqNGjdC9e3eJW7fKqFmzJjgcDhITEyWOJyYmSvU4FjIzM1OovCwvXryAn58fYmJixLud2Nra4vr16/D398e2bdtkXsfj8cDjSX+Rs1gsnftlTnRX4ftF194zDMMoFZMq+r7P3zoIV8dhKqhJd4Tx0gCkYdPxqZjpuV3b4RBCNEhd3++K1Kv0JJAXL15g6tSpaNGiBfz9/ZGVlYUbN27gwYMHePLkCfz9/TFgQOlm4XG5XNjb2yM0NFR8TCQSITQ0FI6OjjKvcXR0lCgPABcuXCiyvCxZWVkACsYbfovD4UAkEsldDyFENX5/ulLpZWh03dP0aG2HQAipgJROAPPy8lCpUiUAQN26dcHn87F27VpYW1urLDgA8PHxwc6dO7Fnzx48fvwYEydORGZmJkaPHg2gYLmWb3sap06dinPnzmHdunV48uQJFi1ahPDwcEyePFlcJiUlBVFRUXj06BEAIDY2FlFRUeJxgs2aNUPjxo0xfvx43LlzBy9evMC6detw4cIF8axmQioaEaNsAqaa8a+L9g9VST26hsYHE0K0oVTLwAQFBYknY3A4HFSvXl0lQX1ryJAhWLt2LRYuXAg7OztERUXh3Llz4okecXFx+PDhg7h8x44dERQUhB07dsDW1hbBwcEICQlBy5YtxWVOnTqF1q1bo0+fPgCAoUOHonXr1uJbu/r6+jhz5gxMTEzQr18/tGrVCnv37sWePXvQu3dvlT9HQkjJTkD2xK9vbT9RumEn2nCbn6HtEAghFRCLUfLPTycnJ0RFRSEjIwPVq1eHQCDApEmT0LFjR7Rs2RJNmjSBnp5aVpkpM9LS0mBkZASBQECTQIhcGIaBQCDQyCQQmz2yJzMV5e7Q2+DzDBRup9PuFhBwVLPkaPTIom+Xng87iBlPV6ikHU0r7nkBiv9bEUJ0W0mfeWUpknconaFdvXoVAPDs2TNEREQgMjISkZGR2Lt3L1JTU8HlctGkSRM8ePBA2SYIITpEpOStSlUlfyW59PggwNFIU4QQUuaVuovOysoKVlZWGDr06/icV69eITw8nLZNI4QQQgjRQUr/ab5w4cIit0Vr0KABBg0ahBUryubtGEKINIZmwKvNyat/q6QeNk0oIYTISekE8N27d+jVqxfq1q2LiRMn4uzZs8jNzVVlbISQMs7v2O/aDqFMWPB6I/Lz80pdz4YW0jsREUKILEongLt370ZCQgIOHjyIKlWqYNq0aahZsyY8PDywd+9epKSkqDJOQoiW/XXsVxy96KfQNdsz/lVTNOXPnYcXkZ2TJXVc1rGidG3ngclV+6IB/S1OCClBqUZns9lsdO7cGWvWrEFsbCxu376N9u3bY/v27ahduza6dOmCtWvX4v3796qKlxCiJYdFUVjyfjtint/Wdigy6da+KYobHzUL7Q61R2ZWusTxKXtc5LreIadyQT0DVuLU2Ggccz6g8hiJZq1rMg+u+XW1HQZRA11Y2F6l0/Osra0xa9Ys3LhxA+/evcOoUaNw/fp1HDx4UJXNEEK06HzEPrnK3X0YWnIhFdl8dDpOc15prD112v6/2RKPb/HSiyj51fjKrtj+yzWJY03qt1JpXETzejp6YkqvjdoOg5RTalmoLzIyEuvWrcPnz59hY2MDDw8PdTRDCNFh12OOaaSdhy/CsSProkba0oSsXMUXhrZr1A16evpqiIZoG1vH9gQnqiFiRGBred0qtSzQNWTIEPTp0wfLly9HkyZN8NNPP+H8+fPqaIoQomGRqWFylVPHFmeyJpXEvPxP5e1oU+EC4KqYFDLWoHup6yDaZVzdXNshEDUQ6cCqCmpJAKtXr45hw4ahdevW8Pb2xsWLFzFz5kx1NEUI0bD7fPlmGDAq2gP4W9sz/oXNHhscOr9efIylnq8xLWLBe3sHdN1rhyMXNsl1RVG9RJN/WqfKwIgWGPANschiIrrlmmo7FFLOqOWbs1GjRli3bp24B6BatWrqaIYQoiXHLm0psYw6egALLf8QIP6ZxSpfCeDp3Lu4w89EKoeNpfE7S1UXm0Nbo5RV19wviX/26PYrNo4tP8McCCBiytkkkEI5OTnYunUr6tWrBzc3N7Rs2RIuLi40G5iQcmLR260lllFHD6AszxPDNdKOpmQosXWecVUzNURCtKm6kYnUsRr52r9tSFRD2a01VUktk0COHz8OAMjMzER0dDQePHiABw8ewNPTE/Hx8Xj+/Lk6miWE6BINfcEF5cvekagiadqgtbZDIApolc3Fg2+GUvBFDLLZJU/2mNd8PhY8WSFXWUJKopYEsJChoSE6dOiADh06qLMZQogO0lQPICmeVQ4Lz3j0b6FLdo26jnaH2gMAfhQ1wTzPvXj1/jGG/jcaADDOsIfM61wdh6Fjq97oGNxZY7ES9dCFrTVLlQC+efMGsbGxaNWqFczMpG9BxMfHw9ycZjARUhGl5nxU0yAToogaMASg+NIy2uKWb4Fzem+1HYbajK/sCj7PALcG3cSjV3fQrkXBTO0Wjdpi8+cleJ/8DENcphd5fRXDavAxHoz1n45oKmSiBmV6DODBgwfRuHFjuLm5oWHDhti3r2Bx2Li4OKxatQrt27dHvXr1VBYoIaRseZ8br+0QCABznoW2Q5BbX2EDLB5+VNthqNVkj7UAAEODKuLkr5Bz2wEY7jarxDUdR/aaB3emmdpiJBWD0gng0qVLMWXKFERHR6NHjx6YOHEi/vjjDzRq1AiBgYFo27Ytjh4t3x9kQkjR7vGz1Vp/bm6OWusvL37tU3aWglk+8gQM+IYyzzXL4cCS9jgGUDC7e+ko+v1alokY7d8CVjoBfPHiBaZOnYoWLVrA398fWVlZuHHjBh48eIAnT57A398fAwYMUGWshBAA2TlZmLCjM8bt6ITPgo9YHTQGLn+3wPrDk+Wu48Ktw6WOY8fJ+aWuozRm7+2n1fbLCrOaZacHsHDZmuUNpkmdW9YjAPVRU8MRlZ5NztfevC65NVRat14xE63m1/5FpW0R1dJja3/nHqUTwLy8PFSqVAkAULduXfD5fKxduxbW1tYqC44QIm323v64wUtFGE+AoUe7Yn/ebSTqsxGQfVXuOu69LP2aYptTT2Hj0amlrkdZF/U/aK1tol79u3hLPG6Ww0HTBq3R2MhGSxEpxzhfhOV9gnB7yC3sbLMO/mPl/4zK49v0zym3JkbxOqFnnjm2263B0J7TscVmuUT58ZVd0SSHjSpCEX5Cc5XGQhRjaFBF2yGUbhJIUFAQ3Nzc0KxZM3A4HFSvXl1VcRFCinCJmyj+OV5fu8tB/J11CdpLAXXLAnNvMAyD5R92azuUMqm/qHGR5/5yDwEA/ObxF3btt9NMQCoQMvgqqlUp6LXsYNNTrW2xwcKMoZLrc3Zu0x8/hP+J29zPGM7vhMkea/HtfYIqhyZgT84NtcZFdJfSPYCdO3eGr68vWrRogZo1ayI7OxsbN27EkSNH8OjRI+Tn56syTkKIHGbt7ovsnCxth1EhDXCaiKE9i569qU0/5BhpO4QSLfUKLvJc3VqWAApuEe9q85eGIvrqD/Ox8K7kDLd8C9STYxyioUiE/zy+Jn/q4lXJWfzz4LY+MstsG3cdlwZcxu9Dt0md+33oNkyp5o5+osY45nxAXWESHaV0D+DVqwVd2c+ePUNERAQiIyMRGRmJvXv3IjU1FVwuF02aNMGDBw9UFiwhpHhnOW9QI3gy5gyvGL1Q4Q+vaDsEMS6XBwAYptcWQflfdyeZXmMQ/krR7oD9Be770eus7o6Z/CGnmtS2dSN5P+B05nW48ttJHHewcUG0TTQAwGaP+m4JR4+MRn5+ntSMXJFQCNsieiFts3noVX8w2lm7waiyasf7yTJ5wDoYnlmMmlXrolPrvkWWk7WrSKFxPy4V//ztAtWmeSIk6tM6TuVZqReCtrKygpWVFYYOHSo+9urVK4SHh+PevXulrZ4QoqDIzEhth6Axo8OnaDsEKXOHByDom8Ski+1P+OuydhPAwh40XbXF+4rUsd+HbsPvGo+kgEVuweg6WcuxsDkc3P85Cs/iHqBpg9a4ePsopj9ZguY5elj54xFYmDXUWJxcLg/j3VeorL4GPEs8wFMAQBNRTSQiRWV1E92jVAIYFxdX7Bp/DRo0QIMGDTBo0CAAwPv371GnTh3lIiSEqJwqd+nIzskCn2egsvrKokOdAoo8Z2pcV4ORlD0juI5SvX/yap3NV8tyQ7sHnC32PJvDEW+/59J+EKLbD1J5DNowe9AufNjvhjwI8XvvLbh+aWjJF5EyS6n+3Xbt2mH8+PG4e/dukWUEAgF27tyJli1b4tixY0oHSAhRPUaFa1CN2NNRZXWVRZWFIrRo1FbiWJ28ggTbMrdg54ZqQu2v+aWrfh+8teRCReCoYauZrbary9TSOapUxbAado2/hb3j78LMxFLiXIccyVmrvYWS50nZo1QP4KNHj7B8+XL06NEDfD4f9vb2MDc3B5/Px+fPn/Ho0SM8fPgQbdq0wZo1a9C7d29Vx00IKUIOq+QthiIzIwGeatp7whPieVwMGtdrqZoKy5ipdb2ljp3yisCdhxfQsVUvAMDVkQ+KHDemKQ1ygVdcrYYgk7K9fwCgx1L+2qJ0sqPfVwDA1+ejUS4LL7gMXPJqw6efPzb+8xuiEIeWTG3MGbwbZ0K6aTtMUgpK/flkbGyM9evX48OHD/Dz84OVlRWSk5Px7NkzAMDw4cMRERGBsLAwSv4I+Y5IqN49IF9ygfz8vGLLvNJX7Sx9EVMxZ/03zmFhaE/p2ZdcLg+dWvcVJzdsDge3h9zSdHgS1vbUnVmefi2XwSWvNtY2mVuqeqZ0XS/xeDDbFvd/jkIV6nEtNTaHgy3u/8P82qOx/OcTsKhthbVjzuLimIfYMPYiqhuZoJewvkJ1/qzfXi2x7m67ASP0O6il7vKsVP3nlSpVwsCBA7FhwwacOHEC586dw/79+zFjxgy0bFkxewMIKc7MXX3gHNgKO07+AQD4L+oMpv/dE//8F6jSdi7cPlTsedWNAPz/+kSSNf5377SKW9CcqnImD/Nrj8aJcfKvcmDAN8TZXv/DknqTcG94JAaxWykbolKa1NdMe4W3v4vCFzFwsv8Rf405D1fHYaVqq5VVB6xtMhc18kWwy+ZhxqDtYHM4+G/kA/ys316uJVtI0cxN6mNoT58it+eb6xGAH5kmGM3vggOO23HN/VKx9U0f6A/O/+9e0kfYQOq8WR4D8xLeP1Ix5jFo16I7Zg3biW65pgpdW9HRHG9CNCQ9MxXn9OLwWY+NzakhAIDf783ERf0PmPNCvv1aT1yWXstLltMPdxW7V66qE8DvrQovXc+ONpnn66O7HL9IZPX8laRuLUsM6DoBenr6mOcZqER0uu/0qHsYzLaFfXYlbG21Uup8pWK2L1OGq+MwXPV+iH3jw8WJCpvDwexhf+OfsdEK1VXc1mpEWnUjEywbdQw+Q/zRqklHVDcywZHOgRgAa/zVbKHEuMG+wgbgcnmIGhWD6JHRWPXLKam9nUfWGYEN3faA9f//DlOrl7yd7K9W08Q//2DZXyXPq6KgBJAQDcnLl+6OyGQr9hF8lSRfj9M13if4nSh6UeJctmp3EPl+Y/P3+mX3FykLwPpf/sVGa19MrzFQfPwP87GYUMUNVYQiDNNrW3QFcpK1xEhZ9yPTBHp6+vhjxH4Ejr9T7Np0muLXcpncZcvuu1Z3WDe0x5KRR+DSfhD8R13GHNMR+L3mUPwx7KBU2Xnt/0T9//9a/LVqXwzr+TusG9pji90a+NadgF/6+BbZjjvTDJON+uFHpzHiY4N7/Cb+2SHbEI45VVX3xMqhUq8DSAjRHEaBHoqA7OtQvI+K1OPUBpvDQTeHgRAJBwBnWNDj6Il/uUzCn1qOUHmWucBrNU0EWdVoBvp0GiV1fLheOxzI/7pihFfNgVJl1MnJ/kdE2/8IQL0LRxNpXC4Pw91mFXnesZUbTrdykzr+7UScQSwbXMm7jwFGrqhfqwX2PNkEG541Fo2SPczl/s9RCH90GW2bdwWbw8GBc2uwKnFf6Z9MOUQJICEawmJpd99edboRcxLWDe21HYZK/P7jDvHPbA4Hv/RbqLa2VjWaIfftf1VozrHEa7xWS921qsteOmXO8N2YlpOFC7cPIS0rGZ49ZqilfXkU7m5hn10JG4efw6OXd3D7yRns+nIZANBDWPT6tkQ7FnoF4dtPYP8u0rPuv8XmcOBg4yJ+7NljBnLOfIGIEaG7vSf6XywfazaqAt0CJoSU2sbPJ9Q+u1lTNLkGnKweM3X6fcB2tdVtaW5d5Dk+zwD9uvyC4W6zSrXsS2lt7rYX4wxcsMrjGIwq14BjKzdMHbQRQ1n26CNsgDkeFWMLxYqk4I84X4zpvxj1zawkzo3mO2kpKt1ACSAhmqKCHsDMvDQVBKIeswN1d69ZUsCkunmZrFtVrBvaY8qgv6SS/En9N2Dl6JMwrmampciIJrA5HPxpNQfWORxst1sDnyF+2g5JqygBJERD2CzJj1uKIEnhOk4I76sqHJU7p/dW/HN+Ob7dTQgpu9w6DseRcVHoaFuwSPuBjjvF5w7+8DeiR0Zju90abYWnUTQGkBBN+S4pcgrprnAVup5Yuf3dEsPMS7e2W0XTM88c5/XjNdZeixw9POSpduHuX6tqf7YvIcpoZdUB12oVrF9Y3cgEANDRthcG3t+DYOahNkNTO+oBJERHlIcxdO/1Wfjzo/RyD2VFuxwDjbc50/1vjbZnY2Cr8jonDpBe74+QsqK6kYk4+Svk63UIPzJNtBSRZpSJBNDf3x+Wlpbg8/lo37497ty5U2z5o0ePolmzZuDz+bCxscGZM2ckzh8/fhw9e/aEsbExWCwWoqKiJM6/fv0aLBZL5n9Hjx5V9dMjFURJs4D77rZF9PPbGoqGyMKF5jfL1eSkEwCYMWiLSuurTtuukXJq2ahjuNDnTMkFyyidTwAPHz4MHx8f+Pr6IjIyEra2tnB1dUVSkuzxUzdv3oSnpye8vb1x7949uLu7w93dHTExMeIymZmZ6NSpE1avXi2zDgsLC3z48EHiv8WLF6Ny5cro1auXWp4nIW+5LKy+NEnbYVRoun2DXTX4PAMMZKlmq87RfCf4ddpZckFCyiizmhbwrTsBPfIkJzmd6Xkc4wxciriqbND5BHD9+vUYO3YsRo8ejebNm2Pbtm0wMDDA7t2yp+tv3LgRbm5umDlzJqytrbF06VK0adMGfn5fZ/uMGDECCxcuhIuL7H88DocDMzMzif9OnDiBwYMHo3Llymp5nqT8k2cdwDf62RqIhOgaTS9H4e1W+lu2Ywy6wWeIH1o16aiCiAjRXQO7T8L6Mf+iQ04VcEUMxhp0h0VtK7g6jJQo55RrrKUIlaPTCWBubi4iIiIkEjU2mw0XFxeEhYXJvCYsLEwqsXN1dS2yvDwiIiIQFRUFb+/iF6AkpDgslm583NIzU7Udgtr953EV93+Owu81h0rsR6qrJg/Q3GLQQMGexIpwZ5rh3vBI8QzJuWajMGnAWvUER4iO2jnuJi4NuobfBm34/yOS3+nrR/6r8ZhKQ6dnAScnJ0MoFMLUVHJjdlNTUzx58kTmNQkJCTLLJyQkKB3Hrl27YG1tjY4di/9LNycnBzk5OeLHaWkFa7YxDKPQFl6kfJLnPcCCat8vsurZdGKqSurWVbNMhqOqYXUAgFfvefDCPFyJOIHfHha9r+i3tPFZ1dcv3dhDZWI2y2OQoF9yr/QDr6/7TzMMA8dWbnD8/+27ysP3WuHnrTw8F6J+VQ2ri98r379n9PW5WFZ/KrY834D4Ej5b6nq/KVKvTieAuuDLly8ICgrCH3/8UWLZlStXYvHixVLHBQIBfbkQxLyWb4JHRkYGANVsHScQCCQei4RCHBJGlrpeXdav4wSp5926cTf0vh+IM3qvir22nZmr1LVlgTIxHx5wFU6nnaWOVxaKkMEp6Nn4hdetTL4eimAYRqWfOVJxGOhXE/9ski+CQCCAk91AONkNROeTnYu8bijLXm2fq8KOJ3nodAJYs2ZNcDgcJCYmShxPTEyEmZnsFdvNzMwUKl+S4OBgZGVlwcvLq8Syc+fOhY+Pj/hxWloaLCwsYGRkhKpVqyrVPik/Tt3bUuKgCxaAypUrw8jISOqXUVpGisJtGhkZSTzeearkP2TKKrM8Bl51Rkg950JLfj6CM4faFXn9z/rtMbrffHWFV6IOOVVwi5eu1LVFPWdlXBtxD1nZ6XgW9wD2zZ1VVq+uKvzjXNZnjpDiGBkZYVoND0R8uIwRHeZLfA4X1hmHJe93SF3TNIcNHy9/8HnqWXJKkfewTieAXC4X9vb2CA0Nhbu7OwBAJBIhNDQUkydPlnmNo6MjQkNDMW3aNPGxCxcuwNHRUakYdu3ahf79+8PExKTEsjweDzweT+p44RIyhJSExUBi2aFvrQz+ReFRuzm5X8RfNPEf32Bz6klVhapT/u19GuYm9Ystw+Pxsaz+FCx4s1nm+dnDNLse3/e2/HIVbQ60UepaZb9f/Fouw+SYBeLHg9itoK/PhZG+Mdq26KpUnWVRUZ85Qkri3W8RvLFI6vgglyn4ciYdt9+dxzvWJ7z8/1Eeq9wOoBLfUG3xlJsEEAB8fHwwcuRItG3bFg4ODtiwYQMyMzMxevRoAICXlxfq1KmDlSsLZrVNnToVTk5OWLduHfr06YNDhw4hPDwcO3Z8zcRTUlIQFxeH+PiC1fdjY2MBQDzjt9Dz589x7do1qXUECVFO6X65/MN+qfA12blZ4gRweYgXtLDMndqd7nGsxOSv0I/O4/AjxgEATl7ZUWQyqA36evrgixhksxV7n/QWWirdppP9j7hv1xd/HZ2MzFwBZnhsU7ouQoikwjHIbz88w/Z/56B57Q5oXE81SzCpgm5MSyzGkCFDsHbtWixcuBB2dnaIiorCuXPnxBM94uLi8OHDB3H5jh07IigoCDt27ICtrS2Cg4MREhKCli2/vuinTp1C69at0adPHwDA0KFD0bp1a2zbJvnlt3v3btStWxc9e/bUwDMlBGBBtWNF7z/9T/xzPPNZpXWr0w851TBMry0AwEgogl/LZWiUK5kYVRGKcP/nKNQ3V261/h+dx5U6TlWbYqr4NnpzBwWUqk02h4MZQ7dioVcQDA10f8Y0IWWNRW0rLBt1DMNcZ2o7FAkshmYnqE1aWhqMjIwgEAhoDCDBrIAfcbaEXjzjfCFCPG7KHI9ks8dGqXaPdA6EdUN7uO9shRdc3f+4/2k1B24dhyM/Pw/BlzbDplEXtGjUFlnZmZiw1xnx7Cx00GsG7x7L0aBOs1K19eNOG7zkFiSc28ZdV9EzKB1F/52jR0arKZKKgWEYCAQCGgNIygVF8g6dvwVMSHnxOe8TID1EVMInPY7K9wTeGjoTmxpeglDFvYvqYm7SCACgp6ePoT2/Tqoy4Bti77i7Km1rU59juBB+AAO7TlFpvaXxp9UczHy2StthEELKOZ2/BUxIRXPk6nqV1neZ+xEnr/2N12Vk/F/tmvU01lZ98yYY038xqlWpqbE2S+LiMFjbIRBCKgBKAAnRMQmCFyqvc8GrjSqvUx36ihrBpLp5yQXLMT09fWy0lm/RakIIURYlgIQQnbFydIi2Q9AJ3RwGajsEQkg5RwkgITqnbIzVU7WmOfR19K0zPY+XWGaO6QgNREIIKY/oG5cQohOmtVuu7RB0ikVtK2xusQQ22fpFlhnuNkuDERFCyhNKAAnREHmXmDjGfoRHLyPUHI3u6dS6r7ZD0DnObQdg/5i7mF/7F6xqNEPb4RBCyhFKAAlRo+TUD/Dd64ldpxcrdN2yUN1bpLjQGVfVbyfXM6+OyussL9gcDob2nI4+nUbBNvvrVG6rHFqzjhCiPEoACVGjZcEjcZyJwYZPwfiIdLmvi+HlqTEq5XlxO8LCrCF22K1Rab2mBhYqra+8WusRgqEcewxktcQatyBth0MIKcNoIWhC1ChU/+s2hc+5ql3gWZN+FDXBEq8jYHM4AABH215AlOrGn/H0DVRWV3lmVtMC838O1HYYhJBygHoACdFRO08v0nYIAIAOOVWwbPQxcfJXKLjLXnTKqVbq+mvnMRjbb0Wp6yGEECI/SgAJ0VGbPh3DwQt/aTsMzO8TIPN40watsXXcdUys0qtU9Z8YfgsGfMNS1UEIIUQxlAASosNWxO/WavsRnuGwrNO02DK//rQG93+OUqr+QexWMKxUWalrCSGEKI8SQEKITLXyReByeXKVZXM4uDv0tsJtzPMMVPgaQgghpUcJICFEil02DysdFNs/mM8zQPTIaAxmt5Kr/CmXo9DTK3qRY0IIIepDCSAhSjp9Yw88d7bDpuMz1drO/ac31Vq/LL91Wg0HGxelrv1jxAG0+ma9OlmGsG3RoE4zpeonhBBSepQAEqKkuc/XIoabjZ3p5xB4fg1s9thg7p6BKm9nxjXNLwptaW5dqusPjI/AMed9Ms+N5P6AecP2lKp+QgghpUPrABKiAus+FCQ7pxGLLmEH0MtxuMrqTtTX/I4PJtXNS11Hk/p2MM4X4ZNewd+Zl/tfQM3qZqWulxBCSOlRDyAhKvbozX8qrzM/X3M7gxzuLHvZF2WscvgLPfJqY1n9KZT8EUKIDqEeQELKgL3nVmqknaCOO9G8YVuV1dfBpic62PRUWX2EEEJUg3oACVGxbGG2yuv869NRldf5vUUWE2Fj1UHt7RBCCNE+SgAJUbFD+eHaDoEQQggpFiWAhPy/zJwMnL1/HHkizY230yVWde20HQIhhBANoTGAhPy/kXu7IpabDefIXVj18xEY6les/WlbNemo7RAIIYRoCPUAEvL/YrkFY/eusOPQIagDbPbYIC0rVam6zt85pMLI1E+VM38JIYToPkoACSnGkmPKLcK89sFyFUeiPv1FjVU685cQQojuowSQkGKE5T1U6rp8lkjFkajPnMHU+0cIIRUNJYCEqIHm9+5QXhXDatoOgRBCiIZRAkhIMRglr2Mre6GG1cwvOz2VhBBCVIcSQEKKkc6R/RG59+JWsddxoDuJ1WhuJ6ljhiIR7LK5WNp6lRYiIoQQom20DAwhJUjJ+oQaBsYSx/Ze8S320/NeXw9CYb6aIyvafLNROPomEO0qtYGP51Z0j72JwKvzUYNbC01qt4e780TweJW0Fh8hhBDtogSQkBI4HXXG3s470brh123S8kW5JV6XmPJenWEVqUa+CENdZ2AoZoiP2TbtiL+aXtZKPIQQQnQP3QImRA6zL0+UeCzPEL98keZ7AFtks7DIZonG2yWEEFK2UA9gGXb17VXcjL+J8bbjUYNfQ9vhlGsf9CSTOQGTWeI1Q878pNE/sSxyGRwaH625BgkhhJRZ1ANYRuWL8rEuYh2CngShz/E++Dv6b3zJ/6LtsCqMfJTcu5fB1uxEkB6VO2u0PUIIIWUXJYBlVOznWCRkJgAAMvIysDFyI/qe6IsTz05AKBJqObryL4aXp+0QJAxn2eHXn9ZrOwxCCCFlBCWAZVQL4xb4Z8A/GNhkINisgn/GpKwkLLy5EINOD8J/7/8Dw5SRxehIqc3x2kezegkhhMiNEsAyzMTABL6OvjjR/wScLZzFx599foaJFydi7IWxePTpkfYCJIQQQohOKhMJoL+/PywtLcHn89G+fXvcuXOn2PJHjx5Fs2bNwOfzYWNjgzNnzkicP378OHr27AljY2OwWCxERUXJrCcsLAzdunWDoaEhqlatii5duuDLF90bZ9ewWkNs7rYZu113o6VxS/Hx2x9uY8jpIZh7fS7iM+K1GCFRp63Nl2s7BEIIIWWMzieAhw8fho+PD3x9fREZGQlbW1u4uroiKSlJZvmbN2/C09MT3t7euHfvHtzd3eHu7o6YmBhxmczMTHTq1AmrV68ust2wsDC4ubmhZ8+euHPnDu7evYvJkyeDzdbdl6ydWTsE9QnCn05/om7luuLjp1+eRt8TfbEufB0EOQItRqi7yurt8mEsO3Rq11/bYRBCCCljWIyO/+Zr37492rVrBz8/PwCASCSChYUFpkyZgjlz5kiVHzJkCDIzM3H69GnxsQ4dOsDOzg7btm2TKPv69Ws0aNAA9+7dg52dncS5Dh06oEePHli6dKnSsaelpcHIyAgCgQBVq1ZVuh5l5ApzcST2CLY/2I7UnFTx8arcqhjXahyGNhsKHoen0Zh0GcMwaLW3VbFlokd+XWLFZo+NukOSy5hKXTF18CZth0FImcUwDAQCAYyMjMBisbQdDiGlokjeodPrAObm5iIiIgJz584VH2Oz2XBxcUFYWJjMa8LCwuDj4yNxzNXVFSEhIXK3m5SUhNu3b2P48OHo2LEjXrx4gWbNmmH58uXo1El6X9VCOTk5yMnJET9OS0sDUPAFo+k8W5+tj+HWw9GvUT8ExARg/+P9yBHmIC03DWvD1+LA4wP41e5X9G3QFxw2R6Ox6SJ5tm3Txb+VPLvP1sm4CCkrCr+f6XNEygNF3sc6nQAmJydDKBTC1NRU4ripqSmePHki85qEhASZ5RMSEuRu9+XLlwCARYsWYe3atbCzs8PevXvRvXt3xMTEwMrKSuZ1K1euxOLFi6WOCwQCrX65jGo0Cr3Me+Hvx3/j37f/ggGDD5kf8MeNPxDwIADjmo9DR9OOFfqv3zuJt0ssk5qaqjOv0Xyjn9HAog24HEMIBHRbnxBlMQyDjIwMANCZzzchyirseJKHTieA2iISFSzgO378eIwePRoA0Lp1a4SGhmL37t1YuXKlzOvmzp0r0fuYlpYGCwsLGBkZafwW8PeMjIywxmwNvFO8sfHeRvz3/j8AwMv0l5hzew5a12qNqa2noo1pG63GqQ0xSdGYcev3EsvxKnFRiWegtjisc3LxmMcFAPROA+pWbYsRP65B55PdJMr1SmdhiNcstcVBSEVS+Mc53QIm5YEi72GdTgBr1qwJDoeDxMREieOJiYkwMzOTeY2ZmZlC5WWpXbs2AKB58+YSx62trREXF1fkdTweDzye9Lg6FoulM18szYybYavLVtxNuIu/Iv5CdHLBuLZ7Sfcw6t9R6FynM6a2mYqmNZpqOVLN+fPifLnKpeekITw5Ah1qd1BLHEfGxSIjPRXPn96DbRtn8XsmemQ0GJEI/176Gyyw4eoyRi3tE1JRFX5H68r3NCHKUuQ9rLtTWgFwuVzY29sjNDRUfEwkEiE0NBSOjo4yr3F0dJQoDwAXLlwosrwslpaWMDc3R2xsrMTxp0+fon79+go8A93VzqwdDvQ+gPXO62FZ1VJ8/Pr76xj4v4GYdW0W4tKKTnbLiycpT3Av/7VcZV1CemJS6CT8dmqs2uIxrGwEy0atpY6z2Gy4uYyj5I8QQohK6HQPIAD4+Phg5MiRaNu2LRwcHLBhwwZkZmaKb816eXmhTp064tuyU6dOhZOTE9atW4c+ffrg0KFDCA8Px44dO8R1pqSkIC4uDvHxBWvjFSZ6ZmZmMDMzA4vFwsyZM+Hr6wtbW1vY2dlhz549ePLkCYKDgzX8CqgPi8VCj/o90NWiK/734n/Ycn+LeHu5s6/O4sLrC/jJ6ieMtx2PWga1tByt6r0SvMKg/w1S+LobaZFqiIYQQgjRHJ1PAIcMGYKPHz9i4cKFSEhIgJ2dHc6dOyee6BEXFyexNl/Hjh0RFBSEBQsWYN68ebCyskJISAhatvy6QPKpU6fECSQADB06FADg6+uLRYsWAQCmTZuG7OxsTJ8+HSkpKbC1tcWFCxfQqFEjDTxrzdJj62GA1QD0btgbR2KPYOeDnfic8xn5TD6OPD2Cky9OwrOZJ35p+Quq86trO1wJDMPg4aeHeJ/xHg8+PsCI5iNgZijf7f4113VnAWX7LG1HQAghpCLR+XUAyzJtrgNYGhm5Gdj3aB8CHwYiK/9rZmKob4gRzUfAq7kXqnCraDHCr3Y+2IlN976ug9fKpBUO9D4g17Wjd7RHOE83Mq8OOVWwc9xNWpOMEA2jzxwpTxTJO3R6DCDRjsrcyphoNxFnPc5iRPMR4LILZqZm5mVi2/1tcDvmhl3Ru5CVp/3k6dvkDwAefHwg/8X0pw8hhJAKihJAUqQa/BqY1W4W/vnpHwxqMgh6rIIRA2m5adgQuQG9j/fGgccHkCvM1XKk8hHkCHDi2QnxOEddwgL1PBBCCNEcSgBJicwMzbDQcSFODTiF/o36g80qeNt8yv6EVXdWoc+JPgh+Gow8UZ5G48oR5pRc6Buzr8/GwpsL4XXWCwDwhS1SR1hKsWWXj9nlhBBCygZKAIncLKpYYHmn5Tje/zh61O8hPp6QmYDFYYvR/0R/nHpxCvmikrdVU4Xs/Gy5ywpyBLjx/gYA4EPmBwDAQ67816vbqCFbtB0CIYSQCoQSQKKwRtUaYb3zehzuexid63QWH3+X8Q7z/5uP1vtaw2aPDd5nvNdilF8FPw1Gp0NF7+GsLYOyTLAosxWOt94BwyrVtB0OIYSQCoQSQKK05sbNscVlC/b12of2tdtLnXc75gabPTbIzMtUS/u5ufL14C0Ok96fWRdYVbOHx68HYNVK/kXKCSGEEFWgBJCUml0tO/zd82/sdt0t83yHoA7Y83CPyieLZGWkq7Q+Tfr902cM+NFX22EQQgipoCgBJCrTzqwdokdGyzy3Nnwt7Pfbw2aPjcKTNxR1+MlhHH5yGL9e/BWxKbElX6Bh5nn5GDT+MfgGlbUdCiGEkApK53cCIWVPYRJos8dG5vm2+9sCACJHREKfra90O0wRC/ktu71M/HNkkuxt2wQ5AqXbLa1RGfVgYKgbC2kTQgipmCgBJGoTPTIaDMMg4GEA/or4S+p8m31tpI498HqgwGr8Ja/kXNT4Q21OCnEZslNrbRNCCCEA3QImasZisfBLy1/wwEu+HTpa7W2FN2lv5CpbVA+gruqamYUjbzgwqV1P26EQQgip4CgBJBrBYrEQPTIa0SOjUUmvUrFl+57oiz0P9yBPWPzC0klfklUZolr9/ukzBhkMg/XCKG2HQgghhIDFMEzZ6kYpQxTZlLmiYRgGrfa2UuiaU+6n0MCogfhx672tkc9oZtHp0vgjOQUmVYag67j1JZaljekJ0Sz6zJHyRJG8g3oAiVZ82yMor/4h/ZGWm4Zh/wyDzR6bMpH8AYAZqys6jV6t7TAIIYQQMZoEQrSupFnD3/rh4A/qDkelLsS9h9kf8ie5hBBCiCZQDyDRGYr2CJYFMUaDtR0CIYQQIoV6AInOUaRHUFe1ys7Br6kC/DB7m7ZDIYQQQqRQDyDRWYU9gks6LtF2KAoZmJaOAx8S0TyL/r4ihBCimygBJDpvgNUARI+MxsWBF7UdikwzPn1G5Ks4uGZkwinrC6Z/TkWo0U+I9wjRdmiEEEKITNRFQcoMU0NT8e4iii4how5m+fk4+zZe/CFa+/ETAOAhuym6Tw/QXmCEEEJICagHkJQ5hUvIHK0xBkHvE7QSw5KPn3Dhm+TvussppDEGyGJ4gLu/VmIihBBC5EU9gKTMSn37EB1ycxH9Kg4AsNy4Og5VraL2diNfxUH//3/OYnh40PcfdG7XDq8bPUBOTjZaWFqoPQZCCCGkNKgHkJRZLEYo8Xj+p8+IfhWH/968K3Xd0a/i0C89E3Xy8rEv/msvY8PcPHHyBwAZvz1Gh3btAACWtU3QlJI/QgghZQD1AJIyq57grszjRiKRuFcwg8XCFzYLJkIRgqsYYnFN4xLr9UjLAACsSP4kPnbs3QdcNqyE/umZ4mOJ0+JhWs2wNE+BEEII0QpKAEmZxHxJRe3cNyWWq8wwqCws2O56YHomemdk4X+VDXGsSmU85nGlyttk52DRpxSp403y8tAkNU/8ONhqNQZS8kcIIaSMogSQlEmvTixBQyWuM2AYDEnPwJD0DInjQgA5LBYMGKbEOn7VX4otwyco0TohhBCiG2gMICmT6rw4pNL6OIBcyR8A+M+botK2CSGEEE2jBJCUPYkPwRNmllxOTVgsltbaJoQQQlSBEkBSpjCCd8DWjmpvJ6b/P9iT3wMvRWbY03Qr/DEYn5gq2Gz0u9rbJoQQQtSNxgAS3SASAmyO9PG8bODKCkDfAOgyC2/O/AVLDYRT2bwZuvnsxYN3AgyxroWkNHeciZ2HYa3MNdA6IYQQol6UABKtY7Y4gpX0qODBws8A+2vH9Jcr61HpxkYAgFAQj1cpIo0kgJZmNQEAFjUMAAD1jA0womMDDbRMCCGEqB/dAiZaxXx8+jX5A4Al1YHMr+vvZd8JFP/MubcH6QnPNRgdIYQQUj5RAki06nW89F6+OSenFfzAMKielyhxrj8nTO0x7crvpfY2CCGEEG2iW8BEq9ip0os5856eAkQiMHFh0Mp82x5LtNEqIYQQojGUABKtqnV7pewTS6prJfmbnDsFGztbaaFlQgghRHPoFjDRqrQvudoOAQDwQlQbfXOWocfgieCwaZ0/Qggh5RslgEQ9rq8D1lkDMceLLcbIufuGOs3L80b33HX4Y+ww9LelZV4IIYSUf2UiAfT394elpSX4fD7at2+PO3fuFFv+6NGjaNasGfh8PmxsbHDmzBmJ88ePH0fPnj1hbGwMFouFqKgoqTqcnZ3BYrEk/pswgfZ/lYswHwhdAqTHA8Gjiy1qhmQNBQW0zd4q/vkdUxMNsvfDMjsIp/Rc8XpVH7RvaEy7fBBCCKkQdD4BPHz4MHx8fODr64vIyEjY2trC1dUVSUlJMsvfvHkTnp6e8Pb2xr179+Du7g53d3fExMSIy2RmZqJTp05YvXp1sW2PHTsWHz58EP+3Zs0alT63cinpCbDUWPJYuvRMX017LTLFmlHdYZl9AJbZB9ApZxNOTu6Mlyt6I2axq7bDI4QQQjSKxejCPbhitG/fHu3atYOfnx8AQCQSwcLCAlOmTMGcOXOkyg8ZMgSZmZk4ffq0+FiHDh1gZ2eHbdu2SZR9/fo1GjRogHv37sHOzk7inLOzM+zs7LBhwwalY09LS4ORkREEAgGqVq2qdD06Iy0eyBYAtayLLrPISL669PhA1TpAygvVxFaM68KWGJk3By9X9cOThDTcf5uKzlYmMK9WSe1tK4phGAgEAhgZGVFvJCEaQJ85Up4oknfo9Czg3NxcREREYO7cueJjbDYbLi4uCAuTvR5cWFgYfHx8JI65uroiJCRE4fYPHDiA/fv3w8zMDP369cMff/wBAwODIsvn5OQgJydH/DgtLQ1AwReMWvLsZxcAoaxJFDLayskAGCHANwIYBmD9f+cvi1XQQ5fyCqywzWDsfgZMWwBVzIBPL4AbGwC74WDd2a7a2POz1Zr8xYgssTp/KK6LWgEATvzaEQzDoKlpFTQ1rQJAN8Yffq/wvaKLsRFSHtFnjpQniryPdToBTE5OhlAohKmpqcRxU1NTPHnyROY1CQkJMssnJCh2G3LYsGGoX78+zM3N8eDBA8yePRuxsbE4frzoSQ0rV67E4sWLpY4LBAK1fLlUPTYG7ByBSutkRe2XPqjq5E/F3opM0Dl3AyBj4ZhNA63Rtp4RDLgsCASqfa3UgWEYZGRkAAD1RhCiAfSZI+VJYceTPHQ6AdSmcePGiX+2sbFB7dq10b17d7x48QKNGjWSec3cuXMleh/T0tJgYWEBIyMjtdwCFjJlYBCnGlwT2uCv/IF4ytRFPjjIARcAcHhcB8w+9gB5QgaXf3eCPqfsvTqFfyjQ7ShCNIM+c6Q8UeQ9rNMJYM2aNcHhcJCYKLkdWGJiIszMzGReY2ZmplB5ebVv3x4A8Pz58yITQB6PBx6PJ3W8cBaxqq3PHQAmL1vmue/7G1kAnDlRuCdqjFSmMlgA9JGP1uxniBA1we/6R8Vlo0WWOCXsiB6cCDiwY/FSZIaGbPVN5DiY3xX+Qne8Y2qiEnLAgQgZqARZPXqyBE9wRFvLGgCAKzO7qi1OTfl25jkhRP3oM0fKi3KTAHK5XNjb2yM0NBTu7u4ACiaBhIaGYvLkyTKvcXR0RGhoKKZNmyY+duHCBTg6OpYqlsKlYmrXrl2qelTpAKsPBMI8uctvFfaXPigs+J+fcIDUqZ3CvsqGprQv4Bd5bqZrU2Tm5MOmjhE6NqqJqpUK3r70pU0IIYQoRqcTQADw8fHByJEj0bZtWzg4OGDDhg3IzMzE6NEF68t5eXmhTp06WLmyYEuxqVOnwsnJCevWrUOfPn1w6NAhhIeHY8eOHeI6U1JSEBcXh/j4eABAbGwsgILeQzMzM7x48QJBQUHo3bs3jI2N8eDBA0yfPh1dunRBq1atNPwKFG1cl4bIyRNKHvwuGSp8FJuQjocfBHBrYYZKXD2JcywW8O7zF9x+9QlvU76oN2gZhrWvh6ndrWBalY88oQhCEQOeHpsSO0IIIURNdD4BHDJkCD5+/IiFCxciISEBdnZ2OHfunHiiR1xcHNjsr2O9OnbsiKCgICxYsADz5s2DlZUVQkJC0LJlS3GZU6dOiRNIABg6dCgAwNfXF4sWLQKXy8XFixfFyaaFhQU8PDywYMECDT1r+Uzq2ljbIaicPocNfY62oyCEEELKN51fB7AsK3frABK1ozXJCNEs+syR8kSRvKPsTZMkhBBCCCGlQgkgIYQQQkgFQwkgIYQQQkgFQwkgIYQQQkgFQwkgIYQQQkgFQwkgIYQQQkgFQwkgIYQQQkgFQwkgIYQQQkgFQwkgIYQQQkgFo/NbwZVlhZuspKWlaTkSUlYwDIO0tDSwWCzalYAQDaDPHClPCvMNeTZ5owRQjdLT0wEAFhYWWo6EEEIIIRVFeno6jIyMii1DewGrkUgkQnx8PKpUqaK2vyzbtWuHu3fvqqVuTbaninpLU4cy1ypyjbxl09LSYGFhgbdv31b4/aM1/d5WBH3uVFOHLnzu6DMniT536m9PnZ87hmGQnp4Oc3NzsNnFj/KjHkA1YrPZqFu3rlrb4HA4Gv3SUld7qqi3NHUoc60i1yhaf9WqVSv8LyNNv7cVQZ871dShS587+swVoM+d+ttT9+eupJ6/QjQJpIybNGlSuWhPFfWWpg5lrlXkGk3/O5UHuvya0edONXXQ50736PJrRp871dZBt4AJ0SFpaWkwMjKCQCDQ2b/CCSlP6DNHKirqASREh/B4PPj6+oLH42k7FEIqBPrMkYqKegAJIYQQQioY6gEkhBBCCKlgKAEkhBBCCKlgKAEkhBBCCKlgKAEkhBBCCKlgKAEkhBBCCKlgKAEkpIw4ffo0mjZtCisrK/z999/aDoeQCmHAgAGoXr06Bg4cqO1QCFEpWgaGkDIgPz8fzZs3x+XLl2FkZAR7e3vcvHkTxsbG2g6NkHLtypUrSE9Px549exAcHKztcAhRGeoBJKQMuHPnDlq0aIE6deqgcuXK6NWrF86fP6/tsAgp95ydnVGlShVth0GIylECSIgGXLt2Df369YO5uTlYLBZCQkKkyvj7+8PS0hJ8Ph/t27fHnTt3xOfi4+NRp04d8eM6derg/fv3mgidkDKrtJ87QsozSgAJ0YDMzEzY2trC399f5vnDhw/Dx8cHvr6+iIyMhK2tLVxdXZGUlKThSAkpP+hzR0jRKAEkRAN69eqFZcuWYcCAATLPr1+/HmPHjsXo0aPRvHlzbNu2DQYGBti9ezcAwNzcXKLH7/379zA3N9dI7ISUVaX93BFSnlECSIiW5ebmIiIiAi4uLuJjbDYbLi4uCAsLAwA4ODggJiYG79+/R0ZGBs6ePQtXV1dthUxImSfP546Q8kxP2wEQUtElJydDKBTC1NRU4ripqSmePHkCANDT08O6devQtWtXiEQizJo1i2YAE1IK8nzuAMDFxQX3799HZmYm6tati6NHj8LR0VHT4RKicpQAElJG9O/fH/3799d2GIRUKBcvXtR2CISoBd0CJkTLatasCQ6Hg8TERInjiYmJMDMz01JUhJRv9LkjFR0lgIRoGZfLhb29PUJDQ8XHRCIRQkND6VYTIWpCnztS0dEtYEI0ICMjA8+fPxc/fvXqFaKiolCjRg3Uq1cPPj4+GDlyJNq2bQsHBwds2LABmZmZGD16tBajJqRso88dIUWjreAI0YArV66ga9euUsdHjhyJwMBAAICfnx/+/PNPJCQkwM7ODps2bUL79u01HCkh5Qd97ggpGiWAhBBCCCEVDI0BJIQQQgipYCgBJIQQQgipYCgBJIQQQgipYCgBJIQQQgipYCgBJIQQQgipYCgBJIQQQgipYCgBJIQQQgipYCgBJIQQQgipYCgBJISQUlq0aBHs7Oy0HYYYi8VCSEiIwtfFxsbCzMwM6enpqg/qG8nJyahVqxbevXun1nYIIUWjBJAQUiZs27YNVapUQX5+vvhYRkYG9PX14ezsLFH2ypUrYLFYePHihYaj1CxVJ55z587FlClTUKVKFZXVKUvNmjXh5eUFX19ftbZDCCkaJYCEkDKha9euyMjIQHh4uPjY9evXYWZmhtu3byM7O1t8/PLly6hXrx4aNWqkjVDLpLi4OJw+fRqjRo3SSHujR4/GgQMHkJKSopH2CCGSKAEkhJQJTZs2Re3atXHlyhXxsStXruDHH39EgwYNcOvWLYnjXbt2BQDs27cPbdu2RZUqVWBmZoZhw4YhKSkJACASiVC3bl1s3bpVoq179+6BzWbjzZs3AIDU1FSMGTMGJiYmqFq1Krp164b79+8XG+/ff/8Na2tr8Pl8NGvWDFu2bBGfe/36NVgsFo4fP46uXbvCwMAAtra2CAsLk6hj586dsLCwgIGBAQYMGID169ejWrVqAIDAwEAsXrwY9+/fB4vFAovFQmBgoPja5ORkDBgwAAYGBrCyssKpU6eKjffIkSOwtbVFnTp1xMdk9TBu2LABlpaW4sejRo2Cu7s7VqxYAVNTU1SrVg1LlixBfn4+Zs6ciRo1aqBu3boICAiQqKdFixYwNzfHiRMnio2LEKIelAASQsqMrl274vLly+LHly9fhrOzM5ycnMTHv3z5gtu3b4sTwLy8PCxduhT3799HSEgIXr9+Le7lYrPZ8PT0RFBQkEQ7Bw4cwA8//ID69esDAAYNGoSkpCScPXsWERERaNOmDbp3715k79WBAwewcOFCLF++HI8fP8aKFSvwxx9/YM+ePRLl5s+fj99//x1RUVFo0qQJPD09xbe4b9y4gQkTJmDq1KmIiopCjx49sHz5cvG1Q4YMwYwZM9CiRQt8+PABHz58wJAhQ8TnFy9ejMGDB+PBgwfo3bs3hg8fXmxv2/Xr19G2bdtiX/+iXLp0CfHx8bh27RrWr18PX19f9O3bF9WrV8ft27cxYcIEjB8/XmrMn4ODA65fv65Um4SQUmIIIaSM2LlzJ2NoaMjk5eUxaWlpjJ6eHpOUlMQEBQUxXbp0YRiGYUJDQxkAzJs3b2TWcffuXQYAk56ezjAMw9y7d49hsVji8kKhkKlTpw6zdetWhmEY5vr160zVqlWZ7OxsiXoaNWrEbN++nWEYhvH19WVsbW0lzgUFBUmUX7p0KePo6MgwDMO8evWKAcD8/fff4vMPHz5kADCPHz9mGIZhhgwZwvTp00eijuHDhzNGRkbix9+3WwgAs2DBAvHjjIwMBgBz9uxZma8JwzCMra0ts2TJEoljsur/66+/mPr164sfjxw5kqlfvz4jFArFx5o2bcp07txZ/Dg/P58xNDRkDh48KFHX9OnTGWdn5yJjIoSoD/UAEkLKDGdnZ2RmZuLu3bu4fv06mjRpAhMTEzg5OYnHAV65cgUNGzZEvXr1AAARERHo168f6tWrhypVqsDJyQlAwZg3ALCzs4O1tbW4F/Dq1atISkrCoEGDAAD3799HRkYGjI2NUblyZfF/r169kjnJJDMzEy9evIC3t7dE+WXLlkmVb9Wqlfjn2rVrA4D49nRsbCwcHBwkyn//uDjf1m1oaIiqVauK65bly5cv4PP5ctf/rRYtWoDN/vrrxNTUFDY2NuLHHA4HxsbGUu1XqlQJWVlZSrVJCCkdPW0HQAgh8mrcuDHq1q2Ly5cv4/Pnz+JkztzcHBYWFrh58yYuX76Mbt26AShIxlxdXeHq6ooDBw7AxMQEcXFxcHV1RW5urrje4cOHIygoCHPmzEFQUBDc3NxgbGwMoGCm8fdjDwsVjsf7VkZGBoCC8Xvt27eXOMfhcCQe6+vri39msVgACsYlqsK3dRfWX1zdNWvWxOfPnyWOsdlsMAwjcSwvL0+utuRpPyUlBSYmJkU/CUKI2lACSAgpU7p27YorV67g8+fPmDlzpvh4ly5dcPbsWdy5cwcTJ04EADx58gSfPn3CqlWrYGFhAQASs4gLDRs2DAsWLEBERASCg4Oxbds28bk2bdogISEBenp6EpMfimJqagpzc3O8fPkSw4cPV/p5Nm3aFHfv3pU49v1jLpcLoVCodBvfat26NR49eiRxzMTEBAkJCWAYRpygRkVFqaQ9AIiJiZFawocQohl0C5gQUqZ07doV//33H6KiosQ9gADg5OSE7du3Izc3VzwBpF69euByudi8eTNevnyJU6dOYenSpVJ1WlpaomPHjvD29oZQKET//v3F51xcXODo6Ah3d3ecP38er1+/xs2bNzF//nyZySRQMAFj5cqV2LRpE54+fYro6GgEBARg/fr1cj/PKVOm4MyZM1i/fj2ePXuG7du34+zZs+JErDDuV69eISoqCsnJycjJyZG7/u+5uroiLCxMIqF0dnbGx48fsWbNGrx48QL+/v44e/as0m18KysrCxEREejZs6dK6iOEKIYSQEJImdK1a1d8+fIFjf+vvft3PS2O4zj+Ut9FZtkUfVOUrErIoBhR8mNSxEIpf4ThZCODwmQwMfhRBrPNZrUb1PkD7h1uie+93b739tXVPc/H+B7e73O2V58+533e3+VyuW71eDwu0zRv62KkHydYk8lE8/lcgUBA3W5XhmH8sm+5XNbxeFQmk5Hdbr/VbTabVquVYrGYKpWKfD6fCoWCzufzw/x71WpVo9FI4/FYwWBQ8Xhck8lEHo/n0+8ZiUQ0HA7V6/UUCoW02WzUbrcf7unlcjmlUiklEgk5nU7NZrNP9/8onU7r7e1Nu93uVvP7/RoMBur3+wqFQjocDup0On89495isZDb7VY0Gv2SfgD+jO3bxwseAICXVKvVdDqdnrY6pd/va7lcarvdPqX/vXA4rFarpVKp9PRZAH7GHUAAeFGGYSiZTMrhcGi9Xms6nT4slP5q9Xpd1+tVpmk+9Xdwl8tF2WxWxWLxaTMA/B4ngADwovL5vPb7vUzTlNfrVbPZVKPR+NePBeA/QAAEAACwGD4CAQAAsBgCIAAAgMUQAAEAACyGAAgAAGAxBEAAAACLIQACAABYDAEQAADAYgiAAAAAFkMABAAAsJjvhlcgz4B9Ik0AAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(7, 4))\n", "plt.plot(wlgrid, abs_only_rprs, lw=2, label='Absorption only')\n", "plt.plot(wlgrid, rprs_cia, lw=2, label='+ CIA')\n", "plt.plot(wlgrid, rprs_full, lw=2, label='+ CIA + Rayleigh')\n", "plt.xscale('log')\n", "plt.xlabel('Wavelength (um)')\n", "plt.ylabel('$(R_p/R_s)^2$')\n", "plt.title('Effect of adding contributions')\n", "plt.legend()\n", "plt.grid(alpha=0.2)" ] }, { "cell_type": "markdown", "id": "6f9fbd94", "metadata": {}, "source": [ "## Reading the Spectrum\n", "\n", "Each point on the transit-depth curve encodes the altitude at which the atmosphere becomes opaque at that wavelength. A higher transit depth means the planet is blocking a larger cross-section of the stellar disc — typically because a strong absorber is present at altitude.\n", "\n", "- **Molecular absorption** produces the narrow, prominent features: peaks at wavelengths where individual gases absorb strongly.\n", "- **CIA** adds a smooth continuum offset, most noticeable at the red end of the spectrum.\n", "- **Rayleigh scattering** lifts the spectrum at short wavelengths, producing a characteristic blue slope." ] } ], "metadata": { "kernelspec": { "display_name": ".venv", "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.12.3" } }, "nbformat": 4, "nbformat_minor": 5 }