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)