PHOEBE 2.1 Documentation

2.1 Docs



IPython Notebook | Python Script

Single Star with Spots

Setup

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_star()

Adding Spots

Let’s add one spot to our star. Since there is only one star, the spot will automatically attach without needing to provide component (as is needed in the binary with spots example

b.add_spot(radius=30, colat=80, long=0, relteff=0.9)
<ParameterSet: 4 parameters | qualifiers: colat, radius, long, relteff>

Spot Parameters

A spot is defined by the colatitude and longitude of its center, its angular radius, and the ratio of temperature of the spot to the local intrinsic value.

print b['spot']
ParameterSet: 4 parameters
            colat@spot01@feature: 80.0 deg
             long@spot01@feature: 0.0 deg
           radius@spot01@feature: 30.0 deg
          relteff@spot01@feature: 0.9

The ‘colat’ parameter defines the latitude on the star measured from its North Pole. The ‘long’ parameter measures the longitude of the spot - with longitude = 0 being defined as pointing towards the observer at t0.

times = np.linspace(0, 10, 11)
b.set_value('period', 10)
b.add_dataset('mesh', times=times, columns=['teffs'])
<ParameterSet: 4 parameters | contexts: compute, dataset>
b.run_compute(distortion_method='rotstar', irrad_method='none')
<ParameterSet: 34 parameters | qualifiers: xyz_elements, uvw_elements, teffs, times>
afig, mplfig = b.plot(x='us', y='vs', fc='teffs',
                      animate=True, save='single_spots_1.gif', save_kwargs={'writer': 'imagemagick'})
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
../../_images/single_spots_13_1.png
animation

animation

If we set t0 to 5 instead of zero, then the spot will cross the line-of-sight at t=5 (since the spot’s longitude is 0).

b.set_value('t0', 5)
b.run_compute(distortion_method='rotstar', irrad_method='none')
Mon, 15 Oct 2018 17:38 BUNDLE       WARNING overwriting model: latest
<ParameterSet: 34 parameters | qualifiers: xyz_elements, uvw_elements, teffs, times>
afig, mplfig = b.plot(x='us', y='vs', fc='teffs',
                     animate=True, save='single_spots_2.gif', save_kwargs={'writer': 'imagemagick'})
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
../../_images/single_spots_18_1.png
animation

animation

And if we change the inclination to 0, we’ll be looking at the north pole of the star. This clearly illustrates the right-handed rotation of the star. At time=t0=5 the spot will now be pointing in the negative y-direction.

b.set_value('incl', 0)
b.run_compute(distortion_method='rotstar', irrad_method='none')
Mon, 15 Oct 2018 17:38 BUNDLE       WARNING overwriting model: latest
<ParameterSet: 34 parameters | qualifiers: xyz_elements, uvw_elements, teffs, times>
afig, mplfig = b.plot(x='us', y='vs', fc='teffs',
                     animate=True, save='single_spots_3.gif', save_kwargs={'writer': 'imagemagick'})
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
WARNING: pad_aspect not supported for animations, ignoring
../../_images/single_spots_23_1.png
animation

animation