Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Updates to helio_gr_test notebook
Browse files Browse the repository at this point in the history
  • Loading branch information
daminton committed Dec 4, 2022
1 parent e2282da commit d7357da
Showing 1 changed file with 17 additions and 9 deletions.
26 changes: 17 additions & 9 deletions examples/helio_gr_test/swiftest_relativity.ipynb
Original file line number Diff line number Diff line change
Expand Up @@ -19,7 +19,7 @@
"metadata": {},
"outputs": [],
"source": [
"sim_gr = swiftest.Simulation(param_file=\"grsim/param.gr.in\", output_file_name=\"bin.gr.nc\")\n",
"sim_gr = swiftest.Simulation(simdir=\"gr\")\n",
"sim_gr.add_solar_system_body([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"])"
]
},
Expand All @@ -29,7 +29,7 @@
"metadata": {},
"outputs": [],
"source": [
"sim_nogr = swiftest.Simulation(param_file=\"nogrsim/param.nogr.in\", output_file_name=\"bin.nogr.nc\")\n",
"sim_nogr = swiftest.Simulation(simdir=\"nogr\")\n",
"sim_nogr.add_solar_system_body([\"Sun\",\"Mercury\",\"Venus\",\"Earth\",\"Mars\",\"Jupiter\",\"Saturn\",\"Uranus\",\"Neptune\"])"
]
},
Expand All @@ -39,8 +39,7 @@
"metadata": {},
"outputs": [],
"source": [
"tstep_out = 10.0\n",
"sim_gr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"helio\",general_relativity=True)"
"run_args = {\"tstop\":1000.0, \"dt\":0.005, \"tstep_out\":10.0, \"dump_cadence\": 0,\"integrator\":\"helio\"}"
]
},
{
Expand All @@ -49,7 +48,16 @@
"metadata": {},
"outputs": [],
"source": [
"sim_nogr.run(tstop=1000.0, dt=0.005, tstep_out=tstep_out, integrator=\"helio\",general_relativity=False)"
"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)"
]
},
{
Expand Down Expand Up @@ -108,8 +116,8 @@
"metadata": {},
"outputs": [],
"source": [
"dvarpi_gr = np.diff(varpisim_gr) * 3600 * 100 / tstep_out\n",
"dvarpi_nogr = np.diff(varpisim_nogr) * 3600 * 100 / tstep_out\n",
"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"
]
},
Expand Down Expand Up @@ -146,9 +154,9 @@
],
"metadata": {
"kernelspec": {
"display_name": "swiftest",
"display_name": "Python (My debug_env Kernel)",
"language": "python",
"name": "swiftest"
"name": "debug_env"
},
"language_info": {
"codemirror_mode": {
Expand Down

0 comments on commit d7357da

Please sign in to comment.