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 4 occurred -- stopping simulation. #50

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version cd8184b902964d9134eaae26f33db86a4694e31e (v0.83)

configuration = "[2,2]"
N_bodies = 4
masses = [11.936308470789099, 7.547039803563019, 10.485550503630007, 14.9226023204833]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [0.7373304083674703, 19.468351016464595, 301.10897805970336]
eccentricities = [0.07776443321165889, 0.6847039014680483, 0.6227576939518841]
inclinations = [0.444389012062038, 1.725357029474652, 2.620807358435491]
arguments_of_pericentre = [3.523291155163777, 4.4182808111681044, 0.4903957377732229]
longitudes_of_ascending_node = [1.3104762524702338, 3.235687730517438, 4.081576882875382]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
kick_distribution_sigma_km_s_NS,kick_distribution_sigma_km_s_BH = 0.0,0.0

configuration = "[2,2]"
N_bodies = 4
masses = [1.6337969506029681, 4.2750927137387515, 11.877354130165823, 10.288911109538986]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [24.278312077726156, 15.8970375523649, 1793.9673042722384]
eccentricities = [0.5521380989079157, 0.041501194627405105, 0.03098032451290733]
inclinations = [2.744543967647522, 0.8542249583240987, 2.114534648076657]
arguments_of_pericentre = [1.5503620961209863, 1.3694357099696468, 3.042271831766447]
longitudes_of_ascending_node = [4.519472526396828, 2.0803372971183856, 2.2168615570984365]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
kick_distribution_sigma_km_s_NS,kick_distribution_sigma_km_s_BH = 0.0,0.0

configuration = "[2,2]"
N_bodies = 4
masses = [10.155717437998788, 10.694291523444676, 9.377971441266535, 9.393160984603881]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [326.68621136397473, 55.380787636848346, 9718.700044479629]
eccentricities = [0.7676894566629486, 0.7164696698248222, 0.07778716556753808]
inclinations = [0.8839497330900175, 2.097147399211663, 2.6356261816685924]
arguments_of_pericentre = [3.561688172958032, 4.457275773557397, 4.630435476178859]
longitudes_of_ascending_node = [1.3847490327754812, 2.2650396692046706, 1.9912735600647105]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
kick_distribution_sigma_km_s_NS,kick_distribution_sigma_km_s_BH = 0.0,0.0

These systems end with the error:

WARNING -- mse.py -- Internal error with code  4 occurred -- -- stopping the simulation but saving data/making plots if specified in command line arguments.
pavanvyn commented 3 years ago

Solved in commit 10156e4ba8fc8973128cdea6e98eb9f8054155f3