From 5325c074041a5307ca7e8957d555dc306fdeecd8 Mon Sep 17 00:00:00 2001
From: ftong <95+ftong@noreply.example.org>
Date: Wed, 8 Jul 2026 17:16:21 +0200
Subject: [PATCH] Update src/seismic_hazard_forecasting.py

AOI feature makes exclude_low_fxy redundant as we now have way to save time by selecting small area to forecast
---
 src/seismic_hazard_forecasting.py | 4 ----
 1 file changed, 4 deletions(-)

diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 572321c..100dc14 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -342,9 +342,6 @@ def main(catalog_file, mc_file, pdf_file, m_file, m_select, mag_label, mc, m_max
     else:
         logger.setLevel(logging.INFO)
 
-    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:])
     
@@ -482,7 +479,6 @@ verbose: {verbose}")
         # Handle Area of Interest (AOI)
         if (None not in AOI_lat) and (None not in AOI_lon):
             use_AOI = True
-            exclude_low_fxy = False # don't exclude any points because we need to analyze all grid points in the AOI
             
             # Create an AOI GeoDataFrame and project it to the same UTM CRS
             aoi_gdf_utm = gpd.GeoDataFrame(