bridge_gravity_with_hydro.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "AMUSE tutorial on high-order bridge\n",
    "====================\n",
    "\n",
    "Hierarchical coupling strategies are fundamental parts of AMUSE.\n",
    "It enables us to combined the output of a wide variety of solvers into a homogeneous solution.\n",
    "In this example we will be nesting multiple bridges, to show the power of bridge."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [],
   "source": [
    "import numpy\n",
    "from amuse.units import (units, constants)\n",
    "#from amuse.lab import Particles"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "                 key         mass         name  semimajor_axis           vx           vy           vz            x            y            z\n",
      "                   -         MSun         none  149597870691.0 * m  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  149597870691.0 * m  149597870691.0 * m  149597870691.0 * m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.000e+00         star    0.000e+00    0.000e+00   -3.465e-09   -0.000e+00   -5.009e-03   -0.000e+00   -0.000e+00\n",
      " 5326756165488562113    9.546e-04       planet    0.000e+00    0.000e+00    3.581e-06    0.000e+00    5.195e+00    0.000e+00    0.000e+00\n",
      "10820997325257666469    9.546e-06         moon    3.550e-02    0.000e+00    4.927e-06    0.000e+00    5.230e+00    0.000e+00    0.000e+00\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n"
     ]
    }
   ],
   "source": [
    "def orbital_period(Mtot, a):\n",
    "    return (((4 * numpy.pi**2) * a**3)/(constants.G * Mtot)).sqrt()\n",
    "\n",
    "from amuse.ext.orbital_elements import new_binary_from_orbital_elements\n",
    "from amuse.ext.orbital_elements import orbital_elements_from_binary\n",
    "\n",
    "def get_star_planet_and_moon():\n",
    "    Mstar = 1.0|units.MSun\n",
    "    Mplanet = 1.0|units.MJupiter\n",
    "    a_planet = 5.2 | units.au\n",
    "    e_planet = 0.0\n",
    "    bodies = new_binary_from_orbital_elements(Mstar, Mplanet, a_planet, e_planet,\n",
    "                                             G=constants.G)\n",
    "    bodies[0].name = \"star\"\n",
    "    bodies[1].name = \"planet\"\n",
    "    planet = bodies[bodies.name==\"planet\"]\n",
    "    RH_planet = a_planet*(1.0-e_planet)*(Mplanet/(3*Mstar))**(1./3.)\n",
    "    a_moon = 0.1*RH_planet\n",
    "    e_moon = 0.0\n",
    "    Mmoon = 0.01*Mplanet\n",
    "    Pmoon = orbital_period(Mplanet+Mmoon, a_moon)\n",
    "    moon = new_binary_from_orbital_elements(planet.mass, \n",
    "                                            Mmoon, a_moon, e_moon,\n",
    "                                            G=constants.G)\n",
    "    moon.position -= moon[0].position\n",
    "    moon.velocity -= moon[0].velocity\n",
    "    moon = moon[1].as_set()\n",
    "    moon.position += planet.position\n",
    "    moon.velocity += planet.velocity\n",
    "    moon.semimajor_axis = a_moon\n",
    "    moon.name = \"moon\"\n",
    "    bodies.add_particle(moon)\n",
    "    bodies.move_to_center()\n",
    "    return bodies\n",
    "bodies = get_star_planet_and_moon()\n",
    "print(bodies)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we have the orbits of the three particles, star, planet and moon. We now want to make a disk around the moon."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "metadata": {},
   "outputs": [],
   "source": [
    "from amuse.ext.protodisk import ProtoPlanetaryDisk\n",
    "\n",
    "from amuse.lab import nbody_system\n",
    "\n",
    "\n",
    "def make_disk_around_celestial_body(moon, Mplanet):\n",
    "\n",
    "    R = 1|units.au\n",
    "    a_moon = moon.semimajor_axis\n",
    "    e_moon = 0.0\n",
    "    Mmoon = moon.mass.sum()\n",
    "    RH_moon = a_moon*(1.0-e_moon)*(Mmoon/(3*Mplanet))**(1./3.)\n",
    "    converter = nbody_system.nbody_to_si(Mmoon, R)\n",
    "    Ndisk = 1000\n",
    "    Rin = 0.03*RH_moon\n",
    "    Rout = 0.3*RH_moon\n",
    "    Pinner = orbital_period(Mmoon, Rin)\n",
    "    Mdisk = 0.01 * Mmoon\n",
    "\n",
    "    disk = ProtoPlanetaryDisk(Ndisk,\n",
    "                              convert_nbody=converter,\n",
    "                              Rmin=Rin/R,\n",
    "                              Rmax=Rout/R,\n",
    "                              q_out=10.0,\n",
    "                              discfraction=Mdisk/Mmoon).result\n",
    "    disk.name = \"disk \"\n",
    "    disk.move_to_center()\n",
    "    disk.position += moon.position\n",
    "    disk.velocity += moon.velocity\n",
    "\n",
    "    masses = Mdisk/float(Ndisk)\n",
    "    disk.mass = masses\n",
    "    rho = 3.0| (units.g/units.cm**3)\n",
    "    disk.radius = (disk.mass/(4*rho))**(1./3.)\n",
    "    return disk, Pinner, converter\n",
    "\n",
    "planet = bodies[bodies.name==\"planet\"]\n",
    "moon = bodies[bodies.name==\"moon\"]\n",
    "disk, Pinner, converter = make_disk_around_celestial_body(moon, \n",
    "                                               planet.mass.sum())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXsAAAEMCAYAAAAlGRZyAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO29X4xcRbon+IvKOjZZ7h2n3fJKTQ0GzIyw1vK4Csxt7/KA3Us3q+GCasFgteBxht5HfFu+18xFFxtxB3dbvfDazGuzyGB6Sma4En8WzENfwVygyo0sgfY23ZhNWlq0dnkHV9rOqop9yIysyMj4Ir6Ic05mVmX8JMuVmedPnDgRX3zx/fl9QkqJhISEhIT1jbFBNyAhISEhoXwkYZ+QkJAwAkjCPiEhIWEEkIR9QkJCwgggCfuEhISEEcD4oBtQFoQQNQA72v/elVIuDLhJCQkJCQPDetbs97b/fYmWwE9ISEgYWaxJYS+EOCiEeMfy3b1CiCfaX30M4E4A/wktgZ+QkJAwsliTwl5KeVr/LIQ42P7+3fbnewE8KqX8GYD/GcATPRdJSEhIGCEMpbAXQtwhhLhD+/xE2wZP4S6sau9fArgDwMdtob8DwGnqxISEhIRRwFA6aKWUn7bNMkDL7u5zsJoLwfellJ+W1sCEhISENYah1OyBjqnmZ+2/fTb3BQBbS29UQkJCwhrF0Ar7th3+1+2/fdE0/4RV7X4HgHccxyYkJCSMHIZS2Lft9QtSyk+llC8BuFe32bdt8Xs1x+xpADva39eUozYhISEhoQWRKI4TEhIS1j+GzkErhJgDsA3APw+6LQkJCQlrCP8KwLdSymnbj0Mn7AFs27x58+TU1NTkoBuSkJCQsFYwPz+Py5cvk78Po7D/56mpqcmzZ88Ouh0JCQkJawb79+/HBx98QFpEhtJBm5CQkJBQLJKwT0hISBgBJGGfkJCQMAJIwj4hISFhBJCEfUJCQsIIYBijcRISEgzMztVx8q0v8M1CAzfWqjhy3+2YmU7RyQl8JGGf0BeECKsk2LoxO1fHkdfOobnSynavLzRw5LVzADDS/ZIQhiTsE0rH7FwdT/32MzSaywBawuqp334GoFdYhRybpz1raTE5duZ8R9ArNFckjp05P9TtThguJJt9Quk4+dYXHeGt0Ggu4+RbX7CPffLUPO4+8R5m5+q52qIWk/pCAxKri0ne63Lue/eJ93Dr0TeDn2Oh0Qz6PiHBhqTZJ5SObxYa7O+pY4FitHzXwlOWllzmbmWt7VISBoek2SeUjhtrVfb31LEK1I6Ai5CFpyiE7Gxs2DKRWb/ftKEykF1KwtpEEvYJueEzURy573ZUs0rXd9WsgiP33d5zLduxJvII5pCFpyjkXWCeeWAXsoro+i6rCGSVsahFJI9JKWHtIgn7hFzg2MBnpifx/EO7MVmrQgCYrFXx/EO7reYG/VgKeQQzZ+EpWhjmXWBmpidx8uCerv47eXAPLhM2e9ciMiifRcLgkWz2CbnAtYHPTE+ybcnqWNPWDXQL5hh7tfqdOq8M+/qR+253PocO6pls/XfyrS9Qtwh21yIyCJ9FwnAgCfuEXCjTBu4SzHmEsmvhCRGG3MXGt8Do1wt5ppBFRIF6L/WFBu4+8V5y9K5jJGGfkAs31qrB2mUIKMFclobKXbxCBTNnZxP6TNxFRAf1vkT7GXzPwlngUoTQcGIgwr5dPHxH+99dUsq/GUQ7EvIjRrssAmXtKLiLl08wxwi8mGcKMY8B9vclAJiVqG2LDGeB60dSXEIcBuWgfRTAXinlaQAQQjwxoHYk5MTM9CQevnMSFdGKFqkIgYfvtAugIh2fZUXVuBy4evttCwLQEm7Tz76NJ0/NBztB+xEpZHOWm4JewVxkOCGkecNME8rDQDR7KeVL2scdAN7Rf5+fn8f+/fut56ZyhcOF2bk6Xv+kjmXZEhnLUuL1T+rYe/PWHq3wyOlzaC5r/C6n4/ldytpRUKYRAD33s0EAuLTYGyXDMTH1a5dk7gbuPvEeazfj2nmonYxrEZydqyftfoAYqM1eCLEDwEUp5buDbEdCPLh25uNvnO8IeoXmssTxN+L5XW7IVuPMa9UMxx7c1bPAxNiObaaRu0+8xxL0lJYM+EMiVV9WhMCylJjsk72bu8jUJjLrQjbRTu7y9U8y5wwWg3bQHpRS/sz8cmpqKmnwawRcO7NNSLi+d8EWknltacV5TF7bsUtQC9C2fh2UOcZs67KUHWHbD8HIdfRKYiVbvL7sXOQUUojnYDEwYS+EOCil/GX77zuklJ8Oqi0J8Sg7GscGzm6i6Ggd6jkna1X87uiPANDmEMBtjol19nJ3LrqJhdo12HYz5vUp4jWOoFcok5YiwY1BRePcC+AXQoin2l+laJwBI9bkwTYBVDOrsKhV7bwvLnB2E0VH6xzYuQ0vf3ihS7CZz2nrC8BuYtLhs3Pbdigff3URr39S9+5cbLsG2/Hm+z+wc1vP9SkzlVpAOChTCUhwY1AO2ncB3DaIeyf0Im+CEuA3ARx7cFdXAQ4AyMYEjj24K7i9nN1EkTsO5YTWxZkAeqKOYuLeAVpYVoQgtf5XPvq65xzbzsV2vnk8gJ73by5sQEvQmwK/mlXw8J2TXQsD0Hq3EOjy09hoKVLMfv8waJt9whCAEijHzpxnZ4gWtSjY4NM6Ab6WfeXaUldUCEeY2PpHAnj/82+j+sIEpRUvS0nuRKhzzON9O5lvFhrk89kg0TJdmf219+at1gimPLQUKWa/WCRhn0AKhIVGs2N68W37ubw0oZPUNuFf/6SOh++cxPuff+u0Y2+uZhgTwJXrq4JsodHsPAfQq9HahEnZtMiTDn+AapcJajdg7lx8juMb24I7pK3KR6GDerd5aCkSj0+xSMJ+DaLorS0nkgRwb/vL0rioCf/+59/2CB1zYVhoNNFNDNz9HFeuLbGECdU/Y0Lg1qNv5n4HPr+H7Teb6cTmK6F2OPrxVHy8zWRTVMz/IPwuo45EcbzGUAZFLYdDXoHa9peVJRky4UPMEfWFBhldYl6b6p9lKQt5By4KaOq352Z2s2ijTcpolemsH09lDT+2bzuLlhoIz47mZAsPovbAekbS7NcYytja2uzpi9eXrDHwrm2/jSws7w6E0qolWqGO+jU5uxMFVwSJKUzM/hmznJuXF8dl4ooxf3HPj0nm4kTu+HZ6nCiuQfEurVckYb/GUNbW1hQILi55Do96Uc41lxlCvyZAZ7DazBGubE+bMNH759ajb1rPU7QB/TBxuUIybb4MzjU4yVy2+9oid3wKCMdhn8epn9CLJOyHDD6tsF9JTL6JZguj1IVkUTsQvR2259bNR5Sgf2zf9h4BSF1vy0TmbZ/rHfTLqUjdRxe8voUmpq0hpjKfAsKN4krCvRgkYT9EsGlNR147h+NvnMfCYpMddlgUnBPN8Hw2VyQ+/upi5/gidyCqHbcefdMqWFzXlACem9lt/c22c3nmAX/cv8u8cPjUvPWc0Of2LfrU9UI0bF8hE5sWHfIceRQQM6pKCHTmQNLu45CE/QBhTmhbdEhzRXZs576ww37h5Ftf9JCaAcDLH17osF2WsQPxXZMKUbSxLeYxEbjODS0VaBPqgD/aiRtBBawKaPNem4msZuqeAE2GZiKPAmKLqvK1K8GPJOwHBJsWzwEVdthPuLRKpUX6nGsxTszQEEWgZYemhENs3L/e7hcOTXVdI7TerE2o62yeCqaGzi1CArQWBtu9sopANia6zHGuewI0GVo1G8PWTRsLUUBcWb9UuxL8SMJ+QPANaBcGHWfs0ipV21zab6wTk6ON//zVcywagRhw2h2yY6Bs5tS40N+77T63fL+Kf/zDRWtsvO1ezWWJLRMZJjaMe9+nwmViJ3C1uVKYAsIZ34OeA2sRSdgPCHkG66DjjJVtmtIiFSjNOY8T0xeiSNnMiyiozW23KYiVAzmP/Ruwh4Sa1AIUfw/VLwuLTcz93U/YBUz6ESAQQxedOHT8SElVAwI1ObZMZJ1Ello1Q1bp9oQOQ5zxzPQkHtu3vSc7ldu2MjMjqX5VBbXzJEGFFiP33Y9qq+29ZxXh7Fsff48vQclVjlEH97g8JSh9SX42QrWiEw3XI5KwHxCoSfPMA7vwu6M/wh9P3I/5Z36Ckwf3sLMY+4nnZnbjhUNTUW0rMzPS1q+ugtoh4Labm2FMjYG/3POD3gZ7GIR9C5FPSNsydR++s+Vw1gW2LSNXPZsSrnmF78x0d11jAWDThgoEWsrQxvExHD4132mTr7+LrH28lpHMOAMC17Y7zHHGsW3jOjFjt+Ybx1cdnFsc0SOhOwluu7k7AGoMnHzrix6naXNFOs1cPvMKN4lJfX569jNvzD7lv8iba2DWNZYAVmQrX8KWqevycSTmzFUkYT9AFBENspZsk2bs9A3ZGBk7HTNJbVm/V5srpMAP2UnYaAV0rVZvV4hd25a5zHWW6uAsRNzxNjtX92bFugR6XjNdKIe/iwE0MWeuIplx1hDWsm3SbPtCo4mrzRW8cGgKvzv6o56JF0O2Rp0jJVh2Zk7bgVUuebPqk3oPXLs2dR8KrsXJZoaJNfmdfOsLb1asS6DnNdOFcvgrmgcdqr8Tc+YqkrBfQ4hlmxwGm2Vo22MmKfXb5UYzlyDkhMk2msv4+avneuzaIfdz3YezWMxMT3b8PbYFlAsf/73+v+1322KXVQSuXFtijUHq2sqGb0L1r62/E3PmKpIZpySUYW6JEYDDYrMMbXtMiJ/rHJ8Jw/W+uFqgmcAV2r+u+yhn6eFT8x3ajLKyqF2MoGrBcZmNTP9AbSLDd1eXyEI4tnuEcvhT/Z2YM1eRNPsSUJa5JUZLidkNlLETCG17jCkkr/mEel8hWmAeXn9XKObrn9S72vebDy+UZs5zFQ/Xk8dcuxd9lzGxYbzH4ezqJ+raXA5/27W2TKwWtt84PppiL2n2JaAMp9DsXB1Xri31fO8TZqEadVk7gVANK4a7JpbvxrcgLl7v7XcXuDsBbm1dIcAyIxXldPSVSVTg7l5iCNeoa8dGgF1trnT+1ktTjpKTNgn7ElC0U8gWZQK0wgqfeWCXc8CGmkPKil6IFd5mpIovC7ZI84kvtI8CZycQUluXyn41UZTTsWjThysjth9mxRSR00IS9iWg6JRyynE3sWE8N3mYCddCxfFDuI6JEd76sWX5Hqj3pUIrQ8AViiG1dSkmTRNF1cSdmZ7Ex19d7IQ6VoTo0C7EwFWABihf8KaInBaSsC8BRWtGeQZrqEZNCb7N1YwUtur69YVGV7aqSyBzhLe+cPhKAeYB9b5CBX1FiC4bcozT1/a9T1gqmKGgQNxCaCY1LUuJ1z+pd+irQ6GPwZgcgrxwFYy30V+vV4ymp6JkFBnzDOSnFwgJyaOcnDa7caO5jL96dR5HXjvXmUxcWgJOirvuNKWchjFCwnRAA+hKz1earGmjVtgykVn76FeP7ukhJgt1+lKJV+Z4erxdDJwC11Fsc8bHhvi6oMYg1eYyQyFdBePXSp5KEUiafUkokuagjPAxSuukdgKU3XhFAiuO6A3ALpB9mi2XAjpUSFirgZ0+B0j0aLJUqJ+qZuXaLfnsxDEOa9/uyIRvIaR2V7HX42BQoZC2GgHAaNnuk7BfA4iNMqHgM6HYBAvXbmyDKp6ht5+qeCQBkm7XRDWr4MDObUHUxRSvuwllP3/+od1eoW6jMXY5JPXj8rxT34IoAUw/+zbpxKcWJBf9QF5Qzw0gNwW1Db4FERgd230S9pEom6PGdv2iikP4tum25+LajU0ogdxTJWlMIKsIq6A1bf869O8bzWX85sMLXef5bNUhE/ubhUZn4Zudq+PYmfN48tQ8njw1j4lsDM0V2Wm/eW9XYpLuRM3zTjnPcmmx2dq5gM+nvyxlz7vxUSyHwOaoL8v5ztkhjko2bRL2EeAUBi9a8y7y+txQQ9uks1WCMqEE8qTG4tijTa9I1KoZNm20V0mS6KUmzsYEIOyauIJvWx5Su1XtMg7s3IZT//XrrsSgRS1u23ZvVx/pNnwgXqBxn6W5bGfMpM7f0s547Wl0AEKUoTJDI30LogBwYOe2XPdYK0gO2gAoZ9aTp+atwuvSYrOQjEZKOBZ1fRf3iEvjn5mexK8e3WPlPalVs47z8IVDU/iT5gx2cda4NFtdvtSqGb53w7hT0Cv42CFdhTFMqGxVqk4rdW+X81Qhr9Mz5FmoKB+bo1lKkBTLHIRmkA+imI2CBPD6J/WRcNImYc+EyXzoQ56JzLlHnutTk5wT8aJHhgCtBaK5LLFp4zjJYLm5msEGZcu301t149rSCslLb7uuDXqkCeeeMXBVfrIhj0Azo3QoojAAqE30vgMqaoyqM5uXopgar/0uZmMi76K7VpDMOExwo0N0xIYFUvbqIq4PuItmUCYVPa1dnc+xs87O1XHFQjeQjYnOPTnP2mguY0y0on9coCI7TNNYoFWCBbPyE7Dax7Y8ASC/QNPfx+xcHUdeO2fdhXx3dckaUx7ijM9LUUx9H1rMpr7Q6PhEJj0mIvM9+Kib1zOSsGciZjDETGSu8Iu9vgIVGko5YU1hzrWznnzrC6vp5Xs3tLJ/uVQAgF/Qu+gjYhbrENjubQpim0ALjSZyQZ1nKwZvVrpy2dTzhkeGZpBzIpPM/gtJINPfA7ewug1ruXAQkIQ9GyGOPSA+dpi7qHDZHEO5aPQ0eRO6MHc5eXUNkjpuoW2SCe1XCi8emnI+G7dft0xkuNpcCV4YfNQVNoFmEp8VFYXi015tAQCHT83j468u4rmZ3Szh60LMYuHLSzl25nwhdAuxC9mwUIXnQRL2TKjEIo7WXatmOPagm6CMAiX8VOSKS/PJK0jMNHkblMBwCekjp8/h2JnzuNxoes0Xtslni7qpZhVsHB/rcKLrmGzz1bvAWVTMhClXCKgJLnWF3s67T7xXaBQKt9KVbZcjAbz84YUOJUKepEBqYdP5+EMWj9m5uvW96+Au5rELGbWTPXbmfBL2PgghDgL4mZTyx4NqQwhmpifxpMPkIIBCtnaU5uFaPGxah6+GqHk+N2nKJaQVmsuyMzltgt5l29aTbGzfxZoXfHkCJq+NzdyxuZqRQifGpFZ0FAq30hV1fYlWaC2QX1t1mbBCtWKO8zSk/2MWMnKH2miuGX6dgQl7KeVpIcTPBnX/GLh4votKeDJNKRzGQUpTs8E0s3AyDBVsQtq1AOqoCIEVKa0LIjX5XHHZoeYFX3up3Yzetuln37YeI4Aok13R7KiuRUJfyFy7HLPaVhGgtOInT83j56+ew09/eBOem9lNnu9b/PpBt+Dqs7VCtzCUZpz5+Xns37/f+tvZs2f72haF2bm6tYhF0QMthnEwVBP0OVptEEBn0dG1XVemqI4VKfHHE/cHtdMGPaNVmQVOvvUFS+ir56UmrUvIzc7VydBPSZzjQ794400zl88kGWtKonxErvG5LGUnC5oS+C5Ba+7IuG0KxZH7bicVhbUSyZPi7BlQ2q852WvVLBebpQ0xjIOUJkhFXevX4w5UCeD9z79ls1Fy26gQUgoxT9lHV9y1q59d/c9JoFLQn/PkW190sWuqpLaTb31BPourn2zPZssQnZmexGP7tjtzDeoLDe+70Nsy/ezbHfZT851wdiqvfPQ1+Ztr8VuRkpx/T89+hsOn5qPGiYmZ6cmu0oY61grdwlAK+6mpKZw9e9b6bxCgtN9NG/3FQ4AwQRZjx6WSpB7bt917n5CBqgi/YvlxqD4IFd55KHhVIhGF2P7nwPacL7dryQr0hhOaz+/rp5npSTx852SXEKcyRJ+b2Y0XDk05E7Fc78Jsy6XFJllnlpPY5Kt7GypoZ+fqTr9VDJ55YFdUjeNhwVAK+2FDiAAwBfvTs58FCbKYbEIqE/K5md1e/nBKG6TOCQ2TVD4Hs2C23gehwpv7PqhFdmaa5qqPLepuE4bmvY+/0Rs+KI3/FWzPz+mn9z//li3gKPoLE9y22KDI5J5/aDdqRCY14M7+BcIFrStfJTbUl5pna8FeDwxQ2Ash7gWwtx2VM9TgCgBKcwsRZJSWfuS+2507hJnp7gIlwGoCiTmNTEerOYAf27edbINvUprn/OrRPXj/82+dfRCqTbvoFxRs7+LIa+cw/ezbuPXom1i8vtQK8TTaSwkPF1mWuXhR9+bSPSiYz8/pp9C+NN9/3raYUO9kZnoSmzbSLsKf/vAm53VCBa2rfQKI5sIx59laEfTAYKNx3gWwZVD3DwHXkRYSFeOafOparrDDkJJ/OoOkmV5OObD23rzV+r0v+kYJDP0cKkvWF7NvW2R99AsKLjI5oGV2UARulxt+JtH3P//W+r2C7tSk7h0K8/k5/RQT4ROTYcrNW9DfiUsAu6JxbO30wdU+iXwRNLFO30Fn4A5lNM6gEFq9yXxRIV55Uwv1cde7EnDMti1eX7IuOmaIKKeIiQkq/FRhvCJw8uCernN9AigkKuX4G+ed9AsKnHehCNzmn/mJ91jO9dQxRURn2J6f0095I3wO7NzWY+s2z5+dq+PKNcuCWxHYtGGcXDxd0UJFw5dXQb0jn0COzRkYhgzcZLNvg+P88m3fuFEx+uThOid9HPT6+ZS5wLxGjKPT52xT3Om+c3Sn7eFT89g4PoYtE5lzi+4Kf1wwvs9L2mWCcz11zA1Z3LRS44R6fo4pI49dWYX96oJeD7lVxzz12896ksvGBLzspy4TZdFQ/UCZHaldo28uxgYH5AkqKAojr9m7skdD440prerhOyfx/uffWrUFLqEYVcZvzFIInII5wGMiT/R2Uxq+7fyN46s1QLdMZPgffvDfdWmQSnjUqhm5vXVNjDEhupLFuJW1uIuC73pKaD09+xkalsImNmRjAt+7YTyoIA0nzyCW6oAyQ+omLMoxq6xULo01lqogFuq63J2OTyCHjnnO7/2M0R9pYV9GfUq9sDGHI4c7CK4SbeSagm0DPDaDUwkTjn3X1sffXV3C7/5w0XrthUazi5RLhy85x2aCclXWyioCV64tdZUJpN6VKag2VzMIgR5BregGfBAADv2FO3OUQlkmgTwOYB1FVZkqAuZ7q01kkBLWRZJbwc0G35wpOls6BiNtxuGEj3Ffhi3x6tqSX8Ojrm8Wm+Bqi53zq5l3K09FmOjf+5J4skr3NtmsVRrjrFSkXKYpy/cuzG2xK7Rw04YKIFuLCzfhZma6VYv3xloVlxtNTGzoNVlwk8xMjTkEZZkEOFFntiIoNlBhybHJcHmgTLAvHJrC1eYK+c5DKrjp4Jii+mnCojDSwr5Izo3YCWgTmEArWuTpWZrB0Gyn+fnYg7u8PgZK2KhM2eln38aTvgxEU7YZn2O3qSpiQgcnOce8n82G/eKhKdQmNpBJQBQ4wiokNDWkb/RFN9aU4ANHIDHXMqvg7Lfd2lRUbDTJ+v1DK7gBfJ/IMMToj7QZxxWe5auAYyLWJjczPYljZ85b2RR1ytkthM1eFc2IsYPGbFvNEEOqVimHdCu0fRx/gU3I2GzYvnBQGzj+lZ/+8KYO14sPobvGonahFDg2dapkoQ6KGK6fdmubqYuCuj9lqqP6PZQA0RyHajHqVyjmSGv2lBnj8X3bgxMm8tTRpCaQrt0+88Auq8lEVUeKSfSI3bbWFxqkvR7onrw2bWmMqfxSgvt3R3+EFw9N5doWx7wvzvPuvXlry0TkQd5do+t6IfQcOjhx4JzxTBHD5ZkjoQih9dDvr5t8XHWP85pgBmHSGmlh7zJjhCKPTc5lB9W1jpMH93RtA8149qLazLE72zJzFczJo29ft0xkLFMHp7JRnhBDW5y4LZ5cJ/qioJ5XTeAr191CZstEGIGeLxtUf3afEKEWAq7w4ZjSqLj5ftqtubsFG1Ec4F4sijDBDCIUc6TNOC4zhl5gm4PYsLLZuTq+u9oreBRMwRk6wFzaGtVmbiETPTNXIRsTWLxuj3DhXteM7aZg64/QpBiFLRMZ7v83P+iEM26uZrhyfamTwEVpeLrJgqtN+koYmnAlI5lmBJ8QoaJ4uCHA5rs0379LePcz9JLqs00bKli8vtzFSWSjEKdkgwBYphvfOBxEKOZIC3uXPTkmnC1GGNvs3gohWo9tcAF+igVbmz/+6iLb7qwyc79ZaKCajWGxubr1Vff7+KuLXSUSOdd88/d/JnMTXH3ge15KIF9tLne10VcGT2+rL3TPRH2hgeln3yaLo5vwZcXq795Fz+ES6CGKjz5mQigA+kkXQPVZVhmDhH9RyxMqyRmH1PW50U4xGGkzjm9L2o8MN251IReoLbiNZZHzTCFmLKVdvnBoyhoe2mgu45WPvg6mRb602Ay2Z3K2xlR/NyKKjAPdJosQ2/OlxSaOnD7HstG6TFbmu6fgU2zGHKY1ReRmayvXX9RvGzXVZ5R/zBwXeUxOnHFIReF9d3WptD4Zac0+Nhu0SHCrC5nQtSRbUe9GczmYF4T7u4JpwqCEDTf23AVOoo5ra6z6K39LVsHhpXFBUUtwFnRq12hb0Kl2uhLMfO+ouSJzFdfmmomKhK3PqLluLtR5TE6uRVW/vi0Kz4xmKxIjLew55fXGhGBlWMYihrjK3CaGClPfVpEbLsk1YXBLF/rgW4RcW+Mjr52LYp7U4SP6sgmIAzu3OU1i3LBUndZD9ScVjqsg0M0+yq0XTIFr2rJhGOgCgLD5Fks7QcsSdIVaUv1ZVp+MrLDnCkyzehCALm6SvPbHGA0iJKzMBvWo1DNwNVTThGETXALAvh1b8I9/uOjUqnUK5ktXrmHRYhKqTWTOfqcm8tXmcpSgt3HXAKvvSm3LTYFvvrv3P/+WFOqKW901pqix6hL0Nuetj7G0TAwDXQDQHycxJUtW5OribnNuK5TVJyMr7CmBWRECK1KSphFfVEOswOec9/TsZ9ZSa6FYaDTZzkyVXKJHpgA8E4YA8D/dthWfXrjsbbNOwTx1/G2rsL/aXGZRMpsTOUajtfEamf3Pfe9H7rsdf3VqHjbCC4kWf8+Tp+a7Jr8vUsYHm6YaamYyQZUG5LbHdu8r15a6COz6gViN3QVViSykQI0tmq1MCoWRFfbUVmlFSvzxxP249eib5HmDsD8+PfsZO0LGh4oQ3mewZWJ1XoMAACAASURBVPu5tCFXGCdXuKh3QjnRKAew3u+2iRwj7C83mnhSI8sC4Kxp6nrvH3910SroFZRSQV07ZltP7TrUbzEa/jMP7Ao+x7y3KRAXGk02H3xR2njRUUGzc3UcOX3OWmfBBz2arewIpZEV9r5tpev3QdgfX/no68KutSxlVPm6mLh3ipbABl/fU/D1u8+2bYOpYd+QjQVXHZudq5NUGFwoIRAqnKldh/6Opp99m90vtWqWWwjNTLcoNsx7UkV4KDNWnp10GYyhJ9/6winoJ2tVXLm2ZB0HoZQLeTCyoZe+0CrX7/1M+1bwOThDttiTtWrfnoF7PYHWxLvtqX+wZudWswr5jL572KgmQtBoLjuFoqsQRh5Br67tqn/rgi/M1lbE2wZFrFcEQorwqNDMIrNNy8hc9WU4/+7oj3DswbCC6WVgZIW9L93e9XtoDG4sV4kOF8XAZK1qnbg2DhrVzn6lrrtyGYT2v1rKdJOGWbnJ9oycNiuqiRBGSi4o0q+8TnRg9dliqZABfyEaGyPoi4emSmNnDOFj8pmxYnbSZezKXcqG+s3s61o1ww3ZGA6fmo+WCaEYWTMO4DdNUL+HePSL2jZSbIqVdqFtqk2+dpad0WjaiVVYmmIVddmPbXVzY9s8M00XPgd6HWWcYwSAx/Ztt/o2OGaXx/dt78ksVvdQwi/Wvq7ALURj+74ImPbxAzu39TxzNas4c0KKjOQpIyroyH23W2322ZjolN4MzWwvA0IWEP9cJIQQZ++55557zp49O+imFAKKHTLGVmdGg2zaUMHf/6/lcmJTNAxFLRK3Hn3TG5L5xxP3R13bhIups9ams/XZsF3ONC4VMdB6d+ef/V96+veW71fJKl42jAlgc7Xlk7BFdphauX4/VbXJzBswfQ2KRpv7jvUFz9YmW5lOalFTSoEt0uuxfduDK33Z3pGtn0JhRuPUqhn+cs8PrAvbDdmYdZzltd/v378fH3zwwQdSyv2230dCs+8nJ4cJl40yFM/N7I4qY2eC2x+2XcmR0+cAuVptKq9W4nM+FulDcIUeLjSaqGYVPL5vOxn15JuMIaabxTYzpskzExo5tCKB/6+xBIGW0G8ur3RYNzeOd1tpzfepCxydx+jUf/26Ky9BUTuo9rpg3sMWYfTKR19jRcqesUclO81MT+Ljry52KTo6gRnAVz7KirO37ZDuPvGe1TQVm9meF+te2JdVr5N7b1u8PtCdTNNPhPSHTXjZog64Yae2RcYlgKtZpbMNri80IMRqQpipbboSknSNq5qNkdE5ShBRMG3z5j1DFvDaRNazvY91EqrxZTqDzbBG32LUaC6TCx2X2oGz4FGJiup8mxB+//NvrQvH8TfO46rGa8SZ32XE2dsQKrzLTjBb92acIs0oIeBs6fsZdqUQ0h8+E4sOn7nFtX0G7DZ9m31XR2VM4FeP7AFg1wrv2L7ZahLJxkRURu3jbbPB7Fwdf/ufe3nrKbs/xx/gslsrxGbAqncb8j4pvHhoCgAtlGPuYRt7eRZS6ppFIMRK4DIbmijClDTyZhyf970sEw9Hw+k3L4jrnrbvQyaZTyuhQt6OnTmPTRvH8c1Co6cUpG0brGN5ReI//Pb32LJpo/XalO27uSJJ/hIXj88rH32NvTdvJRNoKCGXVQSyypizqEmjuey895aJjHQE+qDebZ4SkQouMx4AcifLaZ+CbffJWTBd1ywCoVaCENqRfpiW15Wwtwlul/e9TBMPZ7D1mxfEZVbS2+JysNnACX+k+mOh0eyYH5RP4NiZ87jcaLIm92JzBYsRE3tZyh5tWjkPKVPGspTeBBobri9LNJf9tvxlKa27jsqYWM1ejVDN1bvNS5cA0Ga8Y2fO49rSShThnTkPbIoBRS2wcXzMmsuQd27ZZElo5jwnY5lbDKUIrJs4e4ov+8DObWRsdhkJFgq+wdbvhIqnZz/D4VPz1sloFsJQ/Qh0x7tT4Gw/uZOvuSyxwBT0CjHx85O1Kh6+c7JzbkUIPHznpNcBHqsZc55nslbFyUf2oFZdTR7bMpHhV4/s6djcQ81P+rs1Y72LxEKjGbWI2OYBpRioUFw9/j80WYmT80LJkhga9JnpFt8/VarRVLLy5uO4sG6EPSW43//8WzI5qkzaA1sykZkkpByIZb5goDWIKAK1ihBdwprSqij4ePcVOLVLY6G09BAof4Ba/JalxOuf1L39X0ZiFtCdnLVp43hnrOpOaM64rGZjzoQoJXxeODTlfBZXxnJRMMeeAqUYKDu8XijFlxypg1tAhZIlVH9xFBlfEmM/irusGzOOS3BT3vcyaVd9SU6HT83j+Bvn8d3VpcLCGCm4inasSNl1r9CFjrs7sfXH4vWlYM4aGzgJWjoEYDXVqF2dyxG6LCWyiogivXJhYkOlp3yjOR58NneBVbK4Fw5NkWNICRbK5DKRjeE/th3nptknq4gum73CmGiFgXLhckiG1njgRtdwzTDUHLCZ/gSAW75f7YmsMtvjizbqB7niuhH2MYI7pnBICGzMkXoRDSr87/gb59lx8JzjXALc7J9QJ17IQLT1B8eG7PIbZFoG8cz0JKaOv+3lo3HJpG8WGnjh0BQOn5q3HqcWFjOBhpOQ5cKV68teVk0quci0Z/uUBl/wgGzvQV0Ki0nwFiLoKY1eoaxYeO5O3lU97sDObT3x/noggKvvXYtSP8gV142wjxHcZQ0qCsfOnGfZXC8tNnuKduvtBcIiA1yFRcz+CXXi6RrNgZ3bOoU6TFoEDu1EbSLr2umoNjp7TNtZz87VceX6EqvdFG5sm6XMJB6gO8nHfJ6Q7FkKPlZNarzadjQurdC3mFNU1xwqCFUPwvXOfvXoHu8cKyMWnqsQumQJp7RljEbej+Iu60bY2wSHlC2KXcVJTgmcfiU2xTIg2gZPyLbPlW7u2m5yNHy98o5uGtETZ468dg7H3zjfVfFJ3cem7YfEV+vJPjGRMibU4vfczG7svXlrVGZmTKigC/qEt41XKuvWphXOztVZbdPPDSnMsSxb+wIqjLQIquRYcBVClxLIpeyuLzS8ph2F2bk6rlzrVVKKDuJYN8IeWJ0Ig8yaLQvmxA3Z9oXuYFQ/hvCdu9Bckd6dinlvBU5iSn2hUUjCkCmIQhUBmxac1z+hdl+uDOGQ8nYczVQ/N2bHImGn5NapkgdBYRIyD0L9fCYUZTfgHvNU/4byEXGwroS9QtnODlvxZ05iBJWmL9qztTaRYWHRHnbIta1T275QwTU7V8d3V+0mkVo1y8XTHvIuuGYlnwCbyMbw0J3/kszILZKzHbDvWGIqZqnnopQXl/BevL6EW4++2SXUeBE9qxplXqpmZdYpqxhJKPLu5Dnj0bb4UmOe6t+JDeOF9wUZeimEmBZC/HshxL9r//v3ln//TggxFXNjIcRBIcS9Qogn4ptvR5nODjMO3eT5cIVK2YpoZBWBFx6dwguHpnC1uWKduLbtXNl89FRMd62aYf6Zn5Bxw1xw34UKrcsb8thormDvzVs7YXrAahgli7P95ZeBW24BxsZa/7/8ctD9Z6YnrfUFfJisVZ3Ki6sfL7UVB2VKm52re23AAsDDd64KxLxzRpX5VKGSQDkFRPoFW6jn4/u2d30OqWjWz6p3Ls1+q5TyP/kuIIR4GECQyiKEOAgAUsp3hRBPCCHulVK+G3INF8p0drg0HZ/G6tpGUtQAVOSCciK+8tHXWJaykxRUlDZADTZVHzZvNqb5LlzbevV/nvtJtPpdFzpsvPwy8MQTwOJi6/NXX7U+A8Bjj7Ev8z/u2BpEXwy4nakhvDHNFYljZ87j2IO7nP0oAbz5+z93KIg59Acu2t4xIXoI/wZR1tOGWFOSb3dAmR5t8qcfjlkFUrOXUv6f6m8hxL8QQvxI+/e8dtzrEfe9C8CX7b+/BHCH/uP8/Dz2799v/cdBmVqvb0D6fp+ZnuxJDHGdZ8bBK8zO1aOSgrjwlS10aTjAqtZcq2bWkoDKgTU7V7cmlDx5ah7Tz77deR7zfnqWKRfRwuRv/3ZV0CssLra+D8Cf/t9ihZkSUNyEsoVGs9OPrv67tNjsvAufoFdKBlXicFlKHD41j6dnV/lzakSyFvV9GSgziSlE/vSrYhzAt9n/EsAnAJRasiXnfWvG5+/nvF4XygypLIt/PXSFj/VLcLUZKoJHr4fKtX+aRVcU9GLeNk3z0mI3Ra9+v7tPvBfsN7ghi0wYv2DnyiG/J1Ck5po5KpS5xqfSsk++9UXu+rjAqpIBgHyPEsDLH17ocM9TvqDvri71jfa7TL9eqCOYe2xecIX9a4am/0nO+y4A2Er9ODU1hbwUx2WFVPr412NXZF9YGDck0SVQQhxjtlhzvWBESN/auMgVXMUc1O/H3zjfc78YwdloruCWo2+Gswxu394y3di+D0ARrJMdaJslc6y7oqiOvNYqQlLkwtNo2hPCdCgzGtCbfavQXOFx5heBsk1JIfKnX+HfXGFfE0K8BeBTtIbZNID7ctz3n7Cq3e8A8E6Oa/UFurDd3C4WfGmxGRSN44K5wm9uZ2VStAoh4XYKIdrM7Fwdr3z0tTWq4Ml27oJKovpmoYGJDRUsXl/u1E/96Q9v6pCK5Z1AlxabmDr+No49uBqKlkdwBkd//P3fd9vsAWBiovV9AGIpim0wi4mY45M8ry1QXYl2MZz8nCfijIN+2e2p55do7Rp9cfGDqnyXB9x97Q4A/xuAlwD8GsDpPDeVUp4GsEMIcS+AWpHO2TJg2vcWGk1cba7gxUNT+MPz/xZ/OnF/J7suD6GZsue/cGgK15ZWOtEUlxabPdqQjY3St7PgajM+7hRgNYlK9cmVtqAHWlv733x4oWOnLcLZpKouKRu/LQklBEHRH489Brz0EnDzza042Ztvbn0OcM4C7WpRB/d0EYzVqpnV18GBXpPBHJ++8yhb8WP7tpPf543AurFW9Y4F3+9FEQe6fB0u+30/CMvKAlez/0RK+Uf1QQiRWxOXUv6y/Wehgr6MVdenERcZNzw7V8fPXz3H4gXXj6hVsy7N1wauXyBvbLXCKx99jedmdhfCow5086YX0b4gLfKxx4KFuw22xKuXP7yAG2vVThUobl+p9xb6vhQlhDrXNlfMKC+1S7MlAXGzhZW/h6oX4FNWipxnvkzxkLj4ognLygJX2B8VQvwa3Wacf11aqyJRVrKGTyMuagBwNGoK15ZWvMdQfgFV5zW2BBwF9RzmxHJVZPKhCKeiwiCKxyjhWs3GsNhcfWf1hQZJvkbhm8sNPD37WdCipRy7gN1WTEV5KT+NbZHwlY9UeP/zb8nffORo6p5FClr1PFT2tcrM5iSlDaLqXCi4wv4XhoN2uqT25EJZq65PIy5qAOTRqBvNZfz81XM4fGo+iGbVnKiFORDRbZIwNdoiNP28+Pa/XcXU8bdxudHL2VMETDv6letLHXu9LugVfIK+IgDd3C9lS0ueMBYOFw79xU25Baptkdh781Yvf47yNdmw0l5cXHwyeeaZa8fvUnB0U43r2H4rDjFwZdB2MmN1Qd/+PGc7btCIiVDhwBcL64tJ58LVzqwiWlS6jvOX22yDLjuiGef//uffliZ0f/rDm3q+U2Yq6p7VrILH920vtHAGZWu+rlXFKtr2arOj53XMUqc3llbYsfYu7RrIJ1CvehacihDknKhNZF5beOw889nZObkKasHrZ1x80XA5aG/TqBJMuoTO92g5bwcORQhlQ95V15ZApG85ixoAVDsrQuDkwT2Yf+Yn+OOJ+1mOMq4DMnQhpNLEN22odPq/IgQe37e9p8Qfx0x1QzaGvTdvtVJLxEAnpPJB9VkRTsCi/B4cSIme8UnB975jBSrneZelJOeKlPBSKMTOMx89gzm/KahCSNzKWDb0ozIdBdKME5kZOzBQhFA2zvYYuGJhi0qMoBKZVKFrda8j993Osu9yi55zhaEqCxcLjkC4tNjEkdfO4Xs3jBcSohh6BdN2rj5//NVFsj6tjRivn6gI0WMmo8aHrbC8PmY5NMC28zhjTS9haZ5PUQfr142dZ5zdis9hC3Rnj8eY+wbNxrtuWC9dRYr70ZGUsytkYJoDzlV9iMOgyK2NybGfF7FV5e4idEpkF150VJTKA/N6Ei3b+G8+vNCVV3Fg5zb8l3N/7nIalyno//V/vwn/1/9zped701zGUXwowfP8Q7vx/EO7yXFLnVcjGF0VzMLntigXm5CtTWQ9dvxQhYNjZ/f5kTiRQr65PuhInnVTcNxVpHgQiI3HVTZ1G3uevvX0PRdXOLu4UmwF0vOgSCeW4p7vrw7dzXL6mw8vBEUHxbBeAqtmsXf+aj/uvq078fzu27b27Dg4io9P8Nj4m1znSYkeE0vI+LGZaLKKwHdXl3LHtHPMP65dp6/93Lk+6EiedSPsh81xkpfG1TcwDuzcRtoXdQZMjo1wZnoSmzb2bvIkVk03RWgeIaRdQG/SmEI2Jgrlnu8HKkLgf390Ci8emmKTuFWzSidx77mZ3Zidq+PTC5e7jvn0wuWed8pRfGIFj4sN1bRlv3BoCn+yLBg22GzhmzaM9yQTxlAhc+zs1HMJwNt+7lwvKpAjFuvGjNNPQiEO8q7irq2nioWmtFqdnMoMq9RNQfrWM4SDOxYuM5UN1G/NFYknT807TVljIqwIdtnQ2Us5JGRmpSIq2c5mBuCQ2MWGELrOy8vxYp5/69E3rcf5+J9sMsDXtjwhldy5HlMnu0isG2EPFEsolDcTN288rq/osc/G3mgud7Igze+VxsGx1RetdZiOxGNnzheaKAW0BH0Md3yZuCEbw21P/UOQTd+0k1PnmkKFQ2IXK3hCCfvyKFyhcyiPA5SbcGh7Hm47B62QrhszTpEogv8ir1nJtfXkatsu4cBZMMrWOmamJ3HswV3IYo3ZBDZXs8K54/Oi0VwJEvS6s9P3rmzCz8Y06go35PplXOc9PfsZDp+a984bbvhh6BzKYzq1PdfDd07i9U/q3ucJaafOfwW0iA77FYK5rjT7olCE17yIVZzaqYQUPabC71wLhgDY7TWzRIUAFhb5GalU+cM8WFhsYoERzdMvxJqUFLe7611RQoUbbhgSH+4ay7NzdSvNsTlvbNr3k21mV7PAdugc8j2z7xnM/pg6/jZLDoS2c1AhmEnYW1CU17xIs5IObrjkxIYKVoxkFd0UZFswQmLpzUGrm2NsAzg2PjsUnFJ6/cJkDq4hX4q+i0+GOsdFf0yBI5yocE+ge95QuxSzUI1CyBzy+blCBOzsXJ00L9rGLNVO25gfVAhmMuNY0G+v+excHVPH38YtR9/ELUff7CrHZwM34+/K9VZhELOwNgAsXu+lCA412/jMC6rwyN0n3sMtR9+0bvHLKEU3LIJe9SdFW1wRwkkL4UrRz8YE/kV1nDQDHLnvdqt57L9dWwo2GXDMI65FW583ruPyFh13mVNCTTyudnDlAGUOLovWxYek2VvQT6/57FwdR14712XKuLTYxJOn5jsFom2rva5JUAWOFZal7Gq/bVfAoUg2wRmclxabHfuzbYu/cXwM1azCohZQCU26eWoiG8PGrIKFxeZQafRjAtg4PobDp+bb5fp627VvxxY8sne7s6BJvZ2iD6CLaExPPFNZvk+emm/xJ7VNaTYsr0hr5S8d3KpoOisklVRlZrD7TJB5BJ7LnMLJ0OW2gysHqAWGyrIuOwQzCXsL+uk1d9msVcEOvU1A72TkUMzqWoztuE0bx4Ofrwg65MuNJl44NOVMUwdo85Lqi4XF5tAIeqBlo1dmgAZBEPbphcu4ddtFbNowTpoM9F2Bi2hMPTknssmV6Wozd7hCZJXGmo0JZBXRs2hNbKjgcLuyGUXHoKMIHqsQs5ZrF287fstExp4n1IKhlK9+h2AmMw4BVxZhkfAJS3Oradsa/ubDCxCQ8BU5+mahUWgWX2iSlA2b25mwvzv6I6dJw9Y+M/rDh00b8rW1aDSay94sXLWA9YtYzXYfW1U0E80ViU0bxjumxVo1Q1YRnQpmuo2cytguU+CFRvbYzGDZmMAzD3Qn87kii1zJbXnI1GIxMpr9sNaN5BBncRxcHD5zNfiK4uM2d0Cbqxmayyu4cp0vlPQFyhVBY7aPiv6g0GJWdB9dq2aFx/znhcp6LdKe68rgdVEtTLajuKhevNxoYv6ZnwBomRbNvlSKi1Ke+jkno3br5gpnfPY5fV3m4LKCN1wYCWE/aLY5FzimB66DywWfzb6+0MD0s2/3hL/5YCZJqX7lQhfwmwlha2MudUV/dM4TLfpfVQzelXFbzSo49uAu/Iff/p5dCMSGPNE3tjbpNROKuK6PaoK6j25Go3xEnHEaG/pZBELud/KtL3pMUmaRd19UzaCTqEyMhBknL09NmQglNIvRwGvVrLNNVJE8NpPJpcUmjpw+F53gEWNq2FzNOtE6lFb92L7tPROEs+jdMN7iluGY4VT/UPZ1Doqi0wZ4NROE8b8NetGbyVoVJx/Z4+wLG+eSOQY5JpE8hUYGxfeug7NYcRe0fpiDORgJYT9otjkXXJNYn/BqErhKuwEtB5JuC3zx0BTmn/lJT+zyxAb7pk5pLzEI7c9sTODK9SWnxrplIrPyyHMWPXNBp8wXikGTe10KN9aqzr6rZv7pJgDrAmXL8FREY4/t205e79BdN3WK3viEjY1zSQAdUj1XW0ybs87Do4P6Xt0/b+Z6UeAsVoMmNgvFSJhxhrluJGerZ5qhKPNFNauwzTAuwRy7CPpMAKaNdvH6kpe3nrLjcxPL9Gc59uCunjBX06zBLQxjQmm2VIgf4C/bJ2DfxShQZghXqUFfGUIdlHP2//joAv7LuT/31Op1jTPqvq72UDvwY2fO990Uwgm/HjSxWShGQtgP4qWEOJ9Mga/zl6jvKaGmwuImtXtw7u2yAccugr5+Nu37eQqwmH1Gxdjr53MW1plpd2EY1d81ghrClZkM0NFXkzmEWMjCTY2N2bk62TY9jJTr74rZTVO/LTSawffPC+5Y8R0zTBgJYd/vlxKTmu063jVBJFoRPWqR+Piri05aY4Uj993eo+UCLRtv7CKo97Mq0WeSb+nP64NvQbY5h30LOsdJRzlZK0LgV4+6bd6+Bc+2G9m0oYIr15a64tFdeRXqWr6FDnBXY1JjQ40ZLjip/TG7acpBH3P/IsAZKyFO30FHBI6EsAf66/kP4b7g8JT7IjHM6kkmFG1BD0fHI3u6KIZNDvUYqHNdi1foTiXkvkVMJkpgP//Qbnz81cXO+6oIgZ/+8KYun4KvHR9/dREvf3QB+utuharSNQfMvjxy+hwg0VmoKUFvLtzUuDTbw4HP1Bezm/bliYTcv2yECu5hiAgcGWHfT3C3sFyecq592gWdtkANtIfvbFWoMm2xeeFb7FwT9YVDU9Ft4C7oHPZD9Rz6MR9/dbFrMV2WsvPZFPjUDu71T+pewar3la0vKWoFPdPVtnCTMfQRicc+U1/M4hvCVKruPwht2Sa4fUXpB11/FkjCvhRwt7BcnnI1GIos9NFoLnclJRWpafgWO5cj17bzKXIyczUsm8D++avnrNd85aOvyUmuIyQ0VfVVqAb7pxP3k78VFavPLeph9qGKKKPOCaHuPrBz28C0ZcqR/fKHFzrFYUwMQ0TgSIRe9hvc1OwQnvKZ6UnMP/MTPL5vO8miGAobMdmTp+Zxy9E3cdtT/4CnZ8MSpBR8IWnc/rGF4j15at7LCupCnpwLage2LCWrPSETW/VViLPcd2wovcWWiVZ8/paJrCtW/47tm/Hyhxe63svhU/PO8cIJq+S2T1Xd+tv/3Lvb7Uf+jCvLmLr3MIRpJs2+BHC3sKE85coMEEL4NVmr4sq1peAdAWWiUO1wPRsnKgfw908o9zkH1ESlNEr9WV3gtCdUcwXsfZlVRJfNHuBFl9n6nRobWyYyzP3dT3q+n52rW0NTJYDffHgBv/nwgtXfwjFjmO2rZmNoLK1YzUyuHVJRGcwUXO+RGifDEKYpfHwh/YYQ4uw999xzz9mzZwfdFCeKMC9QESQUKRKVpq5oAbq+Qytm+7mZ3db7+Ip9K1SEwB+e/7fBbY7tH04RdIXQQisuZk2BXn+B7Vld2LShgtrEBna+hMKGisB1ww6v96kvGiemQpirTTFj0IR5jVuPvml9nwLAHy2mp9C+12GOWR9inK1ULoZrTJbtX9i/fz8++OCDD6SU+22/J80+AkXZCkOdWKR2KYHH9213Fpg278OhRQZ6TRdcR1NM9FPoBHdp22a5xCvXl0jHJrC6Bfdpoy5cub6MK9dbbbKNiZnpViFwVQheRfO8//m3PQLU5Fix9SUVrRMyFgsbgwbMMREaipmH5TNk5xvTf+o9mkR8IaHCg0AS9hHI61mPXeGpOOQba1VngWlKYOy9eauXR970D5TpaAqd4JSgcJVLdEEvxlFEyURzTJhmuGUpnQsu5/5FRHmECKEQJ6/efiqibPH6UqfWLnVuKFQCG2eexfbfczO7O/NnLSRUAUnYB8FnCuAM0FhNbHaujiuWUoLZmHCm6bvapCY5tcUGgJ/+8Kauz0XWNjURMsFdWlQerVB3HlKLawgVso+eutFcJk1qep9SgqvfUR4hYcC27GUzouzSYrNTZUu39cdGDqncAu48c/lwbIuQjkFr6qFI0ThM6NEEFDie9dhoEBvlKgB874bxzuSIbRN1zEQ21uOcpWqbXrkeXtuU2w5F6MYt9lCEoGs0lyEErFFDxx7c5eSE18Gh/aUW2oVGE9PPvo2nZz8jI1moPiujti/QTYIG9O78FES7nTpz5cx0K6/DhBn+OztXdxIEUtgykeHkwT1kfoJtnrnmx6BI2MpCEvZM+LRFrmc9VhMjeUPaiSihlXh0UOf+x4f+Tc+xM9OT+N4NvRM2D1umrx1K2+NSxRYVzraw2CTZHY89uMsbJlgEPfWlxSZe/vACKbiO3Hd7KzrHwHdX8y++NhQbLQAAGOZJREFUFNS7+NOJ+/GrR/dY+8EmwAH/ONdNKGbfU4uiQCu/YO7vVtldufPMFe6pQpEHSbVcJJIZpw2ffc83SG9g0NcC8QycvvPy0AWEnktlOiqbd20ig5QIzswtivLgwM5tPc6zbEzgezeM49Jik1UdDGj1rcs5arb1wM5teP/zb50hqRzyNxNUS79ZaBUjtyXbNVdkkN0+FhyTWQj1B7A618y+5xRN0b/jHKuu73ovw1TsKA8GJuyFEAcB/ExK+eNBtUGBY9/zDVJu7Dc33pZTVFydZx4bQzlQlJNOorugdUx0iB5uSJGDKZjPfsv3q/jHP1zs4WQ/9Bct38PLH15gCXoOIVyozXZmehLH3zjvpXXmQgmuy4T/IE+8OTeIINRZe2DnNit/kw5K8QmJVQ85Vpl9XM/Sb2qDMjAwM46U8jT12/z8PPbv32/9VxT0ijg/f/Wc177Hye7j2N5t21NbjLppo339kzoevnOy5zwAfS/4EJqJGZrVyC1iYTvud4agB1oL0Ju//3NQzdpNG8ZLmdjPPNBr/uHkQ7uqR1HCUQBR48DWr4fbmdW6SSPk2qqNPn59n+lx4/iqyNoykZG+G84808EZ04MmX8uLkTTjmJq8j4gM6N22u7bWPvg0Qsq5pLITdc397hPv9SUM1Dzu4TsnOyYLjgANmSjccLiQqJtQbZrSlvPCFmu/b8cWfHrhspcJVJmfzAxVquCKLXeAA4r7BejeqR07c551PV2Au8aBi+nUloPhKwajzzN9p+hKQnNp+MNQ7CgPhlLYT01NocwMWq6QsNn3dCFbVvUr14Tgct0XGQZqO+71T+odTYmTVRnSL9xnKlPTsrVXX/Bi/RK2WPtPL1zuWjzVtRcaza6wzGUpuxzWCjPTdMEVbh+FZC6rhdcVfjpZq1oVCF81Mwp58glceRfmmLclqgHDXYGKi9KEfdsmv9X4+ksp5btl3ZMLzgTw2dE3VzNkFdEVDlnUgPD5BzgOrzEhuhKEbBMijwatOPJnpie9sdeh/cJ1roXEYlezMVDFxLMx4eWZMQVArF+C6sv3P/+2R9jZFlFKwFEFVziLbAw1gW8OUYI7liMmj1LjU+58HD2+xXzQRUm4KM1mL6U8LaV8yfg3cEEP0BOgIgTbjr7QaAKymxlw4/gYDkeGauk+hCvXlqzhdDp0rnubrXFZSq8NP68GfWmx2Uk8efjOSavt2WVX1Z9Z7zNuGGmI72BpRfbkBwi0aCZOPrIHW7S4dN0urMAVGD6ECK2QY/OE3sYkoUkAlnQLAOjqSxOhtnSFPHkknAXBPEYpMTe2dygn3/qiaw6psXvL0Tdx+NR8X31msRhkNM69APYKIQ66nLVlwFWJiBp01iISKxITG8bxzAO7cvGT2LaZ2ZjAlomMtDVTIZe2EnW65qJrIZy6reozpUGr69roGgBggnB0ckxIPm2JCn9U9nAdzWWJLRMZJjaM91xzdq7eZf9daPRGVsUIDBs4uxb1jihzCiXgbsjGOv1Zq2Y49mD+4vMATZq3Yvkyqwg888Cu3h80xGSe5mGN5OwAzT51jU+guxKbi6ZkmDAwYd/W8rcM4t4x8dwuLSsvP0nIQgLYue7VfW49+ibZTo5jmtKgfTbh0G22r8+4AsF23MtEaN/CYtNK20u15eevnsPhU/O4sVZFzbHwKnC0TJ/Q8plUOCYmALi2tNJzDDXefXZ0H01IRQisSFmqCSNPDkaMmdGXgevbCQ1j5M5QOmj7gVDtwqWR5eUncZ0fOshdZGnUdp0zWSntTi+yEWIzLpPTpai26LV9szHR46PRwdEyldBsNJfJyBqXSYWKVvEtnL5dFKf+wMw0zaG0IqWVprhoxOwI1HlAGCV03vE5jJE7IyvsQ+GaEJTWw33hnOxYrjkohizNNVmVoLBNcl0ghG6zYzOJOc6wotqio7kiUatm2LRxPCoax7arskXWUMJEoNfpySXm4+yi1HGufo19Z8OAIpU7wJ1MNqyRO0nYM+GbELZqQleuLXkjYoDiqthQZGkbxsecWYKuyeraDeg+jtAdSMwz2zTUI6+dw/E3zpMx07FtsWGh0ewQeU1sGA8yWfjMAqqt1BbKZBXlRNCo98rRUjnCsKhxuhbge1aqGJArV2DQSMI+ANSEMIVLbSLDd1dXy735HLZ57JE6qEl95fpyh0kwdLL6TBw6QrSnmGemfBvKlk7FTPtoF8xEJwqKydF2Lx9cVLpd74S4vUkuGULMF+IUDnWID6tg48D1zJxnXWv9kMoSlgAqySikjF6R99XvHRoT7LqmL4IpBq72uXj3dejPyi2h6NOSKZ8F951S/cglZVNlEzmJT6Z26euH0NKEax2t5/19T97FWn/mVJawZNiEU96s1liNgRM1E2q7dJk4VGJVURqOL9yNChU1wbVVq3v+/NVzzutSCUuqjRxQuypufHttImOZmmyLj09LzRtNtpYwO1fHkdfOdSXRKazXZ1ZIwj4HKOHkioiJuR7AZ4y00d1y7u26JkBTwF5abJJmlFBQQufYmfO4trTCri2qsod9/EWqv32C/ndHf4TbnvoH63FU8Q4TlMD1sS0CrUVBSn+4n8sk51rk+13tapA4+dYXVkGvsB6fWSEJ+xyghNMN2ViP1sZxZBWhYR170B+XT4HaVbicuybyaEdkgRaCg0UAGLeEQ/oWBbXwcTJHVZuoa4YUt6YErktjV45wKpIKaPWDqa3rkTpUmKfCWo6y8cEc00VyOK01pEpVOeCqHhWTEl6EhhWbju6jFT6wcxu7DaHakUo9j/EenTy4p/OsHC2by8CooHYJrkvnocc4fGoeG8fHsGlDL+1DNavgV4+2yuxRQmiyVu2p3mWW0NTzBWyp/HmoFoYZtjHtGiECWPPP7ELS7HOA0hRqE1lUAkgeDSsvGZNvV+HjIQ9tr97uUBIu/T6c7GHArv1yNL2O5u5YiYqgx6hmFTy+bztZ6Sokksq1Y7HtvNZblI2Ci6rZhsf2bV/zz+xCEvY5cOS+23Hk9LkeM4Kq/9kv/o+8tn7Av6vgauuhGmEMCRd1n1D6XFt/m9zxXHDMV64kKIr5UoESyEBrZ8EJEFCw/R6bnVoEymKNdPWDEIB6vVsmMjzzAI9HaC0jmXFyYGZ6Eps2WIpvr/iLb9sYH2NNML6EHQ58rIIubV0xf3LbqyPGIUbdJ9QcYevvFw5N4U8n7sdKREiy61lM00ro+aq9etF1wF6lrOZgnQSGyy7NrUoWA9dzSrnKeqoXKl/PSJp9TlAVjTgTn9LEi0qmChGkvl0FpQU/tm87npvZHdReHaQprJrh2tJKUOx3jDmC6u8Qrnz9HB0chlHX+T5Qi/zG8d4AAYVhs8WXGfbpy4yWaJHm7b15axL2CX7E2NmLHuBFRFNwBOXG8VUKXe7W17dFpxaZYw/u8raHeo4iJi6XQkFvsy5EuaUvFQTCnOAAvZhfbjQ7CVicaJwYFGV6KTPsUx/T1MIdW7pxLSIJ+5wokoIgdoAXxVlCCcqY+p+282y+BN8iM6hJyBEUCgLAw3d2912oL0ICeP2TepCW6Vrky7TBF+EjUigr7NNcjFy1IdZzbL2OZLPPiRg7O6fqDlXFqag2hCDWJ8A9z7RFD4uWpdr14qEpZ0UsCfREK8UIkFA/y6BCJmPGQ96qZCGw+QG+u9rLBqswTD6MMpE0+wJQBAWBq4AFR3MqU5OL3Ymsl8xM1a8uWgXzmSiNVdUO8GX3hrSr3yGTnPdq1mu+cn2pE7UWU5UsBGQxoHYdYr3vh82HUSaSsB8ATBNBRYgeutth4iqJ3Wqv1cxMyh7tymI1C7xTC/rzD7Wc2YdPzTuLwXAxiJBJ33u15RGYiKlKxgW1GDWaK11kcusln4CLZMbpE8xtLLC6hTUzHH3FKPqN2K12ni26bdsfYtqKhc0E8OSpeUw/+7YzpNEs8A6ANK1R9WX1DM5+PGssDuzc1pOJqr9Xrr8iz3h29Q+1YNYmspEV9EDS7PsCyiyjF4hW0EvWmShbI3Zx4wBxkTEx59n6yyRiqy80cOT0ORw7c55dLYqDY2fOWwXVpcWmtzShgtJaqUWNEnISsFIO5yWYKxKzc3W8/km9a7EyHdRll+6LKbOYVQS+u7pUGGnfWkQS9n0AZZahtB9Vsq6fFYE4sf8xkyLmPK5m2FyW7AIxHMzO1UnSNaC3NKErFt8sSqK3z5XpC5RjxisqVJKiINAd1JwcBWo8c9oZU2bxyrWlnne73imNTSQzTh8Qul1VW/6ioms4JoEisnCLQuz2Pm97OedebjS7IocmCe1U+WFs7fOZt4p2bBeZpcppm+35sjHhzbTmtpNbZlF/TzHJj+sNSbPvA0KzRE3zSSh07UiVSFQc3pQGPAyRM6rdeWqn5Wkv59zaRNbFRXPL9+3v1hW14zNvFe3YLnKnwGlbrPmO286Y/lmrwQJFIgn7PqDoLFEXTHOMLZEkzwQqi7QqD/uljjyT12d+sNl9qQXC53dxLeZ5kuSKrpwW27YyaT9i+meUiqVTSMK+D+hnlmhsJAQ1GQ7s3NbRZDnx0mW324W8k/fAzm34zYcXrL9NEnZfaheSx+9SpGNbEaPZFv2YhTG0bebic2DnNpLGmatwxHIghZ6z3pCEfZ/Qr3jo2EgI22Q4sHMbXv+kzo6X7ke7dZjO0ryTl+LsVxTJLr582zmq7GBM+4pybFPEaGYSXxHkcSZsi4++mHKiaKgFMqZ/BpGTMExIwn6dIU8khDkZpp99u/R4aYUYlsm/3PMDFuMmV5j5zAhUGxUHvkIRfpcYcIjRzD4oM8yTs1vzRdFwdzSjrLFzkYT9OsOBndvw8ocXuoRPVhHYtGE8KB59dq5OEkeZCK1MZZuYoSyTAK2Jm/fTC8yo+HygV5j5zAiU5vnwnZOkaaKfiCFGKzNbm6sEmFE0Ifcd5pyEYUMS9usIVMLLobtuCuac54YwhtjJn579rGshsk3M42+cZy8ySki4NLvjb5zvSYJqLkscf+N8jzDwmRGG3e5bJANryC6L6n/ubi2PU33YqEWGGUnYrwO4yt3ZGBk5cGllWyYyLCyGZa3OztV7dhxA78T87hrNTmjixlrVq9lRC4fte44wD9E8+21eiFmMXKYpTmlNV/9zdmt5nerDEDK8VpCE/RoHJ2QxZuC7cgPm/u4nwddzxc+r9tm0cKAleMbHRCdXAFgVEnk0O524TI+MKiOUtF/mBVv7XYvOkftut5KycYt6uPpflU7kRuPEIMXP85GE/RoHxwkWM/B9uQGhcC04qn2UFi4BnHxkTxATpbpfrZqRFAgmcVnZIbCDMC9waDBM3iEFjpLg06zLdlKn+Hk+krBf4/BNyNiBb5oENlczCNGi5lUp/65JbGqTmwmhqzM9+tpju59Pszv24C4cee1c167ARBlCOMS8UKa5h7PoTObQjgetWQ+7H2WYkLhx1jhckyovp47iF3nh0BSuLa3g0mKTxa1i4zi5cn0J2Vg3Ma4qWK7aV6vaKYSp7wE/jfLM9CROPrKnwzNEoWgbL6caGVAsb40NLgdsEdWiiqw0FUvrbPLgJEFvR9Ls1zhcRTL6qR36jm8uS2yZyDCxgU6Csmnh2ZjAX+75QRcfjS37mOtUvfvEe33RRLnmBapvj505X4i26oqIKaJaVFGadQqhLB8DEfZCiBqAHe1/d0kp/2YQ7VgPKHKyUdcIjXigvl9YbDqdu5wsXqpoOfd5+2Xj5b4Xsq8azVz0zXqElpn0pVBUtagi7PLD4uNYzxiUZv8oAEgpXxJC3CWEeEJK+ZL6cX5+Hvv377eeePbs2b40cC0h72TzaVUhdtnZuTrGchRfMZ/l7hPvFSoEqAXl5Ftf4PCp+R6hnMeeznkv3Fh085ld7TLfp4tFdFhCFFMIZfkYiLDXBTta2v07g2hHQgs+rYqrDSshYxP0sdpzGSUadSHsWugAlG5aCMkc1pPIXO0KIZUblhDFQTt6RwEDtdkLIXYAuCilfFf/fmpqKmnwfQQnfA7wmySokn4VIXp8CByNeXauTpogihICvqItZZsWbH27eH3JyVLpW5y5C+EwhSimEMryUZqwF0IcBLDV+PpLQ7AflFL+rKw2JPDALUjhC7Wk4tlXpOwR9ByNmVOYOy9izAdFmxbMvrUlyoVUsqLep89BPkikEMryUZqwl1Kedv0uhDgopfxl++87pJSfltWWBDeK0KpcXDqmFs51xvkKcxcB30JXlmnBtbPxCb5YwrZnHtg11MKz7ASsUcegonHuBfALIcRT7a9SNM4AUYRW5dJ2zUWDq01TQo2q+xoD30JXtGlhdq7eQ/ZWX2jgyVPzOHbmPI49uMsbHbPWCdsSBoNBOWjfBXDbIO6dYAdXqwplONwykfVcl+uM64cdlyMYzd8AkHH/Lvh4jBYaTZYDuGjCtoTRgJBEYeRBQQhx9p577rknOWiHD5Qt+fmHWvTJNkG2ZSLrMR+4rhNC4lUGfPcLabsJKqHLREUIrEiZNPKEIOzfvx8ffPDBB1LK/bbfUwZtAhschsNjZ853OWovLfZqqyFmhn5qqBzHcZ7kH65jV4Wulp1Fmio8jRaSsE9ggxOiefKtL3qicmzCkBPd029BxBHkeZJ/YkovlpVFmugJRg+JCC2BDQ65VxGZkGWTg1HgtJ1LcGaDjTQsT7vywJdfkLD+kIR9AhschsM8wlBhUIKI0/Y8LI8z05N4/qHdHQbOyVoVj+/b3vlcEXZezjKySBM9weghmXES2ODY2ouIoBmUIOK0PW9Yo8t85UumKhKJnmD0kIR9QhB8tvYiYrwHJYi4bS/Lacy5f1G+jERPMHpIwj6hcOQVhoMURIOOTw/R/PM4VVPi1eghCfuEoUMSRHYUzfk+6IUtob9Iwj5hKJEEUS+SUzUhD1I0TkLCGkERkU4Jo4sk7PuM/fv3k1W4RgWpD1oI7Ycii3sPC9JYaKEf/ZDMOAkJawTJl5GQB0nYJySsISRfRkIskhknISEhYQSQhH1CQkLCCCAJ+4SEhIQRwDAWL/m/N2/ePDk1NTXoppSC+fl5AMB6fT4OUh+0kPoh9YFCEf0wPz+Py5cv16WU/9L2+zAK+zkA2wD886DbkpCQkLCG8K8AfCulnLb9OHTCPiEhISGheCSbfUJCQsIIIAn7hISEhBFAEvYFQQhxhxDiD0KIT9r/fmE55qAQ4l4hxBPtz7X2eQf149vH3Gu7xjCjyD7Qjh/ZPhBCXBJCvCOE+Ot+PkMRKLgf1HcH+/kMeVFUH3Cuw0HKoC0OW6WUtwGtlwNgQf9RCHEvgC+llJ+2X+4dAPYCgJTyJSHEXe0X/jGAO6SUvxRC/I0QYoeU8ss+P0ssCukDKeVL2vE7+vsIuVFkHzwipXy33w9QEIrsh6eklI8IIZ4YxfnQPoa8DhtSyvSv4H8ADlq+qwH4A4A7ADxh+f01APcax/960M8yqD5AS8jvAPDaoJ9lgH1wEMCOQT/HIPsBwBMA/nqt90MRMoG6DvdfMuMUjLY2ctr8Xkq5AODXaL3AO41zdgC4KLu1uL1oDYQ1h4L6YC1pcD0oqA+2ArgohPh12e0tCwX0w20Avo92Pwghan1odqEoSiZQ1+EiCfvi8WPbl+0t27uytR1bMOyPB6WUP9OPVwN9rdkp28jVB0KIe+XaNV8o5B4HUsqX2gLBPG4toYj58Id2P3yClqa/1lCITKCuw0US9gXCo3XcIaX8tP3382hpbRBCHJRS/rL99x1CiF8oZw1atrmtpTW4BBTRB2hpcfe2B/+O9ndrBgWNgyfW2nObKGgs/JN2Tg2x9uoBoaA+8F2HhSTsi8VWABf1L4QQ77T/fKk9ge8F8KhsOWDuBfAL5WVvn/9rAF+2f6vJtrNyDSF3H0gpP21r9lvRmuBrDUWMg1fb5x0EgDzb9wGiiLFwGkCt/RtGcT5Q1wlFyqBNSEhIGAEkzT4hISFhBJCEfUJCQsIIIAn7hISEhBFAEvYJCQkJawiaU9f6mUIS9gkJFogWN4k39FFLc09I6Bc+RneUmvnZiiTsExLs2KrFQJNoh4ju7UN7EtY52grGO5bvOkRpeZCEfUJCQsIQwMyl0HIs3m1/9ppqXEjCPiHBgbZW9U77/7/WNK01ydOSMDi0M6Pv0D4/4RlDdwFQ/FBfokWYBrQI4u7SzjU/W5EojhMSHJBSviuE+HX7/y/RYiL9cZuoai+Atc7hk9AnyBaV8UEhBNAeO23OHwqm8P5++zq/NK77SzCQNPuEBD8+tfydK3U9YTTRNtX8rP23j9W1UG6sJOwTEhIS+oS2Hf7X7b99hXn+Cava/Q4A7ziO9SIJ+4QEB9o21jvaE/Ne7e8fIyflbMJoQVWZahP9vQTgXt3O3nbA7jXI73ZopIi5TIaJCC0hwQKhlUcs8tiEhEEhafYJCQkJI4Ak7BMS7LjIzKA9iFYGY0LCUCOZcRISEhJGAEmzT0hISBgBJGGfkJCQMAJIwj4hISFhBJCEfUJCQsIIIAn7hISEhBFAEvYJCQkJI4D/H6VxG3DXswWuAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "from amuse.plot import scatter\n",
    "def plot(moon, disk):\n",
    "    scatter(disk.x, disk.y)\n",
    "    scatter(moon.x, moon.y, c='r')\n",
    "plot(moon, disk)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<amuse.datamodel.particles.ParticlesSubset at 0x7f4ce7eb48d0>"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from amuse.community.huayno.interface import Huayno\n",
    "from amuse.community.fi.interface import Fi\n",
    "gravity = Huayno(converter)\n",
    "gravity.particles.add_particles(bodies-moon)\n",
    "channel = {\"from stars\": bodies.new_channel_to(gravity.particles),\n",
    "            \"to_stars\": gravity.particles.new_channel_to(bodies)}\n",
    "\n",
    "hydro = Fi(converter, mode=\"openmp\")\n",
    "hydro.parameters.use_hydro_flag = True\n",
    "hydro.parameters.radiation_flag = False\n",
    "hydro.parameters.gamma = 1\n",
    "hydro.parameters.isothermal_flag = True\n",
    "hydro.parameters.integrate_entropy_flag = False\n",
    "hydro.parameters.timestep = 0.01*Pinner \n",
    "hydro.parameters.verbosity = 0\n",
    "hydro.parameters.eps_is_h_flag = False    # h_smooth is constant\n",
    "eps = 0.1 | units.au\n",
    "hydro.parameters.gas_epsilon = eps\n",
    "hydro.parameters.sph_h_const = eps\n",
    "\n",
    "hydro.particles.add_particles(disk)\n",
    "hydro.dm_particles.add_particles(moon.as_set())\n",
    "channel.update({\"from_disk\": disk.new_channel_to(hydro.particles)})\n",
    "channel.update({\"to_disk\": hydro.particles.new_channel_to(disk)})\n",
    "channel.update({\"from_moon\": moon.new_channel_to(hydro.dm_particles)})\n",
    "channel.update({\"to_moon\": hydro.dm_particles.new_channel_to(moon)})\n",
    "\n",
    "bodies.add_particles(disk)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [],
   "source": [
    "from amuse.couple import bridge\n",
    "from amuse.ext.composition_methods import *\n",
    "gravhydro = bridge.Bridge(use_threading=False) #, method=SPLIT_4TH_S_M4)\n",
    "gravhydro.add_system(gravity, (hydro,))\n",
    "gravhydro.add_system(hydro, (gravity,))\n",
    "gravhydro.timestep = 0.2*Pinner"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 23,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Planet: ae= 5.2 AU 1.05367121277e-08\n",
      "Moon: ae= 0.0355011649073 AU 0.0\n",
      "Time: 4.382906388 day dE= 1.0\n",
      "S=                  key         mass         name       radius  semimajor_axis            u           vx           vy           vz            x            y            z\n",
      "                   -         MSun         none  (0.001 * m**3)**0.3333333333333333  149597870691.0 * m  m**2 * s**-2  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  149597870691.0 * m  149597870691.0 * m  149597870691.0 * m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.000e+00         star    0.000e+00    0.000e+00    0.000e+00    2.187e-11   -3.465e-09   -2.847e-22   -5.009e-03   -3.185e-05   -6.218e-19\n",
      " 5326756165488562113    9.546e-04       planet    0.000e+00    0.000e+00    0.000e+00   -1.809e-08    3.582e-06   -7.995e-16    5.195e+00    3.291e-02   -1.641e-12\n",
      "10820997325257666469    9.546e-06         moon    2.175e+08    3.550e-02    0.000e+00   -4.779e-07    4.847e-06   -7.730e-12    5.228e+00    4.504e-02   -1.749e-08\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "g=                  key         mass       radius           vx           vy           vz            x            y            z\n",
      "                   -           kg            m    m * s**-1    m * s**-1    m * s**-1            m            m            m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.989e+30    0.000e+00    7.975e-02   -1.264e+01   -1.038e-12   -7.494e+08   -4.765e+06   -9.303e-08\n",
      " 5326756165488562113    1.899e+27    0.000e+00   -6.595e+01    1.306e+04   -2.915e-06    7.771e+11    4.924e+09   -2.454e-01\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "[-3.18507742038e-05, 0.032913852458] au [0.0450402518118] au\n",
      "Planet: ae= 5.20231704673 AU 0.00135812516381\n",
      "Moon: ae= 0.0355087471562 AU 0.000155418305473\n",
      "Time: 8.765812776 day dE= 1.00000783717\n",
      "S=                  key         mass         name       radius  semimajor_axis            u           vx           vy           vz            x            y            z\n",
      "                   -         MSun         none  (0.001 * m**3)**0.3333333333333333  149597870691.0 * m  m**2 * s**-2  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  149597870691.0 * m  149597870691.0 * m  149597870691.0 * m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.000e+00         star    0.000e+00    0.000e+00    0.000e+00    4.398e-11   -3.465e-09   -1.483e-21   -5.009e-03   -6.405e-05   -8.375e-18\n",
      " 5326756165488562113    9.546e-04       planet    0.000e+00    0.000e+00    0.000e+00   -3.693e-08    3.584e-06   -4.413e-15    5.195e+00    6.620e-02   -2.459e-11\n",
      "10820997325257666469    9.546e-06         moon    2.175e+08    3.550e-02    0.000e+00   -9.047e-07    4.612e-06   -1.070e-13    5.222e+00    8.909e-02   -3.041e-08\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "g=                  key         mass       radius           vx           vy           vz            x            y            z\n",
      "                   -           kg            m    m * s**-1    m * s**-1    m * s**-1            m            m            m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.989e+30    0.000e+00    1.604e-01   -1.263e+01   -5.409e-12   -7.493e+08   -9.581e+06   -1.253e-06\n",
      " 5326756165488562113    1.899e+27    0.000e+00   -1.347e+02    1.307e+04   -1.609e-05    7.771e+11    9.903e+09   -3.678e+00\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "[-6.40464486625e-05, 0.066198986002] au [0.0890939822376] au\n",
      "Planet: ae= 5.20910428051 AU 0.00298789843928\n",
      "Moon: ae= 0.0354834963725 AU 0.000619501933722\n",
      "Time: 13.148719164000001 day dE= 1.00003087729\n",
      "S=                  key         mass         name       radius  semimajor_axis            u           vx           vy           vz            x            y            z\n",
      "                   -         MSun         none  (0.001 * m**3)**0.3333333333333333  149597870691.0 * m  m**2 * s**-2  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  149597870691.0 * m  149597870691.0 * m  149597870691.0 * m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.000e+00         star    0.000e+00    0.000e+00    0.000e+00    6.585e-11   -3.465e-09   -3.201e-21   -5.008e-03   -9.589e-05   -2.941e-17\n",
      " 5326756165488562113    9.546e-04       planet    0.000e+00    0.000e+00    0.000e+00   -5.662e-08    3.587e-06   -9.059e-15    5.194e+00    9.915e-02   -8.566e-11\n",
      "10820997325257666469    9.546e-06         moon    2.175e+08    3.550e-02    0.000e+00   -1.223e-06    4.257e-06   -1.164e-11    5.212e+00    1.299e-01   -5.245e-08\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "g=                  key         mass       radius           vx           vy           vz            x            y            z\n",
      "                   -           kg            m    m * s**-1    m * s**-1    m * s**-1            m            m            m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.989e+30    0.000e+00    2.401e-01   -1.263e+01   -1.167e-11   -7.492e+08   -1.435e+07   -4.400e-06\n",
      " 5326756165488562113    1.899e+27    0.000e+00   -2.065e+02    1.308e+04   -3.303e-05    7.770e+11    1.483e+10   -1.281e+01\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "[-9.58933535063e-05, 0.0991511806367] au [0.12993645851] au\n",
      "Planet: ae= 5.21945835485 AU 0.00495209452991\n",
      "Moon: ae= 0.0354797983657 AU 0.000765016152061\n",
      "Time: 17.531625552 day dE= 1.00006626056\n",
      "S=                  key         mass         name       radius  semimajor_axis            u           vx           vy           vz            x            y            z\n",
      "                   -         MSun         none  (0.001 * m**3)**0.3333333333333333  149597870691.0 * m  m**2 * s**-2  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  149597870691.0 * m  149597870691.0 * m  149597870691.0 * m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.000e+00         star    0.000e+00    0.000e+00    0.000e+00    8.796e-11   -3.464e-09   -5.573e-21   -5.008e-03   -1.281e-04   -6.955e-17\n",
      " 5326756165488562113    9.546e-04       planet    0.000e+00    0.000e+00    0.000e+00   -7.795e-08    3.591e-06   -1.467e-14    5.194e+00    1.325e-01   -1.950e-10\n",
      "10820997325257666469    9.546e-06         moon    2.175e+08    3.550e-02    0.000e+00   -1.404e-06    3.818e-06   -6.901e-12    5.200e+00    1.675e-01   -5.556e-08\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "g=                  key         mass       radius           vx           vy           vz            x            y            z\n",
      "                   -           kg            m    m * s**-1    m * s**-1    m * s**-1            m            m            m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.989e+30    0.000e+00    3.207e-01   -1.263e+01   -2.032e-11   -7.491e+08   -1.916e+07   -1.041e-05\n",
      " 5326756165488562113    1.899e+27    0.000e+00   -2.842e+02    1.309e+04   -5.349e-05    7.770e+11    1.982e+10   -2.917e+01\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "[-0.000128082504309, 0.132495141349] au [0.167489722635] au\n",
      "Planet: ae= 5.2324123078 AU 0.00721727044567\n",
      "Moon: ae= 0.035459166036 AU 0.00125723157926\n",
      "Time: 21.914531940000003 day dE= 1.00011087891\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "Exception ignored in: <function CodeInterface.__del__ at 0x7f4ceb284f80>\n",
      "Traceback (most recent call last):\n",
      "  File \"/home/spz/AMUSE/amuse/src/amuse/rfi/core.py\", line 755, in __del__\n",
      "    self._stop()\n",
      "  File \"/home/spz/AMUSE/amuse/src/amuse/rfi/core.py\", line 845, in _stop\n",
      "    self._stop_worker()\n",
      "  File \"/home/spz/AMUSE/amuse/src/amuse/rfi/core.py\", line 123, in __call__\n",
      "    raise exceptions.CodeException(\"Exception when calling function '{0}', of code '{1}', exception was '{2}'\".format(self.specification.name, type(self.interface).__name__, ex))\n",
      "amuse.support.exceptions.CodeException: Exception when calling function '_stop_worker', of code 'HuaynoInterface', exception was '[Errno 32] Broken pipe'\n",
      "Exception ignored in: <function CodeInterface.__del__ at 0x7f4ceb284f80>\n",
      "Traceback (most recent call last):\n",
      "  File \"/home/spz/AMUSE/amuse/src/amuse/rfi/core.py\", line 755, in __del__\n",
      "    self._stop()\n",
      "  File \"/home/spz/AMUSE/amuse/src/amuse/rfi/core.py\", line 845, in _stop\n",
      "    self._stop_worker()\n",
      "  File \"/home/spz/AMUSE/amuse/src/amuse/rfi/core.py\", line 123, in __call__\n",
      "    raise exceptions.CodeException(\"Exception when calling function '{0}', of code '{1}', exception was '{2}'\".format(self.specification.name, type(self.interface).__name__, ex))\n",
      "amuse.support.exceptions.CodeException: Exception when calling function '_stop_worker', of code 'FiInterface', exception was '[Errno 32] Broken pipe'\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "S=                  key         mass         name       radius  semimajor_axis            u           vx           vy           vz            x            y            z\n",
      "                   -         MSun         none  (0.001 * m**3)**0.3333333333333333  149597870691.0 * m  m**2 * s**-2  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  3646245880.3572216 * m * s**-1  149597870691.0 * m  149597870691.0 * m  149597870691.0 * m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.000e+00         star    0.000e+00    0.000e+00    0.000e+00    1.098e-10   -3.464e-09   -8.781e-21   -5.007e-03   -1.599e-04   -1.349e-16\n",
      " 5326756165488562113    9.546e-04       planet    0.000e+00    0.000e+00    0.000e+00   -1.006e-07    3.595e-06   -2.105e-14    5.193e+00    1.655e-01   -3.587e-10\n",
      "10820997325257666469    9.546e-06         moon    2.175e+08    3.550e-02    0.000e+00   -1.425e-06    3.354e-06   -2.367e-11    5.186e+00    2.004e-01   -5.223e-08\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "g=                  key         mass       radius           vx           vy           vz            x            y            z\n",
      "                   -           kg            m    m * s**-1    m * s**-1    m * s**-1            m            m            m\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "11992421139574210381    1.989e+30    0.000e+00    4.004e-01   -1.263e+01   -3.202e-11   -7.490e+08   -2.392e+07   -2.017e-05\n",
      " 5326756165488562113    1.899e+27    0.000e+00   -3.670e+02    1.311e+04   -7.677e-05    7.768e+11    2.476e+10   -5.366e+01\n",
      "====================  ===========  ===========  ===========  ===========  ===========  ===========  ===========  ===========\n",
      "[-0.000159920371663, 0.165517521875] au [0.200445134594] au\n",
      "Planet: ae= 5.24621757939 AU 0.00954797115782\n",
      "Moon: ae= 0.0354518195719 AU 0.00148205052137\n",
      "Time: 26.297438328000002 day dE= 1.00015896666\n"
     ]
    },
    {
     "ename": "KeyboardInterrupt",
     "evalue": "",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m                         Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-23-56f3cbf30150>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m     35\u001b[0m \u001b[0mt_end\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;36m1.0\u001b[0m \u001b[0;34m|\u001b[0m \u001b[0munits\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0myr\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     36\u001b[0m gravity_hydro_bridge(gravity, hydro, gravhydro, \n\u001b[0;32m---> 37\u001b[0;31m                      bodies, t_end)\n\u001b[0m",
      "\u001b[0;32m<ipython-input-23-56f3cbf30150>\u001b[0m in \u001b[0;36mgravity_hydro_bridge\u001b[0;34m(gravity, hydro, gravhydro, bodies, t_end)\u001b[0m\n\u001b[1;32m     21\u001b[0m               \"dE=\", dE_gravity)#, dE_hydro\n\u001b[1;32m     22\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 23\u001b[0;31m         \u001b[0mgravhydro\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mevolve_model\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mmodel_time\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     24\u001b[0m         \u001b[0mchannel\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m\"to_stars\"\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcopy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     25\u001b[0m         \u001b[0mchannel\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0;34m\"to_disk\"\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcopy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/couple/bridge.py\u001b[0m in \u001b[0;36mevolve_model\u001b[0;34m(self, tend, timestep)\u001b[0m\n\u001b[1;32m    554\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    555\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmethod\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 556\u001b[0;31m             \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mevolve_joined_leapfrog\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtend\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mtimestep\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    557\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    558\u001b[0m             \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mevolve_simple_steps\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtend\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mtimestep\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/couple/bridge.py\u001b[0m in \u001b[0;36mevolve_joined_leapfrog\u001b[0;34m(self, tend, timestep)\u001b[0m\n\u001b[1;32m    574\u001b[0m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mkick_codes\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtimestep\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    575\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 576\u001b[0;31m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdrift_codes\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mtime\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mtimestep\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    577\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    578\u001b[0m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchannels\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcopy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/couple/bridge.py\u001b[0m in \u001b[0;36mdrift_codes\u001b[0;34m(self, tend)\u001b[0m\n\u001b[1;32m    703\u001b[0m         \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    704\u001b[0m             \u001b[0;32mfor\u001b[0m \u001b[0mx\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mthreads\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 705\u001b[0;31m                 \u001b[0mx\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrun\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    706\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    707\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0mkick_codes\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0mdt\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/anaconda3/lib/python3.7/threading.py\u001b[0m in \u001b[0;36mrun\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m    868\u001b[0m         \u001b[0;32mtry\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    869\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_target\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 870\u001b[0;31m                 \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_target\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_args\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_kwargs\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    871\u001b[0m         \u001b[0;32mfinally\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    872\u001b[0m             \u001b[0;31m# Avoid a refcycle if the thread is running a function with\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/couple/bridge.py\u001b[0m in \u001b[0;36mdrift\u001b[0;34m(self, tend)\u001b[0m\n\u001b[1;32m    427\u001b[0m             \u001b[0mprint\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcode\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__class__\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__name__\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m\"is evolving to\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtend\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    428\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 429\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcode\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mevolve_model\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mtend\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    430\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    431\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/support/methods.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *list_arguments, **keyword_arguments)\u001b[0m\n\u001b[1;32m    166\u001b[0m         \u001b[0mobject\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mprecall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    167\u001b[0m         \u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeyword_arguments\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconvert_arguments\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 168\u001b[0;31m         \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmethod\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    169\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    170\u001b[0m         \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconvert_result\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mresult\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/support/methods.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *list_arguments, **keyword_arguments)\u001b[0m\n\u001b[1;32m    166\u001b[0m         \u001b[0mobject\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mprecall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    167\u001b[0m         \u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeyword_arguments\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconvert_arguments\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 168\u001b[0;31m         \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmethod\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    169\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    170\u001b[0m         \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconvert_result\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mresult\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/support/methods.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *list_arguments, **keyword_arguments)\u001b[0m\n\u001b[1;32m    166\u001b[0m         \u001b[0mobject\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mprecall\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    167\u001b[0m         \u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeyword_arguments\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconvert_arguments\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 168\u001b[0;31m         \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmethod\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    169\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    170\u001b[0m         \u001b[0mresult\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mconvert_result\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mresult\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/support/methods.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *list_arguments, **keyword_arguments)\u001b[0m\n\u001b[1;32m    265\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    266\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0m__call__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 267\u001b[0;31m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mmethod\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m*\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    268\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    269\u001b[0m     \u001b[0;32mdef\u001b[0m \u001b[0masynchronous\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m*\u001b[0m\u001b[0mlist_arguments\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m**\u001b[0m\u001b[0mkeyword_arguments\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/rfi/core.py\u001b[0m in \u001b[0;36m__call__\u001b[0;34m(self, *arguments_list, **keyword_arguments)\u001b[0m\n\u001b[1;32m    118\u001b[0m             \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minterface\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchannel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msend_message\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcall_id\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mspecification\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdtype_to_arguments\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdtype_to_values\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    119\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 120\u001b[0;31m             \u001b[0mdtype_to_result\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minterface\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mchannel\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrecv_message\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcall_id\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mspecification\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mid\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mhandle_as_array\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    121\u001b[0m         \u001b[0;32mexcept\u001b[0m \u001b[0mException\u001b[0m \u001b[0;32mas\u001b[0m \u001b[0mex\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    122\u001b[0m             \u001b[0mCODE_LOG\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minfo\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m\"Exception when calling function '{0}', of code '{1}', exception was '{2}'\"\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mformat\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mspecification\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mname\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mtype\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0minterface\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__name__\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mex\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/rfi/channel.py\u001b[0m in \u001b[0;36mrecv_message\u001b[0;34m(self, call_id, function_id, handle_as_array, has_units)\u001b[0m\n\u001b[1;32m   1945\u001b[0m         \u001b[0mmessage\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mSocketMessage\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1946\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1947\u001b[0;31m         \u001b[0mmessage\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mreceive\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0msocket\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1948\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1949\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0mmessage\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0merror\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/rfi/channel.py\u001b[0m in \u001b[0;36mreceive\u001b[0;34m(self, socket)\u001b[0m\n\u001b[1;32m   1496\u001b[0m         \u001b[0;31m# logger.debug(\"receiving message\")\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1497\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1498\u001b[0;31m         \u001b[0mheader_bytes\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m_receive_all\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;36m44\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0msocket\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1499\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1500\u001b[0m         \u001b[0mflags\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnumpy\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mfrombuffer\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mheader_bytes\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdtype\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m\"b\"\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mcount\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m4\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0moffset\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/AMUSE/amuse/src/amuse/rfi/channel.py\u001b[0m in \u001b[0;36m_receive_all\u001b[0;34m(self, nbytes, thesocket)\u001b[0m\n\u001b[1;32m   1478\u001b[0m         \u001b[0;32mwhile\u001b[0m \u001b[0mnbytes\u001b[0m \u001b[0;34m>\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1479\u001b[0m             \u001b[0mchunk\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mmin\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mnbytes\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;36m10240\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m-> 1480\u001b[0;31m             \u001b[0mdata_bytes\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mthesocket\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrecv\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mchunk\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m   1481\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m   1482\u001b[0m             \u001b[0;32mif\u001b[0m \u001b[0mlen\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mdata_bytes\u001b[0m\u001b[0;34m)\u001b[0m \u001b[0;34m==\u001b[0m \u001b[0;36m0\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mKeyboardInterrupt\u001b[0m: "
     ]
    }
   ],
   "source": [
    "from amuse.ext.composition_methods import *\n",
    "from amuse.ext.orbital_elements import orbital_elements_from_binary\n",
    "\n",
    "def gravity_hydro_bridge(gravity, hydro, gravhydro, bodies,\n",
    "                         t_end):\n",
    "\n",
    "    gravity_initial_total_energy = gravity.get_total_energy() + hydro.get_total_energy()\n",
    "    model_time = 0 | units.Myr\n",
    "    dt = 0.012|units.yr  #1.0*Pinner\n",
    "    while model_time < t_end:\n",
    "\n",
    "        model_time += dt\n",
    "\n",
    "        orbit_planet = orbital_elements_from_binary(bodies[:2], G=constants.G)\n",
    "        orbit_moon = orbital_elements_from_binary(bodies[1:3], G=constants.G)\n",
    "        print(\"Planet:\", \"ae=\", orbit_planet[2].in_(units.AU), orbit_planet[3])\n",
    "        print(\"Moon:\", \"ae=\", orbit_moon[2].in_(units.AU), orbit_moon[3])\n",
    "        \n",
    "        dE_gravity = gravity_initial_total_energy/(gravity.get_total_energy()+hydro.get_total_energy())\n",
    "        print(\"Time:\", model_time.in_(units.day), \\\n",
    "              \"dE=\", dE_gravity)#, dE_hydro\n",
    "\n",
    "        gravhydro.evolve_model(model_time)\n",
    "        channel[\"to_stars\"].copy()\n",
    "        channel[\"to_disk\"].copy()\n",
    "        channel[\"to_moon\"].copy()\n",
    "        \n",
    "        print(\"S=\", bodies[:3])\n",
    "        print(\"g=\", gravity.particles)\n",
    "        print(gravity.particles.y.in_(units.au), moon.y.in_(units.au))\n",
    "\n",
    "    gravity.stop()\n",
    "    hydro.stop()\n",
    "\n",
    "t_end = 1.0 | units.yr\n",
    "gravity_hydro_bridge(gravity, hydro, gravhydro, \n",
    "                     bodies, t_end)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEMCAYAAADNtWEcAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjMsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+AADFEAAAgAElEQVR4nO29XZCUVbrv+V+V9QJZuDcJBhHT5og2ngiYMGiqEKed7YVgoMSMW6MGUaLDvj06l7CNUpx2GuixDyjRobc6t+10FGDvCtnsMyijeOEO3Q1UlgRxNGaLiic1JoyGZLZUIllVay4yV9bKletZH+9HflSuX4TdxZuZ77ver+dZ6/lknHMEAoFAYPAY6vYAAoFAINAdggIIBAKBASUogEAgEBhQggIIBAKBASUogEAgEBhQhrs9gE7DGHsWwGXO+RnGWAHAjsZHZzjnlS4OLRAIBDrKIK4AzgEoNP5+iXN+AsAZAM92b0iBQCDQeZaMAmCM7WaMva/ZtqMx69exXvr7nuxGFwgEAr3HklEAjZl8E8bY7sb2M41/79D87LL095fZjS4QCAR6j75RAIyxLYyxLdK/n23Y8Cnux6KAvwxA/HYHgPsbv32zoSh2AHgrg2EHAoFAz9I3TmDO+YWGSQcAtsLutFWVw+2N/bwmbaugdRUQCAQCA0PfrACAppnnucbfNsFdAbAm80EFAoFAn9JXCqBhrnmz8fd6y9f/gsVVwHoA7xu+GwgEAgNH3yiAhv2/wjm/wDl/C8AO2QfQcPJulZy/JwCsb2wvCGdwIBAIBOqwUA46EAgEBpO+cAIzxqYBrAXwb90eSyAQCPQR/wHAD5zzMd2HfaEAAKxdtWpVcXR0tNjtgQQCgUC/UCqVcP36dfLzflEA/zY6Olo8e/Zst8cRCAQCfcO2bdvw0UcfkZaTvnECBwKBQCBdggIIBAKBASUogEAgEBhQggIIBAKBASUzJ7BUffMRzvmLms93o16uYT3n/K1GUtf6xn/3634TCAQCgfTIZAXQyNrd0si+3aKWbSBKNT8NYKso62yo4R8ILAmmpst48MgH+Pn+U3jwyAeYmi53e0iBASOTFQDn/AKAC41Z/WVN4bb7AUw2/r6MurKQq3S21e4plUrYtm2b9nghPDTQb0xNl/HSny+iWpsHAJQrVbz054sAgPGxkO4S6AxZ5wFshb7RirZUM9As8nY11O4JZMnUdBlHT3+B7ypV3FHIY2LnhkSCl9oftf3o6S+awl9Qrc3j6OkvggIIdIxMFUCj8fpTjLHdSscuU6nm3Zzz59SNo6OjYaYfSIWp6TImTsygNl+vg1WuVDFxYgZAvNm3bja/b7KEvZOllu/Js/zvKlXtvqjtgUAWZOUDeFWy4euEvbZUc0NRvNb4ewsCgQw4dPJSU/gLavMceydLuHv/KYweeq9pj3ex0x9891LbbJ4qsShm+XcU8trPqe2BQBZkFQb6JoDLUinmtwBANG3XlWpu/P0qY+w8Y+w8QjOXQEZcm60ZP69Ua5g4PoOXpy5i4sQMypUqOBZXCrISmJouo1I170/lu0oVEzs3IB/lWrbnoxwmdm7w2lcgkISsnMCXsdhq8Yy0/RHp79fkzxs2/3uyGE+g+6Rtc8/6WLUFjv/z0ytYUKbytXmOQycvNfd39PQX3uO7o5Bv+X0nrkkgoKNfisEF+pisI15kgV8YifDjzTnUFhbt++qxCvnIadauCn+BvILwtdkzoDnLHx8rdlzgm5RjJ5V0oDcImcCBzDFFvCRFKBdhprk2W2sKf+pYB5+4F9EQS3xswM9mzwA888C6rglV9VoJ5Tg1XdZ+tm+yhJenLnZlrIHOEFYAgVTRzSKzjHjRKRfbscbHijj3zVX86dNvMU90xIuGGIZzDNXaQttnhXzU/Hti54aW1Q0ARDmGlcuGUanWkGMM85yj2AMzapsi1jmy//jJFfzxkys9Mf5A+gQFEEgNytSzijC5pBHx4qpE8tEQHjzyAb6rVLEqH+HGrTlS+BfyEQ4+cS8AYOL4TMuKIhpizc8AaG352zeuxYef/4Dr1Rr+m1UrekZwlhMo4pCotjQJCiCQGtQMc0U0hHyUa/ksrYiXOwp5UrDJzNYWMNv4HmX/Lxby+Hj/w23bbXZx2Zafhb8jDdu8qcxEPhrC6pXLrdcxJKotPYICCHihE0ZAXUhSAqQyW8Pre0YzcTDqTDBx0c2EfR21aWf4pqVQTP6W2doCdm1ci7c/uULmLwhCotrSIiiAgDM6YTRxfAZgaEuskimMRJlFl6gmmKGGzT0OWZqk4gjOqekynj8203Y+st3e9brajv+nT7+1Cn8gJKotNYICCDijm92qETcqUY7hx5tzzdDJLGzJJhOMD9s3rk08Fsok5Ss4xXlQykxcR9eVgc1UJh+Hoe4AFv8vCIlqS48QBhpwYmq67GRrl1k9EmHlsmEyLDOLcsjjY0Uc3rUJxUIeDHW7fj5ye8w//PyHxMdPK8PXFt2UY8wrtHZi5wZEObfQV476dXt9z2jLdTy8a1Ow/y8xwgogQCLs/eVKFXGi5keWDZOmB98ZrA+q3f7n+085/c5kJnF1xPpk+L48dbEZippjDL/65Z14ZXyTdSyqQ93lHMTxf/OPF3HjllvYbDcS1QKdJSiAQAuq0Bdzd8rQEw0x0gwkBKBu5WCawVJCJ240jGukEGWm8XXEugjOl6cu4o+fXGn+e57z5r9fGd9kvG6Hd20ine4mU5MYl3wdKZ9JsPUPBsEEFGgiZ4MCtNCXOfrU5pbEKBkhpHUmEcq2LWawqnno5amLZBarDd0YVExmmiwymd/+9IpxO3Xd/vD0ZoyPFROZmsbHivh4/8P46shj+MPTm0NRugEmrAACTVyzagVFqaiZ6ngVQoRKlKKycO8o5LUzbl2Iomt4pSlZK0kETZKQSCpQSWy3mZLSKibnsp9QI2jpEhRAoImPk1cIeCEcqrV5p7IHs7fmMPmveuEfDTFM7NygVUTUasRVCOvMMrJgE7N5nwiarM0kNlNSWjZ6035C68qlTTABBQDUX3Sbo1d8LiJCALSYjOY5b5v5uxRrUw/gM7MWQtg3oshUGE0li9r9I0RkErW9W2RZyC/QfcIKIACg3tXKZPPXzeofPPKB1ZHrY1aqzfNmtyzdjJuKS6daMp775mozqkbFJ2M3i9r9/2nXL/APx0otJaeHWH17XHz7ErvsI0n9oKQE01P2BAUQsHa1+vrIY9rtLrZxX0HxXaWK1/eMtvkUGIC/u2cNvv5rtU0g6BQRB/D2J1ew9a41WqFhCk998MgHbcIm7ZBIV6XiKgQpU825b67infNlJxOObh+q0hVkbf4KpqfOEBTAgKETKKblfNHworvYxl1DMOXvj48VcfzcFXz85dXmdg7gwpXr2mQkSphzgHQSm8blI2ySzFJtSsVHCFIrGp2zXbfSocpOdCsjmDqfQycvhVVBivSWwTGQKZTd2ySgZ2/NkfZ0F9u4Swim+tup6TL+RRL+Asr2bJqNUsrBNi4XO7fueu6dLGHsd+85hafa8LG/U+dpC7cF7GUnRGZwJzOCqfO5NluLFQoc0BNWAEsEU1apgBIoOUMBtWuzNXLW6WLGEH8ffPdSm5lJNE65Xq21mXR8on4mdm7A3smS9vuUcpDHHtfOTfk3qGvmu1rwCT81JY7ZEr1sfprVI5G2THaWuK4cq7V5PH9sBvsmS2FFEIOwAlgCiKxS8aKLrFK1nZ9plph0NmxifKyI0oFH8YZSW+bo7s0oHXgUXx15DB/vf7j54poELyXQdR0eoxwzmipEQhRl5rLZuU3jrNbmcfDdS81/+0Qd2Y4/xFhbtBO1GvvVL++0rtJsii5mcdVE+Kwc5zkPK4KYBAWwBPjTp986bacECgPw5H1Fo71fJyR8hZqcgSoLfBXTOHUC/ejpL7QN3FcuG3aaDcYN87QpiEq11rwWccIpKSGoE3i6IniHd23CK+ObtNvl62I7j+uGAIGskM/HhxCi6kdQAEsAynyjbp/YuUEb689Rr4TpMxsWTsMsYsR1gs/UUJ2awboKLkp42pSHyyxVXIs42cTquHKs/e7J15tSsDbFazsPl4ifrCq7frz/Ye9ChFQ5kbAyaCf4AJYAlJ1XFRjjY0XSVi5eGqrJ+Y2f5vDz/aeaZRTeOV92cjDGwTfunrIXryJqFFHHdLUdy7b8VfkIDByzmubxwOK1oPoi20pVy+Oiqpq6Xm/KB2Hy07ishLIO2YwTSRbCSN0IK4AlwK9+eafzdtsMX511rh6JAF43Zwizw9ufXDE6DdOIEXc1FwGNWvcaJ8ANQwRTXFSzV/26MKxcpp9Bi2uhmbwDqLdjdB0jdV1dZ+gmcx3lp3FZCWWdLUyZ6H79wDrSdBcymN0ICmAJ8Mr4Jvz6gXXNGX+OMfz6gXXaLFgXe7csfEc0DV1MPsEsY8SpJf34WBG3rWhfzIrM4jShBEuUGzJe18osbY5yHWOSkhSuAtFH8QqyKJanjklnott61xosH14UYatHoqbCynpMS4VgAuoTbCGEr4xvIsseyPiaV3xeGFGrPosltm1JTwnYtF94k7/B1Pg+jbIKSUpSZCkQO1EsTzXR6Vp/3pTMcCaz4INHPgiJZA2CAugD0rZn+ti7feyvf5sfxr7JEo6e/iL1F8tWu6dTFTtNxzFd14mdG7BvspS4rELckhRZXh+d34ghnR7LFNTzIHICCiNRW7OiaIjhxq25pp8j+AWCCagvSNue6RMd4ROPfW22llk8tm0GG9c8YrsW6ufbN66NdZzxsSKeeWBdW0RLp5qvZHV9gPq5PXlfseXcOIB3zpczi7wx5bToqs6uHolw24ph1Ob1ZTEGlbACyJg0KhqmuXyP094QWDQ7UC0EVVybtbhim8HGMY/Yqoi+PHWxpRFNuVLFO+fLePK+onMzGZlXxut26yxq2UxNl3Ho5CVca5jCCvkIB5+4N1EDmanpMiZOzDSFZrlSxcSJmZb9CT78/IfYDXvi4BsZdLO24N1HeRAICiBD0jLdpLl89ymDLJDNDjrbK0WaL5bOzKBzXsvlj23mKKrxzNuN3rxUFzKRM+FK1mWNVUEN1KOTJo63Cmtf89Ghk5faZsy1eY5DJy85F+TLSrjqngcTppIng9z/OCiADIkjbHW4CD9Xkr6oupnk7K255sxTJm0noHpcqnyyq9I1VRH906ffench0wl6AIknATYFcvT0F22CGgBqCzzWDHxquqzNCRBQ97qTXdPirEx1nw96/+OgADIkrVlRkugPlbgvqkkIUauCcqWKu/efwuqRCAcevzfxrNdlBuujdE1mBJMw0V0rSvEsHx6KPQnQCWKdAjE9T77P2tR0GRPHZ+iubQRpTlJcibsyFaT1XPYzQQFkSJqzorjRHyoTOze0mQtsRdNss2pbZc1rszXSdpwWQkH5hFqaInNMFVLlrGihCCnFQwmkcqVqVKqq/0Hdr6xATIrM91k7evoLq/AvaDKs05ykxEF9BqlGNjIjjrWiljJBAWRIN2ZFTvZm9c2wvCkus2qhCB488oFWGImkLJdZb0uZBVZPotKtOuTv3bg1pzWDCHSCcHysiHPfXG0TtPkohyfvK7Z00hIMMWjDCH0ckgK5LIfcweufZr43dmgDWhWaTqkL7r497xX3blsxREMMB5+4V/tZnElKmv4R+fhyeXSKQXb+CjJTAIyxHY0/H+Gcv6j5fDeACoD1nPO3pG3Pcc4fyWpcnUSdleQYayvelSYuvXF1MzybrdjHlGUzR/iYkijTBwDyezpMStcUmaNu1/k6xP00rRhcqdbm8ceGA9qGrNDE9ftf//xZW00iuauaeB72Tpa0PZ7FfilllmMMR5/anNpzm1W9nqnpsrFWlWCQnb+CTBQAY2wLgC2c89cYYy8yxtZzzi9Ln+8GAM75GcbYs4yxHZzzM5zzE4yx57IYU7cQD3InClOZolpEb1xTG0QVIayp10j3ApkESGEkMl4HW2MSWXm62nopQSdDzVzV7aZibJ0uma8mWYnrN2uZ1cohrbpncGLnBq0PIMoxHN2dnvAH0guScNmvyqA7fwWZKADO+QUAFxhjBQCXZeHf4H4Ak42/LwPYAuCMaZ+lUgnbtm3Tfnb27Nkkw82crB50lTi9cSlsTjX1BbLZ4KMcA+ftglu+Di5Lch9TS7GQJ8M145gebD6dOGaguHz4+Q9t23xNGrpnUPwtO5+zcpZmFTpq+j0DjPc765DdXiNrH8BWAF9qtheUf9+e8Ti6SqdipNOoNyMwzaLUWbVNWQgBss9SitolucfFuQeYZ3hxTQ82B3qc6Jm4+LSF9N1PWgEHNrIKHaX2a5oQANmXte5FMi0FwTk/A+AeYfKRqABY47Ov0dFRnD17Vvtfr5OkjK8PVMMXYLEIlgk53Z9SGAxoqxJJKYtiIY+vjzyG6d8+2qzXo0Nst5WdcBX+clVIHYlKaxAO9PGxIo4+tVkbIZMFumvpU7ZD4NMzIW2SVDfNYr+DWEI6Kx/AqwC+bDh3dcL+L1hcBawH8H4W4+gVOhUNREW1qEWwKOQ6PoWRSJvwUxhpr6bousKxXQc1lFCNAnKd3d4kmrNQ47JtF9gc6Gpc+vPHZrSOyEI+wsrlw81rKBrsuPo1oqH2sF1huvCJgwfoPgXqfuWxximDodvXxM4NOLxrU+oml7ghqdT9L1equOelf8Y8504+pX4iKxPQmwDWNyKBClKUz/uc80cazt4XpM/PND7fAWArY2w35/xERmPrOJ2MkRZRLXJdmDnOwc0ysYVqbR7Lh+v17VWBcm221txvuVLFxPEZMKZvHK7OUl2iokzmByrEVDd+k88jrunBR3GMjxVJk9f1ag2lA4+2bJMjjmTFVxiJ8FNtvhndo9b4AeIlQQlMfQp0JhE5SsnHREKZVw7v2uRVVsOVOGYsl8TApWYWysoJfBl15y4gOXfl8E7O+Wuaz88AWJ3FmDqJrfWe+E6WdcnlWXCc6ERR395UEgAAafOmVjhJoqJ86r+YZvNxV2Q2xaHed6oNJEddmcn3PK7d3bTScMGk9FxWFK7BDJ0KhEiC6/Mlj7vfncYhESxlXBxJWTub4pgCVIYatoGVy4etpiMVW2OYuMIgrTpEwlQmEoVyjOHJ+9oVtPpimxSH7p5GOdZWk16Qxj0Xx4wr/G1KzzVwwOV7ptVTWkI06X7U58KEGHe/O40ZT5i80gkYY2cfeuihh3rZ4WsLg5QjEChThi1KwZWf7z+VSly6zgTkSrGQJ19EanwMwFdHHvM6js78kY9yRgVEmUxEtBIAcp+A3pRH3dPVIxFGlg07PRcu5yof+8ZPdr+OiTf2jBoFlavJzeUcqH0V8hF+mlvwun864jwHLvugyDGGv1mhnxyl9R6nwbZt2/DRRx99xDnfpvs8rABSwOXBkWdAWYeFxgkH1GEqoWuCYTEmXjcrcrHBywpV9jGoMelx/CvUCunabM1awI3qk0vdu8psDdO/fZRUeuVKtc0cpEM32zSRj3JYEQ1pV0dAXUjZBKOpxITAVkdK3pdOQDNmzg1xxWdVSa0UfFbO85yTyrefSkwEBZACLg+OLNyyqMgpf+5aDMuFec7bVgJi3wVNDR7dcau1eew7VsKhk5dQma1hVT5ClGMtv5PNEaqwk/XPtdka/uFY3bka135uekFNBdxMv7PdU5NSdjEd+AgnYYID9LkJrkJ7fKxo9QGtdCyoRilqW26IK66TKpPZJs0JWL8QFEAK2B4c1dYaxwlJPbjnvrmKDz//oU3ocywK42KCFYEIezPV75E/o47D+WId+Uq1hmiIYfVIpC30dujkJaOwW+D1TNW4dta4KySqDDSldOV7evft5mNWa/PaRisCV+GkM3vINYIYA/bcf6fztbtuMTHZPpfRKWpKwfjmJ7hOqkwrBWofPqvgfisxERRACpgEii5uOC2zRbU23xLzr8tREvZIV3uuTD7KYfvGtcZxqi+163FqCxwjy4Yx/dvWcMip6TJptpBRhYYsiMULK/5fvQe2aI/VI1FbC0Hdi922UkGr0hXHnJou41+komwU12ZrmJoua58D6hkTPgbTc8Sl9EDO6716RW0oGzZlmXS2S+UhmPITfB30MqaVwut7Rkk/AuXfc7n+vUxQAClAPXwmB1RaZgvbvETYmKlkIyqGX0TGyL9xMVX4hGqKmvhAa2cnX1RBLGZrcuy2WgXz8K5N2tlnPsrhsV/8DKc++765P13sPVCfveqK76mYCurpvqu7ttQzZqvRkzT8cmLnhpay1brPgVahXBiJwHl9dWATilQeArVdtxLeO1nC6pHIqVezaaVgm5jFuf69TlAAKdCJRK8kjl25mfmpz75vmWHrhL886/EVHmL7vmMlp/yDieMzAEPTH+C61F49Elkjr2TUKpiHd21C6cCj2kxXVVFer9awt9FfWGTC2o4Z17ZMfTft7FbbmOTrMsTqZjeV1SNRc4UjC0f5+RIC+tDJS1ph6eMPM+U8XJut4Z3zZWvUj0s2uu73nUzm7CRBAaRE1gW0fJtgq1Rr9WbmI8uGtSYWeSWwfLheIipJYbnhIWaMHhHELZ722C9+Fvt6mJTYP15oXyXJysO1Xr98HB/lTSWV6ZIJdbgmo5lMN2qTeWqSIEJmXRzUIsIKaM2Hmb01p903ZW4zTRDi5pK4CvKs3/FuEBRAB0gj0UUto0Bhiv4xCW75vapUFyNtdNjsvkdP65uUp0UhH+HDz39IlOxGJfKkDWVbBtA2sxY+l9FD71n7AOtwTUazOSoPnbxkvH+qT8V1haNm0OquCWVuc42CchnLUhTkcQkKIANM7QqTZAuKB5dytOYYwx+e3kwqCZ+69dTEnAFG4QGkE05HKbJ8lMPBJ+jS0q7cUcinkjHtchxq1qlu275xLSb/9Vvtqqham8fzx8x9lXXnU5vn3o5KkxNehJDKv/dZ4Yhng7r2K5frw0pdn6lOhWD2ewkIQVAAKWNqayio1uabNuU4D47OHBTlGFYuG8a+yRIZZ79941qc+uz7mGdWh4MWQLYOYrrommiItfgAqG1A6+zQtBKyhe2JGXBSJWLDx7YMAKOH3jOaxOY5104ebL4QkYyWBrV53hau6mOeFALa1zfhomQ6FYK5FEpACDLtBzCI+MwqxYMjImFcGR8r4vCuTSgW8mCoC1bwurLhaCgdXt/OUF+yi4gedXZXyEf13ztSJGZY4qWgXlJhM5bHXchHuG3FMGrz9ZBNsX+xTUWeHVI139/YM4ovD/9PZF1+uU5RlrPFYiHvXc7ApayDWp9+arqMieMzqYZq2noaqM+Q7nnMR+2iRRbQvj0ydPdb5JKIZ9z3esclTt8AUfxR7rnRC4QVQMqk0ZZPRRfjLtthHzzyQdtLKeLsDzx+L46e/oJ0Xq5cPuw1g1N70Qp8OojpbMAi43j7xrXkWOVra3LmTU2XcUPjXIyGWpuau5y3WFlVqjXn7OpfP7AOr4zXM3GzMBXIwv7gu5eMq4Y4s+K/3/wzq7P75/tPtTmnXUsuAP7JkHGct1mZaXxXL728YhhIBZCl/S6ttnwCKsbdJcxQfMdWo0h9uajIEQA49dn32lhrWwcxFVNiG4VaRpkSOlSo4G0rWu3LOqFianjiWnpZ9OvVvfgTJ2Zw8N1L2hj51UQTHhXW2Pf4WNG4aojbvETXb1hFNA/aN1nCuW+uNhWejMnZGkeg+zhvsxS6vqVcerkU9sApgKy1sdY+P8Rw2wp9+CVgTns3zaxdUthts/pCw/yjvlx37z+l/b7aEEZcO9+XIklim3y/1ASkH2/OkQJal1zkI1SohCAVk6OzNr9YREw9lwOP32stvgbUr5GL8JAVr8+kx2cVywG8/cmVlsxi12Op1z7NHhlZCl3f1UuneoLHYeB8AFn3/VTtocVCHkef2ozp3z6KN/aM1h2cCjduzZE2QdtD8l2lStpHXZKqqK9Qtn4Vce18+7Amsb+LY8p+B466cjKZQ0aW+fXL1dlt5ftLYXN06s4FaPQV3r255dmhEGGsVOK07NdRr5PN9+R7b4RCinOsuGO0kaXQ1b3jJv9Dp3qCx2HgFEAntPH4WBEf738YXx15rKV88PhYEcuG2y95bZ6TCsj2kIgwQ9W5SnaHV6CKefk0GBdmpMO7NrU4EFdoHIHy/nXK0JXvKlXvMM7ZW+7fNQkkcX/f2DNqVHquL7jq25CfHUoJcAD7JvXZ1lGONZO0AP9JD6XMTcECplWPywQr7YlZ1kKXesd1+E6OOsnAKYBuauO6c9Kv1LBJEKthhuKBXLlcH0Wjgzpv3SyHig6R9/HT3GIrymuzNUwcn8HY795ri34YHyvithXxLZB3FPLeSpujHmo59rv3cPf+U7jnpX/G3URUhotAss0EXZWo6dmb2LkBUU6vKHV3OMcYju7e3CKQfCc91HkdePxecl4RN7wz7hht9JLQ9V0xdJKB8wHE7QebBqbZjEkQi99SUUAqPqWDTeets9Garp3W5r3AtT6D8bGisSG5DdeaPCqy01R1qIvS2t81Zv061Gvr4+gUPgqfrFwAXo0d5jlvG0+c/hPUeZ375mpLBVqgfdUTp9cF9bshxtoijlRMPoduJ2upY3vd0oWt0wycAujmg2ESzD6C2MTUdBlDRCJUGqVrV0SL3bLUtH0fm7eIw9e99C7hlh9+/gO2b1zbJoyiHAMDcMuzFIVaWpvCd6WoU6K+ZcB96yWpHcbSnPS8Mr4JW+9aQ56D7p6YjmVrYqSLejNNStTvdVPY9nL4p2DgFADQvQejQIT5jURDqYzHVDAraela3exfNvcA7iGwQlFQgunJ+4r4p5nvjSGO5Uq1TdAw1JudvDK+CS9PXXRq7i1j+2YaK0XfZy+OCYQSgrpSFHGibqhzeHnqovaePHmf/vumfgq6bG5dFE8vh1j28tgEA6kAugUli5YN+0Wn6DDFqA+xeoXPfQnKT7g8zK4JZWIWbRJM75y3R3/oGuCIGPatd62JZSKiEP0Rsn5xXat52rAJmnPfXPXu9WAaKzW5ke+Jiu6ZEs10XH0CvRxi2ctjEwQF0EGoiBt1u8lMoPsMgLFU7gIHGXvuiinZTJ5F2mbv6ixaN5t88MgHsYu0fVepameiSSDuoNwAACAASURBVJnnvJmkpkt6ApInGLpW8xSzZJupTNwz3X5118dndmrqAUCNw2e7q08grs+hE1BjK4xEqeU7JCUogA7i8rCa7IYAtJ/JdnkX4ixDTfZ6sV00njGFf25Zt8qavp9k1r4qH6Uu/AUcwB8/uYI/fnKl6YQHoLVhx1G0ul7Ioponl5R4YSRqqcVPXS/xXFEzbR2us1OfEFwq0dEk5F19AmrXOqB3Qiypoo0/3pwjAyM6zcCFgXYTl9A0ytTy/LEZ7J0saT9zKR+gor7otmJVurHrXlDbeP7ly6tkco9QfklgzCtoxppwRVGuVPEPkyVMnFgsxKa7FgffveRUBMzUC/nabK0tvFaUYHDJR/BRqAXHwoA+ZgwqWY0KkxVCXqxyADSLBcqIJke9GmKpC/9cuWy4zamfZiKqL2EF0EFcIpCoF8vHmemCbdWh1njRjT3OTN1UxsB1Vkn1MR6JhrxDS02zTRsLABYs0UaVak1rfgNar6WuM5aMbgavlmCgCuP5nNuPN+fIxvQyPvefuifquHXRay4+gW5H+5hQx/ZzosRKt/wCQQF0GNvDGkewFpSmMy7YVh3C3AGgRQnIY6ca0xTyEa43SlPrSOIcy0c5UknM1hbqs03PuHmg1aYeRxm4Uq3N49DJSy19EeKavDjQrM1PPVdUfwYGEdLbGslVW+BO5sE4PQB0yOM2CcdetvXbaO2vrA/R7tZ5BBOQhqnpcjNb9O79pzB66L2O1e/2KcEA1Gv+HHziXqxc5q7LC/nIadUB1GeZ4txVM9H2jWu1poe/3/wzox/ANxs7x1jL8t5UkiDJQkmEH2bXzLLOtdlaap3Irs3WjM+mqejeTUX4234joys/MmLpAWDD9Fz0UmavD2pJESpEu1vnEVYACmpDbKC+jJ84bm7HlxYuy+IWGqZRKsJIRbRUlDGtOuRCX1QkiZydvH3j2oZTTi9cqId9arqMGz/pG4TLNt2p6TJ+vGk2lyQhbVNbJzDN2Kl7K/weSWbVYva+aEJsveerG85q1+xdU8Jar2T2+kKZNXOMYYHzrp8H433wwDPGzj700EMPnT17NvNjUWYNoP7S6GrbZwnVPFvG9DLLiJ7B6sM2NV2uFxYjfsdgN01FuXoNGlNUClW+gjrH1Uq0i5NCHFDEPVKvL1W+4/CuulmP+sxHIFHvDPW+2MbUi0I+bojvz/efIk1wXx15LPVxqmzbtg0fffTRR5zzbbrPwwpAwSTkRMy7qWlI2oyPFbX1V2S+q1Tx+p5Rq6JY0NSJkY9BdYFyKbwmesVSDj+qMQxgd/7qGuIEWpErlgKLK0mbg9hU2sMVX5+OKanQVlmzGyQp6dDrvougACRcIibKlWqLoOxEHO+Hn/9gHJMoCQ24xYXrEI5eqo6Lab+Ca7M1FA0P/NR0GQffvdSMihEzfEpQXJutaWPjO0khH4Gx+liydA7bUHMBKHQ5HjoHsam0h+9sN61mQL2UISuTpKRDN4tPuhCcwBJUxISNrON4TS+GriT0G3tG20oIRzlmfeheGd+E1/eMamOqXZ3TlLNu+8a1mDg+0yLArs3WMHFixhh7HifHYeWyXOz4fpliIY/SgUcx/dtH8fWRx9qujS9UWWcboo5T6UC9qZDJCQ64CVJKqO07VsLeyZJXYxaqZPWNn/SNjnq5QYqOJAqrl0tBAxmuABhjOxp/PsI5f1Hz+W4AFQDrOedvUds6SZIZSJazF1PLR/Jh0hXKcYAKJ3RZYcjRReoMkqpqWZvniSJ3tOMYWYZLv3sYU9Nl7J0sxdoHQ3uFVtcwWIqVy4abqwlXZD+IOJ5NjbgIUjI6SHMvnGa7mt9VqjXt6jitWXGWvb1lkppxejlPIRMFwBjbAmAL5/w1xtiLjLH1nPPL0ue7AYBzfoYx9mxDWRTUbZzzM1mMjyJucpP4bVYPJPXCUMJfJ2xd47tNyFEfE8dnWo4hwlHl78nsMwji69UaCjGLnukoV6q497f/l1MHsAfvWYMLV663XFsG4JkH1sVe3lMmq0q15tqorcmPN+fwm3+82NJIyKQvXQWp77Muf1d9zmdvzZElqymTFJDM2dvJUsu9bsZJQiYKgHN+AcAFxlgBwGVZ+De4H8Bk4+/LALYAuF2zrakASqUStm3bpj1eWtFBVEN3MBiTrISJI6sH0veFydrGqo5nVT7Crbl57J2smw904X8mgSPOxzWxyAWq85rMrx9Yh1fGN8VW3KbVjims1ofaAkfNsZWlqUmQiq5uvwkGNM056nNuQ/fcJZ0Vd7LUcr+GoLqQtRN4K4AvNdsLyr9vJ7Z1FFN5YnmbLgoo6wfS54VJO/KAEpDUakDY9sW4gbpyVb8HALkh1vIyyU7irGBAS2emJMKI+q2qzLJ2IBfykXOI8tR0Ge+cL3uNR84H8VXSWdj2O+1I7mUzThIyVQANc85TjLHdnPMT0kcVAGuUr+u2NRkdHU1tpq/iMgO03XzKxNGNyIY0l6y2pbbJti8rPyrUVI5CEPvLUgG4mneSoFshZa7UPGxLPpU8ZeKYR7MylfR6eGW/kEkUEGPsVcbYs41/6gT7X7A4418P4H1iW+aoqdouUQ86eimyIc3IA1tzdJOCK1fqtfkFusYgwjchykyk1cBFR44xvL5nlKznnyYiIuv1PaNtndOywKcIXpJJCaVnCvmoGRklKneK5w6AU0VUH/q1NESvkdUK4E0A64VzV4ryeZ9z/gjn/ARj7AXp8zONz9u2ZU1apptecxTFXbKqqyFKIAshYnMmygXlTE1l0rT/64iT4QoAePtt4De/Aa5cAdatA37/e+CZZ5x/7jvbHomGsHrlcu+Vg89EI0mwg65YnigvQpV8yMI3tpTt8p0kKyfwZdQduYDkyOWcPyL9/Zrm87ZtWWOyJfo4B30eyKyihVz3S31P97JStmshcCjbvsyfPv0Wr4xvMoazZin8TTVpjLz9NvDss8DsbP3f33xT/zfgrAR8Z9t/s2K4xZY/eug9qxLwnWgkdbjLJZptz2+WvrGlapfvJAOfCWxq2+Y7c3F5ILOaEaltEKn9mo5PlYXWzfjk5DMAeOnPn5EF4ET5Bp3gYci+vANV9dLKb36zKPwFs7P17Y4KwHe2/f/++y2M/e49VGZruKOQx99v/llbxyuZOMpNnawURiL8eLM9lJPqu6DW+BHmO7EvzuuhvS4ryEB3GfhMYJ0tMcoxVGbbS/aaMn5tHbUENpt6HKamy8Yer67HN5UOVv0JwKJd9+jpL7Dc0Nhe2ITHx4p48r5iix1Z7vqUFbGv7xV9bSRyuwbf8t5APYpK+KPeOV/Gk/cVm/Z11dkbV3cKH8VXRx7DgcfvxW0r2ueCun0L5S+e97v3n8I+KXP42my9AY4YP3Vvg7O2Nxj4FQA1G/JpZuIzq48TvmYz7ZhKWKj7jdOIWzfj84kF/9Uv72z+ratrlHUTFmCxkJ+XuW3durrZR7fdEZ1pcPvGtWThPZVqbb7Zg/jXD6zD5F++TbVUuUu1WbV0MdAa5mq6b7rP5BVkp7J5A3oGXgEAraabB498YEzV181cfOycvuFrtibxtiJt6n5Nxzc5sl26GukYiYZaom5Mq4ys8Ta3/f73rT4AABgZqW/3QGcafHe6jP/vJ3cbvNx/QSVJlreLk3qB85bSxQ8e+SCRz0a8G+e+udpi3up2g/ReoNMKceBNQCqmmXg0xDB7a67NzGOKblHNQr7ha5RyOXTyUjN81YS6X9Px5fBRYNE5e/DdS83m51RXIx35KIf/tOsXLdtMXb90HaXSxssc9MwzwFtvAXfdVbe93HVX/d8eUUA6pqbLiOOWMF31uDZ1l9+p9ywN+71QaGmbQ/uZtELSfSBXAIyxMdQzecVzpzPncQDnOOfxKm71INQMmTX+R6wO5NmK6Tdiu/j+4V2bcHjXpsRlHVwLilEJbdTxF526iyuBuElMaicvsVrRmXvmOcdsrTPFlr1CIJ95JrHAV4mbiGUirk3d5qTWTU6ShJHK+PaMXup0sryFwGQCWsM5/z9sO2CMPQlgySgAygyyIhpqE7ri5lDRLZRT1qfpRZKXjSpZbItWSkNAFaUeBaoZqxdauvx3/9t/xuFdv+iKqSELAbd949pYvzM9u1Rtobtvb38mxW9WS1FAcbu3ddNBLJtgVjV6QYiIrKzNMd3ok0AqAM75/y3+Zoz9LeqrAcEjnPOXGt97J7PRdQFqhmwq9aD7jSn8zdXO9/LURe3Nz0c5LB8eMs7MdeWMXUn6wKmzxqQKZYgBhjSDWFRrC5n1ebbdX1elrgu/ffK+Iv58/r9iVrEhvXO+jK13rfE+F9+EqpenLuLjL6+2bf+7e9bg7f/4P7Rsc3Ewm0KMO406Xvn96oR/ohvlLVydwK8BOA9A3PnV2QynN5BnyOJlpuSPuDnqrJoqa7Aq75Zf8PLURW2kyMplOfz+f9b3cxUkrXcTZ9VhmjUmUSim8spJSaNEtszUdLmtdHO5UsW+yRLOfXO16Qx3ScQSwl7XevTDz3/ArHJNk5gKfBKq/vTpt9rtn1y+pt0vANL0ZzrHbmCbqKjXWKfoxX7inE83qgm4KoDjyorgfEbj6SlsMxjTzaFuJmPt1RR1Ly/1ot24tWh2Er6EcqWKXGO57VMSmMImoBjqM76v/1p1etALI1Gszl4MaAoIF4U0xIC/XeFXeC2t5fXUdBkTJ2a0ZcM56q02ATSFXWEkwvLhoWbClFxhViRTvf3JFdxRyLdULjWNOU1TAbWKoUw61HbdZKoXhL0Ol+snvqOLzps4MQNwNBPqfFcN3Shv4aoACoyx0wAuoP5ejgHYmdmoUiTJQ2eaEdgEbRxTkozJdio7lF1LAPugztxUOICv/1p1PnbcZCWOusCc2LkBzx+fwbzGDiTMQ/L9+Pn+U85+Btfltal8hku/ZA60rOjUHsNb71rT7E1gWyGastdFRm4S4WEaQ46w6+ccypH2eukGl5XvqnxEru51yt93Zdbpa+SqANYD+F+kf++gvthLJC27QM0IGOAk/HQ3kxIWqiCiXjSB6cFKY6Ylxk4JU5/Z5vUEpZDFcYYAqKpYNHRR8TFh3X173io0qedIjWP3RVe2wxYJMjVdxo2f5tr2FeUYfrw5p41SE/t1fR5MY/jVL+/UmiYfWL86FeXTTVxMczduzXlHxSVZmWW9anINvD7POf9K/IcOlWpOiulBdindkEWJZ9c8ADl7lsKUlZxWLHEa12BV3tzE3MQQY9g7WdIWm9OVlwb8yi98/OXVlmu1d7KE0UPvtVwv6jl6+9P2OPa42MpxiOCBl/58sU0ArR6JsHLZcNs1Ejkc6vOgO0f1WDrKlWrbNc8x1myrqR5jw8v/GXfvP4W7958yHq9XEHkwJkydASmSNGLKOi/AVQHsZ4z9P4yxScbYMfSJArCVH7a9FFnUHHet1f/K+CasXGYWYr5ZyXFI4xpQ1oGVy3JWQW0zhd3z0j/jbkWJq9d4JBryapgimpnbEv3SrmFXbszydNxRyJMmyZFlw+Qqq1Jtr2kltlPChBqDnNcC1J+DPzy9GV//tao9htwHQZSs6AclQIVP24hyrN5CViKJvMiibpiKqwnoVcUJPJbaCDLEt/yweCmAVvNN2kswVzufqak59WCl7SBM4xpQzUpmb83j9T2j2GtoGG9DKAiqRAYD2kImXajW5vH8sRnsmyyRVTGzgDJdXbvxE3kettBjCsqM6JvX4vpspR11lRUTOzeQDn0K4YMC0pMXnXD2mzKBR0WGryz8G/+e1n2v16AicXxCvTrplFHtfVT0TI4xsrlJnFjirO2MpjGNjxWdnKguVGvz+M0/XsTsrfmmsEoit4VySVP4x81pMCkxUx0nXQKjDNWwHXDPa/FRPuVKFXfvP4UcY/jVL+/sSHc2X8bHinjpz59pFYAulFV9F9N6dzqRF2BaAdzDGJOTv9QqvoKr6NFMYGr2ahM4Pho2LeHpWmHT1tnKN5bY5ihPo3+BbUxJG5TI3DCsmrpNlGM4unszxseKuHv/qVT2KddxAtqfdYDOFwFoYeKa1yJHt/notXnO8cdPrjQrnfq8N1lPWKamy2RvC59mOEnpRF6AKRN4SWT4UjP4OC+FSprNXUwhp2qSFQAy4sLXZGOLOqE+f/6Yexata/0h+XOTySMtio34+1OffR8rT8GHHFsU/uLfSRvhFBqlCvZNlpq5IR/vf7gpIPdOlprH0a08dJVe1fIHIj+BSuQS9/HcN1fJaqU2fN4b3Tu3d7KEQycvxev6psFmY5+9NdeWm5EFncgLYLxTxs0EMMbOPvTQQw+dPXs2tX1OTZdx6OSlthffp3csNStS6+e7YItdF/ucmi63tWCMhhiOPrU51oNBHZcB+OrIY8Zx2WrGJGFqupzIN2BDvmZT02X8w7FS6uUmBPIz5Zo3YKOQj/DT3ELb7PDJ+4pkaGqUY1jZcBi7rhAodJ3IVEVSm1/wWpG5vDfUOwck6Pus4JJHIq/meplt27bho48++ohzvk33+UD2AxAPamW21lK8ylfDptVPGLDHrotjHXz3Ulu4X22B4+C7l2I9jDY7o2lctvaTSRgfK2aqAIRDEqgLQFfhz1APa/33n+a0iWkCXYKaS20cF0wZ5X/69FtyZVGb51i5fBilA482t8Wt7a9rs0mttl0T83wycXVUa/PYK62G4j6LLpnrtfn+cGjbGLh+AGps7bXZGn6aW8Dre0a9qnQCtKlI9BP2id+1xa6LY1FJKJVqzbktpe24qn3eJabeJTzNd3xxw/Fc+a5S9SpUN8Tqq6LSgUfxN8vNc6cFXp8lq/b5uMI/Hw21hA1TkVU2s5JrhzgbPuGIriZV+XvUs+KyryTx8lPTZfx4sz3RTsdSKFs9cAogzdhaSnhyTtf7MbF8WH87XB0/cZJGqLwEoD473DdZwvLhIQw5xNGXK1VSuMdJaonTT1fGVp7gjoYzz5UFjuZ4XbKbr8265RMA9ev+xp5RfH3kMbyxZ7TlfryxZxT/5X//H/HVkceakxRTYx0Tug5xcXG9di73UW0TST0raU5IdBw9/YU26VDHUuhrPHAmoDRia1Vb54poqKVmuGu9H3l/rjXZVxPL0yHHInM61KW7rixuNMSQG7JnQsovrNg3EK/ZheoEWxENkdEZKqtHIjz2i5+R9vB8lMP2jWuNJhMdtiZAKqJ7m6mirGr7dgk9piJETD4AdSJBlZVwxVUAinMxmfRk273pWRHX6eC7l6wlGeLM0F1/kxtiXStbnSYDpwCSxtbqhGM+yrVEBbjW+xHoHngh/MUDL5SOTvhHOUYK5jgvgW48tQWOQj7CyuXD2ogQFVW4x1W8ajVJ1wSda7M1vHO+3FJNVETEiETAOFErchMg1RlvGgtlU44b1meKENl61xqySiyw6Ei13UOgnkW9bDjXJmx9x23K95CbBwH2Z0U8EzaHepwZuqti/5vlw5nV4eokA6cAksbWusxkfY9he+BdVgi+SscENZ7r1VrTgSg/6JQQKTec4cJkEWd8cZvRA/X78uHnP7Qo0TQ6k8lNgJI4quNGT6lCRg1JpFYQL09dbFF6tvOv93RujV5KItioDGPRO1vs0/VZkRVBWvHyrjkpOhOgLkRV7QXRawycAkgaW+sykzUl5eji920PvMsKAWgP5xNmDtcqja7Nb8Q5iv2YQvOEySSO4lVfqDhx87ISSqsXr9wEKG5Ip+jYdvT0F9g3WTJWIpWfo+0b17aYeGwRWHHCTlnjHLdvXGsdnw/qCllegYh4/oPvXsLfb2433ZmeFfV9W6XJj3Adt7ovQK8oXetwiV4Qrt3aOr2CGDgFACQr7+A7OxGYksaomZHo8xpH6RRGItyszbeU7m02rUC7sLCFKEY52uZpmjWpttukZYnjIK5zWlEbN36aayqV2FnMDC0mJJ0g1z0zOrMV5UuJE3Yq55ukleQoI94LatJQqbaa7lyfFWo1kGTc1HQjIuz/ZNFAwMkXl9U1NzGQCiAJcU1INseWmknJsdjnNf6SuN1hWpvnOHSyPWfAJmxXLtPbPMVxAdocotpuXUmjPhCweJ3jFEzToRYNBIDnj814rVA4B2rK91VBTs0odeiEj68ClZ/jOE57H2zx/LLpzoek43ZRmlGOafdlq5lkI+trrmPgwkCTQoVN2m6QbRb/4ec/kDM735LMthf/2mytLfzS9oCawh7FspWCA855CTIuXaZc+a5S1V5Hpvy/6zjkMMPxsSIWUsqol++Dz4pFZ5Lw+b38HE9NlxMJMpm48fxxV2tJo/xclCZVomRi5wbyOXLxxXWi+qdKWAEouNjg4piQbLN4U+8CAM3+vy5LYpcHxrXNoDpOFVczQ5zlrOuM2qVcs6g8eu6bq83QT7ki5dR0WRtaaKoeKzK+D528RM7Mc4xhgfNYLSqpe0LV5NHty3RPGYBnlK5q4n66jM+Gr9nT9zi6d5XK4i2MuDUlSiJsx8eKOH7uCj7+8mrLdleHdCeqf6qEFYBElh14bLN4000WL83H+x9uSQaicHlg1EQZU4KNXDRMnc35mBl8k3NcMoHzUQ7P/HIdohw9h5fH/875clOxzHOOd86Xm/b80oFHW5KwVjcat1OsykeYODFjDPH8w9Ob8dWRx5zPRRYU1DPzzAPrmvsTIa1HT3+Bl6cuttyf7RvXkiueYiN6SI1OMd1P38gam0nj8K5NWK0RzC7Hod7Vm8TYXRdoLoqCSoqcmi7jwpXrLdsYgCfvc5swZtGAykZQABJZduCxmY5MAth3DK7ZkqoTWYwPWDR7yJnBuhfO16YuH9NWGkJ3HlGO1atgSmN7ZXwTju7erBUm8gvocn/Hx4r4eP/DeH3PKG7WFshko3yUQ21+gcxJUHs2uJ6LGs6p3hNhHxfCXW6I88dPrrTcH+FMlZ+51xvZxtQkwnQ/fQutucTzT//20bbMZ5fjUPeSShR0ydx2LQNBpX5QPhuqdalKXPNyEoIJSCJrG5zJdGRzpKqx0rbjAHRCmkDnRBa/+65SbYlT1xUNq9bmvcsai2O6RDz4hOwKAa/OxuUX0Of+mmbCooy0rjm6YJ7zthBK2YxXcCxCKLa5RAOp+DpTp6bLZHKYmqzlQtyIORd830mXVbFrGQhqNZeG/EgSoRiHsAKQyKIJvA/jY+Z+pD4mKTGTfWPPqNY8ogtlM5nAqId4nnPnej2uUSa683AxfdleQJ/7S+2LoW6Ks83qRIKTWhpDXV3IPanHfvdemxmHMrO5qlwf4UPlgIicBV+yNGlQ93L1SBT7mC7XyrSvbsuPOGSiABhjBcbYFsbYbsbYq8R3XmCM7WCMvWDa1km6YYNT2b5xrTEixdccdPT0F1ozxW0r2sM6TUKZeojFMlVetv5asVHL30taGsKE7QV0vb9T02UMERFINqe9QBfR9fyxGfx8/yk8f2xGu7q4NltrM+PEMbOp43CNwDLFsceZlWZp0qDMnJyjzezlekyboLbtqxfkhy9ZmYCeBgDO+VuMsfsZY89yzt8SHzYEfIVzfqahKLYA2KFu45xfyGh8WpJmCSdFOCltszsfIUl9V1dOmBI0YrupI5RcLuDtT67gjkK9iiV17VzMA75ZkbYcDdv9pSKBdPtyjdCRkZ3PrsQxs6m4RmBR5+RTllu+Z65mrjiI/ahNnUQiWRxFQz0/rvvqtvyIQyYKQBb2ANYDeF/5yu0ALkv/3kFs66gCADpvg5NxjajxWVL6hJbZBA0HXaXUN4vRJqzjZEW6vIDU/VXr5KjonLq67O1nHljXLD6XFsLMZno2GIC/u2cNvv5rVXtsl4SipHWy1HsmC+ZypYqJ4zM4dPJSS+XcpFnFOr9P3OSpNAR4N+VHHDJ1AjPG1gO4yjk/o3z0JoDnGGNnANwD4EtiW5NSqYRt27Zpj5Nmq8huktQGqcPnpXaZZepqEAHuWYy2UtqmLFiXFzvOCzg1XbY6VRc4b1MkVPa2qSRzHOSCf1RRPA7g0nf/jtKBR8kOXOrzpVth+eSbqNgmMLUF3hTWaZU5SNuUGPf56adZv0zWUUC7OefPqRs555cBvNhQEF8CuKzblvHYeg5qti6SieLOSACz2UPXG9mE7uVyeRFdSmn77C8tTAXwBIWRqFm/Ri4rrbP1f/j5D01BmnQloJrZgHqLRR2iK5yreU23wjq8a1OsEgyA/71Jo8xBN5KnZLpRvydNMlMAjLHdnPPXGn+32PMbNv/1nPMTjLH7Oeev6bbJ+xsdHe2Lmb7NBgrQwjipDZKCmtX41NeX0b1cLi+iz6y+ky+2i+CSa/rbbPlyuWifYmwiM9lWBM2U4StKh9hWfVnUnYlTaympQp/YuaHtGTYVLkybbtTvSZNMFABjbAeAVxljLzU2vdjY/j7n/BHO+QXG2HrG2G4AhwFAt63fcLGBgi121VJnC512IlERQoLVIxF+vDnXEhtNmY9chI6p3IVasjqpPdqHtIrEyfsDaJOIzlHso+gndm4wFt5zeY6yWGGZxkXh0w+CfB/Ui5lOWSYnulG/J02ycgILO766/RHp7xOaz9u29RMuNlAVdbbQSSeS6SFlAKZ/+6izfdNF6JgiZ8R22RTx5H3Flro9tpT6uLZYyqHr04JSIOdXmMIqZXzKBQD1a02Z7eReBXEifrI0naid63z7QVCtRtX3qrbAOzYD77YJKikhEzhF0q5gmDWmma+rIJGxfZcStDo7+qGTl3CzttBWt4dqrJHUFrt8eKj529UjEQ48fi/Z29mIlD7gurLwKRcgOPD4vdYVkkkhZrHCMuWniOxpH4XuYl7p9gy8kyvVLAiZwJ6Y6tfE1fqr8pGxJk7a4xRM7NzgnCWcBrrEIGq1fm225lWXKW4dJ6E45Nj/m41Zf5z7WZvnzWO61mQC/AWWLcnKVtgwiyQt0zmITma6Qny++5O3dzv7thv1e9IkrAA8sM0ybSVuoyHW4gMQ227cmmsKoDSiCFxn9qS04AAAGHtJREFUw7pkmkI+wsEn7s3sAVZXCaZ2kjp8Z3w2wWpSHNT9HGJ0QTD5mOI8XZrFxBFYphWXy+zZ9Ps45jRqxbN6JMKHn//g7Sx1Ma9QM3CfVqhJ6bfYf5mwAvDANstUZwOrR6KWao9Hn9qMo7s3t8wWblsx3OaITVqB1Gc2PD5Wr8j49ZHH8PWRx1A68GhHH2Yqfb6Q15flVQWlWOm49DHWYXJM62Z3b+wZxeXD9WtFZciqvZNtzWLSNhkkbegStyw6dS8PPH5vLAXtUlpBd49EHkYWZd2XGmEF4IFrb14XASpmVz4t/lyJ87J1K5mFch4D7SGUok+y3OjcVHrBRbCaHNOiTwB1HaiKoKKXs+0YQHtGdVLSaOhCTSAOvtveSlTGFAhA5USYxuMaFadbVfZqaGavJY0FBeBBGh5/1w5aSWyYvuPsdjILJWR1mbaT//otJv/ybXPVRAl/V8E6sXMD9k2WtJGENoFBOW7/9Om3Lc5qV+d3GqTR0IWsH9VINFOviYtQi+ssjWNeMU2AuimAu/2e6QgmIA/SqPbnUu8nqUnAd5xZNMJxcULbvqPrk1xb4NbENVGy2eWlGh8rkoJY5CdQ4zOVyKYcrmJ84phpmydMq7wVkdvrbpp8yM/E1HQZo4few97JktXc0klnKTX+wkiUWcc/F7JsOBWXsAKwoM4YnryvaM3UNGGLvU+rSBbgnlCWdiidy0zH5Ttxjy+iqlzvUdEjP0Een8m0I7/Y8n1YrelZK8Je05iZUj1xgXpklcuM05ZoBphXspS5pVPO0u0b17bVdspHOXCOrpqGuh2yqiMoAAM6IRW31KzAVHI3bg0WHT4vW9rJLC4RKC7fiZuh6xtVRdnydfkJ8vhsUV/i2PLzQyGXmkhiGrDV83MtqGdLNLOtZLsl1HQl1UWi3dtEB7dOjbUXk8aCCchAFku2XmwakfaYXGY6Lt/R9tGlOnJL+EZV+SRh6foo54jmMaKHbxziPmcuvW9dBN6Bx+/VPhMivNKmmLsl1Ex9ebudM9CL735QAAayWLL1YuKIGJPcVH35sP+j4ROS6fIy6q7Vnv/+Tqxc5pZcJWO6Zz73U9fP9g9Pb9a+2EmauOjG5eJXcRFmLt+xhVea6KZQM72z3RbAvfjuBxOQgayWbL2aOHJTqntTqbrZiwW26Cb1RXONChHXSnTrMjViz0c5rIiGjKYL3bh19fVdzkEeI9Duc6FCHwv5CCuXDze/e+OnOW0XMlv55r2TJRx891JL4p7NLOUj8FzCK1VEGY04kUJpYHpnfX1jWdBr735QAAb6vc6HD0nL2ppswiIkE0CLc9bVoe4SOisfw1Ye2DWPQEbtCKZCvdi650fNtKbOb/bWXDPskrq+qqJWhdyqfATGkEoXLtNKyRR228nwR5e2oL0kgLtNUAAGemHG0CmSmruo74mQzCQOdZfQ2e8qVRw9/UU9CctQHlgdh4vwz0c5PHlfXQjvmyw5PwemJDc1Sunwrk1t/YjlqB3Tfci6oqxQmKacC1MAQydr5g/SO5sGQQFY6McZg89y2/Zyu5q7VuUjrRljVaOkg48QUNtGUk3aZURct661o1we2LXvskBUsZRbPPrMYNXnR6cI902W8MwD67By+XDbuVZr83j+2IxVUfn6pXTPCOCWkS3jsiJOOrnwNR/14zvbLYICWGL4LLd97fYmiECY5nZXIaBrG+mDrbRGHAf+P818n9oMlopSMfk2XPwTJh+HTagLf4KMeG7kMtkqqxsd7/ZNlprF83TXg7LLc9RXQrYJSq9lzy4lQhRQF3GJ6vDFJ3TVZrf3iVCoEMlHYrtrCJ7vDN0VcRxfB365UiWVUBxlkkXMOaWoqaJuh05ecrrG1dq8UQHfrC2gUq1Zs2qpsuOw/A7ozezZpURQAF0ibsVFGz7LbZvd3iT8VeW1ylK90zUEz1VAvrFnlKzGqYoa+Tg+NfptxIkGSyvmXOQemBqrUMKTyhT2Pb5u388fm2mb0IyPFbFyGW1sMAn0XsyeXUoEBdAlsprZ+CS7xE2M0Skv3UxRjb5wiYF2EZDFRkgfpVSeeWAdeRy1Lk9cxLn5ruImdm5oU1C+MMCpsUoaQnKIwSvHYZ5z7YTGlqBGjdU2sQgkI/gAuoSpXvvLUxdbWuf96pd34pXxTU779QldjRvm6mKm0cWDuzjnfOLY40Z8iM9dqrLK5zOybNhqS7fZp8fHim1VTn2gWmju1djhKdt7IR/hp7kFp3Nf4MDhXZuccxzUcQk/ia2sh06gT02XcePWXNv2pB3req0kczcJCqALTE2XyfjzfDTU4hCc57z5bxcl4CMU4wpQl5nlyLJhr5dKfikLIxGWDw/herVmjWM3KRXTi37wXb0dnAEY1jQv1yU3xa07/8r4Jmy9a01LpNONW3PWKqc5S8KaqoCoomgHn7gXANrCTqljAtCGebooUPGsmBQ7Nek4evoL7TW5bYXfsyUTnMqtBAXQBaiwSwagOreg+aReY951FeATBhcnZM6lSJuP+UF9Ka/N1pCPcnh9z6jz2FRhbwrdBMzRRUd3b3ZSikns07rwUNuseoFzsnKpQDYjUkXRZHOYLSlOlLYW35fHDyxOHqhsajGzl79frlSbysyUQEb2JUjgw+hkTkI/EBRAF6AebN78n3aS1pVxJW5zDxUfG23Sl1I3q9OFVbr4WETJAJfjplkqxMUsNcSYU/aySIqjiqKpxxXHnpoua/sXu5R31oUU69o3+gjZLEqxBKdyK8EJ3AWoB7hYyBsrS2aNa2SS6kg1Rd24kKR/LeAXOvpdpWrtQ+tK2sXFbOchBDNH+zWXuaOQjyXoxsfo/sW2e5FFobMsirelVRE0ixDubhBWAF2AahEoesnqZq+/+uWd3sdxdXaZzA+us7+4TjWTP0QthkYdI041T925jkRDXgIr7bIDPufBUXdM36wtaGfdcXrwis/jzrrTzsDNoqxDGvW9lpIfISiALqCLBOFAszYOgNhRQALXh9Sl0JrL7C/ug2/yh4iX0nYuPo1jtm9ci3+a+V77WW2Ba3vemkhT6Pk2wKnM1vD6nlFSQMYRdFkXQOx2WYc0lMpS8iMEBdAlTn32Pdlx6uP9D3sLfBXXh9TFfJJlzLXJHyK/rKZzcfFJCGTHsEptnnf1JaaE7xADbtxqH3NhJCIFZNIQ2SzCJHtl5pxUqSwlP0JQAF1garpMZmOm9RC5PqS242Vd/pqa9TKgORunZsViuyq0GKvHr6sMsfaesCpZvMSus15K+B589xKA9nHb4gLiCrqsiqktlZlzL7Z2jEtQAF3AFImS1kPk+pCazA6mEL20mNi5AfsmS9oKzkIwmOLfhZJQfRK6ngC2OHsg/Zf45amLLaY+26xXJ3z3EQ3aXdo/pkUayVNLZea8lPqEhCigLpBWFIoJ1wgK3ffkCJNz31zNNNphfKxoreBpCoGlopSO7t7cEpFydPdmFIiyAoK0X+Kp6bI249e35Ee3e9mmVbeKGu8QY8199UN0TRYRT90irAC6gClFP62HyNWWqyboyBE5ajx9VjZbKrlJCAxT8pNLlJLg0MlL5BhsHb/iYOqz4DPrjTPjTLPcQVqmG8pXM885Jo7P4Df/eLHF19HL0TVZmck6TVAAXYB6oUWKPoUu29XUUtH1IRXfe/DIB171XdLCJuBsTl5XYWrKIF3gPPUX2jQun9m7r2M2bWeri+nGpHBaso2JBIbaAkdN4+jWPW+hlk96BAXQBeJEWtiyXU1hnq7HSdr+MS66MgGymUR8rstSBRZNCDYhYPJ3ZGFOMTm4fU1NLso8Tj6HCzZ/kknhAK3hqHES2lVF0wuRREuFoAC6hO8S0iVcU33JfV8W1zj0LISlLnZdN17KhOAiBCZ2bsDE8RnUlBAhtWl8WlAJf888sC51YZVGPgeFbYVGmYgohe2L/LylGUkUVhLBCdw3xJmd+/YccGmWkmW0g228wvmmK4vh4lgdHyvi6FOtzuDVIxGO7t6cyYuvcxa+vmc0cY6HjizzOWxOT+rZTEP4q89bWpFEWTVk6jcyWQEwxgoA1jf+u59z/qLmOy8AuABgC+f8tca23QAqANZzzt/KYmz9SpzZue/LYjO1ACC7T5lwnWm5jHd8rEiGRcapwpk1nTpe1vkcpvPwzWBWKRDlsAv5CAefaC3DnVYM/lLJSUhKViuApwFs5ZyfAADG2LPyhw3hX+Gcn2n8ewtjbAeAy41tlxljWzIaW18SZ3YeJ3xwfIwuCAaA7D5F4TLTEqF/1FF1uQsu30sDOSxx7HfvYfTQe5mHKMYJhTSde9ZhiqaevybyUQ5v7BlF6cCjbWG7Yrs65rQKxNmSCweFTFYAyux9PYD3la/cDuCy9O8dAN4CcJ4x9hTqK4Az8g9KpRK2bdumPd7Zs2cTjrj30TmO5SigwkgEzutJQ6IzVNyEFdOMzneWZJtp2WzXVO5CJxJxdH0KBDr/RBo25bhOTuqadCI+fXys6NRcBgDZB8AnYg1IXqqCSi4cyr7obk+RqROYMbYewFVVmAN4E8BzjLEzAO4B8CXnvMIYexPAcQDq9wOgXxJKaBzetUnbzs/2sqQVdmn6rthusl1TmchZ1quRsdnVTYosbnRKXNNEp64JhS0rmQFeDX5MpGFWo0ycCxxOEWVLhayjgHZzzp9TN3LOLwN4saEgvkTd5LMDwBnO+WuMsVcZY7uFCQkARkdHB2KmHweT0Ph4/8PeD7PNF+BjarHZbCkFwaBvQyiPMeuoDxdFZ1Jkvk1tTB3BetG/IWPzA8jF/XoBU3LhIPkBMosCaghw4dzdony2pfH5ZdSdxCdQdwZfaHzlMIA1WY1tqZFFjZXxsSL+8PTmxPZWm802a3v+1HQZo4few937T2HvZMkr6sNlDDZF5nIPZD9JkrG4HIfyLSQtwWDzURUz8M8kwfQM91ttoiRkogAas/lXGWPnGWPn0RDmjLH3AUAI+kbUz+HGz95ijD3b+O3TIQrInayEaBo1T2z7yKLrk2BquoyJ4zOkbdoWOmoTamkpMpupKY3rYXLG+4REUopC3OfVI+31lnqxUNr4WJGsDdWPVT3jkpUTWNj21e2PSH+fUD6roO4IDniSpVM0DbOCaR9Z2q6Pnv6iLelLxdYiEajXEFLLd6sN1pPcA9MY0qrIasuxcDFf2fwc4r9+SbA6+MS9S6aqZ1xCJvASoNsOwKRkZbt2WcrbZnvjY0UcPf1FmwJQG6wnuQeU/bxYyBv9IDZkQRynKJ36maufQ76fYgz7Jks991z2+3uTBkEBLBG66QDsVWyOyaQzdHV73HvgunrwmVm7lIYAzD2SVeXo6+foh7o9g/7ehFIQASv9UKNdN8aJnRvIuO58NOTsz8jaUe3ia/EtXeBSGkIoGaonRLlSbbnfvtfBtxRJoPOEFUDAiO8srhv2X1MexPLhIVRrC22/WbNyufO4OpF4ppuJytdySJO4ZAozNZl2GKC9N1RPCHG/t29c21KBVrB941rtcZZKB7ClTFAAASM+8e3dWvJTYzx08pJW+APuteyB7tiK1WtJJS5RwtTHryCfvy5D1jZrl30hLmMYpCibXicogIARn1lctwpsUWNUHbcyLrXsKcdmJ3Ax4QC0MKWyuWdvzbVkuroqGpMvhbr+LiunfokYWqoEH0DAiI/dN40lf9qF0Chstey7bad2uWY6O71A+BXUWPdrs7UW34GrojFBXX+bbyOUZO4+QQEEjPgkaiV1luoEwsSJGWsFTpdKqTKrRxZ7L/eqnZq6ZqIXgs5Or1MCK5e3L/JlBZf0PG2+kPGxIj7e/zC+OvJYW1mSXlW+g0RQAAEjPtnAlCC+8dOc06xOJxBq8xyVak07QxSrhX2TJSwfdnuU81EOBx5f7L3cydLSPlCK9w9P18smq4YaSnDaFFzS80xSbbRXle8gEXwAASu+pXrVzNlKtebkDHZ58WVBJ9uXK9Vay6xYJscYFjjX2pg7VVraF5Pj2achjs0Ra6v8aqJYyCey1wcncfcJCiCQKlTmrIsz2LWz1HeVqna1wIE2JWCriZ92hI+LU9PV8UkpXh/BaVNwuvO/8dOctbZ/GkqyV5XvIBEUQCB14i7tXWejdxTy5L446jNTH2GeVoSPS0RRGqGyPoJzfKyIc99cxZ8+/RbznCPHWFtbT/X8dVnE0RDDbSuGUZmtpRatE0oxdJ+gAAKpE3dprwqEwkiEH2/OtRR0E4KOqp2ftH5OElzCYNMIlfURnFPTZbxzvtwM75znHO+cL2PrXWu6UqBPd6wg8LtHUACB1EmytNfNRoWwzzHWFJbbN67FO+fLXTcf+BZcS8vx6So4fRP5wmx8sAhRQIHUSaOPgLwvEREjZrHlShXvnC/jyfuKqRwjLmrYKoW88ul01JGrwgkx+YNJWAEEMiHNpT01iz312feY/u2jxt9mOat1Lbi2feNajP3uPTIzOcuVi6s5rltZ3IHuEhRAoOcxlXqQM1pVIZ91bSKXgmvbN67F5F++RW1ev0ZIq+ELhas5LsTkDyZBAQR6HlN46MF3L+GnuQWtkM96VutScO3BIx8YhX/WDmtXh26IyR9Mgg8g0POYzCOVao0U8lnPal3KZPh03MoKUzkGQZa9mQO9S1AAgZ5nfIxu4E0hZrs6OtnIxXSsXppdp+m4D/QPwQQU6AuoBt4roiGtc1WYOrrRyEVmYucGTJyYaTMDRUOs52bXISZ/8AgKINAXULZsAKSQ74VMU119pEI+wsEn7s1sHCGeP+BKUACBvsE0Q6UEXi/Majs5hn5oxB7oHYICCPQ9vSDke4UQzx/wITiBA4ElRIjnD/gQFEAgsITo1QY3gd4kKICM2LZtG7Zt29btYXSFQT33XjjvbsXz98K5d4N+P+/gAwgElhC9EPkU6B+CAggElhjBKR5wJZiAAoFAYEAJCiAQCAQGlKAAAoFAYEBhnJt6GfUGjLH/umrVquLo6Gi3h+JMqVQCAPTTmNNiUM99UM8bGNxz7/XzLpVKuH79eplz/t/qPu8XBTANYC2Af+v2WAKBQKCP+A8AfuCcj+k+7AsFEAgEAoH0CT6AQCAQGFCCAggEAoEBJSgATxhjWxhjXzLGzjf+e9Xlc8bYNcbY+4yxF1z202ukdd6NbbsZYzsYY892+jzikPK5v9A49xfU4/QaCc57S+Me75a+Oyj3nDr39zt9Di6ETGB/1nDO7wHqNxtAxfHzpzjnZzz202ukct6MsR0ALnPOLzQEwhbO+YUOjD8JaZ37CwAqnPMzDUHR6+ce97xf4pw/xRh7ljG2HsAWAGic97OMsR3Ku9CLpHLunPPLnPMTjLHnOjd0d8IKwBPlwV3POb/s+Hmh8TI47afXSOu8AZwDcLzx0qzvcQEIINVzvx3AVenfO9IdabrEOe/GDP8vDeH3VuM39wMQv72MhkLoZVI8954mKICYMMae5Zyf8Ph8DYCrjLE3ffbTayQ9b855BcCbAI4DuC/TwaZMCvf8TQD3M8YKAO7JcKip4nne96Ch6BhjbzbOtaD85PaMhpo6KZx7TxMUQHwe8fm8MSOoAKjI9kGH/fQaic67YQI601g+q9ei10l07g1zwIuoK4YvsTgr7nW8zhvAl43zPg/gWdTNI2uyGFgHSHruPU1QADGwaXb184Y9sG3Z2w8zBJmUzlu2ex9GnwiGNM5dOAiFWaQfVn6+5w3gL9LfBdSF/1+wuApYD6AnHaIqKZ17TxMUQDzWoNWWC8XLr35+rPGd3QAgvfht++lx0jjvt4QjEMDTnPO3sh1yaiQ+d6H4GtsOZzra9PA678Y9LjTur1gFnQCwvrGt0AcOYEHic2/8ZgeArb242g2ZwIFAIDCghBVAIBAIDChBAQQCgcCAEhRAIBAIDChBAQQCgUCfIwVWaP9NERRAIOBII4/BmsUqSlx0YkyBQINzaE24U/+tJSiAQMCdNS6lKxphjls7MJ7AEkdXSC7NwnpBAQQCgUCPoiYLSnklZxr/TlRPKiiAQMCTxuzr/cb/vyDNyPqi/kugdxBVYaV/P2t5hqjCejuwWGdK928toRx0IOBJo6zxm43/vwzgTc75I43Kn1sB9Euma6DLNMqi72aMAY1np1FLiEJbWI9z/pqy39fgQFgBBALxuKD5u5/KegR6hIaZ57nG37YCgakW1gsKIBAIBLpIw67/ZuPv9Zavp1pYLyiAQMCThs12S+Nl3SH9/Qj6r7x3oIuIbmKc8wuN4nE7ZLu9Wkgu7cJ6oRhcIOBIo/mHU/VSn+8GAt0irAACgUBgQAkKIBBw56pjJvBu1DMxA4GeJpiAAoFAYEAJK4BAIBAYUIICCAQCgQElKIBAIBAYUIICCAQCgQElKIBAIBAYUIICCAQCgQHl/wce/Ln3y1ERHQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot(moon, disk)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "You have adopted the star-planet-moon system with disk. Instead of hyrarchical_bridge.py where the disk was integrated using a direct N-body code, the disk is now a gaseous.\n",
    "\n",
    "\n",
    "Assignmnets and questions:\n",
    "---------------\n",
    "\n",
    "### Assignment 1:\n",
    "\n",
    "\n",
    "### Question 1:\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "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.7.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}