« Phoebe V2のインストール | トップページ | 趣味 »

2021年7月 2日 (金)

Phoebe V2のインストール

http://binary.cocolog-nifty.com/blog/2021/06/post-ad8e35.html
これの続きです。

>パラメータの設定が難しいです。

パラメータ phoebe.default_binary() を試してみます。

import phoebe
from phoebe import u # units
import numpy as np
import matplotlib.pyplot as plt
logger = phoebe.logger()
b = phoebe.default_binary()
b.add_dataset('lc', times=np.linspace(0,1,201), dataset='mylc')
b.run_compute()
afig, mplfig = b['mylc@model'].plot(show=True)

これで T1=T2=6000K の光度曲線を計算して表示します。

202107021 

これを T1=6000K T2=5000K にしてみます。

b.set_value('teff@secondary@component', 5000)
b.run_compute()
afig, mplfig = b['mylc@model'].plot(show=True)

202107022 

今のパラメータを見るには

print(b)

で、表示されます。

ダウンロード - 6000_5000_b.txt

上のファイルをダウンロードするか、以下を眺めるか

>>> print(b)
PHOEBE Bundle: 184 parameters
t0@system: 0.0 d
ra@system: 0.0 deg
dec@system: 0.0 deg
distance@system: 1.0 m
vgamma@system: 0.0 km / s
C ebv@system: 0.0
Av@system: 0.0
Rv@system: 3.1
hierarchy@system: orbit:binary(star:primary, star:secondary)
requiv@primary@component: 1.0 solRad
C requiv_max@primary@component: 2.013275176537638 solRad
teff@primary@component: 6000.0 K
abun@primary@component: 0.0
C logg@primary@component: 4.437551877570185
syncpar@primary@component: 1.0
C period@primary@component: 1.0 d
C freq@primary@component: 6.283185 rad / d
pitch@primary@component: 0.0 deg
yaw@primary@component: 0.0 deg
C incl@primary@component: 90.0 deg
C long_an@primary@component: 0.0 deg
gravb_bol@primary@component: 0.32
irrad_frac_refl_bol@primary...: 0.6
C irrad_frac_lost_bol@primary...: 0.4
ld_mode_bol@primary@component: lookup
ld_func_bol@primary@component: logarithmic
ld_coeffs_source_bol@primar...: auto
C mass@primary@component: 0.9988131358058301 solMass
requiv@secondary@component: 1.0 solRad
C requiv_max@secondary@component: 2.013275176537638 solRad
teff@secondary@component: 5000.0 K
abun@secondary@component: 0.0
C logg@secondary@component: 4.437551877570185
syncpar@secondary@component: 1.0
C period@secondary@component: 1.0 d
C freq@secondary@component: 6.283185 rad / d
pitch@secondary@component: 0.0 deg
yaw@secondary@component: 0.0 deg
C incl@secondary@component: 90.0 deg
C long_an@secondary@component: 0.0 deg
gravb_bol@secondary@component: 0.32
irrad_frac_refl_bol@seconda...: 0.6
C irrad_frac_lost_bol@seconda...: 0.4
ld_mode_bol@secondary@compo...: lookup
ld_func_bol@secondary@compo...: logarithmic
ld_coeffs_source_bol@second...: auto
C mass@secondary@component: 0.9988131358058301 solMass
period@binary@component: 1.0 d
C freq@binary@component: 6.283185 rad / d
dpdt@binary@component: 0.0 s / yr
per0@binary@component: 0.0 deg
dperdt@binary@component: 0.0 deg / yr
ecc@binary@component: 0.0
C t0_perpass@binary@component: -0.25 d
t0_supconj@binary@component: 0.0 d
C t0_ref@binary@component: 0.0 d
C mean_anom@binary@component: 89.99999559997653 deg
incl@binary@component: 90.0 deg
q@binary@component: 1.0
sma@binary@component: 5.3 solRad
long_an@binary@component: 0.0 deg
C asini@binary@component: 5.3 solRad
C ecosw@binary@component: 0.0
C esinw@binary@component: 0.0
C sma@primary@component: 2.65 solRad
C asini@primary@component: 2.65 solRad
C sma@secondary@component: 2.65 solRad
C asini@secondary@component: 2.65 solRad
ebv@constraint: {Av@system} / {Rv@system}
freq@primary@constraint: 6.283185 / {period@primary@component}
logg@primary@constraint: log10((({mass@primary@component} / ({requiv@primary@component} ** 2.000000)) * 2942.206218) * 9.319541)
irrad_frac_lost_bol@primary...: 1.000000 - {irrad_frac_refl_bol@primary@component}
freq@secondary@constraint: 6.283185 / {period@secondary@component}
logg@secondary@constraint: log10((({mass@secondary@component} / ({requiv@secondary@component} ** 2.000000)) * 2942.206218) * 9.319541)
irrad_frac_lost_bol@seconda...: 1.000000 - {irrad_frac_refl_bol@secondary@component}
asini@binary@constraint: {sma@binary@component} * (sin({incl@binary@component}))
t0_perpass@constraint: t0_supconj_to_perpass({t0_supconj@binary@component}, {period@binary@component}, {ecc@binary@component}, {per0@binary@component}, {dpdt@binary@component}, {dperdt@binary@component}, {t0@system})
t0_ref@constraint: t0_supconj_to_ref({t0_supconj@binary@component}, {period@binary@component}, {ecc@binary@component}, {per0@binary@component}, {dpdt@binary@component}, {dperdt@binary@component}, {t0@system})
mean_anom@constraint: (6.283185 * ({t0@system} - {t0_perpass@binary@component})) / {period@binary@component}
ecosw@constraint: {ecc@binary@component} * (cos({per0@binary@component}))
esinw@constraint: {ecc@binary@component} * (sin({per0@binary@component}))
freq@binary@constraint: 6.283185 / {period@binary@component}
mass@primary@constraint: (39.478418 * ({sma@binary@component} ** 3.000000)) / ((({period@binary@component} ** 2.000000) * ({q@binary@component} + 1.000000)) * 2942.206217504419328179210424423218)
sma@primary@constraint: {sma@binary@component} / ((1.000000 / {q@binary@component}) + 1.000000)
asini@primary@constraint: ({sma@binary@component} * (sin({incl@binary@component}))) / ((1.000000 / {q@binary@component}) + 1.000000)
period@primary@constraint: {period@binary@component} / {syncpar@primary@component}
incl@primary@constraint: {incl@binary@component} + {pitch@primary@component}
long_an@primary@constraint: {long_an@binary@component} + {yaw@primary@component}
requiv_max@primary@constraint: requiv_L1({q@binary@component}, {syncpar@primary@component}, {ecc@binary@component}, {sma@binary@component}, {incl@primary@component}, {long_an@primary@component}, {incl@binary@component}, {long_an@binary@component}, 1)
mass@secondary@constraint: (39.478418 * ({sma@binary@component} ** 3.000000)) / ((({period@binary@component} ** 2.000000) * ((1.000000 / {q@binary@component}) + 1.000000)) * 2942.206217504419328179210424423218)
sma@secondary@constraint: {sma@binary@component} / ({q@binary@component} + 1.000000)
asini@secondary@constraint: ({sma@binary@component} * (sin({incl@binary@component}))) / ({q@binary@component} + 1.000000)
period@secondary@constraint: {period@binary@component} / {syncpar@secondary@component}
incl@secondary@constraint: {incl@binary@component} + {pitch@secondary@component}
long_an@secondary@constraint: {long_an@binary@component} + {yaw@secondary@component}
requiv_max@secondary@constr...: requiv_L1({q@binary@component}, {syncpar@secondary@component}, {ecc@binary@component}, {sma@binary@component}, {incl@secondary@component}, {long_an@secondary@component}, {incl@binary@component}, {long_an@binary@component}, 2)
sample_from@phoebe01@compute: []
comments@phoebe01@compute:
dynamics_method@phoebe01@co...: keplerian
ltte@phoebe01@compute: False
irrad_method@phoebe01@compute: horvat
boosting_method@phoebe01@co...: none
eclipse_method@phoebe01@com...: native
horizon_method@phoebe01@com...: boolean
mesh_method@primary@phoebe0...: marching
mesh_method@secondary@phoeb...: marching
ntriangles@primary@phoebe01...: 1500
ntriangles@secondary@phoebe...: 1500
distortion_method@primary@p...: roche
distortion_method@secondary...: roche
atm@primary@phoebe01@compute: ck2004
atm@secondary@phoebe01@compute: ck2004
default_time_source@figure: None
latex_repr@primary@figure:
color@primary@figure: blue
marker@primary@figure: .
linestyle@primary@figure: solid
latex_repr@secondary@figure:
color@secondary@figure: orange
marker@secondary@figure: o
linestyle@secondary@figure: dashed
latex_repr@binary@figure:
R phoebe_version@setting: 2.3.44
dict_filter@setting: {}
dict_set_all@setting: False
run_checks_compute@setting: ['*']
run_checks_solver@setting: ['*']
run_checks_solution@setting: []
run_checks_figure@setting: []
auto_add_figure@setting: True
auto_remove_figure@setting: True
web_client@setting: False
web_client_url@setting: ui.phoebe-project.org
times@mylc@dataset: [0. 0.005 0.01 ... 0.99 0.995 1. ] d
fluxes@mylc@dataset: [] W / m2
passband@mylc@dataset: Johnson:V
intens_weighting@mylc@dataset: energy
compute_times@mylc@dataset: [] d
C compute_phases@mylc@dataset: []
phases_t0@mylc@dataset: t0_supconj
mask_enabled@mylc@dataset: True
mask_phases@mylc@dataset: []
solver_times@mylc@dataset: auto
sigmas@mylc@dataset: [] W / m2
pblum_mode@mylc@dataset: component-coupled
pblum_component@mylc@dataset: primary
l3_mode@mylc@dataset: flux
l3@mylc@dataset: 0.0 W / m2
exptime@mylc@dataset: 0.0 s
enabled@phoebe01@compute: True
fti_method@phoebe01@compute: none
ld_mode@primary@mylc@dataset: interp
ld_mode@secondary@mylc@dataset: interp
pblum@primary@mylc@dataset: 12.566370614359172 W
compute_phases@constraint: times_to_phases({compute_times@mylc@dataset}, {phases_period@mylc@dataset}, {period@binary@component}, {period_anom@binary@component}, {phases_dpdt@mylc@dataset}, {dpdt@binary@component}, {phases_t0@mylc@dataset}, {t0_supconj@binary@component}, {t0_perpass@binary@component}, {t0_ref@binary@component})
latex_repr@mylc@figure:
color@mylc@figure: black
marker@mylc@figure: .
linestyle@mylc@figure: solid
contexts@lcfig01@figure: ['*']
datasets@lcfig01@figure: ['*']
models@lcfig01@figure: ['*']
x@lcfig01@figure: times
y@lcfig01@figure: fluxes
xlabel_source@lcfig01@figure: auto
xunit_source@lcfig01@figure: auto
xlim_source@lcfig01@figure: auto
ylabel_source@lcfig01@figure: auto
yunit_source@lcfig01@figure: auto
ylim_source@lcfig01@figure: auto
color_source@lcfig01@figure: model
marker_source@lcfig01@figure: manual
marker@lcfig01@figure: .
linestyle_source@lcfig01@fi...: manual
linestyle@lcfig01@figure: solid
legend@lcfig01@figure: True
uncover@lcfig01@figure: False
highlight@lcfig01@figure: True
time_source@lcfig01@figure: default
color@latest@figure: blue
linestyle@latest@figure: solid
R times@latest@model: [0. 0.005 0.01 ... 0.99 0.995 1. ] d
R fluxes@latest@model: [0.38523668 0.45848033 0.5541843 ... 0.55433698
0.45851849 0.38523668] W / m2
comments@latest@model:
>>>

 

|

« Phoebe V2のインストール | トップページ | 趣味 »

コメント

コメントを書く



(ウェブ上には掲載しません)




« Phoebe V2のインストール | トップページ | 趣味 »