Oblate planet example 1¶

In [1]:

%pylab inline

%pylab is deprecated, use %matplotlib inline and import the required libraries.
Populating the interactive namespace from numpy and matplotlib

In [10]:

from pytransit import OPModel, RRModel

In [3]:

time = linspace(-0.2, 0.2, 300)

In [13]:

om = OPModel('quadratic')
tm = RRModel('quadratic')

In [14]:

om.set_data(time)
tm.set_data(time)

In [22]:

k, f, alpha = 0.1, 0.1, 0.5*pi
ldc = array([0.2, 0.3])

fig1, axs1 = subplots(3, 3, figsize=(13,5), sharex='all', sharey='all')
fig2, axs2 = subplots(3, 3, figsize=(13,5), sharex='all', sharey='all')

for i, f in enumerate((0.05, 0.3, 0.5)):
    for j, alpha in enumerate(linspace(0, 0.5*pi, 3)):
        flux_opm = om.evaluate(k, f, alpha, ldc, t0=0.0, p=4.0, a=5.0, i=0.45*pi)
        flux = tm.evaluate(k, ldc, t0=0.0, p=4.0, a=5.0, i=0.45*pi)
        axs1[i, j].plot(time, flux, c='0.5', ls='--')
        axs1[i, j].plot(time, flux_opm, 'k-')

        axs2[i,j].plot(time, flux_opm-flux, 'k-')

setp(axs1[:,0], ylabel='Flux')
setp(axs1[-1,:], xlabel='Time [d]')
fig1.tight_layout()

setp(axs2[:,0], ylabel='Oblate - spherical')
setp(axs2[-1,:], xlabel='Time [d]')
fig2.tight_layout()

In [ ]: