PHOEBE 2.0 Documentation

2.0 Docs

  • 1.0
  • 2.0a
  • ver: 2.0


IPython Notebook | Python Script

Rossiter-McLaughlin Effect

Setup

%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()
WARNING: Constant u'Gravitational constant' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING: Constant u'Solar mass' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING: Constant u'Solar radius' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING: Constant u'Solar luminosity' is already has a definition in the u'si' system [astropy.constants.constant]
WARNING:astropy:Constant u'Solar luminosity' is already has a definition in the u'si' system
/usr/local/lib/python2.7/dist-packages/astropy/units/quantity.py:732: FutureWarning: comparison to None will result in an elementwise object comparison in the future.
  return super(Quantity, self).__eq__(other)

Let’s make a significant mass ratio and radius ratio...

b['q'] = 0.7
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'mass@primary@component'=0.449019459045 solMass from 'mass@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'sma@primary@component'=3.29411764706 solRad from 'sma@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@primary@component'=8.6945945137 from 'pot@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'mass@secondary@component'=0.314313621332 solMass from 'mass@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'sma@secondary@component'=4.70588235294 solRad from 'sma@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@secondary@component'=6.74227787671 from 'pot@secondary@constraint' constraint
b['rpole@primary'] = 1.0
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@primary@component'=8.6945945137 from 'pot@primary@constraint' constraint
b['rpole@secondary'] = 0.5
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@secondary@component'=12.3480525785 from 'pot@secondary@constraint' constraint
b['teff@secondary@component'] = 5000
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00

Make sure the primary star is spinning quickly...

b['syncpar@primary@component'] = 2
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'freq@primary@component'=4.18879 rad / d from 'freq@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'period@primary@component'=1.5 d from 'period@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@primary@component'=8.6945945137 from 'pot@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00

Adding Datasets

Now we’ll add radial velocity and mesh datasets. We’ll add two identical datasets for RVs so that we can have one computed dynamically and the other computed numerically (these options will need to be set later).

b.add_dataset('rv', times=np.linspace(0,2,201), dataset='dynamicalrvs')
# TODO: can't set rv_method here because compute options don't exist yet... and that's kind of annoying
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'times' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'times' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'rvs' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'rvs' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'sigmas' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'sigmas' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_func' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_func' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_coeffs' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_coeffs' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating mass constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'mass@primary@component'=0.449019459045 solMass from 'mass@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating comp_sma constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'sma@primary@component'=3.29411764706 solRad from 'sma@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating rotation_period constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'freq@primary@component'=4.18879 rad / d from 'freq@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'period@primary@component'=1.5 d from 'period@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating incl_aligned constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'incl@primary@component'=90.0 deg from 'incl@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating potential constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@primary@component'=8.6945945137 from 'pot@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating mass constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'mass@secondary@component'=0.314313621332 solMass from 'mass@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating comp_sma constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'sma@secondary@component'=4.70588235294 solRad from 'sma@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating rotation_period constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'freq@secondary@component'=2.094395 rad / d from 'freq@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'period@secondary@component'=3.0 d from 'period@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating incl_aligned constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'incl@secondary@component'=90.0 deg from 'incl@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating potential constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@secondary@component'=12.3480525785 from 'pot@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
<ParameterSet: 15 parameters | contexts: compute, dataset>
b.add_dataset('rv', times=np.linspace(0,2,201), dataset='numericalrvs')
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'times' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'times' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'rvs' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'rvs' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'sigmas' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'sigmas' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_func' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_func' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_coeffs' parameter for {'component': 'primary'}
Wed, 21 Sep 2016 18:29 PARAMETERS   INFO    copying 'ld_coeffs' parameter for {'component': 'secondary'}
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating mass constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'mass@primary@component'=0.449019459045 solMass from 'mass@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating comp_sma constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'sma@primary@component'=3.29411764706 solRad from 'sma@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating rotation_period constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'freq@primary@component'=4.18879 rad / d from 'freq@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'period@primary@component'=1.5 d from 'period@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating incl_aligned constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'incl@primary@component'=90.0 deg from 'incl@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating potential constraint for primary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@primary@component'=8.6945945137 from 'pot@primary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating mass constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'mass@secondary@component'=0.314313621332 solMass from 'mass@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating comp_sma constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'sma@secondary@component'=4.70588235294 solRad from 'sma@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating rotation_period constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'freq@secondary@component'=2.094395 rad / d from 'freq@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'period@secondary@component'=3.0 d from 'period@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating incl_aligned constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'incl@secondary@component'=90.0 deg from 'incl@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    re-creating potential constraint for secondary
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 BUNDLE       INFO    setting 'pot@secondary@component'=12.3480525785 from 'pot@secondary@constraint' constraint
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:29 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
<ParameterSet: 15 parameters | contexts: compute, dataset>

Storing the mesh at every timestep is overkill, and will be both computationally and memory intensive. So let’s just sample at the times we care about.

times = b.get_value('times@primary@numericalrvs@dataset')
times = times[times<0.1]
print times
[ 0.    0.01  0.02  0.03  0.04  0.05  0.06  0.07  0.08  0.09]
b.add_dataset('mesh', dataset='mesh01', times=times)
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating mass constraint for primary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'mass@primary@component'=0.449019459045 solMass from 'mass@primary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating comp_sma constraint for primary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'sma@primary@component'=3.29411764706 solRad from 'sma@primary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating rotation_period constraint for primary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'freq@primary@component'=4.18879 rad / d from 'freq@primary@constraint' constraint
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'period@primary@component'=1.5 d from 'period@primary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating incl_aligned constraint for primary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'incl@primary@component'=90.0 deg from 'incl@primary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating potential constraint for primary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'pot@primary@component'=8.6945945137 from 'pot@primary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating mass constraint for secondary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'mass@secondary@component'=0.314313621332 solMass from 'mass@secondary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating comp_sma constraint for secondary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'sma@secondary@component'=4.70588235294 solRad from 'sma@secondary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating rotation_period constraint for secondary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'freq@secondary@component'=2.094395 rad / d from 'freq@secondary@constraint' constraint
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'period@secondary@component'=3.0 d from 'period@secondary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating incl_aligned constraint for secondary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'incl@secondary@component'=90.0 deg from 'incl@secondary@constraint' constraint
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    re-creating potential constraint for secondary
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    setting 'pot@secondary@component'=12.3480525785 from 'pot@secondary@constraint' constraint
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
<ParameterSet: 2 parameters | contexts: compute, dataset>

Running Compute

Now let’s set the rv_method so that one dataset uses the dynamical method and the other uses the flux-weighted (numerical) method. Note that here we have to use set_value_all or loop over the components, as technically there are parameters for each component-dataset pair.

b.set_value_all('rv_method@dynamicalrvs@compute', 'dynamical')
b.set_value_all('rv_method@numericalrvs@compute', 'flux-weighted')
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 PARAMETERS   WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00

Let’s check to make sure that rv_method is set as we’d expect.

print b['rv_method']
ParameterSet: 4 parameters
  rv_method@primary@dynamical...: dynamical
  rv_method@secondary@dynamic...: dynamical
  rv_method@primary@numerical...: flux-weighted
  rv_method@secondary@numeric...: flux-weighted
b.run_compute(irrad_method='none')
Wed, 21 Sep 2016 18:30 BUNDLE       WARNING Object probably has a convective atm (teff=5000K<6600K), for which gravb_bol=0.32 might be a better approx than gravb_bol=1.00
Wed, 21 Sep 2016 18:30 BUNDLE       INFO    running phoebe backend to create 'latest' model
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    initializing mesh for primary
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    covered surface with 3348 triangles
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    initializing mesh for secondary
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    covered surface with 3310 triangles
Wed, 21 Sep 2016 18:30 PASSBANDS    INFO    loading Bolometric:1760-40000 passband
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 PASSBANDS    INFO    loading Johnson:V passband
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:30 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:31 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
Wed, 21 Sep 2016 18:32 UNIVERSE     INFO    handling reflection (convex case), method='wilson'
<ParameterSet: 750 parameters | kinds: rv, mesh>

Plotting

Now let’s plot the radial velocities.

First we’ll plot the dynamical RVs. Note that dynamical RVs show the true radial velocity of the center of mass of each star, and so we do not see the Rossiter McLaughlin effect.

axs, artists = b['dynamicalrvs@model'].plot(component='primary', color='b')
axs, artists = b['dynamicalrvs@model'].plot(component='secondary', color='r')
Wed, 21 Sep 2016 18:32 PARAMETERS   INFO    calling 'mpl' plotting backend for s@primary@dynamicalrvs@latest@model
Wed, 21 Sep 2016 18:32 PARAMETERS   INFO    calling 'mpl' plotting backend for s@secondary@dynamicalrvs@latest@model
../../_images/rossiter_mclaughlin_25_1.png

But the numerical method integrates over the visible surface elements, giving us what we’d observe if deriving RVs from observed spectra of the binary. Here we do see the Rossiter McLaughlin effect. You’ll also notice that RVs are not available for the secondary star when its completely occulted (they’re nans in the array).

axs, artists = b['numericalrvs@model'].plot(component='primary', color='b')
axs, artists = b['numericalrvs@model'].plot(component='secondary', color='r')
Wed, 21 Sep 2016 18:32 PARAMETERS   INFO    calling 'mpl' plotting backend for s@primary@numericalrvs@latest@model
Wed, 21 Sep 2016 18:32 PARAMETERS   INFO    calling 'mpl' plotting backend for s@secondary@numericalrvs@latest@model
../../_images/rossiter_mclaughlin_27_1.png

To visualize what is happening, we can plot the radial velocities of each surface element in the mesh at one of these times.

Here just plot on the mesh@model parameterset - the mesh will automatically get coordinates from mesh01 and then we point to rvs@numericalrvs for the facecolors.

fig = plt.figure(figsize=(12,12))
axs, artists = b['mesh@model'].plot(time=0.03, facecolor='rvs@numericalrvs', edgecolor=None)
/home/kyle/.local/lib/python2.7/site-packages/phoebe/frontend/plotting.py:257: FutureWarning: elementwise comparison failed; returning scalar instead, but in the future will perform elementwise comparison
  if pckwargs['edgecolors'] in ['none', 'None', None] and pckwargs['facecolors'] not in ['none', 'None', None]:
/home/kyle/.local/lib/python2.7/site-packages/phoebe/frontend/plotting.py:257: FutureWarning: comparison to None will result in an elementwise object comparison in the future.
  if pckwargs['edgecolors'] in ['none', 'None', None] and pckwargs['facecolors'] not in ['none', 'None', None]:
../../_images/rossiter_mclaughlin_29_1.png

Here you can see that the secondary star is blocking part of the “red” RVs of the primary star.

This is essentially the same as plotting the negative z-component of the velocities (for convention - our system is in a right handed system with +z towards the viewer, but RV convention has negative RVs for blue shifts).

We could also plot the RV per triangle by plotting ‘vzs’. Note that this is actually defaulting to an inverted colormap to show you the same colorscheme (‘RdBu_r’ vs ‘RdBu’).

fig = plt.figure(figsize=(12,12))
axs, artists = b['mesh01@model'].plot(time=0.09, facecolor='vzs', edgecolor=None)
Wed, 21 Sep 2016 18:32 PARAMETERS   INFO    calling 'mpl' plotting backend for None@mesh01
Wed, 21 Sep 2016 18:32 PLOTTING     INFO    defaulting to 'RdBu' colormap for 'vzs'
Wed, 21 Sep 2016 18:32 PLOTTING     INFO    forcing color range to be centered around 0 for 'vzs'
../../_images/rossiter_mclaughlin_31_1.png
Prev: Complete Binary Animation Next: Wilson-Devinney Style Meshing
.
Last update: 06/07/2017 11:30 a.m. (CET)