PHOEBE 2.1 Documentation

2.1 Docs



Prev: Single Star with Spots Next: Comparing PHOEBE 2 vs PHOEBE Legacy
.

IPython Notebook | Python Script

Comparing Spots in PHOEBE 2 vs PHOEBE Legacy

Setup

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

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

logger = phoebe.logger()

b = phoebe.default_binary()

Adding Spots and Compute Options

b.add_spot(component='primary', relteff=0.8, radius=20, colat=45, colon=90, feature='spot01')
<ParameterSet: 4 parameters | qualifiers: colat, radius, long, relteff>
b.add_dataset('lc', times=np.linspace(0,1,101))
<ParameterSet: 15 parameters | contexts: compute, dataset>
b.add_compute('phoebe', irrad_method='none', compute='phoebe2')
<ParameterSet: 17 parameters | components: primary, secondary>
b.add_compute('legacy', irrad_method='none', compute='phoebe1')
<ParameterSet: 7 parameters | components: primary, secondary>

Let’s use the external atmospheres available for both phoebe1 and phoebe2

b.set_value_all('atm', 'extern_planckint')
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
b.set_value_all('ld_func', 'logarithmic')
b.set_value_all('ld_coeffs', [0.0, 0.0])
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@primary or ld_func@primary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
Mon, 29 Oct 2018 10:09 PARAMETERS   WARNING ld_func='interp' only supported by atm='ck2004'.  Either change atm@secondary or ld_func@secondary@lc01  If not addressed, this warning will continue to be raised and will throw an error at run_compute.
b.run_compute('phoebe2', model='phoebe2model')
<ParameterSet: 2 parameters | qualifiers: fluxes, times>
b.run_compute('phoebe1', model='phoebe1model')
Mon, 29 Oct 2018 10:09 IO           WARNING requiv_max@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING freq@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING pitch@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING yaw@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING long_an@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING irrad_frac_lost_bol@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING ld_func_bol@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING mass@primary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING requiv_max@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING freq@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING pitch@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING yaw@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING long_an@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING irrad_frac_lost_bol@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING ld_func_bol@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING mass@secondary@star@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING pblum_ref@secondary@lc01@lc_dep@dataset has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING fti_method@lc01@phoebe01@phoebe@compute has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING Finite integration Time is not fully supported and will be turned off by legacy wrapper before computation
Mon, 29 Oct 2018 10:09 IO           WARNING lc_method@lc01@phoebe2@phoebe@compute has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING lc_method@lc01@phoebe01@phoebe@compute has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING intens_weighting@lc01@lc_dep@dataset has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING fti_method@lc01@phoebe2@phoebe@compute has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING pblum_ref@primary@lc01@lc_dep@dataset has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING freq@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING t0_perpass@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING t0_supconj@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING mean_anom@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING long_an@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING asini@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING ecosw@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING esinw@binary@orbit@component has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING irrad_method@phoebe1@legacy@compute has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING t0@system has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING ra@system has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING dec@system has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING epoch@system has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING distance@system has no phoebe 1 corollary
Mon, 29 Oct 2018 10:09 IO           WARNING hierarchy@system has no phoebe 1 corollary
<ParameterSet: 2 parameters | qualifiers: fluxes, times>

Plotting

afig, mplfig = b.plot(legend=True, ylim=(1.95, 2.05), show=True)
../../_images/legacy_spots_16_0.png
Prev: Single Star with Spots Next: Comparing PHOEBE 2 vs PHOEBE Legacy
.
Last update: 10/29/2018 9:20 a.m. (CET)