diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index ee252b0..c9aa165 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -596,10 +596,13 @@ verbose: {verbose}")
     elif rate_select:
         logger.info(f"Activity rate modeling selected")
                 
-        time, mag_data, lat_dummy, lon_dummy, depth_dummy = read_mat_cat(catalog_file, mag_label=mag_label, output_datenum=True)
-
-        datenum_data = time  # REMEMBER THE DECIMAL DENOTES DAYS
+        datenum_data, mag_data, lat_dummy, lon_dummy, depth_dummy = read_mat_cat(catalog_file, mag_label=mag_label, output_datenum=True)
 
+        if trim_to_mc:
+            indices = np.argwhere(mag_data < mc)
+            mag_data = np.delete(mag_data, indices)
+            datenum_data = np.delete(datenum_data, indices)
+        
         if time_unit == 'hours':
             multiplicator = 24
         elif time_unit == 'days':