Skip to content

Instantly share code, notes, and snippets.

@hannorein
Created June 20, 2021 23:29
Show Gist options
  • Save hannorein/9bbcb4f1c1d6c44079061d18d47ab226 to your computer and use it in GitHub Desktop.
Save hannorein/9bbcb4f1c1d6c44079061d18d47ab226 to your computer and use it in GitHub Desktop.
2014 UN271 with REBOUND
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "3a53075c",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Searching NASA Horizons for 'Sun'... Found: Target body name: Sun (10).\n",
"Searching NASA Horizons for 'Mercury'... Found: Target body name: Mercury Barycenter (199).\n",
"Searching NASA Horizons for 'Venus'... Found: Target body name: Venus Barycenter (299).\n",
"Searching NASA Horizons for 'Earth'... Found: Target body name: Earth-Moon Barycenter (3).\n",
"Searching NASA Horizons for 'Mars'... Found: Target body name: Mars Barycenter (4).\n",
"Searching NASA Horizons for 'Jupiter'... Found: Target body name: Jupiter Barycenter (5).\n",
"Searching NASA Horizons for 'Saturn'... Found: Target body name: Saturn Barycenter (6).\n",
"Searching NASA Horizons for 'Uranus'... Found: Target body name: Uranus Barycenter (7).\n",
"Searching NASA Horizons for 'Neptune'... Found: Target body name: Neptune Barycenter (8).\n"
]
}
],
"source": [
"import rebound\n",
"import numpy as np\n",
"sim=rebound.Simulation()\n",
"sim.add([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"],date=\"JD2459400.5\")"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "0101f42c",
"metadata": {},
"outputs": [],
"source": [
"# Add data from https://www.minorplanetcenter.net/mpec/K21/K21M53.html\n",
"sim.add(P=610989*np.pi*2., # period given in years, convert to year/2pi\n",
" e=0.9984815,\n",
" M=359.99436/180.*np.pi, # degrees to radians\n",
" inc=95.44606/180.*np.pi,\n",
" Omega=189.92932/180.*np.pi,\n",
" omega=326.37976/180.*np.pi)\n",
"# Could also use Horizon data, but it is slightly older\n",
"# sim.add(\"2014 UN271\")"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "ddb657ce",
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAfEAAAHkCAYAAAAuKRZVAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Z1A+gAAAACXBIWXMAAAsTAAALEwEAmpwYAACBk0lEQVR4nO3dd5wb9Z3/8dd3Rl3bvS7r3hsY29hgOgQIxsGE9EAuuZBG6iUhPYRckh/h0i7l0o6QcuFCEkKOJJRgTO9gMLjbuPe63r6SVmXm+/tjpLV2WdtrW9rRaD/Ph/VY7UgrfdYr6T3f73y/31Faa4QQQgjhPYbbBQghhBDi5EiICyGEEB4lIS6EEEJ4lIS4EEII4VES4kIIIYRHSYgLIYQQHuVaiCulxiilnlBKrVdKrVNKfSa7vU4p9YhSanP2a61bNQohhBClTLk1T1wp1QA0aK1fVUpVAq8AbwGuB5q11t9VSn0FqNVaf9mVIoUQQogS5lpLXGu9X2v9avZ6B7ABGAVcA9yRvdsdOMEuhBBCiF5ca4n3KEKp8cDTwOnALq11TXa7Alpy3wshhBDiCJ/bBSilKoB7gM9qrdud3HZorbVSqs+9DKXUDcANANFodN706dMHotySk7BtfErhz/t/E8KLXnnllcNa66Fu1yGEl7ga4kopP06A/1Fr/bfs5oNKqQat9f7scfNDff2s1vp24HaA+fPn6+XLlw9IzaXmmdZW5ldWEjZNt0sR4pQopXa6XYMQXuPm6HQF/BbYoLX+Ud5N9wHvz15/P3DvQNfmFV2WhQYJcCGEGKTcbImfD7wPWKOUWpnddhPwXeBupdSHgJ3Au9wpr/S1WRbVPtePiAghhHCJawmgtX4WONqB3MsGshavastkqJZWuBBCDFqyYpuHdUpLXAghBjUJcQ/rtCwqpCUuhBCDloS4RyVtGwPwG/InFEKIwUoSwKM6LYuotMKFEGJQkxD3qJh0pQshxKAnIe5RMWmJCyHEoCch7lEyqE0IIYSEuAdprYnbNhEZ1CaEEIOapIAHddk2fqXwSYgLIcSgJingQXI8XAghBEiIe5IcDxdCCAES4p4Us22i0pUuhBCDniSBB0l3uhBCCJAQ9xytNXEJcSGEEEiIe07CsoiYJoY62llchRBCDBauhrhS6ndKqUNKqbV5276plNqrlFqZvbzJzRpLTZfWBOV4uBBCCNxvif8euLKP7T/WWs/JXh4c4JpKWtK2CUmICyGEwOUQ11o/DTS7WYPXdNk2QelKF0IIgfst8aP5lFJqdba7vdbtYkpJl7TEhRBCZJViGvw3MAmYA+wHftjXnZRSNyilliulljc2Ng5gee6SEBdCCJFTcmmgtT6otba01jbwa+Dso9zvdq31fK31/KFDhw5skS5K2rYMbBNCCAGUYIgrpRryvn0rsPZo9x1stNYkZXS6EEKILJ+bT66U+jNwCVCvlNoDfAO4RCk1B9DADuCjbtVXapK2TUApmSMuhBACcDnEtdbX9bH5twNeiEcktZbj4UIIIbpJInhIlxwPF0IIkUcSwUNkZLoQQoh8kggekpKWuBBCiDySCB6S1lpWaxNCCNFNQtxDMlrjkxAXQgiRJSHuIWmt8Ut3uhBCiCxJBA+RlrgQQoh8EuIekrZtCXEhhBDdJMQ9wtYaDZgS4kIIIbIkxD0iozV+CXAhhBB5JMQ9Qo6HCyGE6E1C3CNkZLoQQojeJBU8QlriQgghepMQ94iMbcsxcSGEED1IiHuEjfyxhBBC9CS54BEaUNISF0IIkcfVEFdK/U4pdUgptTZvW51S6hGl1Obs11o3aywVGtnjEkII0ZPbufB74Mpe274CPKa1ngI8lv1+0NNaI+1wIYQQ+VwNca3100Bzr83XAHdkr98BvGUgaypVNtKdLoQQoie3W+J9Ga613p+9fgAY3tedlFI3KKWWK6WWNzY2Dlx1LtEgLXEhhBA9lGKId9Naa5z86uu227XW87XW84cOHTrAlQ08rXVp/7GEEEIMuFLMhYNKqQaA7NdDLtdTEmR0uhBCiN5KMcTvA96fvf5+4F4XaykZNtKdLoQQoie3p5j9GXgBmKaU2qOU+hDwXeCNSqnNwOXZ7wc9pfs8qiCEEGIQ87n55Frr645y02UDWogHSFe6EEKI3kqxO130wVCq7xF+QgghBi0JcY8wAFu61IUQQuSREPcIhTO4TQghhMiREPcIQym0tMSFEELkkRD3CANpiQshhOhJQtwjlFJyTFwIIUQPEuIeYXCU9WeFEEIMWhLiHmGCHBMXQgjRg4S4R/gMg7TbRQghhCgpEuIe4VeKjLTEhRBC5JEQ9wifUqRtGZ8uhBDiCAlxj/AphYUcFxdCCHGEhLhHKKUwAUtCXAghRJaEuIf4lCItIS6EECJLQtxD/IYhg9uEEEJ0c/V84seilNoBdAAWkNFaz3e3IvcFZXCbEEKIPCUb4llv0FofdruIUhEwDFLSEhdCCJEl3ekeElCKpIS4EEKIrFIOcQ08rJR6RSl1Q+8blVI3KKWWK6WWNzY2ulDewAsaBknpThdCCJFVyiF+gdb6TGAR8Eml1EX5N2qtb9daz9dazx86dKg7FQ4wCXEhhBD5SjbEtdZ7s18PAX8Hzna3IvdJiAshhMhXkiGulIoqpSpz14ErgLXuVuU+UykMGaEuhBAiq1RHpw8H/q6UAqfGP2mtH3K3pNIQNgy6bBu/UZL7X0IIIQZQSYa41nobMNvtOkpRKNulXul2IUIIIVwnzTmPiZgmcelOF0IIgYS454QNg4SEuBBCCCTEPSdsGCQsy+0yhBBClAAJcY/JDWiTEepCCCEkxD0obJrSpS6EEEJC3IuihkFMutSFEGLQkxD3oKiMUBdCCIGEuCdFTVNa4kIIISTEvShgGBgg66gLIcQgJyHuUVHTpFNa40IIMahJiHtUhYS4EEIMehLiHhU1TTozGbfLEEII4SIJcY+KmCYaWfRFCCEGs5I8i9kJa2+HRx4BpSCdPrJda+eSuw7g80Eq5Vx3TnV65GuOUs790mnn+vEufj9kMmAYzkWpY3/NnUY0dz13Mc3Xb+t9e+4+pklFOk1HPE5dKHTk8Xv/LkIIIcpWeYR4JgOHDzsBl8m8Pszyv/f7IZk8Euo5vcPe74euriM7Ase6RCLQ2nrke9t2Ln19b9tQUwOHDoFlHbnteJf6eti92/mZ7M9V+f20myZ1bW3ONq17hDxjx8KePT3D3+c78r3P13PnIHdbbS10dPTclv/1aNvC4SP/z7lt/b34/X1/lfOmCyHEUZVsiCulrgT+CzCB32itv3vUO9fVwXXXDVRpJaPSttkbj6OjUZRSTojndgxyYZ/7PpPpua33Jf9223Z6K/K3ZTJHLvnfW5azs5PJQCLhXDo7e94/d0mnj/xs3nXLNFly4AArtm1jbkUFV1RX8/ChQ6xob2duRQWLhg7FDASOhHt+0I8c6ew4tLYeuS3/kv9zueuBgLPjZZrOzwYCR27Pv360bX4/BIPOzwshhItKMsSVUibwC+CNwB7gZaXUfVrr9X3+QO9W9QCzLIvOVCc+04fW2glUQKFw/mW/z+sdyG0zlNFjW/fP5j2GOkoXud8wCBgGcdsmmgsk39H/pLYNH/sYHDzo5I/WTqeEaUJFBUSjzn2am52HGTECGhuPNJYDAadDoK0Nhg51sjq3vbYWdKfzPMGwk3E+n/M1EHC+RqNHcjQUci6mafHBDy5kxYplxOMxotEofr+fVCpFPB4nGo1yZl0de7dsYd7s2fhNk3cuXsziiy/GzN/ZSKdff8ndlkpBOo2VTLJkxQpW7NjB3LFjWTRlCmZHx5H75r6mUs59d+xgxeHDzK2sZFFtLWY6fWRnRSln5yX3C+YuweCRkM/fXlOT/SMr57b8/5jc9VCo57b87fnbct+bJpZlsWTJElasWMHcuXNZtGgRpuxcCDFoKN2PAFRKPQb8UGv9YN6227XWNxSlKKXOBb6ptV6Y/f6rAFrr7/R1//njxunlX/hCz9aX3+8kUyBw5Lh1X5ferbTeXbr96M5de2Atv1n5G5KZJJFAhHQmTdAXxLKdKWAhXwjDMMjYGRSKgBHA0lZ3gBuGgU854avRGMrAUAY+w4etbUzD+VA2lYmpTAxloJTCb/pps/yYSlHnd3YITGWCAp/yYSqz+34KhaEMNm82iHUYZFI+MmmDdNIkkzZJpwxSSYPWFpOg36SqwiSVgZDfTzRkYCiTgN/ENBSGYeIzFIcOHWLlqlUoFD5TobWJYQQxlIltg2mGSad8KBUgk/ERCUM67aezI0rG8pFJ+7FSj/Lqshux7fixXhFA/uvUYMLEi3jnux8kGjWJhCEaVVRVg89URCJOQ7uyUhEKQWWFIhSyeM91V7J8+UvEYs7OwoIFC1i6dOnrQs+yLBYuXMiyZcuOft9cr0c29EmlnD2idNr52nu7ZTnh39XlfJ9/yd2nr+vBoLNntXOnsy33811dWIbBwt27WRaPE8vuyC2or2fpJZdg+v3O4Y3c3lIo5HwfDB7ZHg4fuZ7bnrtEIkfukxtzUWRKqVe01vOL/kRClJH+hvg2YDfwuNb6W9ltr2qtzyxKUUq9A7hSa/3h7PfvAxZorT/V1/3nz5+vl7/8co+WVPf13Adtfgutd8tN6x4tth5dwHAk3HMto/xjyHkXbZpYCmzTwDIVlmlimQbaAMsAy1TYhiKjNJbhXLcM5/62UljawrItLG1haxuNJmNnsG27+zZb29jYZOwMGTtDyjbYl04zwkxjaxtDGXRlusjYGXyGj0Q6ga2dn9dak7bT3Y8X9AXpSHZgY2PZFkEjSGe6k1TGJpW2CBmVHI41k05nSFkay7JBmyRSKdCKzo4EB/Y0k9EZ0GBnQKeiZIwU2s5gYJJJO5//2k7jMzWZjI1f+8ioJAYZ2lbvon31zhN/kRiKqtmzMWvGoKnEtg18ZgCdGYJh+rEyYQL+JIlEPXY6TKp5LekdvwI7deR1ZoaYcO4Xqag/l6qKIEOHgqF8HN75Is/d/00yqUT3fcORMN+77XtccsUb8Bs+Z4fJMLt3uHLXFQqf4evz+/xelfweGDjSM5N9vfe5vcftWvPAvfdy3XvfS2cs1n1bJBTkh1/9N648cw6+VBpf2sKfzGCm0phdaXwZC7MriZFMYyZTzmPF4073Szzu7GjkvnZ1OdeTSee1X1HhdLlkMnTvKeVfotEj16urj+xIRyJHunryv4bDPXYOJMSFOHH97U5vBS4DfqqUuh94b9Eq6iel1A3ADQBjx449cmwzECjsE+WOJ+cCv69jxNljwyqTwZf7PpUBnTpyfLn3cee+vg+FnB0CpY7sKLzuEuj5vc/HBttmgi9EyJfdFs4b6d5rRPvrvuaPiO8HrTUa7exMaLt7x0Kjj+xkdN9mk0zbJLoypNKaVNrOXiySaZuMZfP8k4/zo81foitxrJZ4X38XzZjQME6/4HzCFREs28ayUygb0iqGrbuw7BSW3oGpFWufep7N21I9HkJbXUSqnmHiBRD0VdNlxYinbHaseLJHgAMkEgn+9PBdrK/ZhM7+bjWhWg7FD6C17t7xqg7W0BhvxNYWttYoI7vzhd19qGVoZChN8SaUUoR8IdJWuntnwFAGteFa2pJtNFQ0cKDzAL7sToNCMbxiOI2xRoZXDOfhvz5MZzzWo854V5LPP/9Lvh6KdO/w1YXqaEo0MSQ8hAPxA9jaxhfwkTJTDI0MpTXaijHMYGTlSFKZFLFMjBHREWR0hlgqyYjwMEIZ+MmFtzKndvqRkI/FnOu5S+77RMKZMdLc7BxzicWcr/nXYzGYOBG2bYPKSucihDhh/Q1xpbXOAJ9QSl0PPAvUFq0q2AuMyft+dHZbN6317cDt4LTEi1aJYRRn56C3XPfs0S65nYk+vtak07Qlk4Ryo997D27ra1uuWzi3g5E/sr2v6W7Z71UggDJNjNyI/+NNi8tdfCYE/D3vrxRXTr6OZQ/+lWXLlxOLx4lEIgQCAVLJJPFEAp/PRzp/2mBWKBTiP77ybyxcuLjHr5bbh8r/77RteHT0A3zmpeuIxzu7HyMcruCGKz/PxRcv7jEm8OnAA9z84iskEnn3DVXwjjlf5bz6xd0zBQHsmp6zDXubMcNpdDp/Yt29s5PbEcpdP9623CW3w6C1ZsxVY9hw7wbisSM7QOFImFuvu5XzLz+/+2cVqmdPjra760lkEt29OrZtk7Cc77Wt6bK6jvTc2BbDJ5wOlQ2n9jrPl8k4Yd7R4YT+aacV7rGFGCT6253+Ua31r/K+nwd8Umv9waIUpZQP2ITT+t8LvAy8R2u9rq/7z58/Xy9fvrwYpXhC3LLYnUwyLRI5+QfJnwrXO/h7b4P+TYvL7RzkvuZ+Lv+5bBsrnWbJiy+ycvNm5kyaxBVz5/Lw8uWs3LqVR155hec2bMDWmtxr1TQMLp4zh4d/9CPn2O+x5uNnezYs22bhRz7CstWriSUSRMNhFsyZw9I77sDMTWXLJrGlNQvf8x6WrVhBLB4nGomw4MwzWXrPPc59c4mdPzdfKTRHtikjbw2C3pcC6dexew+R7nQhTly/QtwNSqk3AT/BmWL2O631rUe772APcYDXYjHGh0KEPPjhfSxr1qxh2NChvPjCC9z5xz+SSCT4yAc+wOIrr8TMTavL3zHo6ys4OwuZDEsef5yVa9cyZ+ZMFl188ZHHyL+/1s6o76efZuVrrzFn2jQWnX++E4z58//7uuTfB/q+D/QM9fz58PlB3zv4869ne0e669yw4cjvlJuh0Hsxo96HTXqvpZB/vT87HUd73KPdN3f/o3zmqNpaCXEhTlDJhviJkBCHA8kkACOCQZcrEcfVO9BzS+f2Dv6+7pt/vffP9H783rfnB3Jfj9H7+tECvPdnxomG+FEeW1VXS4gLcYJKcp64OHG1fj/bEwmGBwJHnVcuSkTvQCuz3hMhxMCRNS3LRDC78IucnlQIIQYPCfEyUuvz0SynJxVCiEFDQryMVPt8xC1LTk8qhBCDhIR4GTGUotrno7WPedVCCCHKj4R4manz+WixLMph1oEQQohjkxAvMyHTJGgYtMmxcSGEKHsS4mVoiM9Hk4S4EEKUPQnxMlSRXbGrU4JcCCHKmoR4mRri89EsA9yEEKKsSYiXqWqfjzSQkMVfhBCibEmIlymlFEN8Pg5La1wIIcqWhHgZq/b5SGlNXFrjQghRliTEy5hSinq/X1rjQghRpiTEy1yVaZKxbeIyUl0IIcpOyYW4UuqbSqm9SqmV2cub3K7Jy5RSDAsEaMxkZBU3IYQoM6V6PvEfa63/0+0iykWFz0drJkNbJkON3+92OUIIIQqk5FriojiGBgI0ZTLY0hoXQoiyUaoh/iml1Gql1O+UUrVuF1MOgoZBhWnSJIPchBCibLgS4kqpR5VSa/u4XAP8NzAJmAPsB354lMe4QSm1XCm1vLGxceCK97B6v58OyyIp5xsXQoiyoEp5sJNSajzwgNb69GPdb/78+Xr58uUDU5THtWePjY8JhdwuRYgelFKvaK3nu12HEF5Sct3pSqmGvG/fCqx1q5ZyVOXzYShFq3SrCyGE55Xi6PTvK6XmABrYAXzU1WrK0DC/nz3JJFHTxG+U3H6cEEKIfiq5ENdav8/tGsqd3zCozZ7lbHgw6HY5QgghTpI0wwapap8PpZScc1wIITxMQnyQUkpR5/PRblmkZbS6EEJ4koT4IOYzDKp9PpplSVYhhPAkCfFBLmqa+JSiTbrVhRDCcyTEBbU+H0nbJiHnHRdCCE+REBcYSlHn99OWycjxcSGE8BAJcQE4086qfT5a5CQpQgjhGRLiolvYNAkaBq1yfFwIITxBQlz0UGWaaK3pkCAXQoiSJyEuelBKUev3k7Rt4jLQTQghSpqEuHgdIxvknXLaUiGEKGkS4qJPplLU+nwyYl0IIUqYhLg4Kr9hUG2atGUyZCTIhRCi5EiIi2MKmiYVpklrJoMlU8+EEKKkSIiL4wqZJlEJciGEKDmuhLhS6p1KqXVKKVspNb/XbV9VSm1RSm1USi10oz7xemHTJJydQy6LwQghRGlwqyW+Fngb8HT+RqXUTOBa4DTgSuCXSilz4MsTfYlkg7xNWuRCCFESXAlxrfUGrfXGPm66BrhLa53UWm8HtgBnD2x14lgipklIglwIIUpCqR0THwXszvt+T3abKCFh0ySSDXIZtS6EEO7xFeuBlVKPAiP6uOlrWut7C/D4NwA3AIwdO/ZUH06coJBpopSi3bKoAnxGqe0PCiFE+StaiGutLz+JH9sLjMn7fnR2W1+PfztwO8D8+fOlX9cFQcPAADotiwgQkCAXQogBVWqfuvcB1yqlgkqpCcAU4CWXaxLH4DcMKn0+4pZFQtZaF0KIAeXWFLO3KqX2AOcC/1RKLQXQWq8D7gbWAw8Bn9RaSzKUOFMpqnw+UloTkyAXQogBo3QZjDCeP3++Xr58udtlDHpaazqzIV6RPWYuRH8ppV7RWs8//j2FEDml1p0uPEwpRaXPh08pOixLpqAJIUSRSYiLgsut7ianMhVCiOKSEBdF4TcMKk2TlG0TtyzK4bCNEEKUGglxUTSGUlSYzqq5ndK9LoQQBSchLopKKdW9VGvcskhJ97oQQhSMhLgYEH7DIGqapLUmbllyJjQhhCgACXExYAyliJomplLELIu0tMqFEOKUSIiLARc0DCKmSUprumTQmxBCnDQJceEKM9sqN5QibtvSKhdCiJMgIS5cFTAMwoZBRo6VCyHECZMQF64zlCJsmgSUImHbJKWLXQgh+kVCXJQMn2EQMQyUUnTZNhnpYhdCiGOSEBclRSlFwDAIGgYW0CWLxAghxFH53C5AiL4YShFUCktr0rZNRin8SmHImdGEEKKbtMRFSTOVImSamEDStknJ8XIhhOjmSogrpd6plFqnlLKVUvPzto9XSiWUUiuzl9vcqE+UHp9hEM5OSUtlW+cS5kKIwc6t7vS1wNuAX/Vx21at9ZyBLUd4hc8w8AEZ2yalNYbW+JRCSTe7EGIQciXEtdYbAPngFSdNwlwIIUrzmPgEpdQKpdRTSqkL3S5GlDZfdiS7AmcAnG3LgjFCiEGjaC1xpdSjwIg+bvqa1vreo/zYfmCs1rpJKTUP+IdS6jStdXsfj38DcAPA2LFjC1W28Cif4eyP2lpjZS+mjGYXQpS5ooW41vryk/iZJJDMXn9FKbUVmAos7+O+twO3A8yfP1+aXgJwpqYZSmFrja01mVyYI4dvhBDlp6S605VSQ5VSZvb6RGAKsM3dqoQXGUrhMwz82eBOa40lI9qFEGXGrSlmb1VK7QHOBf6plFqavekiYLVSaiXwf8DHtNbNbtQoyoNSCjO7CpwCbOg+bi6BLoTwOrdGp/8d+Hsf2+8B7hn4isRgYGSPm+vsxdIaBSiQY+dCCE+SZVfFoKOU6g5urXV361zJsXMhhMdIiItBTSlFdhBG92C4XOs8d7sQQpQqCXEhsgylIBvaWmun2z0b6iCBLoQoPRLiQvQh1+UOdA+AsyXQhRAlpqSmmAlRilR2Odf8wW+yKpwQohRIS1yIE5Brgfe3HW5ZFkuWLGHFihXMnTuXRYsWYZpm8QoUQgwqEuJCFIllWSxcuJBly5YRi8WIRqMsWLCApUuXSpALIQpCutOFKJIlS5awbNkyOjs70VrT2dnJcy++yJ333ed2aUKIMiEtcSGKZMWKFcRisR7buuJxPvHgg3x9xAjOrqxkQVUVZ1dVMa+iggqfvB2FECdGPjWEKJK5c+cSjUbp7Ozs3lYRjfKnN7+ZmXPm8FJ7Oy91dHDTtm2sjsWYFApxdlWVE+yVlcyMRjFlFLwQ4hhUOawfPX/+fL18+etOdCaEq07kmHjKtlnV2clLHR0sa2/npfZ2DqRSXFhdzdRIhBmRCDOiUaZHIgzx+136jYpLKfWK1nq+23UI4SUS4kIUUW50+sqVK5kzZ84JjU5vTqdZ0dnJqs5O1sdivBaPsyEeJ2QYzIhEmB6JMDMa7b4+IhDw9Px1CXEhTpyEuBAeorVmfyrFhnic1+Jx1sdi3dfTWjM912rPa72PCQY9cYIXCXEhTpwcExfCQ5RSjAwGGRkMclltbY/bDqdS3a31DfE4D7e0sCEWo9bnI2HbjA4GGR0MMioYZEz2a25bvd/v6Va8EIOVhLgQZaI+EOCCQIALamp6bE9YFnuTSfamUuzu6mJvKsW6eJylLS3sTSbZk0zSaVmMDAR6BH3v0B/i93uiRS/EYOJKiCulfgBcDaSArcAHtNat2du+CnwIsIBPa62XulGjEOUibJpMjkSYHIkc9T59Bf36bGt+TzLJ3mSSDsvi0poaDqfTVPt8zsU0qfL5qPH5qDJNanLbs5fctkrTxGfIshRCFJpbLfFHgK9qrTNKqe8BXwW+rJSaCVwLnAaMBB5VSk3VWlsu1SnEoNCfoI9bFgdSKVoyGVozGdpzXy2L1kyG7V1dPb5vz2Roy96nw7KImCZVptljByA/8IUQJ86Vd47W+uG8b18E3pG9fg1wl9Y6CWxXSm0BzgZeGOAShRC9REyTieHwSf2srTWdlvW64M//Xghx4kph9/eDwF+y10fhhHrOnuw2IYSHGUpR5fNR5fMx+ij3+fKAViREeShaiCulHgVG9HHT17TW92bv8zUgA/zxJB7/BuCG3PfRaPQkKy0y23YuhuFchBgoWjsX8MRrb8iQIfPmz5/v/TmvXmDbkEpBKOR8H4uBaR75/hRYWrMvlaLTshgTDFIhJ/vpt1deeeWw1nroifxM0UJca335sW5XSl0PLAYu00cmq+8FxuTdbXR2W1+PfztwO5ToPPHDh+GFF6CuDubNK8ibQ4gTZttw4IDzehw5Eurr3a7oqObPn0/JvY/LVWcn7NsHU6c6369c6XydM+eUHvbJlhY+v3Ur19fW8vVx4+R8ACdIKbXzRH/GrdHpVwJfAi7WWsfzbroP+JNS6kc4A9umAC+5UOKpWb/eeVOcey5MmOB2NWIwMwwnvOvqYPdup8U1erTT6hKDl9aQH7CmCfH40e9/HB2ZDN/asYOn2tr48eTJXNRrmqMoHrf62H4OVAKPKKVWKqVuA9BarwPuBtYDDwGf9NTI9FQKHn8cNm+GxYslwEXpCIVg8mTng3vzZkgk3K5IuKn3fP9AANLpk3qop1pbeUO2Jf/47NkS4APMrdHpk49x263ArQNYTmEcPgxPPAFjxsDFF0tLR5QepWDUKGhthe3bYcQIp4UuBh+tewa533/CId6ZyfD/du7ksZYWfjhpEpf0WkFQDAw5YFEIue7z886D8ePdrkaIY6upgXAYdu50utdHjfLEwDdRQL3PmXGCIf5Mayuf27qVC6ureWLOHKrk2Ldr5H/+VKRS8Oyz0NEBV18NlZVuVyRE/wSDTvf63r2wZQuMG+dsE4PHSbTEY5bFt3fuZGlzM/85aRKXSuvbdbL7fbIOH4Z774VIxDn+LQEuvMYwnMM/9fWwdavTzS4Gp36E+PNtbVy6ciUJy+LJOXMkwEuEtMRPhnSfi3JSV+fsjOa610eOfP3AJ1FeTuCYeNyyuHXnTh5sbuYHEydyuYyjKCnSEj8RmQw89ZQzuvfqqyXARfnIjV7PZJzu9VTK7YrEQPL7+9xxe7GtjctWraI92/qWAC890hLvr64uePRRp9t88WIZfS7Kj2k6x8YPH4Zt25yu9lJdCVGcuvzQ9vmcHbishGXxnV27uL+pie9NnMgVEt4lS0K8Pzo74cEHnXnfZ53ldjVCFFd9vTN6fc8ep2tdxnuUp/wZCT4ftLeD1rzU0cFnt2xhbkUFT8yeTY3f716N4rgkxI+npQWWLIFZs5yLEINBNOq0ynftAstypqWJ8qUUXYEA3922jX80N/PdiRO5csgQt6sS/SAhfiyHDsHSpbBgwZE1hoUYLEIhZ9zHzp3OGuzSpVpe8rrTl7e385mGBs5IJnl8zhzqpPXtGRLiR7Nnj7OE6iWXwNixblcjhDsCgSNBblkw9IROsCRKXJdl8f3du7mnsZH/iMe5atQoZ5Cb8AwZnd6XrVudJVSvuEICXAi/3wnyjg7njGiiLLza1cUbV69mb7b1fZVtn/T66cI90hLvbf16ePVVuOoq6T4UIsfnc4J81y5nlTeZS+5ZSdvmPw8c4C9NTdw6dSpX505P6/f3GKEuvEFa4vleeQVWr4Y3v1kCXIjeDMMZ7GbbzmlNe6+/LUreyo4Orli1iu3JJI+PHn0kwMHZUZOWuOdISxycD6Pnn3e6Ct/8Zmf1KiHE6ynlnI983z7nOPnYsXLyFA9I2TY/3L2bPx06xC3jx3ONz4fqff7wkziTmXCfvPu0dlZha211VmGTABfi2HKnNA2FYMcOZ8CbKFmrOju5YtUqNicSPD57Nm8ZOhSl1OsPh4RC0p3uQdISf/55iMdh4UKnO0kI0T8jRjjTMPfscVZ3kxZ5SUnbNj/as4c7Dx7kW+PH89b6eie8wQnw3n8vw5CWuAe59q5TSo1RSj2hlFqvlFqnlPpMdnudUuoRpdTm7NfinSpn5Urn2N5ll0mAC3Eyhg2DigrnUJQcIy8Zazo7Wbh6NRtiMR6dPZu35VrfOYbx+qWjlZI18z3IzV3nDPB5rfVM4Bzgk0qpmcBXgMe01lOAx7LfF97mzc4gtsWL5TzKQpyKIUOcQGhsdLuSQS9t2/xg1y6uW7+eT44axf9Mn87wQKB/P9xr/XThDa6FuNZ6v9b61ez1DmADMAq4Brgje7c7gLcU/Mn37IFnn3UCvKKi4A8vxKAzbJgTAE1NblcyaK2LxVi0ejVrsq3vt/dufR+PTDHzpJLoQ1ZKjQfmAsuA4Vrr/dmbDgDDj/IzNwA3AIw9kQVZDh+GRx5xjoHLNDIhCkMp5xj53r1Oi6662u2KBo20bfOzvXv53f79/Pv48bzzRMM7R6aYeZLrI1GUUhXAPcBntdbt+bdprTXQ54E2rfXtWuv5Wuv5Q/u7FGRHB/zzn3DRRc5iFUKIwjEMaGhwZnrEYm5XMyhsiMV405o1vNLRwaOzZ/OuYcNOLsBButM9ytWWuFLKjxPgf9Ra/y27+aBSqkFrvV8p1QAcKsiTdXXB/ffD3LkwaVJBHlII0YvP57TIDxxwjpOHQm5XVJYyts0v9u3j9n37uHncOK49lfDOkXninuTm6HQF/BbYoLX+Ud5N9wHvz15/P3DvKT9ZJnPkfOBnnHHKDyeEOIZg0DlRysGDEgpF8FosxlVr1vBCWxuPzJ7NdcOHn3qAg7TEPcrNlvj5wPuANUqpldltNwHfBe5WSn0I2Am865SeRWvnGHhlJZxzzik9lBg4lm1haQvbtrGx0Vqj0ei8aUz5H1yGMrovCoWhDEzD7OuhxUCIRJwxJwcOOIeuek9nEicsY9v8ct8+frVvHzeNG8d7CtH6zqOzx8RlRXxvcS3EtdbPwlFfL5cV7ImeecZpDSxcKCdsKCGWbZGyUqTtNGkr3f3V1ja2tlFKYSoTn+HDNEwUCqVUd0Dr/KESGrIR3/3zAFrr7jDPPVb+pZAfgKIPlZVOy+7AAedYuSwGc9I2xeN8ZssWqkyTh2fPZlQRpsVua99FfWczMiTRW0pidHrRrFnjzF29+mr5AHGRZVt0ZbroynSRtJIkM0ksbRE0gwTMAH7TT0WgAr/h7w7cQgSs1hpLW1i2ha1tLG2RtJLE0jEs28I0TAJmAJ/hw2/4JdiLobbWCfJDh2D4cNmRPkGW1ty2bx+/3LuXr44dy78Uquu8L34/tiXd6V5TviF+8KBzVrK3vhX6u9iBKAhb23QmO0lkEsTTcdJ2mpAvRMQfoTJQSX2knoBZ/L+JUgqfclrdfclYGTI6Q9pK05XpwrIt/IYfv+nvDncJ9QKor3fej21tUFPjdjWesSUe57NbthA2TR464wzGFHKQYB9rpxs+P2RkDIPXlGeIp1KwdCmcf77MVx0gKStFR7KDzpQT3hFfhIpgBQ2hBoJmsCTD0Gf68OEj5HM+HLXWZOwMaTvd3Vr3m36CZhC/6cdQ0ptzUpRyBrodOOCMVpcR68dkac2v9+3jZ3v38qWxY3nf8OEYhX7/9BHiyueHRKKwzyOKrjxD/MknndMlTpnidiVlLZlJ0pZsoz3Zjq1tKgOV1IXriAaingw8pRR+02mJR/wRtNakrTQpO0U8He/ufg8YAQw5PHNiTNNZnrWpSY6PH8O2RILPbtmCXymWnHEGYwdoh0drjQr40WlZO91ryi/EX3vNOf527bVuV1KWLNuiJdFCa7KVjJ2hOljNyIqRRALldwpXpRQBX4AATtd/2ko7PQ6Zju5A9xv+kuxlKEmhkDNqvanJaZmLblpr7m5s5P/t2MHnx4zh+hEjCt/6Pg6fP4QtUwI9p7xCvLXVGY3+1rfKWckKLJ6O0xRvoj3ZTmWwkhHREUQD0UEVYLlWeq7bPWWlSKQTBMwAQV/Qk70PA66mxjk+3tkp5y3I6rIsvrB1K+tiMe47/XQmRdzZITZ8frRtg21LT4mHlE/SWRY89BCcfbYzkEacMlvbtCRaaE40Y2mLIeEhNFQ2HHWg2GCR3+2utSZlpYilYj1Gu4ujUMrpVj940FkUxu93uyJX7enq4gMbNzIpFOKfZ5xBxKX59BqNaZho03BmE8hgYM8on0+bF16AaBRmz3a7Es+zbIvD8cMcjh8m7A8zomIElcFKt8sqSUopgr4gQV+QtJUmmUnSRVf3YDjRB7/faZE3NQ3qaWfPtrby8c2b+eTIkXx05EhXe7VyayqkfUpC3GPKI8RTKdi0Ca67zu1KPM3WNodih2iMNVIVrGJi7UTC/rDbZXlGrnWev5BNbsqa6KWiwjmfQWurM5d8ENFa86t9+/j53r3cNnUqF7g17U73PLeUaZikTFOWXvWY8gjx9na44goIS+CcDK01h+OHOdB5gIpABVOHTCXoK/yKUIOFaZiEjXCPMJdu9j7klmUNhQbNezdhWdy4ZQtbEgkeOuMMRpfQdDuFws51pwvPKI9PlXDYmVImTlhTvIn9nfsJ+8JMrpssLe8CyoW5rW1SVoqMnSFgBmQAXI5hOMfHDx92znxW5nZ1dXH9a68xIxLhgVmzCLm5nnwfXfemYWL7pCXuNeUR4tGo2xV4TjwdZ1fbLrTWjK8ZT0VARgoXi6EMQr4QGTtDMpPEUAZ+wy9zzcEZ3FZR4RwfL2NPtbbyyU2b+Mzo0Xy4oaEkZ3U4LXFTzjznMeUR4qLfLNtiT/seWrtaGVU1ivqIjOQfKLkTr+RO+GJiShc7QFWVs7ZDGdJa88t9+7ht3z5unzaN80p4BUlTmdgKCXGPkU+QQaStq43trdupDlZz+rDT5VSdLslNTUvbzmj2Qb+kq1LO8fEyE7csPrtlCzu7ulh6xhmMLMKZx06aUq+bC66UckLcstypSZwUCfFBIGNn2NW2i45kBxNrJ1IVrHK7pEFPKUXADGDZFhk7g6nMwb1TVWbzxXckElz/2mvMrqjgvlmzCJbaoZOjdOdrGdjmORLiZa61q5UdrTuoDdVK67sEmYaJoQ0sbZG20nLmtDLweEsL/7Z5M1/ILp/qhb+nUgrDMJxj4hLiniIhXqZsbbOrbRctiRYm1U2S1ncJy50y1da2c6x8sLfKPUprzU/37uW3+/fz22nTOKeEj3/3prVGobBMJMQ9RkK8DCUzSTY1bSLoCzJr+CwZPOURuVHrGTsDNhLkHtKZyfCZLVvYm0zyUKkd/+4HhXLGZfh8EuIeU2IHasSpautqY/XB1QwJD2HqkKkS4B6TW5ddKYVlW+heq2qJ0rMtkWDRmjVU+XzcN2uW5wIcAOW89ixDoWV0uqfIJ3wZ2d+xnz3te5g6ZCrVIe905YnXM5SBRmNrGwPDE8dVB6NHmpv5zJYtfGXsWN43fLi3/k597SD6fOhMGg/9FoOehHgZsLXNtpZtdKY6mTV8FiFf6SzlKE6eUsqZu6tt0HgrIMqcrTU/2bOH/z1wgDumT+esKo+NOen1WlK52JbudM+REPe4tJVmY9NG/IafWcNmyXHUMmQoQ4K8hHRkMnxq82YOp9M8dMYZjPBi9/lRKL9futM9Ro6Je1hXpovVB1dTEahgWv00CfAyllsMRo6Ru2tzPM6i1asZHgjwj9NPL6sAB7BksRfPkZa4R8XTcdYeWsuYqjE0VDa4XY4YAEopCXEXLW1u5sYtW/jauHH8y/DhbpdTcEoplN8PSelO9xIJcQ+KpWKsObSGCTUTGF5Rfh8m4uikO33g2Vrzw927+dOhQ9w5YwZnVla6XVJh5O0Q5l5XGUNho5E+Pe+QEPeYXIBPrJ3IsOgwt8sRLrMsiyVLlrBixQrmzp3LokWLMN08xWWZac9k+MSmTbRbFkvPOINhgYDbJRWGUt2D2/JfQ5XBODMnzaW8FsEtbxLiHtKZ6mTtobVMqp3E0OhQt8sRLrMsi4ULF7Js2TJisRjRaJQFCxawdOlSCfIC2BSPc/1rr3FxTQ3/b/x4/KW2/vmpyJ4ApfdrKBQKct/EiTzylrfLa8gjyuhVWd46k52sObiGyXWTJcAFAEuWLGHZsmV0dnaitaazs5Nly5axZMkSt0s7IWnbZldXl9tl9PDPpiauWbuWT48axXcmTiyvAM/T+zWUSHTx0patnnsNDWbl+cosM/F0nNWHVjO5brKc/1t0W7FiBbFYrMe2WCzGypUr3SnoBMQtiy3xOE+1tnJ/UxM7SyTELa35zs6dfH37dv48YwbXluEAtnx9vYbiqaQnXkPCId3pJS5lpVhz0BnEJi1wkW/u3LlEo1E6Ozu7t0UiEebMmeNeUcfQmk6zL5VibzJJzLIYGQwyKRTi/KoqfCXQ0k1YFjdv3862ri4ePuMM6svl+HdflAKt+34NBYIl+xoSrychXsJsbbP64GqGRobKNDLxOosWLWLBggXdxzN9QR+RCRGuWHiF26UBzpz2xnSavckk+1IpAEYFAsyuqGCo319SI+3jlsV7N2ygIRDg7pkzy7b7vLfca+jp554mk8wQCgc5e8xYFi1a5HZpop8kxEvYhsYNhH1hJtZNdLsUUYJM02Tp0qUsWbKElStXctqs0/hZy8+4+cmb+f4bv+9KTRnb5kAqxd5Uiv3JJFHTZFQwyPlVVdT4S3PMc9yyeM/69YwJhfjJ5MmYJbRzUWy519BFN13E+K7xnDl9Jp9IhmRQm4dIiJeoHa07SGQSnNlwptuliBJmmiaLFy9m8eLFAFycuJg33PEG5o6Yy3WzrhuQGrosq7ubvDGdZojfz6hAgFnRKJESD4OYZXHd+vWMD4X48SAL8BzTNJl+/nTePuPtDFFRzPtecrskcQIkxEtQY6yRve17mT9yfvdym0L0R124jrvfcTeX/u+lTKufVrSdwI5Mhr3JJHtTKdozGUYEAowLhTinqsozXdGdmQzXrl/P5HCYH02ejDHYAjxvsZdDsUOYykT7TFl21WMkxEtMZ7KT1w6/xpwRcwj6ymtdZjEwptVP42eLfsZ191zHCx96gbpw3Sk/ptaapryBaRmtGRkMclokwrBAwHMB2JHJ8O7165kRifCDSZM8V/8p6/X7ZuyMc+4F04S0LLvqJRLiJSRjZ1h5YCVTh0ylMlgmSzsKV7xtxtt4fvfzfPDeD/K3d//tpHp0LK05mEqxLzswLWgYjAwEWFBVRa3PV1ID005EezbAZ0WjfHfixMEX4H2wbAuf4UMpA20aTmu8xA+FCIc3+r0GidcOv0Z9pF7WQxcF8Z3LvkNrVys/eO4H/f6ZlG2zs6uL59vauO/wYTbG41T5fFxaU8PCujpmVVRQV2Ijy09EWybDO9etY3Y0yvckwLtZ2jpypjzTkHOKe4i0xEvE/o79tCfbOWf0OW6XIsqE3/Tzp7f/ibff/XbOG3MeF467sM/7xSzLmQaWTNKcyTA8EGBUIMC8ykqCHjm+3R+t6TTvXL+esyoruXXCBM/uiBTDqMpR+A0/lrbQFZVOiJfZaVbLlYR4CUhmkmxs2siZDWfKQDZRUCMrR/KNi77Bh+77EC9/5GWqQ9UAtOTN307YNiMDAaZGIgwPBMpyhHZLOs07163jnKoqbpEA73ECFHBmw5iGiWVbYCppiXuIq4mhlPqdUuqQUmpt3rY6pdQjSqnN2a+1btY4ENY1rmNM1RiqglVulyLK0JVTrmTh5Cv51MNfZXVnJ0uamlje0YEFnFlRwZuHDOGsqipGBoNlGeDN6TRvX7eO86urJcDz5fWyZOwMpnKOgauupIxQ9xC3m32/B67ste0rwGNa6ynAY9nvy9be9r0kM0km1sqCLqKwtNYcTqVY2dHBJXO/zMa0YtmeFzivupo31tUxu6KC+kCgrEOtOZ3mbWvXcklNDd8cP76sf9cTljfFbHjFcGd0OqBNU1riHuJqiGutnwaae22+Brgje/0O4C0DWdNASqQTbGraxKzhs+TDRRSEnR1R/mpHB/9samJlZydBw+CN9cO5dc5b+Mnjn8WvU26XOSCa0mnesnYtb6yr4+vjxsl7LF+v/4u97Xsxlen8H5kyV9xL3G6J92W41np/9voBoGyHam84vIHxNeOpCFS4XYrwMFtr9ieTLG9v54GmJtbFYkRNk0tqari8ro4Z0ShVPh9vnPRG5jXM4/vPubMk60A6nErx1rVrubKujpvGjpUAPw5LW5iGiUI5C75IS9wzSnpgm9ZaK6V0X7cppW4AbgAYO3bsgNZVCPs79tOV6WJ8zXi3SxEeZGntrFGeTLI/laI6u0b5zOMsdfq9y7/Hgt8s4NrTr2V6/fQBrHjgNGYD/Or6er40ZowEeD9YtoWpTGxtgyEtcS8pxZb4QaVUA0D266G+7qS1vl1rPV9rPX/oUG+dotOyLTYc3sDMoTPlA0b0W8a22d3VxYttbTzQ1MTWRIIhfj8La2u5pLaWKZHIcdcqb6hs4KYLb+LTSz7tfGCXmUOpFG9Zu5Zr6uv5srTA+83WNoYy0GhsU0mIe0gphvh9wPuz198P3OtiLUWxo3UHNaEaakI1bpciSlw6b/GVB5qa2NnVxfBAgCvr6riopoZJ4TChE1xZ6yNnfgSF4u8b/l6kqt1xIJnkmrVrefvQoXzRg71zAyp7PvGc0VWju6e3Kp9futM9xNXudKXUn4FLgHql1B7gG8B3gbuVUh8CdgLvcq/CwktZKba3bufc0ee6XYon5BoEuQZVr+mtZSlp2+xLJtmbTHI4nWZYIMCoYJD5lZUECrD4immYfPvSb/Pev7+XRVMWEfFHClC1uzoyGb68bRvvqq/nxjFj3C7Hc3a07sBvOqeKldHp3uJqiGutj3auxMsGtJABtLV5KyMqRhANRN0upSRkMtDUBIcPQ1ubc0kkIB53vlZUQDrt3FdrCIWcEI9EIBx2LjU1zqWuzrndi7osi73ZY9wtmQzD/f7us4L5irBq2lmjzuKc0efwy5d/yRfO+0LBH38gaa359JYtDPH7uVFa4CelrcNm5w4Dsw4Zne4xJT2wrdzE03H2tO/honEXuV2KaywL9u6F7dud8D54EGprYehQqK+HkSN7BrSvj1doMnkk5GMx5+vGjdDcDH4/DBsGw4dDQ0Nph3o8u9zp3mSSNsuiIRBgUjjMiAFaNe0bF32Da+66hhvm3eDphYZ+vX8/O7q6uO2MM9wuxZNSKQhadezYoZhcp5yWuIS4Z0iID6BNTZsYXzN+0J1i1LJg61bYts25DBkCEybAeec54X2iJ0sKBp1Lba+1/LSGzk6nVX/wIKxeDZWVMGYMjBrl7By4rTOT6W5xd1oWIwMBpmWXOx3ok3FMrJvIeWPO41fLf8UXz//igD53obzS0cEPd+9m6RlnlNU67wNGa2IxRZfRSGtzNg5kipmnSIgPkPaudprjzcwaP8vtUgZMayusWgXr1jlhfdppcPHFTgu7GJRyQruy0tlJsG1obIR9++CZZ6C6GiZOdGoZyLxsz2S6W9wJ22ZUMMhp0ShD/X7Xz6L1+fM+z6I/LuLjZ33cc+sVNKfTfOi11/jR5MmML9aLqpxlX3vJJCjTor3VwLYV2jCkJe4hEuIDZHPzZibUTuhe2rCcHTwIL7/sBOjEiXDdda9vNQ8Ew3C61YcPdz6Tdu+GtWudlv/EiTB6dPHCvC2TYU82uNNaMyoQcJY5LbHTeE4dMpWLxl3Eb1/9LZ855zNul9NvttZ8YtMmrq6v56ohQ9wux9NSKdDYDKkz6YzZ6FDQ2QMWniAhPgAS6QSN8UZmD5/tdilF1dTktHh373a6yq+8su9j2m4wTRg/3rkcPux072/dCtOnw4gRhXmOlnS6O7htYHQwyLzKSup8vpIK7t6+cO4XeOtf3soN824g7PdGi/ane/bQmsnw7+PGuV2Kd2WPY6XToLHQtkE8YaORs5h5SYl8xJa37a3bGVM1Bp9Znv/dySQ895wTirNmwVVXOQPMSlV9vXNpbIQNG2DnTifMq6tP7HG01jTntbgNnOBeUFVFbSn/B/Ry2rDTOGvUWfx+5e/5+Fkfd7uc43qurY3b9u3jsTlz8Mtx8JOX3bG0bbCwqKs1SXRlNyeT7tYm+k3eAUVm2Ra72naV7fKqW7bAr37lvOff+14455zSDvB8Q4fChRc6I+JXrXJGuB+vF1FrTVM6zeqODh5sbuaVjg78SnF+dTVXDhnC6RUVngrwnM+f+3ke3PIgGbu0W2CHUilu2LiRX06dyqjg4BogWiy2DTNqzsA0FF0JOQGK15Rn07CE7O3YS22otuzmhSeT8OijzlSxN7/Z6ab2IqWc0evDhzsD8J5/HmbPdgbH5WtNp9mdTLInmcSvFKMCAS6qrqayVI4XnKIzG84kbaVZsnkJV0+72u1y+mRpzUc2buS9w4dzqRuDLMpNdtW2TEaxoXUlVZN9NHWBDptgpd2uTvRTeXwClbDtLduZMXSG22UU1KFDcN99TvB99KPeaXkfSyAAc+fC/v3wyiswaRJUNzhd5XuyXYujg0HOr66mqkyCu7ePz/84S7cuLdkQ/96uXSjgS7KgS0HZNsyqn0cwCF1tOCdAkWPinlGen0YlojnRTMbOMDTirRO0HMv69fDAA7BwodNiLTeVQzNUB5L8Y2OSwGHNOZODnF1ZSU057Kkcx5WTr+TzD3+eve17GVU1yu1yenispYU/HTzIE3PmDMhCOIOJ1rCq8WUis33E46BNmWLmJXJMvIi2tWxjQu2Ekh6Z3F9aw9NPw7PPwvveV14BHrcsNsfjPN7SwlNtbRhBzbvPrOQcow5rcwVRyj/AAYK+IG+e9mb+uv6vbpfSQ0cmw2379nHb1KkMDQTcLqcsaWyiEQPbQo6Je4yEeJGkrBSHY4cZU+X9kzFoDUuWOIPYrrvOWc7U65K2zdZEgqdaW3m8pYUOy+L0aJQ31dUxu6KCYSE/Z5yhqKqCV18dPIN133Xau/jLur+g885w5bb/2rOHoX4/F9TUuF1KWbIsjUYTCioyGVnsxWukO71I9nXsoy5c131mIK/SGu6/3zlW/P73l/Za5MeTyp4dbE/2JCMNgQDTwmGGHWXJU6Vg8mTnmP/KlTBnjrPcazlbMGoByUyS1QdXM3uE+90t+5NJfnfgAM/Mnet2KeUnO7BNAwpFKKRoa7fQZlCOiXuIhHiR7OvY5/lWeC7A29rg+uu9GWAZ22ZfKsWeZJKm7Gk9J4TDnHsCJxkZN875vFu92jmMUM49ukop3n3au7lr7V0lEeLf2bWL948YIdPJiihjWxjKIBh0xrThk+50L5Hu9CLI2BkOxw8zoqJAS4G55KmnYM8eeOc7vRXgltbsTSZZ1t7OkuZm9iSTjAkGWVRXx4KqKkYFgyc8OGrsWOfsaGvXlv/n27tOexd/e+1vpF2eZrQhFmNJczOfHT3a1TrKndY2pmESDGYPG8kxcU+RlngRHOw8yJDwEE93pa9cCcuWwSc+4Y0udFtrDqVS7E4mOZBKUevzMToY5MyKioKt6jVunLNE5WuvwcyZA3sSlYE0oXYCl4y7hOd2P8cl4y9xrY5vbt/O50aPprpMp/SVDMNm1tDZhEJOrxs+n3Sne4i0xItgX8c+RlaOdLuMk7Z7t7OQy/vf//pFT0qJzgb3qx0dPNjUxKZEgiF+P1fU1XFBTQ3jw+GCL8s5aZJzmGHHjoI+bMmZNXwW/9z0T9eef10sRlMmwwfLYRRlKdOaTEbzWtN6TNM5w2BGy8A2L5Fd3AKztc2BzgPMGu7NU47GYvDHPzrrnxfqxCAnQmtNIpEgFouRSqXIZDIYhoHP5yMYDBKJREiYJru7utiTTBI2TUYHg1xWW0v4RE9MfhKUctZZX7PGOeFLuZ5Aa+GkhXzin59Aa+3KFMnf7N/P5bW1co7wYsr+XbWyGFc9Hq1tUilI2aa3jp8NchLiBXY4fpjKYCUhnwf6oHvRGv72N2fw1qwB2gfRWnPgwAH27t1LU1MTra2tGIZBRUUF4XAYf/bUnbFMht3JJPsyGQzTZEI0yowhQ5jQ0EB1JDIwxWb5fDBlirPwTUVFeX7eTa6bzIHYAXa27Rzwdf87MhnuaWxk2ZlnDujzDjpKgVLYtmZfx15+/OMfEw6fSTItZzHzEgnxAtvfsd+zXemvvuosqXrddcV/rra2NjZs2MCmTZuIRqOMHz+e008/nfr6eoLZVEzZNnuzU8LSlsX8YJDRwSARy6K5uZnm5mZefvllfD4fY8aMYcyYMQQGaOh4RYVzPvKtW53j4+VGKcX5Y87nuV3PDXiIP9naypV1dTSU495RKcn2cijDBg3pdJquZIa09kN7u8vFif6SEC+wxngj8xvmu13GCevogMceg3/5l+KeA7ylpYVXXnmFAwcOMGnSJN785jdTk7eIh601+5JJdieTNKZSDA8EmNJ7Lrffz8iRIxk5ciSnn346ra2t7N69m+eee46GhgYmTpw4IGHe0ADNzXDwoLOOfLm5YOwFPLvrWf7ljH8Z0Od9rq2NaQPcuzKYaQ1KGQQCAQI6TTIjLXEvkRAvoJSVoj3ZTnXoBE9MXQIeeMBpUY4q0pLZ6XSaF198kZ07dzJjxgwuvvhi/Nn1yHPn5d7d1cW+VIoq02RMKNTvkeU1NTXU1NSQSqXYtm0by5YtY9y4cYwZM6aox3OVgokTndHqQ4YUd+fHDeePOZ/vP/f9AT8u/mJ7O7dOnDhgzzfYaWwU4Pf7SXakSGuZYuYlZfax467mRDO1oVrPrZW+Y4dz+dznivP4+/fv5+GHH2bkyJG8/e1vJxwOAxCzLHZ3dbErmcRUijHBIG+oqTnpAWqBQIDp06czduxYNm7cyOHDh5k5cyahIs6Ri0Sgvt4Z0T9hQtGexhXja8ZjGiZbW7YyuW7ygDxn3LLYEI9zZkXFgDzfoKc1KBtDGRiGgWFYJFIS4l4iIV5ALYkW6sJ1bpdxQmwb7r0XLrusOAO0Vq1axauvvsoFF1zAlClTSNs2O7u62N3VRcy2GRUIsKCqqqBzgSORCHPmzGHPnj2sWrWKadOm9eiyL7SRI51FYBIJZ4pOuVBKccGYC3hu13MDFuIrOjuZEYkMyEyDQS/b2LCxUajsISgLKzfFzOlnd7dGcVwS4gXUlGhiXPU4t8s4Ia+95nydX4TD+E8//TQ7d+7kmmuuwaqo4NWODg6mUtT7/Uw+xprlhaCUYsyYMVRWVrJp0yYmTJjA0KHFOSWsz+ccH9+3z5lHXk4uGncRaw6uGbDne7G9nXOrvXc4yssUmkk1U/DZPnx+m7RlOIPebNtZvU2UNJmEWUDNiWaGRLwzcVhreOghuOii7oGqBfP444+zefduzrjqKlYrxYZ4nFqfj8tqazmrqooRwWDRAjxfTU0Np512Gjt37qSxsbFozzN0KMTj0NVVtKdwxZwRc3hk+yMD9nz7kknmSVf6gNJodnXsQCmFbVukUkqWXvUQCfECSaQTaK2J+L0zqnbbNifIC3lucFtrHnjxRZ5ubGTYpZfiyy59elFNDRPCYQIuLN4RjUY57bTT2L17N21tbUV5DsNwgvzgwaI8vGsm1EzAVCaxVGxAnm9XdgEfMXAs22JoaCg+nw+faWNlkBD3EAnxAmlKNHnuePiTT8JZZxWmFd6RybA+FuOOtWt5ZNUqrr3sMq4aMYLTolGq/O6vIR8Oh5k6dSpbt26lq0jN5fp6Z+3pcpqd4zN9BH1BdrTuGJDnqzRN2svpP9ADlKFpSbYAYFkZkikkxD1EQrxAvDaorb0d9u6Fs88++cewtGZ3VxfPtbXxYns78ViMxief5HNXXsn0+voB6S4/ERUVFYwePZodO3agtS744/t8UFPjzB0vJ5NqJ7GtZduAPFe1z0e7hMfA6D6fuAblTDHz+WxsCwlxD5EQL5DD8cNUBLxzLO/ll505zicz+6otk2FNZyePtrSwP5ViUijEZTU1bHr0Uc6eNYtx40p3cN+wYcPw+/0cOHCgKI8/ZAi0tBTloV0zoWbCgIW4tMQHnlIaAwOlFCmri0SXLSHuITI6vUDi6TjRQNTtMvpt/XpYuLD/909nl0DdlUyS0ZoxwSAXV1cTyh6/3LBhA4lEggsuuKBIFRdObh75kCFDCr6yWzTqjDMop+lmE2sn8vzu5wv+uFpr9qVSrI/F2BCPAzAqEGB/KlXw5xJHp3EW8/H7/fj82plVJiHuGRLiBRJLxzwzqO3wYWdxksn9mPrbkcmwo6uLxlSKGr+fGZEI9dmTkuTYts3DDz/Mm970JnweWLbM7/czbNgwDhw4wNixYwv++DU1zuGKcgrxP6z+w0n/vNaaQ+k062Mx1sfj3V83xGIEDIOZkQgzo84O8BkVFdy2eTPfGD/ec4smeZVhaBSqe3R6LKbA75cQ94jS/8T1AMu2SFkpwj5vfGqvXevMCz/agDatNQdTKXZ0dRG3bcYGg5xfU3PU00KuW7eOiooKJvdnr6BE1NfXs2nTJpLJZPcJVwqluhr27Cmf9dQn1k5ka8vWfi2/2tRHWK+PxdDQHdazolHePWwYMyMR6vN6QuYD51RVkbAsVsdizJapZsWVOyae/bs6K7bZzvLBPl95jdAsYxLiBRBPxwn7wp5pOaxZA+ef//rtKdtmV3YZ1IhhMC4UYkQgcNzf66WXXuLcc8/1zO8PYBgGQ4YMoaWlhREFPnF6OOx0qafTToPG63IDNlu6jgzebMvORugd1l22zYxIhNOiUWZEIlxTX8/MSITh/XgdgbNIz7XDhvHngwclxIst9/dQzoptpmkCNh0dgGE5i72IkichXgBeOh5uWc6o9GnTjmxry3aZH0qlGBEIML+ykqp+dosfPHiQAwcOMH369KLUm063YVmdmGY1fn9hP9RramrYtm0bw4cPL/gOSCTiLP7i9cXHWtJpXovHqZ3wL9y4eSORYB0Pt7TQlskwPRvWM6NRFtbVcVo0ysh+hvWxXDtsGJ/ZsoW4ZRGROePFpRRaa4ZHRqC1xrJTckzcYyTECyCejnvmePjevc4a6ZGI5kAyxbZEgqTWjAuFmFlT06+zhuVbs2YN8+bNwyjgIi7J5H4aG++ho2MVWqcwjCC23YXfP4TKyrMYMuRq/P7KU34ev99PJBIhFotRUeBWn9dC/HAqxYZ4/MglO9is07KYEYnQERxDhU5ydX09Xx47lrGhUNGmEE6PRmkIBLh5+3Z+5KFDNJ6jlNPaNjRhX4S9e39GOj2VeByoMSTEPUJCvABi6RhRvzda4tt3aKomJXmmrYugYTAxHO53V2df1qxZwzXXXHNKNVmWxZIlS1ixYgUTJnQxe/ZOamoWMHbsFwmHJ2MYBlrbxOObaG5+hO3bv86IEe+jqmreKT0vOHPHixHi4bAzuK2UaK1pTKfZkO36zg/slNbMiESYkT1u/aa6OmZEIowKBlFK8Yl//oYrQ3VcOWRglhX+0eTJnLl8OdfU13NxEU9eI5wpZvs69vDKKzt4dfMhMpMexrrAwJQQ9wQJ8QJIZBJUFLirt9Ayts2uZJKlB7sY1eBjVjRK7SkesI3FYqTT6VOaF25ZFgsXLmTZsmXEYjHCYZOzzprHY4/dgWmaWCmLRGMCNITqpzBmzHTa21dz8OAfsG2LmppTWK0GZ0nWeHZ6UyGFQs4xcTdOBGVrze5kko3xOK/F4+zs6mJFZ2f3NK6Z2bCeEY3ylvp6ZkQixx37EPFHSGQSA/UrUOf3899Tp/LZzZt5ZPbsHgPgRGFZVoa1P1nL8u1tJJN7eO35j3HghXEsvesu5GBG6ZMQL4BkJkltqNbtMvqUtG12dHWxJ5lkqN9PeGclbzjNR20BBlzt2bOH6urq7ICYk7NkyRKWLVtGZ2cnAPF4hldeWccD9z7AuYFz6drahX+4HzNsYidsgiODVM6fyahRH2Xfvt8RDk8gGDz5s5MFAgFs28a27YIeElDKGeCbTkOx8idl22xNJNgQj3cH9oZ4nE3xONU+H9MjEaZHIsytqOCabFj3nh7YXxF/ZMDWT89ZNGQIqzo7Wbh6NQ+dcQZDJciL4uXnn6Jjezt20lnFMJ2KsWzrVpY89RSLZ81yuTpxPCUb4kqpK4H/AkzgN1rr77pc0lFl7AyK0hqZnbZttnd1sT+VYpjfz3lVVYRNk3sTMGxYYZ5j//791NWd2lKzK1asIBbrGQ6xWIw7vncHj415mq43TMc8/VxMn4+pwTCXNipGPN1K9UXjqKo6i5aWRxkx4rpTqiEQCJBOpws+1SwQKEyId2QybMwG9GvZwN6QbWGPCYWYkQ3ry2pr+dSoUUyNRAp6fnZwQjyeLnyPxfF8eexYYrbNG1et4uHZsxkmQV5wW17bgJ3sORI9lkqx8rXXWOxSTaL/SjLElVIm8AvgjcAe4GWl1H1a6/XuVtY3rTWGKo0VbDO2zc5kkl1dXYwIBDinqqp7frdlwc6dcIq5262jo4PhpzgZeu7cuUSj0e6WODjTjB589UGSL6cILglz5rx5JGIxrvy//+NvQzp5e6gS/9o4NXPewN69P8W20xjGyXct5Frjheb3Q1eXxRNPOMf7586dy6JFi/rsubC1Zk+2CzzXqt6YSNCZybA+HmdKOMyMaJTpkQjXDR/O9EiEyeHwUefuF5pbIa6U4v+NH48f+MjGjXx93DjmRqPdYyiO9X8q+mfC1GmggLzTCRhKMWviRNdqEv3XrxBXSv0bcKfWeqBWhT4b2KK13pZ9/ruAa4CSDHFb267Pkc6FwLZEgnq/n3OyLe98bW0wblzhzh1++PDhU55jvWjRIs4++yyWLXuJeDxOMBgklUyRzCQBSMbjvPriiyTTaZ75yU9oiUb51dSpfGPK1ahUFaZZR1fXXiKR8Sddg2EYWEUYxKOUxVvfupAVK5zj/dFolPlnn80P7rmHLckkGxMJXst2f29KJKjx+ZgaDjM9EmFaJMKb6+uZGg4XdSR4f0X8EVoS7iwKr5Ti3ydM4O5Dh7hq5UqqbrqJA6tWdf+fLliwgKVLl0qQnyTb7uO1r7XME/eI/rbEh+O0hl8Ffgcs1cU4DdQRo4Dded/vARYU8flOiUa72p3enE6zMR4nYhjMr6yk4ihdqfG4s6Z3odi2TfgU1xY1TZMlS/7Jgw8+yNq1G1j78lr+ct9fetwnmU4D8P2f/hQAw+dj5bzzWPKPB/H7q7DtUztWW6ylYh9/fAmvvrqMWMzpZejs7OTJF17gXb/7HfMXLmR6JMLiIUOYNno0UyORfs/Nd4NbLfF87xo2jI5nnuGGV1/Fzr6QOzs7WbZsGUuWLGHxYun8PVGWZfH9m7/aoxUOzhkKb/ntb7nmC19wpzDRb/361NBa36yU+jpwBfAB4OdKqbuB32qttxazwKNRSt0A3AAUZf3rE+FWd3rSttkcj9OayTAtEjnuwJ943DlBR6GkUqmCnNLTNA2uuupy3vKWt3PPr+/h/ofuJ546emDYmQyvrnmFhx9/mHPeYKNUaR4nXbt2BfF4zx0M1dXFB+Nxbj7tNJeqOjlRf9T13iaAfevXo3udDz4Wi7Fy5UoJ8ZOwZMkS0tmd5HwK+NgpTh0VR5eyLBJak7As4rZNwrZJnGRvYL+TJ9vyPpC9ZIBa4P+UUt8/qWc+tr3AmLzvR2e35ddzu9Z6vtZ6/tChJz86uRDc6E4/nEqxvL2doGFwbnV1v0buplLOQkyFopQiU5D1lQ0sywntxdct5oxhZxCNOqEROlqvQiLOqtdWkk63EAgUf5FynV1jWmuNnb0cz+zZc4lEeu41RaNR5syZU6QqiydgBkpiQaPcGIp80WiU2bNnu1SRhynFipUryfQR4sNqavjAiZzmUPSQtm06MhkaUyl2dXWxKR5nVWcny9rbeaq1lZc6OtgYj9OYTmNrTa3Px6ST7NXs7zHxzwD/ChwGfgN8UWudVkoZwGbgSyf17Ef3MjBFKTUBJ7yvBd5T4OcomIHsTtdaszWR4FA6zekVFSc0CrnQc5br6uoKMsfaNP2k0xaZTBfBihB/+epfeHTNY6wbup3nH3qIVStXkrYsMnnH6CLhCKef0YDfH8HvrzrlGpwFZY4Ec++Izv9v6+/x6YULF3HmmQt6HBNfsGABixYtOuV6B1pbso1EeuDmiR/NokWLWLBgQfe6Avn/p7bW3X+nUug18IK5c+YQjkRI9HofTx09WuaIH0PatknaNl3ZVnT+167s51TIMAgZBuHspc7n697mK+CA1P4mQB3wNq31zvyNWmtbKVXwPiytdUYp9SlgKc4Us99prdcV+nkKxVRmQbqVj8fWmg3xOJbWnFVZecJLpBpGYU9MFAwGaWxsLMhj+XxVpFIHMM1xNLyngQu+eT6XjbyMi793Of9+05dZ89JL3fc1DIOzzz6Lc84JUV19Tr8eP9eKzu3F6Lztadsm6PNhc6RrqhADyUzT5O67l7J8+RJWrlzJnDlzPDuSuiPZQWXw1Je6PVWmabJ06VKWLOn7/1RrjQ1o28ZQzq61BPrRLbrySqadfhorX345b8/VIGlZWJY1KIPc1ppUXiAntXa+5r63bRRHQjqQ/VqTDelwgUP6ePp7TPwbx7htQ+HK6fG4DwIPFuOxC81QBpYu/hKFWxMJLNvm9IqKkwoZv9+Zt1wow4YNY+XKlQV5LJ+vinS6ha6uXYSqxzLupnHs/cVedt21ms2vrO4xBSwQCPDBD11EKNRAJDIVoHsnSucuvb5X2Ut3SGe/KsPApxQBnw+zwB/2WoPPZ7J48WLPH6/tSHVQGXA/xMEJ8qP9nyqlMDnyd89oDdkWuqGU66P8S4pSmH4/b//X97B6xSvY6dx7zGb9jh0sefFFFl9+uaslFlouoJNak7TtHteTtk06G9gBpQhmwzloGFSaJkP9foKGQVCpAQ3p4ynd4bAeYhomSStZ1OfoyGRoSqc5q6rqpD+IQiHnUihjxozh/vvvL9jjhUJjiMe30tm5jnDdRCZ+ayKHPneIxNM9u3GTySRr1u3gHdd+jaRtdwe2oRQGzgd57v/oeC2xXAu9GK1jrQs3nc9tHckOxla7O4D0RKhsK9zInqXLxhlxncq20A2cHblCrtLnSUqxbeNm7EyvxV66uli5caNnFnvRWpPO/n1Tvb/aNl3ZoM5oTdAwukM6F9TVpkkgG9ABw/BU742EeAGEfKGiHy9sy2QY5vefUmuxshJ27ChcTaNGjaK1tZXW1lZqTvEkFc4HrSIQnkRX1z6a21eBMph+znAikTCx2JFjdpFImHPPfTt+098d0if7prNtu2hTzGy7sAMJ3dSZ6iyJ7vSTkWudm0rhh+6BiWmtsSzLCfNcsA/C1vqUGdMwAyZW8khvYjQUYs6ECS5W5bCyh7tS2b9X7+tded/7lCKQDeHc16hpUuv3d4d2oAx32iTECyDsCxf95BAhw6A1kzmlNb6rqqCiwlm5rRDhopTizDPPZO3atVxwwQV93kdr3d29beN0cdq5bRzp+la5VjQQCDYQDDaQSbdy1ZV1zJ//B5YvX0M8nugeyHT1VVcVpPvbsqyihngJT/0+IaXUnX6qckGd+9PkQj2TN/PAyGvJ516X5insLJaysy8+j5pJtcS3x0gkEvj8Uc6eOZ1FRRjxb2X/nzPZrut09vt0XjBnsttyA8QChoFfKfzZYPYrRdQwCPh83dsCZfq36Y8y+YhxV9gfpqWruKtZDfH76bAs1sTjjAgEqDVNAieYxKbpzBVvajq19dNzwWxrzdQZM3j6+eeZf+652BwJaMAJ7ewHokH2QxBe/wF5lDdfwBwCoSE89thzRx3IdKps2y74muk5llXYwxdu6kiVxsC2Yugd6uC8Liyc13BGayytybVTc6Gu8lvv2W1ebMn7TJNzv3IB5+9dwN1P3U0o8UUe/H/jMBv393l/K/v/kdvxscj7P8ptywvnTN7FUApf9mICPqXwZ4M5ZBj4s8Gc217ocSrlSEK8AAaiO10pxYRwmM5MhsZ0moOpFKZShHMjJHMvepzh/IZhHBm8lfc49cNh936oqad7zlmuNaw5Ery5VnJuPnR+Kzo/gKfPmcMf/vIXDre0UF9bi5l93uMF9Ik41kCmU5HrSi/WHrxllU9LvDPVWTYt8f7If//k6z6+bttopbqDLJNbQ4AjO7LdO6u5x8x+T7bnKb/3CY68V/Jfjd3Xsz9zNLqPr/nv4949YHauVsuiVdukfFGmnXceUyv20b7uUnaaOzAAOx7vEdaW1ph5O+Vm9vcyc8GcbRnnDl/4el0Ga2u5mMrkI8ZdA9GdnlPh83Uvq5q0LBK23T1nMW7bWLlFSZSCvFZxbqBP1SRYvVcxdBr4s286hRP8uQFhSmtMw3hda7nPYPb7uWDePJ577DHe/e53D8j/QaEU+vSjvRXqsEUpGFkxkqrgqc/H97ru4+vH+cP2CPTcdY4EqIUzDa6vgNU4JzKy6DkVMn/omXY29nh/5+6X+773zkB+70HuNkM7q1yMrhyNqTVGOkWVaVCjfIQzGcxgsPv+vuyOjQRxaZEQL4CQL0RXpuv4dyywoGkSPMGUmDMCNm6EsQXs5r3yyiv52te+xuLFi1+3mlYpMwyjaCGeTjtT+srh8y5tpXl297M0VDa4XYpn5A+me51S2rOzbYb4oK1zB9WVELSSJDr9VBkmoXLaCy1j5TdUzwVhXxhbH5nqVMomTYJduwr7mMOHD+fiiy/mH//4R2EfuMAsy+KBBx7glltu4f777y/K6UdzciFeDna07qChooGAWZpr1ItTY9sW+17dxx/+8AcatzQSjVjYynS6kkTJk5Z4AfhMHz7DRywVoyJY4XY5xzRuHOzd65yWtLq6cI979dVXc+ONN3LRRRcxZsyY4//AALMsi4ULF75uuc5incIyk4F+LGfvCdtatjGpdpLbZYgisCyLz1//SVa+uoIXU89i+k3q6t7Gj9//cwlxj5CWeIFE/BGaEk1ul3FchgFnnw0bCrzOXk1NDe973/u48847C3RSlMJasmQJy5Yto7OzE611j1NYFkMqVT4hvrVlK5PqJMTL0ZIlS1i/ai1WMoPWmkwqQ0vTKyx9/lk5n7hHSIgXSH2k3hMhDk6X+ssvF/5x3/CGN2DbNv/zP/9T+Ac/RStWrCAW63la0NwpLAstk3F2lsrlcOLW5q3SEi9TK1aupCvRczxPOh1n9eZN0hL3CAnxAhkSHsLh+GG3y+iXefNg9+7Cv0eVUnzmM59h9+7dPPLII4V98FM0Z86c1y3qUqzTgqZSUKSp567Y1rpNWuJlau6cOYTCPUe5+v0RTp92mrTEPUJCvEDqI/U0xb3REh8xwgnwNWsK/9gVFRV89KMf5X//93955plnCv8EJ0FrzbJly6irq+s+T3lFRUXRTgtaTl3ptrbZ3rKdibUT3S5FFMGiRYuYOed0jKDTbeQL+KgbMp83XvgGaYl7hAxsK5DacC0dqQ4ydgafUfr/rRdfDM89B0VoiDJq1Ci+9a1v8fOf/5x4PM7ChQsL/yT9ZNs2X//619m6dSurV6/mpZdeKuppQXOfe+UyMn1/x36qglVUBEp7wKY4Oabfz0/+8DM+9Z8fp2NpG1OumEJ98i8oo7V8Vioqc9ISLxBDGdSEamhONLtdSr9cfDE8+yz0OkxcMBMnTuRzn/scf/rTn/jDH/5Q1OlcR9Pe3s4HPvABmpqauO222xg2bBiLFy/m5ptvZvHixUUZlZ5MlldX+tbmrZw+7HS3yxBFkJty+b+/+B+CZpj6+nrGnj4W2zadAR2dnW6XKPpBQryAvHRcvLraaYU//njxnmP06NH89Kc/Zc+ePdx4443s39/3WszF8OKLL3L99dczbdo0fvrTn57yWdb6Q2tnfng5hfjKgyslxMtQbsrlddddx29+fDuv/HoZGzduRCsNKCfEpTvdEyTEC2h4dDgHOw+6XUa/XXMN3Htvcd+r1dXVfOUrX+HCCy/kpptu4s477ySVShXt+Q4fPsxXv/pVfvCDH/DhD3+Ym266icAAHaBOp50eyHI62+GyvcuY1zDP7TJEgfWecmklLTo6Oti/YT+mCVoZMrDNI8ro48Z9wyuG0xhvdLuMfps+HerqnG71YlJK8Y53vINbb72VjRs38oUvfIF//OMfJBKFW2/+4MGD3HrrrVx99dUMGTKE3/3ud7zpTW8q2OP3RypVPmctA+hIdrCpaRNzRsxxuxRRYH1NubRtm+Y9zaTToHymhLhHyMiFAqqP1NPW1UZH0junbbz2WvjNb+CCC4o/r3nkyJHccsstrF+/nj/96U8sWbKE008/nXPOOYc5c+bgP8HRYLFYjGeeeYaHH36YBx98kPe+973cfffdrqwYl06X19xwgOX7ljNn+ByCvjI6PiAAmDt3LtFolM68496GYVA3to7KlIGNId3pHiEhXkBKKcbVjGNH6w5mDZ/ldjn9MncuRKPw2GNwxRUD85wzZ87k29/+NocOHeKJJ57gd7/7HZs3b2b69OnMmzePhoYG6uvrqa2t7e4KTyQSNDU1sWvXLjZv3szevXtpamqivr6eCy+8kG9+85tUVblzli2tnVZ4OOzK0xfNC3te4JzR57hdhiiCRYsWsWDBApY9+SQx28IM+AgHwow9YyxNjysMv7TEvUJCvMDGVY9jc/Nmz4S4UvCRj8APfwgXXjiwQTRs2DDe/e538+53v5vW1lZWr17N9u3bWbVqFbt372bXrl0kk0nGjBlDMBjEtm2qqqqYPn06V1xxBbNmzSJcAsmZTjst8HI6Fg5OiH/v8u+5XYYoAtM0Wbp0KUvOP5/HGmp5dVqSq+vfRNeoLnStcl7MEuKeICFeYONqxvHEjiewbAvT8Ebf6tSpMGUK3HmnE+huqKmp4aKLLuKiiy5yp4CTpLWzzGoJ7EsU1MHOgzQnmpleP93tUkSRmKbJ4nHjmHHpXL5et54vvPML3Lb8NlY2KpQpIe4VZdZ2cF/IF6IuXMe+jn1ul3JCPvxhePFF2LTJ7Uq8JZUqn/OG53txz4ssGLUAQ8lHRFkLBFDpdPe3WmsqKnBCXI6Je4K8Q4tgfM14drTucLuME1JTA+99L3z/+04wieOzLCe8y2V1tnxrDq3hvDHnuV2GKLZAAJVKUxuq6d7U0izHxL1EQrwIxteM91xLHOCSS2DWLLj9drcrKX25bvRyDPBkJsmW5i1cMPYCt0sRxRYMYqYzhP0RtIagGaS2TlriXiIhXgTDosOIpWOeWb0tRyn40Ifg+eeLP3fc6zIZZ2GXcutGB3hixxMEzAAjKka4XYootkAAUmka484iVbF0jLZWhemXY+JeISFeJNOHTGfTYe8dYK6ogG99C/7zP2HrVrerKU2WVX5zwvPd+9q9XDPtGrfLEAMhEMBIpSC7Mxr2h2lvV/j8yulu0trd+sRxSYgXyfSh01l9aDW29t7e7JQp8OlPwze+AU3eOLvqgMk1Tso1wA92HmRt41oum3iZ26WIgVBfj+3z5TKcjmQHhlIoQ8GQIdKl7gES4kVSF66jLlzH1mZvNmcvvRSuugq+9KXinenMa3INk3INcID7N93PFROvIOQro/VjxdEFg2BrNE6L29aaVNJwXuMdHdKl7gES4kU0a9gs1hxa43YZJ+3aa+Gss+DmmwdfkOdO03jLLbfwwAMPYFlW2Qe41pp7N97LNdOlK33QSKcxDzexZ/lubrnlFtY+s56JEy1nrEcmI93pHiCLvRTRlCFTeGLHE7Qn26kKurMk6KlQCj76UfjJT+Bzn3NWdauocLuq4sudpnHZsmXEYjGi0SgLFixg6dKlQPmm+MoDKzGVyaxh3lhtUJw6y+/nun8s5YXmJp5IPoE/5Cca2oFlLcWU05F6grTEi8hn+JhRP4M1B73bGlcKPvtZmDcPPv95OOytAfcnpfdpGjs7O1m2bBlLlixxu7SiunejM6BNleOQe9GnJZs28eqhJjJdGbTWpBIpOjqyr3XDkJa4B0iIF9ms4bPY0LjBkwPccnLrq19wgTNy/bXX3K6ouPo6TWMsFmPlypXuFDQAOlOd7GzbyeKpi90uRQygFfv3E89kemzLZLKvdUPminuBhHiR1UfqqQxVsvHwRrdLOSVKwfveB297G9x4I5Rzo3Ts2LkoFe2xLRqNMmfOHHcKGgD3bryXGfUzGBod6nYpYgDNnTSJSK+BHqaZfa1LS9wTJMQHwPyG+by892W3yyiIN7wBfvELJ8T/4z8gkXC7osJ66im4+eZFjB+/gIqKCpRSVFRUsGDBAhYtWuR2eUURT8f5/crf8/YZb3e7FDHAFp19NvNqqzECPkDhDwUYOTL7WpeWuCfIwLYBMLF2Ii/seYGtzVuZVDfJ7XJO2cSJ8J3vwM9/Dv/6r840tLPOcruqU9PRAT/9KfzP/8Af/mBy4YVLWbJkCStXrmTOnDksWrTIGehThv667q/Ma5hXFq9NcWLMcJi/nzmLyyf5sF64gBlXH2Bk5y+c17qcjtQTpCU+AJRSLBi1gGd3PYsuk+6paBS+/GVnUZhvftNplbe3u13VidMa/v53mDsX9uyBl15y1pA3TZPFixdz8803s3jx4rIN8Hg6zp1r7uQjZ7p0DlrhrmAQXybD6PmjCYS+wpgzZtLZmX2tS4h7goT4AJkyZAoAm5s3u1xJYV14Idx9t3MWtOuugz/8Abq63K6qf9asgRtugG9/G+64A/77v6Guzu2qBpa0wge5YBCVzhAIaoYOs8mkFSNHZm+TEPcECfEBdOG4C3lm5zNl0xrPiUbhE5+An/0MVq2Cq6+GP/8Zkkm3K+vbK6/Au9/tnHr1zDPhhRfg/PPdrmrgSStc5NZOD4Vg/wFNoksd2QmXEPcECfEBNLF2IkFfkA2HN7hdSlFMnOicOCU/zH/5S9i71+3KIJ2GRx5xdjbe+1647DJ48UX4+MedEzkNRtIKFwSDqGQaw9AEApp4p8Ho0dnbJMQ9wZUQV0q9Uym1TillK6Xm97rtq0qpLUqpjUqphW7UV0wXjr2QF3e/6Ol548czfTp897tw223OcfLrroMPfxgeemhgl2/V2ml1f/nLMHOmM2jtkkucHYwbboBweOBqKTWxVIwndzwprfDBLhjESKUYEh7CsOE2jQfzxjqbpoS4B7g1On0t8DbgV/kblVIzgWuB04CRwKNKqala67KZ5zCuZhwVwQpe2fcKZ43y+JDu45g4Eb7yFWfJ1ieecM5RfvvtztKt55/vXGbMKOx65Pv3O+dDf+452LcPdu+GN78ZHn4YJkwo3PN43R2r7mB8zXhphQ92wSDpofU0dzUzZapmWWeGUaOyt1VWSoh7gCshrrXeAPS1vOM1wF1a6ySwXSm1BTgbeGFgKyyuyyZcxp2r7+S0YacR8UfcLqfoAgFYuNC5xOOwfLkTsl/+MlRVOS3myZOdU6BOmgQjR0J1tXOb39/zsbR2Bs61tTkBvW0bbN/ujCw/eNBpZZ93nrOD8JGPOL0CsopoT9tatvHg5gf5/Vt+73Ypwm3BIJgGfsPP0KGa5hUGI0Zkb+vokHniHlBq88RHAS/mfb8nu+11lFI3ADcAjB07tviVFdCQyBDmNszlsW2PcfW0q90uZ0BFInDRRc4FoKUFNm2CLVtg82ZYtszZtnGj0xXv88Hs2c5Sr4mEE+CBgBPUTU1O63rCBGcRmpkzYdo051Ce6Jutbf7jmf/gg3M/SH2k3u1yhNtCIUKvbSFlTeD8C2zuXh1k+vTsbXJM3BOKFuJKqUeBEX3c9DWt9b2n+vha69uB2wHmz5/vueHe54w+h1+/8mu2tWxjYu1Et8txTW0tLFjgXHrT2gnuzk7n8yQcdi4S0ifvbxv+hqEM3jL9LW6XIkpBKIRKpQCYeZrNRz+ZorY2e5tSEuIeULQQ11pffhI/thcYk/f96Oy2shMwA7xpypt4cPODfGTeRwiYg3SI9DEo5bTcI+V/xGFAHIod4tev/prbrroNQ8mekCA7Oj2FRqO1xsg/9iRrp3tCqb2T7wOuVUoFlVITgCnASy7XVDQTaicwtnosT+540u1SRJnTWvOD537A22e8nQm1MsJPZPl8KNuGjIVG99y5k+50T3BritlblVJ7gHOBfyqllgJordcBdwPrgYeAT5bTyPS+XD7xcjY0bmBve1l2OIgS8eSOJ9nVvovr51zvdimilCiFDgYw0xlsbaPo1RKXEC95roS41vrvWuvRWuug1nq41nph3m23aq0naa2naa3L+ISXjrA/zMLJC3l468OkrbTb5Ygy1NrVyj0b7uEr539FDtuI19H+AGOCw7Btm7Avb/EEOSbuCaXWnT4oTa+fTl24joe3Pux2KaLMaK357rPfZXr9dOY2zHW7HFGCdChIU+s+NJouK+/EB3JM3BMkxEvElZOvZEfrDjY0lueSrMIdd629i7SV5qPzPup2KaJE6UCAWh1Ga42p8lZeku50T5AQLxFBX5C3zXgb/9z8Tw7HD7tdjigDr+5/lT+t/RNfPP+L+E3/8X9ADE7BAF3xNiztDG7rJiHuCRLiJaShsoHLJ17OX9b+hWSmRE8BJjzhUOwQ33zym/z7Rf/OiIq+lmsQwqFDwe6BbT1Gp8sxcU+QEC8xc0bMYULtBP7x2j/K7pSlYmCkrBQ3PXYT7z7t3WW/Pr84dToQxJ+yXh/ickzcEyTES9DCSQuJp+M8s+sZt0sRHvTjF37M8Ohw3jPrPW6XIjzAmjgOrKOEuLTES56EeAkyDZN3nvZOXjv8GusPrXe7HOEh/9z0T1YfXM3XLvpaXycYEuL1lEml7UdrTdAXPLI9FJKWuAdIiJeoikAFV0+9mgc2P8Dutt1ulyM8YNmeZdz+6u189/LvDoqz44nCMGIxdCKGpa2ea1Ukk5DJuFeY6BcJ8RLWUNnAW6e/lb+s+wuNsUa3yxElbH3jen6/8vfceumtjKkec/wfECInFGKoWdX3wDZpiZc8CfESN2XIFN448Y3876r/pSXR4nY5ogRtatrETY/dxPvnvJ/Th53udjnCa0IhOtoPS4h7lIS4B8weMZsLx13IHavuoK2rze1yRAnZ1rKNLz3yJT5/7uc5Z/Q5bpcjvCgYxJ86ytrpEuIlT0LcI84edTZnjzqb36/8PR3JDrfLESVgZ+tOvvDwF/j0gk9z4bgL3S5HeFU41D3FzDTyVmyTeeKeICHuIeeNOY+5DXMlyAV72/fyuYc/x8fmf4xLJ1zqdjnCy4IhfBmZJ+5VEuIec9G4i5jbMJffrvitdK0PUrvadvHDF37IB+Z8gCsmXeF2OcLjVDjMUFXhdKfnT0uUY+KeICHuQReMvYAFoxbw61d/LeusDzKbmjZx49IbuWLSFSyeutjtckQ5CIWIdTaTsTM9j4lLd7on+NwuQJycc8ecS8gX4u51d7N46mLGVo91uyRRZC/ueZE7V93JjefcyAVjL3C7HFEmVCjMEBVFKUXQzFvsRVriniAh7mFzG+YS8Ue4c/WdXDXlKmaPmO12SaIItNbcs+Ee7l53N9+85JvMHDrT7ZJEGVHhMImOFkJWCktbR26QZVc9QULc46bVT+NDcz/E/676Xw7HD3PphEtluc0ykrEz/HTZT1nfuJ6fLfoZwyuGu12SKDMqGMKftrBsq+f5xKUl7glyTLwMDK8YzsfP+jibmzfz13V/7bl0ovCsjmQHX37kyzTFm/jpop9KgIuiUNU1dIRNWezFoyTEy0RFoIIPzf0QAV+AX778Sxnw5nHbWrbx7We+zZQhU7jl0ltkLXRRNEYwRK2vAsu2eoa43+8EuShpEuJlxG/6uWbaNZwz+hxuW34baw6ucbskcYK01izZvIQvPfIl3jjhjXxs/sd6frAKUWBGMEhk7yHSdrrnYi+ZDKSlV6/UyTHxMqOUYsHoBYyuGs19G+9jU9MmFk9d3PMUg6IktXW18bOXfkZnqpMfL/yxnMhEDIxQiFodotXWBIxAz9tkYFvJk138MjWqahQfnPtBNJqfvPgTtrdsd7skcQwv7H6Bj/3zYwyLDuNbl3xLAlwMnHAYqytOPBMnf5o4hsSDF0hLvIwFfUHeMfMdrG9czx/X/JEzG87kjRPfiN/0u12ayIqlYs6hj0Nr+NqFX5OzkImBFwoRTNmkdfr1h26kJV7yZFdrEJg5dCY3nnMj7cl2fvLiT9javNXtkgSwfN9yPvbAx/Cbfv77qv8e3AEuo6DdEwoRTNtorfEZee06GdTmCdISHySigSjXnn4t6w6t4661dzGpbhKLpy6mIlDhdmmDTmOskV+98iu2tWzjC+d9QRbpAbjnHrcrGLzCYQIpm3g6LlPMPEhCfJA5bdhpTK6bzCPbHuFXr/yKeQ3zuGDsBT33wEVRJNIJ/r7h7zy+43HeMP4NfOn8LxEwA8f/wXLX1AQ33wwVskPpilCIkAUGxutb4hLiJa88PrlTKbcr8JSgL8jiqYs5q/Ms/rn5nzy36zneNOVNzBkxR1Z7K4KMnWHplqXcte4uzhh2Bre84RYaKhvcLqt03HwzvOMd8NBDblcyOIVCmMk07cl2WbHNg8ojxJuanAEYMpryhAyvGM4H536Qrc1buX/T/Tyz8xmumHQF0+qnSZgXgNaa53Y/xx0r72BYdBjfvPibTKqb5HZZpeXxx+HFF+G55yTE3eL3Y2jIpLp6zhOXEPeE8ghxpeDRR+EKObfyyZhUN4nPLPgMqw+u5h8b/0Foa4grJl3BjPoZEuYnQWvNC3te4G/r/wYKPn7Wxzmz4Uy3yyo9sRh8/vPwwx9CRFakc41S2MEAmUSnHFbzoPL4iw0ZAv/8J5x5JtTXu12NJymlmD1iNrOGz2L1wdU8sOkBntz+JGePPps5I+bIm7sfUlaK53Y9x93r7ibij/Cu097F2aPOlh2ho/nFL+Ccc+DSS92uZNDbP6qKQJecAMWLyuOT2e+HxYvhz3+GT3wCTPP4PyP6ZCiDOSPmMHv4bNY3ruepnU9x72v3ct6Y8zh/7PlUBavcLrHkHOg8wENbHuKx7Y8xs34mHz/r48waNkvC+1iefBL+8Q944AG3KxEAfj/xjuaea0j4/RCQgZelrjxCHOCyy2DNGifI3/tet6vxPKUUpw07jdOGncaBzgM8vfNp/uOZ/+D0YadzZsOZTK+fPqjX9LZsixUHVvDg5gfZ1LSJyyZcxg/e+ANGVIxwu7TSt2sXfPKT8LvfQV2d29UIIBRPEczY+I28EE+nIR53ryjRL+UT4oYBH/sYfPObziCZ8893u6KyMaJiBO867V0snrKYV/a/woObH+SPq//I2aPOZv7I+YysHDkoWp1aa7Y0b+GpnU/x7K5nmTZkGuePOZ+vXPAVmSrWX4kEXH893HgjLFjgdjUiJxwm0d4kr2MPKp8QB2dwzKc/Dd/9LowaBePHu11RWYkEIlw47kIuHHchBzoP8OLuF/ntit9i2ZbTBT9iNhNqJpRVoFu2xfrG9by872VWH1xNxs5w4dgL+c5l35FpYidKa2cg28yZ8KEPuV2NyGOHQoTS9AxxOSbuCeUV4gCjR8O//iv87GdOq7yy0u2KytKIihG8ZcZbuGb6Nexp38Oqg6v485o/k7EzjK8Zz/T66Uyvn051qNrtUk9YU7yJ9Y3rWb5vOSsOrGBk5UjOGnkWN55zI2Orx5bVTsqA+s1vYONG5zi4/B+WlFTAgK5Ez+504QnlF+IAZ58NO3bAL38JX/iCDHQrIqUUY6rHMKZ6DIunLuZw7DDrGtex6uAq/rr+r1T4K5gyZArja8YzvmY8DZUNJXUs3bIt9rTvYXvrdjY0bmBd4zri6ThzR8xl1vBZXD/nemrDtW6X6X0vvAD/9V/w4IMQDrtdjehFh0IYiSQ+szwjoZyV71/sHe9w5p/+9a9w7bVuVzNo1EfruTh6MRePvxhb2+xp28OWli1sbNrIQ1seoi3ZxqTaSQyNDqWhooGGygZGVIygOlhd1BaurW2aE83s79jPodghtrVsY1vLNvZ07GFoZChT6qYwoXYCi6cuZnTVaGltF9K+ffDRj8LPfw5jx7pdjeiDHQriT2V6bpT3gCeUb4gbBnz8406X+vjxznxUMaAMZTC2Zixja458cMfTcXa17WJP+x72tO9h+b7l7O/cD0DIF6I2VEttuJbqYDWVgUoigQhhX5iwP0zQDGIoA6UUCoVSCsu2SNtpkpkkSStJykrR2tVKa1crbV1ttHa1olC81vQalYFKhlcMZ1z1OMZUj+Hi8RczrnocQV/Qrf+i8pdKOce/P/xhuOQSt6sRR6FDISLpRB83yDHxUudKiCulfgBcDaSArcAHtNat2du+CnwIsIBPa62XnvQTVVQ4A92+/30YOVJaASUg4o90Hy/P15XpoiXRQmtXKy1dLbQmnK97O/aSyCRIpJ3jdW3JNjQarTUaTUWggmQmSdAXJGgGqQvXYRomdeE6JtRMoCZUQ3WomqGRoRLWbrjpJhgxAv7t39yuRBzD6BFTuaCujziQEC95brXEHwG+qrXOKKW+B3wV+LJSaiZwLXAaMBJ4VCk1VWttnfQzjR0L//Iv8NOfwre+BdFoIeoXBRbyhWiobJAR3+Xkj3901kVfskS6ZktcTV0DNb3PIid/M09wZYSR1vphrXXuAMyLwOjs9WuAu7TWSa31dmALcPYpP+G558K8ec5AN9s+5YcTQhzHq6/Ct78N//M/MkPEC8JhZw6/8JxSGCb8QWBJ9vooYHfebXuy207du94FlgV//3tBHk4IcRQ7dzoLL/385zBlitvViP6orn59A8fnc5ZeFSWtaCGulHpUKbW2j8s1eff5GpAB/ngSj3+DUmq5Ump5Y2Pj8X/ANOFTn4KXX5YgF6JY9u93ZoN8/vPOUsjCG0Ih6OrquS2dfv02UXKKdkxca335sW5XSl0PLAYu07p79MReYEze3UZnt/X1+LcDtwPMnz+/f6MvKirgq1+FW245ctIUIURhHDrk9Hj967/Cu9/tdjXiRASD0Nrac5sMavMEV7rTlVJXAl8C3qy1zl9h/z7gWqVUUCk1AZgCvFTQJ6+uhq99DR5/HB56qKAPLcSg1dzstMDf/nZnaqfwlnC475OdyOC2kufWMfGfA5XAI0qplUqp2wC01uuAu4H1wEPAJ09pZPrR1NY6Qf7QQ/DIIwV/eCEGlbY2J8AXLoTPfMbtasTJiETkjGUe5coUM6315GPcditwa9GLGDLECfJbbnEGcLzhDUV/SiHKTkcHvOc9zlkDv/Qlabl5lYS4Z5XC6HT3DB3qBPk998Azz7hdjRDeEovB+94Hs2fDv/+7BLiXHS3E5W9a8gZ3iAMMH+6sKnXXXc5JGoQQx5dIwPvfD5MnO/PB5cPe2/oKcRnY5gkS4uAsyfrVr8If/gAvFXYcnRBlJ5l01kNvaIDvfc85T4HwNulO9yx59+WMHg1f/rKzwtSrr7pdjRClKZ12zkhWUQE//rGc5rdcSEvcsyTE840bB1/8Ivz617BihdvVCFFaUimnxyoQcJYw9pXvSRAHnUgEpk59/XY5TFLyJMR7mzgRbrwR7r4bnnrK7WqEKA0tLc5SqlrDL34hAV5u/H7nb5pKHdlmGLLsqgdIiPdl6lRnidZ77oE//Um6lcTgtm0bvOUtzjroP/iBfLCXqy1benapp9M9Q12UJAnxoxk1Cm69FTZtgh/+UNYQFoPT88/DO94Bn/iE05Uug9jKVzTaM8S1lu50D5B35LFUVsLNNzsv7m98A5qa3K5IiIFz111Oj9QvfiFroQ8GkQh0dh75XmvZafMA+Qsdj8/nHAs8/3xnYZgtW9yuSIjisixn7vdtt8H//R+ce67bFYmBEI06C/jk2La0xD1ARqf0h1Lw5jc788m/+11njqx8sIlyFIs56593dMA//gE1NW5XJAaKdKd7krTET8T8+U73+p13Oi0UGfAmysm+fc5ZyGprnde4BPjg0rslLiHuCRLiJ2r8eGfA24oV8LOfyehNUR5WrXJGoL/tbfD978sI9MGo94IvEuKeICF+MmpqnIFuWsO3vgWtrW5XJMTJe+ghuP5654x+N9wgH9yDVUWFcxglR0LcEyTET1YgAJ/+NMydC//1X7BundsVCXFiEgn4z/+E73zHWQ9h4UK3KxJuGjoUMpkj3/v9EAy6V4/oFwnxU6GUM4d28WL4yU+cddele114wapVcPXVsHcv3H8/zJjhdkXCbT5fz17FWMxZ8EWUNAnxQpg3D370I+cN8PnPw+bNblckRN/SaefEJTfcAJ/7nHO9qsrtqkQpqKjoOU/ctmWeuAfIFLNCqax01lx//nmne/KNb4R3vlPWmBalY/NmZydz6FCn9T1smNsViVIiIe5J8hcqtPPOc1rlO3Y4pzbdudPtisRgZ9vw29/CddfBe94Dv/mNBLh4PQlxT1K6DOY6K6UagWKmZT1wuIiPPxC8/jtI/e4r9u9wJvBqER+/2ORv7D6v1z9Na115Ij9QFn29WuuhxXx8pdRyrfX8Yj5HsXn9d5D63Vfs38Hr/0derx+8/zuUQ/0n+jPSVyKEEEJ4lIS4EEII4VES4v1zu9sFFIDXfwep333F/h28/n/k9frB+7/DoKu/LAa2CSGEEIORtMSFEEIIj5IQPwal1A+UUq8ppVYrpf6ulKrJu+2rSqktSqmNSqmSXHRaKfVOpdQ6pZStlJrf67aSrz9HKXVlts4tSqmvuF3P8SilfqeUOqSUWpu3rU4p9YhSanP2a62bNR6LUmqMUuoJpdT67OvnM9ntRfkd5H3mPq+9x0DeZ9201nI5ygW4AvBlr38P+F72+kxgFRAEJgBbAdPtevuofwYwDXgSmJ+33RP1Z2s1s/VNBALZume6Xddxar4IZ87z2rxt3we+kr3+ldxrqRQvQANwZvZ6JbAp+5opyu8g7zPX6/fceyxbt7zPtJaW+LForR/WWudO6/MiMDp7/RrgLq11Umu9HdgCnO1Gjceitd6gtd7Yx02eqD/rbGCL1nqb1joF3IVTf8nSWj8NNPfafA1wR/b6HcBbBrKmE6G13q+1fjV7vQPYAIyiSL+DvM9c57n3GMj7LEdCvP8+CCzJXh8F7M67bU92m1d4qX4v1Xosw7XW+7PXDwDD3Symv5RS44G5wDIG5neQ99nA80qd/THo3mdlsWLbqVBKPQqM6OOmr2mt783e52tABvjjQNbWH/2pX5QWrbVWSpX8tBClVAVwD/BZrXW7Uqr7thP9HeR9JgbaYHmfDfoQ11pffqzblVLXA4uBy3T2IAWwFxiTd7fR2W0D7nj1H0XJ1N8PXqr1WA4qpRq01vuVUg3AIbcLOhallB/ng+WPWuu/ZTef9O8g7zOgdF+7XqmzPwbd+0y6049BKXUl8CXgzVrreN5N9wHXKqWCSqkJwBTgJTdqPEleqv9lYIpSaoJSKgBci1O/19wHvD97/f1AybbelNMU+C2wQWv9o7ybivI7yPvMdeXyHoPB+D5ze4ReKV9wBqLsBlZmL7fl3fY1nBGdG4FFbtd6lPrfinN8KwkcBJZ6qf68Wt+EM3JzK073pes1HafePwP7gXT2//9DwBDgMWAz8ChQ53adx6j/AkADq/Ne+28q1u8g7zP3L157j2VrlveZ1rJimxBCCOFV0p0uhBBCeJSEuBBCCOFREuJCCCGER0mICyGEEB4lIS6EEEJ4lIS4EEII4VES4kIIIYRHSYiLkqWUOit7jumQUiqaPefu6W7XJUQ5UUp9TCm1MnvZrpR6wu2aRP/JYi+ipCmlvg2EgDCwR2v9HZdLEqIsZdfxfhz4vtb6frfrEf0jIS5KWnYt55eBLuA8rbXlcklClCWl1C+BRq31N9yuRfTfoD+LmSh5Q4AKwI/TIo+5W44Q5Sd7FrlxwKdcLkWcIGmJi5KmlLoPuAuYADRoreVDRogCUkrNA+4ALtRat7hdjzgx0hIXJUsp9a9AWmv9J6WUCTyvlLpUa/2427UJUUY+BdQBTzhnx2S51vrD7pYk+kta4kIIIYRHyRQzIYQQwqMkxIUQQgiPkhAXQgghPEpCXAghhPAoCXEhhBDCoyTEhRBCCI+SEBdCCCE8SkJcCCGE8Kj/D0Vn39im41VpAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 576x576 with 3 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig, ax_main, ax_top, ax_right = rebound.OrbitPlot(sim,slices=0.8,color=True)\n",
"ax_right.set_xlim([-25,25])\n",
"ax_top.set_ylim([-25,25])\n",
"ax_main.set_xlim([-25,25])\n",
"ax_main.set_ylim([-25,25]);"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "ff3edb4e",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.8.2"
}
},
"nbformat": 4,
"nbformat_minor": 5
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment