search

hamers / mse

Multiple Stellar Evolution -- A Population Synthesis Code for Multiple-Star Systems
GNU General Public License v3.0
9 stars 4 forks source link

mse.py -- Internal error with code 24 occurred -- stopping simulation. #37

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 3c1483876dfc7045519b13ad8568378a7a36b9a6

configuration = "[2,2]"
N_bodies = 4
masses = [11.545724854578776, 1.0881783040872224, 3.3115343779347834, 6.010473778620222]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [218.78017428049276, 153.5405019638936, 703.9382128960965]
eccentricities = [0.8217645445206639, 0.6461807039125479, 0.07027389562515694]
inclinations = [1.897725151813237, 1.482110956942365, 1.672984517289217]
arguments_of_pericentre = [5.215497922485309, 5.6057896113067995, 1.5135481989561066]
longitudes_of_ascending_node = [2.855244987772052, 3.9616874484784317, 0.1728142629395301]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]

(Sorry for posting this after 5pm! Found an error after a while)

hamers commented 3 years ago

This is a rare case of a marginally stable system in which orbit crossing occurred before secular breakdown/MA11 instability was triggered, which was not taken into account in the code. Fixed in dd220f5.