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

Update src/seismic_hazard_forecasting.py

Browse Source
This commit is contained in:
ftong 2025-07-01 11:55:34 +02:00
parent 52d2c5b0f7
commit 09e21014d1
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
src/seismic_hazard_forecasting.py
View File

@ -116,7 +116,7 @@ def calc_bins2(dates, unit, bin_dur, dates_calc, rate_forecast, rate_unc_high, r
rt, idx, u_e = av_rates(act_rate, pprs) rt, idx, u_e = av_rates(act_rate, pprs)
# --------------------- Line Plot ------------------------- # --------------------- Line Plot -------------------------
fig, ax = plt.subplots(dpi= 1000, figsize=(14, 5)) fig, ax = plt.subplots(figsize=(14, 5))
# Plot activiry rate for all catalog # Plot activiry rate for all catalog
ax.plot(bin_edges_dt[1:len(bin_edges)], act_rate, '-o', linewidth=2, markersize=6) ax.plot(bin_edges_dt[1:len(bin_edges)], act_rate, '-o', linewidth=2, markersize=6)
@ -177,7 +177,7 @@ def calc_bins2(dates, unit, bin_dur, dates_calc, rate_forecast, rate_unc_high, r
[rt[i]] * (int(u_e[idx[i]]) - int(u_e[idx[i - 1]])), linewidth=2) [rt[i]] * (int(u_e[idx[i]]) - int(u_e[idx[i - 1]])), linewidth=2)
plt.legend(loc='best') plt.legend(loc='best')
plt.savefig(filename, dpi=600) plt.savefig(filename)
plt.show() plt.show()
return act_rate, bin_counts, bin_edges, out, pprs, rt, idx, u_e return act_rate, bin_counts, bin_edges, out, pprs, rt, idx, u_e