{ "cells": [ { "cell_type": "markdown", "id": "6cc75a97-e0fb-458a-951b-fe97ba2b656a", "metadata": {}, "source": [ "# Installation and Quick Start\n", "\n", "## Installation\n", "\n", "To install Nimbus, clone the git repository and install it:" ] }, { "cell_type": "code", "execution_count": null, "id": "ef9706ee-bbb0-476d-8161-c6b5f8cc0aad", "metadata": {}, "outputs": [], "source": [ "git clone https://github.com/Kiefersv/Nimbus.git\n", "cd Nimbus\n", "pip install -e ." ] }, { "cell_type": "markdown", "id": "31eddb97-a05d-42f2-a903-02da37dffd24", "metadata": {}, "source": [ "## Quick Start\n", "\n", "To create a Nimbus cloud profile, you will need to provide a temperature, pressure, and K$_\\mathrm{zz}$ profile in addition to the gravity and mean molecular weight. For the cloud structure, you need to define the mass mixing ratio from which the cloud is replenished (here called the deep MMR) and a guess on the particle size, given through f$_\\mathrm{sed}$ (usually 1 is a good guess)." ] }, { "cell_type": "code", "id": "780db206-826a-4db2-aef4-3e946a09ade3", "metadata": { "ExecuteTime": { "end_time": "2026-03-16T01:43:57.066158Z", "start_time": "2026-03-16T01:43:45.823464Z" } }, "source": [ "# import nimbus\n", "import numpy as np\n", "from nimbus import Nimbus\n", "\n", "# define temperature [K] and pressure [bar] structure\n", "temperature = np.asarray([554, 572, 607, 653, 775, 951, 1073, 1111, 1540, 2654, 3000])\n", "pressure = np.asarray([1e-6, 1e-5, 1e-4, 1e-3, 1e-2, 1e-1, 1e0, 1e1, 1e2, 1e3, 1e4])\n", "\n", "# set up Nimbus object\n", "obj = Nimbus()\n", "obj.set_up_atmosphere(\n", " temperature = temperature,\n", " pressure = pressure,\n", " kzz = np.ones_like(pressure) * 1e9,\n", " mmw = 2.34,\n", " gravity = 10**2.49,\n", " species = 'SiO',\n", " deep_mmr = 10**-3,\n", ")\n", "obj.set_up_solver()\n", "\n", "# compute the cloud structure\n", "ds = obj.compute()" ], "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "===========================================================\n", " Welcome to Nimbus \n", "===========================================================\n", "[INFO] For questions contact: kiefersv.mail@gmail.com\n", "[INFO] Settings selected:\n", " -> working directory: .\n", " -> verbose: False\n", " -> analytic plots: False\n", "[INFO] Atmosphere set up with:\n", " -> pressure range: 1.00e+04 - 1.00e-06 bar\n", " -> temperature range: 3.00e+03 - 5.54e+02 K\n", " -> Kzz range: 1.00e+09 - 1.00e+09 cm2/s\n", " -> Mean molecular weight: 2.34e+00 amu\n", " -> Gravity: 3.09e+02 cm/s2\n", " -> SiO deep MMR: 1.00e-03 g/g\n", "[INFO] Solver set up.\n", "[INFO] Max itterations set to 50\n", "[INFO] Cloud structures completed in 7.85s (21 iterations).\n", "[INFO] Saved run under tag: last_run\n" ] } ], "execution_count": 1 }, { "cell_type": "markdown", "id": "1e1231be-262f-457d-8d84-5bdf66de5250", "metadata": {}, "source": [ "The resulting ds is an xarray dataset containing the cloud structure. It can be plotted like this:" ] }, { "cell_type": "code", "id": "d6e4d10b-ead0-4ac6-8d9a-e389d8237bbb", "metadata": { "ExecuteTime": { "end_time": "2026-03-16T01:45:16.184133Z", "start_time": "2026-03-16T01:45:16.012624Z" } }, "source": [ "import matplotlib.pyplot as plt\n", "\n", "# plot results\n", "fig, ax = plt.subplots(1, 3, figsize=(6,4))\n", "ax[2].semilogy(ds['qc'], ds['pressure'])\n", "ax[1].semilogy(ds['qv'], ds['pressure'])\n", "ax[0].semilogy(ds['rg'][:-2], ds['pressure'][:-2])\n", "\n", "# make a pretty plot\n", "for axi in ax:\n", " axi.set_ylim(bottom=ds['pressure'][-1], top=ds['pressure'][0])\n", "plt.subplots_adjust(wspace=0)\n", "ax[0].set_ylabel('pressure [bar]')\n", "ax[2].set_xlabel('Cloud MMR [g/g]')\n", "ax[1].set_xlabel('Gas MMR [g/g]')\n", "ax[0].set_xlabel('Particle radius [cm]')\n", "ax[2].get_xaxis().set_ticks([5e-6, 1e-5])\n", "ax[2].get_xaxis().set_ticklabels(['5e-6', '1e-5'])\n", "ax[1].get_yaxis().set_ticklabels([])\n", "ax[1].get_xaxis().set_ticks([0, 5e-4, 1e-3])\n", "ax[1].get_xaxis().set_ticklabels(['0', '5e-4', '1e-3'])\n", "ax[2].get_yaxis().set_ticklabels([])\n", "plt.show()" ], "outputs": [ { "data": { "text/plain": [ "
" ], "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiMAAAF7CAYAAAAAH4u2AAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjYsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvq6yFwwAAAAlwSFlzAAAPYQAAD2EBqD+naQAAXzNJREFUeJzt3XlYVPX+B/D3DAPDDgKyKYiKoKICQiDmWpiS4VK3rMy1q7ZoFvmzrOuSldrV1EzLm6WWdq9mFi2Wmgu5oSaCqSiCIqCyqMgu65zfHzijI9sMzsyZgffreeapOfOdcz5n5jB+zneVCIIggIiIiEgkUrEDICIiotaNyQgRERGJiskIERERiYrJCBEREYmKyQgRERGJiskIERERiYrJCBEREYmKyQgRERGJiskIERERiYrJCBEREYmKyQgRERGJqlUkI+np6Rg8eDC6d++Onj17orS0VOyQiIiI6A5Ja1gob+DAgfjggw/Qv39/5Ofnw97eHjKZTOywiIiICECL/xf57NmzMDc3R//+/QEATk5OIkdERERE9xK9mebAgQOIjo6Gp6cnJBIJYmNj65RZs2YNfHx8YGlpifDwcBw/flzj/aempsLW1hbR0dHo3bs3Fi1apMPoiYiI6EGJXjNSWlqKwMBATJ48GU8++WSd17du3YqYmBisXbsW4eHhWLlyJYYOHYqUlBS4uroCAIKCglBdXV3nvbt370Z1dTUOHjyIpKQkuLq6YtiwYXjooYcwZMgQvZ8bERERNc2o+oxIJBL8+OOPGDVqlGpbeHg4HnroIaxevRoAoFAo4OXlhRkzZuDtt99ucp/x8fFYsGABdu3aBQBYunQpAOD//u//6i1fUVGBiooK1XOFQoH8/Hw4OztDIpE099SI6rj/WqupqUFWVhYCAgJgZmYmYmTU0vBaI7EIgoDi4mJ4enpCKm24MUb0mpHGVFZWIiEhAXPmzFFtk0qliIyMRHx8vEb7eOihh5CXl4dbt27BwcEBBw4cwLRp0xosv3jxYrz33nsPHDsRERHVysrKQvv27Rt83aiTkRs3bqCmpgZubm5q293c3HD+/HmN9iGTybBo0SIMGDAAgiDgsccewxNPPNFg+Tlz5iAmJkb1vLCwEN7e3sjKyoK9vb1a2YH/3o+bpZXY/M8wBHm10eLMiOrerV67dg3h4eF1rrXEzFsY99VxeDlZ4feZA8QIlUwcrzUSS1FREby8vGBnZ9doOaNORnQlKioKUVFRGpWVy+WQy+V1ttvb29dJRnp18sCfF64joxgYcN9rRM11/7Vma1cDqdwaMkvrOtcg0YPgtUaG0lQ3B9FH0zTGxcUFZmZmyM3NVduem5sLd3d3kaK6K8Cz9o81+VqhyJEQERGZLqNORiwsLBASEoK9e/eqtikUCuzduxcREREiRlYrwNMBAHD2WpHIkRAREZku0ZtpSkpKkJaWpnqenp6OpKQkODk5wdvbGzExMZgwYQJCQ0MRFhaGlStXorS0FJMmTRIx6lrKmpHzOcWoqlHA3MyoczsiIiKjJHoycuLECQwePFj1XNl5dMKECdi4cSPGjBmD69evY968ecjJyUFQUBB27txZp1OrGLydrGErl6GkohqXrpfC373xDjpERERUl+jJyKBBg9DUVCfTp0/H9OnTDRSR5qRSCbp52OGvy7dw9lohkxEiIqJmYLvCA2K/ESIiogfDZOQBdb/Tb+QsR9QQERE1C5ORB3R3eG9Rk81NREREVBeTkQfUxdUO5mYSFJVX48qt22KHQ0REZHKYjDwgC5kUfm61HVfZVENERKQ9JiM6cG9TDREREWmHyYgOdPdQdmJlMkJERKQtJiM6ENCOw3uJiIiai8mIDnTzsIdEAuQUleNmSUXTbyAiIiIVJiM6YCuXwcfZBgBrR4iIiLTFZERH7k5+xmSEiIhIG0xGdEQ1oiabyQgREZE2mIzoyN01ajjXCBERkTaYjOiIcnhv+o1SlFZUixwNERGR6WAyoiNt7eRwtZNDEIDzOWyqISIi0hSTER0KYCdWIiIirTEZ0SFVv5GrTEaIiIg0xWREh1Q1I9nsxEpERKQpJiM6pKwZuZBTgqoahcjREBERmQYmIzrk5WQFO7kMlTUKpOWViB0OERGRSWAyokMSiQTd2ImViIhIK0xGdOzuiBr2GyEiItIEkxEdU/YbSWbNCBGZKEEQOwJqbZiM6Ji/mx0AsM8IEZkciaT2vwKYjZBhMRnRMV9XW0gkwM3SStwsqRA7HCIijd3JRaDgYEAyMCYjOmZlYYb2bawAAKmsHSEiEyJVVo0QGRiTET3o4lrbVMNkhIhMiaqZhp1GyMCYjOhBFzdbAEBqbrHIkRARaU5ZM6JgLkIGxmRED1Q1I7msGSEi08MOrGRoTEb0wE9ZM5LHmhEiMh3KZhrWjJChMRnRg85ta5ORGyWVyC+tFDkaIiLNKJtp2GWEDI3JiB7YyGWqETWcb4SITAU7sJJYmIzoSRfX2tqRC+zESkQmQlUzInIc1PowGdGTLpyJlYhMjGrSM9aMkIExGdET1owQkamRsM8IiYTJiJ4oa0Y48RkRmQr2GSGxMBnRE987NSPXiytQUMYRNURk/DiahsTCZERPbOUytHPkGjVEZDqUfUaYi5ChMRnRI+W08Ow3QkSm4O508ExHyLCYjOiRshMrp4UnIlNwt8+IuHFQ68NkRI+Ua9RweC8RmRLWjJChMRnRIzbTEJEpkUo56RmJg8mIHilH1OQVV6CwrErkaIiIGqfswMpshAyNyYge2Vmaw9PBEgBX8CUi48cOrCQWJiN65svJz4jIRKg6sIobBrVCTEb0zI/TwhORiVAmI6wZIUNjMqJnyk6sHFFDRMZOAs7ASuJoNclIWVkZOnTogFmzZhn0uL53hveyZoSIjJ1Ucvf/uT4NGVKrSUY+/PBD9OnTx+DHVdaM5BZVoPA2R9QQkfFSrtoLsHaEDKtVJCOpqak4f/48oqKiDH5se0tzuNvXjqhhUw0RGTO1mhHxwqBWSPRk5MCBA4iOjoanpyckEgliY2PrlFmzZg18fHxgaWmJ8PBwHD9+XKtjzJo1C4sXL9ZRxNpT1o6ksqmGiIyY5O5MI+zESgYlejJSWlqKwMBArFmzpt7Xt27dipiYGMyfPx8nT55EYGAghg4diry8PFWZoKAg9OjRo87j2rVr+Omnn+Dn5wc/Pz9DnVIdymnhObyXiIyZ5J5/EZiLkCHJxA4gKiqq0eaT5cuXY8qUKZg0aRIAYO3atdixYwfWr1+Pt99+GwCQlJTU4PuPHj2KLVu2YNu2bSgpKUFVVRXs7e0xb968estXVFSgoqJC9byoqKgZZ6XOj9PCUz3uv9aKi3l9kH5oeq3d00rDmhEyKNFrRhpTWVmJhIQEREZGqrZJpVJERkYiPj5eo30sXrwYWVlZuHz5MpYtW4YpU6Y0mIgoyzs4OKgeXl5eD3wene/MNXL5ZukD74tajvuvte7du4sdErVQml5rFrK7/yTcrqwxVHhExp2M3LhxAzU1NXBzc1Pb7ubmhpycHL0cc86cOSgsLFQ9srKyHnifyg6suUUVHC5HKvdfa8nJyWKHRC2UpteaXGYGN3s5ACAzv8yQIVIrJ3ozjSFNnDixyTJyuRxyuVynx3W988ddWa1AQVkV2thY6HT/ZJruv9Z00SRIVB9trrUOzjbILarA5ZulCPRyNEB0REZeM+Li4gIzMzPk5uaqbc/NzYW7u7tIUWlPLjOD050EJLe4XORoiIga5uNsDQDIuMmaETIco05GLCwsEBISgr1796q2KRQK7N27FxERESJGpj1Xu9q7kpxCJiNEZLw6ONsAYB83MizRm2lKSkqQlpamep6eno6kpCQ4OTnB29sbMTExmDBhAkJDQxEWFoaVK1eitLRUNbrGVLjZW+J8TjHyiiqaLkxEJJIOd2pGMlkzQgYkejJy4sQJDB48WPU8JiYGADBhwgRs3LgRY8aMwfXr1zFv3jzk5OQgKCgIO3furNOp1dgpO7HmFLFmhIiMVwcnZc0IkxEyHNGTkUGDBjU5wmT69OmYPn26gSLSD2UP9VwmI0RkxLzv1IzcKKlASUU1bOWi/zNBrYBR9xlpSdwc7g7vJSIyVg5W5qoO92yqIUNhMmIgbnbKZIQ1I0Rk3LydlCNq2ImVDIPJiIG4OzAZISLToBzey34jZChMRgxEOfHZjZIKVNcoRI6GiKhhyuG9mfmsGSHDYDJiIM42cphJJVAIwI2SSrHDISJqkHJ47+UbrBkhw2AyYiBmUolq4jM21RCRMVPWjLDPCBkKkxEDcuVcI0RkApR9RrKLylFexdV7Sf+YjBiQ+51+I3lMRojIiDnZWMBWLoMgAFdusamG9I/JiAG5sWaEiEyARCJhvxEyKCYjBqRMRjjxGREZOx8umEcGxGTEgO4mI6wZISLjppwWPjOfNSOkf0xGDIjr0xCRqeDEZ2RITEYMyJ3NNERkIji8lwyJyYgBKYf2Ft6u4nA5IjJqyg6sV2/dRhVnjSY9YzJiQPaWMliZmwFgUw0RGTc3O0vIZVJUKwRcK7gtdjjUwjEZMSCJRHJPvxE21RCR8ZJKJarVe9lvhPSNyYiBcRZWIjIVqgXz2G+E9IzJiIEpO7FyFlYiMnYcUUOGwmTEwJTNNDmFTEaIyLgpO7FyRA3pG5MRA1NNfFbMPiNEZNw6qGZhZc0I6ReTEQPjLKxEZCqUU8Jn5pdBoRBEjoZaMiYjBubuwGSEiEyDp6MlZFIJKqsV7HRPesVkxMDc7O6Mpiks550GERk1mZkU7dtYAQAuXWe/EdIfJiMG5u5gCVu5DBXVChxMuyF2OEREjQr0cgQAbDySLm4g1KIxGTEwC5kUT4e2BwBsOMw/biIybq892gUyqQR7zuXhMG+gSE+YjIhgYl8fSCRAXMp1XLxeInY4REQN6tzWFi/06QAAeP/XZNSweZn0gMmICDo42+DRrq4AgK+PXBY3GCKiJsx8tAvsLWU4n1OMbSeyxA6HWiAmIyKZ9HBHAMD3CVdQeLtK5GiIiBrWxsYCrz3aBQCwbPcFlFRUixwRtTRMRkTSt7Mz/N3sUFZZwzsNIjJ64yN84ONsjRslFVgbd1HscKiFYTIiEolEgokP+wAANh65zHZYIjJqFjIp5jzeDQCw7uAlXC24LXJE1JIwGRHRqKB2cLQ2x5Vbt7HnXK7Y4RARNeqx7m4I7+iEimoF/r3zvNjhUAvCZEREVhZmeC7MGwCH+RKR8ZNIJJj7RHdIJMBPSdeQmHlL7JCohWAyIrJxfTrATCrB0Uv5OJddJHY4RESN6tHOAU/1rp0r6f1fkyEIbGKmB8dkRGSejlYY1sMdAGtHiMg0/N9Qf1iZm+FkZgF2nM4WOxxqAZiMGIHJdzqyxiZdw82SCnGDISJqgpu9JV4a2BkAsOT38yivqhE5IjJ1TEaMQG/vNujV3gGV1Qr873im2OEQETVpyoCOcLe3xJVbt7Hh8GWxwyETx2TECEgkEky6Uzuy6WgGqmoU4gZERNQEawsZZg/zBwCs2Z+G68Ws1aXmk2lS6Oeff9Z6x0OGDIGVlZXW72utHu/pgQ93nEduUQV+O52NkUHtxA6JiKhRo4LaYeORy/j7SiFW7LmARaN7ih0SmSiNkpFRo0ZptVOJRILU1FR06tSpOTG1SnKZGV7o442Ve1Kx4fBlJiNEZPSkUgn+Nbw7nvlPPLYcz8SECB/4u9uJHRaZII2baXJycqBQKDR6WFtb6zPmFmtseAdYmEmRlFXA8ftEZBLCOjohqoc7FALwwQ4O9aXm0SgZmTBhglZNLi+88ALs7e2bHVRr1dZOjicCPQCAHcKIyGS8HdUVFmZSHEy9gbgL18UOh0yQRsnIhg0bYGenedXb559/DhcXl2YH1ZpNvrOa72+ns5FTWC5yNERETevgbKNaa+vDHefYCZ+0ptVomqqqKshkMpw5c0Zf8bR6Pdo54CGfNqhWCNh8NEPscIiINPLqYF842VggLa8EWzhFAWlJq2TE3Nwc3t7eqKnhBDf6NOlO7ch/j2dyMiEiMgkOVuZ4I7ILAGDFnlQU3q4SOSIyJVrPM/Luu+/inXfeQX5+vj7iIdSujNnO0Qr5pZX4Oema2OEQEWnkuTBv+LraIr+0Emv2p4kdDpkQrZOR1atX48CBA/D09IS/vz969+6t9qAHJzOTYlxEBwDAhiOX2TudiEyCzEyKd4d3AwBsPHwZGTdLRY6ITIVG84zcS9s5R6h5nn3ICyv3XMC57CIcS89Hn07OYodERNSkQX5t0b+LCw6m3sCS38/j8xdCxA6JTIDWycj8+fP1EQfdx9HaAk/2bo//HsvEhsPpTEaIyCRIJLUToUV9cgC/n8nB8fR8hHV0EjssMnItfm2arKwsDBo0CN27d0evXr2wbds2sUPS2KS+PgCAP5JzkZVfJm4wREQa8ne3w7Nh3gBqJ0JTKNjUTI3TOhmpqanBsmXLEBYWBnd3dzg5Oak9jI1MJsPKlSuRnJyM3bt34/XXX0dpqWm0Y3Zxs0P/Li5QCMA38ZfFDoeISGMxQ/xgK5fh7yuFiE26KnY4ZOS0Tkbee+89LF++HGPGjEFhYSFiYmLw5JNPQiqVYsGCBXoI8cF4eHggKCgIAODu7g4XFxeTGgmkXM13y19ZKK2oFjcYIiINudjK8epgXwDAv3em4HYlpymghmmdjHz77bdYt24d3nzzTchkMjz33HP48ssvMW/ePBw9elTrAA4cOIDo6Gh4enpCIpEgNja2Tpk1a9bAx8cHlpaWCA8Px/Hjx7U+DgAkJCSgpqYGXl5ezXq/GAb5uaKjiw2Ky6vxw8krYodDRKSxSQ/7oH0bK+QUleOLA5fEDoeMmNbJSE5ODnr2rF0m2tbWFoWFhQCAJ554Ajt27NA6gNLSUgQGBmLNmjX1vr5161bExMRg/vz5OHnyJAIDAzF06FDk5eWpygQFBaFHjx51Hteu3Z2jIz8/H+PHj8cXX3yhdYxikkolmHDPMF+2vRKRqbA0N8PbUV0BAGv/vIjcIi5xQfXTejRN+/btkZ2dDW9vb3Tu3Bm7d+9G79698ddff0Eul2sdQFRUFKKiohp8ffny5ZgyZQomTZoEAFi7di127NiB9evX4+233wYAJCUlNXqMiooKjBo1Cm+//Tb69u3bZNmKigrV86KiIg3PRH/+EeqFj3dfwKXrpYi/dBMP+3Ldn5bg/mutuLhYxGioJRPzWhve0wMbOlxGQsYtfLDjHD59LthgxybToXXNyOjRo7F3714AwIwZMzB37lx06dIF48ePx+TJk3UaXGVlJRISEhAZGanaJpVKERkZifj4eI32IQgCJk6ciEceeQTjxo1rsvzixYvh4OCgehhDk46tXIZHurkCAJKyCsQNhnTm/mute/fuYodELZSY15pEIsH86O4wk0rwy6lriE1kZ1aqS+tkZMmSJXjnnXcAAGPGjMGBAwfw8ssv4/vvv8eSJUt0GtyNGzdQU1MDNzc3te1ubm7IycnRaB+HDx/G1q1bERsbi6CgIAQFBeH06dMNlp8zZw4KCwtVj6ysrAc6B13p4moLAEjLKxE5EtKV+6+15ORksUOiFkrsa61Xe0e89kjtujVzY89wqgKqQ+tmmvtFREQgIiJCF7HoRb9+/aBQaL6ctVwub1Zzk775Mhlpce6/1oyhSZBaJmO41l4d3BkHU6/jRMYtvL41CVun9oHMrMVPdUUaataVkJKSgunTp+PRRx/Fo48+iunTpyMlJUXXscHFxQVmZmbIzc1V256bmwt3d3edH8+YKZORi9dLuFYNEZkcmZkUK8YEwU4uQ0LGLazmQnp0D62Tke3bt6NHjx5ISEhAYGAgAgMDcfLkSfTo0QPbt2/XaXAWFhYICQlR9VEBAIVCgb179xp1bYw+dHC2gUwqQVllDbIL2SOdiEyPl5M1PhjdAwCwam8qEjJMZ84n0i+tm2lmz56NOXPmYOHChWrb58+fj9mzZ+Opp57San8lJSVIS7ubIaenpyMpKQlOTk7w9vZGTEwMJkyYgNDQUISFhWHlypUoLS1Vja5pLczNpPB2tsal66VIyyuBp6OV2CEREWltZFA7xKVcx4+JVzFzSxJ+m9kf9pbmYodFItO6ZiQ7Oxvjx4+vs/2FF15Adna21gGcOHECwcHBCA6uHe4VExOD4OBgzJs3D0BtJ9lly5Zh3rx5CAoKQlJSEnbu3FmnU2tr4NuW/UaIyPQtHBkALycrXLl1G/N/Oit2OGQEtE5GBg0ahIMHD9bZfujQIfTv31/rAAYNGgRBEOo8Nm7cqCozffp0ZGRkoKKiAseOHUN4eLjWx2kJ7u03QkRkquwszbFyTBDMpBL8mHgVP3HtmlZPo2aan3/+WfX/I0aMwFtvvYWEhAT06dMHAHD06FFs27YN7733nn6iJAAcUUNELUdIByfMeMQXK/ek4l8/nkFv7zbwcrIWOywSiUbJyKhRo+ps++yzz/DZZ5+pbXv11Vfx0ksv6SQwqqtzW9aMEFHLMX2wLw6l3uBwX9KsmUahUGj0qKnhqoz61PlOzciNkkoUlFWKHA0R0YO5f7jvmv0XxQ6JRMIU1ITYymXwcLAEwNoRImoZvJys8f6oO8N996UiIeOWyBGRGDRKRlatWoXycs3ntli7di0X/dIT9hshopZmVHA7jAryRI1CwOtbE1FcXiV2SGRgGiUjb7zxhlbJxezZs3H9+vVmB0UN68zhvUTUAi0c1QPt21ghK/825nG4b6ujUQdWQRDw6KOPQibTbI6027dvP1BQ1LDOrBkhohbI3tIcnzwbhKfXxuPHxKsY5N8WI4PaiR0WGYhG2cX8+fO12unIkSPh5OTUrICocb6qETWlIkdCRKRbtcN9u+CTvRzu29roJRkh/VH2Gcm6VYbyqhpYmpuJHBERke7MeMQXh9JuICHjFt7YmoQtHO7bKvAbNjEuthawt5RBEIBLrB0hohZGZibFyjvDfU9wuG+rwWTExEgkkrsjaji8l4haIA73bX2YjJgg1Ro17MRKRC3UqOB2GMnhvq0GkxETxJoRImoN3r9nuC9X923Zmp2MVFZWIiUlBdXV1bqMhzSgWqOGNSNE1ILZ31ndVyoBfuDqvi2a1slIWVkZXnzxRVhbWyMgIACZmZkAgBkzZmDJkiU6D5DqUtaMXLpRihqFIHI0RET6E+pTO9wXAP714xlk5ZeJHBHpg9bJyJw5c3Dq1CnExcXB0tJStT0yMhJbt27VaXBUv/ZtrGEhk6KyWoErt/iHSUQt24xHfNHb2xHFFdV4Y2sSqmsUYodEOqZ1MhIbG4vVq1ejX79+kEgkqu0BAQG4eJFDsAzBTCpBJxcbAJyJlYhaPpmZFJ88GwzbO8N9P4vjvzUtjdbJyPXr1+Hq6lpne2lpqVpyQvrFaeGJqDWpHe4bAAD4ZC+H+7Y0WicjoaGh2LFjh+q5MgH58ssvERERobvIqFG+XDCPiFqZ0cHtOdy3hdJs5bt7LFq0CFFRUUhOTkZ1dTU++eQTJCcn48iRI/jzzz/1ESPVQzXXCIf3ElEr8v6oHjhx+ZZquO/yMUFih0Q6oHXNSL9+/XDq1ClUV1ejZ8+e2L17N1xdXREfH4+QkBB9xEj18L2nmUYQOKKGiFoH5eq+HO7bsmiVjFRVVWHy5MmQSCRYt24djh8/juTkZGzevBk9e/bUV4xUj44uNpBIgKLyalwvqRA7HCIigwn1ccJ05XDf2DMcVdgCaJWMmJubY/v27fqKhbRgaW4Grza1S2unc8E8ImplXnvEF8Hejigur8Zb2/9mDbGJ07qZZtSoUYiNjdVDKKQtBytzAEBZZY3IkRARGZZydV+5TIrDaTfx86lrYodED0DrDqxdunTBwoULcfjwYYSEhMDGxkbt9ddee01nwVHjZGa1I5kqOQEQEbVCHZxtMOMRXyzbfQHv/5qMQf6uqps0Mi1aJyNfffUVHB0dkZCQgISEBLXXJBIJkxEDMjerrdiqrmH1JBG1TlMGdMIPiVdx6Xoplu1KwfujeogdEjWD1slIenq6PuKgZjC/UzNSrWDNCBG1TnKZGT4Y1QPPrzuGzccy8FRIewR5OYodFmmp2av2kviUNSOV1UxGiKj16tvZBaOD20EQgHd/PM21a0yQ1jUjkydPbvT19evXNzsY0o5MeqeZhiv3ElEr987j3bD3XC7OXivCpqMZmPRwR7FDIi1oXTNy69YttUdeXh727duHH374AQUFBXoIkRqiaqbhXQARtXJt7eSYPawrAODj3ReQW1QuckSkDa1rRn788cc62xQKBV5++WV07txZJ0GRZlTNNOzASkSE58O88X3CFSRlFWDhr8lY83xvsUMiDemkz4hUKkVMTAxWrFihi92RhmSsGSEiUpFKJfhgVA9IJcCOv7Nx4MJ1sUMiDemsA+vFixdRXV2tq92RBizu1IxUMRkhIgIA9GjngIl9a/uLzP3pDMqrOCmkKdC6mSYmJkbtuSAIyM7Oxo4dOzBhwgSdBUZNU9aMVLGZhohIJeYxP/x2OhsZN8vwWdxFxAzxEzskaoLWyUhiYqLac6lUirZt2+Ljjz9ucqQN6dbd0TSsGSEiUrKVyzAvujte+fYk1sZdxKggT3Rqayt2WNQIrZOR/fv36yMOagYLmbKZhjUjRET3iurhjkH+bRGXch1zfzqDzS+GQyKRiB0WNUDrPiO3b99GWdnd5ZozMjKwcuVK7N69W6eBUdNkUmUzDWtGiIjuJZFIsHBEDy6kZyK0TkZGjhyJb775BgBQUFCAsLAwfPzxxxg5ciQ+//xznQdIDTNnB1YiogZ5O1tjxiO+AID3fz2HwttVIkdEDdE6GTl58iT69+8PAPj+++/h7u6OjIwMfPPNN1i1apXOA6SG3Z30jM00RET1mTKgEzq1tcGNkgp8vDtF7HCoAVonI2VlZbCzswMA7N69G08++SSkUin69OmDjIwMnQdIDZOZsc8IEVFjlAvpAcCmoxk4lVUgbkBUL62TEV9fX8TGxiIrKwu7du3CY489BgDIy8uDvb29zgOkhrGZhoioaWoL6cWeRg3X8zI6Wicj8+bNw6xZs+Dj44Pw8HBEREQAqK0lCQ4O1nmA1DBVMw2H9hIRNeqdx7vB3lKGM1eLsCn+stjh0H20Tkb+8Y9/IDMzEydOnMDOnTtV2x999FFOB29g5mymISLSyL0L6S3jQnpGp1nTwbu7uyM4OBhSqRRFRUWIjY2FnZ0dunbtquv4qBEc2ktEpLnnw7wR5OWIkopqvP9rstjh0D20TkaeeeYZrF69GkDtnCOhoaF45pln0KtXL2zfvl3nAVLDlDUjHE1DRNS0exfS+5UL6RkVrZORAwcOqIb2/vjjjxAEAQUFBVi1ahU++OADnQdIDVMmI5WsGSEi0si9C+nN40J6RkPrZKSwsBBOTk4AgJ07d+Kpp56CtbU1hg8fjtTUVJ0HSA2TqeYZYTJCRKSpmMf84GYvx+WbZfg87qLY4RCakYx4eXkhPj4epaWl2Llzp2po761bt2BpaanzAHXh119/hb+/P7p06YIvv/xS7HB0xkLZTMNhakREGrOVyzA/OgAA8HncRVy6XiJyRKR1MvL6669j7NixaN++PTw8PDBo0CAAtc03PXv21HV8D6y6uhoxMTHYt28fEhMTsXTpUty8eVPssHRCWTNSWc2aESIibUT1cMdAv7aorFFg3k9nIQi8qROT1snIK6+8gvj4eKxfvx6HDx+G9M4y9p06dTLKPiPHjx9HQEAA2rVrB1tbW0RFRbWYRf1kUtaMEBE1h0QiwcKRAZDLpDiUdoML6YmsWUN7Q0NDMXz4cFy9ehXV1dUAgOHDh+Phhx/WaXBAbY1LdHQ0PD09IZFIEBsbW6fMmjVr4OPjA0tLS4SHh+P48eOq165du4Z27dqpnrdr1w5Xr17VeZxisJBxaC8RUXN1cLbB9MG1C+l9sOMcisq5kJ5YmrU2zYsvvghra2sEBAQgMzMTADBjxgwsWbJE5wGWlpYiMDAQa9asqff1rVu3IiYmBvPnz8fJkycRGBiIoUOHIi8vT+exGBtVzQiH9hIRNcvUgbUL6V0vrsDHu7iQnli0TkbmzJmDU6dOIS4uTq3DamRkJLZu3arT4AAgKioKH3zwAUaPHl3v68uXL8eUKVMwadIkdO/eHWvXroW1tTXWr18PAPD09FSrCbl69So8PT11HqcYuDYNEdGDuXchvW+4kJ5otE5GYmNjsXr1avTr1w8SiUS1PSAgABcvGnaIVGVlJRISEhAZGanaJpVKERkZifj4eABAWFgYzpw5g6tXr6KkpAS///47hg4d2uA+KyoqUFRUpPYwVsq1aZiMmKb7r7Xi4mKxQ6IWitda4+5dSO+1LYkoKKsUO6RWR+tk5Pr163B1da2zvbS0VC05MYQbN26gpqYGbm5uatvd3NyQk5MDAJDJZPj4448xePBgBAUF4c0334Szs3OD+1y8eDEcHBxUDy8vL72ew4Mw8MdNOnb/tda9e3exQ6IWitda0+Y+0R3t21gh42YZpv83kfM3GZjWyUhoaCh27Niheq5MQL788kvVCr7GZsSIEbhw4QLS0tIwderURsvOmTMHhYWFqkdWVpaBomw+9hgxTfdfa8nJXCuD9IPXWtOcbCywbnworC3McCjtBhb9dl7skFoVmbZvWLRoEaKiopCcnIzq6mp88sknSE5OxpEjR/Dnn3/qI8YGubi4wMzMDLm5uWrbc3Nz4e7u3qx9yuVyyOVyXYRnAKwaMWX3X2vG3CRIpo3Xmma6edhj+TNBeGlzAtYfTkdXDzs8E2q8teMtidY1I/369cOpU6dQXV2Nnj17Yvfu3XB1dUV8fDxCQkL0EWODLCwsEBISgr1796q2KRQK7N2712hrafSBc/UQEenGsB7ueD2yCwDgXz+eQULGLZEjah20qhmpqqrCtGnTMHfuXKxbt05fMakpKSlBWlqa6nl6ejqSkpLg5OQEb29vxMTEYMKECQgNDUVYWBhWrlyJ0tJSTJo0ySDxiYl9RoiIdO+1R7ogJacYv5/JwbRNCfhlxsPwcLASO6wWTauaEXNzc2zfvl1fsdTrxIkTCA4ORnBwMAAgJiYGwcHBmDdvHgBgzJgxWLZsGebNm4egoCAkJSVh586ddTq1tmScxpiISHekUgmWPR2Iru52uFFSganfJHB1Xz3Tuplm1KhR9c6Cqi+DBg2CIAh1Hhs3blSVmT59OjIyMlBRUYFjx44hPDzcYPEREVHLYyOXYd34UDjZWOD01ULM/v5v3vjpkdYdWLt06YKFCxfi8OHDCAkJgY2Njdrrr732ms6Co8axlYaISH+8nKzx2djeeOHLY/j51DV09bDDK4N8xQ6rRdI6Gfnqq6/g6OiIhIQEJCQkqL0mkUiYjIiAuToRkX706eSMBSMC8K/YM1i6KwX+bnZ4tFvr6QZgKFonI+np6fqIg5rB0JPMERG1Ri/06YBz2UX49lgmZm5JQuyrfeHraid2WC1Ks1btVVL23yCR8SsgItKr+dEBCO/ohJKKavzz6xMoLOMKv7rUrGTkq6++Qo8ePWBpaQlLS0v06NEDX375pa5joyawXoSIyDAsZFJ8NrY32jla4fLNMkz/30lOGa9DWicj8+bNw8yZMxEdHY1t27Zh27ZtiI6OxhtvvKEabkuGxYoRIiL9c7aVY934UFiZm+Fg6g0s/p1TxuuK1n1GPv/8c6xbtw7PPfecatuIESPQq1cvzJgxAwsXLtRpgNQwdhkhIjKs7p72WP5MIF7+9iS+OpSOru52eJpTxj8wrWtGqqqqEBoaWmd7SEgIqqurdRIUaYf9doiIDCeqpwdmPlo7Zfy7nDJeJ7RORsaNG4fPP/+8zvYvvvgCY8eO1UlQpBkJe40QEYli5qNdMDTADZU1Cry0OQHZhbfFDsmkad1MA9R2YN29ezf69OkDADh27BgyMzMxfvx4xMTEqMotX75cN1FSo1gvQkRkWFKpBMufCcJTnx/B+ZxiTP0mAdteioCluZnYoZkkrZORM2fOoHfv3gCAixcvAgBcXFzg4uKCM2fOqMpxDgwiImrJlFPGj1h9CKevFuKt7X9j5Zgg/vvXDFonI/v379dHHNQMyuudXUaIiMRRO2V8CMZ9dQw/JV1DV3d7vDyos9hhmZwHmvSMiIiotYvo7Iz5IwIAAP/edR77zueKHJHpYTLSAgjsNUJEJKpxfTrg+XBvCAIw839JSMsrFjskk8JkhIiISAcWRAcgzMcJxZwyXmtMRkwY+0gRERkPC5kUn7/AKeObg8lIC8AOrERExuH+KeOXcMp4jTAZMWEcPkZEZHyUU8YDwJeH0vF9whWRIzJ+TEZaAFaMEBEZl6ieHnjtzpTx7/xwGiczOWV8Y5iMmDDWixARGa/X75kyftqmBOQUlosdktFiMtISsGqEiMjoKKeM7+puh+vFFZi66QTKq2rEDssoMRkhIiLSE+WU8W2szfH3lUK8vf1vrrReDyYjJkzZf7VaocDOM9koLueYdiIiY+PlZI01Y3vDTCpBbNI1LPn9PGtI7sNkxITZymUwN5NAIQAvbT6J4IV/4Nkv4rH2z4s4n1PE7JuIyEj07eyCBdHdAQD/OXAJj604gN1nc/g7fYfWC+WR8bCzNMe2l/oiNvEq/rxwHek3SnH0Uj6OXsrHkt/Pw8PBEgP92mKQvyse9nWGnaW52CHTA+BvFpFpGxfhA3src3y44xwy88swdVMC+ndxwbwnuqOLm53Y4YmKyYiJC/JyRJCXIwDg8o1SxKXkIe7CdcRfvInswnJs+SsLW/7KgkwqQahPGwzyd8Ug/7bwd7PjPCUmQrU6M3sqE5m8kUHtENnNDWv2p+HLg+k4mHoDwz45iPERHfB6pB8crFrnTaNEYB1Ro4qKiuDg4IDCwkLY29uLHY7GyqtqcPTSTcSlXEdcSh4u3yxTe521JsbnypUr8PLyqnOtJWbewujPjqB9GysceusRESOklqKha40MK+NmKT7YcQ5/JNeu8utkY4FZj/ljzENeMJO2jJtFTf8NZc1IC2VpbnanFsQVQABrTUyY8rvgbQNRy9LB2QbrxofiYOp1vPdLMtLySvDOj6fx7bEMzI8OQFhHJ7FDNBjWjDTBVGtGGqNJrckg/7YY6MdaE0Nq6G41KasAo9YcRjtHKxx+mzUj9OBYM2J8qmoU2BSfgRV7LqC4vBoAEB3oiTlRXeHpaCVydM3HmhFqkCa1Jv87noX/HWetiTHgp03U8pmbSTG5X0eMDPLEst0XsOWvTPxy6hr+SM7BK4N8MXVAJ1iam4kdpt6wZqQJLbFmpDGsNRFPQ3erf18pwIjVh+HpYIkjcx4VMUJqKVgzYvzOXC3Ee7+cxV+Xa9e0aedohXeHd0NUD3eTuiFkzQg1C2tNjI/kTt0I7xqIWo8e7Rzw3bQI/PJ3Nhb/dg5XC27jlW9PIqKTM+aP6I6u7i0riWQyQo3ycbHBRJeOmPhwx3prTe6f14S1JrqnGtrLbISoVZFIJBgR6InIbq5YG3cR/zlwCfGXbuLxTw5ibHgHxAzxQxsbC7HD1AkmI6Qx1pqIi/OMELVO1hYyxDzmj6dDvbDot3P4/UwONh3NwC9/X0PMED88H+YNmZlpT6jOPiNNaG19RpqLfU0eXEPt+GevFWL4qkNwtZPj+LuRIkZILQX7jJi2Ixdv4L2fk5GSWwwA6Opuh3nR3dG3s4vIkdXFPiNkUKw10T/eNRARULvOzY7X+uF/xzPx8R8XcD6nGM+vO4aoHu545/Fu8HKyFjtErTEZIb1gXxPdkXBwLxHdR2YmxbgIHzzRyxMr9lzA5qMZ+P1MDvadz8O0AZ3w0qDOsLYwnX/iTSdSMlmsNXkw7MBKRA1pY2OBhSN74Plwb7z3czLiL93Eqn1p2JZwBXMe74boXh4m8RvKPiNNYJ8R/WJfk7saasc/n1OEYSsPwsXWAif+NUTECKmlYJ+RlkkQBOw8k4MPdtQOBQaAMB8nzIvujh7tHESJiX1GyCSw1qRpqnlGeNtARI2QSCSI6umBwV1d8cWBS/gsLg3HL+cjevUhPPuQN2Y95gdnW7nYYdaLNSNNYM2IeMqrahB/6Sb+bCW1Jg3drV7ILcZjKw7AycYCJ+eyZoQeHGtGWodrBbex+Pfz+OXUNQCAnaUMr0f6YXxEB5gbaCgwa0bI5Fmam2GwvysGa1lrMtjfFX5uti2u1oT3DUSkDU9HK3z6XDDG9emA9345i7PXivD+r8n43/FMzHuiOwb4tRU7RBXWjDSBNSPGqSXWmjR0t5qaW4whKw7A0docSfMeEzFCailYM9L61CgEfHciC0t3pSC/tBIAENnNDXOf6IYOzjZ6Oy5rRqhFa06tyeA7fVNMrdbEhEIlIiNlJpXguTBvPN7TA5/sScU38Zex51wuDly4jhf7d8T0wb6wkYuXErBmpAmsGTE9plpr0tDdalpeCSKX/wkHK3Ocms+aEXpwrBmhtLxivPdLMg6m3gAAuNrJ8XZUV4wKagepVHd3QKwZoVarpdWa3J1nhPcNRKQbvq52+GZyGPacy8MHO5KRcbMMMd+dwqajGVgQHYBAL0eDxsNkhFq8+2eDvb/WRDkb7GIjnQ1WmRoxFSEiXZJIJBjS3Q0D/Fzw1aF0rN6XhsTMAoxccxhPh7TH/w3zh6udpUFiYTJCrYpJ15owGyEiPZDLzPDKIF881bs9Pvr9PH5IvIptCVfw+5kcvPaoLyb27QgLmX6HArf4PiNZWVkYN24c8vLyIJPJMHfuXDz99NMav599RloPsfuaNNSOn36jFIOXxcFWLsOZ94bq9JjUOrHPCDXmZOYtvPfzWZy6UggA6ORig7lPdMfgrq5a70vTf0NbfDKSnZ2N3NxcBAUFIScnByEhIbhw4QJsbDQbysRkpPW6v9akolqhek0ftSYN/QNx+UYpBi2Lg42FGc4uHPZAxyACmIxQ0xQKAd+fvIJ/70zBjZIKAMBg/7aY+0R3dGprq/F+2IH1Dg8PD3h4eAAA3N3d4eLigvz8fI2TEWq9jKWviaoDq872SETUOKlUgmdCvRDVwx2r96Vh/eF07E+5jkNpBzDp4Y6Y8YivTn/nDDMf7AM4cOAAoqOj4enpCYlEgtjY2Dpl1qxZAx8fH1haWiI8PBzHjx+vd18JCQmoqamBl5eXnqOmlkbZ12TBiADE/d9gxM0ahAXR3THIvy3kMqmqr8lLmxMQvPAPPPtFPP7z50Wk5BQ/8CgYCYxrdA8RtR52luaY83g37Hp9AB7p6oqqGgFfHLiEwcvi8N1fWVAodHObZPQ1I6WlpQgMDMTkyZPx5JNP1nl969atiImJwdq1axEeHo6VK1di6NChSElJgavr3fat/Px8jB8/HuvWrTNk+NRCGbLW5O7QXj2cCBGRBjq1tcX6iQ9hf0oe3v8lGZdulGL29r+x+VgG5kcHIKRDmwfav0n1GZFIJPjxxx8xatQo1bbw8HA89NBDWL16NQBAoVDAy8sLM2bMwNtvvw0AqKiowJAhQzBlyhSMGzeu0WNUVFSgoqJC9byoqIhtq6SVpvqaDPBri68mhKKyslLtWrt69Sq6d+9e51rLyi9D/3/vh0QCPPuQF4K92iDY2xGd29rqdHIiarnu/11r6Foj0kTh7Sr837ZT2J2cq9o294nueLFfxzplW0WfkcrKSiQkJGDOnDmqbVKpFJGRkYiPjwdQO1HUxIkT8cgjjzSZiADA4sWL8d577+ktZmr5mqo1AWoTa02vNRdbOZxsLJBfWqkadgwAdnIZgrwdEezliGDvNgjyckQbGwu9nhuZJv6uUXMJgoD0G6VIzCxAYtYtJGYW4HxOMWrua57JvFn6QMcx6ZqRa9euoV27djhy5AgiIiJU5WbPno0///wTx44dw6FDhzBgwAD06tVL9fqmTZvQs2fPeo/BmhHSp8s3SnG7qgbdPOy1ulstKKvEsfT82h+EzFv4+0ohblfV1Nl/RxebO8lJbYLi725nsKXCyXixZoQ0VXi7CqeyClTJR1JWAQrKquqUc7OXq2ppg73bILRDm3praltFzYgm+vXrB4VC0XTBO+RyOeRyuR4jotbMx+XuKK77r7WioqIG3+dobYGhAe4YGuAOAKiuUSAlt/hOclL7o3HpeinSb9Q+fki8CgCwNJeiVztlclL7o+Fmb5gZFcl4aHOtUetRoxBwQfU7cguJWQVIyyupU85CJkWvdg6q35Bgb0d4OFjpNBaTTkZcXFxgZmaG3Nxcte25ublwd3cXKSoi/ZOZSRHg6YAATwe80KcDgNrakyTVHU3tj0txeTWOX87H8cv5qvd6OliqflCCvR0R4OkAS3MzsU6FiAzkenGFKulIyizAqSsFKKusW8Pawdla1fwb7O2Iru72ep+B1aSTEQsLC4SEhGDv3r2qphuFQoG9e/di+vTp4gZHZGCO1hYYdGcSNqB20qJLN0pVPz6JmQVIySnCtcJyXDudjR2nswEA5mYSdPewv5ugeLWBl5OV0S0YSESaq6iuQfK1IrWbkyu3btcpZyuXIdDLQdXkEuTlCGdbw7cOGH0yUlJSgrS0NNXz9PR0JCUlwcnJCd7e3oiJicGECRMQGhqKsLAwrFy5EqWlpZg0aZKIUROJTyqVwNfVFr6utng6tHZundKKavx9pVDVES0x8xZulFTi1JVCnLpSiI1Hat/rbGNxt0rWyxG9vBxhKzf6nwuiVkkQBFwtuK3WbHv2ahEqa9S7KEgkQBdXW7W+Hr6utjAzglF5Rv/rcuLECQwePFj1PCYmBgAwYcIEbNy4EWPGjMH169cxb9485OTkICgoCDt37oSbm5tYIRMZLRu5DBGdnRHR2RlA7Y/YlVu3VXdOiZkFOHutEDdLK7HnXB72nMsDUPsj5u9mp6o54dBiIvGUVd65qbinr8f14oo65ZxsLFQd2oO82qCXlwPsjWAl8vqY1GgaMXBtGjIUY1kvpLyqBsnZRXd/6DILcLWgbvUuhxabLmO51qhpCoWA9Julan+PKbl1h9bKpBJ097RX6+vh7WQtenMrR9MQUbNYmpuht3cb9PZuA6B2EqO8onJVvxPl0OLiimocTL2Bg6k3VO/l0GKiB1NYVoWkK3cTj6SsAhTerju01sPBUq2mskc70+6IzmSEiJrkam/JocVEOqb8O0q6J9G/eL3u5GFymRS92juo+nAF6WFordiYjBCR1ji0mEh7ecXldxP4RiYv9HG2Vhvd1tWj5dcwMhkhIp3g0GKiuyqqa3D2mnZ9r4LudDR1aoV9r5iMEJFecGgxtRb1jUpLvlb/0FqOSqsf/7qJyGA4tJhagrpJdQFulNQdWnt/Ut2zvQPsjHRordiYjBCRaCQSCbycrOHlZI0RgZ4AGh5afD6nGOdzirlqMRlUbXNjyT19oWqbG+8bWQuZVIIATzY3NheTESIyKhxaTGIqKKtUu9aSsgpQXF5dpxw7YusWkxEiMnocWkz6UF2jwPmcYlUzYVJmAS7dqDu01tJcil7tHdWaCXkd6RaTESIyOZoMLU7KvIUiDi2me+QVlePkneQ1MbMApxsYWtvJxaa2CfBOXw/WsOkfkxEiahE4tJjuVV6lHFp7605y2sDQWksZgu6ZQj2oPfseiYHJCBG1SBxa3HooR2WdvNPhOTGrAMnXClFVo97LVCoB/NzsVIlHb29HdHLhqCxjwL8uImo1mju0+O4/YhxabAxKKqrx95UCVZ+hpKzapPJ+LrYWCPK62yTXqz2TSmPFb4WIWi0OLTZ+CoWAi9dLVB2VEzMLcCG3uM7QWnMzCQI8HdRqtNq3YXObqWAyQkR0Dw4tFtetUmVH5Dt9PRoYWtvO0epu4uHtiO4e9uyIbMKYjBARNYFDi/WjqkaBlJxiVQ1UYlYB0usZWmtlbnZ31do7tVGu/BxbFCYjRERa4tDi5sktKr+beGQW4O+rBSivUtQp16mtjapvTrC3I/zd7CBjDVOLxmSEiEgHOLRYXe3Q2sK7tUeZt3CtsLxOOXtLGYLu9PEI9nZEkJcjHK3Z96a1YTJCRKQHrWlosSAIyMwvu9vpN6sA57KL6h1a29XdXq3jbycXG45KIiYjRESG0lKGFheXV9UmVff09cgvrW9orRy97+lk2rOdA2yMOKki8fCqICISyQMNLVbOHKrnocUKhYC06yVqfT0u5BVDuG9orYWZFAHt7NX6erRzNP3mJjIMJiNEREZE46HF5U0PLQ7wtNc6GSgsq0JCZr4q8TiVVYDiirpDa9u3sVI1IwV7O6K7pz3kstbREZd0j8kIEZGRa+7Q4gkRHfDeyB4aH+d6cQUeWRZXJ/mwtrhnaK2XI4K8HeFqx6G1pDtMRoiITExTQ4vjLlzHqawCnLlWpNV+M/PLUFxRDblMiuhAT1UfFT83Ww6tJb1iMkJE1ALcO7S4u6c9pm1KaPa+3B0ssezpQB1GR9Q4prpEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJComI0RERCQqJiNEREQkKiYjREREJKpWk4yUlZWhQ4cOmDVrltihEBER0T1aTTLy4Ycfok+fPmKHQURERPdpFclIamoqzp8/j6ioKLFDISIiovsYfTJy4MABREdHw9PTExKJBLGxsXXKrFmzBj4+PrC0tER4eDiOHz+u9vqsWbOwePFiA0VMRERE2jD6ZKS0tBSBgYFYs2ZNva9v3boVMTExmD9/Pk6ePInAwEAMHToUeXl5AICffvoJfn5+8PPzM2TYREREpCGZ2AE0JSoqqtHmleXLl2PKlCmYNGkSAGDt2rXYsWMH1q9fj7fffhtHjx7Fli1bsG3bNpSUlKCqqgr29vaYN29evfurqKhARUWF6nlhYSEAoKioSIdnRVT3Wrt27RoAXmv04EpLiqGoKEPVbQsUFRVpfK2VFBdBUVGG6nKB1yHphPI6EgSh8YKCCQEg/Pjjj6rnFRUVgpmZmdo2QRCE8ePHCyNGjKjz/g0bNghvvvlmo8eYP3++AIAPPvjggw8++NDRIysrq9F/e42+ZqQxN27cQE1NDdzc3NS2u7m54fz5883a55w5cxATE6N6rlAokJ+fD2dnZ0gkEgC1mZ6XlxeysrJgb2/f/BMwMjwvw7r/brWmpgZZWVkICAiAmZmZWlljPQex8XNp2L2fjVwu1/hae5DjaPod8HtrPQRBQHFxMTw9PRstZ9LJiLYmTpzYZBm5XA65XK62zdHRsd6y9vb2LfIPieclng4dOjT6uimcgxj4uTSsoc+mqWtNV8fR9XvI9Dg4ODRZxug7sDbGxcUFZmZmyM3NVduem5sLd3d3kaIiIiIibZh0MmJhYYGQkBDs3btXtU2hUGDv3r2IiIgQMTIiIiLSlNE305SUlCAtLU31PD09HUlJSXBycoK3tzdiYmIwYcIEhIaGIiwsDCtXrkRpaalqdI0+yOVyzJ8/v05zjqnjeRmvlnAO+sDPpWGG+myacxx+b3Q/yZ1RKkYrLi4OgwcPrrN9woQJ2LhxIwBg9erVWLp0KXJychAUFIRVq1YhPDzcwJESERFRcxh9MkJEREQtm0n3GSEiIiLTx2SEiIiIRMVkhIiIiETVKpORplb5vd+2bdvQtWtXWFpaomfPnvjtt9/UXhcEAfPmzYOHhwesrKwQGRmJ1NRUtTL5+fkYO3Ys7O3t4ejoiBdffBElJSUmf14+Pj6QSCRqjyVLlhj1ef3www947LHHVLPqJiUl1dlHeXk5Xn31VTg7O8PW1hZPPfVUnflsDEXb8zdVCxYsqHMtde3aVWf737JlCyQSCUaNGqWzfRqKJquXa+vQoUN4+OGH4ezsDCsrK3Tt2hVDhw7V+jvQJLaCggK8+uqr8PDwgFwuh5+fX52/y3vFxcXViUMikSAnJ+dBT5uMVdMrwrQsW7ZsESwsLIT169cLZ8+eFaZMmSI4OjoKubm59ZY/fPiwYGZmJvz73/8WkpOThX/961+Cubm5cPr0aVWZJUuWCA4ODkJsbKxw6tQpYcSIEULHjh2F27dvq8oMGzZMCAwMFI4ePSocPHhQ8PX1FZ577jmTP68OHToICxcuFLKzs1WPkpISoz6vb775RnjvvfeEdevWCQCExMTEOvt56aWXBC8vL2Hv3r3CiRMnhD59+gh9+/bV2XlpStvzN2Xz588XAgIC1K6l69ev62Tf6enpQrt27YT+/fsLI0eO1Mk+Dem3334T3n33XeGHH34QANRZj6s5Tp48Kfz3v/8Vzpw5I6SnpwubNm0SZDKZ4OnpqdV30FRsFRUVQmhoqPD4448Lhw4dEtLT04W4uDghKSmpwX3u379fACCkpKSoxVJTU/PA503GqdUlI2FhYcKrr76qel5TUyN4enoKixcvrrf8M888IwwfPlxtW3h4uDBt2jRBEARBoVAI7u7uwtKlS1WvFxQUCHK5XPjf//4nCIIgJCcnCwCEv/76S1Xm999/FyQSiXD16lWTPS9BqE1GVqxYoZNzqI+uz+te6enp9SYjBQUFgrm5ubBt2zbVtnPnzgkAhPj4+Ac4G+1pe/6mbP78+UJgYGCDr9+6dUt48cUXBRcXF8HOzk4YPHhwo/+gKVVXVwt9+/YVvvzyS2HChAkmmYzcq75/8MvLy4U333xT8PT0FKytrYWwsDBh//79Wu+7a9euQps2bRp8vanvoL7YVq1aJTg4OGgVmzIZuXXrltbnQKapVTXTVFZWIiEhAZGRkaptUqkUkZGRiI+Pr/c98fHxauUBYOjQoary6enpyMnJUSvj4OCA8PBwVZn4+Hg4OjoiNDRUVSYyMhJSqRTHjh0z2fNSWrJkCZydnREcHIylS5eiurr6gc9JX+eliYSEBFRVVantp2vXrvD29tZqPw+qOedv6lJTU+Hp6YlOnTph7NixyMzMVL329NNPIy8vD7///jsSEhLQu3dvPProo8jPz290nwsXLoSrqytefPFFfYcvmunTpyM+Ph5btmzB33//jaeffhrDhg2r06zamMTERGRlZaG0tFSn38G///1vmJubo3fv3rC0tERGRgaGDBmi0WKmQUFB8PDwwJAhQ3D48GGNz4VMT6tKRhpb5behtsicnJxGyyv/21QZV1dXtddlMhmcnJx00gYq1nkBwGuvvYYtW7Zg//79mDZtGhYtWoTZs2c/8DkB+jkvTeTk5MDCwqLOAona7udBNef8TVl4eDg2btyInTt34vPPP0d6ejr69++P4uJiHDp0CMePH8e2bdsQGhqKLl26YNmyZXB0dMT333/f4D4PHTqEr776CuvWrTPgmRhWZmYmNmzYgG3btqF///7o3LkzZs2ahX79+mHDhg1Nvr99+/aQy+UIDQ3F6NGjsXnzZp19B5mZmbhy5QqKiorg6uqK3bt345NPPgEATJkypcGYPDw8sHbtWmzfvh3bt2+Hl5cXBg0ahJMnTzbvQyKjZ/TTwZNxi4mJUf1/r169YGFhgWnTpmHx4sWc6pm0EhUVpfr/Xr16ITw8HB06dMB3332H8vJylJSUwNnZWe09t2/fxsWLF5GZmYnu3burtr/zzjuYMWMGxo0bh3Xr1sHFxcVg52Fop0+fRk1NDfz8/NS2V1RUqD4vW1tb1fYXXngBa9euVT0/ePAgSkpKcPToUbz99tt4/PHH0atXr2Z9BwDw7LPPQiaT4Z133kFgYCAAoKqqClu3bsXWrVsB1K4hduLEiQZj8/f3h7+/v2p73759cfHiRaxYsQKbNm164M+MjE+rSkaas8qvu7t7o+WV/83NzYWHh4damaCgIFWZvLw8tX1UV1cjPz9fJ6sLi3Ve9QkPD0d1dTUuX76s9mPSHPo4L024u7ujsrISBQUFarUjhl4NurWvSu3o6Ag/Pz+kpaXB0dERHh4eiIuLq7eco6Oj2qgoJycnXLx4EZcvX0Z0dLRqu0KhAFBbM5mSkoLOnTvr+zT0rqSkBGZmZkhISICZmZnaa8p/6O/9bOzt7dXKdOzYEQDQs2dP5ObmYsGCBXjuuecAaP8dAMCKFSswZMgQODk54Y8//gBQ+7twbxLx559/4p///CcqKysbje1eYWFhOHToUOMfBpmsVtVM05xVfiMiItTKA8Aff/yhKt+xY0e4u7urlSkqKsKxY8dUZSIiIlBQUICEhARVmX379kGhUOhkDR2xzqs+SUlJkEqldZqlmkMf56WJkJAQmJubq+0nJSUFmZmZBl0NurWvSl1SUoKLFy/Cw8MDvXv3Rk5ODmQyGXx9fdUeLi4udbY7OTmha9euOH36NJKSklSPESNGYPDgwUhKSoKXl5fYp6gTwcHBqKmpQV5eXp3PRpm03rutsb9NhUKBiooK1XNtvwOgtolF+R0EBwcDqO2D1qlTJ1X5kpISeHh4wMLCQuPYkpKS1G6MqIURuwetoW3ZskWQy+XCxo0bheTkZGHq1KmCo6OjkJOTIwiCIIwbN054++23VeUPHz4syGQyYdmyZcK5c+eE+fPn1zsE1tHRUfjpp5+Ev//+Wxg5cmS9Q3uDg4OFY8eOCYcOHRK6dOmi86G9hj6vI0eOCCtWrBCSkpKEixcvCps3bxbatm0rjB8/3qjP6+bNm0JiYqKwY8cOAYCwZcsWITExUcjOzlaVeemllwRvb29h3759wokTJ4SIiAghIiJCZ+elqabOvyV58803hbi4OCE9PV04fPiwEBkZKbi4uAh5eXmCQqEQ+vXrJwQGBgq7du1SlXnnnXfURqk1xVRH0xQXFwuJiYlCYmKiAEBYvny5kJiYKGRkZAiCIAhjx44VfHx8hO3btwuXLl0Sjh07JixatEj49ddfG9zn6tWrhZ9//lm4cOGCcOHCBeHLL78UzM3NhRdeeEGr72DWrFnC5s2bG4xt9OjRgkQiEaKiooQ9e/YIH3/8sWBjYyOMGzeuwdhWrFghxMbGCqmpqcLp06eFmTNnClKpVNizZ49uP1gyGq0uGREEQfj0008Fb29vwcLCQggLCxOOHj2qem3gwIHChAkT1Mp/9913gp+fn2BhYSEEBAQIO3bsUHtdoVAIc+fOFdzc3AS5XC48+uijQkpKilqZmzdvCs8995xga2sr2NvbC5MmTRKKi4tN+rwSEhKE8PBwwcHBQbC0tBS6desmLFq0SCgvLzfq89qwYYMAoM5j/vz5qjK3b98WXnnlFaFNmzaCtbW1MHr0aLVkxZAaO/+WZMyYMYKHh4dgYWEhtGvXThgzZoyQlpamer2oqEiYMWOG4OnpKZibmwteXl7C2LFjhczMTI2PYarJiHKo6/0P5bVfWVkpzJs3T/Dx8RHMzc0FDw8PYfTo0cLff//d4D5XrVolBAQECNbW1oK9vb0QHBwshISEaP0dREZGNhnb5MmTBblcLgAQzMzMhO7du9c7v4/SRx99JHTu3FmwtLQUnJychEGDBgn79u3TxUdJRoqr9hIREZGoWlWfESIiIjI+TEaIiIhIVExGiIiISFRMRoiIiEhUTEaIiIhIVExGiIiISFRMRoiIiEhUTEZM2OXLlyGRSNTWdmjMxIkTMWrUKL3GpI0FCxaorXOj7/gWLFgAiUQCiUSClStX6u04ADBo0CDVsTT9fkh8Pj4+qu+toKCgxR3PWEgkEsTGxur9OD4+Pnr/W9eHiRMnqq4LQ3xOhj5efZiM6Ni9X6py3YWFCxeiurr6gfd7/z/UXl5eyM7ORo8ePR5o38bik08+wcaNG/V6jICAAGRnZ2Pq1Kl6Pc4PP/yA48eP6/UYupCTk4OZM2fC19cXlpaWcHNzw8MPP4zPP/8cZWVlejuuMllbsmRJndeGDx8OiUSCBQsWPHB5iUQCS0tL+Pn5YfHixdBkjseFCxciOzsbDg4OWp/X119/jX79+mlc/q+//sL27du1Po4xy8nJwYwZM9CpUyfI5XJ4eXkhOjq6zppRxkB5gzJs2LA6ry1duhQSiQSDBg164PISiQRmZmbw8vLC1KlTkZ+f32Rsw4YNQ3Z2ttpq1prKyMiAlZUVSkpKNCr/ySefIDs7W+vj6BKTET1QXkSpqal48803sWDBAixdurRZ+6qpqVGtNHo/MzMzuLu7qxaoEkNlZaXO9uXg4KC2Sq4+yGQyuLu7w9raWq/HcXJyQtu2bfV6jAd16dIlBAcHY/fu3Vi0aBESExMRHx+P2bNn49dff8WePXv0enwvL686yefVq1exd+/eehdE07b8lClTkJ2djZSUFMyZMwfz5s3D2rVrm4zLzs4O7u7ukEgkWp0PAPz0008YMWKExuXbtm0LJycnrY9jrC5fvoyQkBDs27cPS5cuxenTp7Fz504MHjwYr776qtjh1cvDwwP79+/HlStX1LavX78e3t7eD1xeeQOUmZmJDRs2YOfOnXj55ZebjEsul8Pd3R1yuVzLM6q9DgcPHqxatbkpDg4Ooq8EzmRED5QXUYcOHfDyyy8jMjISP//8MwBg+fLl6NmzJ2xsbODl5YVXXnlFLXvduHEjHB0d8fPPP6N79+6Qy+WYPHkyvv76a/z000+qLDsuLq7eZpqzZ8/iiSeegL29Pezs7NC/f39cvHix3jgVCgUWL16Mjh07wsrKCoGBgfj+++8bPTcfHx+8//77GD9+POzt7VU1DG+99Rb8/PxgbW2NTp06Ye7cuaiqqlJ775IlS+Dm5gY7Ozu8+OKLKC8vV3v9/tqf+qpYg4KCVHfAgiBgwYIF8Pb2hlwuh6enJ1577bVG469PQUEBpk2bBjc3N1haWqJHjx749ddfAdz9Pn799Vf4+/vD2toa//jHP1BWVoavv/4aPj4+aNOmDV577TXU1NRofWwxvfLKK5DJZDhx4gSeeeYZdOvWDZ06dcLIkSOxY8cOREdHq8o2dd1mZGQgOjoabdq0gY2NDQICAvDbb781evwnnngCN27cwOHDh1Xbvv76azz22GP1rt6qbXlra2vV3+GkSZPQq1cv1ZL22lq3bh28vLxgbW2N0aNHY/ny5XUS5/LycuzevVuVjGRnZ2P48OGwsrJCx44d8d///tdkmw009corr0AikeD48eN46qmn4Ofnh4CAAMTExODo0aMNvu/06dN45JFHYGVlBWdnZ0ydOlXt+ho0aBBef/11tfeMGjUKEydOVD3Py8tDdHS06vP+9ttvNYrZ1dUVjz32GL7++mvVtiNHjuDGjRsYPnz4A5dX3gC1a9cOkZGRePrpp5t9HR45cgRBQUGwtLREaGgoYmNj620Kvjcprq6uxmuvvQZHR0c4OzvjrbfewoQJE4yqyR5gMmIQVlZWqhoEqVSKVatW4ezZs/j666+xb98+zJ49W618WVkZPvroI3z55Zc4e/YsVq1ahWeeeUZV45KdnY2+ffvWOc7Vq1cxYMAAyOVy7Nu3DwkJCZg8eXKDTUSLFy/GN998g7Vr1+Ls2bN444038MILL+DPP/9s9HyWLVuGwMBAJCYmYu7cuQBq7yY3btyI5ORkfPLJJ1i3bh1WrFihes93332HBQsWYNGiRThx4gQ8PDzw2WefafU53m/79u1YsWIF/vOf/yA1NRWxsbHo2bOnVvtQKBSIiorC4cOHsXnzZiQnJ2PJkiUwMzNTlSkrK8OqVauwZcsW7Ny5E3FxcRg9ejR+++03/Pbbb9i0aRP+85//NJnIGZObN29i9+7dePXVV2FjY1NvmXtrBpq6bl999VVUVFTgwIEDOH36ND766KMm78osLCwwduxYbNiwQbVt48aNmDx5sk7KKwmCgIMHD+L8+fOwsLBotGx9Dh8+jJdeegkzZ85EUlIShgwZgg8//LBOub1796Jdu3bo2rUrAGD8+PG4du0a4uLisH37dnzxxRfIy8vT+vimIj8/Hzt37mzwmmqo1rO0tBRDhw5FmzZt8Ndff2Hbtm3Ys2cPpk+frtXxJ06ciKysLOzfvx/ff/89PvvsM40/78mTJ6vVuq1fvx5jx45t8HrRtrzS5cuXsWvXrmZdh0VFRYiOjkbPnj1x8uRJvP/++3jrrbfqlCsoKMChQ4dUychHH32Eb7/9Fhs2bMDhw4dRVFQkWr+QRom6TF8LdO+qoAqFQvjjjz8EuVwuzJo1q97y27ZtE5ydnVXPlSvKJiUlNbhfpfT0dAGAavXLOXPmCB07dhQqKyubjK28vFywtrYWjhw5olbmxRdfFJ577rkGz69Dhw7CqFGjGnxdaenSpUJISIjqeUREhPDKK6+olQkPDxcCAwPrjU95rBUrVqi9JzAwULW67scffyz4+fk1eL73mz9/vtrxBEEQdu3aJUil0jqrLCspv497Vy6dNm2aYG1trbbq8tChQ4Vp06apvff+78eYHD16VAAg/PDDD2rbnZ2dBRsbG8HGxkaYPXt2g++//7rt2bOnsGDBAo2PP3DgQGHmzJlCUlKSYGdnJ5SUlAh//vmn4OrqKlRVVal9z80tb25uLtjY2Ajm5uYCAMHS0lI4fPhwo3HVd82NGTNGGD58uNq2sWPHCg4ODmrbpkyZovo7P3funABA+Ouvv1Svp6amCgDq7F+5Iu+tW7cajc3YHTt2rN5rqj4AhB9//FEQBEH44osvhDZt2gglJSWq13fs2CFIpVIhJydHEIS73/+9Ro4cqVoZOCUlRQAgHD9+XPW68ju4//O+l/I3obKyUnB1dRX+/PNPoaSkRLCzsxNOnTolzJw5Uxg4cOADlZdKpYKNjY1gaWmpWtF4+fLljX4+9f3ef/7554Kzs7Nw+/Zt1bZ169bV+Y359ttvhdDQUNVzNzc3YenSparn1dXVgre3d72rV9/7vRiaeJ0NWrBff/0Vtra2qKqqgkKhwPPPP69qWtizZw8WL16M8+fPo6ioCNXV1SgvL0dZWZmqH4OFhQV69eql9XGTkpLQv39/mJubN1k2LS0NZWVlGDJkiNr2yspKBAcHN/re0NDQOtu2bt2KVatW4eLFiygpKUF1dTXs7e1Vr587dw4vvfSS2nsiIiKwf//+JmNtyNNPP42VK1eiU6dOGDZsGB5//HFER0dr1YcmKSkJ7du3h5+fX4NlrK2t0blzZ9VzNzc3+Pj4qN35u7m5tYi73uPHj0OhUGDs2LGoqKhQbW/qun3ttdfw8ssvY/fu3YiMjMRTTz2l0TUcGBiILl264Pvvv8f+/fsxbty4Rr8/bcqPHTsW7777Lm7duoX58+ejb9++9dYoNiUlJQWjR49W2xYWFqZqygNqa19++eUXfPfdd6r3yGQy9O7dW1XG19cXbdq00fr4pkJo5gLw586dQ2BgoFptysMPPwyFQoGUlBS4ublptA+ZTIaQkBDVtq5du2rcB83c3BwvvPACNmzYgEuXLsHPz6/R61eb8v7+/vj5559RXl6OzZs3IykpCTNmzNAornulpKSgV69esLS0VG0LCwurU+7eJprCwkLk5uaqlTMzM0NISEiDfRHFwmYaPRg8eDCSkpKQmpqK27dv4+uvv4aNjQ0uX76MJ554Ar169cL27duRkJCANWvWAFDvCGplZdWsznNWVlYal1W2x+7YsQNJSUmqR3JycpPNDfdXwcbHx2Ps2LF4/PHH8euvvyIxMRHvvvvuA3dulUqldX7g7u2H4uXlhZSUFHz22WewsrLCK6+8ggEDBtTpq9IYTT6z+5M7iURS7zZj++NujK+vLyQSCVJSUtS2d+rUCb6+vmqfiybX7T//+U9cunQJ48aNw+nTpxEaGopPP/1Uo1gmT56MNWvW4Pvvv2+yyUWb8g4ODvD19cVDDz2E7777DqtXr9Zbp9zjx4+jurq6WclOS9GlSxdIJBKcP39e5/tu6rdAFyZPnoxt27ZhzZo1Gl+HmpRXjqrs0aOHqgn4vffe02XoKpWVldi5c6dWnaiNBZMRPbCxsYGvry+8vb3V7toSEhKgUCjw8ccfo0+fPvDz88O1a9c02qeFhUWTHSR79eqFgwcPavRHquwcm5mZCV9fX7WHl5eXRjEpHTlyBB06dMC7776L0NBQdOnSBRkZGWplunXrhmPHjqlta6xDG1A70uDe4WZFRUVIT09XK2NlZYXo6GisWrUKcXFxiI+Px+nTpzWOvVevXrhy5QouXLig8XtaAmdnZwwZMgSrV69GaWlpo2U1vW69vLzw0ksv4YcffsCbb76JdevWaRTL888/j9OnT6NHjx7o3r27zssDgK2tLWbOnIlZs2ZpfQfv7++Pv/76S23b/c9/+uknDB8+XNXXyN/fH9XV1UhMTFSVSUtLw61bt7Q6tilxcnLC0KFDsWbNmnqvqYbmUenWrRtOnTql9p7Dhw9DKpXC398fQN3fgpqaGpw5c0b1vGvXrqiurkZCQoJqW0pKilZztwQEBCAgIABnzpzB888/r/PySv/617+wbNkyjX/7lfz9/XH69Gm1Gsv7r8O4uDi0adMGgYGBAGoTcjc3N7VyNTU1OHnypFbHNgQmIwbk6+uLqqoqfPrpp7h06RI2bdqk0VBDoHZkyd9//42UlBTcuHGj3oRj+vTpKCoqwrPPPosTJ04gNTUVmzZtqnP3C9R2OJ01axbeeOMNfP3117h48SJOnjyJTz/9VK2XuCa6dOmCzMxMbNmyBRcvXsSqVavw448/qpWZOXMm1q9fjw0bNuDChQuYP38+zp492+h+H3nkEWzatAkHDx7E6dOnMWHCBLWOpRs3bsRXX32FM2fO4NKlS9i8eTOsrKzQoUMHjWMfOHAgBgwYgKeeegp//PEH0tPT8fvvv2Pnzp1afQam6LPPPkN1dTVCQ0OxdetWnDt3DikpKdi8eTPOnz+v+qw1uW5ff/117Nq1C+np6Th58iT279+Pbt26aRRHmzZtkJ2drfE8FNqWV5o2bRouXLig9bweM2bMwG+//Ybly5cjNTUV//nPf/D777+r1V7+/PPPanejXbt2RWRkJKZOnYrjx48jMTERU6dObXatp6lYs2YNampqEBYWhu3btyM1NRXnzp3DqlWrEBERUe97xo4dC0tLS0yYMAFnzpzB/v37MWPGDIwbN07VRPPII49gx44d2LFjB86fP4+XX35ZLdHw9/fHsGHDMG3aNBw7dgwJCQn45z//qVVtMQDs27cP2dnZGjfvaFseqG2e7tWrFxYtWqRVbM8//zwUCgWmTp2Kc+fOYdeuXVi2bBmAu53N778Ogdrrd/Hixfjpp5+QkpKCmTNn4tatW0Z3HTIZMaDAwEAsX74cH330EXr06IFvv/0Wixcv1ui9U6ZMgb+/P0JDQ9G2bVu14Y1Kzs7O2LdvH0pKSjBw4ECEhIRg3bp1DfYhef/99zF37lwsXrwY3bp1w7Bhw7Bjxw507NhRq/MaMWIE3njjDUyfPh1BQUE4cuSIapSN0pgxYzB37lzMnj0bISEhyMjIaHKs/Zw5czBw4EA88cQTGD58OEaNGqXWd8PR0RHr1q3Dww8/jF69emHPnj345Zdf4OzsrFX827dvx0MPPYTnnnsO3bt3x+zZs01umG5zdO7cGYmJiYiMjMScOXMQGBioal6ZNWsW3n//fQCaXbc1NTV49dVXVdeRn5+fVqOlHB0dGxzVo4vyQO2d+/jx47FgwQKtmtQefvhhrF27FsuXL0dgYCB27tyJN954Q9V2f/HiRaSlpWHo0KFq7/vmm2/g5uaGAQMGYPTo0ZgyZQrs7OzU2vxbmk6dOuHkyZMYPHgw3nzzTfTo0QNDhgzB3r178fnnn9f7Hmtra+zatQv5+fl46KGH8I9//AOPPvooVq9erSozefJkTJgwAePHj8fAgQPRqVMnDB48WG0/GzZsgKenJwYOHIgnn3wSU6dOrXfId2NsbGy0Siy0La/0xhtv4Msvv0RWVpbG77G3t8cvv/yCpKQkBAUF4d1338W8efMAQHVN1ZeMvPXWW3juuecwfvx4REREwNbWFkOHDjW+61CUbrNEIqhvNI0+GfNoGqpffaNp6vPPf/5T6NevnyAItaO6oqKimnxPVlaWAEDYs2eP2vaWMpqGdKe+0TT12bx5s2Bubi6UlZUJCQkJgoODQ5OjC2tqagQ/Pz/hX//6V53XIOJoGtaMUKty+vRp2NraPvAcJ02JiopCQECAXo9B+vHWW2/B1tYWhYWFqm3Lli3DqVOnkJaWpmrKnDBhAgCgffv2mDNnTp397Nu3Dz///DPS09Nx5MgRPPvss/Dx8cGAAQNUZQICApo13Te1fMpRmfeO2vrmm29w6NAhpKenIzY2Fm+99RaeeeYZWFlZobq6Gp9++mmdmvCMjAysW7cOFy5cwOnTp/Hyyy8jPT1drZ/LSy+9pPFsrfoiuZMNEbV4+fn5qjUh2rZt26y1RzR19epV3L59GwDg7e3drEmOyPAyMjJU/bE6deoEqbT2fu2ZZ55BXFwciouL0alTJ8yYMaPOUPX77dq1C2+++SYuXboEOzs79O3bFytXrlTr09TQ8ah1y8vLQ1FREYDa6eeVTZL//ve/8dlnnyEnJwceHh4YNWoUPvzww0aXt8jKysKzzz6LM2fOQBAE1aiee5Piho5nSExGiIiISFRMw4mIiEhUTEaIiIhIVExGiIiISFRMRoiIiEhUTEaIiIhIVExGiIiISFRMRoiIiEhUTEaIiIhIVExGiIiISFT/D1fgXWAvg5d/AAAAAElFTkSuQmCC" }, "metadata": {}, "output_type": "display_data", "jetTransient": { "display_id": null } } ], "execution_count": 2 } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "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.13.2" } }, "nbformat": 4, "nbformat_minor": 5 }