https://github.com/austince/planet-simulations

Planet Simulations for Introduction to Astrophysics PEP 336
https://github.com/austince/planet-simulations

astrophysics data-processing data-science python simulation

Last synced: over 1 year ago
JSON representation

Planet Simulations for Introduction to Astrophysics PEP 336

• Host: GitHub
• URL: https://github.com/austince/planet-simulations
• Owner: austince
• Created: 2018-02-07T22:07:42.000Z (over 8 years ago)
• Default Branch: master
• Last Pushed: 2018-02-26T01:40:10.000Z (over 8 years ago)
• Last Synced: 2025-03-08T18:06:49.232Z (over 1 year ago)
• Topics: astrophysics, data-processing, data-science, python, simulation
• Language: Python
• Homepage:
• Size: 16.8 MB
• Stars: 1
• Watchers: 2
• Forks: 0
• Open Issues: 0
• Metadata Files:
  • Readme: README.md

Awesome Lists containing this project

README

          
# Earth Orbit Simulations for PEP 336

## Requirements

**Only runs with Python 3**  

3rd party packages can be install either with [Anaconda](https://www.anaconda.com/) 

by `conda env create -f env.yml` or with [pip](https://packaging.python.org/tutorials/installing-packages/)

by `pip install -r requirements.txt`.

## Running

The `cli.py` script runs the simulations and writes the 

data to a csv file. All the running options are explained

by running `python cli.py --help` but here are a few good ones

to remember:

- `--plot` will plot the results in a polar projection.

- `--output ` will save the results to `filename`. (Required!)

- `--method {euler,rk2,rk4}` will run one of those three simulations.

- `-dt ` will set the time delta to `days`.

Here's a full example that runs a *RK4* simulation for *5 years* with a 

time delta of *0.5 days* and plots it. 

`python cli.py -dt 0.5 --end-time 1825 --method rk4 --output rk4_dt=0.5_end=1825.csv --plot`

## Examples

### Euler  

![](./examples/euler_dt=0.1_end=1825.csvpolyplot.png)  

![](./examples/euler_dt=0.5_end=1825.csvpolyplot.png)  

![](./examples/euler_dt=1_end=1825.csvpolyplot.png)  

![](./examples/euler_dt=2.5_end=1825.csvpolyplot.png)

![](./examples/euler_dt=5_end=1825.csvpolyplot.png)

### RK2  

![](./examples/rk2_dt=0.1_end=1825.csvpolyplot.png)  

![](./examples/rk2_dt=0.5_end=1825.csvpolyplot.png)

![](./examples/rk2_dt=1_end=1825.csvpolyplot.png)  

![](./examples/rk2_dt=2.5_end=1825.csvpolyplot.png)

![](./examples/rk2_dt=5_end=1825.csvpolyplot.png)

### RK4  

![](./examples/rk4_dt=0.1_end=1825.csvpolyplot.png)  

![](./examples/rk4_dt=0.5_end=1825.csvpolyplot.png)

![](./examples/rk4_dt=1_end=1825.csvpolyplot.png)  

![](./examples/rk4_dt=2.5_end=1825.csvpolyplot.png)

![](./examples/rk4_dt=5_end=1825.csvpolyplot.png)