SimCenter Forum

Research in Natural Hazards Engineering => Regional Hazard Simulation (R2D, rWhale) => Topic started by: Anne Hulsey on September 21, 2021, 12:24:58 AM

Title: openquake poe values
Post by: Anne Hulsey on September 21, 2021, 12:24:58 AM

Hi Kuanshi,

I just discovered that the poe terms in the oq .ini files is for the investigation time, not annual. Could you update line 238 in FetchOpenQuake.py to be (something like):

1 - np.exp(-n_years * (1/RP))

Thanks,
Anne

Title: Re: openquake poe values
Post by: kuanshi on September 21, 2021, 12:28:46 AM

Hi Anne,

Sure - thank you so much for catching this!  Yes, I'll modify the FetchOpenQuake with your correction.

Thanks,
Kuanshi

Title: Re: openquake poe values
Post by: Anne Hulsey on September 21, 2021, 05:53:23 PM

Great, thanks!

After a little more trial and error, I'm still not sure what the best option is. Here is the line I used for what I described above:

'poes': 1 - np.exp(-scen_info['EqRupture']['TimeSpan']*(1.0/scen_info['EqRupture'].get('ReturnPeriod', 100)))

But it gives a lot of significant digits for 500yr RP in 50 years instead of 0.1, as I would prefer. I can't think of a good way to generically round it to the appropriate number of decimal places (perhaps to one significant digit, e.g., 0.1 and 0.002)? For now I am hardcoding 0.1 into the backend script.