PHOEBE 2.1 Documentation

2.1 Docs

Prev: Wilson-Devinney Style Meshing Next: Binary with Spots

IPython Notebook | Python Script

Detached Binary: Roche vs Rotstar


%matplotlib inline

As always, let’s do imports and initialize a logger and a new bundle. See Building a System for more details.

import phoebe
from phoebe import u # units
import numpy as np
import matplotlib.pyplot as plt

logger = phoebe.logger()

b = phoebe.default_binary()

Adding Datasets

Now we’ll create an empty mesh dataset at quarter-phase so we can compare the difference between using roche and rotstar for deformation potentials:

b.add_dataset('mesh', times=[0.75], dataset='mesh01')
<ParameterSet: 4 parameters | contexts: compute, dataset>

Running Compute

Let’s set the radius of the primary component to be large enough to start to show some distortion when using the roche potentials.

b['requiv@primary@component'] = 1.8

Now we’ll compute synthetics at the times provided using the default options

b.run_compute(irrad_method='none', distortion_method='roche', model='rochemodel')
<ParameterSet: 6 parameters | components: primary, secondary>
b.run_compute(irrad_method='none', distortion_method='rotstar', model='rotstarmodel')
<ParameterSet: 6 parameters | components: primary, secondary>


afig, mplfig = b.plot(model='rochemodel',show=True)
afig, mplfig = b.plot(model='rotstarmodel',show=True)
Prev: Wilson-Devinney Style Meshing Next: Binary with Spots
Last update: 10/29/2018 9:20 a.m. (CET)