From 144f7e1fcdb623c103062dbfcf89369c03933152 Mon Sep 17 00:00:00 2001 From: ftong <95+ftong@noreply.example.org> Date: Mon, 22 Jun 2026 23:04:07 +0200 Subject: [PATCH] Update src/seismic_hazard_forecasting.py remove events below mc for activity rate --- src/seismic_hazard_forecasting.py | 9 ++++++--- 1 file changed, 6 insertions(+), 3 deletions(-) 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':