butterflies: more information in first year

src/modules/butterflies/butterflies_actions.c

@@ -42,6 +42,7 @@ enum action_status bf_action_stats_reset(struct cats_grid *grid, struct cats_con
                 cats_dt_rates egg_fraction_source = module_conf->egg_fraction_source;
                 cats_dt_rates eggs_per_female = calculate_rate_for_matrix(&module_conf->eggs_per_female, &l_param, print_rate);
                 cats_dt_rates reproduction_rate = calculate_rate_for_matrix(&module_conf->reproduction_rate, &l_param, print_rate);
+                cats_dt_rates K = calculate_rate_for_matrix(&conf->param[grid->id].carrying_capacity, &l_param, print_rate) * conf->param->max_adult_cc_fraction;
                 cats_dt_rates local_eggs =  (stationary + mobile * egg_fraction_source ) * eggs_per_female;
                 // female -> female
                 // to achieve the target reproduction rate, the number of eggs per female laid in the cell
@@ -56,9 +57,26 @@ enum action_status bf_action_stats_reset(struct cats_grid *grid, struct cats_con
                 log_message(LOG_INFO, "SUMMARY: stationary females at OT: %Lf", stationary);
                 log_message(LOG_INFO, "SUMMARY: eggs to adult rate at OT: %Lf", eggs_to_adults_rate);
                 log_message(LOG_INFO, "SUMMARY: egg fraction source at OT: %Lf", egg_fraction_source);
-                log_message(LOG_INFO, "SUMMARY: carrying capacity at OT: %d", l_param.K);
+                log_message(LOG_INFO, "SUMMARY: carrying capacity at OT: %Lf", K);
                 log_message(LOG_INFO, "SUMMARY: effective female to female rate at OT: %Lf\n", result);
 
+                l_param.suitability = 1.0;
+                eggs_per_female = calculate_rate_for_matrix(&module_conf->eggs_per_female, &l_param, print_rate);
+                reproduction_rate = calculate_rate_for_matrix(&module_conf->reproduction_rate, &l_param, print_rate);
+                K = calculate_rate_for_matrix(&conf->param[grid->id].carrying_capacity, &l_param, print_rate) * conf->param->max_adult_cc_fraction;
+                local_eggs =  (stationary + mobile * egg_fraction_source ) * eggs_per_female;
+                eggs_to_adults_rate = bf_egg_to_adult_survival_rate(reproduction_rate, local_eggs, female_fraction);
+                result =  local_eggs * eggs_to_adults_rate * female_fraction;
+
+                log_message(LOG_INFO, "SUMMARY: reproduction rate at suitability 1: %Lf", reproduction_rate);
+                log_message(LOG_INFO, "SUMMARY: eggs per female at suitability 1: %Lf", eggs_per_female);
+                log_message(LOG_INFO, "SUMMARY: local eggs at suitability 1: %Lf", local_eggs);
+                log_message(LOG_INFO, "SUMMARY: stationary females at suitability 1: %Lf", stationary);
+                log_message(LOG_INFO, "SUMMARY: eggs to adult rate at suitability 1: %Lf", eggs_to_adults_rate);
+                log_message(LOG_INFO, "SUMMARY: egg fraction source at suitability 1: %Lf", egg_fraction_source);
+                log_message(LOG_INFO, "SUMMARY: carrying capacity at suitability 1: %Lf", K);
+                log_message(LOG_INFO, "SUMMARY: effective female to female rate at suitability 1: %Lf\n", result);
+
         }
 
         return action_grid_stats_reset(grid, conf);