From ecbaa9301be95b7d27de6b617a1c476a8ff2f94b Mon Sep 17 00:00:00 2001
From: David A Minton <daminton@purdue.edu>
Date: Thu, 17 Nov 2022 22:50:27 -0500
Subject: [PATCH] Updated the GR test Notebook

---
 .../whm_gr_test/swiftest_relativity.ipynb     | 941 +-----------------
 1 file changed, 34 insertions(+), 907 deletions(-)

diff --git a/examples/whm_gr_test/swiftest_relativity.ipynb b/examples/whm_gr_test/swiftest_relativity.ipynb
index 0cd73753f..6abb9524a 100644
--- a/examples/whm_gr_test/swiftest_relativity.ipynb
+++ b/examples/whm_gr_test/swiftest_relativity.ipynb
@@ -2,7 +2,7 @@
  "cells": [
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -15,913 +15,48 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Creating the Sun as a central body\n",
-      "Fetching ephemerides data for Mercury from JPL/Horizons\n",
-      "Fetching ephemerides data for Venus from JPL/Horizons\n",
-      "Fetching ephemerides data for Earth from JPL/Horizons\n",
-      "Fetching ephemerides data for Mars from JPL/Horizons\n",
-      "Fetching ephemerides data for Jupiter from JPL/Horizons\n",
-      "Fetching ephemerides data for Saturn from JPL/Horizons\n",
-      "Fetching ephemerides data for Uranus from JPL/Horizons\n",
-      "Fetching ephemerides data for Neptune from JPL/Horizons\n",
-      "Writing initial conditions to file init_cond.nc\n",
-      "Writing parameter inputs to file /home/daminton/git_debug/swiftest/examples/whm_gr_test/param.gr.in\n"
-     ]
-    },
-    {
-     "data": {
-      "text/html": [
-       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
-       "<defs>\n",
-       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
-       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
-       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
-       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
-       "</symbol>\n",
-       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
-       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
-       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
-       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
-       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
-       "</symbol>\n",
-       "</defs>\n",
-       "</svg>\n",
-       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
-       " *\n",
-       " */\n",
-       "\n",
-       ":root {\n",
-       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
-       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
-       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
-       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
-       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
-       "  --xr-background-color: var(--jp-layout-color0, white);\n",
-       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
-       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
-       "}\n",
-       "\n",
-       "html[theme=dark],\n",
-       "body.vscode-dark {\n",
-       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
-       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
-       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
-       "  --xr-border-color: #1F1F1F;\n",
-       "  --xr-disabled-color: #515151;\n",
-       "  --xr-background-color: #111111;\n",
-       "  --xr-background-color-row-even: #111111;\n",
-       "  --xr-background-color-row-odd: #313131;\n",
-       "}\n",
-       "\n",
-       ".xr-wrap {\n",
-       "  display: block;\n",
-       "  min-width: 300px;\n",
-       "  max-width: 700px;\n",
-       "}\n",
-       "\n",
-       ".xr-text-repr-fallback {\n",
-       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-header {\n",
-       "  padding-top: 6px;\n",
-       "  padding-bottom: 6px;\n",
-       "  margin-bottom: 4px;\n",
-       "  border-bottom: solid 1px var(--xr-border-color);\n",
-       "}\n",
-       "\n",
-       ".xr-header > div,\n",
-       ".xr-header > ul {\n",
-       "  display: inline;\n",
-       "  margin-top: 0;\n",
-       "  margin-bottom: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-obj-type,\n",
-       ".xr-array-name {\n",
-       "  margin-left: 2px;\n",
-       "  margin-right: 10px;\n",
-       "}\n",
-       "\n",
-       ".xr-obj-type {\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-sections {\n",
-       "  padding-left: 0 !important;\n",
-       "  display: grid;\n",
-       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
-       "}\n",
-       "\n",
-       ".xr-section-item {\n",
-       "  display: contents;\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input {\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input + label {\n",
-       "  color: var(--xr-disabled-color);\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input:enabled + label {\n",
-       "  cursor: pointer;\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input:enabled + label:hover {\n",
-       "  color: var(--xr-font-color0);\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary {\n",
-       "  grid-column: 1;\n",
-       "  color: var(--xr-font-color2);\n",
-       "  font-weight: 500;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary > span {\n",
-       "  display: inline-block;\n",
-       "  padding-left: 0.5em;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:disabled + label {\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in + label:before {\n",
-       "  display: inline-block;\n",
-       "  content: '►';\n",
-       "  font-size: 11px;\n",
-       "  width: 15px;\n",
-       "  text-align: center;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:disabled + label:before {\n",
-       "  color: var(--xr-disabled-color);\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:checked + label:before {\n",
-       "  content: '▼';\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:checked + label > span {\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary,\n",
-       ".xr-section-inline-details {\n",
-       "  padding-top: 4px;\n",
-       "  padding-bottom: 4px;\n",
-       "}\n",
-       "\n",
-       ".xr-section-inline-details {\n",
-       "  grid-column: 2 / -1;\n",
-       "}\n",
-       "\n",
-       ".xr-section-details {\n",
-       "  display: none;\n",
-       "  grid-column: 1 / -1;\n",
-       "  margin-bottom: 5px;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
-       "  display: contents;\n",
-       "}\n",
-       "\n",
-       ".xr-array-wrap {\n",
-       "  grid-column: 1 / -1;\n",
-       "  display: grid;\n",
-       "  grid-template-columns: 20px auto;\n",
-       "}\n",
-       "\n",
-       ".xr-array-wrap > label {\n",
-       "  grid-column: 1;\n",
-       "  vertical-align: top;\n",
-       "}\n",
-       "\n",
-       ".xr-preview {\n",
-       "  color: var(--xr-font-color3);\n",
-       "}\n",
-       "\n",
-       ".xr-array-preview,\n",
-       ".xr-array-data {\n",
-       "  padding: 0 5px !important;\n",
-       "  grid-column: 2;\n",
-       "}\n",
-       "\n",
-       ".xr-array-data,\n",
-       ".xr-array-in:checked ~ .xr-array-preview {\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-array-in:checked ~ .xr-array-data,\n",
-       ".xr-array-preview {\n",
-       "  display: inline-block;\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list {\n",
-       "  display: inline-block !important;\n",
-       "  list-style: none;\n",
-       "  padding: 0 !important;\n",
-       "  margin: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list li {\n",
-       "  display: inline-block;\n",
-       "  padding: 0;\n",
-       "  margin: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list:before {\n",
-       "  content: '(';\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list:after {\n",
-       "  content: ')';\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list li:not(:last-child):after {\n",
-       "  content: ',';\n",
-       "  padding-right: 5px;\n",
-       "}\n",
-       "\n",
-       ".xr-has-index {\n",
-       "  font-weight: bold;\n",
-       "}\n",
-       "\n",
-       ".xr-var-list,\n",
-       ".xr-var-item {\n",
-       "  display: contents;\n",
-       "}\n",
-       "\n",
-       ".xr-var-item > div,\n",
-       ".xr-var-item label,\n",
-       ".xr-var-item > .xr-var-name span {\n",
-       "  background-color: var(--xr-background-color-row-even);\n",
-       "  margin-bottom: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-var-item > .xr-var-name:hover span {\n",
-       "  padding-right: 5px;\n",
-       "}\n",
-       "\n",
-       ".xr-var-list > li:nth-child(odd) > div,\n",
-       ".xr-var-list > li:nth-child(odd) > label,\n",
-       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
-       "  background-color: var(--xr-background-color-row-odd);\n",
-       "}\n",
-       "\n",
-       ".xr-var-name {\n",
-       "  grid-column: 1;\n",
-       "}\n",
-       "\n",
-       ".xr-var-dims {\n",
-       "  grid-column: 2;\n",
-       "}\n",
-       "\n",
-       ".xr-var-dtype {\n",
-       "  grid-column: 3;\n",
-       "  text-align: right;\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-var-preview {\n",
-       "  grid-column: 4;\n",
-       "}\n",
-       "\n",
-       ".xr-var-name,\n",
-       ".xr-var-dims,\n",
-       ".xr-var-dtype,\n",
-       ".xr-preview,\n",
-       ".xr-attrs dt {\n",
-       "  white-space: nowrap;\n",
-       "  overflow: hidden;\n",
-       "  text-overflow: ellipsis;\n",
-       "  padding-right: 10px;\n",
-       "}\n",
-       "\n",
-       ".xr-var-name:hover,\n",
-       ".xr-var-dims:hover,\n",
-       ".xr-var-dtype:hover,\n",
-       ".xr-attrs dt:hover {\n",
-       "  overflow: visible;\n",
-       "  width: auto;\n",
-       "  z-index: 1;\n",
-       "}\n",
-       "\n",
-       ".xr-var-attrs,\n",
-       ".xr-var-data {\n",
-       "  display: none;\n",
-       "  background-color: var(--xr-background-color) !important;\n",
-       "  padding-bottom: 5px !important;\n",
-       "}\n",
-       "\n",
-       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
-       ".xr-var-data-in:checked ~ .xr-var-data {\n",
-       "  display: block;\n",
-       "}\n",
-       "\n",
-       ".xr-var-data > table {\n",
-       "  float: right;\n",
-       "}\n",
-       "\n",
-       ".xr-var-name span,\n",
-       ".xr-var-data,\n",
-       ".xr-attrs {\n",
-       "  padding-left: 25px !important;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs,\n",
-       ".xr-var-attrs,\n",
-       ".xr-var-data {\n",
-       "  grid-column: 1 / -1;\n",
-       "}\n",
-       "\n",
-       "dl.xr-attrs {\n",
-       "  padding: 0;\n",
-       "  margin: 0;\n",
-       "  display: grid;\n",
-       "  grid-template-columns: 125px auto;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dt,\n",
-       ".xr-attrs dd {\n",
-       "  padding: 0;\n",
-       "  margin: 0;\n",
-       "  float: left;\n",
-       "  padding-right: 10px;\n",
-       "  width: auto;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dt {\n",
-       "  font-weight: normal;\n",
-       "  grid-column: 1;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dt:hover span {\n",
-       "  display: inline-block;\n",
-       "  background: var(--xr-background-color);\n",
-       "  padding-right: 10px;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dd {\n",
-       "  grid-column: 2;\n",
-       "  white-space: pre-wrap;\n",
-       "  word-break: break-all;\n",
-       "}\n",
-       "\n",
-       ".xr-icon-database,\n",
-       ".xr-icon-file-text2 {\n",
-       "  display: inline-block;\n",
-       "  vertical-align: middle;\n",
-       "  width: 1em;\n",
-       "  height: 1.5em !important;\n",
-       "  stroke-width: 0;\n",
-       "  stroke: currentColor;\n",
-       "  fill: currentColor;\n",
-       "}\n",
-       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
-       "Dimensions:        (name: 9, time: 1)\n",
-       "Coordinates:\n",
-       "  * name           (name) &lt;U32 &#x27;Sun&#x27; &#x27;Mercury&#x27; &#x27;Venus&#x27; ... &#x27;Uranus&#x27; &#x27;Neptune&#x27;\n",
-       "  * time           (time) float64 0.0\n",
-       "Data variables: (12/14)\n",
-       "    particle_type  (name) &lt;U32 &#x27;Central Body&#x27; &#x27;Massive Body&#x27; ... &#x27;Massive Body&#x27;\n",
-       "    id             (name) int64 0 1 2 3 4 5 6 7 8\n",
-       "    Gmass          (time, name) float64 39.48 6.554e-06 ... 0.001724 0.002034\n",
-       "    radius         (time, name) float64 0.00465 1.631e-05 ... 0.0001646\n",
-       "    j2rp2          (time, name) float64 4.754e-12 nan nan nan ... nan nan nan\n",
-       "    j4rp4          (time, name) float64 -2.247e-18 nan nan nan ... nan nan nan\n",
-       "    ...             ...\n",
-       "    xhz            (time, name) float64 nan -0.002529 ... -0.0407 -0.7287\n",
-       "    vhx            (time, name) float64 nan -9.241 3.454 ... -1.315 -0.08755\n",
-       "    vhy            (time, name) float64 nan 7.755 6.477 ... 1.932 0.5372 1.149\n",
-       "    vhz            (time, name) float64 nan 1.481 -0.1103 ... 0.01905 -0.02168\n",
-       "    ntp            (time) int64 0\n",
-       "    npl            (time) int64 8</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-24d5c828-505c-4d47-a614-260f11ab60e8' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-24d5c828-505c-4d47-a614-260f11ab60e8' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>name</span>: 9</li><li><span class='xr-has-index'>time</span>: 1</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-e65fdf23-7fa6-463b-8c33-4d685235e3c1' class='xr-section-summary-in' type='checkbox'  checked><label for='section-e65fdf23-7fa6-463b-8c33-4d685235e3c1' class='xr-section-summary' >Coordinates: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>name</span></div><div class='xr-var-dims'>(name)</div><div class='xr-var-dtype'>&lt;U32</div><div class='xr-var-preview xr-preview'>&#x27;Sun&#x27; &#x27;Mercury&#x27; ... &#x27;Neptune&#x27;</div><input id='attrs-3dd7d645-9a3a-465a-9aec-e33a094581f5' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3dd7d645-9a3a-465a-9aec-e33a094581f5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7e27bd85-b58d-4960-bbc7-470a6aecc028' class='xr-var-data-in' type='checkbox'><label for='data-7e27bd85-b58d-4960-bbc7-470a6aecc028' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;Sun&#x27;, &#x27;Mercury&#x27;, &#x27;Venus&#x27;, &#x27;Earth&#x27;, &#x27;Mars&#x27;, &#x27;Jupiter&#x27;, &#x27;Saturn&#x27;,\n",
-       "       &#x27;Uranus&#x27;, &#x27;Neptune&#x27;], dtype=&#x27;&lt;U32&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0</div><input id='attrs-b4a3d6b1-2b75-4297-bf70-48da66a4b964' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b4a3d6b1-2b75-4297-bf70-48da66a4b964' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a71efd79-b820-49ee-8910-7abffc2245d4' class='xr-var-data-in' type='checkbox'><label for='data-a71efd79-b820-49ee-8910-7abffc2245d4' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([0.])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-3e8ed869-d3c1-46b5-936f-5868089ab757' class='xr-section-summary-in' type='checkbox'  checked><label for='section-3e8ed869-d3c1-46b5-936f-5868089ab757' class='xr-section-summary' >Data variables: <span>(14)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>particle_type</span></div><div class='xr-var-dims'>(name)</div><div class='xr-var-dtype'>&lt;U32</div><div class='xr-var-preview xr-preview'>&#x27;Central Body&#x27; ... &#x27;Massive Body&#x27;</div><input id='attrs-5c29c659-755d-4b25-9ccf-5ddca9477b90' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-5c29c659-755d-4b25-9ccf-5ddca9477b90' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b362f92d-189f-4cec-9a4b-e14e53e1fd71' class='xr-var-data-in' type='checkbox'><label for='data-b362f92d-189f-4cec-9a4b-e14e53e1fd71' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;Central Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;,\n",
-       "       &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;,\n",
-       "       &#x27;Massive Body&#x27;], dtype=&#x27;&lt;U32&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>id</span></div><div class='xr-var-dims'>(name)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0 1 2 3 4 5 6 7 8</div><input id='attrs-17057a19-fc7b-4f57-a907-307688a3d827' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-17057a19-fc7b-4f57-a907-307688a3d827' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-88d0c856-39ee-439c-8c3a-936f9e89c31e' class='xr-var-data-in' type='checkbox'><label for='data-88d0c856-39ee-439c-8c3a-936f9e89c31e' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([0, 1, 2, 3, 4, 5, 6, 7, 8])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>Gmass</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>39.48 6.554e-06 ... 0.002034</div><input id='attrs-6d157615-df32-4456-99e1-3025ad081179' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6d157615-df32-4456-99e1-3025ad081179' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-295255e7-c548-4cdc-85a4-68f371490789' class='xr-var-data-in' type='checkbox'><label for='data-295255e7-c548-4cdc-85a4-68f371490789' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[3.94769264e+01, 6.55370981e-06, 9.66331340e-05, 1.20026936e-04,\n",
-       "        1.27398020e-05, 3.76922511e-02, 1.12858998e-02, 1.72365895e-03,\n",
-       "        2.03361005e-03]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>radius</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.00465 1.631e-05 ... 0.0001646</div><input id='attrs-f0d7810e-7d03-44bd-8ed8-fc7457626c50' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f0d7810e-7d03-44bd-8ed8-fc7457626c50' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8001a5f4-8347-4cf3-a8c0-95abd2255899' class='xr-var-data-in' type='checkbox'><label for='data-8001a5f4-8347-4cf3-a8c0-95abd2255899' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[4.65046726e-03, 1.63063818e-05, 4.04537843e-05, 4.25875607e-05,\n",
-       "        2.26574081e-05, 4.67326170e-04, 3.89256877e-04, 1.69534499e-04,\n",
-       "        1.64587904e-04]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>j2rp2</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>4.754e-12 nan nan ... nan nan nan</div><input id='attrs-33b47d88-61a7-4217-8ea5-d53d8d634250' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-33b47d88-61a7-4217-8ea5-d53d8d634250' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-63782b0e-891a-4ba3-8c91-6c14e0ec41f0' class='xr-var-data-in' type='checkbox'><label for='data-63782b0e-891a-4ba3-8c91-6c14e0ec41f0' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[4.75358069e-12,            nan,            nan,            nan,\n",
-       "                   nan,            nan,            nan,            nan,\n",
-       "                   nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>j4rp4</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-2.247e-18 nan nan ... nan nan nan</div><input id='attrs-61761cd0-dc1b-4d13-92b8-024d302fa9a4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-61761cd0-dc1b-4d13-92b8-024d302fa9a4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-72df4131-dc97-4743-8ada-7685342dfb6d' class='xr-var-data-in' type='checkbox'><label for='data-72df4131-dc97-4743-8ada-7685342dfb6d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[-2.2473968e-18,            nan,            nan,            nan,\n",
-       "                   nan,            nan,            nan,            nan,\n",
-       "                   nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>xhx</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 0.2253 0.6401 ... 8.179 29.79</div><input id='attrs-f8306230-bcf9-4567-aa23-4c5325265873' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f8306230-bcf9-4567-aa23-4c5325265873' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ea9f1018-9788-43ae-8173-9f3a7211d02c' class='xr-var-data-in' type='checkbox'><label for='data-ea9f1018-9788-43ae-8173-9f3a7211d02c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan,  0.22527007,  0.64013386, -0.78195044, -1.65001718,\n",
-       "        -4.54078501,  8.901068  ,  8.17855335, 29.79419863]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>xhy</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 0.2219 -0.344 ... 17.59 2.045</div><input id='attrs-7fb81775-bdd5-4e9f-bb12-d24f02d7ab5b' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-7fb81775-bdd5-4e9f-bb12-d24f02d7ab5b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-399a26fb-56ce-4398-9079-ac72c4f0b820' class='xr-var-data-in' type='checkbox'><label for='data-399a26fb-56ce-4398-9079-ac72c4f0b820' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan,  0.22185516, -0.34396282, -0.63468545, -0.02334143,\n",
-       "         2.86428737,  2.90284887, 17.59478942,  2.04463099]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>xhz</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan -0.002529 ... -0.0407 -0.7287</div><input id='attrs-33a79169-b3ba-44d4-87ec-56af9d188c8f' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-33a79169-b3ba-44d4-87ec-56af9d188c8f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4a1f3a63-a2c5-4ce3-bf66-edc66e22eadf' class='xr-var-data-in' type='checkbox'><label for='data-4a1f3a63-a2c5-4ce3-bf66-edc66e22eadf' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[            nan, -2.52922505e-03, -4.16625374e-02,\n",
-       "         4.44514635e-05,  3.99643397e-02,  8.96961900e-02,\n",
-       "        -4.04874151e-01, -4.06951825e-02, -7.28719479e-01]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vhx</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan -9.241 ... -1.315 -0.08755</div><input id='attrs-d5a383c3-7d7c-400a-9d70-36fadc5ad33e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-d5a383c3-7d7c-400a-9d70-36fadc5ad33e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-98014e22-9102-4da4-a8d7-8dc2ceb91b20' class='xr-var-data-in' type='checkbox'><label for='data-98014e22-9102-4da4-a8d7-8dc2ceb91b20' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan, -9.24105838,  3.45369085,  3.8578492 ,  0.26241127,\n",
-       "        -1.50370126, -0.74393964, -1.31492768, -0.0875473 ]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vhy</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 7.755 6.477 ... 0.5372 1.149</div><input id='attrs-ca1983fe-0960-4511-ac3b-7057303cc7e4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ca1983fe-0960-4511-ac3b-7057303cc7e4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-74bdf30c-f36b-49ac-b4f0-f67a13abfa65' class='xr-var-data-in' type='checkbox'><label for='data-74bdf30c-f36b-49ac-b4f0-f67a13abfa65' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan,  7.754686  ,  6.47714891, -4.90267373, -4.67368878,\n",
-       "        -2.20463162,  1.93168513,  0.5372287 ,  1.14922255]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vhz</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 1.481 ... 0.01905 -0.02168</div><input id='attrs-60504126-9290-4dcc-b5bf-a6014885a51c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-60504126-9290-4dcc-b5bf-a6014885a51c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-070b3756-82db-4382-afc5-c352f1bd5d42' class='xr-var-data-in' type='checkbox'><label for='data-070b3756-82db-4382-afc5-c352f1bd5d42' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[            nan,  1.48128538e+00, -1.10257596e-01,\n",
-       "         2.88870085e-04, -1.04392391e-01,  4.28033174e-02,\n",
-       "        -4.02921477e-03,  1.90517453e-02, -2.16786846e-02]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ntp</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-a68d3c53-d50b-44bb-930a-e064397ec812' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a68d3c53-d50b-44bb-930a-e064397ec812' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cc7e56cd-1660-4789-be77-c29e8d9ebdaf' class='xr-var-data-in' type='checkbox'><label for='data-cc7e56cd-1660-4789-be77-c29e8d9ebdaf' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([0])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>npl</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>8</div><input id='attrs-4938babb-5d6e-40b1-9df6-1e9430d96b39' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4938babb-5d6e-40b1-9df6-1e9430d96b39' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a11f9f42-6b12-4529-ac26-740858d89b42' class='xr-var-data-in' type='checkbox'><label for='data-a11f9f42-6b12-4529-ac26-740858d89b42' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([8])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-1e288b47-016d-4588-adcc-ecd025bf5071' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-1e288b47-016d-4588-adcc-ecd025bf5071' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
-      ],
-      "text/plain": [
-       "<xarray.Dataset>\n",
-       "Dimensions:        (name: 9, time: 1)\n",
-       "Coordinates:\n",
-       "  * name           (name) <U32 'Sun' 'Mercury' 'Venus' ... 'Uranus' 'Neptune'\n",
-       "  * time           (time) float64 0.0\n",
-       "Data variables: (12/14)\n",
-       "    particle_type  (name) <U32 'Central Body' 'Massive Body' ... 'Massive Body'\n",
-       "    id             (name) int64 0 1 2 3 4 5 6 7 8\n",
-       "    Gmass          (time, name) float64 39.48 6.554e-06 ... 0.001724 0.002034\n",
-       "    radius         (time, name) float64 0.00465 1.631e-05 ... 0.0001646\n",
-       "    j2rp2          (time, name) float64 4.754e-12 nan nan nan ... nan nan nan\n",
-       "    j4rp4          (time, name) float64 -2.247e-18 nan nan nan ... nan nan nan\n",
-       "    ...             ...\n",
-       "    xhz            (time, name) float64 nan -0.002529 ... -0.0407 -0.7287\n",
-       "    vhx            (time, name) float64 nan -9.241 3.454 ... -1.315 -0.08755\n",
-       "    vhy            (time, name) float64 nan 7.755 6.477 ... 1.932 0.5372 1.149\n",
-       "    vhz            (time, name) float64 nan 1.481 -0.1103 ... 0.01905 -0.02168\n",
-       "    ntp            (time) int64 0\n",
-       "    npl            (time) int64 8"
-      ]
-     },
-     "execution_count": 2,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
-    "sim_gr = swiftest.Simulation(tstop=1000.0, dt=0.005, tstep_out=1.0, integrator=\"whm\", init_cond_format=\"XV\", general_relativity=True, param_file=\"param.gr.in\", output_file_name=\"bin.gr.nc\")\n",
+    "sim_gr = swiftest.Simulation(param_file=\"param.gr.in\", output_file_name=\"bin.gr.nc\")\n",
     "sim_gr.add_solar_system_body([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"])"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stdout",
-     "output_type": "stream",
-     "text": [
-      "Creating the Sun as a central body\n",
-      "Fetching ephemerides data for Mercury from JPL/Horizons\n",
-      "Fetching ephemerides data for Venus from JPL/Horizons\n",
-      "Fetching ephemerides data for Earth from JPL/Horizons\n",
-      "Fetching ephemerides data for Mars from JPL/Horizons\n",
-      "Fetching ephemerides data for Jupiter from JPL/Horizons\n",
-      "Fetching ephemerides data for Saturn from JPL/Horizons\n",
-      "Fetching ephemerides data for Uranus from JPL/Horizons\n",
-      "Fetching ephemerides data for Neptune from JPL/Horizons\n",
-      "Writing initial conditions to file init_cond.nc\n",
-      "Writing parameter inputs to file /home/daminton/git_debug/swiftest/examples/whm_gr_test/param.nogr.in\n"
-     ]
-    },
-    {
-     "data": {
-      "text/html": [
-       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
-       "<defs>\n",
-       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
-       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
-       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
-       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
-       "</symbol>\n",
-       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
-       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
-       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
-       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
-       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
-       "</symbol>\n",
-       "</defs>\n",
-       "</svg>\n",
-       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
-       " *\n",
-       " */\n",
-       "\n",
-       ":root {\n",
-       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
-       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
-       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
-       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
-       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
-       "  --xr-background-color: var(--jp-layout-color0, white);\n",
-       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
-       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
-       "}\n",
-       "\n",
-       "html[theme=dark],\n",
-       "body.vscode-dark {\n",
-       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
-       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
-       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
-       "  --xr-border-color: #1F1F1F;\n",
-       "  --xr-disabled-color: #515151;\n",
-       "  --xr-background-color: #111111;\n",
-       "  --xr-background-color-row-even: #111111;\n",
-       "  --xr-background-color-row-odd: #313131;\n",
-       "}\n",
-       "\n",
-       ".xr-wrap {\n",
-       "  display: block;\n",
-       "  min-width: 300px;\n",
-       "  max-width: 700px;\n",
-       "}\n",
-       "\n",
-       ".xr-text-repr-fallback {\n",
-       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-header {\n",
-       "  padding-top: 6px;\n",
-       "  padding-bottom: 6px;\n",
-       "  margin-bottom: 4px;\n",
-       "  border-bottom: solid 1px var(--xr-border-color);\n",
-       "}\n",
-       "\n",
-       ".xr-header > div,\n",
-       ".xr-header > ul {\n",
-       "  display: inline;\n",
-       "  margin-top: 0;\n",
-       "  margin-bottom: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-obj-type,\n",
-       ".xr-array-name {\n",
-       "  margin-left: 2px;\n",
-       "  margin-right: 10px;\n",
-       "}\n",
-       "\n",
-       ".xr-obj-type {\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-sections {\n",
-       "  padding-left: 0 !important;\n",
-       "  display: grid;\n",
-       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
-       "}\n",
-       "\n",
-       ".xr-section-item {\n",
-       "  display: contents;\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input {\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input + label {\n",
-       "  color: var(--xr-disabled-color);\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input:enabled + label {\n",
-       "  cursor: pointer;\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-section-item input:enabled + label:hover {\n",
-       "  color: var(--xr-font-color0);\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary {\n",
-       "  grid-column: 1;\n",
-       "  color: var(--xr-font-color2);\n",
-       "  font-weight: 500;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary > span {\n",
-       "  display: inline-block;\n",
-       "  padding-left: 0.5em;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:disabled + label {\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in + label:before {\n",
-       "  display: inline-block;\n",
-       "  content: '►';\n",
-       "  font-size: 11px;\n",
-       "  width: 15px;\n",
-       "  text-align: center;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:disabled + label:before {\n",
-       "  color: var(--xr-disabled-color);\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:checked + label:before {\n",
-       "  content: '▼';\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:checked + label > span {\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary,\n",
-       ".xr-section-inline-details {\n",
-       "  padding-top: 4px;\n",
-       "  padding-bottom: 4px;\n",
-       "}\n",
-       "\n",
-       ".xr-section-inline-details {\n",
-       "  grid-column: 2 / -1;\n",
-       "}\n",
-       "\n",
-       ".xr-section-details {\n",
-       "  display: none;\n",
-       "  grid-column: 1 / -1;\n",
-       "  margin-bottom: 5px;\n",
-       "}\n",
-       "\n",
-       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
-       "  display: contents;\n",
-       "}\n",
-       "\n",
-       ".xr-array-wrap {\n",
-       "  grid-column: 1 / -1;\n",
-       "  display: grid;\n",
-       "  grid-template-columns: 20px auto;\n",
-       "}\n",
-       "\n",
-       ".xr-array-wrap > label {\n",
-       "  grid-column: 1;\n",
-       "  vertical-align: top;\n",
-       "}\n",
-       "\n",
-       ".xr-preview {\n",
-       "  color: var(--xr-font-color3);\n",
-       "}\n",
-       "\n",
-       ".xr-array-preview,\n",
-       ".xr-array-data {\n",
-       "  padding: 0 5px !important;\n",
-       "  grid-column: 2;\n",
-       "}\n",
-       "\n",
-       ".xr-array-data,\n",
-       ".xr-array-in:checked ~ .xr-array-preview {\n",
-       "  display: none;\n",
-       "}\n",
-       "\n",
-       ".xr-array-in:checked ~ .xr-array-data,\n",
-       ".xr-array-preview {\n",
-       "  display: inline-block;\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list {\n",
-       "  display: inline-block !important;\n",
-       "  list-style: none;\n",
-       "  padding: 0 !important;\n",
-       "  margin: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list li {\n",
-       "  display: inline-block;\n",
-       "  padding: 0;\n",
-       "  margin: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list:before {\n",
-       "  content: '(';\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list:after {\n",
-       "  content: ')';\n",
-       "}\n",
-       "\n",
-       ".xr-dim-list li:not(:last-child):after {\n",
-       "  content: ',';\n",
-       "  padding-right: 5px;\n",
-       "}\n",
-       "\n",
-       ".xr-has-index {\n",
-       "  font-weight: bold;\n",
-       "}\n",
-       "\n",
-       ".xr-var-list,\n",
-       ".xr-var-item {\n",
-       "  display: contents;\n",
-       "}\n",
-       "\n",
-       ".xr-var-item > div,\n",
-       ".xr-var-item label,\n",
-       ".xr-var-item > .xr-var-name span {\n",
-       "  background-color: var(--xr-background-color-row-even);\n",
-       "  margin-bottom: 0;\n",
-       "}\n",
-       "\n",
-       ".xr-var-item > .xr-var-name:hover span {\n",
-       "  padding-right: 5px;\n",
-       "}\n",
-       "\n",
-       ".xr-var-list > li:nth-child(odd) > div,\n",
-       ".xr-var-list > li:nth-child(odd) > label,\n",
-       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
-       "  background-color: var(--xr-background-color-row-odd);\n",
-       "}\n",
-       "\n",
-       ".xr-var-name {\n",
-       "  grid-column: 1;\n",
-       "}\n",
-       "\n",
-       ".xr-var-dims {\n",
-       "  grid-column: 2;\n",
-       "}\n",
-       "\n",
-       ".xr-var-dtype {\n",
-       "  grid-column: 3;\n",
-       "  text-align: right;\n",
-       "  color: var(--xr-font-color2);\n",
-       "}\n",
-       "\n",
-       ".xr-var-preview {\n",
-       "  grid-column: 4;\n",
-       "}\n",
-       "\n",
-       ".xr-var-name,\n",
-       ".xr-var-dims,\n",
-       ".xr-var-dtype,\n",
-       ".xr-preview,\n",
-       ".xr-attrs dt {\n",
-       "  white-space: nowrap;\n",
-       "  overflow: hidden;\n",
-       "  text-overflow: ellipsis;\n",
-       "  padding-right: 10px;\n",
-       "}\n",
-       "\n",
-       ".xr-var-name:hover,\n",
-       ".xr-var-dims:hover,\n",
-       ".xr-var-dtype:hover,\n",
-       ".xr-attrs dt:hover {\n",
-       "  overflow: visible;\n",
-       "  width: auto;\n",
-       "  z-index: 1;\n",
-       "}\n",
-       "\n",
-       ".xr-var-attrs,\n",
-       ".xr-var-data {\n",
-       "  display: none;\n",
-       "  background-color: var(--xr-background-color) !important;\n",
-       "  padding-bottom: 5px !important;\n",
-       "}\n",
-       "\n",
-       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
-       ".xr-var-data-in:checked ~ .xr-var-data {\n",
-       "  display: block;\n",
-       "}\n",
-       "\n",
-       ".xr-var-data > table {\n",
-       "  float: right;\n",
-       "}\n",
-       "\n",
-       ".xr-var-name span,\n",
-       ".xr-var-data,\n",
-       ".xr-attrs {\n",
-       "  padding-left: 25px !important;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs,\n",
-       ".xr-var-attrs,\n",
-       ".xr-var-data {\n",
-       "  grid-column: 1 / -1;\n",
-       "}\n",
-       "\n",
-       "dl.xr-attrs {\n",
-       "  padding: 0;\n",
-       "  margin: 0;\n",
-       "  display: grid;\n",
-       "  grid-template-columns: 125px auto;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dt,\n",
-       ".xr-attrs dd {\n",
-       "  padding: 0;\n",
-       "  margin: 0;\n",
-       "  float: left;\n",
-       "  padding-right: 10px;\n",
-       "  width: auto;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dt {\n",
-       "  font-weight: normal;\n",
-       "  grid-column: 1;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dt:hover span {\n",
-       "  display: inline-block;\n",
-       "  background: var(--xr-background-color);\n",
-       "  padding-right: 10px;\n",
-       "}\n",
-       "\n",
-       ".xr-attrs dd {\n",
-       "  grid-column: 2;\n",
-       "  white-space: pre-wrap;\n",
-       "  word-break: break-all;\n",
-       "}\n",
-       "\n",
-       ".xr-icon-database,\n",
-       ".xr-icon-file-text2 {\n",
-       "  display: inline-block;\n",
-       "  vertical-align: middle;\n",
-       "  width: 1em;\n",
-       "  height: 1.5em !important;\n",
-       "  stroke-width: 0;\n",
-       "  stroke: currentColor;\n",
-       "  fill: currentColor;\n",
-       "}\n",
-       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
-       "Dimensions:        (name: 9, time: 1)\n",
-       "Coordinates:\n",
-       "  * name           (name) &lt;U32 &#x27;Sun&#x27; &#x27;Mercury&#x27; &#x27;Venus&#x27; ... &#x27;Uranus&#x27; &#x27;Neptune&#x27;\n",
-       "  * time           (time) float64 0.0\n",
-       "Data variables: (12/14)\n",
-       "    particle_type  (name) &lt;U32 &#x27;Central Body&#x27; &#x27;Massive Body&#x27; ... &#x27;Massive Body&#x27;\n",
-       "    id             (name) int64 0 1 2 3 4 5 6 7 8\n",
-       "    Gmass          (time, name) float64 39.48 6.554e-06 ... 0.001724 0.002034\n",
-       "    radius         (time, name) float64 0.00465 1.631e-05 ... 0.0001646\n",
-       "    j2rp2          (time, name) float64 4.754e-12 nan nan nan ... nan nan nan\n",
-       "    j4rp4          (time, name) float64 -2.247e-18 nan nan nan ... nan nan nan\n",
-       "    ...             ...\n",
-       "    xhz            (time, name) float64 nan -0.002529 ... -0.0407 -0.7287\n",
-       "    vhx            (time, name) float64 nan -9.241 3.454 ... -1.315 -0.08755\n",
-       "    vhy            (time, name) float64 nan 7.755 6.477 ... 1.932 0.5372 1.149\n",
-       "    vhz            (time, name) float64 nan 1.481 -0.1103 ... 0.01905 -0.02168\n",
-       "    ntp            (time) int64 0\n",
-       "    npl            (time) int64 8</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-d769bb96-7cd9-40d9-ace7-2e7db8872d8d' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-d769bb96-7cd9-40d9-ace7-2e7db8872d8d' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>name</span>: 9</li><li><span class='xr-has-index'>time</span>: 1</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-d6545026-ba8c-4baa-b923-625ca0773cfb' class='xr-section-summary-in' type='checkbox'  checked><label for='section-d6545026-ba8c-4baa-b923-625ca0773cfb' class='xr-section-summary' >Coordinates: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>name</span></div><div class='xr-var-dims'>(name)</div><div class='xr-var-dtype'>&lt;U32</div><div class='xr-var-preview xr-preview'>&#x27;Sun&#x27; &#x27;Mercury&#x27; ... &#x27;Neptune&#x27;</div><input id='attrs-22fde872-9ec9-4387-b835-005227e39975' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-22fde872-9ec9-4387-b835-005227e39975' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-21b41f73-2440-4af6-a1cb-9b4ee0861917' class='xr-var-data-in' type='checkbox'><label for='data-21b41f73-2440-4af6-a1cb-9b4ee0861917' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;Sun&#x27;, &#x27;Mercury&#x27;, &#x27;Venus&#x27;, &#x27;Earth&#x27;, &#x27;Mars&#x27;, &#x27;Jupiter&#x27;, &#x27;Saturn&#x27;,\n",
-       "       &#x27;Uranus&#x27;, &#x27;Neptune&#x27;], dtype=&#x27;&lt;U32&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0</div><input id='attrs-c44ddd8d-965e-47ff-abd6-85dce4a3d967' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c44ddd8d-965e-47ff-abd6-85dce4a3d967' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a953f397-3b70-4dd9-9522-11a1634a17ec' class='xr-var-data-in' type='checkbox'><label for='data-a953f397-3b70-4dd9-9522-11a1634a17ec' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([0.])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-5511f2b8-b9a0-44b2-b3fa-2326b2a6f54e' class='xr-section-summary-in' type='checkbox'  checked><label for='section-5511f2b8-b9a0-44b2-b3fa-2326b2a6f54e' class='xr-section-summary' >Data variables: <span>(14)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>particle_type</span></div><div class='xr-var-dims'>(name)</div><div class='xr-var-dtype'>&lt;U32</div><div class='xr-var-preview xr-preview'>&#x27;Central Body&#x27; ... &#x27;Massive Body&#x27;</div><input id='attrs-19f1454e-150f-426c-88f3-a87452c2e254' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-19f1454e-150f-426c-88f3-a87452c2e254' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5b520aa1-ec07-43b1-bbf1-454ca1bad47a' class='xr-var-data-in' type='checkbox'><label for='data-5b520aa1-ec07-43b1-bbf1-454ca1bad47a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;Central Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;,\n",
-       "       &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;, &#x27;Massive Body&#x27;,\n",
-       "       &#x27;Massive Body&#x27;], dtype=&#x27;&lt;U32&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>id</span></div><div class='xr-var-dims'>(name)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0 1 2 3 4 5 6 7 8</div><input id='attrs-65489649-75ff-480d-b219-68db344d9e79' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-65489649-75ff-480d-b219-68db344d9e79' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0c72e67e-d075-4975-9f1d-af473eca1c75' class='xr-var-data-in' type='checkbox'><label for='data-0c72e67e-d075-4975-9f1d-af473eca1c75' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([0, 1, 2, 3, 4, 5, 6, 7, 8])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>Gmass</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>39.48 6.554e-06 ... 0.002034</div><input id='attrs-14bef83c-9120-4d67-a651-084feadef303' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-14bef83c-9120-4d67-a651-084feadef303' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ffb3617e-192c-4318-81e3-8fec2cc3c6d0' class='xr-var-data-in' type='checkbox'><label for='data-ffb3617e-192c-4318-81e3-8fec2cc3c6d0' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[3.94769264e+01, 6.55370981e-06, 9.66331340e-05, 1.20026936e-04,\n",
-       "        1.27398020e-05, 3.76922511e-02, 1.12858998e-02, 1.72365895e-03,\n",
-       "        2.03361005e-03]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>radius</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.00465 1.631e-05 ... 0.0001646</div><input id='attrs-b5014045-573a-46e1-a41c-777c9f36915e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b5014045-573a-46e1-a41c-777c9f36915e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d814da83-4203-4116-a01c-984dde892660' class='xr-var-data-in' type='checkbox'><label for='data-d814da83-4203-4116-a01c-984dde892660' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[4.65046726e-03, 1.63063818e-05, 4.04537843e-05, 4.25875607e-05,\n",
-       "        2.26574081e-05, 4.67326170e-04, 3.89256877e-04, 1.69534499e-04,\n",
-       "        1.64587904e-04]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>j2rp2</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>4.754e-12 nan nan ... nan nan nan</div><input id='attrs-c86d4393-405c-4b5d-8485-cc396d181eba' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c86d4393-405c-4b5d-8485-cc396d181eba' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-185a8e9e-f944-4754-ae8f-567702fab0ed' class='xr-var-data-in' type='checkbox'><label for='data-185a8e9e-f944-4754-ae8f-567702fab0ed' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[4.75358069e-12,            nan,            nan,            nan,\n",
-       "                   nan,            nan,            nan,            nan,\n",
-       "                   nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>j4rp4</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-2.247e-18 nan nan ... nan nan nan</div><input id='attrs-aed81728-097c-4713-a51f-8753945814e4' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-aed81728-097c-4713-a51f-8753945814e4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-21e0b433-6fad-4838-aff5-bcd405de3348' class='xr-var-data-in' type='checkbox'><label for='data-21e0b433-6fad-4838-aff5-bcd405de3348' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[-2.2473968e-18,            nan,            nan,            nan,\n",
-       "                   nan,            nan,            nan,            nan,\n",
-       "                   nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>xhx</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 0.2253 0.6401 ... 8.179 29.79</div><input id='attrs-e7868fbf-d7b7-4468-8072-a49631e37c17' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e7868fbf-d7b7-4468-8072-a49631e37c17' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8b36c8f6-afe6-4e4c-b232-1728f99923cc' class='xr-var-data-in' type='checkbox'><label for='data-8b36c8f6-afe6-4e4c-b232-1728f99923cc' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan,  0.22527007,  0.64013386, -0.78195044, -1.65001718,\n",
-       "        -4.54078501,  8.901068  ,  8.17855335, 29.79419863]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>xhy</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 0.2219 -0.344 ... 17.59 2.045</div><input id='attrs-4357a954-e884-466f-9182-ee6a624e598a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-4357a954-e884-466f-9182-ee6a624e598a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-69a46900-227f-4164-a4de-8b5d75aba708' class='xr-var-data-in' type='checkbox'><label for='data-69a46900-227f-4164-a4de-8b5d75aba708' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan,  0.22185516, -0.34396282, -0.63468545, -0.02334143,\n",
-       "         2.86428737,  2.90284887, 17.59478942,  2.04463099]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>xhz</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan -0.002529 ... -0.0407 -0.7287</div><input id='attrs-f75fd41a-e565-493f-871e-cd0c72d775d1' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f75fd41a-e565-493f-871e-cd0c72d775d1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f0c4e898-8283-4341-bd38-1c3a027b5cee' class='xr-var-data-in' type='checkbox'><label for='data-f0c4e898-8283-4341-bd38-1c3a027b5cee' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[            nan, -2.52922505e-03, -4.16625374e-02,\n",
-       "         4.44514635e-05,  3.99643397e-02,  8.96961900e-02,\n",
-       "        -4.04874151e-01, -4.06951825e-02, -7.28719479e-01]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vhx</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan -9.241 ... -1.315 -0.08755</div><input id='attrs-08109c76-cabf-4685-948e-5daea05d8d48' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-08109c76-cabf-4685-948e-5daea05d8d48' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-300489a9-057c-4a22-aaf7-73c009b8309b' class='xr-var-data-in' type='checkbox'><label for='data-300489a9-057c-4a22-aaf7-73c009b8309b' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan, -9.24105838,  3.45369085,  3.8578492 ,  0.26241127,\n",
-       "        -1.50370126, -0.74393964, -1.31492768, -0.0875473 ]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vhy</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 7.755 6.477 ... 0.5372 1.149</div><input id='attrs-c196ec7b-eb87-417b-bbcc-3f9c2eb0322a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c196ec7b-eb87-417b-bbcc-3f9c2eb0322a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-dcd832da-96e1-450b-88af-6f8d4d3f625b' class='xr-var-data-in' type='checkbox'><label for='data-dcd832da-96e1-450b-88af-6f8d4d3f625b' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[        nan,  7.754686  ,  6.47714891, -4.90267373, -4.67368878,\n",
-       "        -2.20463162,  1.93168513,  0.5372287 ,  1.14922255]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vhz</span></div><div class='xr-var-dims'>(time, name)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan 1.481 ... 0.01905 -0.02168</div><input id='attrs-faf7ace0-e940-4bd4-a9be-9889a0d83d14' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-faf7ace0-e940-4bd4-a9be-9889a0d83d14' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-58cadf78-c580-4381-9c72-9ac771a58acd' class='xr-var-data-in' type='checkbox'><label for='data-58cadf78-c580-4381-9c72-9ac771a58acd' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([[            nan,  1.48128538e+00, -1.10257596e-01,\n",
-       "         2.88870085e-04, -1.04392391e-01,  4.28033174e-02,\n",
-       "        -4.02921477e-03,  1.90517453e-02, -2.16786846e-02]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ntp</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-7dfc0453-f75b-4f50-80c3-46f8caa8161e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-7dfc0453-f75b-4f50-80c3-46f8caa8161e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-59609ae2-a58b-4815-94d1-538f74371c75' class='xr-var-data-in' type='checkbox'><label for='data-59609ae2-a58b-4815-94d1-538f74371c75' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([0])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>npl</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>8</div><input id='attrs-8a5beeb7-7143-41e4-8aaf-456a76a3ab36' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-8a5beeb7-7143-41e4-8aaf-456a76a3ab36' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a24595ea-29de-428d-8953-f77aafee4bbe' class='xr-var-data-in' type='checkbox'><label for='data-a24595ea-29de-428d-8953-f77aafee4bbe' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([8])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-4e62ee77-dfb0-4f29-aa3a-89c485de33e0' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-4e62ee77-dfb0-4f29-aa3a-89c485de33e0' class='xr-section-summary'  title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
-      ],
-      "text/plain": [
-       "<xarray.Dataset>\n",
-       "Dimensions:        (name: 9, time: 1)\n",
-       "Coordinates:\n",
-       "  * name           (name) <U32 'Sun' 'Mercury' 'Venus' ... 'Uranus' 'Neptune'\n",
-       "  * time           (time) float64 0.0\n",
-       "Data variables: (12/14)\n",
-       "    particle_type  (name) <U32 'Central Body' 'Massive Body' ... 'Massive Body'\n",
-       "    id             (name) int64 0 1 2 3 4 5 6 7 8\n",
-       "    Gmass          (time, name) float64 39.48 6.554e-06 ... 0.001724 0.002034\n",
-       "    radius         (time, name) float64 0.00465 1.631e-05 ... 0.0001646\n",
-       "    j2rp2          (time, name) float64 4.754e-12 nan nan nan ... nan nan nan\n",
-       "    j4rp4          (time, name) float64 -2.247e-18 nan nan nan ... nan nan nan\n",
-       "    ...             ...\n",
-       "    xhz            (time, name) float64 nan -0.002529 ... -0.0407 -0.7287\n",
-       "    vhx            (time, name) float64 nan -9.241 3.454 ... -1.315 -0.08755\n",
-       "    vhy            (time, name) float64 nan 7.755 6.477 ... 1.932 0.5372 1.149\n",
-       "    vhz            (time, name) float64 nan 1.481 -0.1103 ... 0.01905 -0.02168\n",
-       "    ntp            (time) int64 0\n",
-       "    npl            (time) int64 8"
-      ]
-     },
-     "execution_count": 3,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
-    "sim_nogr = swiftest.Simulation(tstop=1000.0, dt=0.005, tstep_out=1.0, integrator=\"whm\", init_cond_format=\"XV\", general_relativity=False,param_file=\"param.nogr.in\", output_file_name=\"bin.nogr.nc\")\n",
+    "sim_nogr = swiftest.Simulation(param_file=\"param.nogr.in\", output_file_name=\"bin.nogr.nc\")\n",
     "sim_nogr.add_solar_system_body([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"])"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
     "%%capture\n",
-    "#Capture the output of the run so it doesn't take up the whole Notebook\n",
-    "sim_gr.run()\n",
-    "sim_nogr.run()"
+    "tstep_out = 10.0\n",
+    "sim_gr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"whm\",general_relativity=True)"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": null,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "%%capture\n",
+    "sim_nogr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"whm\",general_relativity=False)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -934,14 +69,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
     "#Get the ephemerides of Mercury for the same timeframe as the simulation\n",
     "obj = Horizons(id='1', location='@sun',\n",
     "               epochs={'start':start_date, 'stop':stop_date,\n",
-    "                       'step':'1y'})\n",
+    "                       'step':'10y'})\n",
     "el = obj.elements()\n",
     "t = (el['datetime_jd']-el['datetime_jd'][0]) / 365.25\n",
     "varpi_obs = el['w'] + el['Omega']"
@@ -949,18 +84,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 7,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "name": "stderr",
-     "output_type": "stream",
-     "text": [
-      "/home/daminton/.conda/envs/cent7/2020.11-py38/swiftest/lib/python3.8/site-packages/xarray/core/computation.py:733: RuntimeWarning: invalid value encountered in remainder\n",
-      "  result_data = func(*input_data)\n"
-     ]
-    }
-   ],
+   "outputs": [],
    "source": [
     "# Compute the longitude of the periapsis\n",
     "sim_gr.data['varpi'] = np.mod(sim_gr.data['omega'] + sim_gr.data['capom'],360)\n",
@@ -969,7 +95,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 12,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -980,13 +106,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 13,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
-    "dvarpi_gr = np.diff(varpisim_gr)  * 3600 * 100\n",
-    "dvarpi_nogr = np.diff(varpisim_nogr)  * 3600 * 100 \n",
-    "dvarpi_obs = np.diff(varpi_obs) / np.diff(t) * 3600 * 100 "
+    "dvarpi_gr = np.diff(varpisim_gr)  * 3600 * 100 / tstep_out\n",
+    "dvarpi_nogr = np.diff(varpisim_nogr)  * 3600 * 100 / tstep_out\n",
+    "dvarpi_obs = np.diff(varpi_obs) / np.diff(t) * 3600 * 100"
    ]
   },
   {
@@ -997,12 +123,13 @@
    "source": [
     "fig, ax = plt.subplots()\n",
     "\n",
-    "ax.plot(t, varpi_obs, label=\"JPL Horizons\")\n",
-    "ax.plot(tsim, varpisim_gr, label=\"Swiftest WHM GR\")\n",
-    "ax.plot(tsim, varpisim_nogr, label=\"Swiftest WHM No GR\")\n",
+    "ax.plot(t, varpi_obs, label=\"JPL Horizons\",linewidth=2.5)\n",
+    "ax.plot(tsim, varpisim_gr, label=\"Swiftest WHM GR\",linewidth=1.5)\n",
+    "ax.plot(tsim, varpisim_nogr, label=\"Swiftest WHM No GR\",linewidth=1.5)\n",
     "ax.set_xlabel('Time (y)')\n",
     "ax.set_ylabel('Mercury $\\\\varpi$ (deg)')\n",
     "ax.legend()\n",
+    "plt.savefig(\"whm_gr_mercury_precession.png\",dpi=300)\n",
     "print('Mean precession rate for Mercury long. peri. (arcsec/100 y)')\n",
     "print(f'JPL Horizons         : {np.mean(dvarpi_obs)}')\n",
     "print(f'Swiftest No GR       : {np.mean(dvarpi_nogr)}')\n",