Fixed up helio_gr_test scripts and notebooks

MintonGroup · Dec 4, 2022 · d972a21 · d972a21
1 parent eb34415
commit d972a21
Show file tree

Hide file tree

Showing 3 changed files with 0 additions and 23 deletions.
diff --git a/examples/helio_gr_test/helio_gr_test.ipynb b/examples/helio_gr_test/helio_gr_test.ipynb
@@ -88,17 +88,6 @@
     "varpi_obs = el['w'] + el['Omega']"
    ]
   },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "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",
-    "sim_nogr.data['varpi'] = np.mod(sim_nogr.data['omega'] + sim_nogr.data['capom'],360)"
-   ]
-  },
   {
    "cell_type": "code",
    "execution_count": null,
@@ -143,13 +132,6 @@
     "print(f'Obs - Swiftest GR    : {np.mean(dvarpi_obs - dvarpi_gr)}')\n",
     "print(f'Obs - Swiftest No GR : {np.mean(dvarpi_obs - dvarpi_nogr)}')"
    ]
-  },
-  {
-   "cell_type": "code",
-   "execution_count": null,
-   "metadata": {},
-   "outputs": [],
-   "source": []
   }
  ],
  "metadata": {

diff --git a/examples/helio_gr_test/helio_gr_test.py b/examples/helio_gr_test/helio_gr_test.py
@@ -30,10 +30,6 @@
 t = (el['datetime_jd']-el['datetime_jd'][0]) / 365.25
 varpi_obs = el['w'] + el['Omega']
 
-# Compute the longitude of the periapsis
-sim_gr.data['varpi'] = np.mod(sim_gr.data['omega'] + sim_gr.data['capom'],360)
-sim_nogr.data['varpi'] = np.mod(sim_nogr.data['omega'] + sim_nogr.data['capom'],360)
-
 varpisim_gr= sim_gr.data['varpi'].sel(name="Mercury")
 varpisim_nogr= sim_nogr.data['varpi'].sel(name="Mercury")
 tsim = sim_gr.data['time']

diff --git a/python/swiftest/swiftest/io.py b/python/swiftest/swiftest/io.py
@@ -847,7 +847,6 @@ def swifter2xr(param, verbose=True):
         if verbose: print(f"Successfully converted {ds.sizes['time']} output frames.")
     return ds
 
-
 def swiftest2xr(param, verbose=True):
     """
     Converts a Swiftest binary data file into an xarray DataSet.