From 64c61fa1421100e1852ee078e2b4a222dcda9f66 Mon Sep 17 00:00:00 2001
From: Carlisle Wishard <cwishard@purdue.edu>
Date: Tue, 6 Dec 2022 11:06:33 -0500
Subject: [PATCH] deleted notebook

---
 examples/whm_gr_test/whm_gr_test.ipynb | 158 -------------------------
 1 file changed, 158 deletions(-)
 delete mode 100644 examples/whm_gr_test/whm_gr_test.ipynb

diff --git a/examples/whm_gr_test/whm_gr_test.ipynb b/examples/whm_gr_test/whm_gr_test.ipynb
deleted file mode 100644
index a3dac65b2..000000000
--- a/examples/whm_gr_test/whm_gr_test.ipynb
+++ /dev/null
@@ -1,158 +0,0 @@
-{
- "cells": [
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "import swiftest\n",
-    "from astroquery.jplhorizons import Horizons\n",
-    "import datetime\n",
-    "import numpy as np\n",
-    "import matplotlib.pyplot as plt"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "sim_gr = swiftest.Simulation(simdir=\"gr\")\n",
-    "sim_gr.add_solar_system_body([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"])"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "sim_nogr = swiftest.Simulation(simdir=\"nogr\")\n",
-    "sim_nogr.add_solar_system_body([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"])"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "run_args = {\"tstop\":1000.0, \"dt\":0.005, \"tstep_out\":10.0, \"dump_cadence\": 0,\"integrator\":\"whm\"}"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "sim_gr.run(**run_args,general_relativity=True)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "sim_nogr.run(**run_args,general_relativity=False)"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "# Get the start and end date of the simulation so we can compare with the real solar system\n",
-    "start_date = sim_gr.ephemeris_date\n",
-    "tstop_d = sim_gr.param['TSTOP'] * sim_gr.param['TU2S'] / swiftest.JD2S\n",
-    "\n",
-    "stop_date = (datetime.datetime.fromisoformat(start_date) + datetime.timedelta(days=tstop_d)).isoformat()"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "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':'10y'})\n",
-    "el = obj.elements()\n",
-    "t = (el['datetime_jd']-el['datetime_jd'][0]) / 365.25\n",
-    "varpi_obs = el['w'] + el['Omega']"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "varpisim_gr= sim_gr.data['varpi'].sel(name=\"Mercury\")\n",
-    "varpisim_nogr= sim_nogr.data['varpi'].sel(name=\"Mercury\")\n",
-    "tsim = sim_gr.data['time']"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "dvarpi_gr = np.diff(varpisim_gr)  * 3600 * 100 / run_args['tstep_out']\n",
-    "dvarpi_nogr = np.diff(varpisim_nogr)  * 3600 * 100 / run_args['tstep_out']\n",
-    "dvarpi_obs = np.diff(varpi_obs) / np.diff(t) * 3600 * 100"
-   ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": [
-    "fig, ax = plt.subplots()\n",
-    "\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",
-    "print(f'Swiftest GR          : {np.mean(dvarpi_gr)}')\n",
-    "print(f'Obs - Swiftest GR    : {np.mean(dvarpi_obs - dvarpi_gr)}')\n",
-    "print(f'Obs - Swiftest No GR : {np.mean(dvarpi_obs - dvarpi_nogr)}')"
-   ]
-  }
- ],
- "metadata": {
-  "kernelspec": {
-   "display_name": "Python (My debug_env Kernel)",
-   "language": "python",
-   "name": "debug_env"
-  },
-  "language_info": {
-   "codemirror_mode": {
-    "name": "ipython",
-    "version": 3
-   },
-   "file_extension": ".py",
-   "mimetype": "text/x-python",
-   "name": "python",
-   "nbconvert_exporter": "python",
-   "pygments_lexer": "ipython3",
-   "version": "3.8.5"
-  }
- },
- "nbformat": 4,
- "nbformat_minor": 4
-}