Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

AELO year 3 updates to run all IMTs #9648

Open
wants to merge 32 commits into
base: master
Choose a base branch
from
Open

AELO year 3 updates to run all IMTs #9648

wants to merge 32 commits into from

Conversation

AnnaRood
Copy link

@AnnaRood AnnaRood commented May 8, 2024

Updated AELO Year 1 IMTs [PGA, SA(0.2), and SA(1.0)] to run Year 3 [PGA, SA(0.02) -> SA(10)].
[Note: no SA(0) yet]

Updated to ASCE7-22 parameters.

Plots updated to show:

  • One plot for each IMT for the disagg by src (for us)
  • 3-piece disagg plot kept as PGA, 0.2 s , 1.0 s
  • Mean hazard curves for a subset of IMTs, and moved from maximum component to geometric mean hazard curves.
  • MCE results for all periods, and two subplots for a log and linear scale

@micheles micheles added this to the Engine 3.20.0 milestone May 9, 2024
IMTS = ['PGA', 'SA(0.02)', 'SA(0.03)', 'SA(0.05)', 'SA(0.075)', 'SA(0.1)',
'SA(0.15)', 'SA(0.2)', 'SA(0.25)', 'SA(0.3)', 'SA(0.4)', 'SA(0.5)',
'SA(0.75)', 'SA(1.0)', 'SA(1.5)', 'SA(2.0)', 'SA(3.0)', 'SA(4.0)',
'SA(5.0)', 'SA(7.5)', 'SA(10.0)']
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I thought most GMPEs could not be extended to SA(10.0) and we decided to discard it (?)

Copy link
Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

We have now updated so that the IMTS are defined in the job file, instead of using defaults

@AnnaRood AnnaRood changed the title [WIP] AELO year 3 updates to run all IMTs AELO year 3 updates to run all IMTs May 14, 2024
Copy link
Contributor

@kejohnso kejohnso left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

LGTM!

@@ -115,7 +110,9 @@ def calc_rtgm_df(hcurves, site, site_idx, oq):
:returns: pandas dataframe with RTGM and related parameters
:returns: numpy array of conversion factors to max component
"""
IMTS = list(oq.imtls)
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The conventions is to use all caps for constants, so IMTS is a bad name. You can use imts or job_imts or orig_imts or whatever you like more.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants