diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 4c7821d..a144433 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -649,9 +649,9 @@ verbose: {verbose}")
         logger.info(f"95% BCa confidence interval: [{rate_unc_low:.4f}, {rate_unc_high:.4f}]")
         logger.info("------------------------------------------------------")
         
-        lambdas = np.atleast_1d(np.array(rate_forecast/bin_dur, dtype='d'))
-        lambdas_perc = np.array(1.0, dtype='d')
-
+        lambdas = np.array([rate_forecast/bin_dur], dtype='d')
+        lambdas_perc = np.array([1], dtype='d')
+        
         np.savetxt('activity_rate.csv', lambdas, header=f"Activity Rate (Events per {time_unit[:-1]})",
                    delimiter=',', fmt='%1.4f')