PHOEBE 2.1 Documentation

2.1 Docs

Source code for phoebe.distortions.rotstar

import numpy as np

from phoebe import c, u
import libphoebe

[docs]def rotfreq_to_omega(rotfreq,, solar_units=False): """ TODO: add documentation NOTE: everything MUST be in consistent units according to solar_units bool """ if solar_units: omega = rotfreq / (2*np.pi) / np.sqrt(**3/u.d**2).value/scale**3) else: # then SI units omega = rotfreq / (2*np.pi) / np.sqrt(c.GM_sun.value/scale**3) # print "*** rotstar.rotfreq_to_omega", rotfreq, scale, solar_units, omega return omega
[docs]def rpole2potential(rpole, rotfreq, solar_units=False): """ Transforms polar radius to surface potential """ if not solar_units: rpole = rpole/ rpole_ = np.array([0., 0., rpole]) omega = rotfreq_to_omega(rotfreq, solar_units=solar_units) pot = libphoebe.rotstar_Omega(omega, rpole_) # print "*** rotstar.rpole2potential", rpole, rotfreq, solar_units, omega, pot return pot
[docs]def potential2rpole(pot, rotfreq, solar_units=False): """ Transforms surface potential to polar radius """ omega = rotfreq_to_omega(rotfreq, scale=1.0, solar_units=solar_units) # print "*** rotstar.potential2rpole", pot, rotfreq, solar_units, omega rpole = libphoebe.rotstar_pole(omega, pot) if solar_units: return rpole else: return rpole*
Last update: 10/29/2018 9:20 a.m. (CET)