diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index d02200d..0625ae3 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
-        np.save('colorbar_0.svg', np.array(datenum_data))
+        np.save('KDE_magnitude_PDF.csv', np.array(datenum_data))
         continue
         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,