diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 429fb37..4c7821d 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -649,10 +649,10 @@ verbose: {verbose}")
         logger.info(f"95% BCa confidence interval: [{rate_unc_low:.4f}, {rate_unc_high:.4f}]")
         logger.info("------------------------------------------------------")
         
-        lambdas = np.array(rate_forecast/bin_dur, dtype='d')
+        lambdas = np.atleast_1d(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[:-1]})",
+        np.savetxt('activity_rate.csv', lambdas, header=f"Activity Rate (Events per {time_unit[:-1]})",
                    delimiter=',', fmt='%1.4f')
 
     if forecast_select: