Here we'll examine how well jktebop agrees with PHOEBE with increased distortion.
In [1]:
import phoebe
In [2]:
b = phoebe.default_binary()
In [3]:
b.add_dataset('lc', compute_phases=phoebe.linspace(0,1,101))
Out[3]:
In order to allow jktebop to compute models, we'll set requiv_max_limit=1.0, effectively disabling the error that would otherwise be raised at a default factor of 0.5 by b.run_checks_compute.
In [4]:
b.add_compute('jktebop', requiv_max_limit=1.0)
Out[4]:
And to avoid any issues with falling outside the atmosphere grids, we'll set a simple flat limb-darkening model and disable irradiation.
In [5]:
b.set_value_all('ld_mode', 'manual')
b.set_value_all('ld_func', 'linear')
b.set_value_all('ld_coeffs', [0.5])
b.set_value_all('irrad_method', 'none')
For PHOEBE, we'll use blackbody atmospheres (again to avoid any issues of falling out of the grid). For jktebop, we'll keep 'ck2004' - this will only be used to compute the flux-scaling factor based on mean stellar values, so should not fall outside the grid.
At a quick glance, we can see the jktebop agrees quite well at a factor of 0.6, but noticeable differences appear by 0.7 (keep in mind, the default value before an error will be raise within PHOEBE is 0.5, but this can be adjusted as necessary, with caution).
In [6]:
requiv_max = b.get_value('requiv_max', component='primary', context='component')
for requiv_max_factor in [0.6, 0.7, 0.8, 0.9, 1.0]:
b.set_value('requiv', component='primary', value=requiv_max_factor*requiv_max)
b.run_compute(kind='phoebe', atm='blackbody', model='phoebe2_model', overwrite=True)
b.run_compute(kind='jktebop', model='jktebop_model', overwrite=True)
_ = b.plot(context='model',
title='requiv = {:0.1f} / requiv_max'.format(requiv_max_factor),
draw_title=True,
legend=True, show=True)
In [ ]: