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 6 occurred -- stopping simulation. #42

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 6e960c32d6deba8b3420b8483fef1ede37819fef

These are the systems I mentioned in our last meeting, the ones with the collision error.

configuration = "2"
N_bodies = 2
masses = [3.0570774984253286, 3.4335505060992015]
metallicities = [0.02, 0.02]
semimajor_axes = [0.10308355683896057]
eccentricities = [0.4772310671310243]
inclinations = [1.2968326028476609]
arguments_of_pericentre = [6.203836066580958]
longitudes_of_ascending_node = [4.0540468980697595]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]

configuration = "2"
N_bodies = 2
masses = [1.465998881349993, 1.2309546317228084]
metallicities = [0.02, 0.02]
semimajor_axes = [0.20439504179243606]
eccentricities = [0.887249092698013]
inclinations = [1.8833950553687608]
arguments_of_pericentre = [3.260033210427437]
longitudes_of_ascending_node = [5.828964710110773]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1]
object_types = [1, 1]
hamers commented 3 years ago

Contact evolution occurred in these systems between a giant and non-giant, and handle_mass_transfer_cases() did not handle this correctly (invoking collision_product instead of CE). Should be fixed in 9fc01dd.