diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 9a58f6a..965b8ec 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -361,7 +361,7 @@ verbose: {verbose}")
         rate_forecast = multiplicator / np.median(forecasts)  # [per time unit]
 
         # Plot of forecasted activity rate with previous binned activity rate
-        logger.debug(f"{datenum_data}")
+        np.save('datenum_array.npy', np.array(datenum_data))
         act_rate, bin_counts, bin_edges, out, pprs, rt, idx, u_e = calc_bins(np.array(datenum_data), time_unit,
                                                                              time_win_duration, dates_calc,
                                                                              rate_forecast, rate_unc_high, rate_unc_low,