Skip to content

Commit

Permalink
Fixed BA19 year.
Browse files Browse the repository at this point in the history
  • Loading branch information
arkottke committed Jan 22, 2022
1 parent aeea2ea commit e308998
Show file tree
Hide file tree
Showing 3 changed files with 3 additions and 4 deletions.
Binary file removed examples/example-ba18.png
Binary file not shown.
Binary file added examples/example-ba19.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
7 changes: 3 additions & 4 deletions examples/example-ba18.py → examples/example-ba19.py
Original file line number Diff line number Diff line change
Expand Up @@ -8,17 +8,16 @@
scenario = pygmm.Scenario(mag=6.5, dist_rup=30, v_s30=500, mechanism="SS", depth_tor=0)

models = [
pygmm.BaylessAbrahamson18(scenario),
pygmm.BaylessAbrahamson18(
scenario.copy_with(v_s30=pygmm.BaylessAbrahamson18.V_REF)
pygmm.BaylessAbrahamson2019(scenario),
pygmm.BaylessAbrahamson2019(
scenario.copy_with(v_s30=pygmm.BaylessAbrahamson2019.V_REF)
),
]

fig, ax = plt.subplots()

for m in models:
ax.plot(m.freqs, m.eas, label=f"{m.scenario.v_s30:.0f}")
print(m.freqs, m.eas)

ax.legend(title="$V_{s30}$ m/s")

Expand Down

0 comments on commit e308998

Please sign in to comment.