diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 369ed95..e1fa212 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -636,24 +636,15 @@ verbose: {verbose}")
                       rate_forecast, rate_unc_high, rate_unc_low,
                       datenum_data, mag_data)
 
-        print("\n----------------- Forecast Summary -----------------")
-        print(f"Forecasted activity rate (next {bin_dur} {time_unit}(s)): {rate_forecast:.4f}")
-        print(f"95% BCa confidence interval: [{rate_unc_low:.4f}, {rate_unc_high:.4f}]")
-        print("------------------------------------------------------")
+        logger.info("\n----------------- Forecast Summary -----------------")
+        logger.info(f"Forecasted activity rate (next {bin_dur} {time_unit}(s)): {rate_forecast:.4f}")
+        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_perc = np.array(1.0, dtype='d')
 
-
-
-
-
-
-
-        
-        
-        logger.info(f"Forecasted activity rates: {lambdas} events per {time_unit}")
-        np.savetxt('activity_rate.csv', np.vstack((lambdas, lambdas_perc)).T, header="lambda, percentage",
+        np.savetxt('activity_rate.csv', np.vstack((lambdas, lambdas_perc)).T, header=f"Activity Rate (Events per {time_unit})", "percentage",
                    delimiter=',', fmt='%1.4f')
 
     if forecast_select: