Preserve stellar properties when restarting with SeBa (#1105)

src/amuse/community/seba/interface.cc

@@ -464,10 +464,10 @@ int new_particle(int * index_of_the_star, double mass){
 }
 
 
-int new_advanced_particle(int * index_of_the_star, double mass,  double relative_mass, int type_number,  double age, double core_mass, double COcore_mass,  double radius){
+int new_advanced_particle(int * index_of_the_star, double mass,  double relative_mass, int type_number,  double relative_age, double core_mass, double COcore_mass,  double radius){
 
     if (relative_mass == 0) return new_particle(index_of_the_star, mass);
-    if (age < 0) return -1;
+    if (relative_age < 0) return -1;
 
     node * new_node = new node();
     new_node->set_label(next_seba_id);
@@ -484,20 +484,19 @@ int new_advanced_particle(int * index_of_the_star, double mass,  double relative
         seba_insertion_point = new_node;
     }
 
-    stellar_type seba_stellar_type = translate_int_to_stellar_type(type_number);    
+    stellar_type seba_stellar_type = translate_int_to_stellar_type(type_number);
 
-    addstar(new_node, seba_time, seba_stellar_type, seba_metallicity, 0, false);
+    addstar(new_node, seba_time, seba_stellar_type, seba_metallicity, 0, false, seba_stellar_type, relative_mass, mass - core_mass, core_mass, COcore_mass, relative_age);
     new_node->get_starbase()->set_time_offset(seba_time);
     *index_of_the_star = next_seba_id;
 
     next_seba_id++;
 
-    new_node->get_starbase()->set_relative_age(age);
+    new_node->get_starbase()->set_relative_age(relative_age);
     new_node->get_starbase()->set_core_mass(core_mass);
     new_node->get_starbase()->set_COcore_mass(COcore_mass);
     new_node->get_starbase()->set_effective_radius(radius);
 
-
     return 0;
 }
 

src/amuse/community/seba/interface.py

@@ -69,8 +69,8 @@ def new_advanced_particle():
             'stellar_type', dtype='int32', direction=function.IN, default=0,
             description="The initial stellar type of the star")
         function.addParameter(
-            'age', dtype='float64', direction=function.IN, default=0,
-            description="The initial age of the star")
+            'relative_age', dtype='float64', direction=function.IN, default=0,
+            description="The relative initial age of the star")
         function.addParameter(
             'core_mass', dtype='float64', direction=function.IN, default=0, 
             description="The initial core mass of the star")