{ "cells": [ { "cell_type": "markdown", "id": "832d4ebf", "metadata": {}, "source": [ "# 1. Getting Started with PA3Py\n", "\n", "In this notebook, we will learn how to initialize `PA3Py`, run a growth simulation for an embryo, and plot the results." ] }, { "cell_type": "code", "execution_count": 7, "id": "b21de4bb", "metadata": { "execution": { "iopub.execute_input": "2026-07-13T22:23:17.440816Z", "iopub.status.busy": "2026-07-13T22:23:17.440567Z", "iopub.status.idle": "2026-07-13T22:23:19.096545Z", "shell.execute_reply": "2026-07-13T22:23:19.094914Z" } }, "outputs": [], "source": [ "import sys\n", "import os\n", "import matplotlib.pyplot as plt\n", "\n", "# Add package path (if running from the local repository)\n", "sys.path.insert(0, os.path.abspath('../../src'))\n", "from pa3py import PA3Py\n", "from pa3py import constants as c" ] }, { "cell_type": "markdown", "id": "493e2e07", "metadata": {}, "source": [ "## Initializing the Simulation\n", "\n", "We load the simulation simply by passing the path of our HDF5 disk output to the main `PA3Py` class." ] }, { "cell_type": "code", "execution_count": 8, "id": "1782a077", "metadata": { "execution": { "iopub.execute_input": "2026-07-13T22:23:19.099347Z", "iopub.status.busy": "2026-07-13T22:23:19.098934Z", "iopub.status.idle": "2026-07-13T22:23:19.655057Z", "shell.execute_reply": "2026-07-13T22:23:19.654076Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[load_tripodpy_hdf5] Reading 100 snapshots from ../../tests/test_data/run_smooth_a0.001_v10...\n" ] } ], "source": [ "# Path to your TripodPy HDF5 data\n", "data_path = '../../tests/test_data/run_smooth_a0.001_v10'\n", "\n", "# Initialize the engine\n", "sim = PA3Py(data_path)" ] }, { "cell_type": "markdown", "id": "1620c608", "metadata": {}, "source": [ "## Running Growth\n", "\n", "We will tell `PA3Py` to compute the growth for an embryo located at 3 AU, with an initial seed mass of $10^{-3} M_\\oplus$." ] }, { "cell_type": "code", "execution_count": 9, "id": "978708aa", "metadata": { "execution": { "iopub.execute_input": "2026-07-13T22:23:19.658884Z", "iopub.status.busy": "2026-07-13T22:23:19.658554Z", "iopub.status.idle": "2026-07-13T22:23:19.668727Z", "shell.execute_reply": "2026-07-13T22:23:19.667697Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "-------------------------------------------------------------\n", " r [AU] M_tot [ME] M_iso [ME] f_silicates[%] f_H2O[%]\n", "-------------------------------------------------------------\n", " 3.00 9.016 9.02 50.0 50.0\n", "-------------------------------------------------------------\n", "\n" ] } ], "source": [ "# Embryos to simulate\n", "embryos = [3.0]\n", "# Execute the module\n", "results = sim.run_growth(embryos, m_seed_me=1e-3)" ] }, { "cell_type": "markdown", "id": "af1d23df", "metadata": {}, "source": [ "## Plotting Mass Evolution\n", "\n", "The `results` dictionary contains arrays where the first column is time (in seconds) and the second is the total mass of the embryo." ] }, { "cell_type": "code", "execution_count": 10, "id": "9dc7ee84", "metadata": { "execution": { "iopub.execute_input": "2026-07-13T22:23:19.671281Z", "iopub.status.busy": "2026-07-13T22:23:19.671052Z", "iopub.status.idle": "2026-07-13T22:23:20.034296Z", "shell.execute_reply": "2026-07-13T22:23:20.033439Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAqgAAAHbCAYAAAAH5miuAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAAU5RJREFUeJzt3Qd4U2X7x/G7e7e0TAEBB1NQpIAICoqi4kZU3OBW/LvXiwvEvVBRcIugoujr5AVxoiAOpjIVkL3L7G5pyf+6n5KYlq60aXLOyffDda6cJCfJk3OS8OuzTpiIuAQAAACwiPBgFwAAAADwRkAFAACApRBQAQAAYCkEVAAAAFgKARUAAACWQkAFAACApRBQAQAAYCkEVAAAAFgKARUAAACWQkAFYAmrV6+WyZMnB7sYIWPcuHGSlZUV7GIAQLkIqAAqNHjwYHG5XBUuxxxzTMjuvaSkJLnvvvtkzpw5snv3bsnPz5c1a9bIhx9+KKeffrpYQVxcnAwfPlz69OkTtDJcfPHFcuutt1Z7+2HDhsmvv/4q27Ztk7y8PFm+fLk8//zz0qBBA59eNyUlxTxeP6ft2rUrd5vp06fLokWLyr2vfv365rG6/wAEXmQQXhOAzTz44IOmhrOslStXSig67LDD5Ouvv5aWLVvKZ599JhMmTJDs7Gw5+OCDTTidMmWKXH755fLee+8FtZzx8fEyYsQIs/z0009BKcMll1wiHTt2lBdffLFa26enp8sff/xhgr7W8LZv316uvfZaOeOMM6Rz586Sm5tbree54IILTMDcvHmzXHrppeYzDMA+CKgAqvTVV1/JvHnzLLunNIhVN7jUVkREhAmljRs3NjWTv/zyS6n7R44cKf369TPbWaXMdnL++ecfcJvWqH7yySdy1llnyaRJk6r1PJdddplMnTpV1q5da0IyARWwF5r4AdSa1iRqbdWdd94pQ4cOlX/++UdycnJMLWPz5s3NNg888ICsX7/ehLLPP/9cUlNTy30uDXcLFiwwzbNLliyRAQMGlNvtoHfv3jJmzBjZunWrbNiwQU444QRz+7nnnltuM7Pe16NHD89tJ554osyYMcPUfO7atcuUqaKm4LI1c506dZJHHnnkgHDq9u2338q0adOqLLPbjTfeKIsXLzbdBDZu3Cgvv/yyaaJ2u/nmm6WoqKjUbXfccYd5zueee85zW3h4uGRmZsqTTz5pjsn27dvN7VqD6u6WUbbJumnTpiZwa22lNqs/88wz5nmqcvbZZ8v//vc/U14tt9am6zH2fqw2oZ955pnSqlUrz+uXVxNfFe06oerVq1et7bUm+/jjjze1sLoceuihcuyxx/r8ugCChxpUAFXSYKR98rxp2Ni5c2ep27QpNTo6Wl566SVJS0uTe+65Rz766CP54YcfTIB86qmn5PDDDzeB69lnn5Wrr7661ONbt25tasheffVVGT9+vFx55ZXy8ccfy2mnnSbfffddqW3Hjh0rGRkZpsYyISFBfvzxR1m3bp0pg4bNsuXSAPXbb7+Z6yeddJKpFV61apUJb9pXU8s0a9Ys6dKli6l1q4jW4qmaNN+XLbPSwKhl0FD7yiuvSNu2bU1g7datm/Tq1csE05kzZ5oa2eOOO850H1AawIqLi82l29FHH236xmrw1te54YYbzL789NNPzaIWLlzo2V6fU/+I+P333+Wuu+6Sk08+2VzqHxj6uMoMGTLEhPtRo0aZy759+5rQnpycbI67euyxx8xnR/9Iuf32281tum116OctMjLSfCY0cOt+0GNcHfoHif6BpAHaHZ71M6A1sQDsw8XCPuAzwGegvM/A4MGDXRXJy8vzbNeyZUtz29atW13Jycme2x977DFz+4IFC1wRERGe299//31Xfn6+Kzo62nPb6tWrzbYDBgzw3JaUlOTauHGja968eQeUacaMGa7w8PBS5dXX03J5l6FBgwauwsJC1/Dhwz23zZ8/37VlyxZXamqq57ZOnTq5ioqKXO+8806l3wcty86dOw+4PT4+3lW/fn3PomWvqsxaNt0P06ZNc4WFhXluHzp0qNl+yJAh5rret3v3bteTTz7p2SYjI8M1adIk1969e10JCQnmtttuu828h5SUFHNdy6G837t7GTdunLnvgQceOOD9zZkzp8rfhNjY2ANue+WVV1zZ2dmljuvkyZPNsfXlN6Zx48alPmvr1q1zXXDBBdV+/J9//ul69913PdcfffRR17Zt20p9BnWZPn26a9GiReU+R2X7joV9wGdA6nwf0MQPoErabK+1a95L//79D9hOazu1idlNa+bctY1a2+d9e0xMjDRr1qzU47W5WJub3bTZWQcgaa2m9vn09sYbb8i+fftK3abbxsbGlurHOGjQIImKivLUeDZp0sTUNL7zzjumad9NR3NrLWZVI/C1hrC8WkCtLdQmdfcyceLEA7YpW2bdj7ofXnjhBVMj7b3dnj17zMAgpfdpdwLtIqB04JCOateaxbCwME/ztdamalcBfWx1la0p1dpabRKvitZMuiUmJpoaT32s1gxXp6tEZbRmXveNdg/QvqO6P/U1qkO7Xxx55JHywQcfeG7T9YYNG8qpp55aq3IBCBya+AFUafbs2dUaJKVN7N7cQUn7npZ3u/ZD9e6TWN6sADrNkNJ+jNp30628vox///23Kas257799tvmNnfTrjZbK+2b6d62rGXLlpnuBJUNYNLQXLa7g7v5XpuUK2v+L1vmisqyd+9e0/3Afb/S8KddATSAaxDdtGmT6av7559/muvaBUK7AGiXiurSfr7ufqpuGtq1e0ZVOnToII8++qhp2vfuG6vKXveVvv/vv//erGuXBl3XgK59ZN1dHCobHKV/QOj+09kW3GFa971+FnTglC+8/3AAEDgEVAB+411LWp3btfavpjRclUdrUXVKI62d1dpJrV286aabxF/++usvUwOrg4s0JLqtWLHCLGVrF6tT5ur4+eefTf9efT8aSDWwKr3U69p3tVGjRp7bq6Oi41IVDaA6bZXWlj/00EMm/Ot71prup59+ulqDrHyhf2DovtaAWVVA1f6nWtuqf2yUpftHa3i1f6rSMmv/4/LoHynubQAEHk38ACxDB1CV1aZNm1Ijuauio7Y1eGlQ0UBTWFhYamoi9wAoDXRladO0Di6qbPondy2pPndtVVQW7ZJwyCGHlBqspTXDBQUFJox6B1QdEKUnTNCBX+7rdV37pwPetIuBDpQaPXq0p5bTu8uEv8ugNcdV1czqtF86gl+7BWg3D+9F51LVcOo9y4PuX91en7ss9zGpbMAcgLpDQAVgGVrr6T2tlI5Iv+KKK0xTtnfzfmV27NhhRuhrU6+GSJ3uSW9z27Jli3k+nfrJO/AcccQRcsopp1TZBKxN6Dr9lYagis6kVd2aYW2W19B5yy23lLpdZzfQKZW8awt1Oz1rlQZvbfr3rkHV2j59Du0ioe/PzR20qzs9k681r97vU0O19lUuS2srq9vkr++jvBrN8847z3Q7mDt3brWa93WqLJ031Xt58803TXcR7z8s9FhrrfT1119f6nn0felMCrrP3V0NAAQWTfwAqqQDosob+KL9Amsyr2VFtC/mW2+9ZaZY0kB61VVXmcFROt2UL7SZX0OJKm+C9rvvvtuEWG061tdzTzOlfWO1n2dldLojDdE6PZM2u+v0TRoSNYhpwNb5QTVAVtUUrbT/5xNPPGFeU4P0l19+aWruNOhpjWnZvqz6OnoqUD21qvsUnVrjq90O9PiMGzeu1PbaPK1hWgeKaTjTwUc6iEpvqw097vpcOhWY1qBqLameOau8YK59ly+66CIzX6sGbA2Q7lrosnRKKQ3tWuOt70kHlHXt2tUET/2cVXY2Kg2aAwcONAPdNFiWR/evnnZVB0zpfps8ebI5jjpIrXv37uZ9aUjWY6j9ee+///4D+ugCCBymS2Af8BngM+DzNFNK7/eeZurOO+8s9fg+ffqY2wcOHFju86anp3tu06mIdEqifv36uf744w8zXdTSpUur9diyS1RUlGvHjh2uXbt2uWJiYsrdpm/fvq6ZM2e6cnJyzBROX3zxhatdu3bV/i7oVFY6RZNOy5SZmWmmi1q7dq3ro48+cp1xxhk+lVmnldL3WlBQ4Nq8ebNrzJgxnqmivJf+/fub55kyZUqp219//XVz+5VXXnnAY3r06GGmjdLyeU+bpNNMZWVlHbC93q+qev/HHnus65dffjH7b8OGDWYKLD12So+79/Rb7733npmaS1U25ZRO7fTqq6+afaFl0zL//fffrlGjRpn7KiuPTk9W0T5wL7179zbb3HzzzZ7bdEqshx56yLymfub0dfV9XXLJJfwu8rvIZ0CCtw/0z13/dBACAIvQCeh1UI3WkF1zzTXBLg4AwEf0QQXgODoQRkdsa1M/AMB+qEEF4Bjaj1AnaXdP7p6enh7sIgEAaoAaVACOoSOv9Xz2OqG7jv4HANgTNagAAACwFGpQAQAAYCkEVAAAAFiKoybq13NjZ2VlBbsYAAAAqICeJVCnAgyJgKrhdOPGjcEuBgAAAKqgZ96rLKQ6JqC6a071DQeiFlXTvwbiQL0e/I9jaG8cP/vjGNofx9D+kgKcZ9yvV9VrOSaguukbDmRgDPTrwf84hvbG8bM/jqH9cQztL8tieYZBUgAAALAUAioAAAAshYAKAAAASyGgAgAAwFIIqAAAALAUAioAAAAshYAKAAAASyGgAgAAwFIIqAAAALAUAioAAAAshYAKAAAASyGgAgAAwFIIqAAAALCUyGAXwM4KCwslMjJSoqKiDrjP5XJJUVGR53p529T1tmrv3r012lbfV1hYmG22jYiIkPDwcJ+21ecv7xjq/tX9XJ3nDcS2up1uX5Hi4mLZt2+fZbbVY6b7tCK6nW5f223LHr+aPm9V3w1ftrXC995OvxHuY1j2Nn//RnTs21t6XTSw8venz+N+Lv3u7f/+sW3l+yEiMlI+XLVErhr7nBR7fe7YvxXuCJFwr8/s/t/NYG4bERZmjmF0XKxIVpZYhZawkk+ffSQlJUlmZqYkJydLVgB2sL7enXfeWeH9y5cvl4kTJ3qu33fffRIdHV3utmvWrJF33nnHc/3uu++WhISEcrfduHGjvPHGG57rt912m9SrV6/cbbdt2yZjx471XB86dKg0atSo3G13794tL7zwguf6tddeK82aNSt325ycHHnmmWc814cMGSKtWrUqd1v9z+fxxx/3XL/kkkukTZs2UpERI0Z41i+44AI54ogjKtz2scce8/wHdO6550rnzp0r3Pbpp5+W3Nxcs3766adL9+7dK9xW94PuD9WvXz/p1atXhduOGTNGMjIyzPoJJ5xgloq8/vrrsmnTJrPes2dPOeWUUyrcVj8P+rlQ3bp1kzPOOKPCbd9//31ZsWKFWdd9oPuiIh999JEsXbrUrHfo0EEuvPDCCrf9/PPP5Y8//jDrrVu3lksvvbTCbadMmSJz5swx6/pZ0M9ERb755hv55ZdfzHrTpk3luuuuq3DbH3/80SyqYcOGctNNN1W47axZs+Tbb7816/qd0O9GRWbPni1Tp0416/Hx8XLPPfdUuK3uA90XSgPO/fffX+G2S5YskY8//rjcz3NZ/EaIZ58++eSTnt9tf/9GaFgfOXOahK1YK0VrS75/5Uk480QJjyn5jc5fsFT2rlpf8ban9ZbwhLiSbRf+LXtXlHxXyxPfr5dEJCea9YKlK6Vw2T8Vb3tiD4lISzHrhX+vloLFyyvcNq53N4lsmFay7T/rpOCPZRVv27OLRB7U0KzvXbNR8uctrnDb2GOOkqjmTUq23bBF8n//s+Jt0ztKVKuS/yeKNmdI3i/zK9w2pnN7iT6sRcm2GTslb8acirft2Eai2x5i1ot37pHc6b9VuG10+8MkpsPhJdtmZkvut7Mq3DaqdSuJPbKtWd+Xkyc502ZUvO2hB0vs0R1Kti0olJz/Ta9w28iWTSWuayez7ioqkuwvvq9422aNJa7Hv/9XZX3ydYXbRjRpIPG90v/d9vPvtDag/G0bpEp8n3//X8ue/IO4Cv/9Y85beGqyJPQ99t9tv/pJEvv3kSf6nSvbt2wVq+Q1mvgBAI4VXy9F4lOSg10MAD6iBrUWfwFs375dmjRpItnZ2QfcT/Od9Zv4ExMTZcuWLQccQ5r47dHEX/b40cRvvyZ+9zFs0KCBpybF378RB7U5TO765D1xFe+TuV9OlU8fe9avvydVfTd82daXrj2B2Laq76d+N7W1b9euXZKWliZ5eXl++Y2o6bZVfYZ92bauugFZsXtRamqqaQmsl5ISsBbo6tSg0ge1FrTJXg+y9w9QRaqzjZW29f7w2mFb/RJ7f5Grs60+f1XHsCbP6+9t9UfK+0fV6tvqj191P2u12bay4+fL8yq29X0/+OO77D6G/n5eb8kNGpjLsIhw2Z2RIYUFBVU/bzU/66pYQ0Q1y2y3bTWmFlbxO+Xat0/Cw8JkX3FxtfZtdZ+3Jtuq6paBbf+1t7DQHEOroYkfAOBYyQ3re9azMnYEtSwAqo+ACgBwrOSGJYODVGbG9qCWBUD1EVABACFRg7qHgArYBgEVAOBYSQ1o4gfsiIAKAHCs5IYlg6RU5nb6oAJ2QUAFADg+oObs3iPFPsxSACC4CKgAAMf3QWWAFGAvBFQAgCPFJSdJVEyMWSegAvZCQAUAOFKy1wCpTOZABWyFgAoAcKTkRv/OgZq1nTlQATshoAIAHF+DumcbARWwEwIqAMDxk/QzxRRgLwRUAIDjT3OaxVmkAFshoAIAHInTnAL2RUAFADhSkncTP6P4AVshoAIAHCm5QclZpPIys6SooCDYxQHgAwIqAMDRpzllgBRgPwRUAIDjxCTES0x8nFnPZIopwHYIqAAAx9aeqkwm6Qdsh4AKAHAcTnMK2BsBFQDg6NOcUoMK2A8BFQDg7BpU+qACtkNABQA4vA/qjqCWBYDvCKgAAEefRSqT05wCtkNABQA4ugY1i7NIAbZDQAUAODag5ufkSEFubrCLA8BHBFQAgOMk7R8kRe0pYE8EVACAo0THxUpcUqJZZ4AUYE8EVACAoyQ18BrBvy0jqGUBUDMEVACAc0fwM8UUYEsEVACAc+dAZQQ/YEsEVACAgyfp3x7UsgCwcUANDw+XkSNHyqpVqyQ3N1dWrlwpDzzwQLCLBQCwexM/pzkFbClSLODee++VG2+8UQYPHixLliyRrl27yrhx42TPnj3y0ksvBbt4AAAbSfYeJMVZpABbskRA7dmzp3zxxRcydepUc33t2rVy8cUXS/fu3YNdNACAzTBICrA/SwTUX375Ra677jpp3bq1rFixQo488kg57rjj5I477qjwMdHR0RITE+O5npSUVOqyrgX69eB/HEN74/jZX10dw3qNG5nLvQUFEiVhEsXvdJ3he2h/SUHKT1UJExFXnZemqkKEhcnjjz8u99xzjxQXF0tERITcf//98uSTT1b4mOHDh8uIESMCWk4AgPWNWTZHCoqLJSUqRq5ue3SwiwOgHMnJyZKVlSWWrkG98MIL5dJLL5VLLrnE9EHt3LmzvPDCC7Jp0yaZMGFCuY954oknZNSoUaUS+caNG6VZs2aVvmF/CfTrwf84hvbG8bO/ujiGkdFR8sCPU8z6ornzJLlbH788L8rH99D+koKUn6piiYD6zDPPmNrSSZMmmeuLFy+Wli1byrBhwyoMqIWFhWYpS3duIANjoF8P/scxtDeOn/358ximNm3iWd+1dRu/zwHC99D+siyWZywxzVR8fLzs27ev1G3a1K/TTwEAUF0pDRt61vdwmlPAtixRgzp58mTT53TdunWmif/oo482A6TefvvtYBcNAGAjSV5zoGZxmlPAtiwRUG+++WZ55JFHZOzYsdKoUSPT9/S1114zk/cDAFBdnOYUcAZLBNTs7Gy5/fbbzQIAQE0lN/A6ixST9AO2RSdPAIBjJDfiLFKAExBQAQCOwWlOAWcgoAIAHHea06LCQsndkxns4gCoIQIqAMBxg6QyGcEP2BoBFQDgCBGRkZKYlmrWGSAF2BsBFQDgCEn10zzrmRk7gloWALVDQAUAOELS/uZ9xST9gL0RUAEAjpDiNcUUpzkF7I2ACgBwhCSvSfqzaOIHbI2ACgBw3mlOt28PalkA1A4BFQDgwNOcMkgKsDMCKgDAERLrl0wxpbJ27AxqWQDUDgEVAOAIsQkJnvW8rOyglgVA7RBQAQCOEJMQby6Li4qkqKAg2MUBUAsEVACAo2pQC3Jyg10UALVEQAUAOEJMYklAzc/OCXZRANQSARUA4Kga1PwcAipgdwRUAIDthUdGSHRcrFkvoAYVsD0CKgDA9mLi/x3Bn59LH1TA7gioAADbi90/gl9RgwrYHwEVAOCYAVKKQVKA/RFQAQCOmqSfQVKA/RFQAQC2F5tIEz/gJARUAICzalAZJAXYHgEVAOCY05wqBkkB9kdABQA4a5AUpzoFbI+ACgBwVhM/E/UDtkdABQDYXqxXDSpN/ID9EVABALbHICnAWQioAABH9UGlBhWwPwIqAMBRo/iZqB+wPwIqAMBRTfwFjOIHbI+ACgBwzCCpgtw82VdcHOziAKglAioAwDFN/AU5OcEuCgA/IKACAJxTg0rzPuAIBFQAgO3FxJfUoDJACnAGAioAwNaiYmMkIjLSrHMWKcAZCKgAAAeN4KcPKuAEBFQAgGMm6c/Pzg1qWQD4BwEVAOCcGtRcAirgBARUAIBzziKVTRM/4AQEVACAI6aYUgRUwBkIqAAAW2OQFOA8BFQAgK1Rgwo4DwEVAGBrMQySAhyHgAoAsDVqUAHnIaACAGyNUfyA8xBQAQCOGSSVz5mkAEcgoAIAbC028d95UAuYBxVwBAIqAMDWGCQFOA8BFQDgiEFS+4qLpTAvP9jFAeAHBFQAgK3FxJc08dP/FHAOAioAwBE1qJzmFHAOAioAwBF9UAtycoNdFAB+QkAFANhWeESExMTHmXVqUAHnIKACABwxST81qIBzEFABALbFJP2AMxFQAQDOqEFlkn7AMQioAADbogYVcCYCKgDAtmL2TzGlGCQFOAcBFQBg+zlQFYOkAOcgoAIAbCvWqw8qZ5ICnIOACgCw/ST9ikFSgHMQUAEAtkUNKuBMBFQAgEMGSXGqU8ApCKgAAEdMM1WQmxPUsgBwYEBt2rSpvPvuu7J9+3bJzc2VhQsXSnp6erCLBQCwySh+ppkCnCNSLKBevXoya9YsmT59uvTv318yMjKkdevWsmvXrmAXDQBgkzNJ0cQPOIclAuq9994r69evl6uuuspz25o1a4JaJgCA9dHEDziTJZr4zz77bJk7d6589NFHsnXrVpk/f75cc801wS4WAMAmg6QK8/JlX1FxsIsDwEk1qIceeqjceOONMmrUKHn88celW7duMnr0aCksLJQJEyaU+5jo6GiJiYnxXE9KSip1WdcC/XrwP46hvXH87M8fxzA+KdFcFubl8XscBHwP7S8pSPmpKmEi4pIgKygoMDWovXr18tz24osvmqDas2fPch8zfPhwGTFiRABLCQCwmjHL5khBcbHUi46Vq9p0DnZxAFRTcnKyZGVlWbsGdfPmzbJ06dJSty1btkwGDhxY4WOeeOIJU+Pqncg3btwozZo1q/QN+0ugXw/+xzG0N46f/fnjGD404ysJj4yQpQsXSXLX3n4vIyrH99D+koKUn6piiYCqI/jbtm1b6rY2bdrI2rVrK3yMNv/rUpbu3EAGxkC/HvyPY2hvHL/QPYaRMTEmnKrczEx+i4OI76H9ZVksz1hikNTzzz8vPXr0kGHDhslhhx0mF198sVx33XUyZsyYYBcNAGBRsYleU0zlMEk/4CSWCKja/3TAgAEmmC5evFgefPBBue2222TixInBLhoAwA5TTOVwmlPASSzRxK+mTJliFgAAqoOzSAHOZYkaVAAAfBVTqgaVJn7ASQioAABbiuU0p4BjEVABALY+i5RikBTgLARUAIDtB0nlZ9PEDzgJARUAYPtBUoziB5yFgAoAsCUGSQHORUAFANgS00wBzkVABQDYUkw8Z5ICnIqACgCwJWpQAecioAIAbIlTnQLORUAFANhSTGJJE/++ffukMC8v2MUB4EcEVACArWtQC5gDFXAcAioAwJZi9p/qlLNIAc5DQAUA2HqQFGeRApyHgAoAsJ2w8HDPNFMFubnBLg4APyOgAgBs27yv6IMKOE+kLxufddZZPr/At99+K/n5+T4/DgCA6kwxlZ9DDSoQ0gH1888/9+nJXS6XtG7dWlavXu1ruQAAqFYNKn1QgRAPqKpJkyaSkZFRrW0zMzNrUiYAAKp/FqmcHPYWEMp9UMePHy95PkyG/N577xFSAQB+x1mkAGfzqQb1qquu8unJhw4d6mt5AACoUoxXDSqDpADnYRQ/AMDmg6Ro4gechoAKALB3H1ROdQo4DgEVAGDvUfzUoAKO49eAGh8fL507d5b69ev782kBAKh4kBQ1qIDj+DzNVEWGDRsmJ5xwgmRlZUlxcbFkZ2fLTTfdxCT9AAC/i0n0rkFlon7AafxSg3rllVdKamqqnHrqqWZi/kGDBslXX30lL7zwgj+eHgCASqaZYpAU4DR+CagaSJcsWSIPPfSQNGrUSB588EFJS0uTbt26+ePpAQAohUFSgLP5pYm/sLBQtmzZItu3bzdN+nPnzjVN/WFhYf54egAASomJ/7eJv4AmfsBx/FKD+uabb8pll11mmvU1mOpl165dZdq0af54egAAyq1B3ZtfIMVFRewdwGH8UoP65ZdfSkJCgnz//feyYMEC+frrr+WXX36RBx54wB9PDwBAKTH7+6AyxRTgTH4bxf/BBx+YRQdL7dq1y19PCwDAAdw1qDTvA87k94n6CacAgECN4ucsUoAzcSYpAICtRMbESERUSQMgTfyAM/kcUO+//37p37+/mU4KAIBAi/U6zSlnkQKcyec+qCNHjhSXy2XWdWqp+fPny7x58zyXmzZtqotyAgBQaoCUKsjlLFKAE/kcUOfMmSMHHXSQjBs3zsx72qVLFznvvPPkvvvuk4iICMnIyDBh9YwzzqibEgMAQlqs92lOszmLFOBEPgfUHj16yODBg+Xxxx83YfWOO+6QVatWSXR0tHTu3NkE1qOPPrpuSgsACHnepzmlDyrgTDUaJDV+/Hhp06aN/P3336a2VMOq1p7Onj1bXn31Vbn++uv9X1IAADjNKRASajyKPycnR+69915zxqiOHTvKypUr5fLLL/dv6QAAKCNm/xyoinlQAWeq1TRTWmsaExNjJujfsGGD6ZeqE/UDAFBXElJSPOu5mZnsaMCBfO6DqrWmnTp1Mku7du0kPz9fFi5caJr3X3vtNdmzZ0/dlBQAAA2oafU8+yFn5272CeBAPgdU7W+6Zs0a0w9Va05XrFhRNyUDAKAciV4tddmcWhtwJJ+b+GfOnCn169eX4cOHm3lPf/75Zxk9erQMGTJEjjzySAkP5+RUAIC6k5BKDSrgdD7XoJ5wwgnm8vDDD5f09HQzrZQul156qdSrV08KCgpk0aJFcswxx9RFeQEAIS7RK6Bm76KJH3AinwOqm47a12XSpEme21q1amVG9TMPKgCgriSmpXrmQC0qLGRHA6EeUHVg1OLFiz2nOi1L+6bq8t///tdc79Chg5krtbi42D+lBQCEvIR6JaP4c6g9BRzLpw6jCxYsMP1Pq+vXX3+VFi1a1KRcAAAcICw8XOJSks16NiP4AcfyqQY1LCxMHnnkEcnNza3W9nr6UwAA/Fl76h6MSw0q4Fw+BdQZM2ZI27ZtfapBzcvLq0m5AACofAT/bgZIAU7lU0A98cQT664kAAD4MoKfJn7AsZi0FABgGwn7R/CrHCbpBxyLgAoAsA3mQAVCAwEVAGAbNPEDoYGACgCwDZr4gdBAQAUA2AZN/EBoqHFAjY2Nlbi4OM91nZD/1ltvlX79+vmrbAAAlJJQz2uaKUbxA45V44D6xRdfyBVXXGHWU1JS5Pfff5c777zT3H7DDTf4s4wAABgJaSUBdW9BgRRU86QxAEIooHbp0kVmzpxp1s8//3zZunWrtGzZ0oTWW265xZ9lBACgVBM/Z5ECnK3GATU+Pl6ysrLM+imnnCKffvqpuFwu+e2330xQBQCgrpr4maQfcLYaB9SVK1fKueeeK82bN5dTTz1VvvnmG3N7o0aNJDMz059lBABA4pKTJCKq5ASInOYUcLYaB9SRI0fKs88+K2vWrDH9T7Xm1F2bumDBAn+WEQAASfA6zSlN/ICzlfwpWgOffPKJGbl/0EEHyZ9//um5/fvvv5fPPvvMX+UDAMBITP33NKc08QPOFlmbaaa0KV8HRykNqwMGDJBly5bJnDlz/FlGAAAkcf8IfpW9axd7BHAwv08z9fnnnzPNFADA72jiB0IH00wBAGyBJn4gdDDNFADAFhJSUzzrOTTxA45myWmm7r33XjOn6vPPP1+r5wEAOLOJP3vX7qCWBUCITTPVtWtXuf7660vNDAAAAE38QOgIr+00UxooTzvttFLTTN1+++01es6EhAR5//335dprr5VdNN8AALz/j9g/ir+4qEjy95/JEIAz1XiaKaVTTLmnmXKrzRRTY8aMkSlTppiQ+8ADD1S6bXR0tMTExHiuJyUllbqsa4F+Pfgfx9DeOH6hdwyT09LMZd6eTElMTKzTsqF6+B7aX1KQ8lOdBlTVvn17U5OqgdHb5MmTfXqeQYMGSZcuXaRbt27V2n7YsGEyYsSIA27fuHGjBFKgXw/+xzG0N45faBxDHZcweulsKXa5pGXTZpxS22L4HtrfRovlmRoH1EMOOcScMapTp07mhyMsLMzcruvmiSOr/9Q60OrFF1+Ufv36SUFBQbUe88QTT8ioUaNKJXLduc2aNZOsADT9BPr14H8cQ3vj+IXWMYyOi5X7vv/SrM/95Ve5Of34AJUSleF7aH9JQcpPdRZQNVCuXr1aTjrpJHPZvXt3qV+/vjz33HNy1113+fRc6enp0rhxY5k/f/6/BYuMlN69e8v//d//mab8ffv2lXpMYWGhWcrSnRvIwBjo14P/cQztjeMXGscwLfnfJv0923fwu2sxfA/tL8tieabGAfXYY4+Vvn37yo4dO0x41GXWrFmm6X306NGmub66tM9px44dS902btw4+euvv+Spp546IJwCAEJLQmqqZz2HKaYAx6txQI2IiPAk7e3bt0vTpk1l+fLlsnbtWmnbtq1Pz5WdnS1LliwpdVtOTo4Jv2VvBwCEnkTvOVB37gpqWQBYOKAuXrxYjjrqKM88qPfcc49pcr/uuutk1apV/i0lACCkeU/STw0q4Hw1DqiPPvqombdUDR8+3Izanzlzpqn1vOiii2pdsBNPPLHWzwEAcGANKk38gOPVOKC6T23qPu2pTjeVmprKBPsAgDqbpF/RxA84n88B9a233qrWdldffXVNygMAQKWnOaWJH3A+nwPqkCFDzECoBQsWeOY+BQAgUH1QszkVNuB4PgfUV155RS6++GIzUb9OBfXee+/RrA8ACFgf1NzdmextwOHCfX2ATpx/0EEHydNPPy1nnXWWrF+/XiZNmiSnnHJK3ZQQABDy3H1Qc3bvkX3FxSG/PwCn8zmgKp1O6sMPPzShtEOHDmau0rFjx5opp9wj+wEA8HcfVPqfAqGhRgHVm57lyeVymf6oOnk/AAD+FBEVJbGJJZUfBFQgNNQooEZHR5u5TnWqKT17VKdOnUzTf4sWLcwZoAAAqJsBUrvZsUAI8HmQ1JgxY0w41b6nb7/9thkwpZPzAwBQ95P0c5pTIBT4HFBvuOEGWbdunTmdaZ8+fcxSnoEDB/qjfACAEJfoNUl/zk5qUIFQ4HNAnTBhgulzCgBAICR4TdJPEz8QGnwOqFdeeWXdlAQAgCqa+HNo4gdCQq1H8QMAEIg5UFU2TfxASCCgAgBsMQeqYpopIDQQUAEANppmilH8QCggoAIAbNQHdU9QywIgMAioAABLS0wraeLPz8mRosLCYBcHQAAQUAEAlpZQL8Vc0v8UCB0EVACAZYWFh0tcSrJZZwQ/EDoIqAAAS9eehoeX/FdFDSoQOgioAADLYgQ/EJoIqAAAe4zgZ5J+IGQQUAEAlpWwfwS/yt61O6hlARA4BFQAgE3mQGWSfiBUEFABALYIqIziB0IHARUAYFk08QOhiYAKALAsmviB0ERABQBYVkI97z6oe4JaFgCBQ0AFAFhWQlpJQN1bUCAFubnBLg6AACGgAgAs38TPWaSA0EJABQBYvomfEfxAaCGgAgAsKT4lWSKiIs06c6ACoYWACgCwpLRmB3nWd27eEtSyAAgsAioAwJLSmjX1rO/aSEAFQgkBFQBgSWlNvWpQN24KalkABBYBFQBgSWnN/61B3UFABUIKARUAYPk+qLs2bg5qWQAEFgEVAGBJqfub+Avz8iVrx85gFwdAABFQAQCW7oO6cxO1p0CoIaACACwnqX6aRMfFmvVdBFQg5BBQAQCWk+o9Byr9T4GQQ0AFAFhOfa85UHduYIopINQQUAEAlp6knz6oQOghoAIALCe1WRPPOpP0A6GHgAoAsHQT/44NjOIHQg0BFQBg2Sb+/OwcycvMDHZxAAQYARUAYClh4eFS76DGZp3mfSA0EVABAJaS0rCBREZFmXUCKhCaCKgAAEtJa+41gn/jlqCWBUBwEFABAJaSuv8Up4oaVCA0EVABAJZSv9RZpJikHwhFBFQAgHWb+DcxxRQQigioAADrNvEzByoQkgioAABLSdvfxJ+za7cU5OYGuzgAgoCACgCwjPDICKnXuJFZp3kfCF0EVACAZdRr0ljCIyLM+s6N9D8FQhUBFQBgGfX3n+JU7dzACH4gVBFQAQCW63+qaOIHQhcBFQBgGanMgQqAgAoAsGwTP31QgZBFDSoAwDLSvAPqpi1BLQuA4CGgAgAsI7VpE3OZmbFdigoKgl0cAEFCQAUAWEJkTIykNGpo1mneB0IbARUAYAlp+2tPFSP4gdBGQAUAWEJqU68pppgDFQhpBFQAgPXmQN3IJP1AKCOgAgCsN8UUI/iBkGaJgPqf//xHZs+eLZmZmbJ161b57LPPpE2bNsEuFgAgWJP008QPhDRLBNQ+ffrImDFjpEePHtKvXz+JioqSb775RuLj44NdNABAgJv49xUXy+4tW9nvQAiLFAvo379/qetDhgyRjIwMSU9Pl5kzZwatXACAwDfx79mWIcVFRex6IIRZIqCWlZKSYi537txZ4TbR0dESExPjuZ6UlFTqsq4F+vXgfxxDe+P4OesYRsfHSUJqPXN9z5Zt/LbaBN9D+0sKUn6qSpiIuMRCwsLC5Msvv5R69erJ8ccfX+F2w4cPlxEjRgS0bACAupGRnyvvrlxo1jvUayCnNT+cXQ04WHJysmRlZdknoI4dO9Y0+R933HGyceNGn2pQdftmzZpV+ob9JdCvB//jGNobx89Zx7DpUUfIxU+PNLf/+Na7ZoH18T20v6Qg5aeqAqqlmvhfeuklOfPMM6V3796VhlNVWFholrL0zQYyMAb69eB/HEN74/g54xjG7m/eV5tXreF31Wb4HtpflsXyTKSVwumAAQPkhBNOkDVr1gS7OACAAEpr7jUHKpP0AyHPEgFVp5i65JJL5JxzzjHpvXHjxub2PXv2SH5+frCLBwCoY/VLnUVqM/sbCHGWmAd16NChZlDUTz/9JFu2bPEsgwYNCnbRAAAB0OTww8zl3vwCM80UgNAWaZWR+wCA0BSblCgNWjQ36xv/Wi6uffuCXSQAQWaJGlQAQOhq2ra1Z339kmVBLQsAayCgAgCCqmn7tp719Uv+CmpZAFgDARUAEFRN23nVoC5eGtSyALAGAioAIKiatm9jLvNzciRj7XqOBgACKgAgeHKL9kq9JiVTC25Y+jcDpAAY1KACAIJmS162Z30D/U8B7EdABQAEzda8HM86I/gBuBFQAQBBs9WrBnX9YqaYAlCCgAoACHoNam5mpuzYsJEjAcAgoAIAgiK5YQPJKdpr1ul/CsAbARUAENTppdQ6mvcBeCGgAgCComm7fwPqBk5xCsALARUAEPQaVE5xCsAbARUAENQa1Oydu2T3lq0cBQAeBFQAQMClNTtI4lOSzfqmv5ZzBACUQkAFAATcwR07eNY3LSOgAiiNgAoACLiDO7TzrBNQAZRFQAUABFzzI7wC6t8rOAIASiGgAgACKiwsTJrvr0FNjIyW7B07OQIASiGgAgACqkHLgyUuKdGsN45LYO8DOAABFQAQUAd3bO9ZbxJXElQBwBsBFQAQUAd3+DegUoMKoDwEVABAQB3sNUCKgAqgPARUAEDAhEdESLP2bc36ro2bJS4yir0P4AAEVABAwDQ6tJVEx8Wadc4gBaAiBFQAQMC08DqD1EbOIAWgAgRUAEDAdDzxeM/6+kVL2PMAykVABQAERGxigrTtdYxZ37MtQzYsXsaeB1AuAioAICCOOOF4iYyONusLv50uLpeLPQ+gXARUAEBAHHVKX8/6n9/8wF4HUCECKgAg4M37axYsZK8DqBABFQBQ52jeB+ALAioAoM7RvA/AFwRUAECdonkfgK8ifX4EAAA+CJXm/fj4eGnQoIGEhYVJqEhISJD8/Hw5+OCDJScnJ9jFQZCPoX63t2/fLrm5uVJbBFQAQOCa97/+3nF7WwPplVdeKSeccIKEmvDwcJk9e7bcd999sm/fvmAXBxY5hj/++KOMGzeuVn+MElABAIFp3t+aIWv+WOS4va3htE+fPjJp0iT566+/pKioSEIp3LRv316WLVtGQLWpcD8ew8jISGnXrp1ceOGF5vrbb79d8+eqVUkAAKhEhxOO+7d5/zvnNe9r86jWnGo4nTJlioRiuElLS5O1a9cSUG0q3M/H8J9//jGXgwYNkg8//LDGzf0MkgIA1JnOp5zk6Ob9+vXrm0utOQUgpb4P2ie7pgioAIA6EQrN++4BUaHUrA9Uxf19qM2AQQIqAKBOOL15H0DdIaACAOqE05v3AdQdAioAwO8SUus5vnnfCRo3biyjR482A1t0Lsx169bJl19+KX37/js1WCBpLbsuxxxT8tlxi46ONvNr6n06YwKcj4AKAPC7ftdf6Wne/+Pr72jet6CWLVvKvHnzTBi9++67pVOnTnLaaafJ9OnTZcyYMTV+3qioqFqVS0OyTt3lbcCAAZKdnV2r54W9EFABAH5V/+Dm0vPC88x6QW6eTB/3PnvYgsaOHWv+cOjevbt8+umnsmLFClm6dKk8//zz0qNHD892eoahzz//XLKysmTPnj1mSq1GjRp57r/uuutM0L366qtl1apVpiZWpaSkyBtvvCHbtm0zj/v+++/lyCOPrLJc48ePl4suukhiY2M9t1111VXm9rKefPJJ+fvvv80ZkLQWeOTIkWYuTjd9vR9++EEyMzNNGebOnSvp6enmvhYtWpja4p07d5rwu3jxYunfv38t9ij8iYAKAPCr02+9QSKiSkLCT+MnStb2Hexhi0lNTTW1pVpTWt48lRrm3KOwv/jiCzNPpjat9+vXTw499FATUr0dfvjhMnDgQDnvvPOkc+fO5raPP/7YBFkNfRoK58+fb0KqvnZlNOyuWbPGPJ87IPfu3VvefffdA7bV0DxkyBDp0KGD3HrrrXLttdfK7bff7rn//ffflw0bNki3bt1MGTTQ7t2719yn7z0mJsY8t9Ye33vvvdTSWggT9QMA/KbFkUdI51NLBkdl7dgpP74zMST37m0fvi1JDUrmSA0U/UPghYuuqta2Gih1gvaq5m896aSTTHg75JBDTNBTV1xxhalp7dq1qwmd7j6ierv2E1W9evUyNbMaUAsLC81t2o3g3HPPlfPPP9/UrFZGz0CktaYaMDWATp06VTIyMg7Y7rHHHvOs60Tzzz77rKl9feaZZzy1pLqutaxq5cqVnu31vk8++cTUnKrVq1dXa98hMAioAAC/OeuO//Osfz3mTSmo4Vlk7E7Dab3G/zaDW01156fUU2CuX7/eE06VnhJz165d5j53QNVw6A6n6qijjpLExETZsaN07XlcXJwcdthhVb7ue++9Z2o7NRhrQL3lllvK3U5Pqan36XPq62nzvjbnu40aNUrefPNNufzyy+W7774ztbraDUHp4LBXXnlFTjnlFHOfhtVFixjMZxUEVACAXxxx4vFyaHpJ8+621Wvl98++DNk9G4xuDb68pvY31dNa6nnT/UH7gHrTsLh582ZzGtiydu/eXeXzab/Q//3vf/LWW2+ZvqhfffWVJCUlldpG+8lqDevw4cPl66+/Nt0StPb0zjvv9Gzz8MMPy8SJE+WMM84wXQ30um6jfWr1ufVxep+G1GHDhpnHvvzyy7XaF/AflxOWpKQkl9JLJ74eC8eQzwDfQSt/BsIjI1z3fvmh67lFv5qlY9/eIfE72rJlS9eECRPMZbDL4usydepU1/r1613x8fEH3JeSkmIuTz75ZNfevXtdzZs399zXvn17c9zS09Nd4eHhrtdee821YMGCUo93P87X/aLOOeccs37aaaeZ60888YSnTKpPnz7m+h133OFauXJlqce/8cYbrl27dlX4/BMnTnR98cUX5d73+OOPu/7888+gH5dAL+Hh4Z5jGYjvRXW/9wySAgDU2jEDzpZGh7Q066vm/SGLf5jBXrW4m266SSIiImT27NlmcJP2S9Ua1Ztvvll+/fVXs402fWuzt9ZUHn300Waw0YQJE+THH380g5kqoo/T59CaSh1YpVNaHXvssfLoo496RtFXZdq0aeZc7g899FCFtcDaj3TQoEFm4JaWW6ejctOa15deeskM7tLtevbsacqvXRSUzlagNaetWrUy7+3EE0/03IfgI6ACAGolJj5eTr3pGs/1yaNoIrUDHRTUpUsXM+/pc889ZwYLffvtt2Zg1I033ujZ7pxzzjF9TmfMmGGCp/bh1FBYldNPP908Zty4cbJ8+XL58MMPTVDdunVrtcuofVjdo+7Lmjx5sgmZ2iT/xx9/mAD6yCOPeO4vLi6W+vXrm0Ctr//RRx+ZrgLaJUBpONeR/BpKNQzrNkOHDq122VD3XE5YaOIP/jGw2+KE5sVQXjh+1ln633KDp2n/8mcfDaljaOcmfqs2D7PY/xjSxA8ACKouZ54qJ1872KwX7y2SqS++yhEBUGs08QMAaqRtz2PkopEPeK5/PfZN2bH+3+mIAKCmCKgAAJ8dfER7Gfz8454zRs368BP5/s0DT0UJADVBQAUA+KRBy4PlmrHPmcFR6s9vfpDPnhjFXgTgNwRUAIBPZ0i67tUXJDGt5HzqK+fMl4nDHhbXvn3sRQB+Q0AFAFRLQmo9ufaVUVK/eVNzfdPylTLu1nulaP+51gHAXzjVKQCgSh379pHzH7pHkuqnmes7N26WN264XfKzstl7APyOgAoAqFBccpIMGHaHpJ95mue2rB075fUbbpPMjO3sOQB1goAKAChXu+OPlQtHDJOURg09ty3+4Sf578inTUgFgLpCH1QAQCnN2reRS58cIdeOHeUJp3mZWWYw1Lhb/0M4xQH0FKYul0uOOuooSzyPE0yYMEGGDRsW0P2mp73V08dW5tRTT5UFCxZIWFiY1CUCKgBAIqKizFmhbn7vdbnjo/HS5YxTPXvlr59/k2fOu1Tm/W8ae8pBxo0bJ5999pmlXn/9+vXSpEkTWbx4cZ2+9vDhw02g++qrrw6476677jL3aVgLliOPPFJOP/10GT16dLXCo7/223nnnScPPvig5/rq1avl1ltvLbXN119/LXv37pVLL71U6hJN/AAQorQGpGnb1nLUqSfJMeed5Zk6yk1rTSePell+/+TLoJURoWXfvn2ydevWgLzWpk2b5MQTT5RmzZrJxo0bPbdfddVVsnbtWgmmm2++WT7++GPJyckJ6H7btWtXtbZ755135JZbbpH33ntP6go1qAAQQhq0aC7HXjBArnjuMXn4p6lyx8fj5aRrrigVTnX6qI9HPiUjTz6HcBpCBg4cKAsXLpTc3FzZvn27fPvttxK//2QM+seM1qxpTV1+fr5p4tWm3oqEh4fLm2++KatWrTLP99dff5lA412DOWTIEDn33HNNbaUuffr0Kbepunfv3vL777+b19VQ+cQTT0hERESpmsUXX3xRnnrqKdmxY4ds3rzZPH9Vtm3bJt98840MHjzYc9uxxx4rDRo0kClTppTatmvXrmbbjIwM2b17t/z4449y9NFHl9pGX1ODrZZTA6+Wye3GG2+U5cuXS15enmzZssWEz8r23fnnny+TJ0+W6iq733Rf6vVTTjlF5s+fb47B999/Lw0bNpTTTjtNli5dKnv27JH3339f4uLiyq2l1fVWrVrJCy+84DlGblq2bt26yaGHHiohUYM6dOhQufvuu0019Z9//mn+gpgzZ06wiwUAtqNneWp8WCtpfNgh0uSwQ6Xx4YfIQa0Pk3qNG5W7ffHeIln47Q8ya9Knsnr+nwEvrxNFRUVVeJ/+Z19UVOS3bbXJtTb0/90PPvhA7rnnHtPsnpSUJMcff7ynn6E28955551y/fXXm3CqtYxffvmldOrUqcKQtWHDBrngggtMaOzZs6e8/vrrJjxqOHv22Welffv2kpycLFdeeaV5zM6dO6Vp05I5dt30+tSpU02N3RVXXCHt2rWTN954w4TAhx9+2LOdhsxRo0bJMcccY0Kmbj9r1iz57rvvKn3fb7/9tjz99NPy+OOPm+v6vjS0laX7Y/z48SaX6D7RfaHlat26tWRnZ5twf/vtt8tFF10kS5YsMfvTHRbT09NNU/3ll18uv/zyi6SlpZl9W1nzfr169WTu3LlSWyNGjJD/+7//MwH1o48+MktBQYFccsklkpiYaI613l/eftLmfs1ietx0n3vTP1Q0aOv70D9CHB1QL7zwQvPhuuGGG8xfSrfddpvp59C2bVvzFwsAQCQ2MUHiU5IlPiVFEuqlSEJaPTOQSYNnillK1pMbNqhyd2kTvp4JauXsufLnN9Mla/sOdrEf3X///RXep7VpEydO9FzXypno6Ohyt12zZo0JXG76/2NCQsIBQaQ2DjroIBN8P/30U1m3bp25zbs/o/bL1BrKSZMmmev/+c9/TPO4BlcNeWVpoPYuk74HDY76f7276VprE2NiYiptmtaKKw1DGqLU33//bUKrlmXkyJGeWj2t+dXrauXKlWb7k046qcqA+r///U9effVVE7TmzZtnynfccceZoOqtbH/U6667ztSkak2l1ra2aNHCBDZ9PX3vWmZ3BZvep+9XX0vDrO7fP/74o9La0KKiIlPDW1sPPPCACcXqrbfekieffNLUemrfUvXf//7XHMfy9pM29xcXF0tWVla5x0hrs7WsdcUyAfWOO+4wCd39JdSgesYZZ5gPiX4QAcCfTM1QWJiEhYeJ/itZD5dwvR4WXnJ7eLjZruT2cAmLiJBwvR6h1yP2X4ZLeGRkyaVej4w0zY/mMtJ9uX+JipLIqEiJiI42l5FR0RIZGyNRMbpEl1zGxkh0XJypAY2Jj5OYhPiSJT7ezEmqz1NTGkjXL1kmK36fKyt+mysblv3NKUphaE2ZhpRFixaZyiFtztbwoiFMaw+1n6bWSHrT65WNGtdwqf+Ha0DTZmQN4JUFs/JoLeuvv/56wOtqmZo3b26CoDugetOa2kaNym8t8KZBUPtRai2uBjf9w0H3QVn6XI8++qiccMIJZl2/49r9Qd+b0tCtfzhobeK0adNM7ao2g2vA064S2vTvvk8XrbnUgF4e3VcFBQXiD977RUOmBmV3OHXf1r179xo9t5bf3QXEsQFV/2rTKnDtV+KmfxXpl0X/4rKaNsd2l/43XSsfrloiV70yynwAYT/6A+PzMfRhWg2fp+AoZ3MTnHwoh+c1wyooh9e693OX3sb7eTSsldlu/3UT3LzWPdvuf7yGPF3X4Oa+z/t5y31tr+tlt/V+byZShofLy0vnyLBvP9//fsuUy/0+9LrZvmTdlMfB9hUXS+b2HbJr0xbZ+s9q2fLPatn6zyrZsnI1E+sH2GOPPVbhfd79+dQzzzxT7W21T6C/6SCbfv36maZ47beoTdlafm0y1yZ6Xw0aNMg042tTuAZMrYXTWmJ9vrpQtouD7rPqfte1Blhbbjt27FhubbDS5v369eubGmMNmxog9X25a721O4O2+J588slmP44dO9a8X61h1VrTLl26mHCr+1ZrerV2Wftwaj/QsrT/b0JCgslGte264f143Se12U9laVeFumzhtkRA1Q7JkZGRB1Qh63Xtb1Ie/VBo04Cb/jXlfVmXGjQ9SFoc1VE25WaZS9gXx9DeCvcVm9pFJysqKJSC3FwpzM2TvKxsUwual5lpLnP3ZErunj2SlbFDMrdlSOa27ZK9c6fsK95X7nMF4vfRF4H83a4rGiRM7fn+xVtVf/h6b1/bbWsSMtx/zHk/9rfffjOL1hZqTZv2Q9RArIN+tOn7559/9mzbq1cvTz9J96Al937QbbVp+bXXXvNsf9hhh5Uqq4Yl/b/f+/Xd6+7n0cFVWgbvbfS5MzMzTROz+/ay76O891b2vbtfR19D+41q388PP/ywpLXE6373e9VuA1q7rLT2Vgcceb9GYWGhqTnV5ZVXXpFly5aZGmbts+uetkqXRx55xPS31TBb3jRf7lrPjh07mprtsuUu7z2V3W9lr5fdpux+8j6GZd9T2WOkNH/p8dTyVVQeXfT7Ufb7Xd3vuyUCak3o5LXl9bnxniqirizbvV2+2rCyzl8HqIny6lZL1cSWv1rpY0pXsHrXwv674r49zPuHz+sxB1bulrdW+oayNcje10qer/QWnmslFahej9l/j1cZ3GXyfgbTfO/eOsx7u/23e+7f/yOuy/77w8vcZi69rkeEh0lEWLhEhOml3l6yHhkeLpFh4QdcRkdESFS4buPsGt9A/W7XFR2sM3v2bNMUrTVKdqLl1bCgo9GPOOII09Sr4VTDk4YjbcrWEKn3a3DT/3e1plWbwc866yxzu7vl012ZpJcaTLT5V2tLdfS6Bkmd07NHjx7mWLtHv2v40dZTHcmvXQm0ptHdLO9+npkzZ5rBRzqAS/u/6qhy7Vep5encubPZVt+DhkXvUfUpKSnm8WVH2nv3udWmdPf9WtOrQcw9Kl2fz71vlJZb34s2kWvo0ppUPfYaVHWbM8880wQ87bert+v+0Usth9ZGaxcJDaoarDXsugNcReVbtmyZGXDlHf60PG3atDG102VrXN0Vdu79poO3lIZu3a9KuyNoGb1fU/dDbGysZ9uy+1I/C/peNMDr8XLX+Opx08+GDr4q7z00btxYDj74YNOv1/38vrJEQNWdq/1A9A150+va6bg8+qXQQVVuulP1A6QfAm1KqEv6n1BScrLp6KwHvK5fD3VDPzM1OoZlmttqomyTHXwXyO886oYTjqH+J3zfffeZQBHsuTN9peFDA6cGJw2UOrJbR93ryHp9LxratKlaad9RPUY33XSTCZE6TdE555xjuuJpsNFaSKWXWqumgUaDjtbE6u+dBsqXX37ZTHGkr6e0/6gGLh17op+Fvn37msd5P4/ScKtjUXRQmZZZR5XrlFfummQtlzY1u59XaZDS0Ot9W9k+qvqeK7pfn08HY7nv10npdTDVu+++a8qtIVlH/2vTvm6j/4/oDAg6lZaGQO3HqsFOp6PSIKch/OqrrzZhbcWKFeb5PvnkkwqPzcsvvyyXXXZZqTNJ6fvs37+/WbzpvnDPPODeb3oM3bWx7lCpx0n3mfd71v2gQdq9bdl9qeODtDZYa3q17O5aVh0kpn139Q+a8ujgKd1PGurd/YTLfu+rw2WF5bfffnONHj3acz0sLMy1fv1617333lutxyclJbmUXgaivIF+PRaOIZ8BvoNO+ww44Xe0ZcuWrgkTJpjLYJclGEt4eLgrPT3dXAa7LE5aYmNjXWvXrnX16NHDcsewfv36ru3bt7tatWpVo+9Fdb/3lqhBVVobqp2QtT+LNpe4p9HQU6EBAACEivz8fDPvq47RsRrtZqEzNOjUYXXJMgFVJ4/V5gAd3aYT3GpzgjYF+GMeMAAAADv56aefxIq0X6kudc0yAVWNGTPGLAAAAAhdzh8eCgAAAFshoAIAUMsZOXSKIgBS6vtQmxlrCKgAANSQ+yxLFZ1UBghF7fZ/H3Qa0ZriTz4AAGpIJ27XuS4vvPBCzzyUOq93qNBJ4XXOcp33UudURWgfw8jISBNO9fug3wudyL/Gz1WrkgAAEOLc0yGWPcNPqIQbPVmBTsZOQLWn8Do4hhpOaztNKAEVAIBa0H52b7/9tjlbks5b6T63eSjQ+cp1yiH3aUAR2sfQ5XKZZv3a1Jy6EVABAPAD/U9ZT58cSvS0lXoKTK19s+vpakNdkkWPIYOkAAAAYCkEVAAAAFgKARUAAACWEunEvhSBfJ1AvR78j2Nobxw/++MY2h/H0P6SApxnqvs6OtSw5tP8W0jTpk1l48aNwS4GAAAAqtCsWTPZtGmT8wOqO6S6R6DNnj1bunfvXuG2Fd1f3dv1LwANxLqDgz3qrar3Gqjn8+Vx1dm2sm1qch/H0D/Hpbrb1ub+8u7j+NX+mPi6Lb+jNcfvaODxf6Hvv6PByjP6mpWFU8c18Xu/WZ1strIdXdH9vt6utwU7oFb1XgP1fL48rjrbVrZNTe7jGPrnuFR329rcX959HL/aHxNft+V3tOb4HQ08/i+s+e9ooPNMdV7HsYOkxowZU6P7fb3dCvxdtpo+ny+Pq862lW1Tk/s4hv45LtXdtjb3l3cfx6/2x8TXbfkdrTl+RwOP/wvt/Tvq6Cb+QNLq6czMTElOTg56DSpqhmNobxw/++MY2h/H0P6SLJpnHFuDWtcKCgpkxIgR5hL2xDG0N46f/XEM7Y9jaH8FFs0z1KACAADAUqhBBQAAgKUQUAEAAGApBFQAAABYCgEVAAAAlkJABQAAgKUQUOtA8+bNZfr06bJkyRL5888/5fzzz6+Ll0EdSUlJkTlz5siCBQtk0aJFcs0117CvbSouLk7WrFkjzzzzTLCLghpYvXq1+Q3V7+IPP/zAPrSZVq1ameOm/xcuXLhQ4uPjg10k+KBNmzbmu+decnNz5ZxzzpFA0on6Wfy4D5o0aeI66qijzHrjxo1dGzZscMXHx7OPbfI5Cw8Pd8XFxZl1PW6rVq1ypaWlBb1cLL7vg0cffdT14Ycfup555hn2nw0/Q6tXr3YlJCQEvRwsNdsHP/74o+u4444z66mpqa6IiAj2pU0/TwkJCa6MjIyAZhlqUOvAli1bzF/9auvWrbJ9+3ZJS0uri5dCHdBzFefl5Zn1mJgYCQsLMwvs5fDDD5d27drJV199FeyiACGnQ4cOsnfvXvn555/N9V27dklxcXGwi4UaOvvss+X77783taiBQkAtx/HHHy9ffvmlbNy4UVwuV7lV2kOHDjXNTxpkfvvtN+nWrVu5O7hLly4SEREhGzZs8P/RQ50dP23m/+OPP8xx0+bhHTt2sLdtdgyfffZZGTZsWABLDX8fQ33cTz/9JLNnz5ZLLrmEHWyj49e6dWvJzs42zzFv3jy+izbPMhdeeKFMmjRJAomAWo6EhARTA3rTTTdVeKBGjRolDz/8sAmguu3XX38tDRs2LLVdamqqTJgwQa677rq6OXqos+O3Z88e6dy5sxxyyCHmP8ZGjRqxt210DPWv/eXLl8uKFSs4bjb+Hh533HHStWtXczzvu+8+6dSpUwDfQWir7fGLjIw0AUkD0LHHHiv9+vWTk08+OcDvIrQl+CnLJCUlSc+ePWXq1KkSaEHv22DlRZ1zzjmlbvvtt99cL730kud6WFiY6Wd67733em6Ljo52/fTTT67LLrss6O8hlJeaHj/vZcyYMa6BAwcG/b2E6lKTY/j444+71q1bZ/owar+p3bt3ux588MGgv5dQXfzxPXz66addgwcPDvp7CcWlJsevR48ermnTpnnuv+uuu8wS7PcSqkttvoOaY959992Al5kaVB9FRUVJenq6fPfdd/8mfJfLXNe/Et3eeecdM3rxvffe89+fEgjI8dPa0sTERLOenJwsvXv3lr///pu9b6NjqLVtLVq0MDXgd911l7zxxhvyyCOPBLHU8PUY6ohv9/dQa4L69u1rRoPDHsdPZ0LR39J69eqZPvz6O7ps2bIglho1yTLBat5XkQF/RZtr0KCBabrQwU/e9LoOyFC9evWSQYMGmWk1zj33XHPb5ZdfLosXLw5KmeHb8WvZsqW8/vrrnsFRL730EsfOZscQ9j+GjRs3ls8++8ysaz9+/SNj7ty5QSkvfD9+OiBK/1CcMWOG+R395ptvZMqUKexKm/2OJicnS/fu3WXgwIEBLyMBtQ7MmjXL/KDCnvQv/6OPPjrYxYCfjB8/nn1pQzpwQ/uBw76mTZtmFthXZmamNGnSJCivTRO/j3TKqKKiIvPXvTe9rtNLwdo4fvbHMbQ/jqG9cfzsb7sNsgwB1Uc6r5tOmXHSSSd5btPmC73+66+/+vv4wM84fvbHMbQ/jqG9cfzsb69NskzQR5dZ8YwJeiYoXdRtt91m1g8++GBz/4UXXujKy8tzXXHFFa527dq5Xn31VdfOnTtdjRo1CnrZWTh+TvgM8B0M/jHgGAZ/H3L8gr8f7Lwk2D/LBL0Allv69OnjKs+4ceM829x0002uNWvWuPLz881UDd27dw96uVk4fk75DPAdDP4x4BgGfx9y/IK/H+y89LF5lgnbvwIAAABYAn1QAQAAYCkEVAAAAFgKARUAAACWQkAFAACApRBQAQAAYCkEVAAAAFgKARUAAACWQkAFAACApRBQAaAc48aNk88++yzg+2bw4MHicrnM8vzzz1v22PTp08dTzmDsJwDOFhnsAgBAoGmoqsyIESPk1ltvlbAwPdle4O3Zs0fatm0rOTk5Eh0dLfPnz5dZs2bJ9ddfX2q7p556Si644AI58sgjJTs7O6Bl/OWXX6RJkyby4osvSkxMTEBfG4DzEVABhBwNVm6DBg2SkSNHmkDopmFPw2EwA/TWrVs916+44gr59ddf5ZNPPpFvvvnG3HbMMcfI7bffLieffHKdhNOIiAgpLi6u8P69e/eaMubl5RFQAfgdTfwAQo4GK/eitZXuQOheNJyWbeKfPn26jB492jS779y5U7Zs2SLXXHONxMfHy9tvvy2ZmZmyYsUKOe2000q91hFHHCFTp06VrKws85gJEyZI/fr1fSqv1qA+9thj8tZbb0lKSooJhFq+l156SWbMmCG9evUyl7m5ubJu3TpTq6nlcrvssstkzpw5poybN2+W999/Xxo2bHhAc72Wfe7cuVJQUCDHHXecqZn94YcfzON0P+l96enptdr3AFAdBFQA8KF/6Pbt26V79+4mHL7yyivy8ccfm+buLl26mNrNd999V+Li4sz2GiY14C1YsEC6du1qAmDjxo3lo48+8nmfa0DVgKsh+dFHHzWB8r777pNDDz1Upk2bZmpXNVBqjbCGy5dfftnz2KioKHnwwQflqKOOknPPPVdatWol77zzzgGv8eSTT8p//vMfad++vSxcuNAE2Q0bNki3bt1MMNX7teYUAAJBO2OxsA/4DPAZCMnPwODBg127du064PZx48a5PvvsM8/16dOnu2bMmOG5Hh4e7srKynKNHz/ec1vjxo1d6phjjjHX77//fte0adNKPW+zZs3MNq1bt/apPLq0b9/elZub68rPz3elp6eb29544w3Xq6++Wmq7Xr16uYqKilwxMTHlPo8+ViUkJJjrffr0MdfPPvvsUtvt2bPHdcUVV1S6/8ruJxb2AZ8BPgPih31AH1QAqCatVXTbt2+f7NixQxYtWuS5zd1vtFGjRuZSayxPPPFE07xf1mGHHWa6BPhi2bJlpqa0Xr16Mm/ePM9raM3ppZde6tlOB3dpH9JDDjlE/vrrL1O7qwO/dNvU1FQJDy9pPGvRooV5Tjdtwvc2atQoefPNN+Xyyy+X7777ztQWr1q1yqcyA0BN0MQPANVUtnlbm9nLa/J2B8DExESZPHmydO7cudRy+OGHmz6jNVFUVGQWN32N1157rdTzaxDV1/jnn39MX9Svv/7a9CPVEKvN9QMGDDCP1RkCvJUdGPbwww+bPrRTpkyRvn37ytKlS00XAQCoa9SgAkAd0cFNAwcOlDVr1lQ6Ir62r9GhQwcTRsvTqVMnadCggelbqv1JlfaHrS6t5X3hhRfMMnHiRLnyyivl888/91v5AaA81KACQB0ZM2aMpKWlyQcffGBCoQ5oOuWUU8yof3cta23pXKg9e/Y0g7bcNadnn322ua50VL+Oyr/55ptNk/9ZZ51lBkxVJTY21jyHjvDXrgD6Glr76t0lAADqCgEVAOqITumkU0Bpf1Ad4a/9VbUmcvfu3aYPqz/oc2qIbNOmjcycOdPMGKDzum7atMncr7MODBkyxEzor030WpN61113Vfm8WuOr02HptFjLly83Mw989dVXMnz4cL+UGwAqE7Z/tBQAwCJTWWmI1cFMdqDzseqgLXe/VgDwB2pQAcBiNPDpyH+dd9SqdK5VLaP37AEA4C/UoAKAheiofJ3MX2lXAJ3Kyoq0j2qzZs3Mup5q1fvUrABQWwRUAAAAWApN/AAAALAUAioAAAAshYAKAAAASyGgAgAAwFIIqAAAALAUAioAAAAshYAKAAAASyGgAgAAwFIIqAAAABAr+X86JgmabVmLUQAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Extract data for the embryo at 3 AU\n", "history_3au = results[3.0]\n", "\n", "time_years = history_3au[:, 0] / (3.154e7) # Convert s to years\n", "total_mass = history_3au[:, 1] / c.M_EARTH # Convert to Earth Masses\n", "iso_mass = history_3au[:, 2] / c.M_EARTH\n", "\n", "\n", "plt.figure(figsize=(8,5))\n", "plt.plot(time_years, total_mass, lw=2, label=\"Core Mass\")\n", "plt.plot(time_years, iso_mass, '--', color='gray', label=\"Isolation Mass (Limit)\")\n", "plt.xlabel(\"Time [Years]\")\n", "plt.ylabel(r\"Mass [$M_\\oplus$]\") \n", "plt.xscale(\"log\")\n", "plt.title(\"Embryo Growth at 3 AU\")\n", "plt.legend()\n", "plt.grid(True)\n", "plt.show()" ] }, { "cell_type": "markdown", "id": "bae8c214", "metadata": {}, "source": [ "## Water Mass Fraction\n", "\n", "By default, PA3Py runs the classic `silicates` and `H2O` model by tracking the dynamic snowline position.\n", "We can extract the water mass (column 4) and see how the fraction increased." ] }, { "cell_type": "code", "execution_count": 11, "id": "cbf4fb86", "metadata": { "execution": { "iopub.execute_input": "2026-07-13T22:23:20.036602Z", "iopub.status.busy": "2026-07-13T22:23:20.036386Z", "iopub.status.idle": "2026-07-13T22:23:20.279814Z", "shell.execute_reply": "2026-07-13T22:23:20.278648Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAq4AAAHbCAYAAAAK+BjpAAAAOnRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjEwLjgsIGh0dHBzOi8vbWF0cGxvdGxpYi5vcmcvwVt1zgAAAAlwSFlzAAAPYQAAD2EBqD+naQAATfpJREFUeJzt3QmYFMX5x/Hfci0CyyUCgoooAl7RAIJcIiCKIl4o/EWDeESNmkSNiWI0gvcVENHEC0WNBypiMCinKAooigp4gASQS0C578Vl+/+8Zc9mdtld9piZnp75fp6nmZ6eo2uqepZ3qt+uypDkCQAAAEhyFYIuAAAAAFASBK4AAAAIBQJXAAAAhAKBKwAAAEKBwBUAAAChQOAKAACAUCBwBQAAQCgQuAIAACAUCFwBAAAQCgSuQIx4nqcRI0YkVX0uXbpUzz33XFw/8x133KFkNW3aNLfEUpMmTdznvuSSS5TMTjvtNH3xxRfauXOnK2+tWrVCf0zHoz0BhAuBK7APhx12mJ544gktXrzYBQGbN2/WRx99pD/84Q+qWrUq9VfCQK+o5eabb07KOrzwwgv1xz/+UWFUt25dvfbaa+54vfbaa3XxxRdr+/bthT63ffv27sdHWQLbWP7Aij4m1q5dq+nTp+ucc85RsjryyCNdvdnxHUulbQ+rowkTJmjVqlXatWuXVqxYoddff11HH310qff9ySefuPq/+uqrC33cymWP77///oU+Pn/+fH5YIO4qxX8XQHidccYZ7j+B7OxsvfDCC/rqq69UpUoVderUSQ899JD7z+Gqq65SsmrRooVyc3OVDF5++WW98847e223XsFk1L9/fx1zzDEaPnx4vu3Lli1zP1h+/vlnJasTTjhBNWvW1O23366pU6cW+9wOHTpo8ODBGjVqlPtRFhQ7Dv7+97+79UaNGrnv1dixY10Q9eSTTyrZHHXUUa7e3n//fXdMxEpp2+PYY4/Vxo0b3XG6bt06NWzYUJdddplmz57tguB58+aVaL/NmjVT27Zt3Y+Iiy66yP1YB5IRgStQhEMPPVSvvvqq+0+pW7duWrNmTd5j//jHP3T44YerV69eSV1/u3fvVrL4/PPP9dJLLykV2A+ZZFa/fn13u2nTJoWF9RhGHx/2Q/G///2vbrjhhqQMXJPFXXfdtde2Z555RitXrtTvfvc7t5SE9cpbT/ef/vQnvfHGG64nOZYBORArpAoARfjLX/6irKwsXX755fmC1ghLHXj00Uf32n722We7U2Z22s56aC3XsCDrURo5cqR738jzLr300nzP6dKlizstd8EFF+hvf/ub+49oy5YtrgfYetOs53fYsGHuP5utW7fq2Wefddv2leNqpyCHDh3qHoucWnz++efzTv9VrlxZQ4YM0WeffeYCn23btrnTtieffHJcj5W3337b1WlhZs6cqU8//TTvfsWKFXXbbbe5wMY+g32We+65Z6/PX5DlpVqdFjy9G6lruzWWR3nmmWe6Hy+R09e2j+JyXLt27erqyerLesDeeusttWzZstBTrfajx9rFnmd1bG233377laiezj//fNc2O3bs0E8//aQXX3zRHU8RVnYL+ow9z/ZXVJ6zlefhhx92699//33eZy1YP7E6pkvDjutvv/1WTZs2LfI5JT1WI21mQdlvf/vbvOPGeiXbtGlT6JkK+56tX7/epVvYsde7d++8x63tLbgz1uMaqbfI8VNUz6i1QyTlaPXq1a6+LK2jtO2xLz/++KM7PmrXrl2qMwz2mf7zn/+4nl67DyQjelyBIth/VPafzKxZs0pcR5ZCcN5557keWQsmLQ92zJgxOuSQQ7Rhw4a83rCPP/7Y/Yf02GOPueDj9NNPd8GLBaQFT00PGjTI/Ud3//33u9N5v//9791paksBqFOnjjuteOKJJ7ogwYKrwnpgIqpXr64PP/zQ5efZ/qwXtF69ejrrrLN00EEHuf+orQxXXHGFXnnlFT399NN5wfvEiRPdqcS5c+eW6ZipVq1aoblxFnDs2bNHo0ePdkGYBRIWiERY3dkpz5tuuilfj9LAgQNdcGGnl9u1a6dbb73VfS6r//KyINgCfKsT6/EzFhQVpXv37nr33Xe1ZMkS1x4WhFo7zZgxQ61atdqr58ryT62trG3tcQumLNi45ZZbii2XBUx2CtkCLnttgwYNXB5ux44d9etf/9oFHFb2hQsXulPtlipg+ynqB8Gbb76p5s2buyDl+uuvd6eajR2T8TymS6JSpUo6+OCD3TFZlNIeq/Y57TnWg2tltR+nVgeWx56Tk5OXAmDtZj3A9p2z3OC+ffu6HyJ9+vRxtxYc22eyurf6tgDbRG4L06NHD7cfC14tuLc0oyuvvNLd2ve3pO1RFDteLZC3VAF7rd3fV5pIhNXVEUcc4f6G2N8WK4elC9x3330lej2QaB4LdcAxkP8YyMrK8szYsWNLXDdm165d3mGHHZa37dhjj3Xbr7322rxtTz/9tLdq1Sqvbt26+V7/8ssvexs3bvSqVq3q7nfp0sW9dt68eV6lSpXynvfSSy95e/bs8caPH5/v9TNmzPCWLl2ab5vdf+655/LuDx482L3nOeecU+TnqFChgle5cuV822rVquWtXr3ae+aZZ/b6zHfccUex9dKkSROvOO3atcur8507d3oPPfRQvtffdNNN7vMefPDB7v6vfvUr97qnnnoq3/MefPBBt/3kk0/O2zZt2jS3RO5fcskl7jlWpujXRurabiPb3n777b3qM/rz2HtFtn3++efemjVrvDp16uRr+5ycHG/UqFF526yuTMF6HDNmjPfTTz8VW492DNg+7HjIzMzM237GGWe497S2Lfg5W7duvc/j9k9/+lOhdRKPY7qoxep5woQJ3v777+8W24e91gwfPrzI9izpsRppM6vj2rVr523v3bu3296rV6+8bZMnT/bmzp3rValSJd/7fvTRR97ChQvz7vfp02evY6a4pbA66Nevn3uPTp06lag9ilu+/fbbvO/Uli1bvDvvvNPLyMgo0WsfffRRb9myZXn3TznlFPc+xx13XL7nRY5fa6PC3mf+/Pn52oeFOlAc6oBUAaCInhxjPUylMWXKFNfrFmGnV60XzHpaIqzXxk6LZ2RkuB7IyGK9RHZqz3rgotlp30hvUOTK3woVKrjerGi23Xqo7DR6UWzfX375pes1Kor15EYuPLIyWq+u9X5ZL2jBspWG9XKdcsopey3ffPNNXl1br6X1bkXr16+f682zlIbIBXPG0h2iRS7sSXTesfVwWW+n9YTaqf/otp88eXJeeaMVvPDFesGt59t6A4tiPdHWw2o9n9E5tnbBm/X0xetzx+OYLoylH1gPoy12QZGlyNixX9yoE6U9Vq1XPzrv1+rdRD6Lvd7y2a1H3Nqi4Gex3tDotIzSsNSEiMzMTPeedlyb8nyvIqy31OrQclrteLBe/+L+FkTYc+w7ZnUT8d5777lUDet1BZINqQJAISyX1BQXSBRm+fLle22zYMb+QzQHHHCAW7fTuEWNRhC5sKao94xcaRwJ5KK3239Cdoowcgq3IMuttNO8+zJgwACXD2g5mtF5o9EBTGktWrRon6cu7T/Pc88916UGWIqGBRQWsEUPS2X5fpZaYHmK0ew/WqvrWA9PtC+R/dnp+YIsgOjZs6dLk7Ccw6LaNBLw2rFR1I+l4vazYMECd0o/HuJxTBfGgjjLW7ZT+FZXVncluaq+NMdqwc8SCWIjn8VScexH4d133+2Woj7LDz/8oNKyfVgO6//93/+5HyDRYjEUWSQINnZRaSRt4c9//nOxrzv11FPdZ7L0E/v7EJ0rbUPC2Q8Ha5OSKs1zgbIgcAUKYcGD5bjZcEilYQFVYaw3yNh/isZyOe2CqMIUHL6mqPfc177KynpZrGw2FJEN+WW5l7Yvy6mM/o8tHqzXLpJTaIGr3dq+LZc1Fv9BFvWakvRMxVK82i4e4nFMF8Z6Wkuak1nWY7Wkn8Xey3pYC1PwB1NJWS+uDXVl721nPSxn2vZn+4nsN1YsILdeU6uffQWukV7Vwr5jxi44swvQonuNi7qQ0H6gRfcsA/FA4AoUwa6utR4ku3AiujejPOwiC+vNtUCptP9Jx4JdpLOvYNyuWrfnFbzIya7ejjfrabN6t9PEN954ozuFaadz7QrsCLvQyerPLiaxnsYI6zWyXq3ihvCJ9Gza6evo5xXWS1vSwDjyPnYlekHWC2htHt3bWlbR+yk4e5RtK+vQReXtIQvymI71sRrppbX0g319ltLUmx1vlhZjo4NEXzxpPbzled/iWHC5r55cCzRtxAjroY2MkhDNRk2xwDYSuEYfgzbKScH9WarSpEmTYlJ+oCjkuAJFePDBB12viF3BXtipTjuNbVdYl4bl5NmpessJLGxmG8tzjCfb9/HHH1/sjESRXqno3j+76thO3yeCpQs0btzYXS1uZY3OvTORSQzsyuloFuia8ePHF/nekavrTzrppLxt1ttlV3cXZD2/JTmFa1eI2+D5dsV/9POtfe00bGGTLpSF5W1aOoQNyB99StxSEexK+OI+d3EiM2qVZuikZDmmY32sWhBuPwrsB6vlLhf3WUpTb4WVs7BjuLTvG0nVKMh+iNlIF9GjcxTG0nJq1Kihxx9/3LVhwcV+RFq7Ro43C+Ytv9ryaAt+FvsO2agGlqcOxBM9rkAxvS82LI0FTpYvFj1zlp3ys15BuyCntGzIIxvz0y6msiF87OIkG8vRLtCwXpmiplOMBTtNab1UdlrQLu6aM2eO27cNh2UBkZ3SjfxnZadfLRiycTTtMSun/SdXVvb5CrvYw4LJ6B5tC/SsB8/Gs7SL0grm5FoZrd4tuLD/3D/44AMXrNjwWFbmSO9QYewzWAqCDfNjn9tygS3n0C7oKcjqxh6zi75sHE/7EWN1Uxg7HWv/Ydt729ickeGwLEfThseKBasLyze0z26f2YaAigyHZUNe2Zi+ZWGf09iwTtbzZr2NlrJRml7ioI7peByrNkWuTelsF6HZZ7G/A1bPFgzb8Gj2Y8rY6f5Im9gPFgvo7PR8YUNXWeqRtZkNv2XBnaUh2Y+awsaoLW17WDktoLTy2BkFOxNhQ4LZfvY1vJp9Hy1Fw8ZJLsy4ceNcQGoX/lkd22e78847XdlsSDB73Mplfw/tb6WlPVhZgXhjyArqgGOgmGOgWbNm3pNPPuktWbLEDQ20efNm78MPP3TDAUUPmWNGjBhR6FA/0UNS2XLAAQe459oQNNnZ2d4PP/zghuG54oor9hqiyYbdiX5tUUMdFTZUTWH7tiGbbPibFStWuM+zfPly95zooYxuueUW91obnmrOnDluyCV7TsHhoWIxHFbB8tny4osvuscmTZpU6HtWrFjRu/32273Fixe7+rN6vOeee/Yawqjg8Em2NG3a1L2vfTYbNunuu+/2unfvvtfQRtWqVfP+9a9/eRs2bHCPRT57YcNh2dKtWzd3XGzfvt3btGmT9+9//9tr2bJliYYTKmqYrsKWCy64wLWJlX/dunWurho1alSiY6So5a9//as7Hmz4ruhyxPqYLmqx97Phx/b1vMLasyTHaqTNbKipgu9Z2DFsx4gNY2afwT6L1c24ceO88847L9/zLr/8cu+///2v9/PPP+9zaCxrIxv2zI4nGyJs9OjRXsOGDQvdf1HtUdhir509e7a3fv16b/fu3d7KlSvdUGLHHHNMsXVp7WXPf/7554sdwmvbtm2u3NHb+/fv782cOdPbunWrq/dvvvnGfR8Lfv9YqAPFoQ4y/BUAAAAgqZHjCgAAgFAgcAUAAEAoELgCAAAgFAIPXG36PBu42q5stKsT7Yrh1q1b7zUmn81UYo/bFIqFjX0HAACA1BZo4GpD2cyYMcMN93H66ae7sQht6r7o+b5t+BAbK9OGOGnXrp0b486G3LC5ngEAAJBeAhuy4r777vOmT59e7HNsOJLoIUxq1qzpht/o168fQ20w1AbHAMcAxwDHAMcAxwDHgNKnDgIdDuvrr792vac2qLPNh2yDMv/jH/9wMxUZG5zZBn+2AZ/nzp2b9zobYNwGWy5s1hEbHL5gb2xkoHEAAAAkp6ysLJcamrQzZ9mUmTZ13NChQ3XvvffqhBNOcHMj7969281SFJlyz6Y5jGb3C5uOzwwaNChmM9UAAAAgcWzK7+KC10ADV5sj3OZS/utf/+ruWy/qMccc4/JZLXAtC5vK0QLh6OjdenJbtGjhpmyMN5tmcOHChQnbH2KPNgw/2jDcaL/wow3Dr0aC45nI/myK5OIEGriuXr3azSkdzeaEt7mnzZo1a9ytzRMdWY/ctyC3MNZba0th+9pXZcSCBcqJ3B9ijzYMP9ow3Gi/8KMNwy8rwfFMZH9JPaqAjShgkXy05s2ba9myZW596dKlrsK6d++e74PZ6AKzZs1KeHkBAAAQnEB7XIcNG6aZM2e6vNTXXntNbdu21ZVXXumWiEceeUS33XabFi1a5ALZu+66y+U+vPXWW0EWHQAAAAEIdFiDXr16efPmzXNDXH3zzTfeFVdcsddzhgwZ4q1evdo9Z/Lkyd4RRxxR4vfPysryjN0m4vMken8stCHHAN/DVDsG+DsafBvQhsHXYdBLVpLGT4H2uJrx48e7pTh33HGHWwAAAJC+Ap/yFQAAACgJAlcAAACEAoErAAAAQoHAFQAAAKFA4AoAAIBQIHAFAABAKBC4AgAAIBQIXAEAABAKBK4AAAAIBQJXAAAAhAKBKwAAAEKBwBUAAAChQOAKAACAUCBwBQAAQCgQuAIAACAUCFwBAAAQCgSuAAAACAUCVwAAAIQCgSsAAABCgcAVAAAAoUDgCgAAgFAgcAUAAEAoELgCAAAgFAhcAQAAEAoErgAAAAgFAlcAAACEAoErAAAAQoHAFQAAAKFA4AoAAIBQIHAFAABAKBC4AgAAIBQIXAEAABAKBK4AAAAIBQJXAAAAhAKBKwAAAEKBwBUAAAChQOAKAACAUCBwBQAAQCgQuAIAACAUCFwBAAAQCgSuAAAACAUCVwAAAIQCgSsAAABCgcAVAAAAoUDgCgAAgFAgcAUAAEAoELgCAAAgFAhcAQAAEAoErgAAAAgFAlcAAACEAoErAAAAQoHAFQAAAKFA4AoAAIBQCDRwveOOO+R5Xr7l22+/zXs8MzNTjz32mNatW6etW7fqjTfeUP369YMsMgAAANK1x/Wrr75Sw4YN85ZOnTrlPTZs2DD17t1bF1xwgbp06aJGjRrpzTffDLS8AAAACEYlBSwnJ0dr167da3vNmjV1+eWXq3///po2bZrbdumll2rBggVq166dPvnkkwBKCwBIzj6YyJJRyG1kQaJ4XpY2bbLbWsnQR4ZytWFyfXcCD1yPOOIIrVq1Srt27dKsWbM0aNAgrVixQq1bt1aVKlU0ZcqUvOcuXLhQy5YtU/v27YsMXO01lmIQkZWVle823hK9P8QebRh+tGEwPK+CPK+OPK+uW6S6/v0seV51SdXleTXcIlWzhDB5XlX/1v5u/7Ls2JGppk2l7dsXKiOjojyvsh/8VIxaKiTDf2EowrZtUp06traCOgp5GzZocIQyMr6L+/5KGjcF+q234HPgwIEuID3wwANdzuuHH36oY445xqUNZGdna/PmzfleY72z9lhRLPAdPHjwXtstOE6kRO8PsUcbhh9tGFs//SQtWCAtXiytXLn3sn597Pb1/ff274Gxe0MAZfLZZ5/poIOUNAINXCdMmJC3Pn/+fBfIWo9q3759tXPnzjK953333aehQ4fmi+DtP6/GjRu7C7ziLdH7Q+zRhuFHG5aP5+2nPXt+rdxcW1r4S3N53v6Kr2xJu5WRkaMDDqijdetWy/NsW44yMvZIil5y8y0ZGZF1z1+i16MXJEKlShXVtWs3TZv2nnJyrL0Q1jZs0+ZY7dixLGF/t/dZLiUR61397rvv1KxZM02ePNmd8q9Vq1a+XtcGDRpozZo1Rb7H7t273VKQBZGJDCQTvT/EHm0YfrRhSTWV1FHSif5yXCn/e7C/ufYfjl2vsL7AskHSFjvxaC3i39qyXdKuqGV3XmBZo0aW1q7dopo1W+T9HfWIOUOlWrUsvfuuteE5ys7m/8Jwt+GypIpnkipwrV69ug4//HC9+OKLmjNnjgtAu3fvnjeSQPPmzdWkSROXCwsAKCvLGe0s6UxJveyvawles1LSAkk2ZOEiScv9bbb8SG8mgNQPXB966CG9/fbbLj3AhroaMmSI9uzZo1deeUVbtmzRyJEj3Wn/DRs2uPsjRozQzJkzGVEAAEptP0nnSTpH0qk2dksRz7NT7F9Lsg4Cuwh2vh+wJk+PC4D0FWjgetBBB7kgdf/999dPP/2kjz76SCeeeKKbcMDccMMNys3N1ZgxY1zawMSJE3XNNdcEWWQACJkOkgZK6ldEsJojaYakyX6w+ilBKoCkFWjgeuGFFxb7uI0qcN1117kFAFBSNsPgpf7SopDHf5L0rqTxkiZJ2kTVAgiFpMpxBQCUx6GSbpJ0mZ8aEM1O9Y+W9LykmX5KAACEC4ErAITe0ZJukfR/hfxZt5kHn5M0RtKOgMoHALFB4AoAoXWspLslnVVguw039aSkxyUtDahsABB7BK4AEMoc1rskXe5PfxphF7YO9wPWjQGWDwDig8AVAEIjU9L1km4tMEKAzQf/sKRnSAcAkNIIXAEgFM6X9KA/y1XEFj9V4FF/ulQASG0ErgCQ1A6U9ESBPFab+/1pSX/zh7YCgPRA4AoAScsmDhgmqXbUNht39U+SvgqwXAAQDAJXAEg6B0t6SlLPqG1rJdnMgW8GWC4ACFaFgPcPAMjnMr83NTpofVHSUQStANIePa4AkBSy/LFXo6fCXiXpKn9qVgAAgSsABK6VPx1rs6htI/1c1s0BlgsAkgupAgAQqN9LmhUVtG72h766gqAVAAqgxxUAAmEjBTwr6dyobbMl/R/TtAJAEehxBYCEO9IPUqOD1r9L6kTQCgDFoMcVABKqt6SX/IuxzHpJl3ABFgCUAD2uAJAQGZJukzQuKmj9wr8wi1EDAKAk6HEFgLirLmmUf9FVxGh/zNYd1D8AlBA9rgAQ91mwZkQFrbmSBvkXYRG0AkBp0OMKAHFjaQD/kXRg1FBX/SW9Q50DQBkQuAJAXJwp6VU/TcAs8i/MWkh9A0AZkSoAADF3jaS3ooLWjySdSNAKAOVE4AoAMR054GFJj0uq6G+zXtdTJG2gngGgnAhcASAmMiW9LulPUdvu83Nas6ljAIgBclwBICbTt/5b0kn+/RxJv5P0DHULADFE4AoA5dJY0ruSjvXvb/OHvppIvQJAjBG4AkCZtfQD1EP8+z9KOkPSHOoUAOKAHFcAKJP2/sQCkaB1saQOBK0AEEcErgBQar0kTZFU178/xw9aLXgFAMQLgSsAlMoAf4zWav79yZJO9tMEAADxROAKACVmQ109H3V5wMt+76tdkAUAiDcCVwAokQf8yQUiHpV0saSfqT8ASBBGFQCAYtkMWE9Juixq222S7qHeACDBCFwBoEhV/Slbz/bv5/oTC1ggCwBINAJXAChULUnjombDsmlbL5I0hvoCgIAQuALAXg6UNEHSr/z7W/1e12nUFQAEiMAVAPI5wp8Nq6l/34a5Ol3S59QTAASMwBUA8rSS9K6k+v79pZJOlfRf6ggAkgDDYQGA003S+1FB61xJHQlaASCJELgCgC70e1qz/Lr4QFIXSaupGwBIIgSuANLcTf4MWFX8+2Ml9ZS0OeByAQAKInAFkMZ//oZLeihq25OSLpC0K8ByAQCKwsVZANJQpqR/STo/ahuzYQFAsiNwBZBm6kh6K2pigRxJv5U0KuByAQD2hcAVQBo5TNJ/JB3p39/mpwbYZAMAgGRH4AogTXTwe1oP8O+vldRL0pyAywUAKCkuzgKQBvpLei8qaP1GUnuCVgAIGQJXACnL8+zfOyS95F+QZSb7va82KxYAIExIFQCQknbtsuUZSX0LDHd1nX9BFgAgbOhxBZBycnMb6aSTpJycSNCaK+lGSVcTtAJAiNHjCiDFdNKOHWP06aeR+9v9HNdxwRYLAFBu9LgCSCG/cxdheV59dy8j43s/n5WgFQBSAYErgBRgF16NlPQPSZXdllNOkapXP1nSvKALBwBItcD15ptvlud5GjZsWN62zMxMPfbYY1q3bp22bt2qN954Q/Xr/9KTAgC/aCJpuqTL8iqkcuXhevdd63HdQCUBQApJisC1TZs2uuqqqzR37tx82y2I7d27ty644AJ16dJFjRo10ptvvhlYOQEkm3MlfSGprX9/h6QLVbXq7apEBj8ApJzAA9fq1avrpZde0m9/+1tt3Lgxb3vNmjV1+eWX68Ybb9S0adP0+eef69JLL1XHjh3Vrl27QMsMIBlSAx6VZD9k6/jbFvv5rK8GXDYAQLwE3ifx+OOPa/z48Zo6dapuu+22vO2tW7dWlSpVNGXKlLxtCxcu1LJly9S+fXt98sknhb6fvcZSDCKysrLy3cZboveH2KMNk1tu7mHauXOUcnOPz9tWqdKbqlr1D8rI2GItSBuGHN/B8KMNwy8roPgpqQPXfv36qVWrVjrhhBP2eqxhw4bKzs7W5s2b821fu3ate6wogwYN0uDBg/favmrVKiVSoveH2KMNk28WrJdfln73Owtef9lmv1GHD5euvPI8ZWSct9draMNwo/3CjzYMv1VJFs8EFrgedNBBGj58uHr06OEC1Fi57777NHTo0HwRvFV648aN3QVe8Zbo/SH2aMPkk5u7v7Kzhykn55y8bRUqfKeKFQfqz3/+Sn/+c/7n04bhRvuFH20YflkBxU8xCVz//ve/l7oAd999d76c1YIsFaBBgwYudzWvMJUq6aSTTtJ1112n0047zZ3yr1WrVr5eV3vNmjVrinzf3bt3u6Ugq/REBpKJ3h9ijzZMFuf4U7VGjyjygnJzr9GOHTa5QNFow3Cj/cKPNgy/rUkWz5QocL3++us1a9asQgPCwnTq1MkNY1Vc4Go5rcccc0y+bc8995wWLFigBx54QCtWrHD76969e95IAs2bN1eTJk1cWQCkutr+BVi/idq2zp9k4I0AywUASPpUgXPPPVc//fRTiZ67ZYtdIFG8bdu26euvv863bfv27Vq/fn3e9pEjR7rT/hs2bHDvOWLECM2cObPIC7MApIqz7NJNSyqK2vZvSVdK+jHAcgEAkj5wtWGoCl4kVRwbk9UuoiqvG264Qbm5uRozZoxLG5g4caKuueaacr8vgGSeTOBRP3CN2CTpD5JeDLBcAIBk4aXykpWV5Rm7TcX9sdCGqXEMVPakmz1pu/fL+AGRZYInHcT3MM0W/o4G3wa0YfB1GPSSlaTxU7lGFTj66KPdjFYVK1bUjBkz8l1oBQAl00nSE/YXJWrbD3bORdJrVCIAoPwzZ9kpe7vAygLXrl276r333tOtt95a1rcDkJYsWJ0SFbTukTRc0pEErQCAvVQqzbirK1euzLtvQ1ZZj6tdTGVOPPFEjRs3Tvfee29J3xJA2rvUn77VfOKPGPBF2tcKAKCcPa429eof/mAXSPzCAtaePXu6KVZr1KihU045pcSjDgDAL872b3Mk9SRoBQDEJnC1aVlbtGjhhqI67rjjdOWVV7qr/nfu3KlNmza56VsvueSSkr4dgLRn6QHN/Fp43x89AACAGKQK2KwJ1157rdq3b69Ro0a5nNbOnTu7C7NsKc1wWQDwy4xY0WO0AgAQ44uzbNaqNm3auFmxvvjiCzdFK0ErgNIjcAUAxKnH1XpVLT3gyCOP1Ny5c91FWKNHj9YTTzyhgQMHuou1fvyRGW0AlITNiNXGX58jaQXVBgCIXY+rTb9qwalNy2ozaQ0bNkyLFi1S9+7dNWHCBNcTe/XVV5f07QCkteiZsd4KsBwAgJQMXM8++2z16dNHgwYNUo8ePdSrV6+8x5599lk3HJblvALAvpEmAACIY+C6du1anXrqqapcubK6deuWN35rhA2FddFFF5WhCADSS21JJ/vrSyTND7g8AICUy3G1NIGXXnpJQ4cO1erVq9W3b9/4lgxAijpDUmV/nTQBAEAcAlebgKBBgwaqV6+e1q1bRx0DKOekA4bAFQAQx+GwCFoBlJ1N73p65K+JpJlUJgAgtoHrnDlzVLu25aWVzIcffqhGjRqVvBQA0kQ3SVn++tuS9gRcHgBAyqUKHH/88W6a1w0bNqikz8/MtJ4VAChqNAHSBAAAccpxnTp1qjIyMkr0XM/zSlkMAKkvI2r81h2SJgdcHgBASgauTZs2LfUbr1y5sizlAZCy2klq6K9PkrQz4PIAAFIycF2+fHn8SwIgxZEmAABI8KgCAFC+wNUuyPoPlQgAKDUCVwAJcISkFv76R5Lyz7wHAEBJELgCSNBsWRH0tgIAyobAFUCCA9d3qHEAQHyHwyqocuXKql+/vipUyB/7rlixoqxvCSAlVZfUxV9fJumbgMsDAEibwLVZs2Z69tln1aFDh3zbbYxXG7+1UqUyx8IAUna2rMiEJOMDLgsAIMxKHWWOGjVKOTk5OvPMM7V69WomGwCwD6QJAAACClxtOtfWrVtr4cKFMSoCgPQIXHdJmhZwWQAAaXVx1jfffKN69erFpzQAUswxkg7x19/3p3oFACBBgevNN9+sBx98UF26dFHdunWVlZWVbwGA/yFNAAAQYKrAlClT3O3UqVPzbefiLAB7I3AFAAQYuHbt2jWGuweQumpJ6uivW0784oDLAwBIu8B1+vTp8SkJgBTTI+pPDJMOAADKr0yDrtaqVUuXX365jjzySHf/66+/dmO7btmyJQZFApAaSBMAAAR8cZYNhbV48WLdcMMN7uIsW2688Ua37de//nWMiwcgnDIkne6vb7dzNQGXBwCQlj2uw4YN07hx4/Tb3/5We/bscdsqVqyoZ555Ro888ogbbQBAurMfsQ39dbugc3fA5QEApGXg2qZNm3xBq7F1GyLrs88+i3X5AIRSr6h18lsBAAGlClge6yGHRAYU/5+DDz5YW7dujVGxAKROfuu7AZYDAJDWgevo0aM1cuRI9e3bVwcddJBb+vXr51IFXnnllfiUEkCI2Mx6bf31eZJWBFweAEDapgrcdNNN8jxPL7zwgipV+uXlP//8s/75z3/qlltuiUcZAYTKaVG/iUkTAAAEGLhakHr99ddr0KBBOvzww902G1Fg586dMSwWgPBiGCwAQBKN42osUP3qq69iWxoAIWc9raf665slzQq4PACAtAtcx4wZo4EDB7qLr2y9OH369IlV2QCETis/x9VMlZQTcHkAAGkXuG7evNnltUZGFYisA8De+a0RE6kcAEDiA9fLLrssb/3SSy+NbQkApBACVwBAEg2HNXXqVNWqVWuv7VlZWe4xAOmqpqT2/vpCScsCLg8AQOkeuJ588smqUqXKXturVq2qzp07x6pcAEKnW9RJHNIEAAABjipw7LHH5q0fddRR2rBhQ979ihUrqmfPnlq1alXsSwggJEgTAAAkSeD65ZdfuouybHnvvfcKHR7r97//fazLByB0gWu2pPcDLgsAIK0D16ZNmyojI0NLlixR27Zt9dNPP+U9tnv3bv3444/Kzc2NVzkBJLUj7K+Ev/6RpB0BlwcAkNaB6/Lly/PSAgAgP9IEAABJeHHWLbfcUuiQWLbtL3/5S6zKBSBUIrNlGS7MAgAkSeB61VVXacGCBXtt//rrr3X11VfHqlwAQqOypK7++mpJ8wIuDwAgVZU6cG3YsKFWr7b/nPKznNcDDzwwVuUCEBodJdXw1ycFXBYAQCordeC6YsUKdexo/1HlZ9t++OGHUr2X9dDOnTvXTSlry8yZM92wWhGZmZl67LHHtG7dOm3dulVvvPGG6tevX9oiA4gr8lsBAEkauD799NN65JFHNHDgQB1yyCFusfzWYcOGucdKY+XKlS5ntnXr1mrTpo0bZuvf//63GyfW2Hv27t1bF1xwgbp06aJGjRrpzTffLG2RASQkcLVRRSZT1wCAuPJKu9x///3ejh07vJycHLds27bNu/3220v9PoUt69ev9y677DKvZs2aXnZ2ttenT5+8x1q0aOGZdu3alfj9srKy3GvsNhblS7b9sdCGwR4D9T3J85dPk+Z45HsYfBvQfsHXI20YfD2EeclK0vip1D2uxnpJDzjgAJ144ok67rjjVLduXd11113lip4rVKigfv36qXr16po1a5brhbWpZadMmZL3nIULF2rZsmVq3z4yHzqAYDGaAAAgCcdxLWj79u367LPPyl2AY445xgWqVatW1bZt23Tuuefq22+/1fHHH6/s7GyX+xpt7dq17gKxoliwa7mxEVlZWflu4y3R+0Ps0YYlt3PnmcrJ+WV9v/0+VKVKyXHc04bhRvuFH20YflkBxU9xCVytN7Rv374uv9UCxWh9+vQp1XtZL6oFqbVq1dL555+v559/3uWzltWgQYM0ePDgvbavWrVKiZTo/SH2aMPi2UR59hvSJtGzvzfr109QZRsZK4nQhuFG+4UfbRh+q5Isnil14Gqn81944QVNnDhRp556qiZNmqTmzZurQYMGGjt2bKkL8PPPP2vx4sVu/fPPP9cJJ5ygP/7xjxo9erTrObWANrrX1fazZs2aIt/vvvvu09ChQ/NF8FbpjRs3diMTxFui94fYow1LZs+e47Rjx4dufefO/2j//fsnzeFIG4Yb7Rd+tGH4ZQUUP8U8cL311lt1ww036B//+Ie2bNnigsylS5fqySefLHR817LkulrAOmfOHO3evVvdu3fPG0nAAuQmTZq41IKi2GtsKcgqPZGBZKL3h9ijDfflf8Pi5eSMT8rjnTYMN9ov/GjD8NuaZPFMqQPXww8/XOPHj3frFiDaxVSRoatsOKvCTtMX5d5779W7776r5cuXu0i7f//+Ovnkk3Xaaae5oHjkyJGu93TDhg3u/ogRI9xYr5988klpiw0grhdmMfEAACAJA9eNGzfmJdBal65dXPXVV1+pdu3aqlatWqneyyYTsLQDm3HL0gHmzZvngtbISALWs5ubm6sxY8a4XlhLT7jmmmtKW2QAMWff9U7+uqX6LKGOAQDJF7hOnz5dPXr0cMHq66+/ruHDh6tbt25u29SpU0v1XldccUWxj9uoAtddd51bACSTk2wMD3+d3lYAQJIGrhZE2tBV5p577nEXV3Xo0MH1it59993xKCOApE4TYLYsAEASBq4VK1bUmWee6U7ZG8/z9MADD8SrbACSPnDdI+m9gMsCAEgXpZo5a8+ePXriiSfyelwBpKNGko721+1CyfyThAAAEC+lnvJ19uzZbsIAAOmqR9Q6aQIAgCTOcbXxW22IqoMPPtiNtWpTv0abP39+LMsHIOkwDBYAICSB66uvvupuH3300bxtluuakZHhbitVKtMssgBCIUPSKf76FjsHE3B5AADppNRRZtOmTeNTEgAhcJyNwOyv20VZOQGXBwCQTkocuH7wwQc666yz3CxXpnfv3po8ebJ27doVz/IBSNr8VsZvBQAk6cVZnTp1UpUqkQHHpX/9619uxisA6YTxWwEAIRpVIMJyWgGkk/0kdfbXl0r6b8DlAQCkmzIHrgDSjQWtmf46w2ABAJL84qzTTjtNmzf/Mth4hQoV1L17d61duzbfc95+++3YlhBAkmAYLABAiALX559/Pt/9J598Mt99hsMCUhnTvAIAQhK4VqxYMb4lAZDEGko61l//TNLGgMsDAEhH5LgCKAGGwQIABI/AFUAJkN8KAAgegSuAUkzzulXSx9QYACAQBK4A9uFYP8fVvM80rwCAcASuNgRW586dVatWrfiVCECSIU0AABDCwDU3N1eTJk1SnTp14lciAEmGC7MAACFNFfjqq6902GGHxac0AJJMVUkn+evLJH0XcHkAAOms1IHrbbfdpocffli9evVSw4YNlZWVlW8BkGrTvFrwapjmFQAQopmzzDvvvONux40b52bKisjIyGDmLCDlkCYAAAhx4Nq1a9f4lARAEl+YlStpasBlAQCku1IHrtOnT49PSQAkGRsC67ioaV43BFweAEC6K9M4rp06ddKLL76oGTNmqFGjRm7bxRdfrI4dO8a6fAACE5l0wJDfCgAIYeB63nnnaeLEidq5c6datWqlzMxMt93Gdr311lvjUUYAgSC/FQCQAqMKXH311bryyiv1888/52233lcLZAGkWn7rNkmzAi4LAABlCFxbtGhRaJ7r5s2bVbt2beoUSMlpXv/3IxUAgNAErmvWrFGzZs0KzXtdsmRJrMoFIFCkCQAAUiBwffrppzV8+HC1bdvWjdtqF2f179/fTUrwz3/+Mz6lBBBQmoCZRO0DAMI5HNb999+vChUqaOrUqapWrZpLG8jOznaB62OPPRafUgJIoMyoaV5XSFpI7QMAwhm4mnvvvVcPPfSQSxmoUaOGvvnmG23fvj32pQMQgE6S9vPX6W0FAIQ4VWDkyJEuWLURBb799lt9+umnLmi13ld7DEDYkSYAAEiRwPWSSy7RfvtFemP+x7YNGDAgVuUCEBimeQUAhDxVICsrSxkZGW6x9V27duU9VrFiRZ1xxhn68ccf41VOAAlRX9Lx/vrnktZT7wCA8AWumzZtcqMI2PLdd9/t9bhtv+OOO2JdPgCBTfNKfisAIKSBa9euXV1v63vvvac+ffpow4YNeY/t3r1by5Yt0+rVq+NVTgAJz2+dTJ0DAMIZuEZmy2ratKlWrFjhelgBpOrEAzZKyMyAywIAQDmHw1q+fHnexViHHHKIqlSpku/x+fPnl/YtASTNNK+NoqZ53R1weQAAKGfgWq9ePT333HM6/fTTC328UqUyDQ0LIHCnRa1PDLAcAADEaDisRx55RLVr11a7du20c+dO9ezZ0w2RtWjRIp111lmlfTsASYPAFQCQ3ErdPdqtWzedffbZmjNnjnJzc91FWVOmTNGWLVs0aNAgvfPOO/EpKYA4qiaps7/+vaS9Rw4BACB0Pa7Vq1fPG69148aNOuCAA/JyW1u1ahX7EgJIgC6SMv110gQAACkSuC5cuFAtWrRw63PnztVVV12lRo0a6eqrr2Y4LCC0SBMAAKRgqsDw4cN14IEHuvUhQ4ZowoQJuuiii9xYrgMHDoxHGQEkLHDNkTSV+gYAhDtwPfTQQ/X999/rpZdeytv2+eefq0mTJmrZsqUbJmv9eqaHBMLnEEkt/fWPJW0JuDwAAJQzcF28eLG7EGvatGlu9qz3339fq1atciMLfPHFFyV9GwBJhzQBAECKBa42msDJJ5/slgsvvNBNPLBkyRIXxFowa0vkoi0AYULgCgBIscD1gw8+cIvltWZmZqpDhw55gayN41q5cmUtWLBAxxxzTHxLDCCGKkrq7q9bqs8cahcAkLTKNM1Vdna262H96KOP3K3NomWjC1iuK4AwaSeptr8+WVJuwOUBACBGgav1qp544onq2rWr62m12bNWrFih6dOn67rrrnM9sgDChDQBAEAKBq5Tp051gerSpUtdgPrkk0+qf//+WrNmTXxLCCBBgeskahoAkBoTEHTu3NkNd2UXY1kQO3ny5HIHrbfccotmz57tpotdu3atxo4dq+bNm+d7juXTPvbYY1q3bp22bt2qN954Q/Xr1y/XfgGYupJO8KtivqQfqBYAQGoErrVr19aVV16pHTt26Oabb9YPP/ygefPmacSIEerTp4/q1atX6p136dJFjz/+uEs/6NGjh0tFmDRpkqpVs3nTfzFs2DD17t1bF1xwgXu+zdL15ptvlnpfAArqHvUngGleAQDh4JVlqVGjhtezZ0/vgQce8D755BNv165d3vz588v0XpGlXr16nuncubO7X7NmTS87O9vr06dP3nNatGjhntOuXbsSvWdWVpZ7vt2Wp2wlXRK9PxbasOzHwDOe5PnLKSl1LPE9DL4NaL/g65E2DL4ewrxkJWn8VOIe14K2b9+uDRs2uGXjxo3KycnRkUceWa4IulatWu7W3tO0bt3ajRc7ZcqUvOcsXLjQTYTQvn37cu0LQCS/dYekD6kOAEDqXJyVkZGhNm3auNEEbFSBjh07qnr16m72LBsS69prr3W3ZWXv/8gjj7ghtr7++mu3rWHDhm7orc2bN+d7ruXD2mOFsUDX8mIjsrKy8t3GW6L3h9hLhzbcs6elduw4yK1XrDhD1apVsW+PUkU6tGEqo/3CjzYMv6yA4qeYBa6bNm1ygapdkGUB6g033OCmfbXZs2LBcl1t8oJOnTqV630GDRqkwYMH77XdAuxESvT+EHup3IZDh0p/+tMv6w8/3EPXX79FqSiV2zAd0H7hRxuG36ok+zua4ecM7JNdmGUB66JFi2JeCLvA6+yzz9ZJJ52k77//Pm+79ezaKAZ2YVh0r6s9x3pnbSlJj6tVeuPGjd2oBPGW6P0h9tKhDXfseEt79nRz69WqtVHFit8plaRDG6Yy2i/8aMPwywoofqpZs+Y+9xdo8u+IESO8lStXes2aNdvrscjFWeedd17etubNm3NxVhIkbafykvoX9lT3pF3+RVlLkqA8tGHQ9ZVsS+p/B1N/oQ2DbwOFrA1Lur8yTfkaK5YeYJMYWG+rRdcNGjRw2613ddeuXW5815EjR2ro0KHugi27b72zM2fO1CeffBJk0YEQs57WyFmJdwMuCwAAJRdo4HrNNde424JTxQ4cOFDPP/+8W7dc2tzcXI0ZM8alAEycODHvdQDK4vSodQJXAEB4BBq42kgC+2KjClx33XVuARDLwDVb0ntUKQAgNMo8jiuAMGop6VB/fbo/hisAAOFA4AqkFdIEAADhReAKpBUCVwBAeBG4AmmjuqST/PWlkhYEXB4AAEqHwBVIGwyDBQAINwJXIG2QJgAACDcCVyBtMAwWACDcCFyBtMAwWACA8CNwBdLCGVHrzJYFAAgnAlcgLZDfCgAIPwJXIC2Gwersr3/PMFgAgNAicAXSahisdwIuCwAAZUfgCqQ80gQAAKmBwBVIeQyDBQBIDQSuQEo7UtKh/vp0STsCLg8AAGVH4AqktDOj1slvBQCEG4ErkDaB638CLAcAAOVH4AqkrDqSOvrrCyT9N+DyAABQPgSuQEpflFXRX6e3FQAQfgSuQMoiTQAAkFoIXIGUVElST399o6QZAZcHAIDyI3AFUlIHP8fVTJCUE3B5AAAoPwJXICX1jlonvxUAkBoIXIGUzm/d4/e4AgAQfgSuQMppJqmlv265rRsCLg8AALFB4AqknF5R66QJAABSB4ErkHLIbwUApCYCVyCl1JR0kr++RNK3AZcHAIDYIXAFUsqpkir7628HXBYAAGKLwBVIKcyWBQBIXQSuQEp9nc/w17dKmh5weQAAiC0CVyBltJN0gL8+UdLugMsDAEBsEbgCKYM0AQBAaiNwBVJuGKxcSe8GXBYAAGKPwBVICYdJOtZfny3px4DLAwBA7BG4AinhnKj1sQGWAwCA+CFwBVIucH0rwHIAABA/BK5A6NlIAh399W8kfRdweQAAiA8CVyAlLsqKfJXpbQUApC4CVyD0SBMAAKQHAlcg1GpI6uGvr5L0WcDlAQAgfghcgVA7TVJVf/3fkryAywMAQPwQuAKhRpoAACB9ELgCoVUpaprXTZLeD7g8AADEF4ErEFpdJNX218dL+jng8gAAEF8ErkBonRu1zjBYAIDUR+AKhFKGpLP99WxJEwIuDwAA8UfgCoRSa0kH+etTJG0LuDwAAMQfgSsQ+tEExgZYDgAAEofAFQh1fmuupLcDLgsAAIlB4AqETnNJR/nrMyX9GHB5AABIDAJXIHQiF2UZRhMAAKQPAlcgdM6PWidwBQCkDwJXIFSaSGrrr38haXHA5QEAIE0C186dO2vcuHFatWqVPM/T2WdHnwL9xZAhQ/TDDz9ox44dmjx5spo1axZIWYHk6219PcByAACQZoFr9erVNXfuXF177bWFPv6Xv/xFf/jDH3T11VerXbt22r59uyZOnKjMzMyElxVIDn2j1glcAQDppVKQO58wYYJbinL99dfr7rvvdr2yZsCAAVq7dq3OOeccjR49OoElBZItTeBLSf8NuDwAAKRR4Fqcpk2b6sADD9SUKTYr0C+2bNmiTz75RO3bty8ycK1SpUq+HtmsrKx8t/GW6P0hfdpw9+6LlJ0dOc7HKTMzucqXTJK1DVEytF/40YbhlxVQ/BTawLVhw4bu1npYo9n9yGOFGTRokAYPHrzXdsujTaRE7w+p34bt2kmzZ/+yPn/+39S8+d+CLlLSS7Y2ROnQfuFHG4bfqiT7O5q0gWtZ3XfffRo6dGi+CN4qvXHjxtq6dWvc95/o/SE92jA392Bt3/61W69QYb7atOkYdJGSWjK2IUqO9gs/2jD8sgKKn0IbuK5Zs8bdNmjQIG89cv/LLy2/r3C7d+92S0FW6Yn8DyzR+0Oqt2HPvLXc3FeTqFzJLbnaEKVF+4UfbRh+W5Ps72jSjuO6dOlSrV69Wt27d88XjdvoArNmzQq0bEDiXRC1zmgCAID0VCno4bCix2W1C7KOO+44bdiwQStWrNAjjzyi2267TYsWLXKB7F133eXGdH3rLWYLQjo5WNKJ/vpcSYsCLg8AAGkYuLZp00bvv/9+3v1hw4a521GjRunSSy/Vgw8+6ILbp556SrVr19ZHH32knj17KjtyaTWQFph0AACAwAPXDz74QBkZGcU+54477nALkL5IEwAAIKlzXAGYgyS196tinqTvqBYAQNoicAWSGmkCAABEELgCSY00AQAAIghcgaTVWFIHf32+pIUBlwcAgGARuAJJ6/+i1hm7FQAAAlcgaV0ctf5ygOUAACA5ELgCSeloScf76x9LWhxweQAACB6BK5CULopafynAcgAAkDwIXIGkY5Ny9PfXcySNDrg8AAAkBwJXIOl0lNTEX58k6aeAywMAQHIgcAWS+qIs0gQAAIggcAWSSuWoSQe2S/p3wOUBACB5ELgCSeV0SXX99bF+8AoAAAyBK5BUSBMAAKAoBK5A0qgpqbe//qOkyQGXBwCA5ELgCiSNPpKq+uuvStoTcHkAAEguBK5A0mDSAQAAikPgCiSFRpK6+uv/lTQ74PIAAJB8CFyBpHBh1NfxXwGXBQCA5ETgCiQF0gQAANgXAlcgcK0k/dpfn+2nCgAAgIIIXIHA/TZq/ZkAywEAQHIjcAUCVV1Sf399m6RXaA8AAIpA4AoEqp8/8YD8oNWCVwAAUBgCVyBp0gSeDrAcAAAkPwJXIDDHSjrRX58r6VPaAgCAYhC4AknR2/oU7QAAwD4QuAKB2E/Sxf76Dkkv0Q4AAOwDgSsQiPMl1fHXX5O0mXYAAGAfCFyBQHBRFgAApUXgCiTckZI6++tfS5pJGwAAUAIErkDCXRG1zhBYAACUFIErkFCZkgb469mSXqT+AQAoIQJXIKHOlVTPX39D0gbqHwCAEiJwBRLqhqh10gQAACgNAlcgYU6T1DZqpqwPqHsAAEqBwBVImDui1u+k3gEAKCUCVyAhTpHU3l//StJY6h0AgFIicAUS3tt6lySPegcAoJQIXIG4O1lSJ3/9G380AQAAUFoErkDc/S1q/W5JudQ5AABlQOAKxJVN7drVX18oaTT1DQBAGRG4AnF1e9T6PfS2AgBQDgSuQNzYKAI9/PXFkl6mrgEAKAcCVyAhua3W27qHugYAoBwIXIG46CWpp7++VNKL1DMAAOVE4ArEXJsCF2FZb2sO9QwAQDkRuAIxdbik8ZKq+/dfk/QsdQwAQAwQuAIxc4CkCZLq+/fflzSAWbIAAIgRAlcgJqr7Pa3N/PtfSTpHUjb1CwBAjBC4AuVWyU8JOMG/v1LS6ZI2U7cAAMQQgStQLkdIel3SGf79TX7QasErAACIJQJXoEyOk/SqpAV+SoD8tIBz/DQBAACQloHrNddco6VLl2rnzp36+OOPdcIJkVOyQKJ1kPS2pC8l9Yv6Cm30739AkwAAkK6Ba9++fTV06FANGTJErVq10ty5czVx4kQdcIBdwQ3Eh+dJubmHSOoj6T5JkyVtkDRD0plRz1wr6S+Smkj6N80BAECcrypJajfeeKOefvppjRo1yt2/+uqr1atXL1122WV64IEHlFzqKSenvaZMkXJyTpa0I+gCJbmMMj634Hr0UiHqtmIhS6akqgVua7i2iyzbth2gGjWkHTuKO+W/TNKD/hitu8rw2QEAQEoFrpUrV1br1q11333W4/ULz/M0ZcoUtW/fvtDXVKlSRZmZFoz8Iisry90eeOCBeevxkpPTQzt3jtJ550mVKo1TnHeHOCvYfhkZq1WhwnxVrvyOKlUaq4wMmw2rLu2QhGrYL48Efe8Re7Rf+NGG4VcjwX9HI/sLdeBar149VapUSWvX2unY/7H7LVu2LPQ1gwYN0uDBg/favnDhwriVE+niQH85VdIjQRcGJcD3Ptxov/CjDcNvYYLjJwuSt27dGs7AtSysd9ZyYqPVrVtXGzZYfqI0e/ZstW3btsjXF/V4Sbdbha9atUqNGzcutuITYV+fNVHvV5rXleS5xT2nLI/RhrFpl6DasLBtyfI9TMfvYHGP8x2kDYOQjt/D2SH9v9D2+cMPPxT7nKQOXNetW6ecnBw1aNAg33a7v2bNmkJfs3v3brdEi67w3NzcYhugqMdLu922BR247uuzJur9SvO6kjy3uOeU5THaMDbtElQbFvf8oL+H6fgdLO5xvoO0YRDS8XuYG9L/C0uyn6QeVeDnn3/WnDlz1L1797xtGRkZ7v6sWbPK9J6PP/54mR4v7fZkEOuylfX9SvO6kjy3uOeU5THaMDbtElQb0n6xaZPSPJe/o6n9d7S4x/k7Go42fDzF/i8syEvmpW/fvt7OnTu9AQMGeC1btvSeeOIJb8OGDV79+vUDL1thS1ZWlmfsNuiysNCG6XoM8D0Mvg1ov+DrkTYMvh7CvGQlaTyT1KkC5rXXXnNjtt55551q2LChvvzyS/Xs2VM//vijklF2dra7OMxuEU60YfjRhuFG+4UfbRh+2Ukaz2T4ESwAAACQ1JI6xxUAAACIIHAFAABAKBC4AgAAIBQIXAEAABAKBK4AAAAIBQLXBDrooIM0bdo0ff3115o7d67OP//8RO4e5VSrVi19+umn+uKLLzR//nxdccUV1GlI7bfffvr+++/10EMPBV0UlMHSpUvd31D7Lr733nvUYcgceuihrt3s/8J58+apWrVqQRcJpdC8eXP33YssO3bs0Nlnn61ECnww2XRZGjZs6B133HFuvUGDBt7KlSu9atWqBV4ulpLVQYUKFbz99tvPrVu7LVmyxKtbty71F8Jj6O677/ZeffVV76GHHgq8LCylr4OlS5d61atXp+5Cevy8//77XqdOndx6nTp1vIoVKwZeJhaVqQ7se/jTTz8lNJahxzWB1qxZ43oJzNq1a7Vu3TrVrVs3kUVAOdg8zjt37nTrmZmZbvphWxAuzZo1U8uWLfXuu+8GXRQg7Rx11FFuOvePPvrI3d+4caP27NkTdLFQRmeddZamTp3qel0ThcC1FDp37qxx48Zp1apV8jyv0K7xa665xp3GsgDn448/1gknnFDoe7Vq1UoVK1bUypUry956SHj7WbqAzd5m7WanmdevX08rhKwNH374YQ0aNCiBpUas29Be98EHH2j27Nnq378/FRyi9jviiCO0bds29x5z5szhuxjyWKZv374aPXq0EonAtRSqV6/uekyvvfbaIhtw6NChGjJkiAtM7bkTJ050U9ZGq1Onjl544QVdeeWV5Ws9JLz9Nm/erOOPP15NmzZ1/2HWr1+fVghRG1rvwHfffadFixbRbiH+Hnbq1Elt2rRx7Xnrrbfq2GOPTeAnSG/lbb9KlSq5wMkCo/bt26tHjx465ZRTEvwp0lv1GMUyWVlZ6tChg9555x0lGrkdZagDc/bZZ+fb9vHHH3sjRozIu5+RkeHyWG+++ea8bVWqVPE++OAD7+KLL6beQ9h+0cvjjz/u9enTh3YMURvee++93vLly12OpOVlbdq0ybv99ttpwxC1YcHlwQcf9C655BLaMCTtd+KJJ3oTJkzIe/ymm25yS1DHYLov5fkOWhzz4osvJrzM9LjGSOXKldW6dWtNmTLlf78IPM/dt1+VEaNGjXJXU/7rX/+K1a6RoPaz3tUaNWq49Zo1a+qkk07SwoULqf8QtaH1zh1yyCGux/ymm27S008/rbvuuivAUqO0bWhXoEe+h9Zz1K1bN3d1OsLRfjYyi/0trV27trtGwP6OfvvttwGWGmWJZYJKEzCVEr7HFFWvXj13CsQuuopm9+1CENOxY0f169fPDf9xzjnnuG2/+c1v9NVXXwVSZpSu/Zo0aaKnnnoq76KsESNG0HYha0OEvw0bNGigsWPHunW7TsB+fHz22WeBlBelbz+7EMt+QE6fPt39HZ00aZLGjx9PVYbs72jNmjXVtm1b9enTJ+FlJHBNoBkzZrg/tAgn6yn49a9/HXQxECPPP/88dRlCdsGI5ZkjvCZMmOAWhNeWLVvUsGHDQPZNqkCM2NBWOTk5rjcgmt23YbCQ3Gi/8KMNw482DDfaL/zWhSCWIXCNERuXzob26N69e942Ow1i92fNmhWr3SBOaL/wow3DjzYMN9ov/H4OSSwT+FVtYZohwma+ssVcf/31bv3ggw92j/ft29fbuXOnN2DAAK9ly5beE0884W3YsMGrX79+4GVnof1S4RjgOxh8G9CGwdch7Rd8PYR5qR7+WCbwAoRm6dKli1eY5557Lu851157rff99997u3btckNKtG3bNvBys9B+qXIM8B0Mvg1ow+DrkPYLvh7CvHQJeSyT4a8AAAAASY0cVwAAAIQCgSsAAABCgcAVAAAAoUDgCgAAgFAgcAUAAEAoELgCAAAgFAhcAQAAEAoErgAAAAgFAlcAKIXnnntOY8eOTXidXXLJJfI8zy3Dhg1TsurSpUteOYOoJwCprVLQBQCAZGHBVnEGDx6sP/7xj8rIsEkHE2/z5s1q0aKFtm/fripVqujzzz/XjBkzdNVVV+V73gMPPKALLrhAv/rVr7Rt27aElnHmzJlq2LChhg8frszMzITuG0DqI3AFAJ8FXBH9+vXTnXfe6QLFCAsCLWgMigXWa9euzbs/YMAAzZo1S2PGjNGkSZPctnbt2umGG27QKaecEpegtWLFitqzZ0+Rj//888+ujDt37iRwBRBzpAoAgM8CrshivZuRQDGyWNBaMFVg2rRpevTRR93p+w0bNmjNmjW64oorVK1aNT377LPasmWLFi1apJ49e+ar56OPPlrvvPOOtm7d6l7zwgsvaP/99y9VW1iP6z333KORI0eqVq1aLlC08o0YMULTp09Xx44d3e2OHTu0fPly1wtq5Yq4+OKL9emnn7oyrl69Wi+99JIOOOCAvU77W9k/++wzZWdnq1OnTq4n97333nOvs3qyx1q3bs1xBCDuCFwBIAb5p+vWrVPbtm1d0PjPf/5Tr7/+ujtt3qpVK9cb+uKLL2q//fZzz7cg0wK/L774Qm3atHGBYYMGDfTaa6+Vet8WuFrga8Hz3Xff7QLNW2+9VYcddpgmTJjgemMt0LQeZAs6H3vssbzXVq5cWbfffruOO+44nXPOOTr00EM1atSovfZx//3365ZbbtGRRx6pefPmuQB35cqVOuGEE1zAao9bTysAJIIldbFQBxwDHAMcA1HHwCWXXOJt3Lhxrzp57rnnvLFjx+bdnzZtmjd9+vS8+xUqVPC2bt3qPf/883nbGjRo4Jl27dq5+3/961+9CRMm5Hvfxo0bu+ccccQRpSqPLUceeaS3Y8cOb9euXV7r1q3dtqefftp74okn8j2vY8eOXk5OjpeZmVno+9hrTfXq1d39Ll26uPtnnXVWvudt3rzZGzBgQLHHS8F6YqEOOAY4BhSDOiDHFQDKyXohI3Jzc7V+/XrNnz8/b1skL7V+/fru1no4u3bt6tIECjr88MNdakFpfPvtt65ntXbt2pozZ07ePqyn9aKLLsp7nl1UZjmqTZs21YIFC1xvsF1wZs+tU6eOKlT45STcIYcc4t4zwlIBog0dOlTPPPOMfvOb32jKlCmud3nJkiWlKjMAlAWpAgBQTgVPk9vp+sJOnUcCwxo1aujtt9/W8ccfn29p1qyZy0kti5ycHLdE2D6efPLJfO9vAartY/HixS7XdeLEiS5P1YJbO+1/7rnnutfaiAXRCl6QNmTIEJejO378eHXr1k3ffPONSzUAgHijxxUAEswuqurTp4++//77Yq/QL+8+jjrqKBekFubYY49VvXr1XO6q5asay7ctKesVfuSRR9zy8ssv69JLL9Vbb70Vs/IDQGHocQWABHv88cdVt25dvfLKKy5YtAupTj31VDcKQaRXtrxsLNcOHTq4i8UiPa1nnXWWu29slAEbJeD3v/+9Sx3o3bu3u1BrX6pWrerew0YcsJQC24f11kanFgBAvBC4AkCC2dBTNlSV5ZvaiAOWD2s9l5s2bXI5srFg72nBZfPmzfXhhx+6EQxsXNoffvjBPW6jIAwcONBNVGCn+q3n9aabbtrn+1oPsQ3bZcN3fffdd24khHfffVd33HFHTMoNAMWx6V+KnyoGAJAUQ25ZcGsXUYWBjSdrF4tF8mYBIBbocQWAkLBA0EYisHFTk5WNFWtljB7NAABihR5XAAgBGyXAJikwllJgQ24lI8uBbdy4sVu3KWejp6gFgPIicAUAAEAokCoAAACAUCBwBQAAQCgQuAIAACAUCFwBAAAQCgSuAAAACAUCVwAAAIQCgSsAAABCgcAVAAAAoUDgCgAAAIXB/wM8GnA29E2p8wAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "silicate_mass = history_3au[:, 3]\n", "water_mass = history_3au[:, 4]\n", "total_mass_sum = silicate_mass + water_mass\n", "\n", "water_fraction = 100 * water_mass / total_mass_sum\n", "\n", "plt.figure(figsize=(8,5))\n", "plt.plot(time_years, water_fraction, color='blue', lw=2)\n", "plt.xscale(\"log\")\n", "plt.xlabel(\"Time [Years]\")\n", "plt.ylabel(\"Water Fraction [%]\")\n", "plt.title(\"Chemical Evolution of the Planet at 3 AU\")\n", "plt.ylim(0, 60)\n", "plt.grid(True)\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.14", "language": "python", "name": "python3.14" }, "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.14.0" } }, "nbformat": 4, "nbformat_minor": 5 }