From 73fd54cc769c147c45ad92ba19f5a92d07c48e1f Mon Sep 17 00:00:00 2001 From: ftong <95+ftong@noreply.example.org> Date: Tue, 23 Jun 2026 13:09:39 +0200 Subject: [PATCH] Update src/seismic_hazard_forecasting.py add more logging --- src/seismic_hazard_forecasting.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py index b6d2daf..1385f43 100644 --- a/src/seismic_hazard_forecasting.py +++ b/src/seismic_hazard_forecasting.py @@ -524,6 +524,8 @@ verbose: {verbose}") x_range = np.linspace(grid_x_min, grid_x_max, nx) y_range = np.linspace(grid_y_min, grid_y_max, ny) + logger.debug(f"Grid X range: {x_range}, Y range: {y_range}") + t_windowed = time r_windowed = [[x, y]] @@ -797,7 +799,7 @@ verbose: {verbose}") end = timer() logger.info(f"Ground motion exceedance computation time: {round(end - start, 1)} seconds") - logger.info(f"IMT values: {iml_grid_raw[0]}") + logger.debug(f"IMT values: {iml_grid_raw[0]}") if np.isnan(iml_grid_raw).all(): msg = "No valid ground motion intensity measures were forecasted. Try a different ground motion model." @@ -868,6 +870,8 @@ verbose: {verbose}") east, west = grid_lon_max, grid_lon_min # Longitude range bounds = [[south, west], [north, east]] + logger.debug(f"Grid extent North {north}, South {south}, East {east}, West {west}") + map_center = [np.mean([north, south]), np.mean([east, west])] # Create an image from the grid