Update src/seismic_hazard_forecasting.py

remove events below mc for activity rate
2026-06-22 23:04:07 +02:00
parent 337457d49c
commit 144f7e1fcd
1 changed files with 6 additions and 3 deletions
@@ -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':