From d5e5435d8348b5d9c9db2312b6e799b8df8541e4 Mon Sep 17 00:00:00 2001 From: ftong <95+ftong@noreply.example.org> Date: Mon, 22 Jun 2026 23:36:33 +0200 Subject: [PATCH] Update src/seismic_hazard_forecasting.py --- src/seismic_hazard_forecasting.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py index c9aa165..cef5655 100644 --- a/src/seismic_hazard_forecasting.py +++ b/src/seismic_hazard_forecasting.py @@ -749,8 +749,6 @@ verbose: {verbose}") # For each receiver compute estimated ground motion values logger.info(f"Estimating ground motion intensity at {len(indices_filtered)} grid points...") - PGA = np.zeros(shape=(nx * ny)) - start = timer() use_pp = True @@ -796,6 +794,8 @@ 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]}") + if np.isnan(iml_grid_raw).all(): msg = "No valid ground motion intensity measures were forecasted. Try a different ground motion model." logger.error(msg)