Ejemplo n.º 1
0
split_path = os.path.dirname(__file__).split('/')[:-2]
path = ''
for p in split_path:
    path += '/' + p
path += '/'
sys.path.append(path)

import numpy as np
import pandas as pd
import glob
import matplotlib.pyplot as plt
from solar_system import solar_System
from exoplanet_simulation import simulate

try:
    star_sys = solar_System('star.csv', 'planets.csv')
except:
    star_sys = solar_System('Exoplanets_data/GJ_832/star.csv',
                            'Exoplanets_data/GJ_832/planets.csv')

times = np.linspace(0, 1 * 10**(7), 12345) + 0j
# times = np.linspace(-0, .1, 500)+0j
# times = np.linspace(10**6, 10**10, 10000)+0j
# times = np.logspace(6, 10, 10000)+0j
eccs = simulate(star_sys,
                times,
                plot=False,
                plot_orbit=False,
                save=False,
                folder_name='Exoplanets_data/GJ_832')
Ejemplo n.º 2
0
import sys
sys.path.append('/Users/Sahl/Desktop/University/Year_5/Mphys_project/')

import numpy as np
import pandas as pd
import glob
import matplotlib.pyplot as plt
from solar_system import solar_System
from exoplanet_simulation import simulate

try:
    star_sys = solar_System('star.csv', 'planets.csv')
except:
    star_sys = solar_System('Exoplanets_data/HD_12661/star.csv', 'Exoplanets_data/HD_12661/planets.csv')

times = np.linspace(0, 1*10**(5), 1234)+0j
    # times = np.linspace(-0, .1, 500)+0j
    # times = np.linspace(10**6, 10**10, 10000)+0j
    # times = np.logspace(6, 10, 10000)+0j
eccs = simulate(star_sys, times, plot=False, plot_orbit=False, save=False, folder_name='Exoplanets_data/HD_12661')

times = np.real(times)
try:
    df = pd.read_csv('Exoplanets_data/HD_12661'+'/b_nbody.csv')
except:
    df = pd.read_csv('b_nbody.csv')

f, (ax1, ax2) = plt.subplots(2, sharex=True, figsize=(8,7))
ax0 = f.add_subplot(111, frame_on=False)   # creating a single axes
ax0.set_xticks([])
ax0.set_yticks([])
Ejemplo n.º 3
0
split_path = os.path.dirname(__file__).split('/')[:-2]
path = ''
for p in split_path:
    path += '/' + p
path += '/'
sys.path.append(path)

import numpy as np
import pandas as pd
import glob
import matplotlib.pyplot as plt
from solar_system import solar_System
from exoplanet_simulation import simulate

try:
    star_sys = solar_System('star.csv', 'planets.csv')
except:
    star_sys = solar_System('Exoplanets_data/HD_113538/star.csv',
                            'Exoplanets_data/HD_113538/planets.csv')

times = np.linspace(0, 4 * 10**(4), 1234) + 0j
# times = np.linspace(-0, .1, 500)+0j
# times = np.linspace(10**6, 10**10, 10000)+0j
# times = np.logspace(6, 10, 10000)+0j
eccs = simulate(star_sys,
                times,
                plot=False,
                plot_orbit=False,
                save=False,
                folder_name='Exoplanets_data/HD_113538')
Ejemplo n.º 4
0
import sys
sys.path.append('/Users/Sahl/Desktop/University/Year_5/Mphys_project/')

import numpy as np
import pandas as pd
import glob
import matplotlib.pyplot as plt
from solar_system import solar_System
from exoplanet_simulation import simulate

try:
    star_sys = solar_System('star.csv', 'planets.csv')
except:
    star_sys = solar_System('Exoplanets_data/61_Vir/star.csv',
                            'Exoplanets_data/61_Vir/planets.csv')

times = np.linspace(0, 1 * 10**(5), 1234) + 0j
# times = np.linspace(-0, .1, 500)+0j
# times = np.linspace(10**6, 10**10, 10000)+0j
# times = np.logspace(6, 10, 10000)+0j
eccs = simulate(star_sys,
                times,
                plot=False,
                plot_orbit=False,
                save=False,
                folder_name='Exoplanets_data/61_Vir')

times = np.real(times)
try:
    df = pd.read_csv('Exoplanets_data/61_Vir' + '/b_nbody.csv')
except:
Ejemplo n.º 5
0
split_path = os.path.dirname(__file__).split('/')[:-2]
path = ''
for p in split_path:
    path += '/' + p
path += '/'
sys.path.append(path)

import numpy as np
import pandas as pd
import glob
import matplotlib.pyplot as plt
from solar_system import solar_System
from exoplanet_simulation import simulate

try:
    star_sys = solar_System('star.csv', 'planets.csv')
except:
    star_sys = solar_System('Exoplanets_data/HD_108874/star.csv',
                            'Exoplanets_data/HD_108874/planets.csv')

times = np.linspace(0, 1 * 10**(5), 1234) + 0j
# times = np.linspace(-0, .1, 500)+0j
# times = np.linspace(10**6, 10**10, 10000)+0j
# times = np.logspace(6, 10, 10000)+0j
eccs = simulate(star_sys,
                times,
                plot=False,
                plot_orbit=False,
                save=False,
                folder_name='Exoplanets_data/HD_108874')
Ejemplo n.º 6
0
split_path = os.path.dirname(__file__).split('/')[:-2]
path = ''
for p in split_path:
    path += '/'+p
path += '/'
sys.path.append(path)

import numpy as np
import pandas as pd
import glob
import matplotlib.pyplot as plt
from solar_system import solar_System
from exoplanet_simulation import simulate

try:
    star_sys = solar_System('star.csv', 'planets.csv')
except:
    star_sys = solar_System('Exoplanets_data/47_UMa/star.csv', 'Exoplanets_data/47_UMa/planets.csv')

times = np.linspace(0, 7*10**(4), 1234)+0j
    # times = np.linspace(-0, .1, 500)+0j
    # times = np.linspace(10**6, 10**10, 10000)+0j
    # times = np.logspace(6, 10, 10000)+0j
eccs = simulate(star_sys, times, plot=False, plot_orbit=False, save=False, folder_name='Exoplanets_data/47_UMa')

times = np.real(times)
try:
    df = pd.read_csv('Exoplanets_data/47_UMa'+'/b_nbody.csv')
except:
    df = pd.read_csv('b_nbody.csv')