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 1 occurred -- stopping simulation. #7 #41

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 8c5e2fb63a6a13396af2611c29dcc0fa6f88c024

configuration = "[2,2]"
N_bodies = 4
masses = [3.548986035679437, 17.021172984446476, 8.12233857680884, 4.428527949432436]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.2015380813653181, 5.604271111382169, 53.586871789366896]
eccentricities = [0.3527075005652079, 0.9102901776980962, 0.47208484125535405]
inclinations = [1.3308787041010495, 1.9977131616727142, 1.6271690533278615]
arguments_of_pericentre = [6.227825330298463, 4.316756296360188, 4.147507183112485]
longitudes_of_ascending_node = [0.9923061086882445, 3.28162665237633, 0.7969508437343856]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
hamers commented 3 years ago

A SNe event triggered direct collision but the RLOF routines were invoked, causing too large mass transfer rates and ultimately NaNs. This is now checked for and fixed in 5412814.