diff --git a/examples/main.c b/examples/main.c index db7df3b4a129b5469758485c7a042af173e63123..7b00fb4eadef4f7a1819b5306a80e81040af4c15 100644 --- a/examples/main.c +++ b/examples/main.c @@ -90,6 +90,7 @@ int main(int argc, char *argv[]) { struct cooling_function_data cooling_func; struct cosmology cosmo; struct external_potential potential; + struct star_formation starform; struct pm_mesh mesh; struct gpart *gparts = NULL; struct gravity_props gravity_properties; @@ -855,6 +856,11 @@ int main(int argc, char *argv[]) { if (with_cooling) cooling_init(params, &us, &prog_const, &cooling_func); if (myrank == 0) cooling_print(&cooling_func); + /* Initialise the star formation law and its properties */ + if (with_star_formation) + starformation_init(params, &prog_const, &us, &s, &starform); + if (myrank ==0) starformation_init(&starform); + /* Initialise the chemistry */ chemistry_init(params, &us, &prog_const, &chemistry); if (myrank == 0) chemistry_print(&chemistry); @@ -881,7 +887,7 @@ int main(int argc, char *argv[]) { engine_init(&e, &s, params, N_total[0], N_total[1], N_total[2], engine_policies, talking, &reparttype, &us, &prog_const, &cosmo, &hydro_properties, &gravity_properties, &stars_properties, - &mesh, &potential, &cooling_func, &chemistry); + &mesh, &potential, &cooling_func, &starform, &chemistry); engine_config(0, &e, params, nr_nodes, myrank, nr_threads, with_aff, talking, restart_file); diff --git a/src/engine.h b/src/engine.h index d7d916c1bc25fc0d2cbac26c95de406accec5a67..f14f8810fe866f8b3ea079061fdc6cd199969bf4 100644 --- a/src/engine.h +++ b/src/engine.h @@ -355,6 +355,9 @@ struct engine { /* Properties of the cooling scheme */ struct cooling_function_data *cooling_func; + /* Properties of the starformation law */ + const struct star_formation *star_formation; + /* Properties of the chemistry model */ const struct chemistry_global_data *chemistry;