diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 3de67d9..b980171 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -652,7 +652,7 @@ verbose: {verbose}")
         lambdas = np.array(rate_forecast/bin_dur, dtype='d')
         lambdas_perc = np.array(1.0, dtype='d')
 
-        np.savetxt('activity_rate.csv', np.atleast_1d(lambdas), header=f"Activity Rate (Events per {time_unit})",
+        np.savetxt('activity_rate.csv', np.atleast_1d(lambdas), header=f"Activity Rate (Events per {time_unit[:-1]})",
                    delimiter=',', fmt='%1.4f')
 
     if forecast_select: