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

STOP: EVOLV1 FATAL ERROR #39

Closed pavanvyn closed 3 years ago

pavanvyn commented 3 years ago

Version 3ae9006280819f1bc190ce06e5191bb669a90876

configuration = "[2,2]"
N_bodies = 4
masses = [4.844524218671435, 8.11891899747903, 7.833078005839005, 0.46753945289745924]
metallicities = [0.02, 0.02, 0.02, 0.02]
semimajor_axes = [12.62371561001661, 6.004199751613622, 6736.579264887057]
eccentricities = [0.764259851004382, 0.383047422054649, 0.8409918894045305]
inclinations = [1.011816183280177, 2.0220865953864893, 2.022218711344736]
arguments_of_pericentre = [2.0056590053492602, 2.7923321965459107, 2.240399302614785]
longitudes_of_ascending_node = [3.362167998015483, 0.7757501917277166, 3.640222539487186]
tend = 1.400000E+10
N_steps = 100000
stellar_types = [1, 1, 1, 1]
object_types = [1, 1, 1, 1]
hamers commented 3 years ago

For a merged star, SSE failed to converge when finding an internal time-step such that the radius increase would be less than 10%. I believe this is a very rare occurrence, and to get around it, I changed evolv1.f to now stop the search if unable to converge in the new time-step (5007b4e). This seems to not produce unexpected behaviour in the subsequent evolution, at least in this instance.