Re: [PHOEBE-discuss] Example of Simple Binary Script

[Sara M. C. <sar...@gm...>]

# Phoebe 2 Alpha #

import numpy as np
import matplotlib.pyplot as plt
import phoebe; #phoebe.download_atm()


eb = phoebe.Bundle()

# Orbit parameters
eb['period@orbit'] = (1.44631, 'd')
eb['incl@orbit'] = (54.5, 'deg')
eb['sma@orbit'] = (6.54, 'Rsol')
eb['q@orbit'] = (0.125)
eb['ecc@orbit'] = (0)

# Primary component
eb['teff@primary'] = (7150, 'K')
eb['abun@primary'] = (0.30)
eb['syncpar@primary'] = (1.0)
eb['pot@primary'] = (3.2)

# Secondary component
eb['teff@secondary'] = (6500, 'K')
eb['abun@secondary'] = (0.30)
eb['syncpar@secondary'] = (1.0)
eb['pot@secondary'] = (2.14)

# Synthetic light curve
lc01 = eb.lc_fromarrays(phase = np.linspace(-0.6, 0.6, 201))

# Observations
lcData = np.loadtxt('ModTime_1.txt').T
eb.lc_fromarrays(time=lcData[0],flux=lcData[1], # errors give the scatter from the bin in the mean light curve
                dataref = 'keplerlc', # this is the name to which we have to refer when saying twigs, e.g.'scale@keplerlc'
                passband='JOHNSON.V') # this sets the pass band
#eb['lc01@lcobs']

# Compute and plot
eb.run_compute('detailed')
#eb.plot_syn('lc01', fmt='r-')
#eb.plot_obs('keplerlc', fmt='bo')
plt.plot(eb['time@keplerlc@lcsyn'],eb['flux@keplerlc@lcsyn'],'r-', ms=10, mew=2)
plt.plot(eb['time@keplerlc@lcobs'],eb['flux@keplerlc@lcobs'],'bo', ms=10, mew=2)
plt.show()