official-apps/SeismicHazardForecasting
11
0
Fork 2
Code Issues Pull Requests 1 Activity

custom AOI feature, changes to plotting and activity rate calculation #24

Merged
ymlesni merged 50 commits from July2026 into master 2026-07-22 12:43:46 +02:00
Conversation 0 Commits 50 Files Changed 2
+5 -1
1 changed files with 5 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files
Showing only changes of commit 73fd54cc76 - Show all commits
src/seismic_hazard_forecasting.py
+5 -1
View File
@@ -524,6 +524,8 @@ verbose: {verbose}")
x_range = np.linspace(grid_x_min, grid_x_max, nx) x_range = np.linspace(grid_x_min, grid_x_max, nx)
y_range = np.linspace(grid_y_min, grid_y_max, ny) 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 t_windowed = time
r_windowed = [[x, y]] r_windowed = [[x, y]]
@@ -797,7 +799,7 @@ verbose: {verbose}")
end = timer() end = timer()
logger.info(f"Ground motion exceedance computation time: {round(end - start, 1)} seconds") 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(): if np.isnan(iml_grid_raw).all():
msg = "No valid ground motion intensity measures were forecasted. Try a different ground motion model." 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 east, west = grid_lon_max, grid_lon_min # Longitude range
bounds = [[south, west], [north, east]] 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])] map_center = [np.mean([north, south]), np.mean([east, west])]
# Create an image from the grid # Create an image from the grid