def main():
o, arguments = new_option_parser().parse_args()
double_star, stars = make_binary_star(
o.mprim, o.msec, o.semimajor_axis, o.eccentricity,
)
end_time = 1000.0 | units.yr
evolve_model(end_time, double_star, stars)
def main():
o, arguments = new_option_parser().parse_args()
# override binary star parameters
o.semimajor_axis = 1 | units.au
o.eccentricity = 0.5
o.mprim = o.msec = 30 | units.MSun
double_star, stars = make_binary_star(
o.mprim,
o.msec,
o.semimajor_axis,
o.eccentricity,
)
stars.kaps = o.kaps
stars.taulag = o.taulag
# override star parameters
stars.radius = 25 | units.RSun
stars[0].kaps = 0.15
stars[0].taulag = 1e4 | units.s
stars[1].kaps = 0.0
stars[1].taulag = 0.0 | units.s
# this circularizes in 5e4 yr
end_time = 1e4 | units.yr
evolve_model(end_time, double_star, stars)
def main():
m1, m2 = 80 | units.MSun, 55 | units.MSun
a0 = 4000 | units.RSun
e0 = 0.9
o, arguments = new_option_parser(M_default=m1,
m_default=m2,
a_default=a0,
e_default=e0).parse_args()
double_star, stars = make_binary_star(
o.mprim,
o.msec,
o.semimajor_axis,
o.eccentricity,
)
end_time = 1000.0 | units.yr
evolve_model(end_time, double_star, stars)