{ "cells": [ { "cell_type": "markdown", "id": "431b1fc9-0685-4c50-999a-2e8026be504a", "metadata": {}, "source": [ "## gwEccEvNSv2 : analytical model for eccentricity evolution in time domain for non-spinning binaries\n", "Presented by Islam et al. (https://arxiv.org/abs/2604.17868)" ] }, { "cell_type": "code", "execution_count": 1, "id": "fbd9c81b-743b-4940-b6b4-51f0517cd711", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "lal.MSUN_SI != Msun\n" ] } ], "source": [ "import sys\n", "# install gwModels in editable mode from the parent directory\n", "!{sys.executable} -m pip install -e ../ --no-deps --quiet\n", "\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import warnings\n", "warnings.filterwarnings(\"ignore\", \"Wswiglal-redir-stdio\")\n", "\n", "import gwModels\n", "gwModels.utils.set_rcparams()" ] }, { "cell_type": "markdown", "id": "ba5044c7-00b0-458d-8f74-de9741646cb6", "metadata": {}, "source": [ "### Example evolution" ] }, { "cell_type": "code", "execution_count": 2, "id": "dbef91a5-531d-452a-9fd3-c3da0c5c9d82", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "... gwModels eccentricity at t_ref=-2638.56 : 0.07996\n" ] } ], "source": [ "# load eccentric NR data from SXS:BBH:1355 and circular NR data from SXS:BBH:0180\n", "wfs = np.load('../gwModels/data/SXSBBH1355.npy', allow_pickle=True)[()]\n", "\n", "t_ecc = wfs['t_ecc'] \n", "hecc_dict = wfs['hecc_dict'] \n", "t_cir = wfs['t_cir'] \n", "hcir_dict = wfs['hcir_dict']\n", "\n", "# compute eccentricity from waveform data\n", "obj = gwModels.ecc_measures.ComputeEccentricity(t_ecc = t_ecc,\n", " h_ecc_dict = {'h_l2m2': hecc_dict['h_l2m2']},\n", " t_cir = t_cir,\n", " h_cir_dict = {'h_l2m2': hcir_dict['h_l2m2']},\n", " q=1, \n", " distance_btw_peaks=2000,\n", " fit_funcs_orders=['3PN_m1over8', '3PN_m1over8'],\n", " ecc_prefactor=2/3)\n", "\n", "# calculate eccentricity from gwEccEvNSv2 model using q\n", "et_model = gwModels.dynamics.gwEccEvNSv2(t=obj.time_xi, \n", " q=1, \n", " e0=obj.ecc_xi[0])" ] }, { "cell_type": "code", "execution_count": 3, "id": "2088dc18-3413-4083-ad4f-eb08c9411328", "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAFwCAYAAAC1lgFqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAABk4ElEQVR4nO3dd1zTd/4H8FcS9gwJiOIkAfcigHtDtL221yGjy9ohpF7bGx1GO+6u16Hh7PrdXa9Be73etdcTaHujS4m2V23VAlHbWluV4KoTwhdwMJPfH5hIJIywvozX8/Hg8ZDv9/v5ft/f75fA28+U2O12O4iIiIj6EanYARARERF1NSY4RERE1O8wwSEiIqJ+hwkOERER9TteYgcAABaLBQaDAWq1GgAgl8uRmZnZJeXMZjNMJhMAoKysDEqlEqtWreriOyAiIqLeRCL2KCqLxYL4+HiUlJRALpcDAPR6fZuJSHvKWSwW5OXluZzHbDZj7dq1yM3N7bZ7IiIiInGJnuDodDrI5XIYDAbnNkEQEBYWhtZCa085nU4HvV4PlUrlUlar1SI/P7+L74SIiIh6C9H74OTk5DibmBwcNTKOpqWOlrNarS4JkIPVau1ExERERNTbiZrgCIIAQRCa1bAAjcmK2WzuVDmdTofs7GykpqZCEAQAQFZWFnQ6XdfdBBEREfU6oiY4FoulxX0KhQJlZWWdKpecnAyDwYC8vDyEhYUhNTUVycnJ7erATERERH1XrxhF1RJHrUtnyqWkpKCgoMDZ4RgANmzY4GzOcsdms8HX1xcAIJFIXPaVlpZ2KCYiIiLqnPDw8Gbb7HY7SkpKEBUVBan0Sr2NqAlOa0lGa/1k2lvu6hFTWVlZ0Ov1MJvNKC4ubvEcJ0+eRH19vdt9oaGhLZYjIiKinjd8+HAcP34cw4YNc24TNcFRKBQA3NfUCILQYiLT3nIZGRkoKipy7lu1ahVSUlIQHx+P7OzsFpuqgoODATQ2hbWWTFH3aWhogMVigUqlgkwmEzucAYfPX3x8B+LjOxBXe59/ZWUlRo0a5fzb7SB6DY5cLm+xtkar1Xa4nMVicSZCTalUKqxZs8Yl8bmao1lKLpcjLCysrdugbtDQ0ICQkBCEhYXxF4sI+PzFx3cgPr4DcbX3+Tv2Xd2lRPQ+OGlpac2aixydiJOTkztVrqXOyHK5HPHx8R2OmYiIiHo30efB0ev1zs6/DkajEUaj0fm9IAjQarUuw8bbUy4lJQVZWVkuxwiCgPz8/HaNpHp56yF8+sNZVFXXeXRPREREJC7Ra3BUKhVyc3Oh1+uRmJgIi8UCpVLpkoBYrVYUFha6NEm1p5zBYEB2drZz1mMAUCqV7V6m4f/+7w/IHhUHv8FqTByuxPRoBaarlEgcFQZ5gE/XPAAiIiLqcqInOACg0Wig0Wha3K9SqVBeXu5xOQCdmvNG2P4WsP0tSLx9cTpqLHaMnILQmWmQSIAxkcGYoVJiWrQC06IVCA/y7fB1iIiIqGv1igSntwq/5UnUnS5G9Yn9qPnxe9jrahA6Mw12O1CY/z52njsK3+ET4DtsPMaMjMK0aEVjLU+0EoND/cQOn4iIaMBigtOKWbPn4UDZLNQ22GBvqEPDhQrnvouHduHSoV1AwfsAgDPKEdg9fDzemJgM36FjMVIZgOnRCkyLbmzaGq4IEOs2iIiIBhwmOK346z0J8A8KwZ5jAnaXlGG3xQrzsXLU1NsQcfPjqCs9hpoT3zXW8Bz/Duf3fgLfqLHwHToWR0ovwPzOC3g9YiR8h03AKPVoTFeHO/vxjFIGNBvSRkRERF2DCU4b/LxlmKlWYqZaCQCorbfh6xMCdpdYsbskEkVH1LgQ9xMAQH3FWUj9AgEADefLcH7Ph87znPEPwd6h4/DXYRMQHPcTDFaGNjZpqZSYEa1AzKAgJjxERN3IZDI5B5lotVqkpKSIHFHXslgsMBqNyM7OhkKhaHXGfrVaDavViszMTKxZs6ZDk9rqdDrk5OQgNze31WldgMaVBfR6PSwWS6txdSUmOB7y8ZIiYZQCCaMUeGAhUN9gw/6TlZdreAbhqyNWVFXXwys4HMMeevtKDc+J/bhUXIDqI3sRknADzlbV4N1PPsNblkL4Dp+IobGTMHPsUMxUNSZTqvBAJjxERF1Iq9WivLwchYWFYofSLVQqFQwGA5RKJdauXQuz2ex2II7JZIJGo4FCoYDBYOjw9YxGY7ufpUajgV6vh06n6/D1PMUEp5O8ZFJMGS7HlOFyZM5To8Fmx/enK7HbYsXukjJ8FR6O8tEzAQC22kuos/4IicwbAHDx8Feo3LkJAHBW6oXvhozGWyMmIkCdiOFjp2LG5WRnpkqJkWzSIiLqMLPZDJVKBblc3mZtQ18nl8uRlpbWbG44h44uZN1Z7lYX6E5McLqYTCrBhKhQTIgKxb1zomGz2XH43HnstpRdbtaS41xVDQBAPvs2BMRMQ/Xxb1Fz/FtUH9+Pmh+/g636As4OHYf/7DuJf7z9FmT+wRgxXoO540dgxuWEh52WiYg8M5DWFtTpdEhKSmqW4AiCAJVKJVJUPYsJTjeTSiUYHRmM0ZHBWDZzFOx2O46UXbyS8FiCcDJqDDB9Key2BtSeLYHUu3GIud1uR/mnr8N2sQJnJVJ8HanCxuET4Td8EtRTpmP2uGHOJq0oub/Id0pEA4XNZkf5xVqxw2gmLMAHUmnzmm6z2Qyj0QiLxYKsrCxnTY5er4dKpYJOp0N+fj4AOJtssrKyMGrUKJw6dQr/+te/sHr1agCNzTt6vR4AsGHDBlgsFlitVhQVFbn0f9m0aRPWrFnT5lxtgiBAr9cjPj4eRUVFUKvV0Gg0MBqNSE9Pd17LEb/BYIBcLkdubi4UCgWSkpIgCAJyc3NdruVogsrLy3Ppa1RYWNhqDZbj+QCNfXpWrVrVbL9cLm+1NsZxDsdCma31dcrOzoZKpYIgCM4Frjszf11TTHB6mEQiQXR4IKLDA3HrtBGw2+04br2EnZZS7Cwuw05LAM5U1jiPj0x7BtXHv7lcy7Mftaf/haqCf8Er5BWcugjkfPE9qo99jZjJiZg3SYWZaiVmqJSIDOE8PETUPcov1iL+WZPYYTRT9GQylG4mXdVoNNDpdDCZTC5/sNesWYO1a9cCANLT07FpU2OXgdTUVOh0OixcuBCHDh2CTCaDVqtFfn4+kpOTYTAYoNPpYLVanX+81Wo19Hq9S5+WjIyMVhd2dhyTnp6OlJQUCIKA+Ph4FBcXO5MQq9Xq0olXEASUlZU5kxBHEuWuViYlJQWbNm1ySTBaa55y3HfT9Rwd9w00LpGkVqudCYggCEhNTW12Dsf9AI39nlQqldtEz7HcUtPrmUxd93PFBEdkEokEI5QBGKEcgfTExoSnpPQCdlrKsLO4DLssfiiNVCEk4UbY7Tbn0HTviFEAgOpj3+Dc+8/j3PtAYcQo+A2fCN/hEzFm6jTMmxKDmarGhCcimDMtExE1JZfLYTabnX9gNRoNzGazc7RVQ0MDgMbOu1arFSaTCcnJyVAoFLBYLC41IVcnGBqNpsUFn5vKy8tzJkVyuRwWi8VZ8wE0Lix9dcfcpmWsVmuLTU46nQ5qtRqCIDjP3VKNUtP7bnpPjvtOSEhAVlYW7Ha7y/Nrej6LxYK8vDyXc6SmprbYFwgAcnNzkZaWBrlcDpVKhYSEhBaflaeY4PQyEokEqoggqCKCcMf0kbDb7Th89jx2Wcqw09KY8FgvJzcA4DNkNBRLHkD1sW9Qc/xbVJk/QJX5A1TujsWJ5S/hH7uPoa7sONRDB2FB3GjMUodjpkqJ0ABv8W6SiKiXuDo5KCwsdJswqFQqZw2Ou3JyuRxqtdrj62s0mma1Kk3P7UgiHE1Njn5ETZOgljhqTnJycpCZmQmz2dxic1Fb9+1IklpjMpkgl8tdamGKi4tbTPRSUlJgNBoRFhYGjUaD9PT0Zk1incEEp5eTSCSIjQxG7OU+PDabHQfPVjU2ZxWXYXeJNyqCr0Xw1Gtht9tRX34S1ce/dfbjAQDrto04aSnCTsUw+I2YBP+Rk6GZPhsL42IxWx2OhFFh8POWiXiXRETiuPqPdk+PMNLpdFi7di02bNiA7OxstzUdjuYzjUaDhIQE6HQ65OXlITk5uc0RYTqdDkajEZmZma0mKO2577ZGQTk6MDeNqa348vPznbVHjnvvqiSHCU4fI5VKMHZwCMYODsE9sxtHaX13qhK7LGXYZWlMeKoUQ13KBE1YBFmAHNXHvsH5vR/j/N6PsfnfQNENj+HP4+fDx0uKyRHeWDh5FGaplZg0NBReMqlId0hEvV1YgA+Knux9Q63DAnw6fY7k5GRnv5ymLBYL0tPTO31+dzZs2ACTyYSUlJQW+9Lo9XpotVpnoqLVatvVIdfRxOVIiFrS1n23p8lNo9G4PUdLtT/Z2dnIzMx0LpydmZmJpKQkJjjUSCqVYOLQUEwcGooVc1VosNmx/2TF5Q7LZSgosQLj5yNw/HwAQJ1wGtVHv0bNsa/hM2Q0AOBi+Vm8+9w9+G+kCn4jJiMsNg7z583DwokjMDsmnLMsE5ELqVTitjNvf6DRaJCcnAyTyYSFCxcCaOyfAqDV0UAdrfkpKipCZmZmq+d2NDU5ZgB2JEEtXbPpTMGOeX+u7mx8tab37UiErr7vzMxMZ1LiuL7ZbHbGkZycjISEhGYjtxxNZFcTBMHlfE3vrSswwelnZFIJJg+TY/IwOXTz1ahrsOHbHyvwZXEZviwuRcERKbzlgxE8ZbGzjL2uBgFjZqH62DeoLHgflQXv42/vSLFp2HhE3rYWkSF+mB0TjllqJWbHhHNIOhH1KWazGWvXroXFYoFer0d6ejqsVisMBoNz6HjT2pPc3Fzo9XocPnwYZ8+eRVVVlXM0VNNzZWVlYdWqVcjKynLO6KtSqaBQKLB27VrnEPDWlkLQarWQSCTOoddyuRw6na5ZQpCenu5SA6PT6ZolLI5h5NnZ2c5rO4bBOziSCpPJBIVCgaysLGetkOO+HTU1xcXFLqPAjEYjsrKykJeX57wfR62NI5HKz8+HXq+H1Wp1Nmk5+v80fQeO4e7AldFUFosFGzZsaPd7bYvE3rRLNAEAKisrERoaCqvVirCwMLHD6VLVdQ0wHy3HjsOl+KK4DN+cEGC7/BNgt9tQd+4oqo99jeqjXwMABi19CgBQsSsXlyxF8BsxGarJ07B44RzMHxuFmWol5F1QLXy1hoYGHDp0CLGxsZDJ2D+op/H5i4/vQHzd/Q4EQUBGRgY2bNgAuVzunAtGr9cjNTW1y+aD6ava+/wdf7MrKioQEhLi3M4anAHGz1uGWTHhmBUTDgCouFSH3ZYyfFlchh2HS3FYEg2fQdEISbjRpZyt+gJqTx9GzfFvseeLf2DvBl+8MnQcQqbdiGlzkzErRonZ6nAkjlLA34e/jImI2pKdnQ2dTuesyXCMmDIYDC0Oq6b2Y4IzwIX6e2PxhMFYPGEwAOBMZTW+LC7FF4fL8MXhUpyqqAYAhC24G/K5d6Lm1CHUHPsa1cf2oebH79BwKQnf/FiBr08IeO6xB+A/WI3E2fNw/YIZmDd6EMYPCXE7sygR0UCn0WjcrsRtMpmg1WpFiqr/YIJDLiJD/HBz3DDcHDfMOengF8Vl+PJwKb4sLkPFsHHwGzYOobPSYa+/MlV7Q1UZLhzYjgv7P8XHWzdi8/pQ+I2cgvDR8bj+ljQsGD8Uc2MjMDiUMywTEQFXhlBfvTxCW8sbUPswwaEWNZ10cNmMkWiw2fHdyUp8UVyKLw6XouCIFdV1NgCAV0g4hv/yn6g58R2qj+xF9dF9uHjgcxwvLsAHsfPx4XdlqD13FBENpbh2cTKWxMdgerQCAT78ESSigas9c9lQx/CvC7WbTCrBpGGhmDQsFPfPV6OmvgHmo8LlJq1S7DtRAWm0Bv7RjVN3N1ysQF3ZcUhkjT9mF77dilNfvYev3/wtXhocg8DoqZg6fQ5+es0iJE8czuYsIiLqMkxwqMN8vWSYqW5czfyRxWNQWV2HncVl2HGoFNsPncMRALKAUOfxwZrr4BU25HINz9co35mLT3fmonD7DXg5WQdFoA8mBlfjJzMnYU5shGj3RUREfR8THOoyIX7eWDJhMJZc7rB8rOwith8+hx2HGmt4KhGJ4KmOZSVsqD1jQfXRvfAdOgEAUFZ1CW8/fyfehh1+IyYjalw8rr32Gtw0Lx4z1Eo2ZxFRr2CxWGA0GpGdnQ2FQuEyz0xxcTFMJpNzDafeco2rz9d0MkCgcaXw7OxsJCQkwGAwOBfRzMrKch5TVlbmnBW5K9eM6i6cB8eN/jwPjljqG2z4+scKZ+2O+ZiABpvrj56tthoVOzeh+uhe1J46DKBxvyxkEIYv/z2mT4zFvNERmD86AuOHhHB25W7EOVjEx3cgvrbeQXx8PBISEpoN6RYEAampqZ1KcLrrGllZWVi7di0yMzOdK5I7OCbgc9DpdNDpdC4rhqempgKAy4rh7b0ucGWW5fYMg+c8ONQneMmk0IwIg2ZEGH6eFIsqR3PW4VJsP1SKktILkPr4IWz+cgDL0XCpCjXHvsGlo3tRe6YYNv8w7LJYse2TD1CxMweK0QmYtygZ6dclYdH4qG6ZbJCIqDUtLT7pWCuqN17DMWOxVquFTqdzWRrh6tXQc3JymiUiGzZsQEZGhkfXdJc4abXaLkkAW8MEh0QR7Oc6/85x68XLyU5jk1YlghEwZhYCxsxyKddwsRL1whmc/vyfyPn8n8h93h/+IyZh4vwbcPut6VgwJgKThoayszIR9ThBEGC1Wp1rR7W0yKTY13CM3EpNTXVZisHdtRzD1h3kcjkSExM9itexXpUjTp1Oh/j4+Gbn7mpcMpp6heGKANw2bQRevSMee369GO/dPwN3xSmQOCoMXk2SleCp12DYQ29h8F0vInTunfAZFI2LxYX4fv/XeMl0EDf+6QuMvOFBXP/IC3jr8+9Qer5GxLsiot5CEATodDrn7MFZWVkwmUxITU1FXl4e1Go11Go1TCYTsrOzoVarMW3aNBw/fhyCICA+Ph5qtdq5AKU7FovFZeHJpolHdna2y1dTre1rzzXy8vIQHx+PsLAwmEwm53FqtRpardbtKuC5ubkwm82tXk+j0UCr1TrP6eDof9Pe6xYWFrrE0NZioV2FNTjU68ikEkwZLkdAdRhiY2Nxqd6OncVl+PzgOXx28CyOWy/Bd8ho+A4ZDcy6Fbaai7A31AEA7PV1+HHL6zjxUQ0+fGkVfKPGYuTkGVhyzRLccd1CxI9SwkvGvJ6oq0yYMMHt9n//+9+IiYnB4cOHceONN7o9Zv/+/QCAzZs34+GHH262X61W4z//+Q8A4E9/+hNeffXVZscsXrwYL730UptxZmRkID09HSkpKc6Epbi42DkHjdVqdZlVWBAEnDt3DsOHD4dcLseaNWug0Wia1TgUFhYiKysLZWVlyMvLc9s3xbHf0UyTl5fnXHG7tX3tvUZKSopzkU5H/I5FNlvqDCyXy2EwGKDX65GWlua2FsjRlOVoCktOToZer3deo73XLS8vdzmvIxnqztobgAkO9QFBvl7Qjo+Ednwk7HY7jpRdxGc/nMX/Dp7DzuIy1CDAeazEyxtD7v0jqkvMuHRkD6qP7sPBj7/DwY//gnf3b4ByyHBMH+qPuME+uGXeVM6sTDRA5OXlOZMIuVwOi8Xi0kSSlpbmMlIJAN577z3cd999AOBsFrpaQkJCqyOKHKt6N/0jv2nTJiQmJra6z5NrAI3Jh9VqhdlsdnYKbqvpatWqVTAajcjIyHCbmKlUKueIrfz8fOcSErm5uc4ErCPXXbt2LYxGY5c3312NCQ71KRKJBNHhgYgOj8Y9s6NRXdeA3SVW/O+Hc/jfwbMoPncB3mFD4B12HYI118HeUI+ak9+j5uT38A4bgsrqerz/r/ex8aOX8QvFMERNmI6FyclYdvO1mDtuOHy8WLtD5AlHLUxLYmJi2jxmyZIlbR7zwAMP4IEHHvA4PgdHf5Wmru5botFonLUnjj++x48fR2xsbLuuodPpXK5hNpthtVohl8td/pg7kgmTydTiPk+u4UgsMjMzYTQaYTQaYTKZkJaW1mbMubm5iI+Pb7Xprelsy3q9HhkZGS41TJ5cV6/XIz09vUdWSudvc+rT/LxlmD86Ar++YTy2PrIA21ctxLM3TYR2fCQCfWSQyLzgN3wiQqdf+TB6hQ1F4KRk2Gsv4tj2d/Hmb1YiOS4GQ+el4WdvFyG38DjOVbHvDlF/otPpsHbtWgiCgKysLLfDlNPT07Fp0yZYLBYkJCQgIyMDmzdvhtlsbtdyCo6Ovw6FhYWt9jPpSB8Ud9dw0Ol0yMnJAdDYF6Y9NSQajQaZmZnO4d9NY8vLy2t2vMFggCAILrG397qOvk49NYcOa3CoXxmuCMCdM0bizhkjUVtvQ+FRK/538Bz+98M5fH+6CgDgd3nBULvdjrrSY6g+sgeXLEVo8A/DR9+cxkffnMa5//4eg5QKJC9egvvSrse0mCiOzCLq4zZs2ACTyYSUlBS3zU0pKSnQ6/XQarXIzMxEcHAwkpKSEBMTg/vvv9+ja5nNZqhUKqhUKreJjCAIbmuVHPvak5w4ruGgUqmgUCiQl5fX4vBydwwGA6Kjo2EwGFxqVgoKCtwu+qlSqVzia891Hf1uHOdvOhqsu7AGh/otHy8pZqnDsebacfjkl/Owa00SspZOxnWThiDYzwsSiQQ+ESMRkngTItOfQUjiTQAAe0M9ak8dwpHt72PjU/dj5vhoyEcnYOHyR5Gz4ztUVdeJe2NE5LGioiLI5fIWkxvgSu2IYzK6tkb7WK3WFq+n1+udCY6jM7GDIAjIyclpdZ8n12hKp9MhIyOjxRqnq2cwBhqb5zZs2NDsWtnZ2c1GUDkSxKu1dl2z2exsSnP0fXLMqNydWINDA8bgUD+kJQ5HWuJw1DfYsPe4gP8dPIfPfjiHb36scB4nkXkhKsOIeuuPuGQpxCVLEaqOfI3Pis04HDIVT3x8BBPkDRhhP4vM9BswYUQEZ1Um6uW0Wi0kEgnkcjkUCoVz9M/VfUHS09Nd/kinpaVh6dKlLsc4lj1wzO/SNEFxdMpt2oE5NzcXer0eWVlZUKlUsFqtzuu2tM/TazhkZmaiuLi4WQ2QxWKBXq9HXl4eBEGAwWBwOSYlJQWbNm1yKePolN302k23t+e6giAgKSnJ2aG6qe5uquJSDW5wqQbx9fQ09Wcrq/HZD+ew7fuz2H7oHC7UNrjst9VeQu3pQ/AbMRkAUFn4H5RvzQZk3ghTT8G0uUm4M+VGpCQlws+77/+/gcsEiI/voOsIgoCMjAxs2LABcrncOYGdXq9Hampqix1e+Q7ExaUaiLrAoJArtTu19TYUHLFi2/dn8en3Z2EpvQCpj78zuQEA/5hpgN2GS5YilBfvxeaDhdj8ugEPTl2M1F89h0VjB2F+bDiiwgJauSoR9QTH5H6O2gXHiCmDwdCuNZGob2KCQ3QVHy8pZseEY3ZMOJ66fjyOlF5oTHZ+OIvdFitqG2zwlg+Gd+JNCEm8CbbaS6g+9jUuWYrgPWw88r87g/zvzuDU3x5BqFyOeUlLkHlnCrTTJrIpi0gEGo3GZRI/B8e8LtQ/McEhasOo8EDcOyca986JxoWaeuw4XIpPLyc8ZyprIPXxR0DMdATETHeWsdfXQerjj3M/FOHdA7vx7h9/h4DB0Zg6exFWPfoIrkkYDV8vVnkT9QRHYuPo5wLA2X/FXYdZ6h+Y4BB5INDXC0smDMaSCYNht9ux/2QlPv3+LLb9cBZ7jwtw9GiTeHkj8tZnYau5iEtH9uDS4QJcshTgy/ffxEPDlyD4o6OYElqDIXUn8dBdKRgVFSHujRH1c00nq6OBgQkOUQdJJBJMHBqKiUND8VBSLMrO1+B/Bxs7Kn9+8Bwqq+sh9Q1A4JjZCBwzG3ZbA+qsP0LqG4ALtQ346F+5qNjxNl584iFExE7F/OQluH9ZKpKmTxX71oiI+jwmOERdRBnki1s0w3CLZhjqGmwoOloO03dnkH/gDI6WXYREKoNP+Ajn8UGTtZAFynHp8FcoLd6HvB+KkPen5zFszlKsXPMMksdFIm54KLzZlEVE5DEmOETdwFsmxQyVEjNUSjxx3TgUnzuP/O/OwnTgDMzHymG3A17B4Qieei2Cp14LW101qo9+jUvFX6EmcjyyP7fA+L9inPv7rzB8lAo33HADHlqeipFRkWLfGhFRn8AEh6ibSSQSxAwKRsygYKxcoEbp+Rps+/4sTN+dwfZDpbhU1wCptx8CYqYhIGaas5yt5gIaIMXBnVvwws4teOGJXyBCPRGLlvwEax9/BNFR4SLeFRFR78YEh6iHhQf5Ii1hONIShqO6rgFfFpci/7uz2HrgDM42WeRT5heEIXe9iIbz5Y0zKhcXoLTEjJwNh7DLfzo00RHQyGswTm7HzdcshFTKlVeIiByY4BCJyM9bhkVjI7FobCRston49mTF5X47Z3HgVCUAQBYUhqDJWgRN1sJeX4u6suOAzBvmYwK2/vV1VBa8D59gBeJmJ2HZrUtxT+oNCAjgBINENLDxv3xEvYRUKsHkYXI8vHgMPv7FXOzQL8TTP52AubHh8JY1ThAo8fKBT6TaWSZwwkKEzEiF3ScIuz/JxYN334pguRLXr/w1dlnK0GDjSixENDCxBoeolxoWFoDls0Zh+axRqKyuw/9+OIct353Bp9+fxfmaegCAT6QKPpEqhM1fjjrrj7h0eDcuHtqNglIZbs3ehTB/Gc5/aEDyooX4+T23YdyYWJHvioioZzDBIeoDQvy8ccOUKNwwJQo19Q34srgMW/afRv53Z1B6vhYA4K0YCu9ptyBk2i3OcmdPHMGpPTtxuOAzvGb4DRTD1Ei65jqsXH4rFsyewaUjiKjfYoJD1Mf4esmwcMwgLBwzCM/eZIf5WDk2f3sam787jePWSy7HeiuHY9jP30b1kX24eGgnhMNfIXfj/+G9f7yJ2175GD+ZPBTThvhg5JBweHnx1wER9R/8jUbUh8mkEiSOUiBxlAJPXDcO35+uwub9p7F5/xlnJ2Wptx8CYqcjIHY67LYG1Jz8Hg1VVmw/3Phl3fxH1Bz6EtMXLEHmXbci9cbrmOwQUZ/H32JE/YREIsG4ISEYNyQEv0wejWNlF7Hlu9PYvP80Co82Ti4okcrgN2yCSzkv5QhUn/wB2z/MxfYPc3GPXwA0sxfh/nvvQmws++wQUd/EBIeonxqhDMCKuSqsmKvCuaoabD1wBpv3n8YXh8tQ22BzHheS8FOEJPwUdeWncPHgl7h48Et8tfUDHB08Fx9U7IJ2jBK1B7/AsvRbEBYWJuIdERG1HxMcogEgItgXt04bgVunjUBVdR0+++EcPtl/GtsOnMWlugYAgHfYEIROX4rQ6UtRX1kKWVAYio4K2LHNhLN5v8UvH9BhXPxM3HlrGu65PRWDBw8W+a6IiFrGeXCIBpjgyyOy/nS7BuantHjtTg1+OiUKgT5XFvX0CgmHRNr4vU/UGCiv/Tn8RsXhu6JdePyRn2NIVBSmJd+II6UXxLoNIqJWsQaHaADz95HhmolDcM3EIaiua8D2Q6X4+JtT2PLdGedcOzL/YARNXoygyYthq7mAS8WFuHjwSxy86I8F6z/D+CEhgDkHk0ZEQHfPnVCpVCLfFRERExwiuszPWwbt+Ehox0fiYk0d8rZ/ja+tMuQfOIuKS3UAAKlvIALHz0fg+PnOcvtPCvjx/U34uOocsp79DUaOmYQ7bktHxvI7MGrUKJHuhogGOjZREVEzvl5STBsWCMPSSSh8Mhlv3jsNtyYOR1iAd7NjJRIphmZmI2LprxE4YSGOlRTj+d8+iejoaDzzVj5OV1TDbueSEUTUs1iDQ0St8pZJMX90BOaPjsCzN03E7hIrPvrmFDbvP+2cRVni5Y2AmGkIiJkGe30dLh0xo/ro19j4TQ3+sn8rBpftxdkvcrH8jtuw/I5bMXLkSJHvioj6O9bgEFG7ecmkmB0TjudunoTdjyfjn5kzcNfMkQgP8nEe05jsTIciKQMSiQR2O/BDyXGUHD6I3zyxGqNGjULsRA2eW5eF48ePi3g3RNSfsQaHiDpEJpVghkqJGSolfnPDBOy2lOG/X5/Ex9+ehnCxzuXYkPgbEDxlCS6VmHHh++0oPvQVnlyjR+7Og/j1b34L7fhI1F2sgkKhEOluiKi/YYJDRJ0mk0owKyYcs2LC8bsbJ2LH4VJ8sO8Utuw/jarLo7EkXj4IiJ2BgNgZsNXVoLrEjHODovFwzj7IGqpx9JU7EDt+EjKW34k7br8VkZGRIt8VEfVlbKIioi7lLZNi4ZhBeCFtCgqeTEb2snjcMCUK/t5X5tmRevsiYPRMeMsbJwusuVAF35jp+OG7b/HIw7/EkCFRSJyzABtffx02m62lSxERtahX1OBYLBYYDAao1WoAgFwuR2ZmZpeVs1gsMBqNUCqVKCsrQ2JiIlJSUrr2JoioGT9vGRZPGIzFEwbjYm09tn1/Fh/sO4VtP5xFbf2VxMUrJAIRP10FW201Lh3ejQsHPkfhzi+wv/g4TkTMwE1xQxFmq8CgQYPg7+8v4h0RUV8heoJjsVgQHx+PkpISyOVyAIBer0dWVhZWrVrV6XImkwlGoxG5ubkAAEEQkJSUxASHqIcF+Hjh+slRuH5yFKqq62A6cAb/3XcK2w+dQ11D4zByqY+fc56dhurzaKg8h79+eQR//fIIqv71O1w49i2uvf6nyLh7GbRaLby9mw9bJyICAIld5AkqdDod5HI5DAaDc5sgCAgLC2t17oz2lBMEAdHR0S5JkMlkQmpqKsrLy1s8d2VlJUJDQ2G1Wrm4oEgaGhpw6NAhxMbGQiaTtV2AulRPPv+Ki3XYvP80/rPvJL4sLoWthY99lfkDnP92G2pPHQQABIWG4ZalS/H0U0/0ywkF+RkQH9+BuNr7/B1/sysqKhASEuLcLnofnJycHGcTk0PTZKQz5dauXYuEhATndgBITk5uNbkhop4VGuCNtMTheGvFdOxak4QnrxuHSUNDmx0XrLkeQ+56EVGZGyCfuww13sH421824leb9uI98wmcOFMKs9nMSQWJCIDICY4gCBAEwe3aNXK5HGazuVPl8vLyoNVqATQmPS2dj4h6h0EhflgxV4X/PjQHWx+Zj58visEIRYDLMd5hQxA6Kx1R972KqAwj9pR74+GcfUi452nEx8djpHo0fvfMszh69KhId0FEvYGofXAsFkuL+xQKBcrKyjpVznFcdnY20tLSYLFYoNVqYTAYoNFo2owvMjISEonEZdvFixfbLEed19DQAJvNhoaGBrFDGZB6w/MfpfDHL5Ji8PNFauw9XoF/7zuJD785DeuFWucx3oqhzn9LB49BsOZ6/Hjgc/zm10/hN79+CvHTZ+HJ1Y/hhhtuEOMWOqU3vIOBju9AXO6ef0BAQLPjWqq1Fb2TcWsEQehwOUdyk5+fj/z8fACARqOBXq9HUlJSu5qp3D20Q4cOdSgm8ozNZoPVasXhw4chlYrekjrg9LbnHwTgjrHeuHX0MOw5eQnbLFX48tgFVNdf+Yz6DIqGQns/whataJxQcP+nKCrchcde34L9DYOxSBWEI99/g4kTJ/aJzsm97R0MRHwH4nL3/N39Xe6VCU7TvjFXs1qtXVLu6pqa5ORkCIKA7OzsNoeinzp1ip2MRdLQ0IDDhw8jJiaGnftE0Juf/9gxwG0LgYu19TAdOIv/7DuF7YdKUX+5d7JE5uVcF8tWcwE1kGBDQRmyPz2IY6/ejZDQMNx+azruumsZEhMTm9XS9ha9+R0MFHwH4nL3/C9dutTsuMrKSrezoIua4DgCcldTIwhCi4lMe8o5jrm6I7JDUVFRm/HJZDL+UItIKpXyHYiotz//YH8ZbtYMx82a4Sg7X4OPvjmF9/b8iD3HBOcxUt9A579tMm/I5y7D+W+34bXX/ozXXvszRqjUyLjnbjzxxBO9MtHp7e9gIOA7EFd7nn9L+0SvwZHL5S3W1jg6CHeknOOYlpq5Wkp8iKjvUQb5YtnMUVg2cxSKz53H++Yf8Z75BE5WVDuPkfmHIHRmGkJmpKL29GFc2L8NJw58jhfefA+hM9Nxs2YoKk6WYOjQoa3WEhNR3yB6H5y0tDQUFxe7bHP0n0lOTu5UubS0NBQUFLgc40h4Wjs3EfVd6oggPLpkDB7WjsYuSxnyzCfwybencbG2saOiRCKB75BY+A6JRdjC+9BwUcAL+QfxQv5BVG16DFWnLLjpppuRcd+9WLRoEf/nTtRHid5rSq/XIy8vz2Wb0WiE0Wh0fi8IArRarcsw7/aUMxgMMJvNLqOu9Ho9UlJS2jWKioj6LunlBUBfTJuKgieS8ULqFMxSK9G0JUoi84JXcDiAxo6KknFaSCPUyN30TyxevBhDR4zEk08+idLSUpHugog6SvQaHJVKhdzcXOj1eiQmJsJisUCpVLp0ALZarSgsLHRpkmpPOblcjqKiIuj1emeVs1qtdkmCiKj/C/T1wtL4YVgaPww/CpfwvvkE3jX/iJLSC85jJBIJgiYvRtDkxaiz/ojz325F6bfb8Pw6A5TTb8KyhcHwbqiGTCZDUFCQiHdDRO0h+lINvRGXahAfp0gX10B4/na7HXuOC3i36AT+u+8kKqvrmx9ja0DduaPwiVTBSyqB4tAH+OajN5GWmor77r0Xc+fO7bbOyQPhHfR2fAfi6vNLNRARiUEikUAzIgzP3TwJXz2RjD/drsHCMRGQNm3CksrgE9k4Y3q9zY6S2hDYQofhb2++ifnz5yNaHYNnn30WP/74o0h3QUQtEb2JiohIbH7eMlw3eQiumzwEpyou4d2iE8gpPIFjVteZywPHzUXguLmoPXcUF77diuP7t+Gpp57COYRg3aoH4O8jQ11dXZ+YSJCov2MNDhFRE0NC/fHgolh89ugCvJMxAzfHDYWvl+uvSp+IkQhbeC+GrfwrIpY+hffLh2LacyZkvPI+BkUOxsMPP4wDBw6IdAdEBDDBISJySyqVYKZaiZfSp+KrJ5Lx7E0TMXmY6yrnjbMmT4fU2w9VNfX4767vcV7ij5deegnjx4/HrFmz8eabb3INOyIRsImKiKgNof7euHPGSNw5YyQOnKpETuFxvL/nRwgX61yO8xs5GVErXkPNif04v28zdn21Azt3foncf3+ID97LESl6ooGJNThERB4YNyQEv7lhAnY/noQ/3h6HubHhrnPrSCTwGz4R4dc/gqEP/h1hyToUBiTipj99gU0Fx7A0JRUbNmxAVVWVeDdBNAAwwSEi6gBfLxmunxyFv983HTv0i/CLpFgMDvFzOUbmF4SQ+BvgN2w89h4X8PCGT/D+f/6LzMxMDB4ShczMTBQUFLS4GjIRdRwTHCKiThoq98evtKOxQ78QG+5KwKKxg+BuehxvxVAMe+DvUCx5APVBkdiwYQOmTZuGufPmMckh6mIe9cHZuHEjVqxY0V2xEBH1aV4yKbTjI6EdH4kfhUvY9NUxbCo8jjOVNc5jpL4BCJ56LYKnXoua04dxft9m/GALxDMfHMCymSPx9RcmxMTEYOzYsSLeCVHf59FMxkqlEkeOHEFwcHB3xiQ6zmQsPs4gKi4+/65T32DD1u/P4h+7j+HzQ+fQ2m9ce30dzhjvQc15AfPmzcONN96IlStXwt/fv+cCJid+DsTVozMZl5eXQ6PR4L333ut4xEREA4iXTIolEwbjzXun4fPHFuKBhWpEBPu6P1jmBfmNjyNw/AJs/2InHnnkEYwcNQpPPvkkzp8/37OBE/VxHiU4crkchYWFiI6OxurVq7F+/XpUVlZ2V2xERP3KcEUAHlsyFl+uXoTX7tRgbmy4y36JRAK/YeMRfsOjGLryDcjnL0d5jQSv/HkDDpZWA2j836rNZhMjfKI+xaM+OI7VvOPi4hAXF4eKigps2rQJxcXFWLx4MRYtWtQtQRIR9SfeMimumTgE10wcgqNlF/D27mPYVHAcFZeuzKsjC5QjdEYqQqbdgvqKs7jltd2YOlyOik9exrEDe/CzlStxzz33QKFQiHgnRL1Xp0ZRhYaGIiMjA+vWrUNRURGUSiUef/xx1uoQEbXTSGUgHv/JOOxak4SspZMxISrEZb9EKoN32BAAwN7jAvaU++DYqVI8+uijGDpsGDIzM/Htt9+KETpRr9bpYeIbN25EbGws9Ho97HY7bDYbNm3ahDVr1mDbtm1dESMRUb/n7yNDWuJwfPDQHLy7chZunDIEXm5+Q8vn3I6olX+F8vpHYJcPx4YNGzBp0iTs27ev54Mm6sU8aqJKTExEQUEB9u7dC6PRiOzsbNjtdiQnJ2PdunVYunSpy/Fbt27F+vXrkZmZ6dKzmYiI3JNIJIgfGYapw0Jw61gfFJT54J2C4zhVUX3lGC9vBE1YiKAJC1Fz8gdcPLQLT/yvEvfZTiCw7Ad8+/U+3HvvvZDL5eLdCJHIPBomLpVKERYWBkEQEBoaiszMTOh0OkRHR7dabv369Xj00Uc7HWxP4TBx8XF4prj4/MXX9B3YIYHpwFn8fdcRfHG4rNVylR8YUL5/OwIDA7F8+XI89NBDnFOng/g5EFePDhMHgOjoaOTk5MBqtWLdunVtJjdERNQ5XjIprpk4GG+vmAHTw/OwfOZIBPq4/4UffM2voPzJr1AfNBivvvoqxo0bhyVLlrCfDg04HiU4Go0GhYWFzZqiWpKWlgalUomystb/x0FERO0TMygYT984ETsfT8JT14/HcIXrJIASLx8ETUpC5PKXEXmHAQFj5sC0dRv2nrwAm82O6upqXLhwQaToiXqOR31wdDqdRydfs2YNEhMT8dhjj3lUjoiIWhfi54375kTj7lmjYDpwBn/ZUYLdJVbn/sY5dSbAb9gENFyqxJPbSvGXr/+HEWe/wH82voCVK1fiwQcfxJAhQ0S8C6Lu41ENTnp6erNtJSUlzkn/rh41FRcXx+SGiKgbyaQSLJkwGJt0M/HBQ3Nwi2YovGWuK33K/Bv7JVhKL+DDg+dxQRqA559/HqNGjcK9996L/fv3ixE6UbfyKMHJyMhoti06Ohrr1q3Do48+CrvdjvXr13dZcERE1H4Th4bixbSp+GL1Ivw8KRbKQJ9mxwSOnYPB976KiJufgCwyFm+88QYmTpyITZs2iRAxUffxKMFpa8BVUlIS+9sQEYlsULAfHtaOxherFyErZTLGDnZdIFkilSFg9EwMut2AwcteQOCEBfikPBLmY+U4ffo0/va3v6G2tlak6Im6hkcJjkQiaXX/3r17YTKZOhUQERF1DT9vGdIShuPjX8zFPzKmY+GYiGbH+EaNQfj1j+LTIxdwy6tfYuG9q7F8+XJER0fDYDBAEISeD5yoC7TYyXjPnj0wGo0oLCxESUkJAEAQBCiVSrfHOz4ERqOx66MkIqIOk0gkmKUOxyx1OA6eqUL25xb8e++PqGtoXit/fvS1kF+UoHTPf7F69Wo899xzWLlyJX75y1+yQzL1KS0mOHFxcXjttdcAACaTCampqZDL5UhNTXV7vFwuh1arRVJSUvdESkREnTY6MhjrU6fg0cVj8MaXJfjHrmOoqql37pf5BSF0RgpCEm/Ehe8+x4WCd5GVlYXYseOx4p7lIkZO5Jl2DRNPTk5GUVERVq9e7Ux6iIio7xoc6oc1147Dgwtj8M+vjuP1HSU4XdlkOQiZN4ImJSFw4kJcKi7EKyURqMw/iGkhVXj1lRewevVqTJ48WcQ7IGpdu/vgqFQqt8PEiYio7wr280bGPBU+X7UQL6ROwZjIqzokS6QIiJmGqlo7/m/rIdzy+J/wzjvvYMqUKbj++uuxY8cOkSInap1HnYzbM4Pxxo0bOxwMERGJw8dLiqXxw/DJL+fir/ckYpbafX/LgMSlGLz8ZQSNnYMPP/oIc+fOxdy5c7F79+4ejpiodR6vRdUWdjImIuq7JBIJFowZhH9kzMB/HpyNJRMimx3jOzgGyhtXI+q+PyN4shZf7twFy5kKAI3Tidhstp4Om6gZtwlOWloaRo8e3fxgqRQymazVL7PZ3O1BExFR95s8TA7jsgRs+dU83Bw3FDKp61Qh3sphUFz7C0StfANrvqyD7u+F+PPf8zBp0iT885//RENDg0iRE7XQybioqAglJSWorKx0WXpcpVIhOTkZarXa7cnsdjsMBkP3REpERKIYHRmMl9Kn4lfJo/Ha58XIKzyB2oYrtTSywDAAwOb9Z5C393NUFpfgtttuw9NPP40nn3wSt956K2Qy96ufE3UXtwlOcXExKioqXJIboHEoeFujqAoLC7suOiIi6jVGKAPw/M2T8IukWGz43IK3dx/DpTrXWprgqdcgYPRMVBb+G8V7PsSdd96J3/3ud3jjjTcwa9YskSKngajFPjihoaHNtm3durXNE7IGh4iof4sM8cOT149vXPNqUQxC/Fz/rywLCEXYvLsQmbkRobNuxdGTZ1AhCQIA1NXVoa6uToywaYDxqJOxI+mprKxsts+xknh0dHQXhEVERL2dItAHDy8egy9WL4L+mrFQXLW4p8w/GPK5dyJS91es/Pdx3PfXAvw26xWMGTMGGzduZKJD3crjUVT3338/wsLCkJiY6LI9Ojoaa9as6bLAiIiobwj288bKBWpsX7UQa65tnuhIvLwBAFu/P4vXPi/Bj2fLkJGRgXHjxuGtt95iZ2TqFh4lOOvXr4dOp8Ojjz6KzMxMl33R0dFYu3Yt1q9f36UBEhFR3xDo6wXd/MZEZ7WbRAcAguN+gsiMjQideyeOnzyDZcuWYfLkydi/f78IEVN/1q6lGhxKS0sRFxeHuLi4Fo9x13eHiIgGjkBfL9w/X41lM0bi77uOIvtzC6wXap37pb4BkM+6FcFx16Gy4H0cP/Q56v0bR2LV1tbC29sbEomkpdMTtYtHNTjt+YGzWCwdDoaIiPoPR6KzfdVC6K8Zi7AAb5f9Mv9ghM27C2F3v4ZbNu7Fmve+xq9WrcHcuXPx2WefiRM09RseJTjl5eVuOxg3xQSHiIiaCvT1wsoFauzQL3Kb6EhkXmiw2fHOV8fxj50W7Ny1GwsXLkRycjJ27dolUtTU13mU4KSmpiIpKQn79u1rtq+yshLp6enNOh8TEREBVxKd7fpFeEQ7GsG+zXtJhC5cgcErXkPo5CRs+/RTzJw5EzfddBPKyspEiJj6Mo/64CQlJWHRokWIi4uDWq2GSqWCSqVCYWEhzGYzNBoNNm3a1F2xEhFRPxDk64WHkmJx54yReO3zYvz1iyOoqb8yM7K3fDDk1/4KgYkpuLDzbeze8y18A4JEjJj6Io+HiRsMBmzZsgUhISHIz8+H0WhEUVERMjIyUFBQ0B0xEhFRPxQW6IM1147D56sW4s4ZI+B19VpX4cMhv2E1vG9Zi8Uv78B75hO477778NRTT7XZXYKoQ6uJJycno6ioCDabDcXFxbDZbG0u4UBEROROZIgfnr1pErY9sgC3xA3F1eNZpL4BOFlRjV++vRs5H32KZ599FjExMfjDH/6A2tpa9yelAa9DCU5TV89cvHHjxs6ekoiIBqARygC8mD4Vn/xiHhaPj2y2X+rtB8Wyl6G45uc4Xwf8/Oc/x/jx47Fp0ybY7XYRIqberNMJztWMRmNXn5KIiAaQMYODkX1XAv71wGzMUClc9kmkMgRPWYzwe/+MsPl34fjJM/jTn/l3h5pzm+CkpaVh9OjRzQ+WSiGTyVr9MpvN3R40ERH1f1OHy/FOxgz85e4ExA5y7WQs9fZDyIw0DFphxImJd+LF/IOoqq7Db3/7Wxw9elSkiKk3cTuKqqioCCUlJaisrERISIhzu0qlQnJyMtRqtduT2e12riZORERdRiKRYNHYSMyLjUBu0Qm8mH8Q56pqnPtlAaGwIRR/2HYY2Tkf4uDrT2PdunV49NFHsXr1agQFcfTVQOU2wSkuLkZFRYVLcgMAcrm8zc7EhYWFXRcdERERAC+ZFLdNG4GfTonChu0WZH9uwcVa10U6a8LHYFD6s7j4+et47rnn8Je//AVr167FsmXLIJV2eY8M6uVafOPu1pTaunVrmydkDQ4REXWXQF8v/DJ5ND57bAFunz4CV40sh/+oqVDc+TIUi3+G8vOXcPfdd+Ppp58WJ1gSlUcpbWhoaJtzD1w9qoqIiKirDQr2w/M3T8KWX81D8rhBLvskUhmC436CiHtfQ9iMpcCYJFysrUddXR2OHz8uUsTU0zxKcNLS0hAWFoYjR450UzhERETtFzMoGBuXJ+Kt+6ZjdORVHZH9ghAy/x68+XUVFq3/HzJXP4uxY8di7dq1qKmpaeGM1F943Cj52GOPYdSoUd0QChERUcfMiQ3HRz+fi9/dOAGh/t7N9p+urMZHP/pAEqTA448/jsmTJ2PLli0iREo9xaMER6VSYd26da0ek56e3qmAiIiIOsJLJsVdM0fhs0cXYPnMkZBd1UHHf9RUKJf9AfL5d8Fy5BiWLFmClJQUnD17VqSIqTt5lOCsWbMGa9aswaefftriMRaLpdNBERERdVRYoA+evnEiPv7FXMyJCXfZJ/HyRuiMNETe9yrk4+fg0893wM/PT6RIqTt5tJp4WloarFYrsrKyADQOG1corswyyeSGiIh6i9GRwfj7fdNgOnAWz374HY6WXXTu8woZhNAbVqPhYgVWvLMfz9w4ER+9sxGTJk1CcnKyiFFTV/EowSkoKEBCQgIyMjLc7l+0aBHy8vK6JDAiIqLOkkgk0I6PxLzR4di4vQR/2HYI1XU2535ZQCi+KrHimqxPcDr7N6i+UIU777wTL774ost/4Knv8SjBUalUbXbKKikp6VRAREREXc3XS4YHFsbgp1Oi8PR/v4PpwBmX/XafACiXvYILnxrx1ltv4aOPPsLvf/97zJw5U6SIqbM86oOzYcOGNo/hRH9ERNRbDVcEYOPyBGy8KwHDwvxd9nmFDkLIjU8i/MbVuFQP3Hffffjtb38rTqDUaR7V4MTFxXXJMURERGJKHh+J2THhePWzwzD+z4LahsZmK4lEgsCxc+A3airOb38TismL0GCzQyYDbDYbl3zoQzx6U+5mMS4pKcHq1auxfv16bNu2rcsCIyIi6k7+PjI8sngMPvnlXMyNdR1tJfMLQqj2AXxcNQzp2bvw9/c+wpQpU7Bz506RoiVPeZTguOtcHB0d7Vy51W63Y/369R4HYbFYoNPpkJWVhaysLGRnZ3dbOa1W63F8RETUf6kigvC3e6fhj7fHITzIt9n+Pccr8MiGj3Hghx8we/ZsPPzww7h48aKbM1Fv4lETld1ub3V/UlISTCaTRwFYLBbEx8ejpKQEcrkcAKDX65GVlYVVq1Z1abmsrCyP4yMiov5PIpHg+slRmBsTgec/OoBNha5rVgVMuQbewybgfP4f8dJLL+GDDz7AG2+8gdmzZ4sUMbXFoxociUTS6v69e/d6nEAYDAZkZmY6kxSgcUJBvV7fpeUsFgsKCgo8io2IiAaW0ABvGFIm452MGRilDHDZ560cDnna8whbeC9Kjh7DokWLcOrUKZEipba0mODs2bMH999/PxISEqBUKqFUKpGXl+f899VfMpkM8fHx0Ol0HgWQk5MDtVrtss2RtLSWLHlaLi8vj8tIEBFRu8xUK/HhQ7ORPkkOryZLPkikMoRMuwWD7noFo278BY5XNzZpVVVViRUqtaDFJqq4uDi89tprABoThtTUVMjlcqSmpro9Xi6XQ6vVIikpqd0XFwQBgiBApVK5PZ/ZbHY7o6Sn5fLy8pCSkgKz2dzu2IiIaGDz85bhnnglls0fjyf+tR/7TlQ493krh6FGOQzp2Ttx2xQF/vZYGlJTUvDcc88hICCglbNST2lXH5zk5GQUFRVh9erVzqSnK7S2tINCoUBZWVmnywmCAKvVCpVK5XGCExkZ2axZjh3LekZDQwNsNhsaGhrEDmVA4vMXH9+B+BzvYPSgQOTqZuDNnUfxQv5Bl5mQ7Xbgb59/jypZCF5++WV8/PHHePPNN5GQkCBi5P2Du8+Au+Sxpf7B7e5krFKperyJRxCETpfLzs5utbNya9w9tEOHDnXoXOQZm80Gq9WKw4cPc94JEfD5i4/vQHxXv4O5gwDVDUPxwo5z+O5stfM4r5BBkKethazwXzi8/S3MmTMHDz74IFasWAGZTCbiHfRt7j4D7v4udzrBAYClS5e2eUxlZSVCQkLadb6mHYSvZrVaO13OZDJ1atG0U6dOISwsrMPlqeMaGhpw+PBhxMTE8BeECPj8xcd3ID537yAWwNw4O9748gheyD+E2vrLEwRKZQiZthR+o+JQ+fGLePnllzFp0iQsX75cxDvo29w9/0uXLjU7rrKy0u26YR4lOO2RmpqKzZs3t+tYR0DuamoEQWgxkWlvObPZ3OHaGwCQyWT8xSIiqVTKdyAiPn/x8R2Iz907kMkA3fwYJI0bjEdy92HfccG5z2eQCso7XoT/vk9QGz0HEokUtbU18PX1bXMkMjXXns9AS/s6lOBs27bNbXJhtVpRWFjY7vPI5XLI5fIWa2tampSvPeWys7NRXFzsMmzc0QdHr9dDqVR2KvkhIqKBLWZQEN69fyayt1vwcv6hK8s9ePkgMP6nWLf5ED4/XAbhP+sQ6OcDo9HIFcp7kEcJzp49exAfHw+JRNJim5enGWpaWhqKi4tdtjk6EbfWvNSRctnZ2TCZTFwQlIiIuoSXTIqfLYhB0thIPJK7F9/+6Lqk0Rc/nEZFsRUVB3Zg586deOeddzB37lyRoh1YPOq5ptfrYTAYsGXLFhQXFzf7KiwsdDt0u61z5uXluWwzGo0wGo3O7wVBgFardRkF1Z5yV+top2UiIqLWjBkcjPdWzsbKBWo0/X++xMsHoTfoobz2FzhbWoaFCxfi+eefh81ma/lk1CU8qsHRaDR47LHHWtwfHR3t8WriKpUKubm50Ov1SExMhMVigVKpRGZmpvMYR9NX0yap9pRzsFgsMBqNzoQoNTUVWq3W7bFEREQd4eMlhf6asZgXG4GHc/biVEXjSCuJRIKgyVr4Ro2F8GEWnnjiCQiCgKysLJEj7t88SnBiYmLaPCYnJ8fjIDQaDTQaTYv7VSoVysvLPS7XtLzBYGDTFBERdbuZaiU++cU8PP6vb/Dh11eWcvAOHw7l7evhvf1viJxxI+x2OzsedyOPmqjsdjsqKytbPWbNmjWdCoiIiKivCw3wxh9vi8MLqVMQ6HNllI/U2xfyRRn4w1cCfva2GX8ybsCvf/1rTujYDTxKcDIyMpCdnY333nuvxUSHq3UTERE1Nk0tjR+Gj38xD3Ej5M32f/TNSTy+7v/wzDPPQKvV4ty5cz0fZD/mURNVbGwsrFYrO+sSERG10whlAHJ0M7F+yw8w/u/KUkMSiRTypb8DPt2ITz/9CBqNBu+++y6mTZsmYrT9h0cJTllZGZKTk6FSqaBUKpvtLy0txcaNG7ssOCIiov7AWybFmmvHYXq0Ag/n7INwsQ5A4yirMO3P4D14NE7lv4q5c+fCaDTi7rvvFjfgfsCjBEelUrXZibikpKRTAREREfVXi8ZG4qOfz8VD7+xB0dErg2eCJiXDO2IUhP+sQ1VDly8yMCB51AcnNze3zWM4UomIiKhlUXJ//DNzBnTzXeeN8x0cg4h7/gTjsQhs3n8aFosFx44dEynKvs+jBCc6OhoAcOTIEbz33ns4cuSIc9+2bdtcjiEiIiL3HE1Wf7k7AfIAb+d2qbcvztfUI/PNrzBHex0SEhKwY8cOESPtuzxKcADg/vvvh0qlQmpqqstMwna7nUPEiYiIPLBobCT+++AcjB8S4rJdIpWhfspSlFeeR1JSEt58802RIuy7PEpwfv/738NisaCoqAgNDQ0u61ElJSVh9erVWL9+fZcHSURE1F8NVwTg3ZWzcHPcUJftAaNnIuJ2AyT+ctx9993Q6/WcL8cDHiU4BQUF2LJli3M5hqtnYAwNDUVoaGjXRUdERDQA+PvI8GLaFPz2hvHwkl752+ozSIWIO1+A/9CxyMrKwiuvvCJilH2LRwlOYmJim8dUVFR0OBgiIqKBSiKR4O7Z0Xh7xXSEB/k4t8uCwhBx6/OQz7kD/lOudWk9oZZ5lOA0XewSgNuHXFxc3LmIiIiIBrDpKiU+eGgupg6XO7dJvHwQOvs2GPIteOJf3+L5tetQVFQkXpB9gMdrUb3++uvO75s2UVVWViI9PR3x8fFdFx0REdEANDjUD//MnIGfTolqtu9vWwrw69/8BvPnz8cnn3wiQnR9g0cJzrp167B27VrExsZi5cqVyM/Px/r165Geno7o6GgIgoAVK1Z0V6xEREQDhp+3DK/cOhUPa0e7bPcKjUR42rOosctw/fXX44033hApwt7N42HiRUVFmDp1KoxGI/Lz87Fq1Srk5uZi9erV2Lx5c3fESERENCBJJBL8PCkWf7w9Dr5eV/5k+w0bj0G3GyALDse9996LZ555hn1zruJxghMaGorc3FyUl5cjPz8fxcXFsNlseOyxx7ojPiIiogHv+slR2KSbiYhgX+c2b+VwRNzxe/gOVuO559fCYrG0coaBx+MEB2jsbxMaGoqkpCTnzMWOmYyJiIio600dLsd/HpyNcU0mBfQKUmDQrWuhTP0dzOU+rZQeeDo0k3FYWFizIePR0dGcyZiIiKgbDQn1R45uBubEhDu3SX0D4B01Do/lfY3fvrUVd911Fy5cuCBilL2DRwnO+vXrodPp8OijjyIzM9NlX3R0NNauXcuZjImIiLpRsJ83/nJ3Im6a2nyE1ct/fh1///vfsWTJEgiC0PPB9SIercleWlqKuLg450zG7nAmYyIiou7l4yXFi2lTERniB+PnV/rehM65Hfb6GnzxxXtYsHAh8rdsQUREhIiRisejGpyrl2Zwh52ciIiIup9UKsGan4zDr68fD8efZ4lEAvmCexA6907s27sXc+fNw4kTJ8QNVCQeJTjl5eWorKxs9RgmOERERD3n3jnR+ONtGnjLGrMciUQC+axbEZaUiR++/x7PrzWIHKE4PEpwUlNTkZSUhH379jXb55jJuD3rVREREVHXuW7yEGy4K8FlrpyQhJ9iUOrTODsuFVXVdSJGJw6PEpykpCQsWrQIcXFxiI2NxZIlS7By5UokJiYiLCwMFosFjz76aHfFSkRERC1YMGYQ/nbvNAT5Xule66+KR+GJKty+YReeW/d7nDx5UsQIe5bHw8QNBgO2bNmCkJAQ5Ofnw2g0oqioCBkZGSgoKOiOGImIiKgdpquUeHvFdIT6e7tsL/xqN55cswrzFywYMElOhyb6S05ORlFREWw2m3Mm49dee62rYyMiIiIPTRkuxybdDIQHXZn12HfoOCi09+PwoUOYv2ABTp8+LWKEPaPDMxk7cCZjIiKi3mXs4BDk6GZgcIifc1uw5nqEJetw+NAhaBcvRnl5uYgRdj/OZExERNQPqSKC8E7mDAxqsn5VSPwNkM9fjm+/+Qa3LE3p1wt0ejTRn2Mm49DQUMTExLjsazqTMTsaExERiS86PBDvZM7Ardm7cK6qBgAQMj0F9vpaIG4uquts8PeRiRxl9+BMxkRERP2YOiII72RMx63Zu1B6vrZxnpw5d6AEwAP/MOM3SUMwLGoIZLL+lehwJmMiIqJ+LmZQMP6RMQPKQNcVxzfv+hoTJsfhZz/7Wb9rruJMxkRERAPA6MhgvLViOkL8rjTeeAUp0RA6DNnZ2Vi3bp2I0XU9zmRMREQ0QIwbEoI37kmEn3fjn3+Jlzcibl4D74hRePzxx/HPf/5T5Ai7DmcyJiIiGkDiRyrw2p3x8JI2djuR+gZiUMpvIAtSYNldy7Fjxw6RI+wanMmYiIhogFkwZhBeTJ/qXIXcKyQCg1J+iwaJDM+//Kq4wXURj0ZROThmMgaAkpIS52R/RERE1Df8dEoUKi7V4al/fQsA8IlUYfCd63E0ahQOnKrEuCEhIkfYOe1KcCorKxES0vxGV69eDYlEArVaDZVKhYSEBLfHERERUe+zbMZInKmoxh8/PQwA8IkYiQt1dtzzl91YNugEdPcsg5dXh+pCROc26pUrV6KwsBCCICA5ORlqtRqZmZnNkhdHj+s9e/ZAr9dj69atkEgkqK+v7/7IiYiIqNMeWTwax6wX8Z99VxbhPLT933hw859w6MC3ePnFF0SMruPc9sExGo1QKBQoLCzEn//8Zzz66KOt1szExcVhy5YtWLFiRb8bR09ERNSfSSQSZKVMRsLIMOe2oInJ8I0ai1deehFvvfW2iNF1XIudjHNzcz2eldhoNDLBISIi6mP8vGXIvisBo5QBABqHj4ff/DhkQQrce999MJvNIkfoObcJjkqlalZjc+TIkWZf7ib902g03RMpERERdRtFoA/euGca5AHeAACvIAUibn4CdQ0N+MkNN8JqtYocoWdaTHCaKikpQX5+PjIzM6FWq6HVapGXl4fCwsI2yxIREVHfEB0eiD/epsHlKXLgGzUGiuT7UVZxAdvN+8UNzkNuOxnL5XKX76Ojo5GRkYGMjAzExMQgJyenxQU327NeFREREfVOc2LDsebacXjuowMAgKApSxAwdg7+b18DFs2pQ7Cft8gRto/bGpzWkhSVStXqauJERETUt62YG40bp0YBaMwJZH5BKD53Aff8flOfmdTXbYIjCEKLBa6u3SEiIqL+RSKRYN0tkzG+yWR/tpoLeO9ZHX5y4y0oLy8XMbr2cZvgWK1WVFVVuS3QWu1ORUUFVxMnIiLqB/x9ZMi+K97Z6VjqG4iwhfeg9NQJ3H73CpGja5vbPjhFRUWt1tTIZLLuioeIiIh6iWFhAVifMgUr/tY4qChoyjW4VGLGJ/95D3/9+9u4e9kdIkfYshbnXw4NDYVCofDoZGVlZW6HjhMREVHflDw+EivmRGPjjhJIJBIolzyIkycO4P6VK6FdtABDhw4VO0S33CY4Go3G7RDw9khISOhUQERERNS7rLpmLAqOWLHvRAVkAaFQXvMQzr33DJ5/7W386ZlVYofnlts+OJ7W3DTFeXCIiIj6Fx8vKf5wmwbBvo31IgGx0xF135+xw2sqzlRWixyde24TnC1btnT4hDk5OR0uS0RERL3TCGUAnr9lkvN77/DhqLhUh4f/8RXOnj0rYmTutbgWFREREVFTN0yJwg1Topzf22ouIOfxW7Hk5lt73VqUTHCIiIio3Z65cQIGBfsCaBw67hs1Fnu//BSvvv53kSNzxQSHiIiI2k0e4ANDymTn92GLVkDqH4JHH/5lr1qQkwkOEREReWThmEG4bdoIAIAsIBRhi+5DdVU50lc8KHJkVzDBISIiIo89cd04DAn1AwAETlgEv5FTYHr/HWzb/qXIkTVigkNEREQeC/L1wtM/nQCgcRmnsGQdQmfdhv8e6x2pRe+IgoiIiPqcxRMGY/H4SACAT/gIyOfegfe/KcNuS5nIkTHBISIiok54+sYJCPK9sjBCw6VK3JbxECoq3S/a3VOY4BAREVGHDQn1x6OLRzu/v/jDlyjZ9g6W/2KNiFExwSEiIqJOunPGSIyJDAYABE3Wwjt8JP7z9kYcLDkmWkwtribekywWCwwGA9RqNQBALpcjMzOzS8qZTCbk5+dDEARYLBakpqa269xERETUPl4yKX59w3jcsXE3JFIZ5PPuwrn3nsFdD63Grg/+IU5Moly1CYvFgvj4eJSUlEAulwMA9Ho9srKysGpVyyuUtqecyWSC2WyGwWAAAAiCgPj4eBQVFcFoNHbrfREREQ0ks2PCoR0fifzvzsA/Zhp8o8Zi98c52F70OObGT+zxeERvojIYDMjMzHQmKQCwZs0a6PX6TpczGo0uSZJcLoder0d2djYsFkuX3QMREREBj/9kHLxlEkgkEsjnLwdsDfjV0+tFiUX0BCcnJ8fZxOTgSFpMJlOnyuXl5TVLlBISEto8NxEREXkuOjwQd88aBQDwGzEJg9J+h9Jxt+CH0z0/okrUBEcQBAiCAJVK1WyfXC6H2WzuVLmUlJRmSRARERF1n58tiHEOG/eP1gASGV7M/6HH4xC1D05rzUQKhQJlZe4nCmpvudzc3Gb7CwsLAQDJycltxhcZGQmJROKy7eLFi22Wo85raGiAzWZDQ0OD2KEMSHz+4uM7EB/fQceE+Mlw3+xReGXbYQBAXelxvPW8AdeNfAnXzZ7a7vO4e/4BAQHNjrPb7W7Li97JuDWCIHR5OYPBAIPB4Lb252ruHtqhQ4c6FBN5xmazwWq14vDhw5BKRW9JHXD4/MXHdyA+voOOmzfYhjd8paissaHhYjkufr8dv3z8txi98fftPoe75+/u73KvTHCadhC+WmtLrne0XGpqKpKTk1sdndXUqVOnEBYW1q5jqWs1NDTg8OHDiImJgUwmEzucAYfPX3x8B+LjO+icn1l9sO6TH+A7fBJ8o8aiZNcnOF61Dos0Y9tV3t3zv3TpUrPjKisroVAomm0XNcFxBOSuxkUQhBYTmY6Uy87OhkKh8Gh4uEwm4w+1iKRSKd+BiPj8xcd3ID6+g45bPisaG7aXoOxCLUJmpuHcu7/Dqt+uw54P/97uc7Tn+be0T9Q6N7lcDrlc3mKti1ar7ZJyeXl5EATBJbnpaPMXERERtc3fR4Z750Q3/ludCG/lCHy99X3sP3KmR64veqNiWloaiouLXbY5OhG31hG4veXMZjOsVqtLs5QgCBwmTkRE1M3unDESQb5ekEgkCI6/HraaC3jyj2/1yLVFT3D0ej3y8vJcthmNxma1LVqt1mXYeHvKWSwWrF27FgqFAnl5ec4vvV7frk7GRERE1HGh/t64Y8YIAEDghIUYvPxl7Pcdj3NVNd1+bdFHUalUKuTm5kKv1yMxMREWiwVKpdJlvSir1YrCwkKXJqn2lIuPj4cgCM0SIQBcqoGIiKgH3Dc7Gm/sOIJaH3/4Do5BbYMNmwqO4cFFsd16XdETHADQaDTQaDQt7lepVCgvL/e4nLsyRERE1HMGhfjhp1OjkFd0AgBw4YcvsK5gE+6f/1d4ybqvIUn0JioiIiLq3+6aOdL57/N7P8GPn/0DeTu+69ZrMsEhIiKibjV5mBxThssBAEGTkgFbA7L+tKFbr8kEh4iIiLrdXTMaa3ECRs+E1DcQ3376bxSf7b5FOJngEBERUbe7bvIQKAJ9IPHyQcD4BagrPYr/27S5267HBIeIiIi6nZ+3DDdNHQrgcjMVgJx/vAWbzf1aUp3FBIeIiIh6xC2axgTHZ3AMFNr7IZ18PXZZyrrlWkxwiIiIqEdMiArB2MHBjTMba66Hd9gQ5JlPdMu1mOAQERFRj5BIJFiqGeb8vuFSJd775FNU1zV0+bWY4BAREVGPuXFqFKSSxn+f+ccaHM/5HbYdON3l12GCQ0RERD1mUIgfEkcpAAABsTNguyDgr+9+3OXXYYJDREREPeraiYMBAP6jZwIAPjVtRm29rUuvwQSHiIiIetSSywmOT6QassAwVBz8Cju7eDQVExwiIiLqUUNC/TF1uBwSiQR+qnjUlR5DzqfmLr0GExwiIiLqcddcrsUJUE+Db9RY7Pi2BHZ71036xwSHiIiIetyCMREAgIAxszB42XpUBA6HpfRCl52fCQ4RERH1uDGRwYgM8XXZ9vnBc112fiY4RERE1OMkEgnmxjbW4lQf+xpn/vkE3s/f0WXnZ4JDREREopg3ujHBsTc0oProPuz+Yjtq6rtmVmMmOERERCSKuTHhAADfoWMBiRRVR77Gtz9Wdsm5meAQERGRKMICfTAmMhhSH3/4DI5FzYn9KCgp7ZJzM8EhIiIi0cSPCgMA+A2fAFv1eeR/UdQl52WCQ0RERKJJvJzg+ESNAQAUFhV1yXw4THCIiIhINAkjGxfe9Bs2AcrrHkbDkIkoPtf5+XCY4BAREZFohoX5IzLEF7JAOYImLoJXcDi+PiF0+rxMcIiIiEg0EokEk4bKAQB2ux31FWfxzXFrp8/LBIeIiIhENT4qBABQ8eU/8eNr92Jn4d5On5MJDhEREYlqwuUExzt8BABg/7ffdLqjMRMcIiIiEpUjwfFRNiY4laeO4KRQ3alzMsEhIiIiUQ2V+yPU3xteYUMAqQx1Zcfx3anOzWjMBIeIiIhEJZFIMCYyGBKZF7zDolBXdgyW0s4NFWeCQ0RERKJTRQQCALyVw9FwvhzFpys6dT6vrgiKiIiIqDOiwxsTHOW1P4fExx9HhdpOnY81OERERCQ6VUQQAEDqFwSJVIaSTjZRsQaHiIiIROeowbHVXMSlEjNq5INxviaqw+djDQ4RERGJboQiAFIJYKs+j9J/r8OF/Z/iVFV9h8/HBIeIiIhE5+MlRWSIH2TBSkAiRX3FGZy7wASHiIiI+rghoX6QSGWQhUSgvuIszl1kgkNERER93BC5PwDAKyQCDRVnUMoaHCIiIurrokL9AACyICVsNRdwWuj4SCqOoiIiIqJeYUhoYw2OT8RINFSOw2nr+Q6fiwkOERER9QpDLtfghM5MQ+jMNFRJO56msImKiIiIegVlkK/L95U1DR0+FxMcIiIi6hUUgT4AgDrhNKxbN8BavA819bYOnYtNVERERNQrOBIcW/V5VBX+G7KAUAgXaxHg6+3xuViDQ0RERL1CqL83JJLG9agAwHapCtYLHVt0kwkOERER9QoyqQRhAT6Q+l5el6r2IqwX6zp0LiY4RERE1GuEBXhD6hsAoHHhzXLW4BAREVFfFxbgA4lUBom3L2w1F1BxqWM1OOxkTERERL1GoG9jahI05Rp4hQ7C+ZqOLdfABIeIiIh6jaDLCY4iKQMAcKGDc+GwiYqIiIh6jUBfmcv3F2pZg0NERER9nKOJqnzb66iznsCF+OwOnYcJDhEREfUajiaq2nNHUHPyezZRERERUd/nqMGReHnDXl/X4U7GTHCIiIio13AmODJvwFaP89WcB4eIiIj6OH/vxk7GElnj+lMXq2s6dB4mOERERNRr+HhdTk0uJzg1HUxw2MmYiIiIeg0fWWOCEzhhAXwHx6Be0rFUhQkOERER9Rq+l2tw/EdOAUZO6XCCwyYqIiIi6jWcTVSX1dbbOnQeJjhERETUazgSHOGLd3D8/27H+bKTHTpPr2iislgsMBgMUKvVAAC5XI7MzMwuKdeRc4eHh3fkNqgLBQQEwG6349KlS2KHMiDx+YuP70B8fAficPTBsdfXwHapErW1fXQ1cYvFgvj4eJSUlEAulwMA9Ho9srKysGrVqk6V6+i5iYiISByOGpzgqT9BQMx02AIUsNvtkEgkHp1HYrfb7d0RYHvpdDrI5XIYDAbnNkEQEBYWhtZCa0+5jp7bx8cHdXV1sFqtCAsL68ztUQf5+vo6/+ckk8naLkBdis9ffHwH4uM7EIfl3HkseuF/LtsOPXctvGXue9VUVlYiNDQUFRUVCAkJcW4XvQ9OTk6Os/nIwVHbYjKZOlWuo+cmIiIicVzdyRjoWEdjURMcQRAgCAJUKlWzfXK5HGazucPlOnpuIiIiEo+vV/Pasr9+ecTj84jaB8disbS4T6FQoKysrMPlOnpuAM7mq0GDBjVr8ztz5kyL5ajryGQy2O12lJeXs2pYBHz+4uM7EB/fgTjsDTbYai66bDt1rgzl5UpERkY2P/7y3+yru56I3sm4NYIgdFu51o4pKSnB8OHDUV/ffAVThULRoZioYyIiIsQOYUDj8xcf34H4+A7E99zLwHNtHFNVVYXQ0FDn96ImOI7+MO5YrdZOlevouQEgKioKx48fR3BwsMe9tomIiKjn2O12VFVVISoqymW7qAmOozbEXW2KIAgtJintKdfRcwOAVCrFsGHDWo2diIiIeoemNTcOonYylsvlkMvlLdaoaLXaDpfr6LmJiIio7xN9mHhaWhqKi4tdtjk6CCcnJ3eqXEfPTURERH2b6BP9WSwWaLVal0REr9dDrVY7l1QQBAGpqakwGAzQaDTtLnf1MSaTCU888QR8fX3h6+uL1NTUZss2xMfHY82aNc4EKDs7GwBcZj7uriUiBgKTyYT8/HwIggCLxeL2HXTV8+U7aF12djaKi4tdJsJ04Oeg+7X2/PkZ6Hn8me9duuQ52nuBoqIi+6pVq+y5ubl2g8FgNxgMLvuLi4vtcrncnp+f71G5psc89dRT9muvvdZ5THl5uV2lUtkzMzNdjgfg8nX1fkcs5eXlzm2rVq1yuXZ7jhmI8vPzXZ6Bu3fQVc+X78C94uJie2Zmpj0zM9Mul8vtq1atcnscPwfdoz3Pn58BcfBnvvfoqufYKxKcnpCSktJsm9FotAOwFxcXO7dlZmbajUaj3Wg0umxvuv/qX0rl5eX2prlie44ZiNrzDrrq+fIdtE2j0bSY4PBz0P1aev78DIiDP/O9R1c9R9H74PSUvLw86PV6l20JCQkAXJdtcDRxZWZmup0FmUtEdFx73kFXPV++g87h50A8/AyIgz/zvUdXPccBk+CkpKQ0e2AtEQQBJpOp2XIOXCKic9p6B131fPkOugY/Bz2PnwFx8WdefF35HAdMgpObm9usg1JhYSEA1xFV+fn5MJlMzpoFrVbrfKDdvUREf9fWO+iq58t30Hn8HIiDnwHx8Ge+d+jK59irl2robgaDAQaDwSVTNBqNzu81Gg10Oh1SU1ObDTd3p7NLRAxETd9BW5l5Vz1fvoO28XPQO/Ez0H34M983ePIcB0wNztVSU1ORnJzsMgQQQLNqMY1GA4vFApPJ1O1LRAw0V7+Drnq+fAedx8+BOPgZEA9/5nuHrnyOfaoGx2w2IyMjo93Hb9iwwTlvTlPZ2dlQKBQwGo0u2/V6PdLT013KOJZ8sFgszqrL7loioi/oznfQVUtw8B24aukdtISfg9Z15/PnZ6DjOvNe+DPfe3Tlc+xTCY5Go0FRUVGnzpGXlwdBEFz+sDoeWlZWFtRqtcsPuSNjVKlUkHOJiG5/B13xfPkOOoefg9Z15/PnZ6DjOvNe+DPfe3TlcxxQTVRmsxlWq9WlWcrRax5o7A9ydSdYR/WkJ8s/cImIlrX1Drrq+fIddBw/B+LiZ6Dn8We+d+my59iVk/P0ZsXFxfaUlBR7bm6uy1dmZqa9qKjIbrc3zrSbm5vrLOOYabfptuLiYrtKpXI596pVq+xGo9GjYwai9ryDrnq+fAdtczeTt93Oz0FPaen58zPQ8/gz37t01XMUfS2qnhIWFtZi7+umj8CxVhLQmDHqdLpmGaPZbMamTZuQmJjozCqv7qzcnmMGmva+g656vnwHzQmCgLVr10IQBGRnZ0MulyMtLQ1qtdrl2fBz0D3a+/z5Geh5/JnvXbriOQ6YBIeIiIgGjgHVB4eIiIgGBiY4RERE1O8wwSEiIqJ+hwkOERER9TtMcIiIiKjfYYJDRERE/Q4THCIiIup3mOAQERFRv8MEh4j6NbPZLHYIRCQCzmRMRP2aWq1utnAfEfV/rMEhon4rLy+vxdWHU1NTER8fD4lEAolE4lzvpi06nQ4SiQRhYWHQarXIzs7uypCJqIswwSGiPsNiscBkMrX7eKPRCJ1O53Zfbm4ucnNznQlQexIck8kEq9UKANi6dSvy8/ORmZnZ7niIqOcwwSGiPkOn08FgMLTrWEEQYLFYoNFoWjzGZDI5E6D2JDgWiwWCIEAul7d6XiISHxMcIuozTCYTtFptu47NyclBSkpKq8cUFxc7j2mrn052djbS0tJgMplabPYiot6DCQ4R9QmOpqn2Jhe5ubktNk9dTS6XtzraymKxQKVSOWt52ptkEZF4mOAQUZ+Qn58PAO1qGhIEAVarFSqVqtVj1Go1ALgkL+44Oit7mmQRkXi8xA6AiKglJpMJRqMRQGOSIZfLkZqaCqCxP05LiUZ2djbS09PbPLejvEqlQl5entvj8vLynB2J8/PzIZfLW02ciKh3YIJDRL1WcnKyMwmRSCRYs2YNVq1a1Wa5TZs2YevWra0eU1BQ4Ox/k5iYiLy8PGcHYgdBEADAuc1kMrXZr4eIegc2URFRr+dJ05CjqalpotIWR43M1c1U2dnZzoTG0UeH/W+I+gYmOETU63nS/6a1uW8cBEGAUql0fu8uwbl6tBT73xD1LWyiIqJez5Oh2Xl5eSgqKvLofO4SHIvF4nIM+98Q9S2swSGiXk0QBJjNZmfn4taYzWaoVKo2m6cKCgpcaoMcxzvmwsnOzm42QzHnvyHqW5jgEFGv5q5pqKXlGtrTPNUSx1BxR5LUFPvfEPU9THCIqFcrKChwaRoSBMHZJ+dq7RnlZLFYXPrfOKhUKhQWFrqtqWH/G6K+hwkOEfVqSqXSpUZl7dq1WLNmTbPjTCZTuzoh6/V6t01YKpUKgiC4XTxz06ZNzmOIqG+Q2O12u9hBEBG1RBAEpKamOpOXNWvWuE1QdDodUlNTW6xl0ev1yMvLg8VigVwuR0JCgktNUHZ2NhQKhbMGSBAEZGRkQBAElxocrVbbrrl4iEhcTHCIqF9Qq9VtLphJRAMHm6iIqM9zrBVFROTABIeI+rxNmzZ1ePQUEfVPbKIioj6PzVNEdDXW4BBRn9Z0vSgiIgcmOETUpxUVFbF5ioiaYRMVERER9TuswSEiIqJ+hwkOERER9TtMcIiIiKjfYYJDRERE/c7/A5fDU5Wx8WiRAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.figure(figsize=(6,4))\n", "plt.plot(obj.time_xi, obj.ecc_xi, c='C0', lw=3, label='from gwModels')\n", "plt.plot(obj.time_xi, et_model, c='k', ls='--', label='$\\\\rm gwEccEvNSv2$')\n", "plt.xlabel('$t/M$')\n", "plt.ylabel('$\\\\rm Eccentricity$')\n", "plt.ylim(0.00,0.08)\n", "plt.xlim(xmin=-2500,xmax=50)\n", "plt.legend(frameon=False)\n", "plt.tight_layout()\n", "plt.show()" ] }, { "cell_type": "code", "execution_count": null, "id": "cf923fe2-7343-4d39-acff-7ca6da819ed2", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "kitp-py310", "language": "python", "name": "kitp-py310" }, "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.10.18" } }, "nbformat": 4, "nbformat_minor": 5 }