diff --git a/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb b/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb index 2c1c7d294..fbe6d8246 100644 --- a/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb +++ b/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb @@ -81,8 +81,8 @@ { "data": { "text/plain": [ - "[,\n", - " ]" + "[,\n", + " ]" ] }, "execution_count": 6, @@ -91,7 +91,7 @@ }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEGCAYAAABy53LJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAvVUlEQVR4nO3dd3hUZdrH8e+dAqGEFhIIhCZF6RBCQEAQEaWoCChNkSqLiorIuuzq+trWuq66KxZkUVQwIgqCIFhAERElIIo0pUogEDqEEjKZ+/1jRjZigEzIyZTcn+vKlTnnPM/ML2Vy57TnEVXFGGOMya8wfwcwxhgTXKxwGGOM8YkVDmOMMT6xwmGMMcYnVjiMMcb4JMLfAQpT5cqVtXbt2v6OYYwxQWPlypX7VDXWlz4hVThq165Namqqv2MYY0zQEJHtvvaxQ1XGGGN8YoXDGGOMT6xwGGOM8Ymj5zhEpBvwAhAOTFbVJ8/Y3gt4FHADLmCsqi71btsGHAVyAJeqJhUkQ3Z2NmlpaZw8ebLAX0cgioqKIiEhgcjISH9HMcYUM44VDhEJByYCXYE0YIWIzFHVdbmafQ7MUVUVkWbADOCSXNs7q+q+C8mRlpZGdHQ0tWvXRkQu5KkChqqyf/9+0tLSqFOnjr/jGGOKGScPVSUDm1R1i6qeAlKAXrkbqGqm/m+UxTJAoY+4ePLkSWJiYkKmaACICDExMSG3F2WMCQ5OFo7qwI5cy2nedb8jIr1FZAMwDxiea5MCn4jIShEZdbYXEZFRIpIqIql79+49W5uC5A9oofg1GWOCg5OFI6+/bH/Yo1DVWap6CXA9nvMdv2mvqolAd+AOEemY14uo6iRVTVLVpNhYn+5hMcaY4JXjggV/hcNpRf7SThaONKBGruUEYNfZGqvqEqCuiFT2Lu/yfs4AZuE59BVw2rVrl+f6oUOHMnPmzCJOY4wpNr54HJa/BNuXFflLO1k4VgD1RaSOiJQABgBzcjcQkXriPeYiIolACWC/iJQRkWjv+jLAVcBPDmYtsGXLiv6HZowp5n75FL56FhJvgWb9ivzlHbuqSlVdIjIGWIjnctwpqrpWREZ7t78C9AVuEZFs4ATQ33uFVRVglremRADTVXWBU1kvRNmyZcnMzERVufPOO1m0aBF16tTBZlY0xjjicBp8MAqqNIHuT+PKcXPS5aZsyaIbQcrRV1LV+cD8M9a9kuvxU8BTefTbAjR3MlthmzVrFhs3bmTNmjXs2bOHRo0aMXz48PN3NMaY/MrJhveGeT7fOBUiS/H43HUs3bSXWbe3p0wRFQ+7c7yQLFmyhIEDBxIeHk61atW44oor/B3JGBNqFt4Pad/Bdf+GyvWY/u2vTPl6Kx3qxRZZ0QArHIXKLpE1xjjmhxT47lW4dAw06cOyTft48MOfuPziWP7W45Lz9y9EVjgKSceOHUlJSSEnJ4f09HQWL17s70jGmFCR/gPMvRtqXwZXPszWfce4bdoq6lQuw78HtiQivGj/lIfUfBz+1Lt3bxYtWkTTpk1p0KABnTp18nckY0woOH4A3r0ZSsfADa9zKMvNiDdWECbw3yGtKRdV9OPVWeG4QJmZmYDnMNWLL77o5zTGmJDizoGZw+Hobhi2gFNRMYye8i1pB0/w9sg21Iwp7ZdYdqjKGGMC1aJHYcti6PksWj2RB2avYfmWAzx1Q1OS61TyWywrHMYYE4jWfQhLn4NWQyHxFl75cgszUtO4q0t9erdM8Gs0KxzGGBNoMjbA7NshoTV0f5r5a9J5asEGrm1ejXuurO/vdFY4jDEmoBw/AO8MgMjS0O9NVqef4J53V5NYswLP3NAsIC77t5PjxhgTKHKy4b0hcGQnDJ3HTndFRk79mtjokky6JYmoyHB/JwRsj8MYYwLHwvth6xK49gWOxrZkxBsryMrO4fWhralctqS/051mhcNhO3bsoHPnzjRs2JDGjRvzwgsv+DuSMSYQrXzj9J3hrqYDuPOd7/klI5OXbk6kfpVof6f7HTtU5bCIiAieffZZEhMTOXr0KK1ataJr1640atTI39GMMYFi+zKYNx7qdoErH+axeev5YuNe/tG7CZfVD7wJ6myPw2Hx8fEkJiYCEB0dTcOGDdm5c6efUxljAsahXz13hlesBTdMYeq3abyxbBsjO9Thpja1/J0uT8Vqj+PhuWtZt+tIoT5no2rl+L9rG+er7bZt2/j+++9p06ZNoWYwxgSprEx4Z6BnGtiBKSzefoqH567lyoZV+GuPhv5Od1a2x1FEMjMz6du3L88//zzlypXzdxxjjL+53TD7NshYBzdOYX12FcZMX0XD+HK8MKAF4WH+v+z2bIrVHkd+9wwKW3Z2Nn379uWmm26iT58+fslgjAkwXz4F6+fAVf8gPbY9wyYuo2xUBP8d0rpI59YoCNvjcJiqMmLECBo2bMi4ceP8HccYEwjWzoYvn4QWN3G05SiGvb6CzCwXrw9Npmr5KH+nOy8rHA77+uuveeutt1i0aBEtWrSgRYsWzJ8///wdjTGhKf1HzyGqhGSyuz/L7dO/Z1NGJi/dlEijasFxGNvR/SER6Qa8AIQDk1X1yTO29wIeBdyACxirqkvz0zdYdOjQAVX1dwxjTCA4ku4ZTqRURbT/W/x1zs989cs+nr6hGR0bBN5lt2fj2B6HiIQDE4HuQCNgoIicefPC50BzVW0BDAcm+9DXGGOCx6ljnqJx4hAMTOGFb48wc2Uad3epT7+kGv5O5xMnD1UlA5tUdYuqngJSgF65G6hqpv7v3/EygOa3rzHGBA23Gz4YBbt/hBum8N7Oijz/2S/0TUxgbACMdusrJwtHdWBHruU077rfEZHeIrIBmIdnryPffY0xJih8/hBs+Aiufpyvwlrx1w/W0KFeZZ7o0zQgRrv1lZOFI6/vxh8O9qvqLFW9BLgez/mOfPcFEJFRIpIqIql79+4taFZjjHHGyqnw9QvQeiTraw7itrdXUS+uLC/dnEiJiOC8PsnJ1GlA7gN3CcCuszVW1SVAXRGp7EtfVZ2kqkmqmhQbGzwnl4wxxcCWL2DeOKjbhfR2DzHsjVTKlozg9WGtKRcV6e90BeZk4VgB1BeROiJSAhgAzMndQETqiXc/TUQSgRLA/vz0NcaYgLZ3I7x7C8TU58i1rzFs6veeezWGtSa+fCl/p7sgjhUOVXUBY4CFwHpghqquFZHRIjLa26wv8JOIrMZzFVV/9cizr1NZnTZ8+HDi4uJo0qTJ6XUHDhyga9eu1K9fn65du3Lw4MHT25544gnq1avHxRdfzMKFC/0R2RhzIY7tg+n9IKIE2QNSuH3mJjZlZPLyzYk0jA+OezXOxdEDbKo6X1UbqGpdVf2Hd90rqvqK9/FTqtpYVVuo6qW/3cNxtr7BaujQoSxYsOB365588km6dOnCL7/8QpcuXXjySc9tKuvWrSMlJYW1a9eyYMECbr/9dnJycvwR2xhTEK4sSLkJju5GB0znr4sOs3TTPp7o0zQgh0gviOA8MxNkOnbsSKVKlX637sMPP2TIkCEADBkyhNmzZ59eP2DAAEqWLEmdOnWoV68e3333XVFHNsYUhCp8OAZ2LIfrX+bZdeVP36txY5Ddq3EugT2SVmH7eALsXlO4z1m1KXT3/ab2PXv2EB8fD3jm7MjIyABg586dtG3b9nS7hIQEm7/DmGDx5VOwZgZc8XemHknkxcVrGZhcIyjv1TgX2+MIMHkNTxKM13kbU+x8Pw2+eAKaD2Je+UE85J1X49FeTULuPVy89jgKsGfglCpVqpCenk58fDzp6enExcUBnj2MHTv+d+9jWloa1apV81dMY0x+bPoc5t4FF13O8sYPcs/U1bSqWZEXB7UkIjz0/j8Pva8oSFx33XVMnToVgKlTp9KrV6/T61NSUsjKymLr1q388ssvJCcn+zOqMeZcdq+BGUOg8sVs6DiRW6f9SK2Y0kwekkRUZLi/0zmieO1x+MnAgQP54osv2LdvHwkJCTz88MNMmDCBfv368d///peaNWvy3nvvAdC4cWP69etHo0aNiIiIYOLEiYSHh+YvnzFB73AaTLsRSkazq+ebDH57A2WjIpg6PJkKpUv4O51jJJSG/E5KStLU1NTfrVu/fj0NGwbu3L0XIpS/NmMC3olDMKUbHNnJoQFz6PP+YfZlZjHztnY0qBLt73T5JiIrVTXJlz62x2GMMb5yZcG7N8P+TZzs/y5D5h9n56ETTBvZJqiKRkHZOQ5jjPHFb/dqbPsK17X/ZvSyaNakHeLFQYkk1a50/v4hoFgUjlA6HPebUPyajAkKnz8Ca2agnR/gvl8a8cXGvTzeuyldG1Xxd7IiE/KFIyoqiv3794fUH1pVZf/+/URFBf6k9saElNQpsPRf0GooTx7ryQerdjKuawMGJNf0d7IiFfLnOBISEkhLSyPU5uqIiooiISHB3zGMKT42LoB590L9q5lSfgyvzv+ZwW1rcecV9fydrMiFfOGIjIykTp06/o5hjAlmaakwcxhUbcaH9R/lkQ9+pnuTqjx0XeOQuys8P0L+UJUxxlyQvT977tUoG8eXSRMZN3sz7evF8Fz/FoSHFb+iAVY4jDHm7A7vhLf7QFg4qzq9zq2zdtCkenleHRy6d4XnR8gfqjLGmAI5cRDe7gsnDvFzjxRumbWXWpVK88bQ1pQtWbz/dBbvr94YY/KSfQKmD4ADm9nZ8y36zzlOhdKRvDWiDRXLhO5QIvllhcMYY3LLccF7w2DHt+zv8Sp9F0QSHqa8PaINVcvbJfBg5ziMMeZ/VOGju+Hnjzna5QluWFKF46dcvDUimdqVy/g7XcCwwmGMMb/5/BH4/m1OtruX/quasPvwSV4f1pqG8eX8nSygOFo4RKSbiGwUkU0iMiGP7TeJyI/ej2Ui0jzXtm0iskZEVotI6pl9jTGmUC1/GZb+C1eLIQze3IVfMo7yyuBWtKpVPMaf8oVj5zhEJByYCHQF0oAVIjJHVdflarYV6KSqB0WkOzAJaJNre2dV3edURmOMAWDNTFgwAffF1zDqwEBSfz3AiwMT6dQg1t/JApKTexzJwCZV3aKqp4AUoFfuBqq6TFUPeheXAzaGhjGmaG36DGaNRmu2Y6xrDIt+PsA/rm9Kz2bx/k4WsJwsHNWBHbmW07zrzmYE8HGuZQU+EZGVIjLqbJ1EZJSIpIpIaqiNR2WMcdj2byDlZjTuEu6P+htz1h7ggZ4NGdSmeA1a6CsnL8fN6178PIeoFZHOeApHh1yr26vqLhGJAz4VkQ2quuQPT6g6Cc8hLpKSkkJnCFxjjLPSf4Dp/dDy1Xmi0uNMX3WEcV0bMPKyi/ydLOA5uceRBtTItZwA7DqzkYg0AyYDvVR1/2/rVXWX93MGMAvPoS9jjLlw+36Bt/qgJcvxn2rPMGnVUUZ3qlssR7otCCcLxwqgvojUEZESwABgTu4GIlIT+AAYrKo/51pfRkSif3sMXAX85GBWY0xxcWgHvHk9iPB63ef414oTDG1Xm790u7hYjnRbEI4dqlJVl4iMARYC4cAUVV0rIqO9218BHgRigJe8PzCXd9L0KsAs77oIYLqqLnAqqzGmmMjMgDd7QdZRZjR9lUeWZtMvKYEHr2lkRcMHEkoz4yUlJWlqqt3yYYzJw4mD8Ma1cGAz81q8xB1fleC65tWK9fDoACKy0vsPe77ZWFXGmNB36hhM6wf7NvJFq/9wx5ISdG1UhWf7NS/WRaOgrHAYY0KbKwtSboKdqXzX+jmGfRVNxwaxvDioJZHhNupSQVjhMMaErhwXvD8CtixmTdITDFwaR+vaFXn15laUjCi+EzFdKCu3xpjQ5M6BD++A9XPZ0OJv9F5Wm2YJ5ZkytDWlSljRuBBWOIwxocftho/Gwo8pbGpyD9euaErj6uWZOjy52M/eVxjsO2iMCS2q8PGfYdWbbGt0Oz2+b0PD+GjeHJ5MuahIf6cLCbbHYYwJHaqw8H5YMZkdDUdy1Q+X0aBqWd4a3obypaxoFBYrHMaY0KDqmYhp+UR2XnwLXdZ0oV5cNG+PaEP50lY0CpMVDmNMaPjyaVj6L3bXG8AVa7tzUeWyTBvZhgqlS/g7WcixwmGMCX5Ln4MvHiejbl86b7yO2jGeolGxjBUNJ9jJcWNMcPvmJfjsIfbVvpbOP99AQsUyTLu1DTFlS/o7WciyPQ5jTPBaMRkW/pUDtbrRectAqlYozfRb21LZioajrHAYY4LTqrdg3r0cSOhC562DiS1XhndubUtstBUNp1nhMMYEn1VvwZw7ORh/GZdvH0JshWhS/tSWuHJR/k5WLFjhMMYEl1VvwpwxHIjvQMcdtxJfqSIpo9oSF21Fo6jYyXFjTPBY+QbMvZv98R3ptGMkNWIrMW1kGyrZ1VNFyvY4jDHBIfV1mHs3+6p2pNOvI6kVV4npVjT8wgqHMSbwpU6Bj8aSUbUTHX8dQd34GKaPbGv3afiJHaoyxgS2Ff+FeePIqNKJTr+O4JKEyky1AQv9ytE9DhHpJiIbRWSTiEzIY/tNIvKj92OZiDTPb19jTDGwYjLMG8fuKp3otGMEjWvE2ii3AcCxwiEi4cBEoDvQCBgoIo3OaLYV6KSqzYBHgUk+9DXGhLLvXoN595JepROX/zqcZrXimDo8mWgrGn7n5B5HMrBJVbeo6ikgBeiVu4GqLlPVg97F5UBCfvsaY0LYd6/B/PGkxXqKRuJFVXl9WGvK2CRMAcHJwlEd2JFrOc277mxGAB8XsK8xJlQsexHmj2db5U503jGCDhdXZ8rQ1pQuYUUjUDj5k5A81mmeDUU64ykcHQrQdxQwCqBmzZq+pzTGBAZVWPJPWPwYG2O6cE3aEK5uVoPn+rcgMtwuAA0kTv400oAauZYTgF1nNhKRZsBkoJeq7velL4CqTlLVJFVNio2NLZTgxpgipgqfPwyLH+OHSt3osXMovZNq88KAllY0ApCTP5EVQH0RqSMiJYABwJzcDUSkJvABMFhVf/alrzEmRKjCggmw9Dm+rXQd1++6mVva1+XJPs0ID8vr4IPxN8cOVamqS0TGAAuBcGCKqq4VkdHe7a8ADwIxwEsiAuDy7j3k2deprMYYP3HnwEf3wKqpLK54A8N29ebOK+ozrmsDvH8TTAAS1TxPHQSlpKQkTU1N9XcMY0x+5Lhg9m2wZgYflR/EmD09mdC9IaM71fV3smJFRFaqapIvfewyBWNM0XOdgvdHwPo5zCg3hPv2XM2j1zdhcNta/k5m8sEKhzGmaGWfhBm3wC8LmVJ2FI/tvZx/9WtOn8SE8/c1AcEKhzGm6GRlQsogdOsS/h11OxMPduSlm1rQrUm8v5MZH1jhMMYUjeMHYNqN6K5VPBZ5JynH2/PGsCTa1avs72TGR1Y4jDHOO7IL3uqN+8BWxst4vnQnkzIqmaYJ5f2dzBSAFQ5jjLP2b4Y3r8d1/AAjs//CpjItmTmiDXUql/F3MlNAVjiMMc5J/wHe7ktWtosBx//G8cpNeX9EMlXK2fzgwSxfd46LyIgzlsNF5P+ciWSMCQnbvoY3riEzJ4IemfcTntCSGX+61IpGCMjvkCNdRGS+iMSLSBM8Q6BHO5jLGBPMNn6Mvt2H/WExdD30N2o3aMFbI9pQvrTNpREK8nWoSlUHiUh/YA1wHBioql87mswYE5xWv4N+eAc7SzXg2gN3c0ViI57s29QGKwwh+T1UVR+4G3gf2AYMFpHSDuYyxgSjb16C2aPZGNWcqw+M58aOLfjnjc2saISY/J4cnwvcoaqfi2fksXvwjGDb2LFkxpjgoQqfPwJL/8U3Jdsz5OAo7uvZjJGXXeTvZMYB+S0cycBIEbkDz4RKS/EMdW6MKe5cp2DOnfBjCnMirua+Y0N4blAiPZvZ3eChKr+FYzJwFPiPd3kgcCnQz4lQxpggcfIIzBgMW75gYthAJrl689bI1rSuXcnfyYyD8ls4LlbV5rmWF4vID04EMsYEiSPpMO1G3Bnr+Zv7dpaW6sr7w5KpF1fW38mMw/JbOL4XkbaquhxARNoAdlWVMcXV3o3wdl+yM/czIms8B+I78MHQ1sRF2z0axUF+C0cb4BYR+dW7XBNYLyJrAFXVZo6kM8YEnu3foO8M4FhOOP2P309sg2TeHZRImZI2EEVxkd+fdDdHUxhjgsO6D9H3byUjvAp9M++lQ1Iij13fhAi73LZYye8NgNudDmKMCXDfvop+/Bd+jmxI/yN3M6JrK8ZcUc/mBi+GbN/SGHNu7hz49EH45kWWhrfh9uN38Ej/VvRuaTP2FVeO7l+KSDcR2Sgim0RkQh7bLxGRb0QkS0TGn7Ftm4isEZHVIpLqZE5jzFmcOgbvDoZvXuQdunGPjuP1Wy+zolHMObbHISLhwESgK5AGrBCROaq6LlezA8BdwPVneZrOqrrPqYzGmHM4sgveGYA7fQ2PuobwVaW+fDCkNTVjbLSh4s7JQ1XJwCZV3QIgIilAL+B04VDVDCBDRHo6mMMY46v0H9DpAzh17BB/OnUvOXW78v6gRMqXstFtjbOHqqoDO3Itp3nX5ZcCn4jIShEZdbZGIjJKRFJFJHXv3r0FjGqMOW3jx+iU7hw87qLXiQdJSO7F60NbW9Ewpzm5x5HXpRbqQ//2qrpLROKAT0Vkg6ou+cMTqk4CJgEkJSX58vzGmNxUYflL6ML7+Tm8HrccH8tt17RnSLvaduWU+R0nC0caUCPXcgKwK7+dVXWX93OGiMzCc+jrD4XDGFMIcrLh4/sgdQqLpS33ue7gmSFt6XxJnL+TmQDkZOFYAdQXkTrATjyj6Q7KT0cRKQOEqepR7+OrgEccS2pMcXbyMLw3FDYvYpL7Ot4sNYS3hyVzSdVy/k5mApRjhUNVXSIyBlgIhANTVHWtiIz2bn9FRKoCqUA5wC0iY4FGQGVglnf3OAKYrqoLnMpqTLG1fzP6zkDc+zcxIXsU22v0YfbNiVQuW9LfyUwAc/QGQFWdD8w/Y90ruR7vxnMI60xHgOZ5rDfGFJbNi9AZQ8nMVkZlTaBucnfevqYxJSJs+BBzbnbnuDHFzW8nwT95gK0kMOLUvYy6/goGJtf0dzITJKxwGFOcuLLgo3tg9TQ+12QejriL525tT5JNvGR8YIXDmOLi6G405WZk5wqed/VhcZVhzLilNfHlS/k7mQkyVjiMKQ52rsT9ziCyjx3i7lNjKdW8N+/2aUpUZLi/k5kgZIXDmFC3+h3cc+8mw12e4Vn/R5/u3RjRoY7d1GcKzAqHMaHKdQoW/hVWTGalNuK+sHE8MuxyLqsf6+9kJshZ4TAmFB3Zhc64BUlbwWuuHsyvMpq3bm5NQkUb2dZcOCscxoSabV+TM2MI2SeOMv7UnVRI7k/KNY0oGWHnM0zhsMJhTKhQheUvo588wE6N4/acRxl+Qw/6JNqkS6ZwWeEwJhScOobOuQv5aSafuZN4vuw9PHtLRxtvyjjCCocxwW7fJnLeHYzsXc8z2f3YcvGtvNOvJeWibP4M4wwrHMYEszUzyZlzF5muMO7KnkD7q2/k5csuskttjaOscBgTjLJPogv+iqycwmq9mAcj7+WB4Vdyad0YfyczxYAVDmOCzf7N5MwYQvieNbziupbltW/jjf5JxEbbUOimaFjhMCaYrJ1FzuwxHMuGcdl/JrHrAKZ0rEtYmB2aMkXHCocxwcCVhS68H1nxGj+66/NI1HjuH3qVjWpr/MIKhzGBbv9mXDOGEbHnB15z9SC13l283q8VFUqX8HcyU0xZ4TAmUKnC6unkzBvPMVcYf3HdS3K3wbzSvrZdNWX8ygqHMYHo5GFy5owlfN0HrHA35Nmy4/n7oCtpllDB38mMscJhTMDZsYLsGcMIO7qTZ7L7cTDxdqZe25TSJeztagKDo7PSi0g3EdkoIptEZEIe2y8RkW9EJEtExvvS15iQ485Bv3wG95Sr2XPkJMPlEZoOfJTH+7a0omECimO/jSISDkwEugJpwAoRmaOq63I1OwDcBVxfgL7GhI7DO8meeSuRO75mTk475tX8M0/3b0eVclH+TmbMHzj5b0wysElVtwCISArQCzj9x19VM4AMEenpa19jQoIqrHkP19xxZGdn84DrNupfdSuvdLjI7s0wAcvJwlEd2JFrOQ1oU9h9RWQUMAqgZs2avqc0xl+OHyB7zlgiN3zIancD/lNuHH8Z1JNG1WxEWxPYnCwcef27pIXdV1UnAZMAkpKS8vv8xvjXL59y6oPbkRMHeNrVH/eld/HqVQ2JirTJlkzgc7JwpAE1ci0nALuKoK8xgSsrE9fCB4hY9Tpb3DX4Z+mnGT2gt90BboKKk4VjBVBfROoAO4EBwKAi6GtMYNrxHVkzRhJ59FdedfUko9W9/Ltnc7tiygQdx35jVdUlImOAhUA4MEVV14rIaO/2V0SkKpAKlAPcIjIWaKSqR/Lq61RWYxyVfQLXon8Q9s1EMjSGJ0o+yk03D6R9vcr+TmZMgYhq6JwWSEpK0tTUVH/HMOZ/fl3OyZmjiTqylemuzmxo+hf+3CuJaJudzwQIEVmpqkm+9LF9ZGOccOoY2Z88TETqJPZqZZ4p8TB9Bg5i0MVx/k5mzAWzwmFMYdu6hBPv306pzB284bqanYnjebxnImVL2tvNhAb7TTamsGQdJevjByi5+g12u6vwfJnHubn/QIbaFVMmxFjhMKYQ6Ib5nPxwHCVP7GZyTg+OXvoXnura1O7LMCHJCocxF+LILo5/eC+lN89nu7sGkyo8w4gB/Whcrby/kxnjGCscxhSEOwfXt6/h/uwRxJXNszqICl3G8nSH+kSEOzrotDF+Z4XDGF+l/0DmzDsou38NX+Y049M693F7nyupVqGUv5MZUySscBiTX1mZnPjkUUqunMQJjeaZEuPo1P9PPNawqr+TGVOkrHAYcz6quH/6gJMfTaB0VgbTc7qwt80EJlzVilIl7OS3KX6scBhzLnvWcWTWPZTbvZwt7tqkxP6LW/r1o0GVaH8nM8ZvrHAYk5eTh8lc+Bilvp9Mjpbm6cg/0fCaO3m0eQIiNsGSKd6scBiTm9tN9qppZC98kNKnDvKuduFgm/sY07WVjWJrjJe9E4zx0rRUDn9wLxUOrGaNux4Laj7O4D69qFGptL+jGRNQrHAYc+hXDs39OxU2zyZby/N0qbtp12cMf2tgAxIakxcrHKb4OnmEI589RamVrxLlhslhfSl1xb3c074hkXYTnzFnZYXDFD85Lk58OwX3oscp5zrIh+4O7Gp1Hzdd3Y5yNk+GMedlhcMUH6pkb/yEzLkTqHhsC9+5L2Fp3ScZ0KuX3fVtjA+scJhiwb19OQfm3E/l/akcdlfhjdgH6dp7BOMSKvg7mjFBxwqHCWm6ew17P/w7cemLUS3PxNKjaXzNnYxtVN3uxzCmgBwtHCLSDXgBCAcmq+qTZ2wX7/YewHFgqKqu8m7bBhwFcgCXr3PimuJN929h75z/o/L2uZTUUrxWcjDVrx7LbS3rEhZmBcOYC+FY4RCRcGAi0BVIA1aIyBxVXZerWXegvvejDfCy9/NvOqvqPqcymhB0ZBcZ8x6j0sYUojWcaRG9KdvlXoa1aWTDnRtTSJzc40gGNqnqFgARSQF6AbkLRy/gTVVVYLmIVBCReFVNdzCXCUWH09jz8ZNU3JBCBXXzQXhXpON4+nVIpGSEDURoTGFysnBUB3bkWk7j93sTZ2tTHUgHFPhERBR4VVUn5fUiIjIKGAVQs2bNwklugoYe3M6e+U8S88t7VFI3c8M6c6Lt3fS5or2NXGuMQ5wsHHkdSFYf2rRX1V0iEgd8KiIbVHXJHxp7CsokgKSkpDOf34QoPbCV9HmPE7f5fSopzAnvguvSsVzXqa0VDGMc5mThSANq5FpOAHblt42q/vY5Q0Rm4Tn09YfCYYoXd8bPpM9/gqrbZhOjYXwYcRV0uIeeHZKIirSCYUxRcLJwrADqi0gdYCcwABh0Rps5wBjv+Y82wGFVTReRMkCYqh71Pr4KeMTBrCbAZW1bzt6Pn6bankXEaASzIntQotNYrr20FSUi7KS3MUXJscKhqi4RGQMsxHM57hRVXSsio73bXwHm47kUdxOey3GHebtXAWZ5r7OPAKar6gKnspoA5XaT+dM8Dn/2LNWPfE9ZLcO7pfpRodMdXJ/c1K6SMsZPxHNBU2hISkrS1NRUf8cwF8qVxYHl08hZ+gKxJ7eRppX5ouKN1L36NtpeUtNu3DOmEInISl/vk7M7x03A0KN72Pn5y5Rd8yaVcvazXmuxMOHvJPUYxs3VY/wdzxjjZYXD+N3Jbd+R/skLVN+1gARcfE1zdlz8dy7vPoCbbfBBYwKOFQ7jH65T7P32XbKWvUzCsbXEahQLoroR3mYUV3ToQHu7pNaYgGWFwxQp1/5t7Pj8VSpteIdY90G2aDzvxd3JRV1v5dr6dv7CmGBghcM4z3WKvakfcOyb/1Lz8ApqKSwLa0lGo4dof3U/bixvc3obE0yscBjHZKWv59fPXqHK1lnEug+TrTHMqTiYCu2G0b5VC5ue1ZggZYXDFCo9cZC0r9/F/f00ah37kdoazrKI1hxuMpDkLjdyfcUy/o5ojLlAVjjMhXOdYs+quRz+dhq19y+hBtls1ap8GPcnqnUcxmWNL7E5MIwJIVY4TMGocnDjV+z+6k2q71pAFT1KuJZjUdkehLfoT3L7rvQqXcLfKY0xDrDCYfLP7ebAL8vYvXwGsb8uIDZnD1FaguUlL+Vkwxto2ak33SpF+zulMcZhVjjMublzyFi7mL3fvkfVXZ8S495PWQ1nVURzVtS7jYsvH0DnhHh/pzTGFCErHOYPNPsEu1Z/yqFVs6i2exFxeohyGsnKEq04Wu8uLmrflzY1E+yeC2OKKSscBoDjGdvYvnw2smkhtY+kUp1TVNCSfB+VzIl619CgfR/aV4vzd0xjTACwwlFMqSuLXT8tYd/3HxGzczEJru00BHZoHMvK90AaXMUll3anQ0wlf0c1xgQYKxzFRY6LPT9/x54fPqHEjq+pdewHqpNFrIbzU0QT1tfoRUyLa2jcrBU1Iu3XwhhzdvYXIlS5c8jYtIr0Hz4l4tevqHV0NVU4ThVgMwl8W74HclFH6ib3JDG+ir/TGmOCiBWOEHHycAa//riEzE3LKJ2xihon1hPHSeKAX6nKqnKd0VqXUb3lVdStcxF17cS2MaaArHAEoVPHj7BrYyqHtqTCzpXEHvqR6u5dNABcGsaWsNqsqtANd0Jr4ptdSb16F1PT7tw2xhQSKxyBTJXDGTvYufFbjm1fTeTetcRmbiQ+J53a4pnyd7+WY2upRmyKvZ4yddtRq1l7GlSqRAM/RzfGhC4rHAEgO+s4u7eu5eD2dWTt2Uj4gU1EH9tG1ewdlOc45b3t0qjCrlL12RxzDSWqN6dKg9bUqFWPpAib9MgYU3QcLRwi0g14AQgHJqvqk2dsF+/2HsBxYKiqrspP32By8vhR9u/czKH0rZzYt42cgzuIOLqTUifSqXgqnSruvdQQpYa3/W5iyChRgzUVryKncgOiayWScElrEipXJsGvX4kxxjhYOEQkHJgIdAXSgBUiMkdV1+Vq1h2o7/1oA7wMtMln3yKlbjcnj2dy/NhhTmYe5uSxo5w6cZjs40fJPnYI19EM9Nhewo7vJzLrAFGnDlI25xDl3IcpTybVgere58pRYa/EcCAijrQyTdhS4SIi4hpQLqER8Rc1pmqFSlT11xdqjDHn4eQeRzKwSVW3AIhICtALyP3HvxfwpqoqsFxEKohIPFA7H30LzaZHEynpPkEYOYSpm3ByCMNNOG7CyCFCcyhFFqVEKXWO58lR4ZCU40hYeY5FVGRP6frsLFkJd3Q8ERVrUjq2FhXjL6Jy9dpULVHSioMxJig5WTiqAztyLafh2as4X5vqZ1l/Zl8ARGQUMAqgZs2aBQp6qMxFiLpQCUclApUwNCwCJAzCIjyPI8tAiTJIyWjCosoSHhVNZKloIktHE1W2IuVi4ilfKY6Y8HBiCpTCGGOCg5OFI6/rPzWfbfLT17NSdRIwCSApKSnPNueTNG5mQboZY0yx5GThSIPT53sBEoBd+WxTIh99jTHG+EGYg8+9AqgvInVEpAQwAJhzRps5wC3i0RY4rKrp+exrjDHGDxzb41BVl4iMARbiuaR2iqquFZHR3u2vAPPxXIq7Cc/luMPO1deprMYYY/JPPBc0hYakpCRNTU31dwxjjAkaIrJSVZN86ePkoSpjjDEhyAqHMcYYn1jhMMYY4xMrHMYYY3wSUifHRWQvsL2A3SsD+woxTlGwzEUj2DIHW16wzEUlr8y1VDXWlycJqcJxIUQk1dcrC/zNMheNYMscbHnBMheVwspsh6qMMcb4xAqHMcYYn1jh+J9J/g5QAJa5aARb5mDLC5a5qBRKZjvHYYwxxie2x2GMMcYnVjiMMcb4JOQLh4h0E5GNIrJJRCbksV1E5N/e7T+KSGJ++wZaZhGpISKLRWS9iKwVkbsDPXOu7eEi8r2IfBQMmb3THM8UkQ3e7/elQZD5Hu/vxU8i8o6IRAVI5ktE5BsRyRKR8b70DbTMAf4ePOv32bs9/+9BVQ3ZDzxDsm8GLsIzOdQPQKMz2vQAPsYz62Bb4Nv89g3AzPFAovdxNPBzoGfOtX0cMB34KNB/N7zbpgIjvY9LABUCOTOe6Zi3AqW8yzOAoQGSOQ5oDfwDGO9L3wDMHMjvwTwz59qe7/dgqO9xJAObVHWLqp4CUoBeZ7TpBbypHsuBCiISn8++AZVZVdNVdRWAqh4F1uP5gxGwmQFEJAHoCUwugqwXnFlEygEdgf8CqOopVT0UyJm92yKAUiISAZSmaGbVPG9mVc1Q1RVAtq99Ay1zIL8Hz/F99vk9GOqFozqwI9dyGn/8IZ6tTX76OuFCMp8mIrWBlsC3hR/xDy408/PAfYDboXx5uZDMFwF7gde9u/aTRaSMk2HPk+e8bVR1J/BP4FcgHc9sm584mPWceYqg74UolNcNwPfguTyPD+/BUC8ckse6M68/Plub/PR1woVk9mwUKQu8D4xV1SOFmO1sCpxZRK4BMlR1ZeHHOqcL+T5HAInAy6raEjgGFMXx9wv5PlfE8x9oHaAaUEZEbi7kfHm5kPdRIL8Hz/0EgfkezLtjAd6DoV440oAauZYT+OPu+dna5KevEy4kMyISiecXdpqqfuBgznzlyUeb9sB1IrINz+71FSLytnNRz5snP23SgDRV/e0/yZl4ConTLiTzlcBWVd2rqtnAB0A7B7OeL4/TfS/EBb1uAL8Hz8b396DTJ238+YHnP8MteP7L+u2EUeMz2vTk9ycTv8tv3wDMLMCbwPPB8n0+o83lFN3J8QvKDHwFXOx9/BDwTCBnBtoAa/Gc2xA8J/fvDITMudo+xO9PNAfse/AcmQP2PXi2zGdsy9d7sMi+MH994LnK5Gc8Vxzc7103Ghid6wc90bt9DZB0rr6BnBnogGf39EdgtfejRyBnLsgvbSBkBloAqd7v9WygYhBkfhjYAPwEvAWUDJDMVfH8x3wEOOR9XO5sfQM5c4C/B8/6fc71HPl6D9qQI8YYY3wS6uc4jDHGFDIrHMYYY3xihcMYY4xPrHAYY4zxiRUOY4wxPrHCYUwevKPf3p5ruZqIzHTota4XkQfP0+afInKFE69vjK/sclxj8uAdZ+gjVW1SBK+1DLhOVfedo00t4DVVvcrpPMacj+1xGJO3J4G6IrJaRJ4Rkdoi8hOAiAwVkdkiMldEtorIGBEZ5x3wcLmIVPK2qysiC0RkpYh8JSKXnPkiItIAyFLVfSIS7X2+SO+2ciKyTUQiVXU7ECMiVYvwe2BMnqxwGJO3CcBmVW2hqn/OY3sTYBCe4az/ARxXz4CH3wC3eNtMwjOsRytgPPBSHs/THsg9DPcXeIYNARgAvK+esaXwtmt/gV+XMRcswt8BjAlSi71/6I+KyGFgrnf9GqCZd3TUdsB7IqcHLi2Zx/PE4xmi/TeT8QxvPRsYBtyaa1sGnpFtjfErKxzGFExWrsfuXMtuPO+rMOCQqrY4z/OcAMr/tqCqX3sPi3UCwlX1p1xto7ztjfErO1RlTN6O4pn6s0DUMwfDVhG5EU7PBd48j6brgXpnrHsTeAd4/Yz1DfAMUGiMX1nhMCYPqrof+FpEfhKRZwr4NDcBI0TkBzxDmuc17ekSoKXkOp4FTAMq4ikewOk5HurhGZHXGL+yy3GN8TMReQGYq6qfeZdvAHqp6uBcbXoDiar6dz/FNOY0O8dhjP89jmeiJUTkP0B3PHMr5BYBPFvEuYzJk+1xGGOM8Ymd4zDGGOMTKxzGGGN8YoXDGGOMT6xwGGOM8YkVDmOMMT75fxGpsOxzhiZMAAAAAElFTkSuQmCC\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY4AAAEGCAYAAABy53LJAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAvHElEQVR4nO3dd3hUZdrH8e9NAqQQeicEIh2kCCFRwAKIIgoISBdBQETFtSwq+9p1VSzr2kUsLCiCiqKoCK5YQLGQWOgoUkxogVBDSJu53z9mZEcMmIFMzszk/lxXLnLOeZ7JLyFn7jynPEdUFWOMMaa4yjkdwBhjTGixwmGMMcYvVjiMMcb4xQqHMcYYv1jhMMYY45dIpwOUpJo1a2rjxo2djmGMMSEjLS1tj6rW8qdPWBWOxo0bk5qa6nQMY4wJGSKy1d8+dqjKGGOMXwJaOESkt4hsEJGNIjKliO39RWSliPwoIqki0s1n2xYRWfX7tkDmNMYYU3wBO1QlIhHAs0AvIANYISILVHWtT7MlwAJVVRFpB7wJtPTZ3l1V9wQqozHGGP8F8hxHMrBRVTcBiMhcoD9wtHCoarZP+1igxOc/KSgoICMjg9zc3JJ+aUdFRUURHx9P+fLlnY5ijCljAlk4GgDpPssZQMqxjURkAPAQUBu42GeTAh+LiAIvqOr0kwmRkZFBXFwcjRs3RkRO5iWCjqqSlZVFRkYGiYmJTscxxpQxgTzHUdS79J9GFKo6X1VbApcC9/ts6qqqHYGLgOtE5Jwiv4jIBO/5kdTdu3f/aXtubi41atQIm6IBICLUqFEj7EZRxpjQEMjCkQE09FmOB7Yfr7GqLgWaiEhN7/J277+ZwHw8h76K6jddVZNUNalWraIvRQ6novG7cPyejDGhIZCFYwXQTEQSRaQCMAxY4NtARJqK9x1QRDoCFYAsEYkVkTjv+ljgAmB1ALMaY0xocRXCov+DAxml/qUDVjhUtRCYBCwG1gFvquoaEZkoIhO9zQYBq0XkRzxXYA1VzwNC6gBfishPwHfAh6q6KFBZT0WXLl2KXD9mzBjmzZtXymmMMWXG5w/BN8/C1uWl/qUDeue4qi4EFh6zbprP5w8DDxfRbxPQPpDZSsry5aX/n2aMKeM2fgLL/gVnjIJ2Q0r9y4fVlCNOqFSpEtnZ2agq119/PZ9++imJiYnYkxWNMQFxcDu8MwFqt4KLHsHlVnILXMRWLL23c5typITMnz+fDRs2sGrVKl588UUbiRhjSp6rEOaNg4JcGDwTKsTw0MJ1DHp+OYfzCksthhWOErJ06VKGDx9OREQE9evXp0ePHk5HMsaEm0/uht+WQ98noFZz3ljxGy99uZmUxOo24ghVdomsMSZgVr8NXz8Dna+CdkP4dlMWd7y7mrOb1eTOS1qXahQrHCXknHPOYe7cubhcLnbs2MFnn33mdCRjTLjIXAfvXQ8NU+DCB0nfm8M1s7+nYbUYnhnekciI0n0rt5PjJWTAgAF8+umntG3blubNm3Puuec6HckYEw5yD8DckVAhFgbP5FChMG7mCgpdbl4anUSVmNKfr84KxynKzvbM0ygiPPPMMw6nMcaEFbcb5l8D+7fC6PdxVarLjbNS+XX3YWZemcxptSo5EssOVRljTLD68nHY8CFc8AA06sJDC9exZH0m9/RtTbdmNR2LZYXDGGOC0cZP4NN/QtvBkHI1s7/dyktfbmZMl8aMOquxo9GscBhjTLDZtxXeHg+1W0PfJ1m2cQ93vbeG7i1qccfFrZxOZ4XDGGOCSl42zB0B6oahr7Jxv5trZ39Ps9qVeHpE6V9BVRTnExhjjPFwu+HdayBzLVw2g6yK8Vz5nxVUjIzg5TGdqVSKN/mdiBUOY4wJFksfhXULoNf95DY6jwmvppF5MI+XRifRoGq00+mOssIRYOnp6XTv3p1WrVrRpk0bnnzySacjGWOC0br34fMHof1w9Mxrue3tlaRt3ce/h3agQ8OqTqf7g+AY94SxyMhI/vWvf9GxY0cOHTpEp06d6NWrF61bl+4UAcaYILZrDbxzNTRIgkue4KlPf+W9H7dzy4Ut6NO2ntPp/sRGHAFWr149OnbsCEBcXBytWrVi27ZtDqcyxgSNw1kwZzhUjIOhr/Hemiz+/cnPDOoYz7XnNXE6XZHK1Ijj3vfXsHb7wRJ9zdb1K3N33zbFartlyxZ++OEHUlJSSjSDMSZEuQrgrdFwaCdc+RFp+ypyy7xvSU6szkMD2wbtxKk24igl2dnZDBo0iCeeeILKlSs7HccYEwwW/x9sWQb9niI9phUTZqVRv0oUL1zeiQqRwfv2XKZGHMUdGZS0goICBg0axMiRIxk4cKAjGYwxQSZtJnw3Hc6axMEWgxj73HIK3crLYzpTLbaC0+lOKHhLWphQVcaNG0erVq24+eabnY5jjAkGW7+GD/8OTXpQ0OMerpv9PZv3HOb5yzvSxKGJC/0R0MIhIr1FZIOIbBSRKUVs7y8iK0XkRxFJFZFuxe0bKr766iteffVVPv30Uzp06ECHDh1YuHCh07GMMU45kAFvjoKqCeigV7jjvXUs+2UPDw5oS5cmzk1c6I+AHaoSkQjgWaAXkAGsEJEFqrrWp9kSYIGqqoi0A94EWhazb0jo1q0bqup0DGNMMMjLhjnDPM8MH/Mhz3y9hzdS0/lbj6YM6dzQ6XTFFsgRRzKwUVU3qWo+MBfo79tAVbP1f++qsYAWt68xxoQUtwveucpzz8bg/zA/I5Z//fdnBp7RgJt6NXc6nV8CWTgaAOk+yxnedX8gIgNEZD3wITDWn77e/hO8h7lSd+/eXSLBjTGmxP33LtiwEHo/zHLpwK3zVtKlSQ2mDmoXtJfdHk8gC0dRP4k/HbNR1fmq2hK4FLjfn77e/tNVNUlVk2rVqnWyWY0xJnDS/gNfPwPJE/i58XCufi2NxJqxPB/kl90eTyATZwC+B+3ige3Ha6yqS4EmIlLT377GGBO0Nn3uuYKq6fns6nI3Y175jujyEcy4Mpkq0aX/vPCSEMjCsQJoJiKJIlIBGAYs8G0gIk3FO0YTkY5ABSCrOH2NMSbo7fkF3rwCajTjcN8XGTvrB/YfKeCVMZ2DarZbfwXsqipVLRSRScBiIAJ4RVXXiMhE7/ZpwCDgChEpAI4AQ70ny4vsG6isxhhT4nL2wutDoFx5CofO4bp3NrJ+5yFeGp3E6Q2qOJ3ulAT04JqqLlTV5qraRFUf8K6b5i0aqOrDqtpGVTuo6lmq+uWJ+oaqsWPHUrt2bU4//fSj6/bu3UuvXr1o1qwZvXr1Yt++fUe3PfTQQzRt2pQWLVqwePFiJyIbY05FYT68MQoObEOHzebOpdl8vmE39/c/ne4tajud7pSF3lmZEDRmzBgWLVr0h3VTp06lZ8+e/PLLL/Ts2ZOpU6cCsHbtWubOncuaNWtYtGgR1157LS6Xy4nYxpiToQof3ARbv4T+z/L8pprM+e43rj2vCSNSEpxOVyKscJSCc845h+rVq/9h3Xvvvcfo0aMBGD16NO++++7R9cOGDaNixYokJibStGlTvvvuu9KObIw5WV89AT++BudO4T13Fx5ZtIH+Heoz+YIWTicrMWVqkkM+mgI7V5Xsa9ZtCxdN9bvbrl27qFfP84CWevXqkZmZCcC2bds488wzj7aLj4+353cYEyrWvQ+f3AOnD+KbhKu45ZUVpCRW55HL2lGuXGjdq3EiNuIIMkVNTxJqNwcZUyZlpMLbV0F8Zzac+RBXvZpGQo0Ypo9KomJkhNPpSlTZGnGcxMggUOrUqcOOHTuoV68eO3bsoHZtzwmz+Ph40tP/d9N8RkYG9evXdyqmMaY49m6C14dCXB129HmFK/6zipgKEcwcm0yVmNC8V+NEbMThkH79+jFz5kwAZs6cSf/+/Y+unzt3Lnl5eWzevJlffvmF5ORkJ6MaY04kZy/MHgzq4sDAuVw+ZxM5+S5mjk0O6Xs1TqRsjTgcMnz4cD7//HP27NlDfHw89957L1OmTGHIkCG8/PLLJCQk8NZbbwHQpk0bhgwZQuvWrYmMjOTZZ58lIiK8hrnGhI2CXM/zwvenkzdiPmPe30v6viO8OjaZlnXD90mfEk5TficlJWlqauof1q1bt45WrVo5lCiwwvl7Myboud0w70pY+y6uQTOYkNaQzzZk8tzIjvQ+vZ7T6YpNRNJUNcmfPnaoyhhjTsYnd8Pad9Fe9/OPDU1Ysj6T+/qfHlJF42RZ4TDGGH999yIsfwo6X8Xj2RfwZmoGf+vRlMvPbOR0slJRJgpHOB2O+104fk/GhIQNH8FHt0Lzi3it2jU8/dmvDE9uGHIPYzoVYV84oqKiyMrKCqs3WlUlKyuLqKgop6MYU7Zs+x7mjYV67Vnc6gHufH8957eqw/39Ty9T91uF/VVV8fHxZGRkEG5PB4yKiiI+Pt7pGMaUHfu2eGa7ja1JapdpXD/nZzomVOPp4WcQGRH2f4P/QdgXjvLly5OYmOh0DGNMKPv9Xg1XAb9e/AZXvrGVRjVieHl0EtEVyt7l8mFfOIwx5pTk53hGGvu2knnpHIbP30dsxUhmjk2makwFp9M5omyNr4wxxh+uQs+9GtvSOHDx8wxdHEFugYtZ45KpH6Z3hReHFQ5jjCmKKnxwA/y8iCMXPMzwZbXZeSCXGVd2pnmdOKfTOcoKhzHGFOXT++GH1yjodgujf2rLL5mHmDaqE50aVf/rvmHOCocxxhzr2xdg2b9wdRzNhN96sWLrXv49tAPnNq/ldLKgYIXDGGN8rX4HProNbdGHyYdH8dnPe3jg0rZc0s4eb/C7gBYOEektIhtEZKOITCli+0gRWen9WC4i7X22bRGRVSLyo4ikHtvXGGNK3OalMP9qNOFM7o+azPyfMrmtd8uweVZ4SQnY5bgiEgE8C/QCMoAVIrJAVdf6NNsMnKuq+0TkImA6kOKzvbuq7glURmOMOWrHSpgzAqo34bm69/PK0p1cfc5pXHNeE6eTBZ1AjjiSgY2quklV84G5QH/fBqq6XFX3eRe/AexWaGNM6du3BWZfBlFVmNPiCR5dmsnQpIZMuail08mCUiALRwMg3Wc5w7vueMYBH/ksK/CxiKSJyITjdRKRCSKSKiKp4TatiDGmFGTvhlcHQmEeizs+xz8+yaJP27o8OLBtmZp/yh+BvHO8qJ94kTMNikh3PIWjm8/qrqq6XURqA/8VkfWquvRPL6g6Hc8hLpKSksJnJkNjTODlHoDXBsLB7Xx79itcu/gwZzeryb+HdiCinBWN4wnkiCMDaOizHA9sP7aRiLQDXgL6q2rW7+tVdbv330xgPp5DX8YYUzLyc+D1YZC5llXdnmHUx0LbBlWYdnknKkaWvfmn/BHIwrECaCYiiSJSARgGLPBtICIJwDvAKFX92Wd9rIjE/f45cAGwOoBZjTFlSWE+vDUafvuaX7o+zuAlsTStXYmZVyYTW9Gm8PsrAfsJqWqhiEwCFgMRwCuqukZEJnq3TwPuAmoAz3mPJRZ6n31bB5jvXRcJvK6qiwKV1RhThrhd8O5E+OVjfus6lUu/qEN8tWheHZdMlZjyTqcLCRJODzhKSkrS1FS75cMYcxyq8MFNkDaDXcn/R6/vOlA1pgJvXn0WdauUzQejiUia9w/2YrMxmTGm7FhyH6TNYN8Z19InrSOVKpZj9viUMls0TpZNOWKMKRu+ehK+fJxDp4+i9+oeiAivjU+hYfUYp5OFHCscxpjwl/Yf+O9dHGnenz6/9Cffpcwen8JptSo5nSwkWeEwxoS31e/A+zeSl9iTfttGsT/XzayxKbSoW7afqXEq7ByHMSZ8/fIJvDOBggbJDM6aSMZBF6+OS6ZtfBWnk4U0KxzGmPC06Qt4YySuWi0ZlXMz6/e6mDGmM0mN7UFMp8oOVRljws9v38CcYbiqNmZ0wT9Iy3Qz7fKOdG1a0+lkYcEKhzEmvGxLg9cuwx1Xj/F6B9/uEp4b2YkeLes4nSxsWOEwxoSPnavg1YG4Y6pzdbm7WbYjgmdGdKRXaysaJcnOcRhjwkPmepjVH3f5GCZF3sun28vzzPAzuLBNXaeThR0bcRhjQl/WrzCrHyqR3Bh1P4u3V+TJYR24qG09p5OFJRtxGGNC276tMLMf6irklkoP8UFGNP8e2p5L2tV3OlnYssJhjAldB7d7Rhr5h7i98lTeTq/E40Pa07/DiR42ak6VFQ5jTGjKzvSMNA7v4e6qDzEnvQqPXdaeAWfEO50s7Nk5DmNM6Mne7SkaB7fxYLX7mfVbDR4e1I5BnaxolAYrHMaY0JK9G2b2RfdtYWq1e3lxax0eGtiWIUkN/7qvKRF2qMoYEzp8isaDVe/lpfQGPDKoHUM6W9EoTVY4jDGhwado3F/lHv6zLZ7HLmtvh6ccYIXDGBP8fIrGvXF3M2t7Q/49tINdPeUQKxzGmOCWvdtzye2+Ldxd6S5m72rE08PP4OJ2dnOfUwJ6clxEeovIBhHZKCJTitg+UkRWej+Wi0j74vY1xpQBvxeNvZu5M/ZOXs9szLMjrGg4LWCFQ0QigGeBi4DWwHARaX1Ms83AuaraDrgfmO5HX2NMOPMpGrfH3MEbexJ5/vJO9D7diobTAjniSAY2quomVc0H5gL9fRuo6nJV3edd/AaIL25fY0wY8yka/4i6g3lZTXhhVCeb5TZIBLJwNADSfZYzvOuOZxzwkb99RWSCiKSKSOru3btPIa4xJihkZx4tGlMq3s78/U14cXSSPU8jiASycEgR67TIhiLd8RSO2/ztq6rTVTVJVZNq1ap1UkGNMUHi4HaY0Qf33s1MLv9/vHewKa+M6cy5zW3fDiaBvKoqA/C9Kyce2H5sIxFpB7wEXKSqWf70NcaEkX1bYVY/3Nl7uK7cHXyZ04xXx3Wmsz0jPOgEcsSxAmgmIokiUgEYBizwbSAiCcA7wChV/dmfvsaYMJL1K8zog+vwPsa67+Cbwua8ftWZVjSCVMBGHKpaKCKTgMVABPCKqq4RkYne7dOAu4AawHMiAlDoPexUZN9AZTXGOMj75L7CwnxGFtzO5sgmvDE+heZ14pxOZo5DVIs8dRCSkpKSNDU11ekYxpji2rkKZl1KvgpDcqawO/o0Xr8qhUY1Yp1OVmaISJqqJvnTx+4cN8Y4Y9v38OoA8spF0S/7Ngqqnsa88SnUqxLtdDLzF6xwGGNK32/fwuzLyImszCX7b6Vi7dN4c1wyNStVdDqZKQZ7HocxpnRtXgavDuBQZHV67ZtC5fpNmXvVmVY0QogVDmNM6fn5Y5h9Gfsr1qXH3ltJaNyM18anUCWmvNPJjB+scBhjSseqeejc4eyKSqTHnlto16I5M67sTKWKdsQ81FjhMMYE3ncvom+PZ0t0W3ru+TvdO7Zm2qhORJWPcDqZOQlW6o0xgaMKSx+Dz/7JqtguDM6awJhzWjLlopZ4790yIahYIw4RGXfMcoSI3B2YSMaYsOB2w+L/g8/+ybKY8xmQdQ03X9SOf/RpZUUjxBX3UFVPEVkoIvVE5HQ8U6DbbZ3GmKK5CmHBJPjmORZE9ePK/Vcy9bIzuPrcJk4nMyWgWIeqVHWEiAwFVgE5wHBV/SqgyYwxoakgF94eB+s/YEaF4UzN7se0yztxvj1LI2wUq3CISDPgBuBtoBUwSkR+UNWcQIYzxoSYvEMwZzhsWcZjEeOYlX8hr423GW7DTXFPjr8PXKeqS8RzcPImPDPYtglYMmNMaDmcBbMHoTtWcjvX80m583jzqmRa1q3sdDJTwopbOJKB8SJyHZ4HKn2JZ6pzY4yBvZvhtUG49mdwXeHNbKjSjbfHJtOweozTyUwAFLdwvAQcAp72Lg8HzgKGBCKUMSaEbP8RZg8mLy+XkblTKGyQwrzRSdSwKUTCVnELRwtVbe+z/JmI/BSIQMaYELJxCfrmFRykEgMP38lprTry1LAziK5gN/aFs+JejvuDiJz5+4KIpAB2VZUxZdmPc9DXh7CNOvQ6eCfdzuzCtMs7WdEoA4o74kgBrhCR37zLCcA6EVkFqKq2C0g6Y0zwUYUvH4cl97G6QntGHLyev/XpxPizE+3GvjKiuIWjd0BTGGNCg9sFH90KK15iSeQ5/O3wBKYO70Tf9vWdTmZKUXFvANwa6CDGmCBXcATeHg/rP2CW9OPxwpHMGJ9McqLdo1HWBHR2XBHpLSIbRGSjiEwpYntLEflaRPJEZPIx27aIyCoR+VFE7EHixjgpZy/MuhRd/yEPua/ghYpXMu/arlY0yqiAzY4rIhHAs0AvIANYISILVHWtT7O9wN+AS4/zMt1VdU+gMhpjimHvJpg9GNferdxQcD1b6l7A/DGdqR0X5XQy45BAjjiSgY2quklV84G5QH/fBqqaqaorgIIA5jDGnKytX6Mv9iRnfybDc6eQ27w/b0w4y4pGGRfIwtEASPdZzvCuKy4FPhaRNBGZUKLJjDF/beVb6Kx+7CyIpk/O3XQ4+2JeGNWJWHtiX5kXyN+Aoq7LUz/6d1XV7SJSG/iviKxX1aV/+iKeojIBICEh4eSSGmP+RxW+eBg+f4hVEW0Ym3MDtww8i6Gdbf8yHoEccWQADX2W44Htxe2sqtu9/2YC8/Ec+iqq3XRVTVLVpFq1ap1CXGMMhXkw/2r4/CHel3MZ67qdp8b1tKJh/iCQI44VQDMRSQS24ZkUcURxOopILFBOVQ95P78AuC9gSY0xniun5o6E35bzuGso71cezltXJpNYM9bpZCbIBKxwqGqhiEwCFgMRwCuqukZEJnq3TxORukAqUBlwi8iNQGugJjDfexdqJPC6qi4KVFZjyrw9G9HXB+Pal8FN+ZPY07gv8y/vSNWYCk4nM0EooGe5VHUhsPCYddN8Pt+J5xDWsQ4C7YtYb4wpaVu+ROeOJLtAGZ37fzRP6snjl55O+YiA3uZlQphdHmFMWZY6A104mXTqcvmRv3NFn/MY183mnDInZoXDmLLIVQCL/gErXmQ5HbiFG3hgTDe6t6jtdDITAqxwGFPWHM5C3xqNbFnGi65LeKvqOF67IpnTalVyOpkJEVY4jClLdq3FPWcY7gPbuTV/IgdbDObtoe2JiyrvdDITQqxwGFNWrP8Q99tXsd9VgXG5d3B2jz7c2LMZ5crZ+QzjHyscxoQ7VVj2GHz6T9ZLE651TWbKyO70Pr2e08lMiLLCYUw4yzsE710Ha99jgbsrT8X+jRdGd6VF3Tink5kQZoXDmHC15xfcc0fAno08XDCctYljmDfCbuozp84KhzHhaN0HuOdfzaGCclyTP4W2Z/djxgUtiLSb+kwJsMJhTDhxu+CzB2HZY6yjCTfq35k8sicXtqnrdDITRqxwGBMucvaib49Hfl3CG67zeK369UwfdZZNUmhKnBUOY8LBjpW45l6O+8B27iwYT367Ubw5oC3RFSKcTmbCkBUOY0LdD6/h/uBmslyVuLbwLgb068+I5ASbb8oEjBUOY0JV/mH0w78jP83hW3cb/hk9mQfH9aB9w6pOJzNhzgqHMaFo9wZcb1yB7NnAk4UD+em0Cbw6tBPVY+1SWxN4VjiMCTUr38S14AYOFkZyY8EUul44mJe7nWZTh5hSY4XDmFBRcAT9aAry/X9Ic7fkgehbuHtsTzomVHM6mSljrHAYEwqyfqVw7igid6/hucJ+rGw2iVmDO1Ilxma1NaXPCocxwe6nN3C9fxM5hcLkwlvpctEInu/S2K6aMo6xwmFMsMo9iPvDv1Nu1Zt8727BI7GTuWvkhbSNr+J0MlPGBXTiGhHpLSIbRGSjiEwpYntLEflaRPJEZLI/fY0Ja9vSKHj+bHTVPP5dMIg32jzHKzcMsKJhgkLARhwiEgE8C/QCMoAVIrJAVdf6NNsL/A249CT6GhN+3G50+VPokvvZ7a7KFLmHy4YM5qb29Z1OZsxRgTxUlQxsVNVNACIyF+gPHH3zV9VMIFNELva3rzFh59BOCuZNoPzWL1joSmZevVt4aMTZNKga7XQyY/4gkIWjAZDus5wBpJR0XxGZAEwASEhI8D+lMcFgwyIK3rkGV142dxdeRYMeV/PSeU2JsHszTBAKZOEo6jdeS7qvqk4HpgMkJSUV9/WNCQ55h3B99A8ifnyVX9yNmBp7LzeP6EcHmzbEBLFAFo4MoKHPcjywvRT6GhMatn5N/rwJRB5K5/nCvqS3u5Hn+3cgtqJd7GiCWyB/Q1cAzUQkEdgGDANGlEJfY4JbYR7uTx9Alj/FLq3FvZH3M2LYEK5pWcfpZMYUS8AKh6oWisgkYDEQAbyiqmtEZKJ3+zQRqQukApUBt4jcCLRW1YNF9Q1UVmNKza415L05jopZ65hT2J3vW93Co5d2pppNTmhCiKiGz2mBpKQkTU1NdTqGMX/mduFe/iy65D72uWO4v9xELhhwJRe3q+d0MlPGiUiaqib508cOphoTaLt/JvftiUTtTGORqzOLEm/j9sFnUzsuyulkxpwUKxzGBIqrEPfyp3F/9iC5rvLcyyQ69ruafyc1tHmmTEizwmFMIOxay5F5E4ne/ROLXZ35b+NbufWyc6hbxUYZJvRZ4TCmJLkKKFz6OLL0EXLcMdwTcTNnDxjPv9rVt1GGCRtWOIwpKTtWkvPW1cTsXcsC11l813IKUy7tYldMmbBjhcOYU5V/mPwlDxH57XNkaxz3lb+NC4eN558tazudzJiAsMJhzCnQnxeT++6NROdsZ46rO5va38LtlyQTF2VP5jPhywqHMSfj0E5yFtxCzC8LSHc3YHqVhxkxeCjD7fnfpgywwmGMP9xuCle8jOvju4kozOdJHUpcz78ztVszIiMC+lw0Y4KGFQ5jimvnag7Nm0Tcnh/4xtWGjxvfxsRBF1DfnpdhyhgrHMb8lSP7yFl8P1E/ziBfY7m3wg10G3It97Wu63QyYxxhhcOY43G7KUibSeHH91Ax/wCvu3uxN2Uyt16QRHSFCKfTGeMYKxzGFCUjjYPv3EDlvav4wd2ChfEPM2ZgXxrXjHU6mTGOs8JhjK/s3Rz88E4qr5vDEa3K09E30fXSa7jHnpVhzFFWOIwBKMwjd/k0+OIRoguPMIO+lDvvNm45uw0VIu1qKWN8WeEwZZsqhavf5cjCO4g7ksEXrnasaDGZ0f17UyuuotPpjAlKVjhMmaXpK9j/7q1Uy/qebe6GzK/5AJcMHMXk+CpORzMmqFnhMGXP/t/Yu+B2qm9aQKFW4fHo6zij33VMaWUz2BpTHFY4TNlxZB8H/vsoMT9MJ8YNL0dcRuXzJ/O3lBZ217cxfrDCYcJfXjYHlz5D+a+fJs51mPfpxu7kWxl+fhdiK9ouYIy/ArrXiEhv4EkgAnhJVaces1282/sAOcAYVf3eu20LcAhwAYX+PkzdGArzyF7+Eix9lMqF+1ji7sj61jcwuE9vale2J/EZc7ICVjhEJAJ4FugFZAArRGSBqq71aXYR0Mz7kQI87/33d91VdU+gMpow5XaRk/o6BUseoEreDr5xtyK1yT/p33cgPavHOJ3OmJAXyBFHMrBRVTcBiMhcoD/gWzj6A7NUVYFvRKSqiNRT1R0BzGXCldvNkVXvcmTRfVQ/spmf3KexLOExLuo3jEm145xOZ0zYCGThaACk+yxn8MfRxPHaNAB2AAp8LCIKvKCq04v6IiIyAZgAkJCQUDLJTWhxu8j54S2OLHmYGjmb2Oauz9y6d9P90nFMqm+X1hpT0gJZOIq6rlH9aNNVVbeLSG3gvyKyXlWX/qmxp6BMB0hKSjr29U04cxVyKHUO+Z89So3craS745lX53ZSLhnHtY1qOJ3OmLAVyMKRATT0WY4Hthe3jar+/m+miMzHc+jrT4XDlEGuAg5++xqFXzxK9bxtrHMn8E79e+jSdwxXN7An8BkTaIEsHCuAZiKSCGwDhgEjjmmzAJjkPf+RAhxQ1R0iEguUU9VD3s8vAO4LYFYTCvJz2PfVDPj6aarl72C1O5H5CQ9wbt9RXFXHDkkZU1oCVjhUtVBEJgGL8VyO+4qqrhGRid7t04CFeC7F3Yjnctwrvd3rAPO9d/FGAq+r6qJAZTVB7vAedn7yFJV+mkE190F+dDfl3cSH6dH3csbVrOR0OmPKHPFc0BQekpKSNDU11ekYpoS49/zKto8eo86v86hAPp+RREbr8Zx/QX/qVbXLao0pCSKS5u99cnbbrAk6eVu+ZddHjxC/awm1NYJFkeeRn3wdvc87h+52p7cxjrO90ASHwjyyvp1L3vIXqH94DVU0hnkxg6l63iT6JLW1uaSMCSJWOIyj3PvSSf/kWaqtm0MN9342uuszp/b1NOl1FYObJ9hstcYEISscpvSpkr3hc3YveZqE3Z8Rr8qycklkthpNtwsGMbyanb8wJphZ4TClRrN3k/HFDCqsfI06eVsp0EosiB1IbNerOTelExUjI5yOaIwpBiscJrDcLvavXkzWspdJ2P05DSnkR23GsvgptL5wLAMS6jid0BjjJyscJiAKsrbw25IXqbrhDWq4duPWSiyO7Utk59Gc3eVsOtjVUcaELNt7TYnRnH2kL5+L66c3aXToBxIVvivXjmXNbqJdz2FcUtfmjzImHFjhMKem4Ag7VrxLduocGu39igQK2ax1+bD6FVTvOoaUMzrYpbTGhBkrHMZ/rgJ2r/qErG9m03DnEuqRwy6typK4fpTvMITkrj3pG13B6ZTGmACxwmGKRfNz2PH9Rxz44R3iMz+nlmZTUaP5Jrobha0vo+O5/bioil1Ga0xZYIXDHJf7yAG2fjOf3FXv0XjvcuqTS6zGkBp9FrlN+tDmnIH0rFPd6ZjGmFJmhcP8jyrZO9aT/u0CIn79hMbZ35NIIZlala8qnQ+t+tK268V0r2aPYTWmLLPCUca5c7NJ/34xB1d/RO1dy6jj2kkrYJPW54tqA6jY9lLan9WL82MqOh3VGBMkrHCUNa4Cdm34lsxVn1Dht6UkHv6JRhRyWCuyukJ7fmp4BbXPuJg2p7fjNLsayhhTBCscYU4L89m5/hsyV31CxYyvSTi8kjrkUgfYSEOWVRtI+RYX0jKlFynV7Cl6xpi/ZoUjzOQeyCRj9Zcc+vVronb+QKOcVdQjl3p4CsV3VS6ERmfT4IyeNGmUSNNyNvusMcY/VjhCmDs/l20bvmXP+uWU255GrQOrqe/eQVPApcKmcgmkVu0NjbvRsMP5NGnUmKY2Tbkx5hRZ4QgFquzfuZkdP6dxOP0nInevpXr2L9QvzKChuGkI7NLqbI1pxc+1BhFzWgqN2nSlWa0aNHM6uzEm7FjhCCKFeTns3rqBrPR15O36Gfb+SuyhzTTI+5Wq5FDV224btdgR1YTNtXpQvkF76rXpRqPGzahjh52MMaUgoIVDRHoDTwIRwEuqOvWY7eLd3gfIAcao6vfF6RuKcg8fYO+OrRzYtYUjWekU7ttGuUPbiM7eSo28DGq791BPlHre9ns1jh2RDVlZrRfu2m2IS2hPg+YdqV+zFg3skJMxxiEBKxwiEgE8C/QCMoAVIrJAVdf6NLsIaOb9SAGeB1KK2dcRbpeLvCOHyTuS7fk39zB52fvJPbiH/EO7cWXvxZ2TRbkj+4jM20eFggNUyt9DdXcWlTlMfaC+z+vt0zgyI+vxW6X2bKySSGTNpsQ1aEHtRq2pWas21a1AGGOCTCBHHMnARlXdBCAic4H+gO+bf39glqoq8I2IVBWRekDjYvQtMZvva095zUVQPG/Tiuj/Pi+HmwrkE6V5REkB0UD0CV7PrcJBiSVb4jgcUZm9UQ3ZEd0ZjatHRNUGRNdsSJU6jahetzHV4ipTLRDflDHGBEggC0cDIN1nOQPPqOKv2jQoZl8ARGQCMAEgISHhpILuqdwKcRf8/oqoCOD9EE85cUdGo+WjITIaykcj5WMoVzGGchViiIiKI6pyTWKq1iKuWm0qV61J1fLlj56TMMaYcBLIwlHUMRYtZpvi9PWsVJ0OTAdISkoqss1f6Xzj3JPpZowxZVIgC0cG0NBnOR7YXsw2FYrR1xhjjAMCORnRCqCZiCSKSAVgGLDgmDYLgCvE40zggKruKGZfY4wxDgjYiENVC0VkErAYzyW1r6jqGhGZ6N0+DViI51LcjXgux73yRH0DldUYY0zxieeCpvCQlJSkqampTscwxpiQISJpqprkTx+bN9sYY4xfrHAYY4zxixUOY4wxfrHCYYwxxi9hdXJcRHYDW0+ye01gTwnGKQ2WuXSEWuZQywuWubQUlbmRqtby50XCqnCcChFJ9ffKAqdZ5tIRaplDLS9Y5tJSUpntUJUxxhi/WOEwxhjjFysc/zPd6QAnwTKXjlDLHGp5wTKXlhLJbOc4jDHG+MVGHMYYY/xihcMYY4xfwr5wiEhvEdkgIhtFZEoR20VEnvJuXykiHYvbN9gyi0hDEflMRNaJyBoRuSHYM/tsjxCRH0Tkg1DI7H3M8TwRWe/9eZ8VAplv8v5erBaROSISFSSZW4rI1yKSJyKT/ekbbJmDfB887s/Zu734+6Cqhu0HninZfwVOw/NwqJ+A1se06QN8hOepg2cC3xa3bxBmrgd09H4eB/wc7Jl9tt8MvA58EOy/G95tM4Hx3s8rAFWDOTOexzFvBqK9y28CY4Ikc22gM/AAMNmfvkGYOZj3wSIz+2wv9j4Y7iOOZGCjqm5S1XxgLtD/mDb9gVnq8Q1QVUTqFbNvUGVW1R2q+j2Aqh4C1uF5wwjazAAiEg9cDLxUCllPObOIVAbOAV4GUNV8Vd0fzJm92yKBaBGJBGIonadq/mVmVc1U1RVAgb99gy1zMO+DJ/g5+70PhnvhaACk+yxn8Of/xOO1KU7fQDiVzEeJSGPgDODbko/4J6ea+QngVsAdoHxFOZXMpwG7gRneof1LIhIbyLB/kecv26jqNuAx4DdgB56nbX4cwKwnzFMKfU9FiXzdINwHT+QJ/NgHw71wSBHrjr3++HhtitM3EE4ls2ejSCXgbeBGVT1YgtmO56Qzi8glQKaqppV8rBM6lZ9zJNAReF5VzwAOA6Vx/P1Ufs7V8PwFmgjUB2JF5PISzleUU9mPgnkfPPELBOc+WHTHk9gHw71wZAANfZbj+fPw/HhtitM3EE4lMyJSHs8v7GxVfSeAOYuVpxhtugL9RGQLnuF1DxF5LXBR/zJPcdpkABmq+vtfkvPwFJJAO5XM5wObVXW3qhYA7wBdApj1r/IEuu+pOKWvG8T74PH4vw8G+qSNkx94/jLchOevrN9PGLU5ps3F/PFk4nfF7RuEmQWYBTwRKj/nY9qcR+mdHD+lzMAyoIX383uAR4M5M5ACrMFzbkPwnNy/Phgy+7S9hz+eaA7affAEmYN2Hzxe5mO2FWsfLLVvzKkPPFeZ/IznioPbvesmAhN9/qOf9W5fBSSdqG8wZwa64RmergR+9H70CebMJ/NLGwyZgQ5Aqvdn/S5QLQQy3wusB1YDrwIVgyRzXTx/MR8E9ns/r3y8vsGcOcj3weP+nH1eo1j7oE05Yowxxi/hfo7DGGNMCbPCYYwxxi9WOIwxxvjFCocxxhi/WOEwxhjjFyscxhTBO/vttT7L9UVkXoC+1qUictdftHlMRHoE4usb4y+7HNeYInjnGfpAVU8vha+1HOinqntO0KYR8KKqXhDoPMb8FRtxGFO0qUATEflRRB4VkcYishpARMaIyLsi8r6IbBaRSSJys3fCw29EpLq3XRMRWSQiaSKyTERaHvtFRKQ5kKeqe0Qkzvt65b3bKovIFhEpr6pbgRoiUrcUfwbGFMkKhzFFmwL8qqodVPWWIrafDozAM531A0COeiY8/Bq4wttmOp5pPToBk4HninidroDvNNyf45k2BGAY8LZ65pbC267rKX5fxpyySKcDGBOiPvO+0R8SkQPA+971q4B23tlRuwBviRyduLRiEa9TD88U7b97Cc/01u8CVwJX+WzLxDOzrTGOssJhzMnJ8/nc7bPsxrNflQP2q2qHv3idI0CV3xdU9SvvYbFzgQhVXe3TNsrb3hhH2aEqY4p2CM+jP0+Kep7BsFlEBsPRZ4G3L6LpOqDpMetmAXOAGcesb45ngkJjHGWFw5giqGoW8JWIrBaRR0/yZUYC40TkJzxTmhf12NOlwBniczwLmA1Uw1M8gKPPeGiKZ0ZeYxxll+Ma4zAReRJ4X1U/8S5fBvRX1VE+bQYAHVX1TodiGnOUneMwxnkP4nnQEiLyNHARnmcr+IoE/lXKuYwpko04jDHG+MXOcRhjjPGLFQ5jjDF+scJhjDHGL1Y4jDHG+MUKhzHGGL/8P1Op12cPf+MQAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] diff --git a/src/symba/symba_step.f90 b/src/symba/symba_step.f90 index b56db66d5..d34956f4b 100644 --- a/src/symba/symba_step.f90 +++ b/src/symba/symba_step.f90 @@ -52,6 +52,7 @@ module subroutine symba_step_interp_system(self, param, t, dt) real(DP), intent(in) :: dt !! Current stepsize ! Internals real(DP) :: dth !! Half step size + integer(I4B) :: irec !! Recursion level dth = 0.5_DP * dt associate(system => self) @@ -72,11 +73,11 @@ module subroutine symba_step_interp_system(self, param, t, dt) call pl%kick(dth) call tp%kick(dth) - - call pl%drift(system, param, dt, pl%status(:) == ACTIVE) - call tp%drift(system, param, dt, tp%status(:) == ACTIVE) - - call system%recursive_step(param, 0) + irec = -1 + call pl%drift(system, param, dt, pl%status(:) == ACTIVE .and. pl%levelg(:) == irec) + call tp%drift(system, param, dt, tp%status(:) == ACTIVE .and. tp%levelg(:) == irec) + irec = 0 + call system%recursive_step(param, irec) call pl%set_beg_end(xend = pl%xh) call pl%accel(system, param, t + dt)