From c8245128fe59d0a8e0c90154b726d89b9b6cef54 Mon Sep 17 00:00:00 2001
From: David Hendriks <davidhendriks93@gmail.com>
Date: Sun, 4 Jul 2021 17:58:55 +0100
Subject: [PATCH] added basic example for ipython notebook
---
examples/notebooks/basic_example.ipynb | 128 +++++++++++++++++++++++++
1 file changed, 128 insertions(+)
create mode 100644 examples/notebooks/basic_example.ipynb
diff --git a/examples/notebooks/basic_example.ipynb b/examples/notebooks/basic_example.ipynb
new file mode 100644
index 000000000..ab0bd54d7
--- /dev/null
+++ b/examples/notebooks/basic_example.ipynb
@@ -0,0 +1,128 @@
+{
+ "cells": [
+ {
+ "cell_type": "code",
+ "execution_count": 1,
+ "id": "66ca8e31-7cbf-4484-bb6c-c9235e14bda2",
+ "metadata": {},
+ "outputs": [],
+ "source": [
+ "import pandas as pd\n",
+ "import numpy as np\n",
+ "\n",
+ "import binarycpython\n",
+ "from binarycpython.utils.custom_logging_functions import binary_c_log_code\n",
+ "from binarycpython.utils.run_system_wrapper import run_system\n",
+ "from binarycpython.utils.functions import example_parse_output\n"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 10,
+ "id": "31851c7f-4ccf-4b5a-9f74-8a3e06c6b2d8",
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "<AxesSubplot:xlabel='time'>"
+ ]
+ },
+ "execution_count": 10,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEGCAYAAABmXi5tAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAbuElEQVR4nO3dfXRc9X3n8fdXj2PLtvwkG2PZSCBjYwzGRhAoCe0GU0zAOOlJTqBNqYNjH5qFpqG7PbBkN+np6YaUnD0pB1rqEwiwB0xch6Y2D6ENhJrNssQ2DQ9+wsYPIAN+BFlI1vN3/5greSwkW9LcmXs19/M6R0d3fvfOzFc/zXz00+/eudfcHRERKXxFURcgIiL5ocAXEUkIBb6ISEIo8EVEEkKBLyKSECVRF3AqkydP9pqamqjLEBEZMTZv3nzY3av6WxfrwK+pqWHTpk1RlyEiMmKY2b6B1mlKR0QkIRT4IiIJocAXEUmIWM/hi4j01dHRQUNDA62trVGXEqlUKkV1dTWlpaWDvk8sA9/MlgBL6urqoi5FRGKmoaGBsWPHUlNTg5lFXU4k3J0jR47Q0NBAbW3toO8Xyykdd1/v7isrKyujLkVEYqa1tZVJkyYlNuwBzIxJkyYN+b+cWAa+iMipJDnsewynDxT4Iqex70gz/2fn4ajLEMmaAl/kNH733pf42kOvRl2GjHAvvfQS119/PQDr1q3jnnvuyXsNsdxpKyIyUrg77k5R0eDHzzfccAM33HBDDqvqn0b4IiJDtHfvXmbPns3NN9/MvHnzWL58OfX19Zx//vl897vf7d3uF7/4BXPmzGHhwoU89dRTve2PPPIIt912GwDLli1j7dq1vevGjBkDwAcffMCVV17JRRddxLx583j55ZezrlsjfBEZsf5q/Ra2vn8s1Mece+Y4vrvk/NNut3PnTh599FEuu+wyjh49ysSJE+nq6uKqq67ijTfe4Nxzz2XFihW8+OKL1NXV8dWvfnVIdTzxxBNcc8013H333XR1ddHS0jLcH6lXLEf4ZrbEzFY1NjZGXYqISL/OOussLrvsMgDWrFnDwoULWbBgAVu2bGHr1q1s376d2tpaZs2ahZnxta99bUiPf8kll/CTn/yE733ve7z55puMHTs265pjOcJ39/XA+vr6+hVR1yIi8TWYkXiuVFRUALBnzx5++MMfsnHjRiZMmMCyZcuGdHx8SUkJ3d3dAHR3d9Pe3g7AlVdeyYYNG3jmmWdYtmwZd9xxBzfffHNWNcdyhC8iMlIcO3aMiooKKisrOXDgAM899xwAc+bMYe/evbzzzjsArF69ut/719TUsHnzZiB99E5HRwcA+/btY+rUqaxYsYJvfOMbvPbaa1nXGssRvojISDF//nwWLFjAnDlzmDFjBldccQWQPtfNqlWruO666xg9ejSf+9znaGpq+tT9V6xYwdKlS5k/fz6LFy/u/c/hpZde4t5776W0tJQxY8bw2GOPZV2ruXvWD5Ir9fX1rgugSNRq7nwGgL33XBdxJQKwbds2zjvvvKjLiIX++sLMNrt7fX/ba0pHRCQhFPgiIgmhwBeRESfOU9H5Mpw+UOCLyIiSSqU4cuRIokO/53z4qVRqSPfTUToiMqJUV1fT0NDAoUOHoi4lUj1XvBoKBb7IAA42tVI5avCXj5P8KC0tHdJVnuSEvAW+mX0RuA4YBzzk7v+ar+cWGY5L/+YFrjl/atRliIQmqzl8M3vYzA6a2Vt92heb2Q4z22VmdwK4+8/dfQVwKzC0swiJROT5LQeiLkEkNNnutH0EWJzZYGbFwAPAtcBc4CYzm5uxyXeC9SKxleQdglK4sgp8d98AHO3TfCmwy913u3s78CSw1NJ+ADzn7gOeFMLMVprZJjPblPSdMhId5b0UolwcljkdeC/jdkPQdjuwCPiymd060J3dfZW717t7fVVVVQ7KEzm97iDxda1sKSR522nr7vcB9w1mWzNbAiypq6vLbVEiA+gORvgGaLAvhSIXI/z9wIyM29VB26C5+3p3X1lZWRlqYSKD5UHMF2mILwUkF4G/EZhlZrVmVgbcCKwbygPoilcStZ45fAW+FJJsD8tcDbwCzDazBjNb7u6dwG3A88A2YI27bxnK42qEL1HTHL4Uoqzm8N39pgHanwWezeaxRaLUrRG+FKBYnjxNUzoStZ4RfpHyXgpILANfUzoSNU9fUxrTCF8KSCwDXyRqmsOXQhTLwNeUjkSt59h7zeFLIYll4GtKR6KmOXwpRLEMfJGonQh8Jb4UDgW+SD96PnilnbZSSGIZ+JrDl6hpSkcKUSwDX3P4ErXek6cp8KWAxDLwRaLmmsOXAqTAF+mHTp4mhSiWga85fIlaZzCnUxTLd4jI8MTy5aw5fIlaZ1f63AqlSnwpIHo1i/Sjoys9wi8p1pSOFA4Fvkg/OrvTI/wSjfClgOjVLNKPnhF+qUb4UkAU+CL96JnDLynWW0QKRyxfzTpKR6LWc5ROiT5qKwUkloGvo3Qkah09R+lohC8FRK9mkX506igdKUAKfJF+6CgdKUR6NYv0Q0fpSCFS4Iv0o3eErzl8KSB6NYv0o70zHfjlJXqLSOHQq1mkH8fbuwAYVVoccSUi4Yll4Os4fIlaS0cQ+GUKfCkcsQx8HYcvUWsNRvia0pFColezSD+OByN8HaMjhUSBL9KPnsBvCz5xK1IIFPgi/Tjeng76o5+0R1yJSHgU+CL9aA1G+B8ea424EpHwKPBF+tHc3gnAu0dbIq5EJDwKfJF+NB7vAGDfkXTgV+jwTCkACnyRfvQEfg8djy+FQIEv0o9jfQJfpBAo8EX6cPeTRvjVE0ZFWI1IePIW+GZ2tpk9ZGZr8/WcIsNxvKOLji5n8pgy7v3yhfzuuVVRlyQSiqwC38weNrODZvZWn/bFZrbDzHaZ2Z0A7r7b3Zdn83wi+fBxS3p0/xe/P5uv1M+IuBqR8GQ7wn8EWJzZYGbFwAPAtcBc4CYzm5vl84jkzcGmNgCmjC2PuBKRcGUV+O6+ATjap/lSYFcwom8HngSWDvYxzWylmW0ys02HDh3KpjyRYTkQfNhq6rhUxJWIhCsXc/jTgfcybjcA081skpk9CCwws7sGurO7r3L3enevr6rS3Knk38Eg8KeM0whfCktJvp7I3Y8Atw5mWzNbAiypq6vLbVEi/ThwrI3iImNShQJfCksuRvj7gcw9XdVB26DpfPgSpQ+PtVI1ppziIp0cWQpLLgJ/IzDLzGrNrAy4EVg3lAfQFa8kSu8eaWHmxNFRlyESumwPy1wNvALMNrMGM1vu7p3AbcDzwDZgjbtvGcrjaoQvUdpzpJmayQp8KTxZzeG7+00DtD8LPJvNY4tE4ZO2Tg41tVEzuSLqUkRCF8tTK2hKR6Ky93AzALWTFPhSeGIZ+JrSkai8faAJgLopYyKuRCR8sQx8kai8ub+RUaXFnF2lwJfCE8vA15SOROWt/Y3MPXOcDsmUghTLwNeUjkShq9vZ+v4xLpiu150UplgGvkgUtrzfSHN7Fwtmjo+6FJGcUOCLBH696wgAv3PO5IgrEcmNWAa+5vAlCv/3ncOcO3UMVTotshSoWAa+5vAl3z5p6+TVPUf5bJ3O0CqFK5aBL5JvL2w7QHtnN9decEbUpYjkjAJfBHjmjQ+YOq6ci2dOiLoUkZyJZeBrDl/y6fAnbfxqx0Guu+BMinT8vRSwWAa+5vAln9Zseo+OLucPPzMz6lJEciqWgS+SLx1d3Tz+/97l8rMn6fw5UvAU+JJo//zafvZ/fJzln62NuhSRnFPgS2J1dHVz34s7ubC6kqvOmxJ1OSI5p8CXxFq1YTcNHx3n21efi5l21krhi2Xg6ygdybWt7x/jR798m+svnMZ/mq3RvSRDLANfR+lILrV1dnHHmt9SOaqMv146L+pyRPImq2vaioxE972wk+0fNvHjm+uZUFEWdTkieRPLEb5Irmze9xH/8NI7fOXiahbNnRp1OSJ5pcCXxDja3M7tT7zGtMpR/Pclc6MuRyTvNKUjidDV7Xzryf/g8CftrP3TyxmXKo26JJG8U+BLwXN3/vrprby88zD/80sXcGH1+KhLEomEAl8KWmtHF9/5+Vus3dzALVfUctOlM6IuSSQysQx8M1sCLKmrq4u6FBnB3mxo5C9/9gbbPjjGn101i28vmqUPWEmixTLw3X09sL6+vn5F1LXIyLP3cDP3/2oXP3utgUkV5Ty8rJ7Pz9EROSKxDHyRoero6uaFbQd4/NV3eXnnYcpLilh+RS1/tmiWdtCKBBT4MmI1tXbw728f4t+2HuDF7Qdpau1kWmWKP180iz+8dCZTxqWiLlEkVhT4MmK0dnTx2rsf8Zs9R3l191E27TtKR5czsaKMxeefwbUXnMGVs6ooKdbHS0T6o8CXWGtq7eCRX+9lw85DvP5eI+1d3ZjB3Gnj+PoVtSw6byoXnzWBYl2aUOS0FPgSK40tHew40MSOA01sff8Yv9x2gENNbcyfMZ6vX1HDZ86eyMVnTaRylOblRYZKgS+R+Ki5nT1Hmtl9qJm3DzSx/cMm3v6wiQ+PtfZuMy5VwsKzJvDtRecyf8b46IoVKRAKfMkJd+doczvvfXScvYeb2XO4mb1Hmtl7pIW9h5tpPN7Ru21ZSRF1VWO4/JxJzD5jLLOnjmX2GWOZVpnScfMiIVLgy7C0d3bzYWMr+z8+zv6Pj/N+8JV5u7Wju3d7MzizchQ1k0dz/YXTqJ1cwVmTKqidXEHNpNHa0SqSBwp8OUlLeycHj7VxsKmNg02tJy0famoLbrfyUUvHp+47eUw508enmHPGWD4/ewrTJ4xi+vhR1E6uYMbE0aRKiyP4iUSkR94C38wqgL8H2oGX3P3xfD13UnV1O43HOzja3N779VFLxnJzO0f73G5u7/rU45QWG1Vjyqkal2LmpNHU10xgytgU08anmD5+FGeOH8W0ypQCXSTmsgp8M3sYuB446O7zMtoXA38HFAM/dvd7gD8A1rr7ejP7KaDAP42Orm6a2zppau3kWGsHx473fO/gWGtn8D3d3ti73JHe/ngHTW2dAz52RVkxEyrKmBh81VWNYUJFGZPGlDFlbIopY8uZMq6cKWNTjB9VSpEOexQZ8bId4T8C3A881tNgZsXAA8DVQAOw0czWAdXAm8Fmnx5GjlDd3U5rZxetHd20dnQFX91BWxctbV00t3fS3NZFS/A9fbuTlvYuPmnr7G3PXN/S1kV7V/dpn39sqoRxqVLGjSplXKqEGRNHB7dLGJsqZeLo0pOCfWJFGRNGl2k0LpJAWQW+u28ws5o+zZcCu9x9N4CZPQksJR3+1cBvOcWVtsxsJbASYObMmcOqq/F4B20dXbR1dtN2UhgH3zu7ON7eRWtnN22ZIR2sy9y+rXPg9W0d3YMK5b5SpUVUlJVQUV7C6LJiKspLGJsqYVplitFlJVSUp9sqyooZXVbCmFQJlaNKe4O8J+DHlJfoA0ciMmi5mMOfDryXcbsB+AxwH3C/mV0HrB/ozu6+ClgFUF9f78Mp4EsP/Jrdh5uHdB8zSJUUkyotIlVaTKq0mPKSnuUiJo0p+9T6nnWp0mJSJUUntZWXFpMqKc4I7xJGlxdTUaaQFpFo5G2nrbs3A18fzLbZng//9qvqaGnvorykmLKSok+FcSoI495gLi2irLhIx3yLSEHLReDvBzIvK1QdtA1atufD/9KC6uHcTUSkoOXi0y4bgVlmVmtmZcCNwLocPI+IiAxBVoFvZquBV4DZZtZgZsvdvRO4DXge2AascfctQ3zcJWa2qrGxMZvyREQkQ7ZH6dw0QPuzwLNZPK4ucSgiErJYnsBEI3wRkfDFMvDdfb27r6ysrIy6FBGRghHLwBcRkfDFMvA1pSMiEr5YBr6mdEREwhfLwBcRkfAp8EVEEiKWga85fBGR8MUy8DWHLyISvlgGvoiIhE+BLyKSELEMfM3hi4iEL5aBrzl8EZHwxTLwRUQkfAp8EZGEUOCLiCSEAl9EJCFiGfg6SkdEJHyxDHwdpSMiEr5YBr6IiIRPgS8ikhAKfBGRhFDgi4gkhAJfRCQhFPgiIgkRy8DXcfgiIuGLZeDrOHwRkfDFMvBFRCR8CnwRkYRQ4IuIJIQCX0QkIRT4IiIJocAXEUkIBb6ISEIo8EVEEkKBLyKSEHkLfDM728weMrO1+XpOERE5YVCBb2YPm9lBM3urT/tiM9thZrvM7M5TPYa773b35dkUKyIiw1cyyO0eAe4HHutpMLNi4AHgaqAB2Ghm64Bi4Pt97n+Lux/MuloRERm2QQW+u28ws5o+zZcCu9x9N4CZPQksdffvA9cPtyAzWwmsBJg5c+ZwH0ZERPrIZg5/OvBexu2GoK1fZjbJzB4EFpjZXQNt5+6r3L3e3eurqqqyKE9ERDINdkona+5+BLh1MNua2RJgSV1dXW6LEhFJkGxG+PuBGRm3q4O2rOl8+CIi4csm8DcCs8ys1szKgBuBdeGUJSIiYRvsYZmrgVeA2WbWYGbL3b0TuA14HtgGrHH3LWEUpUscioiEb7BH6dw0QPuzwLOhVpR+3PXA+vr6+hVhP7aISFLF8tQKGuGLiIQvloGvnbYiIuGLZeCLiEj4Yhn4mtIREQlfLANfUzoiIuGLZeCLiEj4FPgiIgkRy8DXHL6ISPhiGfiawxcRCV8sA19ERMKnwBcRSQgFvohIQsQy8LXTVkQkfLEMfO20FREJXywDX0REwqfAFxFJCAW+iEhCxDLwtdNWRCR8sQx87bQVEQlfLANfRETCp8AXEUkIBb6ISEIo8EVEEkKBLyKSEAp8EZGEiGXg6zh8EZHwxTLwdRy+iEj4Yhn4IiISPgW+iEhCKPBFRBJCgS8ikhAKfBGRhFDgi4gkhAJfRCQhSqIuQGQkONrcztX/69+jLkMS5NlvfY7S4nDH5Ap8kdP44oLpfNzSgeNRlyKSlbwFvpl9EbgOGAc85O7/mq/nFsnGJTUTuaRmYtRliGRtUP8vmNnDZnbQzN7q077YzHaY2S4zu/NUj+HuP3f3FcCtwFeHX7KIiAzHYEf4jwD3A4/1NJhZMfAAcDXQAGw0s3VAMfD9Pve/xd0PBsvfCe4nIiJ5NKjAd/cNZlbTp/lSYJe77wYwsyeBpe7+feD6vo9hZgbcAzzn7q8N9FxmthJYCTBz5szBlCciIoOQzS7g6cB7GbcbgraB3A4sAr5sZrcOtJG7r3L3enevr6qqyqI8ERHJlLedtu5+H3DfYLY1syXAkrq6utwWJSKSINmM8PcDMzJuVwdtWdP58EVEwpdN4G8EZplZrZmVATcC68IpS0REwjbYwzJXA68As82swcyWu3sncBvwPLANWOPuW8IoSpc4FBEJn7nH99ODZnYI2DfMu08GDodYTtjiXh+oxjDEvT6If41xrw/iVeNZ7t7vES+xDvxsmNkmd6+Puo6BxL0+UI1hiHt9EP8a414fjIwaQWfLFBFJDAW+iEhCFHLgr4q6gNOIe32gGsMQ9/og/jXGvT4YGTUW7hy+iIicrJBH+CIikkGBLyKSEAUX+EM5R38OnnuGmf3KzLaa2RYz+1bQPtHM/s3MdgbfJwTtZmb3BbW+YWYLMx7rT4Ltd5rZn4RcZ7GZ/YeZPR3crjWzV4M6fhp8chozKw9u7wrW12Q8xl1B+w4zuybk+sab2Voz225m28zs8jj1oZl9O/j9vmVmq80sFXUf9nfNijD7zMwuNrM3g/vcZ2YWUo33Br/nN8zsn81sfMa6fvtnoPf4QL+DbOrLWPcXZuZmNjm4HUkfZs3dC+aL9Ln43wHOBsqA14G5eXz+acDCYHks8DYwF/hb4M6g/U7gB8HyF4DnAAMuA14N2icCu4PvE4LlCSHWeQfwBPB0cHsNcGOw/CDwp8HyN4EHg+UbgZ8Gy3ODvi0HaoM+Lw6xvkeBbwTLZcD4uPQh6TPC7gFGZfTdsqj7ELgSWAi8ldEWWp8Bvwm2teC+14ZU4+8DJcHyDzJq7Ld/OMV7fKDfQTb1Be0zSJ9RYB8wOco+zPr1m+8nzOkPA5cDz2fcvgu4K8J6/oX0BWJ2ANOCtmnAjmD5H4GbMrbfEay/CfjHjPaTtsuypmrgBeDzwNPBi+9wxpuutw+DF/nlwXJJsJ317dfM7UKor5J0oFqf9lj0ISdOCz4x6JOngWvi0IdADSeHaSh9FqzbntF+0nbZ1Nhn3ZeAx4PlfvuHAd7jp3odZ1sfsBaYD+zlROBH1ofZfBXalM5Qz9GfM8G/7guAV4Gp7v5BsOpDYGqwPFC9ufw5fgT8JdAd3J4EfOzpcyP1fa7eOoL1jcH2uayvFjgE/MTS004/NrMKYtKH7r4f+CHwLvAB6T7ZTLz6sEdYfTY9WM5lrQC3kB75DqfGU72Oh83MlgL73f31Pqvi2oenVGiBHwtmNgb4GfDn7n4sc52n/7xHciysmV0PHHT3zVE8/yCVkP63+h/cfQHQTHo6olfEfTgBWEr6D9OZQAWwOIpahiLKPhsMM7sb6AQej7qWHmY2GvhvwP+IupawFFrg5+wc/YNlZqWkw/5xd38qaD5gZtOC9dOAnuv7DlRvrn6OK4AbzGwv8CTpaZ2/A8abWc/FcDKfq7eOYH0lcCSH9UF65NPg7q8Gt9eS/gMQlz5cBOxx90Pu3gE8Rbpf49SHPcLqs/3Bck5qNbNlpC+L+kfBH6bh1HiEgX8Hw3UO6T/srwfvmWrgNTM7Yxj15bQPBy3fc0i5/CI9OtxN+pfUs0Pn/Dw+v5G+0PuP+rTfy8k7z/42WL6Ok3f8/CZon0h6HntC8LUHmBhyrb/HiZ22/8TJO7u+GSz/Z07e4bgmWD6fk3eo7SbcnbYvA7OD5e8F/ReLPgQ+A2wBRgfP+Sjpy3dG3od8eg4/tD7j0zscvxBSjYuBrUBVn+367R9O8R4f6HeQTX191u3lxBx+ZH2Y1Wsk30+Y8x8ovff8bdJ78u/O83N/lvS/zW8Avw2+vkB6fvEFYCfwy4wXgAEPBLW+CdRnPNYtwK7g6+s5qPX3OBH4Zwcvxl3Bm6Y8aE8Ft3cF68/OuP/dQd07CPloA+AiYFPQjz8P3jix6UPgr4DtwFvA/w5CKdI+BFaT3qfQQfq/pOVh9hlQH/y87wD302enehY17iI9593zfnnwdP3DAO/xgX4H2dTXZ/1eTgR+JH2Y7ZdOrSAikhCFNocvIiIDUOCLiCSEAl9EJCEU+CIiCaHAFxFJCAW+SMDSZ+n8ZrB8ppmtjbomkTDpsEyRQHD+o6fdfV7UtYjkQsnpNxFJjHuAc8zst6Q/rHSeu88LPvr/RdLnzZlF+uRpZcAfA22kPzF51MzOIf1hnCqgBVjh7tvz/UOIDERTOiIn3Am84+4XAf+1z7p5wB8AlwB/A7R4+uRurwA3B9usAm5394uB/wL8fT6KFhksjfBFBudX7t4ENJlZI7A+aH8TuDA4Q+rvAP+UcSGj8vyXKTIwBb7I4LRlLHdn3O4m/T4qIn0+9ovyXJfIoGlKR+SEJtKXphwyT1/3YI+ZfQV6r3k6P8ziRLKlwBcJuPsR4NfBRazvHcZD/BGw3MxeJ30K5aVh1ieSLR2WKSKSEBrhi4gkhAJfRCQhFPgiIgmhwBcRSQgFvohIQijwRUQSQoEvIpIQ/x8q6Ro0odz3FAAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ "<Figure size 432x288 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": [
+ "# You can decide to `write` your own logging_line, which allows you to write a more complex logging statement with conditionals.\n",
+ "logging_line = 'Printf(\"MY_STELLAR_DATA time=%g radius=%g\\\\n\",stardata->model.time,stardata->star[0].radius)'\n",
+ "\n",
+ "# Generate entire shared lib code around logging lines\n",
+ "custom_logging_code = binary_c_log_code(logging_line, verbose=-1)\n",
+ "\n",
+ "# Run system. all arguments can be given as optional arguments. the custom_logging_code is \n",
+ "# one of them and will be processed automatically.\n",
+ "output = run_system(\n",
+ " M_1=1, metallicity=0.002, M_2=0.1, separation=0, orbital_period=100000000000, \\\n",
+ " custom_logging_code=custom_logging_code,\n",
+ ")\n",
+ "\n",
+ "df = pd.DataFrame.from_dict(example_parse_output(output, 'MY_STELLAR_DATA'), dtype=np.float64)\n",
+ "\n",
+ "df.plot('time', 'radius', logy=True)"
+ ]
+ },
+ {
+ "cell_type": "code",
+ "execution_count": 7,
+ "id": "db0f3990-8680-4ea7-80d3-cdb430950c7f",
+ "metadata": {},
+ "outputs": [
+ {
+ "data": {
+ "text/plain": [
+ "<AxesSubplot:xlabel='time'>"
+ ]
+ },
+ "execution_count": 7,
+ "metadata": {},
+ "output_type": "execute_result"
+ },
+ {
+ "data": {
+ "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXwAAAEGCAYAAABmXi5tAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAbuElEQVR4nO3dfXRc9X3n8fdXj2PLtvwkG2PZSCBjYwzGRhAoCe0GU0zAOOlJTqBNqYNjH5qFpqG7PbBkN+np6YaUnD0pB1rqEwiwB0xch6Y2D6ENhJrNssQ2DQ9+wsYPIAN+BFlI1vN3/5greSwkW9LcmXs19/M6R0d3fvfOzFc/zXz00+/eudfcHRERKXxFURcgIiL5ocAXEUkIBb6ISEIo8EVEEkKBLyKSECVRF3AqkydP9pqamqjLEBEZMTZv3nzY3av6WxfrwK+pqWHTpk1RlyEiMmKY2b6B1mlKR0QkIRT4IiIJocAXEUmIWM/hi4j01dHRQUNDA62trVGXEqlUKkV1dTWlpaWDvk8sA9/MlgBL6urqoi5FRGKmoaGBsWPHUlNTg5lFXU4k3J0jR47Q0NBAbW3toO8Xyykdd1/v7isrKyujLkVEYqa1tZVJkyYlNuwBzIxJkyYN+b+cWAa+iMipJDnsewynDxT4Iqex70gz/2fn4ajLEMmaAl/kNH733pf42kOvRl2GjHAvvfQS119/PQDr1q3jnnvuyXsNsdxpKyIyUrg77k5R0eDHzzfccAM33HBDDqvqn0b4IiJDtHfvXmbPns3NN9/MvHnzWL58OfX19Zx//vl897vf7d3uF7/4BXPmzGHhwoU89dRTve2PPPIIt912GwDLli1j7dq1vevGjBkDwAcffMCVV17JRRddxLx583j55ZezrlsjfBEZsf5q/Ra2vn8s1Mece+Y4vrvk/NNut3PnTh599FEuu+wyjh49ysSJE+nq6uKqq67ijTfe4Nxzz2XFihW8+OKL1NXV8dWvfnVIdTzxxBNcc8013H333XR1ddHS0jLcH6lXLEf4ZrbEzFY1NjZGXYqISL/OOussLrvsMgDWrFnDwoULWbBgAVu2bGHr1q1s376d2tpaZs2ahZnxta99bUiPf8kll/CTn/yE733ve7z55puMHTs265pjOcJ39/XA+vr6+hVR1yIi8TWYkXiuVFRUALBnzx5++MMfsnHjRiZMmMCyZcuGdHx8SUkJ3d3dAHR3d9Pe3g7AlVdeyYYNG3jmmWdYtmwZd9xxBzfffHNWNcdyhC8iMlIcO3aMiooKKisrOXDgAM899xwAc+bMYe/evbzzzjsArF69ut/719TUsHnzZiB99E5HRwcA+/btY+rUqaxYsYJvfOMbvPbaa1nXGssRvojISDF//nwWLFjAnDlzmDFjBldccQWQPtfNqlWruO666xg9ejSf+9znaGpq+tT9V6xYwdKlS5k/fz6LFy/u/c/hpZde4t5776W0tJQxY8bw2GOPZV2ruXvWD5Ir9fX1rgugSNRq7nwGgL33XBdxJQKwbds2zjvvvKjLiIX++sLMNrt7fX/ba0pHRCQhFPgiIgmhwBeRESfOU9H5Mpw+UOCLyIiSSqU4cuRIokO/53z4qVRqSPfTUToiMqJUV1fT0NDAoUOHoi4lUj1XvBoKBb7IAA42tVI5avCXj5P8KC0tHdJVnuSEvAW+mX0RuA4YBzzk7v+ar+cWGY5L/+YFrjl/atRliIQmqzl8M3vYzA6a2Vt92heb2Q4z22VmdwK4+8/dfQVwKzC0swiJROT5LQeiLkEkNNnutH0EWJzZYGbFwAPAtcBc4CYzm5uxyXeC9SKxleQdglK4sgp8d98AHO3TfCmwy913u3s78CSw1NJ+ADzn7gOeFMLMVprZJjPblPSdMhId5b0UolwcljkdeC/jdkPQdjuwCPiymd060J3dfZW717t7fVVVVQ7KEzm97iDxda1sKSR522nr7vcB9w1mWzNbAiypq6vLbVEiA+gORvgGaLAvhSIXI/z9wIyM29VB26C5+3p3X1lZWRlqYSKD5UHMF2mILwUkF4G/EZhlZrVmVgbcCKwbygPoilcStZ45fAW+FJJsD8tcDbwCzDazBjNb7u6dwG3A88A2YI27bxnK42qEL1HTHL4Uoqzm8N39pgHanwWezeaxRaLUrRG+FKBYnjxNUzoStZ4RfpHyXgpILANfUzoSNU9fUxrTCF8KSCwDXyRqmsOXQhTLwNeUjkSt59h7zeFLIYll4GtKR6KmOXwpRLEMfJGonQh8Jb4UDgW+SD96PnilnbZSSGIZ+JrDl6hpSkcKUSwDX3P4ErXek6cp8KWAxDLwRaLmmsOXAqTAF+mHTp4mhSiWga85fIlaZzCnUxTLd4jI8MTy5aw5fIlaZ1f63AqlSnwpIHo1i/Sjoys9wi8p1pSOFA4Fvkg/OrvTI/wSjfClgOjVLNKPnhF+qUb4UkAU+CL96JnDLynWW0QKRyxfzTpKR6LWc5ROiT5qKwUkloGvo3Qkah09R+lohC8FRK9mkX506igdKUAKfJF+6CgdKUR6NYv0Q0fpSCFS4Iv0o3eErzl8KSB6NYv0o70zHfjlJXqLSOHQq1mkH8fbuwAYVVoccSUi4Yll4Os4fIlaS0cQ+GUKfCkcsQx8HYcvUWsNRvia0pFColezSD+OByN8HaMjhUSBL9KPnsBvCz5xK1IIFPgi/Tjeng76o5+0R1yJSHgU+CL9aA1G+B8ea424EpHwKPBF+tHc3gnAu0dbIq5EJDwKfJF+NB7vAGDfkXTgV+jwTCkACnyRfvQEfg8djy+FQIEv0o9jfQJfpBAo8EX6cPeTRvjVE0ZFWI1IePIW+GZ2tpk9ZGZr8/WcIsNxvKOLji5n8pgy7v3yhfzuuVVRlyQSiqwC38weNrODZvZWn/bFZrbDzHaZ2Z0A7r7b3Zdn83wi+fBxS3p0/xe/P5uv1M+IuBqR8GQ7wn8EWJzZYGbFwAPAtcBc4CYzm5vl84jkzcGmNgCmjC2PuBKRcGUV+O6+ATjap/lSYFcwom8HngSWDvYxzWylmW0ys02HDh3KpjyRYTkQfNhq6rhUxJWIhCsXc/jTgfcybjcA081skpk9CCwws7sGurO7r3L3enevr6rS3Knk38Eg8KeM0whfCktJvp7I3Y8Atw5mWzNbAiypq6vLbVEi/ThwrI3iImNShQJfCksuRvj7gcw9XdVB26DpfPgSpQ+PtVI1ppziIp0cWQpLLgJ/IzDLzGrNrAy4EVg3lAfQFa8kSu8eaWHmxNFRlyESumwPy1wNvALMNrMGM1vu7p3AbcDzwDZgjbtvGcrjaoQvUdpzpJmayQp8KTxZzeG7+00DtD8LPJvNY4tE4ZO2Tg41tVEzuSLqUkRCF8tTK2hKR6Ky93AzALWTFPhSeGIZ+JrSkai8faAJgLopYyKuRCR8sQx8kai8ub+RUaXFnF2lwJfCE8vA15SOROWt/Y3MPXOcDsmUghTLwNeUjkShq9vZ+v4xLpiu150UplgGvkgUtrzfSHN7Fwtmjo+6FJGcUOCLBH696wgAv3PO5IgrEcmNWAa+5vAlCv/3ncOcO3UMVTotshSoWAa+5vAl3z5p6+TVPUf5bJ3O0CqFK5aBL5JvL2w7QHtnN9decEbUpYjkjAJfBHjmjQ+YOq6ci2dOiLoUkZyJZeBrDl/y6fAnbfxqx0Guu+BMinT8vRSwWAa+5vAln9Zseo+OLucPPzMz6lJEciqWgS+SLx1d3Tz+/97l8rMn6fw5UvAU+JJo//zafvZ/fJzln62NuhSRnFPgS2J1dHVz34s7ubC6kqvOmxJ1OSI5p8CXxFq1YTcNHx3n21efi5l21krhi2Xg6ygdybWt7x/jR798m+svnMZ/mq3RvSRDLANfR+lILrV1dnHHmt9SOaqMv146L+pyRPImq2vaioxE972wk+0fNvHjm+uZUFEWdTkieRPLEb5Irmze9xH/8NI7fOXiahbNnRp1OSJ5pcCXxDja3M7tT7zGtMpR/Pclc6MuRyTvNKUjidDV7Xzryf/g8CftrP3TyxmXKo26JJG8U+BLwXN3/vrprby88zD/80sXcGH1+KhLEomEAl8KWmtHF9/5+Vus3dzALVfUctOlM6IuSSQysQx8M1sCLKmrq4u6FBnB3mxo5C9/9gbbPjjGn101i28vmqUPWEmixTLw3X09sL6+vn5F1LXIyLP3cDP3/2oXP3utgUkV5Ty8rJ7Pz9EROSKxDHyRoero6uaFbQd4/NV3eXnnYcpLilh+RS1/tmiWdtCKBBT4MmI1tXbw728f4t+2HuDF7Qdpau1kWmWKP180iz+8dCZTxqWiLlEkVhT4MmK0dnTx2rsf8Zs9R3l191E27TtKR5czsaKMxeefwbUXnMGVs6ooKdbHS0T6o8CXWGtq7eCRX+9lw85DvP5eI+1d3ZjB3Gnj+PoVtSw6byoXnzWBYl2aUOS0FPgSK40tHew40MSOA01sff8Yv9x2gENNbcyfMZ6vX1HDZ86eyMVnTaRylOblRYZKgS+R+Ki5nT1Hmtl9qJm3DzSx/cMm3v6wiQ+PtfZuMy5VwsKzJvDtRecyf8b46IoVKRAKfMkJd+doczvvfXScvYeb2XO4mb1Hmtl7pIW9h5tpPN7Ru21ZSRF1VWO4/JxJzD5jLLOnjmX2GWOZVpnScfMiIVLgy7C0d3bzYWMr+z8+zv6Pj/N+8JV5u7Wju3d7MzizchQ1k0dz/YXTqJ1cwVmTKqidXEHNpNHa0SqSBwp8OUlLeycHj7VxsKmNg02tJy0famoLbrfyUUvHp+47eUw508enmHPGWD4/ewrTJ4xi+vhR1E6uYMbE0aRKiyP4iUSkR94C38wqgL8H2oGX3P3xfD13UnV1O43HOzja3N779VFLxnJzO0f73G5u7/rU45QWG1Vjyqkal2LmpNHU10xgytgU08anmD5+FGeOH8W0ypQCXSTmsgp8M3sYuB446O7zMtoXA38HFAM/dvd7gD8A1rr7ejP7KaDAP42Orm6a2zppau3kWGsHx473fO/gWGtn8D3d3ti73JHe/ngHTW2dAz52RVkxEyrKmBh81VWNYUJFGZPGlDFlbIopY8uZMq6cKWNTjB9VSpEOexQZ8bId4T8C3A881tNgZsXAA8DVQAOw0czWAdXAm8Fmnx5GjlDd3U5rZxetHd20dnQFX91BWxctbV00t3fS3NZFS/A9fbuTlvYuPmnr7G3PXN/S1kV7V/dpn39sqoRxqVLGjSplXKqEGRNHB7dLGJsqZeLo0pOCfWJFGRNGl2k0LpJAWQW+u28ws5o+zZcCu9x9N4CZPQksJR3+1cBvOcWVtsxsJbASYObMmcOqq/F4B20dXbR1dtN2UhgH3zu7ON7eRWtnN22ZIR2sy9y+rXPg9W0d3YMK5b5SpUVUlJVQUV7C6LJiKspLGJsqYVplitFlJVSUp9sqyooZXVbCmFQJlaNKe4O8J+DHlJfoA0ciMmi5mMOfDryXcbsB+AxwH3C/mV0HrB/ozu6+ClgFUF9f78Mp4EsP/Jrdh5uHdB8zSJUUkyotIlVaTKq0mPKSnuUiJo0p+9T6nnWp0mJSJUUntZWXFpMqKc4I7xJGlxdTUaaQFpFo5G2nrbs3A18fzLbZng//9qvqaGnvorykmLKSok+FcSoI495gLi2irLhIx3yLSEHLReDvBzIvK1QdtA1atufD/9KC6uHcTUSkoOXi0y4bgVlmVmtmZcCNwLocPI+IiAxBVoFvZquBV4DZZtZgZsvdvRO4DXge2AascfctQ3zcJWa2qrGxMZvyREQkQ7ZH6dw0QPuzwLNZPK4ucSgiErJYnsBEI3wRkfDFMvDdfb27r6ysrIy6FBGRghHLwBcRkfDFMvA1pSMiEr5YBr6mdEREwhfLwBcRkfAp8EVEEiKWga85fBGR8MUy8DWHLyISvlgGvoiIhE+BLyKSELEMfM3hi4iEL5aBrzl8EZHwxTLwRUQkfAp8EZGEUOCLiCSEAl9EJCFiGfg6SkdEJHyxDHwdpSMiEr5YBr6IiIRPgS8ikhAKfBGRhFDgi4gkhAJfRCQhFPgiIgkRy8DXcfgiIuGLZeDrOHwRkfDFMvBFRCR8CnwRkYRQ4IuIJIQCX0QkIRT4IiIJocAXEUkIBb6ISEIo8EVEEkKBLyKSEHkLfDM728weMrO1+XpOERE5YVCBb2YPm9lBM3urT/tiM9thZrvM7M5TPYa773b35dkUKyIiw1cyyO0eAe4HHutpMLNi4AHgaqAB2Ghm64Bi4Pt97n+Lux/MuloRERm2QQW+u28ws5o+zZcCu9x9N4CZPQksdffvA9cPtyAzWwmsBJg5c+ZwH0ZERPrIZg5/OvBexu2GoK1fZjbJzB4EFpjZXQNt5+6r3L3e3eurqqqyKE9ERDINdkona+5+BLh1MNua2RJgSV1dXW6LEhFJkGxG+PuBGRm3q4O2rOl8+CIi4csm8DcCs8ys1szKgBuBdeGUJSIiYRvsYZmrgVeA2WbWYGbL3b0TuA14HtgGrHH3LWEUpUscioiEb7BH6dw0QPuzwLOhVpR+3PXA+vr6+hVhP7aISFLF8tQKGuGLiIQvloGvnbYiIuGLZeCLiEj4Yhn4mtIREQlfLANfUzoiIuGLZeCLiEj4FPgiIgkRy8DXHL6ISPhiGfiawxcRCV8sA19ERMKnwBcRSQgFvohIQsQy8LXTVkQkfLEMfO20FREJXywDX0REwqfAFxFJCAW+iEhCxDLwtdNWRCR8sQx87bQVEQlfLANfRETCp8AXEUkIBb6ISEIo8EVEEkKBLyKSEAp8EZGEiGXg6zh8EZHwxTLwdRy+iEj4Yhn4IiISPgW+iEhCKPBFRBJCgS8ikhAKfBGRhFDgi4gkhAJfRCQhSqIuQGQkONrcztX/69+jLkMS5NlvfY7S4nDH5Ap8kdP44oLpfNzSgeNRlyKSlbwFvpl9EbgOGAc85O7/mq/nFsnGJTUTuaRmYtRliGRtUP8vmNnDZnbQzN7q077YzHaY2S4zu/NUj+HuP3f3FcCtwFeHX7KIiAzHYEf4jwD3A4/1NJhZMfAAcDXQAGw0s3VAMfD9Pve/xd0PBsvfCe4nIiJ5NKjAd/cNZlbTp/lSYJe77wYwsyeBpe7+feD6vo9hZgbcAzzn7q8N9FxmthJYCTBz5szBlCciIoOQzS7g6cB7GbcbgraB3A4sAr5sZrcOtJG7r3L3enevr6qqyqI8ERHJlLedtu5+H3DfYLY1syXAkrq6utwWJSKSINmM8PcDMzJuVwdtWdP58EVEwpdN4G8EZplZrZmVATcC68IpS0REwjbYwzJXA68As82swcyWu3sncBvwPLANWOPuW8IoSpc4FBEJn7nH99ODZnYI2DfMu08GDodYTtjiXh+oxjDEvT6If41xrw/iVeNZ7t7vES+xDvxsmNkmd6+Puo6BxL0+UI1hiHt9EP8a414fjIwaQWfLFBFJDAW+iEhCFHLgr4q6gNOIe32gGsMQ9/og/jXGvT4YGTUW7hy+iIicrJBH+CIikkGBLyKSEAUX+EM5R38OnnuGmf3KzLaa2RYz+1bQPtHM/s3MdgbfJwTtZmb3BbW+YWYLMx7rT4Ltd5rZn4RcZ7GZ/YeZPR3crjWzV4M6fhp8chozKw9u7wrW12Q8xl1B+w4zuybk+sab2Voz225m28zs8jj1oZl9O/j9vmVmq80sFXUf9nfNijD7zMwuNrM3g/vcZ2YWUo33Br/nN8zsn81sfMa6fvtnoPf4QL+DbOrLWPcXZuZmNjm4HUkfZs3dC+aL9Ln43wHOBsqA14G5eXz+acDCYHks8DYwF/hb4M6g/U7gB8HyF4DnAAMuA14N2icCu4PvE4LlCSHWeQfwBPB0cHsNcGOw/CDwp8HyN4EHg+UbgZ8Gy3ODvi0HaoM+Lw6xvkeBbwTLZcD4uPQh6TPC7gFGZfTdsqj7ELgSWAi8ldEWWp8Bvwm2teC+14ZU4+8DJcHyDzJq7Ld/OMV7fKDfQTb1Be0zSJ9RYB8wOco+zPr1m+8nzOkPA5cDz2fcvgu4K8J6/oX0BWJ2ANOCtmnAjmD5H4GbMrbfEay/CfjHjPaTtsuypmrgBeDzwNPBi+9wxpuutw+DF/nlwXJJsJ317dfM7UKor5J0oFqf9lj0ISdOCz4x6JOngWvi0IdADSeHaSh9FqzbntF+0nbZ1Nhn3ZeAx4PlfvuHAd7jp3odZ1sfsBaYD+zlROBH1ofZfBXalM5Qz9GfM8G/7guAV4Gp7v5BsOpDYGqwPFC9ufw5fgT8JdAd3J4EfOzpcyP1fa7eOoL1jcH2uayvFjgE/MTS004/NrMKYtKH7r4f+CHwLvAB6T7ZTLz6sEdYfTY9WM5lrQC3kB75DqfGU72Oh83MlgL73f31Pqvi2oenVGiBHwtmNgb4GfDn7n4sc52n/7xHciysmV0PHHT3zVE8/yCVkP63+h/cfQHQTHo6olfEfTgBWEr6D9OZQAWwOIpahiLKPhsMM7sb6AQej7qWHmY2GvhvwP+IupawFFrg5+wc/YNlZqWkw/5xd38qaD5gZtOC9dOAnuv7DlRvrn6OK4AbzGwv8CTpaZ2/A8abWc/FcDKfq7eOYH0lcCSH9UF65NPg7q8Gt9eS/gMQlz5cBOxx90Pu3gE8Rbpf49SHPcLqs/3Bck5qNbNlpC+L+kfBH6bh1HiEgX8Hw3UO6T/srwfvmWrgNTM7Yxj15bQPBy3fc0i5/CI9OtxN+pfUs0Pn/Dw+v5G+0PuP+rTfy8k7z/42WL6Ok3f8/CZon0h6HntC8LUHmBhyrb/HiZ22/8TJO7u+GSz/Z07e4bgmWD6fk3eo7SbcnbYvA7OD5e8F/ReLPgQ+A2wBRgfP+Sjpy3dG3od8eg4/tD7j0zscvxBSjYuBrUBVn+367R9O8R4f6HeQTX191u3lxBx+ZH2Y1Wsk30+Y8x8ovff8bdJ78u/O83N/lvS/zW8Avw2+vkB6fvEFYCfwy4wXgAEPBLW+CdRnPNYtwK7g6+s5qPX3OBH4Zwcvxl3Bm6Y8aE8Ft3cF68/OuP/dQd07CPloA+AiYFPQjz8P3jix6UPgr4DtwFvA/w5CKdI+BFaT3qfQQfq/pOVh9hlQH/y87wD302enehY17iI9593zfnnwdP3DAO/xgX4H2dTXZ/1eTgR+JH2Y7ZdOrSAikhCFNocvIiIDUOCLiCSEAl9EJCEU+CIiCaHAFxFJCAW+SMDSZ+n8ZrB8ppmtjbomkTDpsEyRQHD+o6fdfV7UtYjkQsnpNxFJjHuAc8zst6Q/rHSeu88LPvr/RdLnzZlF+uRpZcAfA22kPzF51MzOIf1hnCqgBVjh7tvz/UOIDERTOiIn3Am84+4XAf+1z7p5wB8AlwB/A7R4+uRurwA3B9usAm5394uB/wL8fT6KFhksjfBFBudX7t4ENJlZI7A+aH8TuDA4Q+rvAP+UcSGj8vyXKTIwBb7I4LRlLHdn3O4m/T4qIn0+9ovyXJfIoGlKR+SEJtKXphwyT1/3YI+ZfQV6r3k6P8ziRLKlwBcJuPsR4NfBRazvHcZD/BGw3MxeJ30K5aVh1ieSLR2WKSKSEBrhi4gkhAJfRCQhFPgiIgmhwBcRSQgFvohIQijwRUQSQoEvIpIQ/x8q6Ro0odz3FAAAAABJRU5ErkJggg==\n",
+ "text/plain": [
+ "<Figure size 432x288 with 1 Axes>"
+ ]
+ },
+ "metadata": {
+ "needs_background": "light"
+ },
+ "output_type": "display_data"
+ }
+ ],
+ "source": []
+ },
+ {
+ "cell_type": "code",
+ "execution_count": null,
+ "id": "39509b1c-f5fb-4edd-847a-97757c2197df",
+ "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.6.4"
+ }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}
--
GitLab