{ "cells": [ { "cell_type": "markdown", "id": "431b1fc9-0685-4c50-999a-2e8026be504a", "metadata": {}, "source": [ "## gwEccEvNS : a simple model for eccentricity evolution in time domain for non-spinning binaries\n", "Presented by Islam and Venumadhav (https://arxiv.org/abs/2502.02739)" ] }, { "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, # choice for the distance between peak\n", " fit_funcs_orders=['3PN_m1over8', '3PN_m1over8'], # default choices for the fit orders using PN+pseudo PN expressions\n", " ecc_prefactor=2/3) # default choice for the Newtonian limit\n", "\n", "# calculate eccentricity from gwEccEvNS model using q\n", "et_model = gwModels.dynamics.gwEccEvNS_model(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": "iVBORw0KGgoAAAANSUhEUgAAAjgAAAFwCAYAAAC1lgFqAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8ekN5oAAAACXBIWXMAAA9hAAAPYQGoP6dpAABmPklEQVR4nO3deVhTV/oH8G8IuwIhqCiuBHDfCOBuXSAutW4ti512ugqp02U6rTUynXY6004Vfl2nyzRg97EdgXasdtoq0a7aqhC11h2C+waECyh7kt8fNCmRgOyX5ft5Hp9H7j333nPvJfByznvOkVgsFguIiIiIuhEnsStARERE1NYY4BAREVG3wwCHiIiIuh0GOERERNTtOItdAQAwGAxISkpCUFAQAEAmkyEhIaFNjtPr9dDpdACAwsJC+Pn5Yc2aNW18B0RERNSZSMQeRWUwGBAWFoa8vDzIZDIAgEajuWEg0pTjDAYDMjIy7M6j1+uxbt06pKent9s9ERERkbhED3DUajVkMhmSkpJs2wRBgK+vLxqrWlOOU6vV0Gg0UCgUdseqVCpkZma28Z0QERFRZyF6Dk5aWpqti8nK2iJj7Vpq6XFGo9EuALIyGo2tqDERERF1dqIGOIIgQBCEei0sQG2wotfrW3WcWq1GSkoKYmJiIAgCACA5ORlqtbrtboKIiIg6HVEDHIPB0OA+uVyOwsLCVh0XFRWFpKQkZGRkwNfXFzExMYiKimpSAjMRERF1XZ1iFFVDrK0urTkuOjoa+/btsyUcA0BqaqqtO8sRs9kMNzc3AIBEIrHbV1BQ0KI6ERERUev06dOn3jaLxYK8vDwEBATAyem3dhtRA5zGgozG8mSaetz1I6aSk5Oh0Wig1+uRm5vb4DkuXLiAmpoah/t8fHwaPI6IiIg63uDBg3H27FkMGjTItk3UAEculwNw3FIjCEKDgUxTj4uPj0d2drZt35o1axAdHY2wsDCkpKQ02FXl5eUFoLYrrLFgitqPyWSCwWCAQqGAVCoVuzo9Dp+/+PgOxMd3IK6mPv+SkhIMGzbM9rvbSvQWHJlM1mBrjUqlavFxBoPBFgjVpVAokJiYaBf4XM/aLSWTyeDr63uj26B2YDKZ4O3tDV9fX/5gEQGfv/j4DsTHdyCupj5/677rU0pEz8GJjY2t111kTSKOiopq1XENJSPLZDKEhYW1uM5ERETUuYk+D45Go7El/1pptVpotVrb14IgQKVS2Q0bb8px0dHRSE5OtisjCAIyMzObNJLqlR0n8fXxKyitqG7WPREREZG4RG/BUSgUSE9Ph0ajQUREBAwGA/z8/OwCEKPRiKysLLsuqaYcl5SUhJSUFNusxwDg5+fX5GUaXl73LP41dAI8Bo/G+OAhmBwox2SFHyKG+ULm6do2D4CIiIjanOgBDgAolUoolcoG9ysUChQVFTX7OACtmvOmVP85SvWfAwAuyAfh22ETIVc9AIkEGOHvhSkKP0wKlGNSoBx9eru1+DpERETUtjpFgNNZ+d/5Aqrz81B57ggqzx1BVf4pAIDFAuzd8gG+v3AMboNGw23QGIwaMxZTgvvVtvIE+qG/j7u4lSciIurBGOA0YnLoeBwtHI2qiQsBAObqCtu+6sKzKDvxI8qO7wIAXHZxx48BI+EVsRSeQREY6ueJyYFyTAr0w+RAOQbLPUW5ByIiop6IAU4j3rs3HB69vbH/jIA9eYXYYzBCf6YIlTVm+C18BL5zV6LywjFUnjuCinNHUHn+KHqPrx3BdargGva++Thc+gyB26DRGDZaiRnjg215PMP8POsNaSMiIqK2wQDnBtxdpJga5IepQX4AgKoaM34+J2BPnhF78ozI9vLCtcDaPCCLqaa2/wqA+ZqAqsu5qDi1H6VZn6FgM3DAdwBSB46Bb+RKDOhbm78zWeGHKYFyBPfrzYCHiKgd6XQ62yATlUqF6OhokWvUtgwGA7RaLVJSUiCXyxudsT8oKAhGoxEJCQlITExs0aS2arUaaWlpSE9Pb3RaF6B2ZQGNRgODwdBovdoSA5xmcnV2QvgwOcKHyfHgHKDGZMbhCyW2Fp69p4woraiBtLcvBj28ETXG87WtO+eOoPL8YZSd2AW/hQ/jSmkl0j/7Au8c+ALuA0fDf0QoZk+LwPQQf0wN8oOiTy8GPEREbUilUqGoqAhZWVliV6VdKBQKJCUlwc/PD+vWrYNer3c4EEen00GpVEIulyMpKanF19NqtU1+lkqlEhqNBmq1usXXay4GOK3kLHXChMEyTBgsQ8JNQTCZLTh2qQR7DEbsySvE3jxXFPkNgteEeQAAc8VVSJxqZ12sLjiN8pN7UH7iRxR9DRx/xwPvDRyNXqNugmLazZiiqG05mqrww1B2aRERtZher4dCoYBMJrtha0NXJ5PJEBsbW29uOKuWLmTdWo5WF2hPDHDamNRJgjEBPhgT4IP7ZgTCbLYgJ/8q9hgKbd1a+aWVAADviGXoPX4eKi8cR+XZX1Bx7jAqzvwMF79BuFJaiS0HL+D915IgcXJCwEglIm+ajpmjB2Gqwo9Jy0REzdST1hZUq9WIjIysF+AIggCFQiFSrToWA5x25uQkwXB/Lwz398Lvpw6DxWLBqcKy3wIeQyEuuIXCIzAUAGCpqbKN1rJYLLh2SAfT1UIU7/oYR991Rmr/YLgNGYuQWbdidthoTP21lSdA5iHmbRJRD2I2W1BUViV2Nerx9XSFk1P9lm69Xg+tVguDwYDk5GRbS45Go4FCoYBarUZmZiYA2LpskpOTMWzYMFy8eBGbN2/G2rVrAdR272g0GgBAamoqDAYDjEYjsrOz7fJfNm3ahMTExBvO1SYIAjQaDcLCwpCdnY2goCAolUpotVrExcXZrmWtf1JSEmQyGdLT0yGXyxEZGQlBEJCenm53LWsXVEZGhl2uUVZWVqMtWNbnA9Tm9KxZs6befplM1mhrjPUc1oUyG8t1SklJgUKhgCAItgWuWzN/XV0McDqYRCJBYJ9eCOzTCysmDYHFYsFZYzl+NBTgx9xC/GgoxOUSV1vZgQ+8japLOag49wsqzx6uzee5cAwXxsxFRvY5bPruEIp/ysDgMeGImn0T5k4MwhSFH/y9OQ8PEbWPorIqhD2nE7sa9WT/JQp+DiZdVSqVUKvV0Ol0dr+wExMTsW7dOgBAXFwcNm3aBACIiYmBWq3GnDlzcPLkSUilUqhUKmRmZiIqKgpJSUlQq9UwGo22X95BQUHQaDR2OS3x8fGNLuxsLRMXF4fo6GgIgoCwsDDk5ubaghCj0WiXxCsIAgoLC21BiDWIctQqEx0djU2bNtkFGI11T1nvu+56jtb7BmqXSAoKCrIFIIIgICYmpt45rPcD1OY9KRQKh4GedbmlutfT6dru+4oBjsgkEgmG+HliiN8QxEXUBjx5Bdfwo6EQP+YW4ieDEQXSkXAbOBKYHA2L2YTq/NNw8RsMAKg4exilWZ/hSNZnOPI+8K8+Q+E+ZCyGjp2ExcuWY6rCD1MUfujrxZmWiYjqkslk0Ov1tl+wSqUSer3eNtrKZDIBqE3eNRqN0Ol0iIqKglwuh8FgsGsJuT7AUCqVDS74XFdGRoYtKJLJZDAYDLaWD6B2YenrE3PrHmM0GhvsclKr1QgKCoIgCLZzN9SiVPe+696T9b7Dw8ORnJwMy68jha31rXs+g8GAjIwMu3PExMQ0mAsEAOnp6YiNjYVMJoNCoUB4eHiDz6q5GOB0MhKJBIq+vaHo2xt3TB4Ki8WCnCtX8ZOhtnXnJ4MRRqffvpk9R0xHwP3/QsXZQ6g4exiVZw+hVP8/nLxswEcDwvDRnjMoy9mLvq7ViJwzGzdPHY+pCj/4eLqId5NERJ3E9cFBVlaWw4BBoVDYWnAcHSeTyRAUFNTs6yuVynqtKnXPbQ0irF1N1jyiukFQQ6wtJ2lpaUhISIBer2+wu+hG920Nkhqj0+kgk8nsWmFyc3MbDPSio6Oh1Wrh6+sLpVKJuLi4el1ircEAp5OTSCQI8fdCyK85PGazBSeulNZ2Z+XW5vEUSwbDpc9geIXeDIvFghrhEswVV23nKM3agvzTB3Dk43X4l48/3IeMQ8jEybh54QIsCB+J8GG+cHeRiniXRETiuP6XdkePMFKr1Vi3bh1SU1ORkpLisKXD2n2mVCoRHh4OtVqNjIwMREVF3XBEmFqthlarRUJCQqMBSlPu+0ajoKwJzHXrdKP6ZWZm2lqPrPfeVkEOA5wuxslJgpH9vTGyvzfunV47SuvIxRL8ZCjET78mLpdKBtgd43fzo6g89wsqzhxCxZlDuHpIh/2HdDhTeA3/PhgFV2cn9CvYj4WRs3Dz5NEYN9AHzlInke6QiDo7X09XZP+l8w219vV0bfU5oqKibHk5dRkMBsTFxbX6/I6kpqZCp9MhOjq6wVwajUYDlUplC1RUKlWTEnKtXVzWgKghN7rvpnS5KZVKh+doqPUnJSUFCQkJtoWzExISEBkZyQCHajk5STB2oA/GDvTBypkKmMwWHL5QbEtY3pdnxDXvPnAePRu9Rs8GANSUFKDi7CG4D50AALhWcAG7tE9ilxZ4Rj4QXoETMHHSdCxdEIWbJ4/iLMtEZMfJSeIwmbc7UCqViIqKgk6nw5w5cwDU5qcAaHQ0UEtbfrKzs5GQkNDoua1dTdYZgK1BUEPXrDtTsHXen+uTja9X976tgdD1952QkGALSqzX1+v1tnpERUUhPDy83sgtaxfZ9QRBsDtf3XtrCwxwuhmpkwTjB8kwfpAM6llBqDaZ8cv5YuzOLcTu3ALsO1UEePdB7zFzfjvGwwt+i/6EitOHUHH2EIqyv8DX2V/gmw1/x0t/3AR/WS9M7CPBpEA5bp40kkPSiahL0ev1WLduHQwGAzQaDeLi4mA0GpGUlGQbOl639SQ9PR0ajQY5OTm4cuUKSktLbaOh6p4rOTkZa9asQXJysm1GX4VCAblcjnXr1tmGgDe2FIJKpYJEIrENvZbJZFCr1fUCgri4OLsWGLVaXS9gsQ4jT0lJsV3bOgzeyhpU6HQ6yOVyJCcn21qFrPdtbanJzc21GwWm1WqRnJyMjIwM2/1YW22sgVRmZiY0Gg2MRqOtS8ua/1P3HViHuwO/jaYyGAxITU1t8nu9EYmlbko0AQBKSkrg4+MDo9EIX19fsavTpiqqTdCfLsIPOQXYlVuIQ+cEmK/7DqgpvoyKM4dgKhPgM7n2A1T07fso+SkdLn2GoN9wJabMmInoRfOwIHw4ZG3QLHw9k8mEkydPIiQkBFIp84M6Gp+/+PgOxNfe70AQBMTHxyM1NRUymcw2F4xGo0FMTEybzQfTVTX1+Vt/ZxcXF8Pb29u2nS04PYy7ixTTgvtgWnAfAEBxeTX2GAqxO7cQP+QUIOfKVTj7+KP3OH+749wGjkKv0bNRceZnnN+9GZ/s3oxPkgHZjDswI/YBTAv2Q/jAXpgxciA8XPnDmIjoRlJSUqBWq20tGdYRU0lJSQ0Oq6amY4DTw/l4uGDemP6YN6Y/AOBySQV25xZgV04hduUU4GJx7azKnsGT4Bk8qXaUVtEFW8KyS/9gHDpfjJ/PCTj3xu/h0luOwAmTERkZhbuWzUd4cIDDmUWJiHo6pVLpcCVunU4HlUolUq26DwY4ZMff2x3LQwdheegg26SDu3ILsTunALtzC1FcXg0X+UC4yAfCa+IC23GWqnK4Dx6HijM/49j2j3Fs+8d4I1GKXoNG4fbElzEvfARmhvRFfx/OsExEBPw2hPr65RFutLwBNQ0DHGpQ3UkHfz9lKExmC45cKMGu3ALsyinAvlNGVFSbAQBObp7ou1QDi8WM6vzTqDh9sPZf/mlsN5Qh89TPKM/NgumXLxA29SbctngBfrdwJnq7t33+DhFRV9GUuWyoZRjgUJNJnSQYN8gH4wb54IFZQaisMUF/Wvi1S6sAB88Vw2R2gmu/QLj2C4R3xDJYLBbbEPPqovMoOqHH9mN7sf3dF7DKwxsDR4VhwbJoPHLfHRg9wJvdWURE1CYY4FCLuTlLMTWodjXzx+eNQElFNX7MLcQPJwvw/cl8nCoss5s/xzt8KXqPn4/Kc4dtLTxn9d9gk7Mvtl8bCnkvV/Q+uR0TQgbjntsWoRfnGiQiohZigENtxtvdBfPH9Mf8XxOWzxSW4fucfPxwsraFp6SiBk6u7vBQhMFDEQYAMJWXAqYaAEBhaTkObHoD31eV4/Wn/wjPfkMwNmIGFi+cB/XvlqGvr49o90ZEZGUwGKDVapGSkgK5XG43z0xubi50Op1tDafOco3rz1d3MkCgdqXwlJQUhIeHIykpybaIZnJysq1MYWGhbVbktlwzqr1wHhwHuvM8OGKpMZnx8/liW+uO/owA0/UT8ODXWZZPH0TF6QOoOH0QpqtGAMAQtRZTleMQPsAVfasvY8WiSLi5dc+ZVDsDzsEiPr4D8d3oHYSFhSE8PLzekG5BEBATE9OqAKe9rpGcnIx169YhISHBtiK5lXUCPiu1Wg21Wm23YnhMTAwA2K0Y3l44Dw51Cc5SJyiH+EI5xBePRIag1NqdlVOA708WIK/gWm057z7oPS4SvcdF1g5JN55DxbkjgE8AfjIYodv6NQo/fxH3ubpj2NhJiIyKwsoVSxExcRyXkyCiDtXQ4pPWtaI64zWsMxarVCqo1Wq7pRGuXw09LS2tXmCVmpqK+Pj4Zl9XDMxyIFF4udfOv/P3pWPx9erZ+H7NHKy7dRxuHtcf3u61cbdEIoGL32B4TZhvC17cBgyHz8w74dIvCIYDu5Ca/DQmKydgYOhsvKo7iYNnBVRVVYt4Z0TUU1lnIgZq57hpj5XJ2+Ia1pFb1taYplzLSiaTISIiotnXFANbcKhTGCz3xO2ThuD2SUNgMltw8IwRm/ecwNEiC/afEVDza3eWi3wgZNNWANNWwFxZhoqzh1CRtx+VfYfiZd0JvKw7gcL0p+CJSsyYNQd3xSzBwsjZcHXlcHSinsy6NlNYWBiys7MRFBQEpVIJrVaLuLg4aDQaALXrLVnXdJLJZFi/fj369u2LefPmQRAEpKen23XZ1FU3GLh+6HdKSord13WXYWhsX1OukZGRYVvnyTpxoMFggEqlgkKhgFarrbeIZXp6Onx9festdlmXUqmESqWCVqu1u5+ukH8DMMChTkjqJMGEwTJ4VvgiJCQE5TUW/JhbiO9O5OObE1dw1lgOoHbuHc/gyfAMnmw71mKxwOIhQ0FuFj599w18+u4bkLq6I2TiZPztuedx69wpcJay4ZKorYwZM8bh9s8++wzBwcHIycnB0qVLHZY5fPgwAGDbtm147LHH6u0PCgrCli1bAABvvPEG3nzzzXpl5s2bh5dffvmG9YyPj0dcXByio6MhCALCwsKQm5tr+8VtNBrtZhUWBAH5+fkYPHgwZDIZEhMToVQq6wUKWVlZSE5ORmFhITIyMhzmplj3W/NbMjIybCtuN7avqdeIjo62LdJprb91kc2GghGZTIakpCRoNBrExsY6XAzU2pVl7QqLioqCRqPpMvP2MMChTq+3mzNUo/2hGu0Pi8WCU4Vl+Ob4FXx7Ih8/5haissZsKyuRSNDnlsdhMZtQddmAilP7UZ6nx7GsXXhs80k8t6sUE3yqcHbnh4hbtgjRS25G3759Rbw7IuoIGRkZtiBCJpPBYDDYZg0GgNjYWLuRSgDw6aef4v777wdQGwBdH9wAQHh4eKMtGtaWo6KiItu2TZs2ISIiotF9zbkGUBt8GI1G6PV6WwtTQyuYW61ZswZarRbx8fEOAzOFQmEbsZWZmWlbQiI9Pb1LzLTMAIe6FIlEgsA+vRDYJxD3Tg9ERbUJe/KM+PZ4Pr49cQW5+bXJyhInKdwGhMBtQAh8psbCXF0BJxd3lFTU4Mvsb1H4eTq+/TwdfwAQEDwacyOjcN+K2zBn9k3i3iBRF2NthWlIcHDwDcvMnz//hmUefPBBPPjgg82un5WjfJW6AYt1oUtr64k1ODh79ixCQkKadA21Wm13Db1eD6PRCJlMZhdsWIMJnU7X4L7mXMMa0CQkJECr1UKr1UKn0yE2NvaGdU5PT0dYWBj0en2DZerOtqzRaBAfH98lAhy21VOX5u4ixazhffH04tHY8XhtsvJzy8ZCNdofveqsau7k8tsaWL3HzMHAVe/Cb+Ej8Bw5E5fOn8O/tf/EkpWP4w8bs5GedRbbv9mF06dPi3FLRNQO1Go11q1bB0EQkJyc7HC17ri4OGzatAkGgwHh4eGIj4/Htm3boNfrm9Qto1Ao7PJzsrKyGk0CbkmCsKNrWKnVaqSlpQGozdW5UQsOUBv4JSQk1Es4FgQBGRkZ9conJSVBEIR2SaBua2zBoW5lsNwTd04ZijunDEVVjRlZp4349kQ+vj2ej2OXSm3lnL37ovf4eeg9fh4sFjOqLuUCFjO+OHQJXxy6hIvvP4qqSznwHxIE1bz5uDN6KWbNugnu7lwslKirSk1NhU6nQ3R0tMPupujoaGg0GqhUKiQkJMDLywuRkZEIDg7GAw880Kxr6fV6KBQKKBQKh8GAIAgNjoISBKFJwYn1GlYKhQJyuRwZGRkNDi93JCkpCYGBgUhKSrJLON63b5/DlhqFQtGk+omNLTjUbbk6O2FaUB8kLhyFrx69CT8lRiL5tvFYNG4AvNx/i+0lEqfa7qyAEbZt3pOj0WtcFAqMRfj3hjexYMF8ePn44v/e24zSimpwfkyiriU7OxsymazB4Ab4rXXEOsuvtVxDrRVGo7HB62k0GluAY00mthIEAWlpaY3ua8416lKr1YiPj2+wxen6GYyB2u651NTUetdKSUmBTqez22YNELsCtuBQj9Hfxx2xEYMRGzEYNSYzDpwV8O2JfHxzPB+Hzhfble01cgZ6jZwBi8WC6vxTKM/LRkWeHq8dqETKyUwESwuR/c5fsGDBQsQuuwVz586Fp6enSHdGRDeiUqkgkUggk8kgl8tto46uHyIdFxdnFxzExsbitttusytjXfZAr9fburysrEm5dROY09PTodFokJycDIVCAaPRaLtuQ/uaew2rhIQE5Obm1mthMRgM0Gg0yMjIgCAItmHwVtHR0di0aZPdMdak7LrXrru9s+NSDQ5wqQbxdfQ09VdKKvDN8XzsPHYF35/Mx7UqU6Plyw3ZKNz2Okwl+QAAqYsrQidNxb133o4/PKBu9NiugMsEiI/voO0IgoD4+HikpqZCJpPZJrDTaDSIiYlpcB4YvgNxcakGojbQz/u31p2qGjP2nTJi57Er+PrYFRh+XUaiLg9FGAY+8A5qCs+hPC8b5YZsZP20G8eLLNjvGYa5I/uh5PC38Pf1QlRUFHr37i3CXRERUNvVolarbS0W1hFTSUlJDpONqXtggEN0HVdnJ0wP7oPpwX3w1C2jcargWm2wc/wK9hiMqDLVzrsjkUjg0mcwXPoMhnfEMpirKmCuvIrMI5eReeQyzqesRU3RBUhdXDB52kzELl+CW265pd56L0TUvpRKpd0kflbWeV2oe2KAQ3QDw/r0wn0zAnHfjEBcq6zBDzkF+PrXgOdySaWtnJOrO5xcfxtl1S/mGZQbslGesxe7f/gOu7/diUcffRQvv/8pVv1uCZwlgNlshouLixi3RdRjWAMba54LAFv+SldJmKXmY4BD1Ay93Jwxf0x/zB/THxaLBYcvlODrY1ew8/gVHDgroG5Gm4tvAFzCAuAdthjmqnJUnD6IckM2Xv4ZSM3JRIjlHDJffgxRKhWWL12MhQsXwt/fX7ybI+rG6k5WRz0DAxyiFpJIJBg70AdjB/rg4cgQFF6txLcnahOVvzuRj5KKGltZJ1cPeIZMgWfIFADAtSoTfjx9BTXeAdiy+b/YsvlTSCQSjA8Nw8p778ZDDz0k1m0REXULDHCI2ohfbzfcqhyEW5WDUG0yI/t0EXRHLiPz6GWcLiyrV95j6AR43PUiTNeE2kTlnH34+Rc9nv/wK5QoIhE1yh+Hvt4Mua8voqKi4OXlJcJdERF1TQxwiNqBi9QJUxR+mKLww5OLRiE3/yoyj1yB7uhl6M8U2XVlSXvJ0HtsJHqPjYTFVANzVTlSvjNA+20OLry5GjVXi+Ds4oqbbroJSxbfgsWLFzc4URkREdVigEPUziQSCYL7eSG4nxdWzQ5CwdVK7Dx2Bbojl/H9yQKUV/82545E6gyph9evxzmh3x3/h/LcfSjP3Yevv/kWO3fo8Oijj2Lb1z9g3uzpMJlMkEgkcHLipORERHUxwCHqYH16uyE2fDBiwwejotqE3bkFyDxyBTuOXsaV0kq7si6y/nAJW/xronIFKk4fQMWp/Uj4shDKI7sQUHIUaS//BcuWLsGSJUsQGRkJDw8Pke6MiKjzYIBDJCJ3FynmjvTH3JH+MJvH4pcLxb/m7VzB0YsldmWdXN3tEpX1ZwT8cCIPgskNqampSE1Nhbu7B+bPn4d77rkHy5YtE+GOiIg6BwY4RJ2Ek5ME4wfJMH6QDI/NG4FzRWXYcbQ2b+cnQyGqTfVXVfEcPg2ew6eh2nge5Tl7UZazB59t2Yp8iS/6j5+JiGFybPlsM0aOHImRI0dCIpGIcGdERB2PAQ5RJzXI1xN3TxuGu6cNQ0lFNb49no/tRy7j62NXcLWyxq6si3wgXCYth/ek5TCVl+CMyYQVKT9B5gYcSb4T1RVlCA4OwdJfu7KmTZsGZ2d+/Imo++JPOKIuwNvdBYsnBGDxhABU1piwO7cQ2w9fQuaRyyi4WmVXVurx22JzReUm+C7WoCxnD07l7MGLL76IF198EXI/Pxz6+WcEBAR09K0QEXUIBjhEXYybsxRzRvTDnBH98NwyC/RnirDtl0vYduQSzhrL7cpKnKTwUITBQxEGi2oVqi7novzkHlRcycXaL89i4bgaVOf8hE/TPsKtt96KW265BXK5XKQ7IyJqOwxwiLowqZMEEcPkiBgmx5OLRuHYpVJsO3wJ2w5frpekLJFI4NY/GG79gwEA350swHcnC1C86zMU79qKzz77DFKpFHPmzMHSpUsxfvx4hISEiHFbREStxgCHqJuQSCQYNcAbowZ449Go4ThTWIbtRy5h2+FLyDptP7lgXT7Tf4feoYtQnrMHZSd+xM5vvoVOp8N96j9g+vTpAICLFy9iwIABHXg3REStwwCHqJsa4ueJlTMVWDlTgfzSSuw4ehnbDl/CrpxCVJnMdmWlnj7oPX4eeo+fB3NlGcoN2fgKIYhN+QlzFF54fEkEQkJCcOutt2L58uWYOHEiR2QRUacmsVga+ruu5yopKYGPjw+MRiN8fX3Frk6PZDKZcPLkSYSEhEAqlYpdnW6ltKIa3xzPx1eHL2Hn0St2Myk7UlNaiOIfNqLKsBdVVwUAwLBhw7B8+XIkJydzNFY74WdAfHwH4mrq87f+zi4uLoa392+DLPiTiaiH8aozIqu8yoRvT1zBF4cuYcfRy7hWVT/Ycfbyg9/CR2Axm1B5/ijKju/GZcMepG/5Co/8+VkM6+OMH374AWVlZZgzZw5cXFxEuCsiInsMcIh6MA9XKRaMHYAFYwegotqE708W4MtDF7H9yOV6c+1InKRwHzwW7oPHwmKJh7msGLNf+AajB3jj9Ed/wZG930Eul2PZsmWIiYlBZGQkgx0iEg0DHCICULtshGq0P1Sj/VFWWY2M73/Gz0YpMo9eQXF5tV1ZiUQCaS8ZAODIxRJUj78dPm5DUJ37I9555x2888478PX1xUcffYQFCxaIcDdE1NMxwCGietycnTBpUC/cMScE6yDB7txCfHnoIrYdvoSisup65V38BkM2bQUwbQW8jOdRdnwXzIYfcbC0FyYWV8CpsgQajQbR0dFQqVRwdXUV4a6IqCdhgENEjXKROmHW8L6YNbwvnls2FnvyjPji12Dn+lmUgdplI3ymxgJTY/Gv7BK8pd+B/gV6/PTee3jvvffg4+ODpUuXIjo6GvPmzYObm5sId0VE3Z2T2BUgoq7DWeqE6cF98I/l47Dnz1H4T8IU3DV1KPr0brhFxmIBLvopEZCQAt/Z90AqG4APPvgAS5bUrov1WzkO6CSitsMWHCJqEamTBFMUfpii8MNfF4/BHkMhtv58AV/+cgmCo24s3wC4TI4GJkfDXbiEihO7UakIwqf6c1CN9seK25ZBJpNhxYoVmD9/PruxiKhVGOAQUatJnSSYFtwH04L74O9Lx+KHnAJ8fvAith++hNLrRmMBgIusP1wm3YpTAB5LOwgXiQklJ87iSs4X+OijjyCTyXDbbbdhxYoVmD17NufaIaJm408NImpTLlIn22KgFdVj8d2JfGz9+SJ0Ry43OKlgtUUKj9vWIUC4hJqTu2DJ3YW3334b7777Ls6fP4/+/fujsrISLi4ucHJizzoR3VinCHAMBgOSkpIQFBQEAJDJZEhISGiz4wwGA7RaLfz8/FBYWIiIiAhER0e37U0QUT3uLlLMG9Mf88b0R1lVDXYeu4LPD17EzuNXUFVjrlfeRdYfLhG3ARG3IaDgLFyMOXhrbyGWhbpDt2kDXnvtNcTFxWHFihUIDQ3lchFE1CDRl2owGAwICwtDXl4eZDIZAECj0cDPzw9r1qxp9XE6nQ5arRbp6ekAAEEQEBkZiezs7AbPzaUaxMcp0sXV3s+/tKIauqOXsfXgRXx/Mh/Vphv/GHI9oUPBrjQUXbkAABg+fDhWrFiBu+66y/ZHTnfCz4D4+A7E1dqlGkQPcNRqNWQyGZKSkmzbBEGAr69vo6MqmnKcIAgIDAy0C4J0Oh1iYmJQVFTU4LkZ4IiPP1jE1ZHPv7isGtsOX8KWgxewO7cA5kZ+IlksFlRdOAaPsz+h8NC3KDEWIDU1FStXrgQAnD9/HgMHDmzX+nYUfgbEx3cgrtYGOKJ3ZqelpdX766tuMNKa49atW4fw8HDbdgCIiopqNLghoo7l4+mC2IjB+PfKyfgpMRJ/WTQK4wb6OCwrkUjgNnAUzFPuhez+tzHg9uehqwjEp/pzOHjkGAYNGoSpU6fijTfeQH5+fgffCRF1JqIGOIIgQBAEKBSKevtkMhn0en2rjsvIyIBKpQJQG/Q0dD4i6hz6ebtj5UwFtj48Azsen4VH5gZjiNzTYVmJkxSuQ8bjp3MVeCztIJa//j2GT1+IAwcP4qGHHkJAQAAWL16MTz75pIPvgog6A1GTjA0GQ4P75HI5CgsLW3WctVxKSgpiY2NhMBigUqmQlJQEpVJ5w/r5+/vXS2IsKyu74XHUeiaTCWazGSaT41E31L46w/MfJvfAHyOD8cjcIBw4W4zPDl7A/w5dgvFa/dmTAcDsHYDKGQ+i76T7ID2TBee8H/DFF18AAJYtWwYAOHv2LAICArpEd0NneAc9Hd+BuBw9f0/P+n/wNJTO0ilGUTVEEIQWH2cNbjIzM5GZmQkAUCqV0Gg0iIyMbFI3laOHdvLkyRbViZrHbDbDaDQiJyeHw4JF0Nmef28Ad4x0wYrhg7D/Qjl2Gkqx+8w1VNTU/4w6uXrAEjwT1cEzETC1CFfcq7B+8z7MCeyFu2KXo7S0FLfccguWLFmCESNGdPzNNFFnewc9Ed+BuBw9f0e/lztlgFM3N+Z6RqOxTY67vqUmKioKgiAgJSXlhkPRL168yCRjkZhMJuTk5CA4OLhL/LXd3XTm5z9yBHD7HKCsqga6o1ew5eBFfH+yADUOspOlvX1xGUDqvkK8vfcKvEbdhFK9zrbi+bhx4/C73/0OarXaLjmxM+jM76Cn4DsQl6PnX15eXq9cSUkJ5HJ5ve2iBjjWCjlqqREEocFApinHWcs0NHy0sWHiVlKplN/UInJycuI7EFFnf/5eHlIsVw7GcuVgFF6txBeHLuLT/eex/4zgsLxZIkXxmNvgPvpWBOafgOz8T8jdk4lnn30WDz30EKRSKQoLC+Hs7AwfH8dJzh2ts7+DnoDvQFxNef4N7RO1zU0mk0EmkzXYWmNNEG7JcdYyDXVzdcd5M4h6Kr/ebvj91GH47x+mY8fjs/DQnGAE+Lg7LCuRSGDuNwLG0Lshi38XI1e+hHf3XMS5ojK89NJL6N+/P+644w5kZmYy94KoCxO9UzE2Nha5ubl226z5M1FRUa06LjY2Fvv27bMrYw14Gjs3EXVdQX17Y/X8EfhBMxcfrZyMW5UD4enq+C88idQFhR6D8GLmCcxI+hrfFHhi2Mhx+OijjzBv3jwMGzYMTz75JE6dOtWxN0FErSZ6gKPRaJCRkWG3TavVQqvV2r4WBAEqlcpumHdTjktKSoJer7cbdaXRaBAdHd2kUVRE1HU5/boA6EuxE7HvySi8GDMB04L80NjqDud9J6J8/l8R/NA7mBYdj2qzBc8//zz2798PoDaZsaSkpIPugIhaQ/RRVAqFAunp6dBoNIiIiIDBYICfn59dArDRaERWVpZdl1RTjpPJZMjOzoZGo7Hl8wQFBdkFQUTU/fVyc8ZtYYNwW9ggnBfK8V/9OXyiP4+8gmsOy1f36ofzQUvhplgMZfFJXPYZhYKrlTAcOYhZs2bh1ltvxT333IO5c+cyN4OokxJ9qYbOiEs1iI9TpIurJzx/i8WC/WcFfJJ9DlsPXkBJRU2j5Z2dJBgrvQjDl6k4lL0HADBo0CDcdddduPvuuzF8+PA2rV9PeAedHd+BuFq7VIPoLThERGKQSCRQDvGFcogvnrplNHYcvYKM7LP49kS+w/WwaswWHDD3B6Kewvib8tH30h4c/e5zPP/88zh//jzee+89ALVzd3DOFCLxMcAhoh7P3UWKReMHYNH4AbhYXI5Pss8hLesczhgdz1xe7NoXxUNugfMdN2NW1SlMnDoK5VUmuEqBsWPHYvr06YiPj8ekSZPqzYZORB2Df2YQEdUxwMcDD80NwTerZ+Pj+ClYHjoQbs6Of1RKJE445abAK/pKTPqHDn98Zyec3Tzx9ttvY8qUKRg/fjxeffXVBpedIaL2wwCHiMgBJycJpgb54eW4idj7ZBSeWzYW4wc1PAFgaWUNPs+tQumCv+MmzbuYF3M3zp07h0cffRQTJkyA2WzuwNoTEbuoiIhuwMfDBXdOGYo7pwzF0YslSMs6i//uPw+hrNph+dPoCyhiMGDVrZhVfgQT+7vb8nKee+45ODk54Z577kFAQEBH3gZRj8IWHCKiZhg1wBt/XTwGe/4cidd/F4qZIX0anFunzCzFAbdxeK8oBMve2IX/7D2NDz74EE8++SSGDBmCpUuXYuvWraipaXwEFxE1H1twiIhawM1ZilvGB+CW8QE4L5Qjbd9ZbNp3FpdKKhyWP3BWwIGzAnrHvYTYmhxczvoSn3/+ObZs2YKAgADs2rULw4YN69ibIOrGGOAQEbXSQJkH/qQajofnBuPr4/n4eO8ZfH38ChzNMna1GtiDYCDsYcybrYbPuV24cmI/hgwZAgDYv38/rly5grlz53bwXRB1L83qotqwYUN71YOIqMtzljpBNdof79wTgR80c/HI3GD4e7s1WP5oqSt+8pmDgplr8Nz/jiGv4Br+8Y9/YMGCBRg5ciQ2bNiA/Pz8DrwDou6jWQGORqNBaWlpe9WFiKjbGCjzwGPzRmCXZi60vw/DrOF9G8zVKamowTu78jDnhW9gnnQnVsQ/gqtXr+LFF1/E0KFDcccdd+D48eMdewNEXVyzApyioiIolUp8+umn7VUfIqJuxVnqhPlj+uP9+ybhuyfm4ME5Qejr1XCrjt7ogh/l8xD0yAeIefTvCJ80GR999BGqqqoAAJWVlSguLu6o6hN1Wc0KcGQyGbKyshAYGIi1a9fihRde4Mq6RERNNFjuiSfmj8TutXPx1p1KzAzp02DZC6U12OumRMGsRNz/2uew+A4GAGRkZGDgwIFISEjAgQMHOqjmRF1Ps5KMrat5h4aGIjQ0FMXFxdi0aRNyc3Mxb948JsURETWBi9QJC8YOwIKxA3C68Bo27jmDTfvOori8/rw6lTVm6M4Butd3YeJgGcaYTBgxchRSU1ORmpqKmTNn4uGHH8ayZcvg4uIiwt0QdU6tmgfHx8cH8fHxWL9+PbKzs+Hn54c///nPbNUhImqioX698OebR+GnxEgk3zYeYwK8Gyx74KyAjRf8YFn6PB56JQ3RK+7Anj17EBsbi5deeqkDa03U+bV6or8NGzYgJCQEGo0GFosFZrMZmzZtQmJiInbu3NkWdSQi6vY8XKWIjRiMzx+egU9WTcPSCQPQwBJYMF6rwtaLnsgedjtuf+lzqB9LxF133QUAKCwsxMqVK6HX6zuw9kSdT7MCnIiICADAgQMHsGrVKkilUiQkJCAwMBDp6ekwGo1Yv3494uPjsW7dOlgsFubpEBE1g0QiQdhQX7wUOwEfxgzFY1EhGODj7rCs2QJ8c7YKX7lMR0JGLv67/xw+27IVb7/9NsLCwjBjxgxs2rQJ1dWOl5Qg6s4kFoujqagcc3Jygq+vLwRBgI+PDxISEqBWqxEYGNjocS+88AJWr17d6sp2lJKSEvj4+MBoNMLX11fs6vRIJpMJJ0+eREhICKRSqdjV6XH4/MVX9x1YIIHu6BV8+NMp7MppfGXyfl5umON3FWd3/Refpm9CZWUlAgIC8Nxzz+Hee+/toNp3D/wciKupz9/6O7u4uBje3r918TZ7JuPAwEAkJibitttua1mNiYioWZylTlgwtj8WjO2PnCul+PDH08jIPodrVaZ6Za+UVmJTqQvchq5Awpv3Qnrya6R/+LbdelenT5/G0KFDO/IWiDpcs7qolEolsrKymhzcxMbGws/PD4WFjf/FQURETRPczwt/WzoWP/45Ek/dMhqD5R4Oy1XWmLHlxDX81zIJN/1lE4ZMXgCz2YKTJ08iMDAQKpUKX331FZrRiE/UpTSrBUetVjfr5ImJiYiIiMATTzzRrOOIiKhx3u4uuH9GIO6ZNgy6o5fxzg952JNndFh29ykBu08JUPTJxcKhwG0xsfjvJxnQ6XQYPXo0HnvsMdxxxx1wd3ec60PUFTUrwImLi6u3LS8vD1qtFn369IFSqbSbC8c6Xw4REbUPqZME88f0x/wx/fHL+WK8sysPWw9eQLWpfsuMoeAa3igAZGPuQ+JiNYx7t+Df77+DlStXYtu2bUhLSxPhDojaR7O6qOLj4+ttCwwMxPr167F69WrbqCkiIup4Ywf64KXYidi1di4eiQyBXy9Xh+WEsmp8+EsZdN7zsPLNr/Dks+vx4IMPAgAsFgv+/ve/49ixYx1ZdaI216wA50Z9tZGRkcy3ISISWT8vdzymGo5da+ciOXo8Rvb3cliussaMTw4ZsfHaWHx8pjf0Z4qwb98+/PWvf8WoUaNwyy234Ouvv2aeDnVJzQpwJA0thfurAwcOQKfTtapCRETUNtxdpIgNH4wv/zgTH8VPxpwRfR2Ws1iArw5fwq1v7kZydjVe/vdWLF9+K7744gvMnTsXU6ZMwZYtWzq49kSt02AOzv79+6HVapGVlYW8vDwAgCAI8PPzc1heEAQAgFarbftaEhFRi0kkEkwL6oNpQX1w4nIpUr4z4LMD5x3m6WSdFpAFCYKmP4yX73wUR7Z/hPffexebN2/GkiVLRKg9Ucs0GOCEhobirbfeAgDodDrExMRAJpMhJibGYXmZTAaVSoXIyMj2qSkREbXacH8vvBAzAavnjcC7u/Pw0U9nUFpZU69cbv41vJwP9Bl0G5756G4sn9AfAFBVVYVZs2YhNjYW8fHx6N27d0ffAlGTNGkUVVRUFLKzs7F27Vpb0ENERF1Xfx93JC4chYfmBOM/e8/i7R/ycKmkol65gquVeGtvJT76uRj3TK/CZJ9rOHXqFB577DE8++yzeOihh/Dwww+jb1/H3V9EYmlyDo5CoXA4TJyIiLouL3cXxN+kwHdr5uDFmAkY4e84Ibmkogb/3HES8Z9fhvrNL/DCq2/Az88Pzz77LIYOHYqnn366g2tO1LhmJRk3ZQbjDRs2tLgyREQkDldnJ9wWNghfPToT790bgWlBjvMty6pMeG/PBWgvB+K2del48+0PMWrUKLu1goqLizuq2kQNalaA0xRMMiYi6rokEglmj+iHj+KnYMtD0zF/jL/DclUmM/6TdR4v5Mgx84kNiL7vIQBAQUEBhgwZgjvvvBNHjx7tyKoT2XEY4MTGxmL48OH1Czs5QSqVNvpPr9e3e6WJiKj9jR8kg/b34dj+p5uwPHQgpE71pwoxmS3YfOACFr3xE9QfZmHP0bOYPHkyNm7ciDFjxiAuLg6HDh0SofbU0zlMMs7OzkZeXh5KSkrslh5XKBSIiopCUFCQw5NZLBYkJSW1T02JiEgUw/298HLcRPwpajje+i4XGVnnUGUy1yu37fBlbDsMRP7uH/jwoSfwsfYVpKWlIS0tDX/605/w0ksviVB76qkcBji5ubkoLi62C26A2qHgNxpFlZWV1Xa1IyKiTmOInyeeXz4Of4wMQep3Bmzccwbl1aZ65XYcu4IdACJjnsFHf1iNTamvYsKECbb9ubm5Df6hTNRWGszB8fHxqbdtx44dNzwhW3CIiLo3f293/OWW0bVrXs0Nhre74xlHdhy7gsTvy+C37EmERi4FABw9ehQhISFYuHAh9u7d25HVph6mWUnG1qCnpKSk3r6dO3cCqF18k4iIuj95L1c8Nm8Edq2dC82CkZA3sLjnjmNXsOT1Xbj/vX3IM1ZixYoV2LZtGyZPnowlS5bgwIEDHVtx6hGaPYrqgQcegK+vLyIiIuy2BwYGIjExsc0qRkREXYOXuwtWzQ7C92vmIHFh44HOH7aeh+f8P+FT3Y+IiYnB1q1bERoaimeffbaDa03dXbMCnBdeeAFqtRqrV69GQkKC3b7AwECsW7cOL7zwQptWkIiIuoZebs5Qz6oNdNY2Eujojl7Bo9sL4L3oCWRs+x6LFy/GrFmzANQOVjlz5kxHVpu6qSYt1WBVUFCA0NBQhIaGNljGUe4OERH1HL3cnPHArCD8fspQfPjTaaR8Z4DxWlW9crqjV6A7Ciy5/W8YPLp2apLt27dj0aJFuPfee/GXv/wFQ4cO7ejqUzfRrBYciaT+HAjXMxgMLa4MERF1H9ZA5/s1c6BZMBK+ni4Oy205eAFRL32LxE9/hsWtF6ZNm4YNGzYgJCQEDz74IC5evNjBNafuoFkBTlFRkcME47oY4BARUV293JyxanYQftDMbTDQMZkt+HjvWTyiK4XqiX8h47PPERoaijfffBPBwcHIyMgQoebUlTUrwImJiUFkZCQOHjxYb19JSQni4uLqJR8TEREBvwU632vm4nHVcHi51c+SqKox4+1dp/B0lhRxz36A/6R/guHDh0OpVAIAampqUFlZ2dFVpy6oWQFOZGQk5s6di9DQUISEhGD+/PlYtWoVIiIi4OvrC4PBgNWrV7dXXYmIqBvo7eaMhyND8N2aOVDPUsDNuf6vomtVJvxzZw7WH+mNB15JR8Dg2lyct99+GyNHjsTGjRthNtefTZnIqtnDxJOSkrB9+3Z4e3sjMzMTWq0W2dnZiI+Px759+9qjjkRE1A359nJF4sJR+G7NHNw5ZQicHax1VVRWjXVfHsfcF77Bp/pzsFhq0yXuvPNOhIWFYfv27bBYLCLUnjq7Fq0mHhUVhezsbJjNZuTm5sJsNt9wCQciIiJH/L3d8dyycdj5+GzcGjoQjsazXCiuwGNpB7GlcjQ+ytyH1atX4+jRo5g/fz5UKhVOnTrV4fWmzq1FAU5d189cvGHDhtaekoiIeqAhfp54KW4ivvrjTZg32t9hmSMXS/CHT07gyqgYfPFDNu655x4cPHgQMpmsYytLnV6rA5zrabXatj4lERH1ICP6eyHlrnBsfnA6pijkDst8eyIfKz85hb6LHsWu7F9sAc7zzz+PZ555BteuXevAGlNn5DDAiY2NxfDhw+sXdnKCVCpt9J9er2/3ShMRUfc3cbAMH8dPwTv3hCOkX+96+80WIC3rHJak6vHi9uMoLqtERkYG/va3v2HEiBH48MMPmYjcgzmcyTg7Oxt5eXkoKSmBt7e3bbtCoUBUVFSDy9xbLBauJk5ERG1GIpFg7kh/3BTSF+nZ5/BS5gnkl9oPE6+oNuO1nTn4eO8Z/On1dBRlf4G/PfMM7rrrLvzzn//EK6+8gunTp4t0ByQWhwFObm4uiouL7YIbAJDJZDdMJs7Kymq72hEREQFwljrh9klDsGRCAFK/NyDlOwPKqkx2ZQquVuGpLccwesBE/Ee3F1988Dpef/113HvvvTh69CikUqlItScxNJiD42hNqR07dtzwhGzBISKi9tLLzRmPRg3HN0/Mxu8mD4GDkeU4crEE8ZuOoSri99ixOwvvv/++LbjZvHkzysrKOrjWJIZmJRn7+PjccKmG60dVERERtbV+Xu54fvk4bP/TTYga1c9hmf/9fBHxn13Aj6VylFXVQK/XY/ny5Rg9ejQ2b97M+XO6uWYFOLGxsfD19eV8A0RE1CkE9/PChrsj8O/7J2O4f/1E5MoaM/65MwdzX/gWhho5Xn75ZRQVFWH58uW4+eabceLECRFqTR2h2cPEn3jiCQwbNqwdqkJERNQyM0L64ItHZuLvS8fAx6P+Yp6XSiqw+tPD+M41Apu/3oe7774bX331FcaNG8f527qpZgU4CoUC69evb7RMXFxcqypERETUEs5SJ9w1dRi+WT0bd08dCqmDBB39GQH3pZ/EkOWr8dWObzBu3DiEhoYCqB0JzG6r7qNZAU5iYiISExPx9ddfN1jGYDC0ulJEREQt5dvLFX9bOhZf/nEmZgT3qbffYgE++PE0ntxdhUTtf20rlW/ZsgXz589HTk5OR1eZ2kGzc3B0Oh2ioqIglUrh5+eHkJAQ2z9O9EdERJ3FcH8vfHj/JKTeFY6hfp719hdcrcLj6T8jLuUnHL9UiqysLGRmZmLcuHF4/vnnUVVVJUKtqa1ILM1oj5PL5QgPD4dCoXC432KxICMjA4WFhW1WQTGUlJTAx8cHRqMRvr6+YlenRzKZTDh58qQtcKaOxecvPr6DtlVZY8KG7/Pw2s6TqKiuP7uxs5ME980IRJjrJfzxoVU4fvw4xowZgyeffBKxsbF8ByJo6mfA+jv7+vn7HE701xCFQoHt27c3WiYvL685pyQiImp3bs5SPDgnGEsmBOBvW49Ad/Sy3f4aswUp3xkwwMcdSRu/wv6t72HdunW4//77ERkZiX79HA9Fp86rWV1UqampNyzDif6IiKizGiz3xIa7w7HhrnAM8vWot/9icQX+mH4YZ4ctwlff7sazzz4LPz8/AMCZM2eYhNyFNCvAsWaat7YMERGRmKJG+yPzT7Pw8NxguErr/yrccewKHvoiHybFDJjMFly5cgUTJ07E4sWLce7cORFqTM3VrADH0SzGeXl5WLt2LV544QXs3LmzzSpGRETUnjxcpXh83gh89ehMzAypP9rqWpUJb+4pQFzKT8grLMOiRYvwv//9D2PGjME777zD1pxOrllJxnFxcdi0aVOD+3fs2IH9+/dj9erVzaqEwWBAUlKSbZVymUyGhISEdjlOpVIhMzOz0TJMMhYfEyzFxecvPr6DjmWxWPC/QxfxzJYjKLhaWW+/i1SCVbOCoKg4jgdXrcKFCxcwf/58pKSkYMiQISLUuPvr0CTjG8VCkZGR0Ol0zTklDAYDwsLCkJeXB5lMBgDQaDRITk7GmjVr2vS45OTkZtePiIi6P4lEglvGB2BmcF88/8VRbMo6a7e/2mTBP3fmIKivNz788gf8+5Vn8e6772Lbtm2Ij48XqdbUmGZ1UUkkDpZtrePAgQPNDiCSkpKQkJBgC1KA2gkFNRpNmx5nMBiwb9++ZtWNiIh6Fh9PFyRFj8fH8VMwzMHcObn513DfR0cweOlj2PHNt1i5ciUA4OzZszh9+nRHV5ca0WCAs3//fjzwwAMIDw+Hn58f/Pz8kJGRYfv/9f+kUinCwsKgVqubVYG0tDRbF5OVNWhpLFhq7nEZGRlcRoKIiJpkapAf/vfwdMSNk8HZwZIP7/94Gn/ba8K+U0WwWCy4//77MXbsWKSkpDA3p5NosIsqNDQUb731FoDagCEmJgYymQwxMTEOy8tkMqhUKkRGRjb54oIgQBAEhxMHymQy6PV6REVFtfq4jIwMREdHc5ZlIiJqMncXKe4N88PvZ43Gk5sP4+C5Yrv9pwvLEJfyI+6bHohHHv0Tjh07BrVaja1bt2LDhg3w9/cXqeYENDEHJyoqCtnZ2Vi7dq0t6GkLja1bJZfLG5wRuTnHCYIAo9EIhULR7ADH39+/XrdcWVlZs85BLWMymWA2m2EymcSuSo/E5y8+vgPxWd/B8H69kK6egvd/PI0XM0/YzYRssQBv/5CHwD6eeH/rN3j7/57Gxo0bMW7cOKSkpGDx4sXi3UAX5+gz4OlZv9uwoRazJicZKxSKDu/iEQSh1celpKQ0mqzcGEcP7eTJky06FzWP2WyG0WhETk4OnJyalSpGbYDPX3x8B+K7/h3M7AcoFg/Eiz/k48iVCruyeQVluP8/R3HbglVInqjEP577O15//XWMHDlSpNp3fY4+A45+L7c6wAGA22677YZlSkpK7IZpNaZugvD1jEZjq4+zLgzaUhcvXuQwcZGYTCbk5OQgODiYQ2RFwOcvPr4D8Tl6ByEAZoZa8O7uU3gx8ySqan5rzTFbgPRfBIT0C0NG5i6MHeyHvn37wmQy4fDhwxg/frxId9I1OXr+5eXl9cqVlJRALpfX296sAKcpYmJisG3btiaVtVbIUUuNIAgNBjJNPU6v17e49QYApFIpf7CIyMnJie9ARHz+4uM7EJ+jdyCVAupZwYgc1R+Ppx/EwbOC3TEnr1zFqs3X8MR8d6zs549XX30JiYmJeOqpp/CXv/yF77MZmvIZaGhfiwKcnTt3OgwujEYjsrKymnwemUwGmUzWYGuNSqVq8XEpKSnIzc21GzZuzcHRaDTw8/NrVfBDREQ9W3C/3vjkgalI+d6AVzJPosr0W2tOtcmC5784hm9P5CN+xlyMHTsWzzzzDL755hts3LgRAQEBIta8Z2hWgLN//36EhYVBIpE02Od1o7lyrhcbG4vc3Fy7bdYk4sa6l1pyXEpKCnQ6HRcEJSKiNuEsdcIfZgcjcqQ/Hk8/gF/O2y9ptCunEL+cd8HfN/wXme+9iDfeeAMTJkzAhx9+iAULFohU656hWZlrGo0GSUlJ2L59O3Jzc+v9y8rKcjh0+0bnzMjIsNum1Wqh1WptXwuCAJVKZTcKqinHXa+lSctERESNGdHfC5+umo5Vs4Nw/d/5xeXV+GPGEXjOjsfG/6Shuroa999/v8N8Emo7zWrBUSqVeOKJJxrcHxgY2OzVxBUKBdLT06HRaBAREQGDwQA/Pz+7NaWsXV91u6SacpyVwWCAVqu1BUQxMTFQqVRNWu+KiIioKVydnaBZMBI3hfTFY2kHcLHYfqRVWtY57PHrg4+//A6+kjJ4eHgAqE2ctf6f2k6zFtvcsGGDbVrq7oyLbYqPCw2Ki89ffHwH4mvNOyguq8afNx/C/36+WG+fi1SCJ28ehbunDUN2djYWLVqElJQULF26tK2q3i20drHNZnVRWSwWlJSUNFomMTGxOackIiLqdnw8XfD67aF4MWYCerna/3KuNlnwzNYj+MNGPS4WGGE2m7Fs2TI8+eSTnNixDTUrwImPj0dKSgo+/fTTBgMdrtZNRERUO+jmtrBB+PKPNyF0iKze/i9/uYT/+9kF//58JyIiIvD8889j4cKFKCgo6PjKdkPNysEJCQmB0Whksi4REVETDfHzRJp6Kl7Yfhzab+2XGjpjLMNDW85iTdIHmPjxS0hNTUV4eDgOHToELy8vkWrcPTQrwCksLERUVBQUCgX8/Pzq7S8oKMCGDRvarHJERETdgYvUCYkLR2FyoByPpR2EUFZt21dtsuAf23Kx8CY13giLQP6lCwxu2kCzAhyFQoG0tLRGy+Tl5bWqQkRERN3V3JH++OKRmXj44/3IPl1kt+/LXy7hRN8QaFfdDgC4evUqXnzxRSQmJsLV1VWM6nZpzcrBSU9Pv2EZTqJHRETUsACZB/6TMAXqWfXnjcvNv4Zlb+zCtsOX8Nprr+GZZ56BSqViXk4LNCvACQwMBACcOnUKn376KU6dOmXbt3PnTrsyRERE5Ji1y+qde8Ih83Sx23e1sgbqD7PhNGEJHn74YXz33XeYNGkSDh8+LFJtu6ZmBTgA8MADD0ChUCAmJsZuJmGLxcIh4kRERM0wd6Q/tj40A6MHeNfb96/vTqF44p145bU3cfbsWUydOhWff/65CLXsmpoV4Pzf//0fDAYDsrOzYTKZ7NajioyMxNq1a/HCCy+0eSWJiIi6q8FyT3yyahqWhw6st+/bE/n4pGwk3tm0Ba6urvj73/8Os9ns4Cx0vWYlGe/btw/bt2+3fX39wpo+Pj7w8fFpm5oRERH1EB6uUrwUOwETBvnguf8dRY35twaEM8YyJB10xqsff4FZowLg5FTbNmGxWJq9wHVP0qwWnIiIiBuWKS4ubnFliIiIeiqJRIJ7pgdi48rJ6NPbftTU1coaPLUzH9tOVcNiseDf//43br75ZpSWlopU286vWQFO3cUuAcDRMla5ubmtqxEREVEPNlnhh88fnomJg2V2280W4NnPj+DJzb9g1+7d+OqrrzBr1ixcvFh/vStqwVpUb7/9tu3ruk1jJSUliIuLQ1hYWNvVjoiIqAfq7+OO/yRMwZIJAfX2fbTnDIpDf48/P/VX7N+/H1OnTsXRo0dFqGXn1qwcnPXr1yM4OBjr169HVFQUDIbaKaf37dsHnU6H8PDwHrHaOBERUXtzd5Hi1RUTEdyvN17KPGG3b3euEZf6zEbyq/2R+NhDmD59OrZs2YIZM2aIVNvOp9nDxLOzszFx4kRotVpkZmZizZo1SE9Px9q1a7Ft27b2qCMREVGPJJFI8EhkCF7/XSjcnO1/ZRsKruE/xUF47b1NqKmpwdmzZ0WqZefUrBYcoHakVHp6OoqLi5GVlQWFQsHJ/YiIiNrRLeMDMMjXE/EfZCG/tNK2vfBaFV470QsfbN+DZVNGAQCqq6vh4uLS0Kl6jGa34AC1+TY+Pj6IjIy0BTfWmYyJiIio7U0cLMOWh6Zj1HWTApZVmfD4ljykZ51FUVERJk+ejH/9618i1bLzaNFMxr6+vvWGjAcGBnImYyIionY0wMcDaeopmBHcx267yWzBExk/4/Vth3D16lX84Q9/wPr160WqZefQrADnhRdegFqtxurVq5GQkGC3LzAwEOvWreNMxkRERO3Iy90F79wTgWUT64+wevtAKW7+cwrGjx+PxMRErF271uGULj1Bs3JwCgoKEBoaitDQ0AbLcCZjIiKi9uXq7ISXYifC39sd2u8Mdvs2H7uGhY+8Bo9UDZKSklBeXo5XXnmlx8163KwWnKY8HOvQcSIiImo/Tk4SJN48Ck/fMhrX/3r+8mQpRt+fhJk33YQff/wRZWVl4lRSRM0KcIqKilBSUtJoGQY4REREHee+GYF4/XYlXKT2Uc7O3FIMjPsbPvv8S/Tq1Uuk2omnWQFOTEwMIiMjcfDgwXr7rDMZN2W9KiIiImo7i8YPQOpd4fXmyvnxzDX88dMTKK2oxoYNG/D444/3mJycZuXgREZGYu7cuQgNDUVQUBAUCgUUCgWysrKg1+uhVCqxadOm9qorERERNWD2iH744L5JuP/9LFytrLFt33vKiN+l7EZRxofY9f13AGoHDXX3nJxmDxNPSkrC9u3b4e3tjczMTGi1WmRnZyM+Ph779u1rjzoSERFRE0xW+GHjysnw8bCf6O/QhavwXvokpk6fgZdeeglPP/20SDXsOC2a6C8qKgrZ2dkwm83Izc2F2WzGW2+91dZ1IyIiomaaMFiGTeop6NPbzW77kfxq+C5/GsqwcDz33HN46aWXRKphx2jxTMZWnMmYiIiocxnZ3xtp6ino7+1ut/1wfhX6x/4NI0eNwpo1a3Dy5EmRatj+OJMxERFRN6To2xsfJ0xBPy/7lpzDRgsCf78OGzelIyQkRKTatT/OZExERNRNBfbphY8TpqDvdUHOkWJnfJLfH+VVJhQWFmLXrl0i1bD9cCZjIiKibiyob298HD8ZK1J+QsHVKtv2vaeM+MPGbJx8+3EcPHgAP/zwAyZMmCBiTdsWZzImIiLq5oL7eeGj+Cnw6+Vqt/3r4/kYEHk3qqqqsGjRIpw7d06kGrY9zmRMRETUAwz398K/V06Gt7t9582+ygG45aFncf78eSxatOiGv+e7Cs5kTERE1EOMGuCNd++NgLuL/a//bJcxWHTfn/Dzzz8jNjYW1dXVItWw7XAmYyIioh4kbKgcb90ZhpXvZ6HG/NuyDYf6zMXsJRdRVnQeZWVlXT6nljMZExER9TCzR/TDS3ET7VYhl0gkOD3ydvzl9X93+eAGaGYLjpV1JmMAyMvLs032R0RERF3DkgkBKC6vxlObf7FtM0ukeOyTo0jr44NvNm/E6NGjMXv2bPEq2QpNCnBKSkrg7e1db/vatWshkUhs3VXh4eEOyxEREVHn8/spQ3G5uAKvf51j23a1sga/f+0rnHhjLdzd3JCVlYWhQ4eKWMuWcdhFtWrVKkRERCAkJASrVq1CSkqKw6zq9evXY926dQgLC8P69evh6+sLZ+cWNQoRERGRCB6fNxxLJgTYbSu09MaIuLUoKCjA8uXLUVZWJlLtWs5hgKPVaiGXy5GVlYV//etfWL16daMtM6Ghodi+fTtWrlwJi8XSYDkiIiLqXCQSCZKjxyN8qK/d9gK/CVAuS8D+/fuRkJDQ5X6/N5hknJ6e3uwkI61W2+UeABERUU/n7iJFyl3hGObnabe9YPgtGDs1Ehs3bsTbb78tUu1axmF/kkKhqNdic+rUqXrl5HJ5vXJKpbLtakdEREQdQt7LFe/eOwnL39wFoax2HhyJxAlCxEoEFV1Gv379RK5h8zhswVEoFHZf5+XlITMzEwkJCQgKCoJKpUJGRgaysrJueCwRERF1DYF9euH125VwqjN8XOrhBfeYZIyZOle8irWAwwBHJpPZfR0YGIj4+Hhs374dgYGBSEtLw+rVqzF3bv2bbcp6VURERNQ5zQjpg8SFo+y2Xa0yI+GDLGQd+BlPP/10l0hHcdhF1ViQolAoGl1NnIiIiLq2lTMD8cuFYnx24IJtW27+Ndy+ah1yftoOhUKBe+65R7wKNoHDFhxBEBo84PrWHSIiIupeJBIJ1t86HqMH2OfZVky6F759/PHII4/g9OnTItWuaRwGOEajEaWlpQ4PaKx1p7i4mKuJExERdQMerlKk3BUGmaeLbZvUwwseUQ+htLQU99xzD8xms4g1bJzDACc7OxsymQxSqbTev4yMDIfbpVIp5HI59Hp9R98DERERtYNBvp54IXqC3TaXoaEYOHUpvvnmG7zyyiviVKwJGpx22MfHB3K5vFknKywsdDjjMREREXVNUaP9sXJGIDb8kGfb5jTl9/A9dQB5eXmNHCkuhwGOUql0OAS8KcLDw1tVISIiIupc1iwYiX2njDh4rhgA4OTqjt5x/wfVPdNFrlnDHHZRNbflpi7Og0NERNS9uDo74bXblfBy+61dxMnNE09t/gVn8ovx448/ilg7xxwGONu3b2/xCdPS0lp8LBEREXVOQ/w88fyt4+y2FZdXY1rkzVCpVJ1uVFWDa1ERERER1bV4QgAWX7fyePWohbh27RpWrVrVqSYAZIBDRERETfbs0jHo5+Vm+9ojMBQ+4yLx5Zdf4j//+Y+INbPHAIeIiIiaTObpiqTo8XbbvGbfB7fePnj88ccbnEevozHAISIiomaZM6Ifbp80xPa11NMHvabfiYsXL+LVV18VsWa/aXAeHCIiIqKGPLloFL45fgUXiysAAL3Hz4On1Ix71Q+KXLNabMEhIiKiZuvt5oy/LRlj+1riJIV03CK8+u0ZEWv1GwY4RERE1CLzxvTHvNH+dts+/uEYYu6KR2Zmpki1qsUAh4iIiFrsb0vHoHedCQAtVWX49OMPsGbNGlEX42SAQ0RERC02wMcDq+cNt33t7N0PvZW34MCBA6JO/ssAh4iIiFrlzilDMcLfy/a195RoOLl54sm/PIXq6mpR6tQpRlEZDAYkJSUhKCgIACCTyZCQkNAmx+l0OmRmZkIQBBgMBsTExDTp3ERERNQ0zlInPL14NO7YsAcAIPXwhlfEchh+2Ih33nkHarW64+vU4Ve8jsFgQFhYGPLy8iCTyQAAGo0GycnJWLNmTauO0+l00Ov1SEpKAgAIgoCwsDBkZ2dDq9W2630RERH1JNOD+0A12h+ZRy4DALwjluGq/nN8lJYhSoAjehdVUlISEhISbEEKACQmJkKj0bT6OK1WaxckyWQyaDQapKSkwGAwtNk9EBEREfDnm0fBRSoBADi5esD/jmRMuH+dKHURPcBJS0uzdTFZWYMWnU7XquMyMjLqBUrh4eE3PDcRERE1X2CfXrhn2jDb1y7ygdh66DKOX+r45RtEDXAEQYAgCFAoFPX2yWQy6PX6Vh0XHR1dLwgiIiKi9vOH2cF2w8ari/Ox4JYl+OCDDzq0HqLm4DTWTSSXy1FYWNiq49LT0+vtz8rKAgBERUXdsH7+/v6QSCR228rKym54HLWeyWSC2WyGyWQSuyo9Ep+/+PgOxMd30DLe7lLcP30YXt2ZA6C2q+r0ob3469/OYcWKFZBKpU06j6Pn7+npWa+cxWJxeLzoScaNEQShzY9LSkpCUlKSw9af6zl6aCdPnmxRnah5zGYzjEYjcnJy4OQkek9qj8PnLz6+A/HxHbTcTf3NeNfNCSWVZji594aXchFO/ZSOt956C/PmzWvSORw9f0e/lztlgFM3Qfh6RqOxzY+LiYlBVFRUo6Oz6rp48SJ8fX2bVJbalslkQk5ODoKDg5sc7VPb4fMXH9+B+PgOWucPRles/+o4AMA7fAlK9v0X72/8Dx58sGmLcTp6/uXl5fXKlZSUQC6X19suaoBjrZCjFhdBEBoMZFpyXEpKCuRyebOGh0ulUn5Ti8jJyYnvQER8/uLjOxAf30HL3T0tEKnf56HwWhWkvXzRa9RN0O/diUOHDmHixIlNOkdTnn9D+0Rtc5PJZJDJZA22uqhUqjY5LiMjA4Ig2AU3Le3+IiIiohvzcJXivhmBtq+9wpYAAP7z6dYOub7onYqxsbHIzc2122ZNIm4sEbipx+n1ehiNRrtuKUEQOEyciIiond05ZahtRJVb/2AExGshmbisQ64teoCj0WiQkZFht02r1dZrbVGpVHbDxptynMFgwLp16yCXy5GRkWH7p9FompRkTERERC3n4+GCO6YMsX3tIh+IT7LPIb+0st2vLfooKoVCgfT0dGg0GkRERMBgMMDPz89uvSij0YisrCy7LqmmHBcWFgZBEOoFQgC4VAMREVEHuH96IN794RSqTGYAQNHh73Dz8n9hb+Zn9aZiaUuiBzgAoFQqoVQqG9yvUChQVFTU7OMcHUNEREQdp5+3O5ZMDEBG9jkAQOWFY8jatxXffvc9Zs+6qd2uK3oXFREREXVvd00davt/73G1ebLPv/Kvdr0mAxwiIiJqV+MHyTBhsAwA4Np3GFwHhODrLz/D1atX2+2aDHCIiIio3d01pW4rjgo1leV4891/t9v1GOAQERFRu1s0fgDkvVwBAJ6jbgKcnPH2+xvb7XoMcIiIiKjdubtIsWziQACA1L03+i5dAy/VQzCbHa8l1VoMcIiIiKhD3KocaPu/5/BpKLD0xk+Gwna5FgMcIiIi6hBjArwxsr+X7eua0kKsf+uDdrkWAxwiIiLqEBKJBLcpB9m+Nm57HVtfegLnL11p82sxwCEiIqIOs3RiAJx+ncDYc/hUwGLGq++ltfl1GOAQERFRh+nn7Y6IYXIAgLsiHACwZev/2vw6DHCIiIioQy0c2x8A4NxbDlf/IOTof0B5ZXWbXoMBDhEREXWo+b8GOADgoQiHqeIq3vnv9ja9BgMcIiIi6lADfDww8delGzyCJ8F96ETsaePh4gxwiIiIqMMt+LUVxy1gBPxXPIdcp0GwWNpu0j8GOERERNThZo/oa/f1eaEchoJrbXZ+BjhERETU4Ub4e8Hf2w0AUHHuMC5+8Dj+tfHTNjs/AxwiIiLqcBKJBDNDaltxJM5uqLp4HNt1O9vs/AxwiIiISBQ3Da8NcFz7BULi1gs5B/ehssbUJudmgENERESimBncBwAgcZLCfdBolF84gazcy21ybgY4REREJArfXq4Y4V+7+KZrwAjAXIOtO3e3ybkZ4BAREZFowob5AgDc+ocAAL7bvbdNzssAh4iIiEQTYQ1wBo5E3+V/RmlAeJvMh8MAh4iIiEQTPrR24U0nt17wHD4NpZJeyM1v/Xw4DHCIiIhININ8PWzz4VjMJlTln8a+nIutPi8DHCIiIhKNRCLBuIEyAEDJ3v/i4jsP4uvd+1p9XgY4REREJKrRAd4AABe/wQCAAz//0upzMsAhIiIiUY2xBjh9hgAADCePtTrRmAEOERERicoa4Dj79IPE2RXXLp/GBaGiVedkgENERESiGijzgI+HCyROUjjLB6G64CyOXCxp1TkZ4BAREZGoJBKJbUZjlz6DYakqw7FzBa06p3NbVIyIiIioNRR9e2HvKSP63PwoJFIXXLjKHBwiIiLq4gL79AIASKQuAABDQesm+2OAQ0RERKJT9O0NADBXXsPVnzNxSN+6uXAY4BAREZHorC045soyFH75Ki7u34GrlaYWn48BDhEREYluiNwTThJA2lsOODmjpvgyLpbWtPh8DHCIiIhIdK7OTvD3dq8dKu7dFzXCZeRfY4BDREREXdwAH3cAgNS7D0xXC5FfxgCHiIiIurgBMg8AgNRTBnPFVVwuLm/xuTgPDhEREXUKAb+24Lj2D4G5sgyXiq62+FwMcIiIiKhTGOBT24LjM/lW+Ey+FcUW9xafi11URERE1ClYc3CsCpmDQ0RERF2dX283AEC18TwKv/wnLhzNavG52EVFREREnYK8lysA62zG2+HsOwCVNWZ4SqXNPhdbcIiIiKhTsAY4Tm7WWY2vQSiratG5GOAQERFRp+Dj4QKJpE6AU3ENxmsMcIiIiKgLkzpJ4OvpWqcFpwzGsuoWnYsBDhEREXUavp4ukDi7AFIXWCqvoaiFLThMMiYiIqJOw9fTFcA1eEcsg7OsP4rLW9aCwwCHiIiIOo1ebrWhie+suwEAVytbNhcOu6iIiIio0+jtZt/2cq3S1KLzsAWHiIiIOo1ebrVz3hR+9TpqhEu4Nu3DFp2HAQ4RERF1GtYuqhrhAqou5eIau6iIiIioq7N2UUmkrrCYalrcRcUAh4iIiDoNawuOxNkVlpoqlFZwHhwiIiLq4qwBDpxdAFhwtbyiRedhgENERESdhodLbZKxROoCACgrr2zReZhkTERERJ2Gq3Nt20vv8Sq4DxmPmha2xTDAISIiok7DVVob0LgPGgMMGoMaSctCFXZRERERUafh5mwfmlTVmFt0HgY4RERE1GlYu6iKvvsAZ168FaUFF1t0HgY4RERE1GlYAxyYTbDUVKGqugsvtmkwGJCUlISgoCAAgEwmQ0JCQpsc15Jz9+nTpyW3QW3I09MTFosF5eXlYlelR+LzFx/fgfj4DsRhzcGRuLjDycMb1SZLi84jsVgsLTuyjRgMBoSFhSEvLw8ymQwAoNFo4OfnhzVr1rTquJae29XVFdXV1TAajfD19W2T+6TmcXNzs/1gkUqlYlenx+HzFx/fgfj4DsRx9GIJFr76ve1riQQwPH8zJBKJw/IlJSXw8fFBcXExvL29fztO7ABHrVZDJpMhKSnJtk0QBPj6+qKxqjXluJaemwGO+PiDRVx8/uLjOxAf34E4DPlXMffFb+22nfzHQrhIHWfVNBTgiJ6Dk5aWZus+srK2tuh0ulYd19JzExERkThcneuHJi0ZSSVqgCMIAgRBgEKhqLdPJpNBr9e3+LiWnpuIiIjE4+Zcv7Xsvd2nmn0eUZOMDQZDg/vkcjkKCwtbfFxLzw3A1n3Vr1+/en1+ly9fbvA4ajtSqRQWiwVFRUVsGhYBn7/4+A7Ex3cgDovJDHNlmd22i/mFKCryg7+/f/3yv/7Ovj71pFOMomqIIAjtdlxjZfLy8jB48GDU1NTU2yeXy1tUJ2qZvn37il2FHo3PX3x8B+LjOxDfP14B/nGDMqWlpfDx8bF9LWqAY82HccRoNLbquJaeGwACAgJw9uxZeHl5NZi1TUREROKzWCwoLS1FQECA3XZRAxxra4ij1hRBEBoMUppyXEvPDQBOTk4YNGhQo3UnIiKizqFuy42VqEnGMpkMMpmswRYVlUrV4uNaem4iIiLq+kQfJh4bG4vc3Fy7bdYE4aioqFYd19JzExERUdcm+kR/BoMBKpXKLhDRaDQICgqyLakgCAJiYmKQlJQEpVLZ5OOuL6PT6fDkk0/Czc0Nbm5uiImJqbdsQ1hYGBITE20BUEpKCgDYzXzcXktE9AQ6nQ6ZmZkQBAEGg8HhO2ir58t30LiUlBTk5ubaTYRpxc9B+2vs+fMz0PH4Pd+5tMlztHQC2dnZljVr1ljS09MtSUlJlqSkJLv9ubm5FplMZsnMzGzWcXXLPPXUU5aFCxfayhQVFVkUCoUlISHBrjwAu3/X77fWpaioyLZtzZo1dtduSpmeKDMz0+4ZOHoHbfV8+Q4cy83NtSQkJFgSEhIsMpnMsmbNGofl+DloH015/vwMiIPf851HWz3HThHgdITo6Oh627RarQWAJTc317YtISHBotVqLVqt1m573f3X/1AqKiqy1I0Vm1KmJ2rKO2ir58t3cGNKpbLBAIefg/bX0PPnZ0Ac/J7vPNrqOYqeg9NRMjIyoNFo7LaFh4cDsF+2wdrFlZCQ4HAWZC4R0XJNeQdt9Xz5DlqHnwPx8DMgDn7Pdx5t9Rx7TIATHR1d74E1RBAE6HS6ess5cImI1rnRO2ir58t30Db4Oeh4/AyIi9/z4mvL59hjApz09PR6CUpZWVkA7EdUZWZmQqfT2VoWVCqV7YG29xIR3d2N3kFbPV++g9bj50Ac/AyIh9/znUNbPsdOvVRDe0tKSkJSUpJdpKjVam1fK5VKqNVqxMTE1Btu7khrl4joieq+gxtF5m31fPkOboyfg86Jn4H2w+/5rqE5z7HHtOBcLyYmBlFRUXZDAAHUaxZTKpUwGAzQ6XTtvkRET3P9O2ir58t30Hr8HIiDnwHx8Hu+c2jL59ilWnD0ej3i4+ObXD41NdU2b05dKSkpkMvl0Gq1dts1Gg3i4uLsjrEu+WAwGGxNl+21RERX0J7voK2W4OA7sNfQO2gIPweNa8/nz89Ay7XmvfB7vvNoy+fYpQIcpVKJ7OzsVp0jIyMDgiDY/WK1PrTk5GQEBQXZfZNbI0aFQgEZl4ho93fQFs+X76B1+DloXHs+f34GWq4174Xf851HWz7HHtVFpdfrYTQa7bqlrFnzQG0+yPVJsNbmyeYs/8AlIhp2o3fQVs+X76Dl+DkQFz8DHY/f851Lmz3HtpycpzPLzc21REdHW9LT0+3+JSQkWLKzsy0WS+1Mu+np6bZjrDPt1t2Wm5trUSgUdudes2aNRavVNqtMT9SUd9BWz5fv4MYczeRtsfBz0FEaev78DHQ8fs93Lm31HEVfi6qj+Pr6Nph9XfcRWNdKAmojRrVaXS9i1Ov12LRpEyIiImxR5fXJyk0p09M09R201fPlO6hPEASsW7cOgiAgJSUFMpkMsbGxCAoKsns2/By0j6Y+f34GOh6/5zuXtniOPSbAISIiop6jR+XgEBERUc/AAIeIiIi6HQY4RERE1O0wwCEiIqJuhwEOERERdTsMcIiIiKjbYYBDRERE3Q4DHCIiIup2GOAQUbem1+vFrgIRiYAzGRNRtxYUFFRv4T4i6v7YgkNE3VZGRkaDqw/HxMQgLCwMEokEEonEtt7NjajVakgkEvj6+kKlUiElJaUtq0xEbYQBDhF1GQaDATqdrsnltVot1Gq1w33p6elIT0+3BUBNCXB0Oh2MRiMAYMeOHcjMzERCQkKT60NEHYcBDhF1GWq1GklJSU0qKwgCDAYDlEplg2V0Op0tAGpKgGMwGCAIAmQyWaPnJSLxMcAhoi5Dp9NBpVI1qWxaWhqio6MbLZObm2src6M8nZSUFMTGxkKn0zXY7UVEnQcDHCLqEqxdU00NLtLT0xvsnrqeTCZrdLSVwWCAQqGwtfI0NcgiIvEwwCGiLiEzMxMAmtQ1JAgCjEYjFApFo2WCgoIAwC54ccSarNzcIIuIxOMsdgWIiBqi0+mg1WoB1AYZMpkMMTExAGrzcRoKNFJSUhAXF3fDc1uPVygUyMjIcFguIyPDlkicmZkJmUzWaOBERJ0DAxwi6rSioqJsQYhEIkFiYiLWrFlzw+M2bdqEHTt2NFpm3759tvybiIgIZGRk2BKIrQRBAADbNp1Od8O8HiLqHNhFRUSdXnO6hqxdTXUDlRuxtshc302VkpJiC2isOTrMvyHqGhjgEFGn15z8m8bmvrESBAF+fn62rx0FONePlmL+DVHXwi4qIur0mjM0OyMjA9nZ2c06n6MAx2Aw2JVh/g1R18IWHCLq1ARBgF6vtyUXN0av10OhUNywe2rfvn12rUHW8ta5cFJSUurNUMz5b4i6FgY4RNSpOeoaami5hqZ0TzXEOlTcGiTVxfwboq6HAQ4RdWr79u2z6xoSBMGWk3O9poxyMhgMdvk3VgqFAllZWQ5baph/Q9T1MMAhok7Nz8/PrkVl3bp1SExMrFdOp9M1KQlZo9E47MJSKBQQBMHh4pmbNm2ylSGirkFisVgsYleCiKghgiAgJibGFrwkJiY6DFDUajViYmIabGXRaDTIyMiAwWCATCZDeHi4XUtQSkoK5HK5rQVIEATEx8dDEAS7FhyVStWkuXiISFwMcIioWwgKCrrhgplE1HOwi4qIujzrWlFERFYMcIioy9u0aVOLR08RUffELioi6vLYPUVE12MLDhF1aXXXiyIismKAQ0RdWnZ2NruniKgedlERERFRt8MWHCIiIup2GOAQERFRt8MAh4iIiLodBjhERETU7fw/jEx1ogHlPbQAAAAASUVORK5CYII=", "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 gwEccEvNS$')\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": "00988450-8e6c-4089-a326-129ecb110058", "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 }