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(