From 2b6cd521313867eed32ef42d7445fc669de833f1 Mon Sep 17 00:00:00 2001 From: ftong <95+ftong@noreply.example.org> Date: Wed, 24 Jun 2026 11:39:03 +0200 Subject: [PATCH] Update src/seismic_hazard_forecasting.py --- src/seismic_hazard_forecasting.py | 14 +++++++------- 1 file changed, 7 insertions(+), 7 deletions(-) diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py index 4b768fb..220a193 100644 --- a/src/seismic_hazard_forecasting.py +++ b/src/seismic_hazard_forecasting.py @@ -7,13 +7,13 @@ from matplotlib.ticker import MultipleLocator import matplotlib.pyplot as plt import numpy as np -global ncp_choice, tcp_max, torder_min, torder_max, AOI_lat, AOI_lon +global ncp_choice, tcp_max, torder_min, torder_max ncp_choice = 'default' tcp_max = 5 torder_min = 0 torder_max = 1 -AOI_lat = np.array([51.48, 51.54]) -AOI_lon = np.array([16.15, 16.24]) +#AOI_lat = np.array([51.48, 51.54]) +#AOI_lon = np.array([16.15, 16.24]) #AOI_lat = np.array([None, None]) #AOI_lon = np.array([None, None]) @@ -253,8 +253,8 @@ def bins_and_beast(dates, unit, bin_dur, multiplicator): def main(catalog_file, mc_file, pdf_file, m_file, m_select, mag_label, mc, m_max, m_kde_method, xy_select, grid_dim, xy_win_method, rate_select, time_win_duration, - forecast_select, custom_rate, forecast_len, time_unit, model, products_string, verbose): -# forecast_select, custom_rate, forecast_len, time_unit, AOI_extent, model, products_string, verbose): +# forecast_select, custom_rate, forecast_len, time_unit, model, products_string, verbose): + forecast_select, custom_rate, forecast_len, time_unit, AOI_extent, model, products_string, verbose): """ Python application that reads an earthquake catalog and performs seismic hazard forecasting. Arguments: @@ -344,8 +344,8 @@ def main(catalog_file, mc_file, pdf_file, m_file, m_select, mag_label, mc, m_max exclude_low_fxy = False # skip low probability areas of the map thresh_fxy = 1e-3 # minimum fxy value (location PDF) needed to do PGA estimation (to skip low probability areas); also should scale according to number of grid points -# AOI_lat = np.array(AOI_extent[:2]) -# AOI_lon = np.array(AOI_extent[2:]) + AOI_lat = np.array(AOI_extent[:2]) + AOI_lon = np.array(AOI_extent[2:]) # log user selections logger.debug(f"User input files\n Catalog: {catalog_file}\n Mc: {mc_file}\n Mag_PDF: {pdf_file}\n Mag: {m_file}")