diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 3c0bb9d..3417616 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -17,7 +17,7 @@ def main(catalog_file, mc_file, pdf_file, m_file, m_select, mag_label, mc, m_max
                     then the program looks for a label of 'Mw' for magnitude in the catalog
         mc: The magnitude of completeness (Mc) of the catalog
         m_max:M_max. The magnitude distribution is estimated for the range from Mc to M_max. If no value is provided,
-                    then the program sets M_max to be 3 magnitude units above the maximum magnitude value in the catalog.
+                    then the program sets M_max to be 1 magnitude units above the maximum magnitude value in the catalog.
         m_kde_method: The kernel density estimator to use.
         xy_select: If True, perform an estimation of the magnitude distribution using KDE with the chosen KDE method
         grid_dim: The grid cell size (in metres) of the final ground motion product map. A smaller cell size will
@@ -160,9 +160,10 @@ verbose: {verbose}")
             lon = np.delete(lon, indices)
 
         # if user does not provide a m_max, set m_max to 1 magnitude unit above max magnitude in catalog
+        logger.info(f"mag max of catalog is: {mag.max()}")
         if m_max == None:
             m_max = mag.max() + 1.0
-
+        logger.info(f"m_max is: {m_max}")
         start = timer()
 
         t_windowed, r_windowed = win_CNE(time, [lon, lat, mag], win_size=len(mag), win_overlap=0, min_events=1)